pyiron_lammps

The pyiron_lammps packages provides primarily two functions. A write_lammps_structure() function to write an ase.atoms.Atoms structure to an LAMMPS data file and a parse_lammps_output_files() function to parse the log.lammps, dump.out and dump.h5 files from the LAMMPS thermo and dump commands.

Installation

The pyiron_lammps package is distributed via both pypi: pip install pyiron-lammps and conda-forge: conda install -c conda-forge pyiron_lammps

Write LAMMPS structure

The write_lammps_structure() function is designed to write an ase.atoms.Atoms structure to an LAMMPS data file: ```python from pyironlammps import writelammps_structure

writelammpsstructure( structure, potentialelements, units="metal", filename="lammps.data", workingdirectory=None, ) `` Thestructureparameter refers to thease.atoms.Atomsstructure and thepotentialelementsrefers to the list of elements implemented in the specific interatomic potential. For example the [NiAlH_jea.eam.alloy](https://github.com/lammps/lammps/blob/develop/potentials/NiAlH_jea.eam.alloy) potential implements the elementsNi,AlandH, so when writing a structure for a simulation with this potential thepotentialelements=["Ni", "Al", "H"]. It is important to maintain the order of the elements as LAMMPS internally references the elements based on their index, starting from one. Theunitsparameter refers to the LAMMPS internal [units](https://docs.lammps.org/units.html) to convert thease.atoms.Atomsobject which is defined in Angstrom to the length scale of the LAMMPS simulation. Finally,filenameparameter and the current working directoryworking_directory parameter are designed to select the location where the LAMMPS structure should be written. With the default parameters the LAMMPS structure is written in thelammps.data` file in the current directory.

Parse LAMMPS output

In addition to writing the LAMMPS input structure pyiron_lammps also provide the parse_lammps_output_files() function to parse the LAMMPS output files, namely the log.lammps, dump.out and dump.h5 files: ```python from pyironlammps import parselammpsoutputfiles

parselammpsoutputfiles( workingdirectory, structure, potentialelements, units="metal", dumph5filename="dump.h5", dumpoutfilename="dump.out", loglammpsfilename="log.lammps", ) `` In analogy to thewritelammpsstructure()function theworkingdirectoryparameter refers to the directory which contains the output files. Thestructureparameter reefers to thease.atoms.Atomsobject which should be used as template to parse the structure from the dump files. This structure is again required as LAMMPS internally references elements only by an index, so the template structure is required to map the elements from the interatomic potential back to the elements of thease.atoms.Atomsobject. In the same way thepotentialelementsrefers to the list of elements implemented in the specific interatomic potential. Theunitsparameter refers to the LAMMPS internal [units](https://docs.lammps.org/units.html) to convert thease.atoms.Atomsobject which is defined in Angstrom to the length scale of the LAMMPS simulation. Finally, the parametersdumph5filename,dumpoutfilenameandloglammpsfile_name` refer to the output file names.

For the dump.out file the following LAMMPS dump command should be added to the LAMMPS input file: dump 1 all custom 100 dump.out id type xsu ysu zsu fx fy fz vx vy vz dump_modify 1 sort id format line "%d %d %20.15g %20.15g %20.15g %20.15g %20.15g %20.15g %20.15g %20.15g %20.15g"

For the log.lammps file the following LAMMPS thermo command should be added to the LAMMPS input file: thermo_style custom step temp pe etotal pxx pxy pxz pyy pyz pzz vol thermo_modify format float %20.15g thermo 100

Usage

Currently, the pyiron_lammps parser is primarily used in the pyiron_atomistics package and its successor the atomistics package to provide a simple LAMMPS parser. It only depends on ase, numpy, pandas and scipy and has an optional dependency on h5py to parse the LAMMPS h5md format.