Generation of initial molecular dynamics configurations in arbitrary geometries and in parallel

Macpherson, G.B. and Borg, M.K. and Reese, J.M. (2007) Generation of initial molecular dynamics configurations in arbitrary geometries and in parallel. Molecular Simulation, 33 (15). pp. 1199-1212. ISSN 1029-0435 (https://doi.org/10.1080/08927020701730724)

[thumbnail of strathprints005374]
Preview
 Text. Filename: strathprints005374.pdf
Accepted Author Manuscript
Download (4MB)| Preview

Abstract

A computational pre-processing tool for generating initial configurations of molecules for molecular dynamics simulations in geometries described by a mesh of unstructured arbitrary polyhedra is described. The mesh is divided into separate zones and each can be filled with a single crystal lattice of atoms. Each zone is filled by creating an expanding cube of crystal unit cells, initiated from an anchor point for the lattice. Each unit cell places the appropriate atoms for the user-specified crystal structure and orientation. The cube expands until the entire zone is filled with the lattice; zones with concave and disconnected volumes may be filled. When the mesh is spatially decomposed into portions for distributed parallel processing, each portion may be filled independently, meaning that the entire molecular system never needs to fit onto a single processor, allowing very large systems to be created. The computational time required to fill a zone with molecules scales linearly with the number of cells in the zone for a fixed number of molecules, and better than linearly with the number of molecules for a fixed number of mesh cells. Our tool, molConfig, has been implemented in the open source C++ code OpenFOAM.

CORE (COnnecting REpositories)