From: Francesco Pietra (chiendarret_at_gmail.com)
Date: Sun Jun 02 2013 - 08:45:04 CDT
I have a problem with a water molecule bound to Zn(II) in a metal cluster
of a protein, under parm7 amber ff. As I posted before (unanswered), by
setting "rigidBonds water", NAMD 2.9-CUDA hangs, asking for bond parameters
for the two H-atoms of the water molecule bound to Zn(II). Although the
whole cluster, such water included, is called MOL.
By not freezing any bond, the simulation goes on smoothly, albeit slowly.
Therefore, I have now tried to retain all bond motions by the use of MOLLY,
of which I had no previous experience, with multiple timesteps. Got the
same error message, asking for params as said above. In all cases,
"timestep 1.00" The relevant portion of the conf file:
colvars on # provide a colvarsConfig (and colvarsINput if restart)
## SIMULATION PARAMETERS
## AMBER FF settings
# rigidBonds water
# rigidTolerance 1.0e-8
## Integrator Parameters
nonbondedFreq 2 # ts between nonbonded eval
fullElectFrequency 6 # number ts between full electrostatic eval
Which error related to MOLLY in conf file? If none, I wonder whether MOLLY
really retains bond motions for water, or the way MOLLY does that is
different by allwing all bond motions by the traditional flags.
This archive was generated by hypermail 2.1.6 : Tue Dec 31 2013 - 23:23:17 CST