Solvent Heating with Rigid Bonds

From: ben estimated.com (ben_at_estimated.com)
Date: Sun Jan 26 2020 - 22:07:34 CST

• Next message: Laura Lopes: "problem with TIP4P on nightly version"
• Previous message: Shadi Rahnama: "Why very large force reported at the starting point of constant velocity SMD run in NAMD log file?"
• Messages sorted by: [ date ] [ thread ] [ subject ] [ author ] [ attachment ]

I'm heating my system from 200K to 310K by coupling a heat bath to the only the solvent atoms. I did this by using NAMD's "tCouple" and setting the B value of all the water molecules in the PDB file to 1. I am also using a 2fs timestep and therefore have rigidbonds to "all." If the water molecules are unable to vibrate/stretch because of this, does it still make sense to heat them and expect my simulation to be physiologically realistic? (I assume this problem would still exist even if I set rigidbonds to "water").

--Thank you,

Ben

• Next message: Laura Lopes: "problem with TIP4P on nightly version"
• Previous message: Shadi Rahnama: "Why very large force reported at the starting point of constant velocity SMD run in NAMD log file?"
• Messages sorted by: [ date ] [ thread ] [ subject ] [ author ] [ attachment ]

This archive was generated by hypermail 2.1.6 : Fri Dec 31 2021 - 23:17:08 CST