restraining

From: Christopher Rowley (cnrowley_at_mun.ca)
Date: Thu Aug 23 2012 - 13:55:59 CDT

Hi,

I'm using NAMD for simulations of a model for hydrogen sulphide (H2S).
The model is defined as being rigid, like TIP3P, so I've included the
command "rigidbonds all" Nevertheless, the H-S-H bond is flexible in MD
simulations.

The PSF file contains a bond between the two hydrogens, but this doesn't
seem to make a difference.

Unless I've misunderstood, H-H bonds should be identified in Molecule.C

Molecule.C:
           if ( ! is_water(a2) ) { // H-H but not water
             rigidBondLengths[a1] = ( mode == RIGID_ALL ? x0 : 0. );
             rigidBondLengths[a2] = ( mode == RIGID_ALL ? x0 : 0. );
           }

Is there a way to force NAMD to constrain this distance/angke?

Thanks,
Chris

-- 
Christopher Rowley
Ph.D. (chemistry)
Assistant Professor
Department of Chemistry
Memorial University of Newfoundland
http://www.chem.mun.ca/homes/cnrhome/
This electronic communication is governed by the terms and conditions at
http://www.mun.ca/cc/policies/electronic_communications_disclaimer_2012.php

This archive was generated by hypermail 2.1.6 : Mon Dec 31 2012 - 23:21:56 CST