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.
disposable parameters that can be optimized for a particular calculation. Note thatthe discretization points in continuous space can also be considered to constitute asort of “finite basis set” – despite different statements in the literature – <strong>and</strong> thatthe “infinite basis set limit” is reached as h → 0 for N fixed. A variation on thetheme are Mehrstellen schemes where the discretization of the entire differentialequation <strong>and</strong> not only of the derivative operator is optimized 89 .The first real–space approach devised for ab <strong>initio</strong> <strong>molecular</strong> <strong>dynamics</strong> was basedon the lowest–order finite–difference approximation in conjunction with a equally–spaced cubic mesh in real space 109 . A variety of other implementations of moresophisticated real–space methods followed <strong>and</strong> include e.g. non–uniform meshes,multigrid acceleration, different discretization techniques, <strong>and</strong> finite–element methods686,61,39,130,131,632,633,431,634 . Among the chief advantages of the real–spacemethods is that linear scaling approaches 453,243 can be implemented in a naturalway <strong>and</strong> that the multiple–length scale problem can be coped with by adapting thegrid. However, the extension to such non–uniform meshes induces the (in)famousPulay forces (see Sect. 2.5) if the mesh moves as the nuclei move.3 Basic Techniques: <strong>Implementation</strong> within the CPMD Code3.1 Introduction <strong>and</strong> Basic Def<strong>initio</strong>nsThis section discusses the implementation of the plane wave–pseudopotential <strong>molecular</strong><strong>dynamics</strong> method within the CPMD computer code 142 . It concentrates on thebasics leaving advanced methods to later chapters. In addition all formulas are forthe non-spin polarized case. This allows to show the essential features of a planewave code as well as the reasons for its high performance in detail. The implementationof other versions of the presented algorithms <strong>and</strong> of the more advancedtechniques in Sect. 4 is in most cases very similar.There are many reviews on the pseudopotential plane wave method alone or inconnection with the Car–Parrinello algorithm. Older articles 312,157,487,591 as wellas the book by Singh 578 concentrate on the electronic structure part. Other reviews513,472,223,224 present the plane wave method in connection with the <strong>molecular</strong><strong>dynamics</strong> technique.3.1.1 Unit Cell <strong>and</strong> Plane Wave BasisThe unit cell of a periodically repeated system is defined by the Bravais latticevectors a 1 , a 2 , <strong>and</strong> a 3 . The Bravais vectors can be combined into a three by threematrix h = [a 1 ,a 2 ,a 3 ] 459 . The volume Ω of the cell is calculated as the determinantof hΩ = deth . (104)Further, scaled coordinates s are introduced that are related to r via hr = hs . (105)43