Hello All,
I am trying to load a crystal to fracture under uniaxial tension in LAMMPS version dated 23 Jun 2022. The simulation proceeds without errors until the stress-strain curve flattens and starts to assume a negative slope, at which point it exits due to one of the two following errors:
- ERROR: Lost atoms: original 53125 current 53123
OR - ERROR on proc 217: Out of range atoms - cannot compute PPPMDisp
Due to these errors, I cannot take the system all the way to fracture. Furthermore, at these final stages of the simulation (i.e., once the stress-strain curves flattens and starts to assume a negative slope), in addition to these errors, I sometimes observe that the temperature suddenly drops all the way to zero despite being in the NPT ensemble at 300 K.
Any suggestions for extending the simulation all the way to fracture is appreciated.
INPUT (units metal used here):
Initialization
read_restart restart.file
Potential Parameters
pair_style buck/long/coul/long long long 12
kspace_style pppm/disp 1e-5
kspace_modify force/disp/real 4.335e-6
kspace_modify force/disp/kspace 8.671e-5
“list of pair coefficients here”
Simulation Specifics
variable atm equal 0
timestep 1e-3 # in PICOSECONDS
neigh_modify every 5 delay 0 check yes
restart 20000 restart.file
fix 1 all npt temp 300 300 0.1 x {atm} {atm} 1 y {atm} {atm} 1
fix 2 all deform 1 z erate 0.0004 units box remap x
variable stp equal “step”
variable len_z equal “lz”
variable tmpr equal “temp”
variable plat_x equal “pxx/10000”
variable plat_y equal “pyy/10000”
variable pz equal “pzz/10000”
fix dp all print 200 “{stp} {len_z} {plat_x} {plat_y} {pz} {tmpr}” file opdat screen no
thermo 200
thermo_style custom step lz press temp
compute p all stress/atom NULL
dump data all custom 200 dump.data id type mass x y z vx vy vz c_p[*]
run 480000