Re: Kb, Ktheta values for TIP3P water model

From: Axel Kohlmeyer (akohlmey_at_gmail.com)
Date: Tue Aug 12 2014 - 12:27:25 CDT

On Tue, Aug 12, 2014 at 1:19 PM, Viswanath Pasumarthi
<v.pasumarthi_at_iitg.ernet.in> wrote:
> Yes Kenno, there is a warning message which I have failed to notice when
> force constant values for O-H bond and H-O-H angle are taken to be zero.
>
> Warning: Ignored 400 bonds with zero force constants.
> Warning: Will get H-H distance in rigid H2O from H-O-H angle.
> Warning: Ignored 200 angles with zero force constants.
>
> Moving on, as Maxim suggested, on checking the variation of length of O-H
> bond with force constants, kb = 1,2,5,10,100,450, the bond lengths varied
> from 5% to 15% on average for 10 randomly selected bonds. This only
> reinforces the failure of rigidBonds all constraint. But the second
> warning message posted above indicates that NAMD identified water as a
> rigid molecule.

the only case that i can imagine where a) you do everything correctly
and b) bonds are not fully rigid as expected, is when you run a
minimization and not an MD. but if you get as big a distortion as you
claim, then it is very likely that your overall structure, geometry or
other input has issues.

there is (still) far too little information to say anything definite, tho.

axel.

>
> The following is an excerpt from output log file for a system involving
> ionic liquid+water mixture. 20 IL molecules + 200 H2O molecules were
> considered. It can be seen that 900 rigid bonds are identified. Each IL
> molecule contains 15 covalent bonds involving hydrogen, which makes 300
> rigid bonds corresponding to IL. The rest 600, contributes from 200 water
> molecules, which implies even H-H is considered as an explicit bond for
> water molecule. Is it alright and does it mean the bond parameters for H-H
> should also be provided?
>
> Info: STRUCTURE SUMMARY:
> Info: 1400 ATOMS
> Info: 1180 BONDS
> Info: 1600 ANGLES
> Info: 1660 DIHEDRALS
> Info: 0 IMPROPERS
> Info: 0 CROSSTERMS
> Info: 0 EXCLUSIONS
> Info: 900 RIGID BONDS
> Info: 3297 DEGREES OF FREEDOM
> Info: 700 HYDROGEN GROUPS
> Info: 4 ATOMS IN LARGEST HYDROGEN GROUP
> Info: 700 MIGRATION GROUPS
> Info: 4 ATOMS IN LARGEST MIGRATION GROUP
> Info: TOTAL MASS = 11990.3 amu
> Info: TOTAL CHARGE = 18.0571 e
> Info: MASS DENSITY = 0.73744 g/cm^3
> Info: ATOM DENSITY = 0.0518519 atoms/A^3
>
>
> PSF:
> 801 MOL 1 UNK O OW -0.928231 15.9994 0
> 802 MOL 1 UNK H HT 0.464116 1.00794 0
> 803 MOL 1 UNK H HT 0.464116 1.00794 0
>
>
> BONDS
> !atom type Kb b0
> !TIP3P WATER
> OW HT 450.0000 0.9572 !
>
> ANGLES
> !atom types Ktheta Theta0 Kub S0
> !TIP3P Water
> HT OW HT 55.0000 104.52 !
>
> NONBONDED
> !atom ignored epsilon Rmin/2 ignored eps,1-4 Rmin/2,1-4
> !TIP3P Water
> HT 0.000000 -0.046000 0.224500 !
> OW 0.000000 -0.152100 1.768200 !
>
>
>
>> On 08/11/2014 01:56 PM, Viswanath Pasumarthi wrote:
>>> considering a value of k=0 mean absence of spring (physical
>>> bond) between the respective atoms. It would have been easier to
>>> understand this mistake had there been an error message pointing to the
>>> choice of k value as zero.
>>
>> But there *is* a message, and if you would have noticed/posted it, we
>> wouldn't have needed to have this long, marginally constructive
>> discussion. Look in your NAMD output for the following:
>>
>> Warning: Ignored .* bonds with zero force constants.
>> Warning: Ignored .* angles with zero force constants.
>>
>> Where .* is an integer.
>>
>>> However, using rigidBonds constraints still
>>> required bond constraint values to be specified in the parameter file.
>>> This, as Kenno cautioned has turned into a flexible simulation, as
>>> covalent bonds involving hydrogen became more flexible at lower k values
>>> (1,2,5,10) and functioning decent at the values reported in modified
>>> TIP3P
>>> model.
>>
>> I have to repeat that I tried and failed to reproduce this; my rigidBonds
>> simulations behave correctly for all force constants (even absurd values
>> such as 0.0001 or -10000) except zero. Most probably, there is an issue
>> with your psf and/or rigidBonds is not applied correctly. Which I have
>> been saying from the very beginning!
>>
>> You probably need to try to use more detailed language (for example
>> "distances" or "force constants" rather than "values") when describing
>> your problem, and attach relevant input and output files to your posts,
>> else we'll never get anywhere.
>>
>>
>

-- 
Dr. Axel Kohlmeyer  akohlmey_at_gmail.com  http://goo.gl/1wk0
College of Science & Technology, Temple University, Philadelphia PA, USA
International Centre for Theoretical Physics, Trieste. Italy.

This archive was generated by hypermail 2.1.6 : Wed Dec 31 2014 - 23:22:43 CST