Hi Lammps-users,
I have been trying to quench b-cristobalite from 5000K to 300K using the BKS potential but I keep getting an error of lost atoms. I obtained the atom coordinates from Wyckoff, and generated a system of 648 atoms. There is no such issues at lower temperatures of 3000K and below. I have tried different damping from 10.0 to 0.5 but the same thing occurred.
I tried to use the same methods as mentioned in the subject thread as well, to no avail. If anyone with experience doing this could kindly advise me on what is wrong with my method, it would be greatly appreciated.
Initialise
dimension 3
boundary p p p
units metal
atom_style charge
Create geometry
read_data silica.dat
neighbor 10.0 bin
neigh_modify every 1 delay 0 check yes
Potentials
kspace_style pppm 0.00001
pair_style buck/coul/long 8.0
Si type 1, Oxy type 2
pair_coeff 1 1 0.0 0.10 0.0
pair_coeff 1 2 18003.7572 0.205205 133.5381
pair_coeff 2 2 1388.7730 0.362319 175.0
Initialise velocities
compute new all temp
compute 1 all pe/atom
compute 2 all ke/atom
compute 3 all stress/atom pair
velocity all create 5000.0 4928459 temp new
Thermostat
fix 1 all nvt temp 5000.0 5000.0 1.0
Run
timestep 0.0005
thermo 100
#Dump format: dump ID group-ID style N file args
dump 1 all custom 1000 *.dump x y z
dump_modify 1 sort id
#Number of steps to run
run 100000
fix 1 all nvt temp 5000.0 300.0 1.0
run 3000000
write_restart restart.*.heat
Regards,
Yeo Jingjie
Ph.D. Student
School of Mechanical and Aerospace Engineering
Nanyang Technological University, Singapore