Abstract
A path-integral molecular dynamics technique for strongly interacting atoms using ab initio potentials derived from density functional theory is implemented. This allows the efficient inclusion of nuclear quantum dispersion in ab initio simulations at finite temperatures. We present an application to the quantum cluster H +5 .
References
Car, R., Parrinello, M.: Phys. Rev. Lett.55, 2471 (1985)
Remler, D.K., Madden, P.A.: Mol. Phys.70, 921 (1990); Galli, G., Parrinello, M.: In: Computer Simulations in Materials Science, p. 282. Meyer, M., Pontikis, V. (eds.). Dordrecht: Kluwer 1991; Payne, M.C., Teter, M.P., Allan, D.C., Arias, T.A., Joannopoulos, J.D.: Rev. Mod. Phys.64, 1045 (1992); especially the density functional aspect is reviewed by Jones, R.O., Gunnarson, O.: Rev. Mod. Phys.61, 689 (1989)
Berne, B.J., Thirumalai, D.: Annu. Rev. Phys. Chem.37, 401 (1986); Schmidt, K.E., Ceperley, D.M.: In: Monte Carlo Methods in Condensed Matter Physics. Binder, K. (ed.). Berlin, Heidelberg, New York: Springer 1992
Marx, D., Nielaba, P., Binder, K.: Phys. Rev. Lett.67, 3124 (1991); Phys. Rev. B47, 7788 (1993)
Parrinello, M., Rahman, A.: J. Chem. Phys.80, 860 (1984); De Raedt, B., Sprik, M., Klein, M.L.: J. Chem. Phys.80, 5719 (1984)
Hall, R.W., Berne, B.J.: J. Chem. Phys.81, 3641 (1984)
Martyna, G.J., Klein, M.L., Tuckerman, M.: J. Chem. Phys.97, 2635 (1992); Tuckerman, M.E., Berne, B.J., Martyna, G.J., Klein, M.L.: J. Chem. Phys.99, 2796 (1993)
Blöchl, P.E., Parrinello, M.: Phys. Rev. B45, 9413 (1992); Fois, E.S., Penman, J.I., Madden, P.A.: J. Chem. Phys.98, 6361 (1993); for reversible integration of NH chains within CPMD see Tuckerman, M.E., Parrinello, M.: J. Chem. Phys. (in press)
Yamaguchi, Y., Gaw, J.F., Remington, R.B., Schaefer III, H. F.: J. Phys. Chem.86, 5072 (1987); Špirko, V., Kraemer, W.P.: J. Mol. Spectrosc.159, 521 (1993)
Parameters: box size: 20 a.u., plane wave cutoff: 25 Ry, pseudopotential: exchange gradient-corrected [13] Car type (Giannozzi, P.: priv. comm.). Tests: (i) H2 bond length and harmonic frequency from Morse fit are 1.41(1.401) a.u. and 4430(4401) cm−1, experimental data in brackets; (ii) H +3 bond length is 1.66 a.u. (cf. 1.65000 a.u. from a benchmark study of Röhse, R., Kutzelnigg, W., Jaquet, R., Klopper, W.: subm. to J. Chem. Phys.); (iii) H2 ground-state nuclear density agreed with the one from the Morse wave functions,P: 32. H +5 quantum simulations:M *H : 4M H, time step: 20 a.u., μ: 3000 a.u. All calculations were performed without imposing periocic boundary conditions (Fois, E.S., Hutter, J.: to be published).Analyzed lengths of the quantum and classic runs were 6850 and 20 000 MD steps, respectively
Nuclear exchange can be incorporated by introducing (bosonic or fixed-node fermionic) Monte Carlo permutation moves along the lines of Ceperley, D.M., Pollock, E.L.: Phys. Rev. Lett.56, 351 (1986); Ceperley, D.M.: Phys. Rev. Lett.69, 331 (1992)
Sindzingre, Ph., Ceperley, D.M., Klein, M.L.: Phys. Rev. Lett.67, 1871 (1991)
Becke, A.: Phys. Rev. A38, 3098 (1988); J. Chem. Phys.15 2155 (1992)