AMBER Archive (2005)

Subject: AMBER: About WHAM program

From: Atsutoshi Okabe (okabe_at_cherry.bio.titech.ac.jp)
Date: Sat Nov 05 2005 - 22:16:43 CST


Hello,

 Thank you for advising me on WHAM program.
I understood hist_min hist_max and loc_win_min .
Then I prepared metadatafile(test.out) and timeseriesfile , and I carried
out the next wham command

> ./wham 2 20 180 0 298 0 test.out TEST

Then wham produced output file(TEST) resembling

#Coor Free +/- Prob +/-
2.050000 nan nan nan nan
2.150000 nan nan nan nan
2.250000 nan nan nan nan

Why Free energy column and the others are $B!H(Bnan$B!I(B?

Could you please give me a suggestion?

Thank you
A. Okabe

-----Original Message-----
From: owner-amber_at_scripps.edu [mailto:owner-amber_at_scripps.edu] On Behalf Of
Guanglei Cui
Sent: Wednesday, November 02, 2005 10:47 PM
To: amber_at_scripps.edu
Subject: Re: AMBER: About WHAM program

Hi,

You don't need to manipulate DUMPAVE except perhaps discarding sampling
that's not equilibrated.

hist_min and hist_max define the range of the distance in which the
potential of mean force is calculated. loc_win_min is the center of your
umbrella window, i.e. the value you used for r2 and r3.

Regards,

Guanglei

Atsutoshi Okabe wrote:
> Dear all
>
> I want to calculate potentials of mean force(for distance,2$B"r!A#2#0"r(B 0.1
$B"r(B
> /bin) using umbrella sampling in Amber8.
>
> So I could get a DUMPAVE file.(The first field is each time point, the
> second field is
>
> the values of the reaction coordinate)
>
> But I can't understand how to manipulate the DUMPAVE file in WHAM program.
>
> I read the WHAM manual written by Pro. Alan Grossfield , but I can't
> understand.
>
> I can't also understand "wham" command. What does "hist_min and hist_max"
> mean
>
> in my case ? Doest it mean that "hist_min = 2, hist_max = 20" ?
>
> Adittionally I can't understand metadatafile formats. What does
> "loc_win_min"
>
> mean in my case ?
>
>
>
> Could you please give me a suggestion?
>
>
>
> Thank you
>
> Atsutoshi Okabe

-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu