- Editors' Suggestion
- Letter
- Open Access
Atomic forces by quantum Monte Carlo: Application to phonon dispersion calculations
Phys. Rev. B 103, L121110 – Published 15 March, 2021
DOI: https://doi.org/10.1103/PhysRevB.103.L121110
Abstract
We report a successful application of the ab initio quantum Monte Carlo (QMC) framework to a phonon dispersion calculation. A full phonon dispersion of diamond is successfully calculated at the variational Monte Carlo (VMC) level, based on the frozen-phonon technique. The VMC-phonon dispersion is in good agreement with the experimental results, giving renormalized harmonic optical frequencies very close to the experimental values, and improving upon previous density functional theory estimates. Key to success for the QMC approach is the statistical error reduction in the atomic force evaluation. We show that this can be achieved by using well conditioned atomic basis sets and by explicitly removing the basis-set redundancy, which reduces the statistical error of forces by up to two orders of magnitude by combining it with the so-called space-warp transformation algorithm. This leads to affordable and accurate QMC-phonons calculations, which are up to times more efficient than a bare force treatment, and paves the way to new applications, particularly in correlated materials, where phonons have been poorly reproduced so far.
Physics Subject Headings (PhySH)
Article Text
Supplemental Material
References (67)
- R. M. Martin, Electronic Structure: Basic Theory and Practical Methods (Cambridge University Press, Cambridge, 2004).
- D. Sholl and J. A. Steckel, Density Functional Theory: A Practical Introduction (John Wiley & Sons, Hoboken, NJ, 2011).
- S. Baroni, S. de Gironcoli, A. Dal Corso, and P. Giannozzi, Rev. Mod. Phys. 73, 515 (2001).
- A. Togo and I. Tanaka, Scr. Mater. 108, 1 (2015).
- N. Mounet and N. Marzari, Phys. Rev. B 71, 205214 (2005).
- M. Lazzeri, C. Attaccalite, L. Wirtz, and F. Mauri, Phys. Rev. B 78, 081406(R) (2008).
- M. Krisch, D. L. Farber, R. Xu, D. Antonangeli, C. M. Aracne, A. Beraud, T.-C. Chiang, J. Zarestky, D. Y. Kim, E. I. Isaev, R. Ahuja, and B. Johansson, Proc. Natl. Acad. Sci. U.S.A. 108, 9342 (2011).
- D. Reznik, G. Sangiovanni, O. Gunnarsson, and T. Devereaux, Nature 455, E6 (2008).
- A. Floris, S. de Gironcoli, E. K. U. Gross, and M. Cococcioni, Phys. Rev. B 84, 161102(R) (2011).
- I. Leonov, A. I. Poteryaev, V. I. Anisimov, and D. Vollhardt, Phys. Rev. B 85, 020401(R) (2012).
- I. Leonov, V. I. Anisimov, and D. Vollhardt, Phys. Rev. Lett. 112, 146401 (2014).
- W. Foulkes, L. Mitas, R. Needs, and G. Rajagopal, Rev. Mod. Phys. 73, 33 (2001).
- M. Marchi, S. Azadi, and S. Sorella, Phys. Rev. Lett. 107, 086807 (2011).
- M. Casula and S. Sorella, Phys. Rev. B 88, 155125 (2013).
- S. Sorella, K. Seki, O. O. Brovko, T. Shirakawa, S. Miyakoshi, S. Yunoki, and E. Tosatti, Phys. Rev. Lett. 121, 066402 (2018).
- A. Zen, D. Zhelyazov, and L. Guidoni, J. Chem. Theory Comput. 8, 4204 (2012).
- Y. Luo, A. Zen, and S. Sorella, J. Chem. Phys. 141, 194112 (2014).
- Y. Y. F. Liu, B. Andrews, and G. J. Conduit, J. Chem. Phys. 150, 034104 (2019).
- K. Nakano, R. Maezono, and S. Sorella, J. Chem. Theory Comput. 15, 4044 (2019).
- R. Maezono, A. Ma, M. D. Towler, and R. J. Needs, Phys. Rev. Lett. 98, 025701 (2007).
- S. Azadi, C. Cavazzoni, and S. Sorella, Phys. Rev. B 82, 125112 (2010).
- M. S. Liu, L. A. Bursill, S. Prawer, and R. Beserman, Phys. Rev. B 61, 3391 (2000).
- M. Schwoerer-Böhning, A. T. Macrander, and D. A. Arms, Phys. Rev. Lett. 80, 5572 (1998).
- J. Kulda, H. Kainzmaier, D. Strauch, B. Dorner, M. Lorenzen, and M. Krisch, Phys. Rev. B 66, 241202(R) (2002).
- D. Vanderbilt, S. G. Louie, and M. L. Cohen, Phys. Rev. Lett. 53, 1477 (1984).
- W. Parrish, Acta Cryst. 13, 838 (1960).
- J. L. Warren, R. G. Wenzel, and J. L. Yarnell, in Inelastic Scattering of Neutrons, Vol. 1 (International Atomic Energy Agency, Vienna, 1965), pp. 361–371; J. L. Warren, J. L. Yarnell, G. Dolling, and R. A. Cowley, Phys. Rev. 158, 805 (1967).
- A. Rohatgi, “ Webplotdigitizer: Version 4.3,” https://automeris.io/WebPlotDigitizer (2020).
- M. Casula, C. Filippi, and S. Sorella, Phys. Rev. Lett. 95, 100201 (2005).
- K. Nakano, R. Maezono, and S. Sorella, Phys. Rev. B 101, 155106 (2020).
- M. Casula and S. Sorella, J. Chem. Phys. 119, 6500 (2003).
- K. Nakano, C. Attaccalite, M. Barborini, L. Capriotti, M. Casula, E. Coccia, M. Dagrada, C. Genovese, Y. Luo, G. Mazzola et al., J. Chem. Phys. 152, 204121 (2020).
- M. Burkatzki, C. Filippi, and M. Dolg, J. Chem. Phys. 126, 234105 (2007).
- See Supplemental Material at http://link.aps.org/supplemental/10.1103/PhysRevB.103.L121110 for details of the DFT and QMC calculations, and for the supplemental Figures S-1-S-7 and Tables S-I-S-III, which includes Refs. [47, 65, 66, 67].
- J. P. Perdew and A. Zunger, Phys. Rev. B 23, 5048 (1981).
- S. Sorella, M. Casula, and D. Rocca, J. Chem. Phys. 127, 014105 (2007).
- C. J. Umrigar, J. Toulouse, C. Filippi, S. Sorella, and R. G. Hennig, Phys. Rev. Lett. 98, 110201 (2007).
- R. G. Hennig, A. Wadehra, K. P. Driver, W. D. Parker, C. J. Umrigar, and J. W. Wilkins, Phys. Rev. B 82, 014101 (2010).
- S. Sorella, M. Casula, L. Spanu, and A. Dal Corso, Phys. Rev. B 83, 075119 (2011).
- F. Becca and S. Sorella, Quantum Monte Carlo Approaches for Correlated Systems (Cambridge University Press, Cambridge, 2017).
- P. Giannozzi, S. Baroni, N. Bonini, M. Calandra, R. Car, C. Cavazzoni, D. Ceresoli, G. L. Chiarotti, M. Cococcioni, I. Dabo et al., J. Phys. Condens. Matter. 21, 395502 (2009).
- J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996).
- P. Pavone, K. Karch, O. Schütt, W. Windl, D. Strauch, P. Giannozzi, and S. Baroni, Phys. Rev. B 48, 3156 (1993).
- G. Kresse, J. Furthmüller, and J. Hafner, Europhys. Lett. 32, 729 (1995).
- F. Occelli, P. Loubeyre, and R. LeToullec, Nat. Mater. 2, 151 (2003).
- H. McSkimin and P. Andreatch Jr., J. Appl. Phys. 43, 2944 (1972).
- P. Vinet, J. R. Smith, J. Ferrante, and J. H. Rose, Phys. Rev. B 35, 1945 (1987).
- S. Sorella and L. Capriotti, J. Chem. Phys. 133, 234111 (2010).
- A conventional supercell containing 64 atoms in the simulation cell is large enough for obtaining a phonon dispersion almost consistent with experimental ones, as shown in the Supplemental Material [34]. This was confirmed by comparing the phonon dispersion obtained by the finite displacement method and that obtained by the linear-response method with very dense and grids at the DFT level.
- N. J. Harris, J. Phys. Chem. 99, 14689 (1995).
- F. Mouhat, S. Sorella, R. Vuilleumier, A. M. Saitta, and M. Casula, J. Chem. Theory Comput. 13, 2400 (2017).
- T. Morresi, L. Paulatto, R. Vuilleumier, and M. Casula, arXiv:2103.04094.
- In the Raman mode, two nearest-neighbor carbon atoms are displaced in opposite directions by a distance from their high-symmetry positions [20]. The distortion in the [1 0 0] direction was employed in this study [25]. We calculated VMC forces for two distorted structures ( bohr and bohr) and obtained the Raman frequency by fitting the forces with a linear function. We also calculated VMC and LRDMC energies for the undistorted structure and two distorted structures ( bohr and bohr), then obtained Raman frequencies by fitting the energies with a quadratic function.
- S. Moroni, S. Saccani, and C. Filippi, J. Chem. Theory Comput. 10, 4823 (2014).
- Since smaller displacements were employed in this study (i.e., 0.03 and 0.05 bohr) than the previous one, the error bars in phonon frequencies at point are a bit larger.
- N. D. Drummond, M. D. Towler, and R. J. Needs, Phys. Rev. B 70, 235119 (2004).
- A. V. Krukau, O. A. Vydrov, A. F. Izmaylov, and G. E. Scuseria, J. Chem. Phys. 125, 224106 (2006).
- J. Sun, A. Ruzsinszky, and J. P. Perdew, Phys. Rev. Lett. 115, 036402 (2015).
- P. Hao, Y. Fang, J. Sun, G. I. Csonka, P. H. T. Philipsen, and J. P. Perdew, Phys. Rev. B 85, 014111 (2012).
- M. F. Peintinger, D. V. Oliveira, and T. Bredow, J. Comput. Chem. 34, 451 (2013).
- = the maximum eigenvalue / the minimum eigenvalue of the overlap matrix ().
- C. Attaccalite and S. Sorella, Phys. Rev. Lett. 100, 114501 (2008).
- A. Subedi and L. Boeri, Phys. Rev. B 84, 020508(R) (2011).
- M. Casula, M. Calandra, and F. Mauri, Phys. Rev. B 86, 075445 (2012).
- H. J. Monkhorst and J. D. Pack, Phys. Rev. B 13, 5188 (1976).
- A. Dal Corso, Comput. Mater. Sci. 95, 337 (2014).
- P. Virtanen, R. Gommers, T. E. Oliphant, M. Haberland, T. Reddy, D. Cournapeau, E. Burovski, P. Peterson, W. Weckesser, J. Bright, S. J. van der Walt, M. Brett, J. Wilson, K. J. Millman, N. Mayorov, A. R. J. Nelson, E. Jones, R. Kern, E. Larson, C. J. Carey, İ. Polat, Y. Feng, E. W. Moore, J. VanderPlas, D. Laxalde, J. Perktold, R. Cimrman, I. Henriksen, E. A. Quintero, C. R. Harris, A. M. Archibald, A. H. Ribeiro, F. Pedregosa, P. van Mulbregt, and SciPy 1.0 Contributors, Nat. Methods 17, 261 (2020).