Re: AMBER: problems with periodic boundary conditions

From: Michael Crowley <crowley.scripps.edu>
Date: Fri, 5 Sep 2003 07:32:29 -0700 (PDT)

Dear Monika,

The older versions of AMBER would "wrap" molecules back into the
primary cell before output and your output would look nice".
However, with that output, one could not analyze trajectory data for
transport properties like diffusion. In later versions, in particular
AMBER 7, we do not wrap molecules back into the primary cell unless
the user specifies iwrap=1 (see docs) in the input. The simulation is
correct even though one observes molecules moving out of the box
and the density appears to be decreasing at the faces. All bonded
interactions are the same, and all nonbonded interactions are still
calculated between nearest images (essentially wrapped back into the
primary cell for nonbond calcs).

Our approach is that the trajectory should be the raw MD output so
that as much analysis can be done as possible. It is not inntended that
researchers should visualize the raw data. The output can be wrapped back
into the cell without changing the raw data using ptraj with the
"image" command to produce a "pretty" output for visualization.

Sincerely,
Mike

-----------------------------------------------------------------
Physical mail: Dr. Michael F. Crowley
                 Department of Molecular Biology, TPC6
                 The Scripps Research Institute
                 10550 North Torrey Pines Road
                 La Jolla, California 92037

Electronic mail: crowley.scripps.edu
Telephone: 858/784-9290
Fax: 858/784-8688
-----------------------------------------------------------------



-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo.scripps.edu
Received on Fri Sep 05 2003 - 15:53:01 PDT
Custom Search