A successful phonon calculation within the Quantum Monte Carlo framework
Scientists expand the scope of the quantum Monte Carlo framework by reducing the error in calculation of atomic forces in solids
Numerical calculations predicting electronic properties of materials independent of experimental data, such as the density functional theory, often fail for "unusual" materials. The quantum Monte Carlo (QMC) framework alleviates this issue but incurs heavy computational costs. In a breakthrough study, an international team of scientists reports enhanced QMC computation speeds based on an error reduction approach and successfully perform proof-of-concept calculations for a periodic solid, opening doors to application of QMC in material development research.
Fig. 1. The major scientific concept of the study
Phonon dispersion of diamond calculated at the variational Monte Carlo level by TurboRVB Image courtesy: Kousuke Nakano from JAIST
The focus and ultimate goal of computational research in materials science and condensed matter physics is to solve the Schrödinger equation--the fundamental equation describing how electrons behave inside matter--exactly (without resorting to simplifying approximations). While experiments can certainly provide interesting insights into a material's properties, it is often computations that reveal the underlying physical mechanism. However, computations need not rely on experimental data and can, in fact, be performed independently, an approach known as "ab initio calculations". The density functional theory (DFT) is a popular example of such an approach.
For most material scientists and condensed matter physicists, DFT calculations are the bread and butter of their profession. However, despite being a powerful technique, DFT has had limited success with "strongly correlated materials"--materials with unusual electronic and magnetic properties. These materials, while interesting on their own, also possess technological useful properties, a fact that strongly motivates an ab initio framework suited to describe them.
To that end, a framework known as "ab initio quantum Monte Carlo" (QMC) has shown considerable promise and is expected to be the next generation of electronic structure calculations due to its superiority over DFT. However, even QMC is largely restricted to calculations of energy and atomic forces, limiting its utility in computing useful material properties.
Now, in a breakthrough study published in Physical Review B (Editors' Suggestion), scientists have taken things to the next level based on an approach that allows them to reduce the statistical error in atomic force evaluation by two orders of magnitude and subsequently speeds up the computation by a factor of 104! "The drastic reduction in computational time will greatly expand the range of QMC calculations and enable highly accurate prediction of atomic properties of materials that have been difficult to handle," observes Assistant Professor Kousuke Nakano from Japan Advanced Institute of Science and Technology (JAIST), who, along with his colleagues Prof. Ryo Maezono from JAIST, Prof. Sandro Sorella from International School for advanced Studies (SISSA), Italy, and Dr. Tommaso Morresi and Prof. Michele Casula from Sorbonne Université, France, led this groundbreaking achievement.
The team applied their developed method to calculate the atomic vibrations of diamond, a typical reference material, as a proof-of-concept and showed that the results were consistent with experimental values. To perform these calculations, they used a large computer, Cray-XC40, located at the Research Center for Advanced Computing Infrastructure at Japan Advanced Institute of Science and Technology (JAIST), Japan, along with another located at RIKEN, Japan. The team made use of a QMC software package called "TurboRVB", initially launched by Prof. Sorella and Prof. Casula and developed later by Prof. Nakano along with others, to perform phonon dispersion calculations for diamond that were previously inaccessible, greatly expanding its scope.
Prof. Nakano looks forward to the applications of QMC in materials informatics (MI), a field dedicated to the design and search for novel materials using techniques of information science and computational physics. "While MI is currently governed by DFT, the rapid developments in computer performance, such as the exascale supercomputer, will help QMC gain popularity. In that regard, our developed method is going to be very useful for designing novel materials with real-life applications," concludes an optimistic Dr. Nakano.
|Title of original paper:||Atomic forces by quantum Monte Carlo: application to phonon dispersion calculation|
|Journal:||Physical Review B (Editors' Suggestion)|
Computational resources were provided from the facilities of Research Center for Advanced Computing Infrastructure at Japan Advanced Institute of Science and Technology (JAIST), from PRACE Project No. 2019204934, and from RIKEN (the supercomputer Fugaku) through the HPCI System Re- search Project (Project ID: hp200164). This study was funded by French state funds managed by the ANR within the Investissements d'Avenir programme under reference ANR-11-IDEX-0004-02, and more specifically within the framework of the Cluster of Excellence MATISSE led by Sorbonne University, the Simons Foundation and Grant-in-Aid for Scientific Research on Innovative Areas (No. 16H06439), MEXT-KAKENHI (19H04692 and 16KK0097), Toyota Motor Corporation, Air Force Office of Scientific Research (AFOSR-AOARD/FA2386-17-1-4049; FA2386-19-1-4015), JSPS Bilateral Joint Projects (with India DST), and PRIN 2017BZPKSZ.
March 19, 2021