Hi, all, there is a problem I am facing for one month. I am using CVFF to simulate my system. While the input include the electrostatic force, it is hard to minimize the energy of the system, the final energy is negative.
and there are errors when I fix system in NPT, as
out of range cannot compute PPPM.
I try to use neighbor modify to fix that, however, it did not work.
I am not sure it is because the structure issue or the coefficients for the input.
here is the input for my job:
units real
dimension 3
boundary p p p
atom_style full # include charges
pair_style lj/cut/coul/long 12 # only Lennard-Jones potential
bond_style harmonic
angle_style harmonic
dihedral_style harmonic
improper_style cvff
read_data CNTket_1.dat
kspace_style pppm 1e-6
#kspace_style none
special_bonds lj 0.0 0.0 1.0 coul 0 0 1
#newton on off
bond_coeff 2 372.4000 1.3730
angle_coeff 5 63.3100 119.9500
angle_coeff 10 63.0300 118.5800
improper_coeff 5 0.4000 -1 2
Thanks,
Haoxiang Huang