LUNAR: Automated Input Generation and Analysis for Reactive LAMMPS Simulations

Generating simulation-ready molecular models for the LAMMPS molecular dynamics (MD) simulation software package is a difficult task and impedes the more widespread and efficient use of MD in materials design and development. Fixed-bond force fields generally require manual assignment of atom types, bonded interactions, charges, and simulation domain sizes. A new LAMMPS pre- and postprocessing toolkit (LUNAR) is presented that efficiently builds molecular systems for LAMMPS. LUNAR automatically assigns atom types, generates bonded interactions, assigns charges, and provides initial configuration methods to generate large molecular systems. LUNAR can also incorporate chemical reactivity into simulations by facilitating the use of the REACTER protocol. Additionally, LUNAR provides postprocessing for free volume calculations, cure characterization calculations, and property predictions from LAMMPS thermodynamic outputs. LUNAR has been validated via building and simulation of pure epoxy and cyanate ester polymer systems with a comparison of the corresponding predicted structures and properties to benchmark values, including experimental results from the literature. LUNAR provides the tools for the computationally driven development of next-generation composite materials in the Integrated Computational Materials Engineering (ICME) and Materials Genome Initiative (MGI) frameworks. LUNAR is written in Python with the usage of NumPy and can be used via a graphical user interface, a command line interface, or an integrated design environment. LUNAR is freely available via GitHub.