Re:Re: Re: Dielectric Constant

From: ydhuang2727 (
Date: Thu Apr 29 2010 - 07:39:49 CDT

Dear Jerome,
Thanks for your kind analysis.
The dielectric constant of a protein is about 2 to 4, thus when the protein is polarizable, a proper dielectric constant bigger than 1.0 should be adopted to reduce the intensity bwteen two charged atoms. Should we?

If this is the case, i can abtain the same values as what namd got.
Appreciate your attention!
Good day!

在2010-04-29 15:59:42,"Jérôme Hénin" <> 写道:
>Dear Yangdong,
>> Based on the CHARMM22 force field, i use my own code to roughly minimize
>> a peptide and compare each term of interactions with the results of
>> minimization of namd2.6. My approachment is simply to set velocity zero
>> everay step, which is inferior to annealing algorithm that namd may adpot.
>> All interactions are comparable and consistent with the results of namd
>> except for electrostatic interaction when dielectric constant is 1.0, the
>> absolute value of which by my code is nearly twice as big as that by namd.
>> As to VDW and Electrostatic interactions, all pairs of atoms without bond
>> interaction( bond stress, angle bend, dihedral torsion, improper dihedral)
>> are considered in my code. I did not use cut-off method. Can you find
>> something wrong?
>Maybe. In the CHARMM ff, the rule is to include nonbonded terms for
>all atoms except those linked through one or two bonds (which will
>typically be coupled by bond and angle terms). Atoms connected by a
>torsion term (3 bonds) have their full nonbonded interactions (the
>so-called 1,4-interactions). There are a few exceptions: for some atom
>types, the CHARMM parameter files list modified Lennard-Jones
>parameters for 1,4 interactions. Electrostatics is never modified,

This archive was generated by hypermail 2.1.6 : Wed Feb 29 2012 - 15:55:43 CST