Ab initio molecular dynamics: Theory and Implementation
Ab initio molecular dynamics: Theory and Implementation
Ab initio molecular dynamics: Theory and Implementation
- No tags were found...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
up e.g. in Refs. 132,37,596,597,428,429,333 .A final observation concerning parallel supercomputers might be useful, see alsoSect. 3.9. It is evident from the Lagrangian Eq. (319) <strong>and</strong> the resulting equationsof motion (e.g. Eqs. (331)–(337)) that most of the numerical workload comes fromcalculating the ab <strong>initio</strong> forces on the nuclei. Given a fixed path configurationEq. (316) the P underlying electronic structure problems are independent fromeach other <strong>and</strong> can be solved without communication on P nodes of a distributedmemory machine. Communication is only necessary to send the final result, essentiallythe forces, to a special node that computes the quantum kinetic contributionto the energy <strong>and</strong> integrates finally the equations of motions. It is even conceivableto distribute this task on different supercomputers, i.e. “meta–computing” iswithin reach for such calculations. Thus, the algorithm is “embarrassingly parallel”provided that the memory per node is sufficient to solve the complete Kohn–Shamproblem at a given time slice. If this is not the case the electronic structure calculationitself has to be parallelized on another hierarchical level as outlined inSect. 3.9.4.4.3 <strong>Ab</strong> Initio Path Centroids: DynamicsInitially the <strong>molecular</strong> <strong>dynamics</strong> approach to path integral simulations was inventedmerely as a trick in order to sample configuration space similar to the Monte Carlomethod. This perception changed recently with the introduction of the so–called“centroid <strong>molecular</strong> <strong>dynamics</strong>” technique 102 , see Refs. 103,104,105,665,505,506,507 forbackground information. In a nutshell it is found that the time evolution of thecenters of mass or centroidsR c I(t) = 1 PP∑s ′ =1R (s′ )I(t) (320)of the closed Feynman paths that represent the quantum nuclei contains quasiclassicalinformation about the true quantum <strong>dynamics</strong>. The centroid <strong>molecular</strong><strong>dynamics</strong> approach can be shown to be exact for harmonic potentials <strong>and</strong> to havethe correct classical limit. The path centroids move in an effective potential whichis generated by all the other modes of the paths at the given temperature. Thiseffective potential thus includes the effects of quantum fluctuations on the (quasiclassical)time evolution of the centroid degrees of freedom. Roughly speakingthe trajectory of the path centroids can be regarded as a classical trajectory of thesystem, which is approximately “renormalized” due to quantum effects.The original centroid <strong>molecular</strong> <strong>dynamics</strong> technique 102,103,104,105,665 relies onthe use of model potentials as the st<strong>and</strong>ard time–independent path integral simulations.This limitation was overcome independently in Refs. 469,411 by combiningab <strong>initio</strong> path integrals with centroid <strong>molecular</strong> <strong>dynamics</strong>. The resulting technique,ab <strong>initio</strong> centroid <strong>molecular</strong> <strong>dynamics</strong> can be considered as a quasiclassical generalizationof st<strong>and</strong>ard ab <strong>initio</strong> <strong>molecular</strong> <strong>dynamics</strong>. At the same time, it preservesthe virtues of the ab <strong>initio</strong> path integral technique 395,399,644,404 to generate exacttime–independent quantum equilibrium averages.113