I have an amorphous silica surface system with 576 atoms which was generated using BKS forcefield. I would like to relax the system with Reaxff. To do that I use the fixes in lammps as,

fix             1 mobile nve
fix             2 mobile qeq/reax 1 0.0 10.0 1.0e-6 reax/c
fix             3 mobile press/berendsen z 1.0 1.0 1000 modulus 360000
fix             4 mobile temp/berendsen 298.0 298.0 100.0

can someone point out whether I should use the pressure tensors in all lattice directions or having the z direction specified is enough? so that it physically makes sense!



You must log in to answer this question.

Browse other questions tagged .