| 9:38:19 AM PST - Mon, Nov 16th 2015 |
|
The qmd module does not calculate a Hessian, so it would not be possible to output the Hessian at every step without modifying the code. I believe a better alternative is to run the qmd calculation and then use the trajectory that is output (the xyz file) to calculate the Hessian at the desired times. This would be most conveniently accomplished by making use of the geometry block's ability to read an xyz file to input coordinates:
geometry
load frame n trajectory.xyz
end
where n would be which coordinates in the xyz file you would want (the first set in the file being 1, the second set 2, etc.), and trajectory.xyz being your xyz file from the qmd simulation.
Best,
Sean
|
|
|