This week in pchem we are discussing the energy minimization techniques that are used in computational chemistry. The students will build a mock energy function that models the dihedral rotation of 1,2-dichloroethane. Then three methods (The Monte Carlo Method, Newton minimization, and Metropolis simulated annealing) will be employed to solve for the preferred (lowest energy) dihedral angle (a).
A performance plot will also be generated that shows the lowest energy and its root mean squared deviation RMSD from the known minimum structure (a = 180 degrees). This plot clearly shows Newton’s propensity to get stuck in local minima. It also clearly shows that the Monte Carlo method will always find the global minimum, but with increasing inefficiency. And finally, the Metropolis simulated annealing technique is found to be flexible enough to accurately locate the minimum energy structure every time provided that the step size and temperatures are “tuned”.
Stay “tuned” for a planned video of the spreadsheet in action.
You can participate! Download the Rosetta@home screen saver, and solve protein folding problems in your sleep. (I have no official connection to the Rosetta at home folks, but their work is great! http://boinc.bakerlab.org/rosetta/)