AMBER Archive (2008)Subject: AMBER: Radial Distribution Function
From: prabhakar g (prabha.kaist_at_gmail.com)
Date: Mon Dec 08 2008  20:37:52 CST
Dear Amber experts,
I have carried out a 1ns restrained md simulation of I2(iodide) molecule in
MEOHBOX using ff99 force field,
(force constant of 100 kcal/mol)
Here is my simulation protocol using sander(Amber9)(from tutorial)
Minimization I (restrained minimization)
Minimization II ( whole system minimization)
Heating 20ps(NVT)
equilibration 300ps (NPT)
production 1ns (NVT)
am interested in calculating the radial distribution function(rdf) for
studies related to photodissociation.
I calculated the rdf with ptraj but for some of the atomic pairs i am
getting an unusual rdf
that is the curve does not converges to 1 for some of the atomic pairs (
file attached)
here is my ptraj input
trajin i2_meoh_md1.2ns.mdcrd
radial rdfC1I1.out 0.1 15.0 :2569_at_C1 :1_at_I1 density 0.015036
radial rdfC1I2.out 0.1 15.0 :2569_at_C1 :1_at_I2 density 0.015036
radial rdfO1I1.out 0.1 15.0 :2569_at_O1 :1_at_I1 density 0.015036
radial rdfO1I2.out 0.1 15.0 :2569_at_O1 :1_at_I2 density 0.015036
radial rdf1HCI1.out 0.1 15.0 :2569_at_1HC :1_at_I1 density 0.015036
radial rdf1HCI2.out 0.1 15.0 :2569_at_1HC :1_at_I2 density 0.015036
radial rdf2HCI1.out 0.1 15.0 :2569_at_2HC :1_at_I1 density 0.015036
radial rdf2HCI2.out 0.1 15.0 :2569_at_2HC :1_at_I2 density 0.015036
radial rdf3HCI1.out 0.1 15.0 :2569_at_3HC :1_at_I1 density 0.015036
radial rdf3HCI2.out 0.1 15.0 :2569_at_3HC :1_at_I2 density 0.015036
radial rdf1HOI1.out 0.1 15.0 :2569_at_1HO :1_at_I1 density 0.015036
radial rdf1HOI2.out 0.1 15.0 :2569_at_1HO :1_at_I2 density 0.015036
whereas the same simulation in moldy(http://www.ccp5.ac.uk/moldy/moldy.html)
I get a very gud radial distrubution function. my aim was to scale the
computing time of both the programs ..
could anyone help regarding this..
thanks in advance
prabhakar

The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
to majordomo_at_scripps.edu
