I have a system of Peptide-popc-ions-water to simulate. I'm trying to run MD 
simulation following the Justin's tutorial of kalp15-DPPC. 

The force field which I use, is charmm36. 

In NVT step, I get LINCS-warning as follow:
Step 0, time 0 (ps)  LINCS WARNING
relative constraint deviation after LINCS:
rms 2.460224, max 54.696251 (between atoms 4719 and 4716)
bonds that rotated more than 30 degrees:
 atom 1 atom 2  angle  previous, current, constraint length

Step 0, time 0 (ps)  LINCS WARNING
relative constraint deviation after LINCS:
rms 0.182707, max 15.613569 (between atoms 4725 and 4726)

When I checked the warning, I see the first warning is related to 2 atoms of 
71POPC, the chain of 71POPC residue is broken near these two atoms. And the 
other warning is related to a bond between H and C in the other POPC residue.

Some settings of nvt.mdp is as:
ns_type        = grid        
nstlist        = 5      
rlist        = 1.2        
rlistlong       = 1.4
rcoulomb    = 1.2       
rvdw        = 1.2        
pbc        = xyz   
vdwtype         = switch
rvdw_switch     = 0.8
; Parameters for treating bonded interactions
continuation    = no      
constraint_algorithm = LINCS    
constraints    = all-bonds   
lincs_iter    = 1        
lincs_order    = 4        

Do I need to increase the rlistlong? I have done a simulation of POPC-water 
before with these settings, but everything was ok. So what's the problem? 
Anybody may help me?

Your suggestions would be appreciated. Thanks in advance.

gmx-users mailing list    gmx-users@gromacs.org
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the
www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Reply via email to