AMBER: Radial Distribution functions by ptraj

From: Priti Hansia <>
Date: Fri, 2 Dec 2005 01:14:02 +0530 (IST)

Dear Amber users,

I am trying to calculate rdf function by using ptraj. Here is my input file:

trajin md1.mdcrd.gz
trajin md2.mdcrd.gz
trajin md3.mdcrd.gz
trajin md4.mdcrd.gz
trajin md5.mdcrd.gz
radial rdf 0.1 10.0 :WAT :1-94

It reads the files alright. However, after processing trajin md2.mdcrd.gz
it gives segmentation fault !

I have gone through the amber archive carefully and have noticed that some
other users also have mentioned this. However, I could not find out why it
is happening. Am I wrong while mentioning the solute and solvent mask?
Ptraj with my input gives following:

      spacing is 0.100
      maximum is 10.000
      [This implies 99 bins]
      density is 0.033
      RDF will bin all solute-mask to solvent-mask distances
      Solvent atom selection is :95-4589
      Solute atom selection is :1-94

My system has 94 amino acids and 95 to 4589 are water molecules. Does
'radial rdf 0.1 10.0 :WAT :1-94' mean that RDF by ptraj will represent the
interactions of all the atoms in 94 residues (~1500) with all the water
molecules (~4500) ? What will be my input if I want to calculate RDF for a
single atom or a single residue?

Any help will be of great help to me.


The AMBER Mail Reflector
To post, send mail to
To unsubscribe, send "unsubscribe amber" to
Received on Thu Dec 01 2005 - 19:53:00 PST
Custom Search