RE: AMBER: Calculating dipole moments

From: Ross Walker <ross.rosswalker.co.uk>
Date: Thu, 29 Nov 2007 08:08:50 -0800

Hi Austin,

> the "CHARGE" and "PERT_CHARGE" sections, I get the charges for the
> unperturbed and perturbed residues, respectively. Seems
> likely that the
> charges employed in the dipole moment calculations are the from the
> PERT_CHARGE section ie using perturbed charges for the residue I am
> mutating, and unperturbed charges for the rest of the
> trisaccharide. In
> this case, what will be the best way to proceed? Extract
> structures from
> the TI calculation, generate a new unperturbed prmtop file, and then
> compute the average dipole moment for each structure?

If you are happy with scripting I believe this will be the fastest way to
go. You could write a script to take each structure in the mdcrd file - have
ptraj make a restart / inpcrd file then use this with the regular prmtop
file to run sander set to 1 step of MD and calculate the dipoles and then
use grep to extract the data you need and append it to a results file.

Anything beyond this will require modification to the code and Amber 8 is
not really maintained anymore - especially for TI calculations since the
implementation of this has change significantly in Amber 9. Do you have
access to Amber 9? If so then you could try this, I suspect that the dipoles
printing will have more chance of working correctly in Amber 9 since it uses
regular prmtop files instead of perturbed ones - hence the main charge array
(used by the dipole printing routine) should contain the real charges from
the prmtop file.

Good luck
Ross

/\
\/
|\oss Walker

| HPC Consultant and Staff Scientist |
| San Diego Supercomputer Center |
| Tel: +1 858 822 0854 | EMail:- ross.rosswalker.co.uk |
| http://www.rosswalker.co.uk | PGP Key available on request |

Note: Electronic Mail is not secure, has no guarantee of delivery, may not
be read every day, and should not be used for urgent or sensitive issues.


-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Sun Dec 02 2007 - 06:07:21 PST
Custom Search