Cell relax in LAMMPS


These days I am learning how to use LAMMPS.

To allow the simulation cell to relax, you need to specify box/relax with fix command as, :

fix 1 all box/relax tri 0.0

Here it mentions that the cell is triclinic, which meams x, y, z, yz, zx, and xy are controlled independently. And final value 0.0 means pressure to be achieved.

Then you minimize the system as, :

minimize    1.0e-7 0.001 1000 10000

Here the values indicate etol, ftol, maxiter, and maxeval respectively.

For example, if you use MEAM potential with meam.library and meam.coeff files. The LAMMPS input file can be as follows, :

units metal 
boundary p p p 
atom_modify sort 0 0.0 

read_data data_lammps

### interactions 
pair_style meam 
pair_coeff * * meam.library Al Mg Si meam.coeff Al Mg Si 

### run
fix 1 all box/relax tri 0.0
dump dump_all all custom 1 trj_lammps* id type x y z vx vy vz fx fy fz
thermo_style custom step temp press cpu pxx pyy pzz pxy pxz pyz ke pe etotal vol lx ly lz atoms
thermo_modify flush yes
thermo 1
minimize    1.0e-7 0.001 1000 10000