do you have any confirmation for
that from some published literature?
cheers,
axel.
Hi, Axel,
In Chenoweth, van Duin, and Goddard’s J. Phys. Chem. A 2008, 112, 1040-1053 paper, it was mentioned on page 1042 that “During the equilibration simulations, the C-O and H-O bond parameters were switched off to prevent reactions from occurring during the equilibration of the system. The equilibrated systems were then used in NVT-MD simulations where the temperature was set to 2500 K using a temperature damping constant of 0.5 ps and a MD time step of 0.1 fs.”
Can I do something similar with LAMMPS/ReaxFF or LAMMPS/ReaxFF-c?
do you have any confirmation for
that from some published literature?
cheers,
axel.
Hi, Axel,
In Chenoweth, van Duin, and Goddard's J. Phys. Chem. A 2008, 112, 1040-1053
paper, it was mentioned on page 1042 that "During the equilibration
simulations, the C-O and H-O bond parameters were switched off to prevent
reactions from occurring during the equilibration of the system. The
equilibrated systems were then used in NVT-MD simulations where the
temperature was set to 2500 K using a temperature damping constant of 0.5 ps
and a MD time step of 0.1 fs."
Can I do something similar with LAMMPS/ReaxFF or LAMMPS/ReaxFF-c?
you can certainly modify the parameter files according to those
instructions. this doesn't require any modification of LAMMPS itself.
There is no way to do this in "user mode." You would either have to modify
the code (more hard) or modify the ffield.reax file (less hard). Either way,
you will have to understand the structure of the ReaxFF potential to get the
desired behavior.