Hello everyone, My lammps input script is as following :

units real variable cfac equal 1.01325e-4 variable cunits string GPa

atom_style full boundary p p p

pair_style lj/cut/coul/cut 20

kspace_style pppm 1.0e-4 #kspace_modify mesh 4 4 4 bond_style harmonic angle_style harmonic dihedral_style charmmfsw improper_style cvff

read_data tkx_50_supercell_data.txt

include pair_lj_coul.dat

neighbor 2.0 nsq neigh_modify delay 0 every 1 check yes one 100000 page 1000000

Define minimization parameters

variable etol equal 0.0 variable ftol equal 1.0e-10 variable maxiter equal 1000 variable maxeval equal 1000 variable dmax equal 1.0e-2

Setup minimization style

min_style cg min_modify dmax ${dmax} line quadratic

######################################

EQUILIBRATION

timestep 0.001 velocity all create 300 12345 mom yes rot no fix 1 all npt temp 300 300 1 iso 0 0 1 drag 1

Set thermo output

thermo 1000 thermo_style custom step lx ly lz press pxx pyy pzz pe temp dump trj all atom 1000 Au.lammpstrj

Run for at least 10 picosecond (assuming 1 fs timestep)

run 50000

Whenever I am trying to run the simulation, I am facing the error: must redefine kspace style after changing to triclinic box. What should be done here?

More Md. Rifat Hossain's questions See All
Similar questions and discussions