Re: [AMBER] pressure tensor components

From: case <case.biomaps.rutgers.edu>
Date: Sun, 13 Dec 2009 10:10:02 -0500

On Thu, Dec 10, 2009, Silvia Carlotto wrote:

> I am trying to compute the surface tension of an interface
>
> The formula for the surface tension,
> gamma, that I want to use is the following:
>
> gamma=(Lz/2)[<Pzz>-0.5(<Pxx>+<Pyy>)]
> Pzz is the normal pressure tensor component (average value for the entire
> simulation) and Pxx, Pyy are the lateral
>
> pressure tensor components.
>
> Is it possible to obtained these values directly from the simulation?

The individual pressure components are stored in the pres() array in runmd.f
(in sander). Search for "ntp > 0". You could then add print statements to
write out the information you want.

...dac


_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Sun Dec 13 2009 - 07:30:03 PST
Custom Search