Hi Neha,
PBCAL and the "corrected reaction field energy" are the same but in different units.
PBCAL is given in kcal/mol, the "corrected reaction field energy" in kT/mol.
This is done to be compatible with the DELPHI program since it calculates the reaction field energy in kT/mol as well.
To calculate the "corrected reaction field energy" from PBCAL you simply need to multiply PBCAL with (1000.0 * 4.184 / (T * R)) with the absolute temperature T and the gas constant R. You can reproduce that from the source code if you grep for "kcal2kt" in $AMBERHOME/src/mm_pbsa/mm_pbsa_calceneent.pm and for "kt2kcal" in $AMBERHOME/src/mm_pbsa/mm_pbsa_statistics.pm, respectively.
Cheers!
Hannes
-------- Original-Nachricht --------
> Datum: Fri, 13 Jun 2008 13:42:54 +0800
> Von: "Neha Gandhi" <n.gandhiau.gmail.com>
> An: amber.scripps.edu
> Betreff: AMBER: PCAL and polar contribution in the output file
> Hi All,
>
> In MM-pbsa analysis, the ouput consists of corrected reaction field
> energy. I am able to calculate non-polar contribution which is
> PBSUR=surface area*0.00542+0.92. Where is the electrostatic solvation
> energy in the output file such as snapshot_rec.all.out? I cannot
> correlate PBCAL with corrected reaction field energy. I appreciate if
> somebody can help me to calculate PBSOL for each snapshot. How PBCAL
> and corrected reaction field energy are related?
>
>
> --
> Regards,
> Neha Gandhi,
> School of Biomedical Sciences,
> Curtin University of Technology,
> GPO Box U1987 Perth,
> Western Australia 6845
> -----------------------------------------------------------------------
> The AMBER Mail Reflector
> To post, send mail to amber.scripps.edu
> To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
> to majordomo.scripps.edu
--
GMX startet ShortView.de. Hier findest Du Leute mit Deinen Interessen!
Jetzt dabei sein: http://www.shortview.de/?mc=sv_ext_mf.gmx
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
to majordomo.scripps.edu
Received on Wed Jun 18 2008 - 06:07:09 PDT
