Dear all,
I am trying to set calculations for a dna molecule by treating each nucleobase as a separate rigid body.
For this reason I use the configuration file below:
/********************************/
units real
atom_style hybrid full dipole
bond_style harmonic
angle_style harmonic
dihedral_style opls
improper_style harmonic
pair_style hybrid/overlay lj/class2 10.0 coul/cut 10.0
dielectric 1
read_data 10bp_0.topo
pair_coeff * *12 lj/class2 1.0 1.0
pair_coeff * * coul/cut
neighbor 2.0 bin
timestep 1.0
thermo_style multi
thermo 1
group rigids type 4 5 6 7 8 9 10 11 13
group no-rigids subtract all rigids
#delete_bonds rigids multi
#neigh_modify exclude molecule rigids
fix 1 no-rigids nvt temp 1.0 1.0 5.0
fix 2 rigids rigid/nvt molecule temp 1.0 1.0 5.0
run 0
/*****************************************/
So the types 4 5 6 7 8 9 10 11 13 are the types of the atoms in the nucleobase => rigids group
and the rest of the nucleoside => no-rigids group
I have to fixes. The first one responsible to integrate the no-rigids group works with no problem if it is the only fix. When I add the next fix (2) I take the problem below.
//
LAMMPS (5 May 2012)
Scanning data file …
2 = max bonds/atom
3 = max angles/atom
5 = max dihedrals/atom
1 = max impropers/atom
Reading data file …
orthogonal box = (-236.84 -241.014 -229.768) to (264.212 268.385 257.139)
1 by 1 by 1 MPI processor grid
230 atoms
148 bonds
204 angles
248 dihedrals
10 impropers
Finding 1-2 1-3 1-4 neighbors …
3 = max # of 1-2 neighbors
4 = max # of 1-3 neighbors
8 = max # of 1-4 neighbors
10 = max # of special neighbors
170 atoms in group rigids
60 atoms in group no-rigids
20 rigid bodies with 170 atoms / which is correct since I have 20 nucleosides Also I checked and the number of atoms is correct */
Segmentation fault
/******/
I would appreciate any suggestion about what is going wrong.
Kind regards,
Chris