We are currently investigating the impact resistance of multilayers using molecular dynamics. This project showed us that modelling material interfaces could be challenging even for someone with a background in atomic simulations. Therefore, we would like to share some useful information and the LAMMPS input files to model the aluminum-polyurethane system shown below. A movie of the simulation is vailable here.
The command create_atoms in LAMMPS was used to model aluminum. In this tutorial, Mark Tschopp outlines how to use the create_atoms command to generate a crystal containing a grain boundary. I used the EAM potential by Adams and Ercolessi to model aluminum.
Modelling a polymer can be challenging depending on the complexity of an individual polymer chain. The structure of the polyurethane chain shown in the above video can be depicted as below.
I created the polymer chain using the Enhanced Monte Carlo Package by Pieter J. in ’t Veld. Moltemplate by Andrew Jewett is another useful tool to build polymers. I used the polymer consistent force field (PCFF) for modelling polyurethane.
Mark Tschopp has provided an overview of modelling of a polymer chain and an amorphous polyethylene sample in LAMMPS. Moreover, Appendix A of the reference [1] and the reference [2] provide some useful information about modelling polymers.
The two data files (i.e. polyurethane and aluminum) can be combined using MATLAB or even MS Excel. When combining the two models, careful attention needed to be paid to the image flags of the atoms in the polymer to make sure that they do not have crossed the periodic boundary. The two potentails used for the simulation (e.g., PCFF and EAM) must be in similar LAMMPS units. In this simulation, the non-bonded interactions between aluminum and polyurea were modelled using the van der Waals parameters available in the interface force field.
I used OVITO to visulalize MD tragectories of this simulation. A great introductory tutorial about OVITO is available here. VMD is another very useful visualization tool.
Many simulation questions have been answered in the LAMMPS mail list. If you would like to learn the basics of atomic simulations, the book Modeling Materials: Continuum, Atomistic and Multiscale Techniques could be useful.
I hope that this information and the LAMMPS files will be useful for those who are interested in modelling material interfaces. Good luck!
[1] A.P. Awasthi, D.C. Lagoudas, D.C. Hammerand, Modeling of graphene–polymer interfacial mechanical behavior using molecular dynamics, Model. Simul. Mater. Sci. Eng. 17 (2009) 015002.
[2] D. Hossain, M.A. Tschopp, D.K. Ward, J.L. Bouvard, P. Wang, M.F. Horstemeyer, Molecular dynamics simulations of deformation mechanisms of amorphous polyethylene, Polymer 51 (2010) pp. 6071-6083.