Re: [AMBER] aMD question: How to read the amd.log file

From: diego alonso <diegoarmino.yahoo.com>
Date: Tue, 16 Apr 2013 16:00:12 -0700 (PDT)

Please, disregard my last message. I've found the answer here http://archive.ambermd.org/201303/0027.html
But now I have another question. My amd.log file looks like this,
 
       1000      1000         1  -244983.961652883328     5150.554060385000        0.877386849692        1.000000000000      122.560169373394        0.000000000000
        1000      1000         2  -245258.279528016457     4908.733101503000        0.864339807604        0.000000000000      152.255951286500      230.274951345411
        1000      1000         3  -245595.981875224505     4786.527299220000        0.844405798738        0.000000000000      204.922763049394      435.283846299278
        1000      1000         4  -245589.265768913552     4795.329286301000        0.844216231486        0.000000000000      205.467551961184      420.517889114242

Now, as I understand it, the weighting factor should be the exponential of 1/kT times (column_8+column_9), am I right?
If that's so I'm getting some pretty gigantic numbers for exp(\beta \deltaV). For the last line it would be exp((205+420)/0.59)=10^271 which is a somewhat ridiculous number. Am I calculating the reweighting factor wrong, or should I try to reduce (a lot) the boosting potential by changing my aMD parameters?

Also, if I calculate \DeltaVp=(EthreshP-col4)^2/(alphaP+EthreshP-col4) I don't get column 8, nor do I get column 9 if I calculate \DeltaVd using column 5 and the corresponding EthreshD and alphaD parameters. For example, for row 1:

\DeltaVp= (-243830 + 244983.9)^2/(17073 - 243830 + 244983.9) = 73 != 122.56


What am I missing?
 
PS: My system consists of a 600 aminocids protein with a total 900000 atoms including water molecules. EthreshP=-243830 alphaP=17073,  EthreshD=5046 and alphaD=444 Kcal/mol.


________________________________
 From: Diego Javier Alonso de Armiño <diegoarmino.gmail.com>
To: AMBER Mailing List <amber.ambermd.org>
Sent: Tuesday, April 16, 2013 6:56 PM
Subject: [AMBER] aMD question: How to read the amd.log file
 

Dear Amber community,

I'm trying Accelerated Molecular Dynamics  (dual boost, amd=3) on my
system. I want to calculate (corrected) free energy surfaces for certain
pairs of coordinates.
According to what I've read, the reweighting procedure consists of
obtaining the histogram with the following equation

Hij = \sum_k \delta_ij Exp(\beta*\DeltaV)

where the sum is over all the frames k, \delta_ij is 0 if the frame "k" the
coordinates in question fall into bin (ij) and cero otherwise, \beta is
1/kT and \DeltaV is the biasing potential.

The problem is, the manual doesn't explain how to read the amd.log file.
What's the meaning of each of the nine columns of the amd.log file?
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Tue Apr 16 2013 - 16:30:02 PDT
Custom Search