Accelerated MD: question about dual boost mode

From: patrick wintrode (pat_wde2_at_yahoo.com)
Date: Thu Mar 31 2011 - 15:01:49 CDT

• Next message: flavio seixas: "angles"
• Previous message: Christian Jorgensen: "vmd-l: No probability distribution plot in ParseFEP"
• Next in thread: Jeff Wereszczynski: "Re: Accelerated MD: question about dual boost mode"
• Reply: Jeff Wereszczynski: "Re: Accelerated MD: question about dual boost mode"
• Messages sorted by: [ date ] [ thread ] [ subject ] [ author ] [ attachment ]

Hi.

I've been testing out the new accelerated MD functionality in NAMD. I find that applying the boost just to the dihedrals works fine. However, if I use the dual boost mode, even quite small boost energies for the total-dihedral potential causes the simulation to become unstable and crash due to a constraint failure almost instantly (within ~15 steps).

My system is a 372 residue protein and the total system size protein+solvent is ~45,000 atoms.

The dihedral boost is E = 600 and alpha = 120. The smallest total-dihedral boost I've tried so far is E = 10 and alpha = 2.

Is my total-dihedral boost just too large? From reading the McCammon group's papers, it doesn't appear to be.

Any suggestions would be appreciated.

Thanks.

Patrick L. Wintrode
Assistant Professor
Dept. of Physiology & Biophysics
Case Western Reserve University
Cleveland, Ohio 44106

• Next message: flavio seixas: "angles"
• Previous message: Christian Jorgensen: "vmd-l: No probability distribution plot in ParseFEP"
• Next in thread: Jeff Wereszczynski: "Re: Accelerated MD: question about dual boost mode"
• Reply: Jeff Wereszczynski: "Re: Accelerated MD: question about dual boost mode"
• Messages sorted by: [ date ] [ thread ] [ subject ] [ author ] [ attachment ]

This archive was generated by hypermail 2.1.6 : Mon Dec 31 2012 - 23:20:03 CST