Amber masthead
Filler image AmberTools24 Amber24 Manuals Tutorials Force Fields Contacts History
Filler image

Useful links:

Amber Home
Download Amber
Installation
Amber Citations
GPU Support
Updates
Mailing Lists
For Educators
File Formats
Contributors
Workshops
Sampling Configuration Space
 

(Note: These tutorials are meant to provide illustrative examples of how to use the AMBER software suite to carry out simulations that can be run on a simple workstation in a reasonable period of time. They do not necessarily provide the optimal choice of parameters or methods for the particular application area.)
Copyright Dwight McGee, Hailey Bureau, Caley Allen, Rigoberto Hernandez

Section 3: Analysis of Results


Now that the ASMD simulations are complete, we need to restitch the PMF output from each stage created by ASMD.py to obtain the actual PMF for each of the simulations at different velocities.

Restitch the PMF


Lets first begin with the 100 Å/ns simulations. To do this we can use the python script create.PMF.py. To get a full a list of all the available options use the -h or --help flags. In order for the script to work, all input files being stitchd must have the following naming convention jar.stage#.dat where # is the number of the stage.

Type the following command: to produce the output file jar.PMF.100ang.25tps.dat

>create.PMF.py -s 5 -o jar.100ang.25tps.dat

Next, create the PMF for the 10 Å/ns simulation to produce the output file jar.PMF.10ang.25tps.dat

>create.PMF.py -s 5 -o jar.10ang.25tps.dat

Plotting the PMF


You can use your favorite program to plot the PMFs.

For example, if you were to use the GNU graphing suite, you can type the following from the terminal command line: gnuplot
gnuplot> plot "jar.PMF.100ang.25tps.dat" w l ti "100 Å/ns" lw 3 lc rgb "black", "jar.PMF.10ang.25tps.dat" w l ti "10 Å/ns" lw 3 lc rgb "blue"

As expected, the PMF for the faster velocity is higher than the slower velocity pull at the end of the pull. At faster velocities, the system has less time to adjust to the pertubation being applied, thus resulting in a higher PMF.

Things to consider: What would happen if you: employed a different solvent model, added more trajectories, pulled at a slower velocity, increased or decreased the number of stages, or pulled for a longer distance?



CLICK HERE TO GO TO BACK TO THE INTRODUCTION