AMBER Archive (2002)

Subject: ptraj rdf output

From: Prabal Maiti (maiti_at_wag.caltech.edu)
Date: Fri Sep 20 2002 - 18:24:35 CDT


Dear Amber users:
I tried to calculate radial distribution function between water and Na+
and Na+ and DNA (this may not be so well defined). For each of this I
get two output, one with _standard extension and another with volume
extension.
Both has two columns of output apart from the 1st column with
distance(bin). Can someone explain me what are these numbers in these 2
columns. Also I would have exptected the g(r) eventually goes to 1 if it
is proporely normalized. But the numbers in 2nd and 3rd columns are
different than that.

Can someome please explain me these numbers,

>From the code I see that _standarrd is Allen & Tildesley style histogram
    fprintf(outFile, "%8.3f %8.3f %8.3f\n",
            ((double) i + 0.5) * radialInfo->spacing,
            hist, occupancy);
In this case hist should apporach 1 at large distance which is not
the case. it ranges from 0.8 to 5. What
exactly is the occupnacy means. IS it the number of molecules within
a given distnace?

Here is my ptraj input
trajin md_prod_final.trj.gz 590 690 1
radial rdf_DNA-WAT 0.05 8.0 :WAT :1-196
radial rdf_DNA-NA 0.05 8.0 :Na+ :1-196
go

Thanks
Praabal

-- 
------------------------------------------------------------------------
Prabal K Maiti                             Email:maiti_at_wag.caltech.edu
Beckman Institute (139-74)                       pkmaiti_at_yahoo.com
California Institute of Technology         Tel:626-395-8151 (O)
400 South Wilson Avenue                       :626-744-1363 (H)                    
Pasadena, CA 91125                         Fax:626-585-0918
                                           Home add: 110 S. Michigan Ave.
http://www.wag.caltech.edu/home/maiti      Apt # 19, Pasadena, CA 91106
------------------------------------------------------------------------