Showing 1–50 of 97 results for author: Kresse, G

Density-Based Long-Range Electrostatic Descriptors for Machine Learning Force Fields

Authors: Carolin Faller, Merzuk Kaltak, Georg Kresse

Abstract: This study presents a long-range descriptor for machine learning force fields (MLFFs) that maintains translational and rotational symmetry, similar to short-range descriptors while being able to incorporate long-range electrostatic interactions. The proposed descriptor is based on an atomic density representation and is structurally similar to classical short-range atom-centered descriptors, makin… ▽ More This study presents a long-range descriptor for machine learning force fields (MLFFs) that maintains translational and rotational symmetry, similar to short-range descriptors while being able to incorporate long-range electrostatic interactions. The proposed descriptor is based on an atomic density representation and is structurally similar to classical short-range atom-centered descriptors, making it straightforward to integrate into machine learning schemes. The effectiveness of our model is demonstrated through comparative analysis with the long-distance equivariant (LODE) descriptor. In a toy model with purely electrostatic interactions, our model achieves errors below 0.1%. The application of our descriptors, in combination with local descriptors representing the atomic density, to materials where monopole-monopole interactions are important such as sodium chloride successfully captures long-range interactions, improving predictive accuracy. The study highlights the limitations of the combined LODE method in materials where intermediate-range effects play a significant role. Our work presents a promising approach to addressing the challenge of incorporating long-range interactions into MLFFs, which enhances predictive accuracy for charged materials to the level of state-of-the-art Message Passing Neural Networks. △ Less

Submitted 25 June, 2024; originally announced June 2024.

Stoichiometric reconstruction of the Al$_{2}$O$_{3}$(0001) surface

Authors: Johanna I. Hütner, Andrea Conti, David Kugler, Florian Mittendorfer, Georg Kresse, Michael Schmid, Ulrike Diebold, Jan Balajka

Abstract: Macroscopic properties of materials stem from fundamental atomic-scale details, yet for insulators, resolving surface structures remains a challenge. The basal (0001) plane of $α$-Al$_{2}$O$_{3}$ was imaged with noncontact atomic force microscopy with an atomically-defined tip apex. The surface forms a complex $({\sqrt31} {\times} {\sqrt31})R{\pm}9°$ reconstruction. The lateral positions of the in… ▽ More Macroscopic properties of materials stem from fundamental atomic-scale details, yet for insulators, resolving surface structures remains a challenge. The basal (0001) plane of $α$-Al$_{2}$O$_{3}$ was imaged with noncontact atomic force microscopy with an atomically-defined tip apex. The surface forms a complex $({\sqrt31} {\times} {\sqrt31})R{\pm}9°$ reconstruction. The lateral positions of the individual O and Al surface atoms come directly from experiment; how these connect to the underlying crystal bulk was determined based on computational modeling. Before the restructuring, the surface Al atoms assume an unfavorable, threefold planar coordination; the reconstruction allows a rehybridization with subsurface O that leads to a substantial energy gain. The reconstructed surface remains stoichiometric, Al$_{2}$O$_{3}$. △ Less

Submitted 29 May, 2024; originally announced May 2024.

Derivative learning of tensorial quantities -- Predicting finite temperature infrared spectra from first principles

Authors: Bernhard Schmiedmayer, Georg Kresse

Abstract: We develop a strategy that integrates machine learning and first-principles calculations to achieve technical accurate predictions of infrared spectra. Specifically, the methodology allows to predict infrared spectra for complex systems at finite temperatures. The method's effectiveness is demonstrated in challenging scenarios, such as the analysis of water and the organic-inorganic halide perovsk… ▽ More We develop a strategy that integrates machine learning and first-principles calculations to achieve technical accurate predictions of infrared spectra. Specifically, the methodology allows to predict infrared spectra for complex systems at finite temperatures. The method's effectiveness is demonstrated in challenging scenarios, such as the analysis of water and the organic-inorganic halide perovskite MAPbI$_{3}$, where our results consistently align with experimental data. A distinctive feature of the methodology is the incorporation of derivative learning, which proves indispensable for obtaining accurate polarization data in bulk materials and facilitates the training of a machine learning surrogate model of the polarization adapted to rotational and translational symmetries. We achieve polarisation prediction accuracies of about 1 % by training only on the predicted Born effective charges. △ Less

Submitted 30 April, 2024; originally announced April 2024.

Comments: 9 pages, 4 figures

Towards Large-Scale AFQMC Calculations: Large Time Step Auxiliary-Field Quantum Monte Carlo

Authors: Zoran Sukurma, Martin Schlipf, Moritz Humer, Amir Taheridehkordi, Georg Kresse

Abstract: We report modifications of the ph-AFQMC algorithm that allow the use of large time steps and reliable time step extrapolation. Our modified algorithm eliminates size-consistency errors present in the standard algorithm when large time steps are employed. We investigate various methods to approximate the exponential of the one-body operator within the AFQMC framework, distinctly demonstrating the s… ▽ More We report modifications of the ph-AFQMC algorithm that allow the use of large time steps and reliable time step extrapolation. Our modified algorithm eliminates size-consistency errors present in the standard algorithm when large time steps are employed. We investigate various methods to approximate the exponential of the one-body operator within the AFQMC framework, distinctly demonstrating the superiority of Krylov methods over the conventional Taylor expansion. We assess various propagators within AFQMC and demonstrate that the Split-2 propagator is the optimal method, exhibiting the smallest time-step errors. For the HEAT set molecules, the time-step extrapolated energies deviate on average by only 0.19 kcal/mol from the accurate small time-step energies. For small water clusters, we obtain accurate complete basis-set binding energies using time-step extrapolation with a mean absolute error of 0.07 kcal/mol compared to CCSD(T). Using large time-step ph-AFQMC for the N$_2$ dimer, we show that accurate bond lengths can be obtained while reducing CPU time by an order of magnitude. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: 13 pages, 6 figures

Comparing machine learning potentials for water: Kernel-based regression and Behler-Parrinello neural networks

Authors: Pablo Montero de Hijes, Christoph Dellago, Ryosuke Jinnouchi, Bernhard Schmiedmayer, Georg Kresse

Abstract: In this paper we investigate the performance of different machine learning potentials (MLPs) in predicting key thermodynamic properties of water using RPBE+D3. Specifically, we scrutinize kernel-based regression and high-dimensional neural networks trained on a highly accurate dataset consisting of about 1,500 structures, as well as a smaller data set, about half the size, obtained using only on-t… ▽ More In this paper we investigate the performance of different machine learning potentials (MLPs) in predicting key thermodynamic properties of water using RPBE+D3. Specifically, we scrutinize kernel-based regression and high-dimensional neural networks trained on a highly accurate dataset consisting of about 1,500 structures, as well as a smaller data set, about half the size, obtained using only on-the-fly learning. The study reveals that despite minor differences between the MLPs, their agreement on observables such as the diffusion constant and pair-correlation functions is excellent, especially for the large training dataset. Variations in the predicted density isobars, albeit somewhat larger, are also acceptable, particularly given the errors inherent to approximate density functional theory. Overall, the study emphasizes the relevance of the database over the fitting method. Finally, the study underscores the limitations of root mean square errors and the need for comprehensive testing, advocating the use of multiple MLPs for enhanced certainty, particularly when simulating complex thermodynamic properties that may not be fully captured by simpler tests. △ Less

Submitted 23 December, 2023; originally announced December 2023.

Layer-by-layer phase transformation in Ti$_3$O$_5$ revealed by machine learning molecular dynamics simulations

Authors: Mingfeng Liu, Jiantao Wang, Junwei Hu, Peitao Liu, Haiyang Niu, Xuexi Yan, Jiangxu Li, Haile Yan, Bo Yang, Yan Sun, Chunlin Chen, Georg Kresse, Liang Zuo, Xing-Qiu Chen

Abstract: Reconstructive phase transitions involving breaking and reconstruction of primary chemical bonds are ubiquitous and important for many technological applications. In contrast to displacive phase transitions, the dynamics of reconstructive phase transitions are usually slow due to the large energy barrier. Nevertheless, the reconstructive phase transformation from $β$- to $λ$-Ti$_3$O$_5$ exhibits a… ▽ More Reconstructive phase transitions involving breaking and reconstruction of primary chemical bonds are ubiquitous and important for many technological applications. In contrast to displacive phase transitions, the dynamics of reconstructive phase transitions are usually slow due to the large energy barrier. Nevertheless, the reconstructive phase transformation from $β$- to $λ$-Ti$_3$O$_5$ exhibits an ultrafast and reversible behavior. Despite extensive studies, the underlying microscopic mechanism remains unclear. Here, we discover a kinetically favorable in-plane nucleated layer-by-layer transformation mechanism through metadynamics and large-scale molecular dynamics simulations. This is enabled by developing an efficient machine learning potential with near first-principles accuracy through an on-the-fly active learning method and an advanced sampling technique. Our results reveal that the $β$-$λ$ phase transformation initiates with the formation of two-dimensional nuclei in the $ab$-plane and then proceeds layer-by-layer through a multistep barrier-lowering kinetic process via intermediate metastable phases. Our work not only provides important insight into the ultrafast and reversible nature of the $β$-$λ$ transition, but also presents useful strategies and methods for tackling other complex structural phase transitions. △ Less

Submitted 2 April, 2024; v1 submitted 9 October, 2023; originally announced October 2023.

Comments: 26 pages,23 figures (including Supporting Information)

Journal ref: Nat Commun 15, 3079 (2024)

Machine Learning-Aided First-Principles Calculations of Redox Potentials

Authors: Ryosuke Jinnouchi, Ferenc Karsai, Georg Kresse

Abstract: Redox potentials of electron transfer reactions are of fundamental importance for the performance and description of electrochemical devices. Despite decades of research, accurate computational predictions for the redox potential of even simple metals remain very challenging. Here we use a combination of first principles calculations and machine learning to predict the redox potentials of three re… ▽ More Redox potentials of electron transfer reactions are of fundamental importance for the performance and description of electrochemical devices. Despite decades of research, accurate computational predictions for the redox potential of even simple metals remain very challenging. Here we use a combination of first principles calculations and machine learning to predict the redox potentials of three redox couples, $\mathrm{Fe}^{2+}$/$\mathrm{Fe}^{3+}$, $\mathrm{Cu}^{+}$/$\mathrm{Cu}^{2+}$ and $\mathrm{Ag}^{+}$/$\mathrm{Ag}^{2+}$. Using a hybrid functional with a fraction of 25\% exact exchange (PBE0) the predicted values are 0.92, 0.26 and 1.99 V in good agreement with the best experimental estimates (0.77, 0.15, 1.98 V). We explain in detail, how we combine machine learning, thermodynamic integration from machine learning to semi-local functionals, as well as a combination of thermodynamic perturbation theory and $Δ$-machine learning to determine the redox potentials for computationally expensive hybrid functionals. The combination of these approaches allows one to obtain statistically accurate results. △ Less

Submitted 22 March, 2024; v1 submitted 22 September, 2023; originally announced September 2023.

Machine learning density functionals from the random-phase approximation

Authors: Stefan Riemelmoser, Carla Verdi, Merzuk Kaltak, Georg Kresse

Abstract: Kohn-Sham density functional theory (DFT) is the standard method for first-principles calculations in computational chemistry and materials science. More accurate theories such as the random-phase approximation (RPA) are limited in application due to their large computational cost. Here, we construct a DFT substitute functional for the RPA using supervised and unsupervised machine learning (ML) te… ▽ More Kohn-Sham density functional theory (DFT) is the standard method for first-principles calculations in computational chemistry and materials science. More accurate theories such as the random-phase approximation (RPA) are limited in application due to their large computational cost. Here, we construct a DFT substitute functional for the RPA using supervised and unsupervised machine learning (ML) techniques. Our ML-RPA model can be interpreted as a non-local extension to the standard gradient approximation. We train an ML-RPA functional for diamond surfaces and liquid water and show that ML-RPA can outperform the standard gradient functionals in terms of accuracy. Our work demonstrates how ML-RPA can extend the applicability of the RPA to larger system sizes, time scales and chemical spaces. △ Less

Submitted 1 August, 2023; originally announced August 2023.

Journal ref: Journal of Chemical Theory and Computation 2023 19 (20), 7287-7299

How to verify the precision of density-functional-theory implementations via reproducible and universal workflows

Authors: Emanuele Bosoni, Louis Beal, Marnik Bercx, Peter Blaha, Stefan Blügel, Jens Bröder, Martin Callsen, Stefaan Cottenier, Augustin Degomme, Vladimir Dikan, Kristjan Eimre, Espen Flage-Larsen, Marco Fornari, Alberto Garcia, Luigi Genovese, Matteo Giantomassi, Sebastiaan P. Huber, Henning Janssen, Georg Kastlunger, Matthias Krack, Georg Kresse, Thomas D. Kühne, Kurt Lejaeghere, Georg K. H. Madsen, Martijn Marsman , et al. (20 additional authors not shown)

Abstract: In the past decades many density-functional theory methods and codes adopting periodic boundary conditions have been developed and are now extensively used in condensed matter physics and materials science research. Only in 2016, however, their precision (i.e., to which extent properties computed with different codes agree among each other) was systematically assessed on elemental crystals: a firs… ▽ More In the past decades many density-functional theory methods and codes adopting periodic boundary conditions have been developed and are now extensively used in condensed matter physics and materials science research. Only in 2016, however, their precision (i.e., to which extent properties computed with different codes agree among each other) was systematically assessed on elemental crystals: a first crucial step to evaluate the reliability of such computations. We discuss here general recommendations for verification studies aiming at further testing precision and transferability of density-functional-theory computational approaches and codes. We illustrate such recommendations using a greatly expanded protocol covering the whole periodic table from Z=1 to 96 and characterizing 10 prototypical cubic compounds for each element: 4 unaries and 6 oxides, spanning a wide range of coordination numbers and oxidation states. The primary outcome is a reference dataset of 960 equations of state cross-checked between two all-electron codes, then used to verify and improve nine pseudopotential-based approaches. Such effort is facilitated by deploying AiiDA common workflows that perform automatic input parameter selection, provide identical input/output interfaces across codes, and ensure full reproducibility. Finally, we discuss the extent to which the current results for total energies can be reused for different goals (e.g., obtaining formation energies). △ Less

Submitted 26 May, 2023; originally announced May 2023.

Comments: Main text: 23 pages, 4 figures. Supplementary: 68 pages. Nature Review Physics 2023

Journal ref: Nat. Rev. Phys. 6, 45 (2024)

Phaseless auxiliary field quantum Monte Carlo with projector-augmented wave method for solids

Authors: Amir Taheridehkordi, Martin Schlipf, Zoran Sukurma, Moritz Humer, Andreas Grüneis, Georg Kresse

Abstract: We implement the phaseless auxiliary field quantum Monte Carlo method using the plane-wave based projector augmented wave method and explore the accuracy and the feasibility of applying our implementation to solids. We use a singular value decomposition to compress the two-body Hamiltonian and thus reduce the computational cost. Consistent correlation energies from the primitive-cell sampling and… ▽ More We implement the phaseless auxiliary field quantum Monte Carlo method using the plane-wave based projector augmented wave method and explore the accuracy and the feasibility of applying our implementation to solids. We use a singular value decomposition to compress the two-body Hamiltonian and thus reduce the computational cost. Consistent correlation energies from the primitive-cell sampling and the corresponding supercell calculations numerically verify our implementation. We calculate the equation of state for diamond and the correlation energies for a range of prototypical solid materials. A down-sampling technique along with natural orbitals accelerates the convergence with respect to the number of orbitals and crystal momentum points. We illustrate the competitiveness of our implementation in accuracy and computational cost for dense crystal momentum point meshes comparing to a well-established quantum-chemistry approach, the coupled-cluster ansatz including singles, doubles and perturbative triple particle-hole excitation operators. △ Less

Submitted 27 April, 2023; originally announced April 2023.

Comments: 13 pages, 7 figures

Benchmark Phaseless Auxiliary-Field Quantum Monte Carlo Method for Small Molecules

Authors: Z. Sukurma, M. Schlipf, M. Humer, A. Taheridehkordi, G. Kresse

Abstract: We report a scalable Fortran implementation of the phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) and demonstrate its excellent performance and beneficial scaling with respect to system size. Furthermore, we investigate modifications of the phaseless approximation that can help to reduce the overcorrelation problems common to the ph-AFQMC. We apply the method to the 26 molecules in the H… ▽ More We report a scalable Fortran implementation of the phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) and demonstrate its excellent performance and beneficial scaling with respect to system size. Furthermore, we investigate modifications of the phaseless approximation that can help to reduce the overcorrelation problems common to the ph-AFQMC. We apply the method to the 26 molecules in the HEAT set, the benzene molecule, and water clusters. We observe a mean absolute deviation of the total energy of 1.15 kcal/mol for the molecules in the HEAT set; close to chemical accuracy. For the benzene molecule, the modified algorithm despite using a single-Slater-determinant trial wavefunction yields the same accuracy as the original phaseless scheme with 400 Slater determinants. Despite these improvements, we find systematic errors for the CN, CO$_2$, and O$_2$ molecules that need to be addressed with more accurate trial wavefunctions. For water clusters, we find that the ph-AFQMC yields excellent binding energies that differ from CCSD(T) by typically less than 0.5 kcal/mol. △ Less

Submitted 7 March, 2023; originally announced March 2023.

Comments: 13 pages, 7 figures

Journal ref: J. Chem. Theory Comput. 2023, 19, 15, 4921-4932

Quantum paraelectricity and structural phase transitions in strontium titanate beyond density-functional theory

Authors: Carla Verdi, Luigi Ranalli, Cesare Franchini, Georg Kresse

Abstract: We demonstrate an approach for calculating temperature-dependent quantum and anharmonic effects with beyond density-functional theory accuracy. By combining machine-learned potentials and the stochastic self-consistent harmonic approximation, we investigate the cubic to tetragonal transition in strontium titanate and show that the paraelectric phase is stabilized by anharmonic quantum fluctuations… ▽ More We demonstrate an approach for calculating temperature-dependent quantum and anharmonic effects with beyond density-functional theory accuracy. By combining machine-learned potentials and the stochastic self-consistent harmonic approximation, we investigate the cubic to tetragonal transition in strontium titanate and show that the paraelectric phase is stabilized by anharmonic quantum fluctuations. We find that a quantitative understanding of the quantum paraelectric behavior requires a higher-level treatment of electronic correlation effects via the random phase approximation. This approach enables detailed studies of emergent properties in strongly anharmonic materials beyond density-functional theory. △ Less

Submitted 1 March, 2023; v1 submitted 17 November, 2022; originally announced November 2022.

Temperature-dependent anharmonic phonons in quantum paraelectric KTaO$_3$ by first principles and machine-learned force fields

Authors: Luigi Ranalli, Carla Verdi, Lorenzo Monacelli, Matteo Calandra, Georg Kresse, Cesare Franchini

Abstract: Understanding collective phenomena in quantum materials from first principles is a promising route toward engineering materials properties on demand and designing new functionalities. This work examines the quantum paraelectric state, an elusive state of matter characterized by the smooth saturation of the ferroelectric instability at low temperature due to quantum fluctuations associated with anh… ▽ More Understanding collective phenomena in quantum materials from first principles is a promising route toward engineering materials properties on demand and designing new functionalities. This work examines the quantum paraelectric state, an elusive state of matter characterized by the smooth saturation of the ferroelectric instability at low temperature due to quantum fluctuations associated with anharmonic phonon effects. The temperature-dependent evolution of the soft ferroelectric phonon mode in the quantum paraelectric KTaO$_3$ in the range 0-300 K is modelled by combining density functional theory (DFT) calculations with the stochastic self-consistent harmonic approximation assisted by an on-the-fly machine-learned force field. The calculated data show that including anharmonic terms is essential to stabilize the spurious imaginary ferroelectric phonon predicted by DFT, in agreement with experiments. Augmenting the DFT workflow with machine-learned force fields allows for efficient stochastic sampling of the configurational space using large supercells in a broad and dense temperature range, inaccessible by conventional ab initio protocols. This work proposes a robust computational workflow capable of accounting for collective behaviors involving different degrees of freedom and occurring at large time/length scales, paving the way for precise modeling and control of quantum effects in materials. △ Less

Submitted 24 September, 2022; originally announced September 2022.

Comments: 10 pages, 6 figures. To be published in Advanced Quantum Technologies. Editor: Luigi Ranalli

Approaching the basis-set limit of the dRPA correlation energy with explicitly correlated and Projector Augmented-wave methods

Authors: Moritz Humer, Michael E. Harding, Martin Schlipf, Amir Taheridehkordi, Zoran Sukurma, Wim Klopper, Georg Kresse

Abstract: The direct random-phase approximation (dRPA) is used to calculate and compare atomization energies for the HEAT set and 10 selected molecules of the G2-1 set using both plane waves and Gaussian-type orbitals. We describe detailed procedures to obtain highly accurate and well converged results for the projector augmented-wave (PAW) method as implemented in the Vienna Ab-initio Simulation Package (V… ▽ More The direct random-phase approximation (dRPA) is used to calculate and compare atomization energies for the HEAT set and 10 selected molecules of the G2-1 set using both plane waves and Gaussian-type orbitals. We describe detailed procedures to obtain highly accurate and well converged results for the projector augmented-wave (PAW) method as implemented in the Vienna Ab-initio Simulation Package (VASP) as well as the explicitly correlated dRPA-F12 method as implemented in the TURBOMOLE package. The two approaches agree within chemical accuracy (1 kcal/mol) for the atomization energies of all considered molecules, both for the exact exchange as well as for the dRPA. The root mean-square deviation is 0.41 kcal/mol for the exact exchange (evaluated using density functional theory orbitals) and 0.33 kcal/mol for exact exchange plus the random-phase approximation. △ Less

Submitted 31 August, 2022; originally announced August 2022.

Comments: 16 pages, 2 figures

Combining Machine Learning and Many-Body Calculations: Coverage-Dependent Adsorption of CO on Rh(111)

Authors: Peitao Liu, Jiantao Wang, Noah Avargues, Carla Verdi, Andreas Singraber, Ferenc Karsai, Xing-Qiu Chen, Georg Kresse

Abstract: Adsorption of carbon monoxide (CO) on transition-metal surfaces is a prototypical process in surface sciences and catalysis. Despite its simplicity, it has posed great challenges to theoretical modeling. Pretty much all existing density functionals fail to accurately describe surface energies, CO adsorption site preference, as well as adsorption energies simultaneously. Although the random phase a… ▽ More Adsorption of carbon monoxide (CO) on transition-metal surfaces is a prototypical process in surface sciences and catalysis. Despite its simplicity, it has posed great challenges to theoretical modeling. Pretty much all existing density functionals fail to accurately describe surface energies, CO adsorption site preference, as well as adsorption energies simultaneously. Although the random phase approximation (RPA) cures these density functional theory failures, its large computational cost makes it prohibitive to study the CO adsorption for any but the simplest ordered cases. Here, we address these challenges by developing a machine-learned force field (MLFF) with near RPA accuracy for the prediction of coverage-dependent adsorption of CO on the Rh(111) surface through an efficient on-the-fly active learning procedure and a $Δ$-machine learning approach. We show that the RPA-derived MLFF is capable to accurately predict the Rh(111) surface energy, CO adsorption site preference as well as adsorption energies at different coverages that are all in good agreement with experiments. Moreover, the coverage-dependent ground-state adsorption patterns and adsorption saturation coverage are identified. △ Less

Submitted 17 February, 2023; v1 submitted 20 August, 2022; originally announced August 2022.

Comments: 6 pages, 2 figures

Journal ref: Phys. Rev. Lett. 130, 078001 (2023)

Core-hole excitations using the projector augmented-wave method and the Bethe-Salpeter equation

Authors: Martin Unzog, Alexey Tal, Georg Kresse

Abstract: We present an implementation of the Bethe-Salpeter equation (BSE) for core-conduction band pairs within the framework of the projector augmented-wave method. For validation, the method is applied to the $K$-edges of diamond, graphite, hexagonal boron-nitride, as well as four lithium-halides (LiF, LiCl, LiI, LiBr). We compare our results with experiment, previous theoretical BSE results, and the de… ▽ More We present an implementation of the Bethe-Salpeter equation (BSE) for core-conduction band pairs within the framework of the projector augmented-wave method. For validation, the method is applied to the $K$-edges of diamond, graphite, hexagonal boron-nitride, as well as four lithium-halides (LiF, LiCl, LiI, LiBr). We compare our results with experiment, previous theoretical BSE results, and the density functional theory-based supercell core-hole method. In all considered cases, the agreement with experiment is excellent, in particular for the position of the peaks as well as the fine structure. Comparing BSE to supercell core-hole spectra we find that the latter often qualitatively reproduces the experimental spectrum, however, it sometimes lacks important details. This is shown for the $K$-edges of diamond and nitrogen in hexagonal boron-nitride, where we are capable to resolve within the BSE experimental features that are lacking in the core-hole method. Additionally, we show that in certain systems the supercell core-hole method performs better if the excited electron is added to the background charge. We attribute this improved performance to a reduced self-interaction. △ Less

Submitted 23 June, 2022; originally announced June 2022.

Zero-point Renormalization of the Band Gap of Semiconductors and Insulators Using the PAW Method

Authors: Manuel Engel, Henrique Miranda, Laurent Chaput, Atsushi Togo, Carla Verdi, Martijn Marsman, Georg Kresse

Abstract: We evaluate the zero-point renormalization (ZPR) due to electron-phonon interactions of 28 solids using the projector-augmented-wave (PAW) method. The calculations cover diamond, many zincblende semiconductors, rock-salt and wurtzite oxides, as well as silicate and titania. Particular care is taken to include long-range electrostatic interactions via a generalized Fröhlich model, as discussed in P… ▽ More We evaluate the zero-point renormalization (ZPR) due to electron-phonon interactions of 28 solids using the projector-augmented-wave (PAW) method. The calculations cover diamond, many zincblende semiconductors, rock-salt and wurtzite oxides, as well as silicate and titania. Particular care is taken to include long-range electrostatic interactions via a generalized Fröhlich model, as discussed in Phys. Rev. Lett. 115, 176401 (2015) and Phys. Rev. B 92, 054307 (2015). The data are compared to recent calculations, npj Computational Materials 6, 167 (2020), and generally very good agreement is found. We discuss in detail the evaluation of the electron-phonon matrix elements within the PAW method. We show that two distinct versions can be obtained depending on when the atomic derivatives are taken. If the PAW transformation is applied before taking derivatives with respect to the ionic positions, equations similar to the ones conventionally used in pseudopotential codes are obtained. If the PAW transformation is used after taking the derivatives, the full-potential spirit is largely maintained. We show that both variants yield very similar ZPRs for selected materials when the rigid-ion approximation is employed. In practice, we find however that the pseudo version converges more rapidly with respect to the number of included unoccupied states. △ Less

Submitted 28 June, 2022; v1 submitted 9 May, 2022; originally announced May 2022.

Comments: Added pseudopotential information; revised Sec. II B to avoid possible misconceptions; revised discussion in Sec. III A; added diamond calculation with harder pseudopotential; corrected ZPR comparisons for AlSb and GaP due to inconsistent band-gap transitions; conclusion unchanged

Thermal transport and phase transitions of zirconia by on-the-fly machine-learned interatomic potentials

Authors: Carla Verdi, Ferenc Karsai, Peitao Liu, Ryosuke Jinnouchi, Georg Kresse

Abstract: Machine-learned interatomic potentials enable realistic finite temperature calculations of complex materials properties with first-principles accuracy. It is not yet clear, however, how accurately they describe anharmonic properties, which are crucial for predicting the lattice thermal conductivity and phase transitions in solids and, thus, shape their technological applications. Here we employ a… ▽ More Machine-learned interatomic potentials enable realistic finite temperature calculations of complex materials properties with first-principles accuracy. It is not yet clear, however, how accurately they describe anharmonic properties, which are crucial for predicting the lattice thermal conductivity and phase transitions in solids and, thus, shape their technological applications. Here we employ a recently developed on-the-fly learning technique based on molecular dynamics and Bayesian inference in order to generate an interatomic potential capable to describe the thermodynamic properties of zirconia, an important transition metal oxide. This machine-learned potential accurately captures the temperature-induced phase transitions below the melting point. We further showcase the predictive power of the potential by calculating the heat transport on the basis of Green-Kubo theory, which allows to account for anharmonic effects to all orders. This study indicates that machine-learned potentials trained on the fly offer a routine solution for accurate and efficient simulations of the thermodynamic properties of a vast class of anharmonic materials. △ Less

Submitted 30 September, 2021; originally announced September 2021.

Journal ref: npj Comput. Mater. 7, 156 (2021)

Phase Transitions of Zirconia: Machine-Learned Force Fields Beyond Density Functional Theory

Authors: Peitao Liu, Carla Verdi, Ferenc Karsai, Georg Kresse

Abstract: We present an approach to generate machine-learned force fields (MLFF) with beyond density functional theory (DFT) accuracy. Our approach combines on-the-fly active learning and $Δ$-machine learning in order to generate an MLFF for zirconia based on the random phase approximation (RPA). Specifically, an MLFF trained on-the-fly during DFT based molecular dynamics simulations is corrected by another… ▽ More We present an approach to generate machine-learned force fields (MLFF) with beyond density functional theory (DFT) accuracy. Our approach combines on-the-fly active learning and $Δ$-machine learning in order to generate an MLFF for zirconia based on the random phase approximation (RPA). Specifically, an MLFF trained on-the-fly during DFT based molecular dynamics simulations is corrected by another MLFF that is trained on the differences between RPA and DFT calculated energies, forces and stress tensors. Thanks to the relatively smooth nature of the differences, the expensive RPA calculations are performed only on a small number of representative structures of small unit cells. These structures are determined by a singular value decomposition rank compression of the kernel matrix with low spatial resolution. This dramatically reduces the computational cost and allows us to generate an MLFF fully capable of reproducing high-level quantum-mechanical calculations beyond DFT. We carefully validate our approach and demonstrate its success in studying the phase transitions of zirconia. △ Less

Submitted 13 September, 2021; originally announced September 2021.

Optical and excitonic properties of transition metal oxide perovskites by the Bethe-Salpeter equation

Authors: Lorenzo Varrassi, Peitao Liu, Zeynep Ergönenc Yavas, Menno Bokdam, Georg Kresse, Cesare Franchini

Abstract: We present a systematic investigation of the role and importance of excitonic effects on the optical properties of transitions metal oxide perovskites. A representative set of fourteen compounds has been selected, including 3$d$ (SrTiO$_3$, LaScO$_3$, LaTiO$_3$, LaVO$_3$, LaCrO$_3$, LaMnO$_3$, LaFeO$_3$ and SrMnO$_3$), 4$d$ (SrZrO$_3$, SrTcO$_3$ and Ca$_2$RuO$_4$) and 5$d$ (SrHfO$_3$, KTaO$_3$ and… ▽ More We present a systematic investigation of the role and importance of excitonic effects on the optical properties of transitions metal oxide perovskites. A representative set of fourteen compounds has been selected, including 3$d$ (SrTiO$_3$, LaScO$_3$, LaTiO$_3$, LaVO$_3$, LaCrO$_3$, LaMnO$_3$, LaFeO$_3$ and SrMnO$_3$), 4$d$ (SrZrO$_3$, SrTcO$_3$ and Ca$_2$RuO$_4$) and 5$d$ (SrHfO$_3$, KTaO$_3$ and NaOsO$_3$) perovskites, covering a band gap ranging from 0.1 eV to 6.1 eV and exhibiting different electronic, structural and magnetic properties. Optical conductivities and optical transitions including electron-hole interactions are calculated through the solution of the Bethe-Salpeter equation (BSE) with quasi-particle energies evaluated by single-shot $G_0W_0$ approximation. The exciton binding energies are computed by means of a model-BSE (mBSE), carefully benchmarked against the full BSE method, in order to obtain well-converged results in terms of k-point sampling. The predicted results are compared with available measured data, with an overall satisfactory agreement between theory and experiment. △ Less

Submitted 28 June, 2021; originally announced June 2021.

Electronic state unfolding for plane waves: energy bands, Fermi surfaces and spectral functions

Authors: David Dirnberger, Georg Kresse, Cesare Franchini, Michele Reticcioli

Abstract: Modern computing facilities grant access to first-principles density-functional theory study of complex physical and chemical phenomena in materials, that require large supercell to properly model the system. However, supercells are associated to small Brillouin zones in the reciprocal space, leading to folded electronic eigenstates that make the analysis and interpretation extremely challenging.… ▽ More Modern computing facilities grant access to first-principles density-functional theory study of complex physical and chemical phenomena in materials, that require large supercell to properly model the system. However, supercells are associated to small Brillouin zones in the reciprocal space, leading to folded electronic eigenstates that make the analysis and interpretation extremely challenging. Various techniques have been proposed and developed in order to reconstruct the electronic band structures of super cells, unfolded into the reciprocal space of an ideal primitive cell. Here, we propose an efficient unfolding scheme embedded directly in the Vienna Ab-initio Simulation Package (VASP), that requires modest computational resources and allows for an automatized mapping from the reciprocal space of the supercell to primitive cell Brillouin zone. This algorithm can computes band structures, Fermi surfaces and spectral functions, by using an integrated post-processing tool (bands4vasp). The method is here applied to a selected variety of complex physical situations: the effect of doping on the band dispersion in the BaFe$_{\rm 2(1-x)}$Ru$_{\rm 2x}$As$_2$ superconductor, the interaction between adsorbates and polaronic states on the TiO$_2$(110) surface, and the band splitting induced by non-collinear spin fluctuations in EuCd$_2$As$_2$. △ Less

Submitted 26 May, 2021; v1 submitted 17 March, 2021; originally announced March 2021.

Optimized effective potentials from the random-phase approximation: Accuracy of the quasiparticle approximation

Authors: Stefan Riemelmoser, Merzuk Kaltak, Georg Kresse

Abstract: The optimized effective potential (OEP) method presents an unambiguous way to construct the Kohn-Sham potential corresponding to a given diagrammatic approximation for the exchange-correlation functional. The OEP from the random-phase approximation (RPA) has played an important role ever since the conception of the OEP formalism. However, the solution of the OEP equation is computationally fairly… ▽ More The optimized effective potential (OEP) method presents an unambiguous way to construct the Kohn-Sham potential corresponding to a given diagrammatic approximation for the exchange-correlation functional. The OEP from the random-phase approximation (RPA) has played an important role ever since the conception of the OEP formalism. However, the solution of the OEP equation is computationally fairly expensive and has to be done in a self-consistent way. So far, large scale solid state applications have therefore been performed only using the quasiparticle approximation (QPA), neglecting certain dynamical screening effects. We obtain the exact RPA-OEP for 15 semiconductors and insulators by direct solution of the linearized Sham-Schlüter equation. We investigate the accuracy of the QPA on Kohn-Sham band gaps and dielectric constants, and comment on the issue of self-consistency. △ Less

Submitted 16 July, 2023; v1 submitted 25 January, 2021; originally announced January 2021.

Comments: postprint

Journal ref: J. Chem. Phys. 154, 154103 (2021)

New insights into the 1D carbon chain through the RPA

Authors: Benjamin Ramberger, Georg Kresse

Abstract: We investigated the electronic and structural properties of the infinite linear carbon chain (carbyne) using density functional theory (DFT) and the random phase approximation (RPA) to the correlation energy. The studies are performed in vacuo and for carbyne inside a carbon nano tube (CNT). In the vacuum, semi-local DFT and RPA predict bond length alternations of about 0.04 Å and 0.13 Å, respecti… ▽ More We investigated the electronic and structural properties of the infinite linear carbon chain (carbyne) using density functional theory (DFT) and the random phase approximation (RPA) to the correlation energy. The studies are performed in vacuo and for carbyne inside a carbon nano tube (CNT). In the vacuum, semi-local DFT and RPA predict bond length alternations of about 0.04 Å and 0.13 Å, respectively. The frequency of the highest optical mode at the $Γ$ point is 1219 cm$^{-1}$ and about 2000 cm$^{-1}$ for DFT and the RPA. Agreement of the RPA to previous high level quantum chemistry and diffusion Monte-Carlo results is excellent. For the RPA we calculate the phonon-dispersion in the full Brillouine zone and find marked quantitative differences to DFT calculations not only at the $Γ$ point but also throughout the entire Brillouine zone. To model carbyne inside a carbon nanotube, we considered a (10,0) CNT. Here the DFT calculations are even qualitatively sensitive to the k-points sampling. At the limes of a very dense k-points sampling, semi-local DFT predicts no bond length alternation (BLA), whereas in the RPA a sizeable BLA of 0.09 Å prevails. The reduced BLA leads to a significant red shift of the vibrational frequencies of about 350 cm$^{-1}$, so that they are in good agreement with experimental estimates. Overall, the good agreement between the RPA and previously reported results from correlated wavefunction methods and experimental Raman data suggests that the RPA provides reliable results at moderate computational costs. It hence presents a useful addition to the repertoire of correlated wavefunction methods and its accuracy clearly prevails for low dimensional systems, where semi-local density functionals struggle to yield even qualitatively correct results. △ Less

Submitted 17 March, 2021; v1 submitted 18 December, 2020; originally announced December 2020.

Comments: 7 pages, 5 figures

Journal ref: Phys. Chem. Chem. Phys., 2021, 23, 5254

Local embedding of Coupled Cluster theory into the Random Phase Approximation using plane-waves

Authors: Tobias Schäfer, Florian Libisch, Georg Kresse, Andreas Grüneis

Abstract: We present an embedding approach to treat local electron correlation effects in periodic environments. In a single, consistent framework, our plane-wave based scheme embeds a local high-level correlation calculation (here Coupled Cluster Theory, CC), employing localized orbitals, into a low-level correlation calculation (here the direct Random Phase Approximation, RPA). This choice allows for an a… ▽ More We present an embedding approach to treat local electron correlation effects in periodic environments. In a single, consistent framework, our plane-wave based scheme embeds a local high-level correlation calculation (here Coupled Cluster Theory, CC), employing localized orbitals, into a low-level correlation calculation (here the direct Random Phase Approximation, RPA). This choice allows for an accurate and effcient treatment of long-range dispersion effects. Accelerated convergence with respect to the local fragment size can be observed if the low-level and high-level long-range dispersion are quantitatively similar, as is the case for CC in RPA. To demonstrate the capabilities of the introduced embedding approach, we calculate adsorption energies of molecules on a surface and in a chabazite crystal cage, as well as the formation energy of a lattice impurity in a solid at the level of highly accurate many-electron perturbation theories. The absorption energy of a methane molecule in a zeolite chabazite, for instance, is converged with an error well below 20 meV at the CC level. As our largest periodic benchmark system, we apply our scheme to the adsorption of a water molecule on titania in a supercell containing more than 1000 electrons. △ Less

Submitted 11 December, 2020; originally announced December 2020.

Comments: 6 pages, 2 figures

Journal ref: J. Chem. Phys. 154, 011101 (2021)

$α$-$β$ phase transition of zirconium predicted by on-the-fly machine-learned force field

Authors: Peitao Liu, Carla Verdi, Ferenc Karsai, Georg Kresse

Abstract: The accurate prediction of solid-solid structural phase transitions at finite temperature is a challenging task, since the dynamics is so slow that direct simulations of the phase transitions by first-principles (FP) methods are typically not possible. Here, we study the $α$-$β$ phase transition of Zr at ambient pressure by means of on-the-fly machine-learned force fields. These are automatically… ▽ More The accurate prediction of solid-solid structural phase transitions at finite temperature is a challenging task, since the dynamics is so slow that direct simulations of the phase transitions by first-principles (FP) methods are typically not possible. Here, we study the $α$-$β$ phase transition of Zr at ambient pressure by means of on-the-fly machine-learned force fields. These are automatically generated during FP molecular dynamics (MD) simulations without the need of human intervention, while retaining almost FP accuracy. Our MD simulations successfully reproduce the first-order displacive nature of the phase transition, which is manifested by an abrupt jump of the volume and a cooperative displacement of atoms at the phase transition temperature. The phase transition is further identified by the simulated x-ray powder diffraction, and the predicted phase transition temperature is in reasonable agreement with experiment. Furthermore, we show that using a singular value decomposition and pseudo inversion of the design matrix generally improves the machine-learned force field compared to the usual inversion of the squared matrix in the regularized Bayesian regression. △ Less

Submitted 24 May, 2021; v1 submitted 8 December, 2020; originally announced December 2020.

Comments: 10 pages, 7 figures

Journal ref: Phys. Rev. Materials 5, 053804 (2021)

Nonrad: Computing Nonradiative Capture Coefficients from First Principles

Authors: Mark E. Turiansky, Audrius Alkauskas, Manuel Engel, Georg Kresse, Darshana Wickramaratne, Jimmy-Xuan Shen, Cyrus E. Dreyer, Chris G. Van de Walle

Abstract: Point defects in semiconductor crystals provide a means for carriers to recombine nonradiatively. This recombination process impacts the performance of devices. We present the Nonrad code that implements the first-principles approach of Alkauskas et al. [Phys. Rev. B 90, 075202 (2014)] for the evaluation of nonradiative capture coefficients based on a quantum-mechanical description of the capture… ▽ More Point defects in semiconductor crystals provide a means for carriers to recombine nonradiatively. This recombination process impacts the performance of devices. We present the Nonrad code that implements the first-principles approach of Alkauskas et al. [Phys. Rev. B 90, 075202 (2014)] for the evaluation of nonradiative capture coefficients based on a quantum-mechanical description of the capture process. An approach for evaluating electron-phonon coupling within the projector augmented wave formalism is presented. We also show that the common procedure of replacing Dirac delta functions with Gaussians can introduce errors into the resulting capture rate, and implement an alternative scheme to properly account for vibrational broadening. Lastly, we assess the accuracy of using an analytic approximation to the Sommerfeld parameter by comparing with direct numerical evaluation. △ Less

Submitted 22 November, 2020; v1 submitted 14 November, 2020; originally announced November 2020.

Comparative ab initio study of the structural, electronic, magnetic, and dynamical properties of LiOsO$_3$ and NaOsO$_3$

Authors: Peitao Liu, Jiangang He, Bongjae Kim, Sergii Khmelevskyi, Alessandro Toschi, Georg Kresse, Cesare Franchini

Abstract: Despite similar chemical compositions, LiOsO$_3$ and NaOsO$_3$ exhibit remarkably distinct structural, electronic, magnetic, and spectroscopic properties. At low temperature, LiOsO$_3$ is a polar bad metal with a rhombohedral $R3c$ structure without the presence of long-range magnetic order, whereas NaOsO$_3$ is a $G$-type antiferromagnetic insulator with an orthorhombic $Pnma$ structure. By means… ▽ More Despite similar chemical compositions, LiOsO$_3$ and NaOsO$_3$ exhibit remarkably distinct structural, electronic, magnetic, and spectroscopic properties. At low temperature, LiOsO$_3$ is a polar bad metal with a rhombohedral $R3c$ structure without the presence of long-range magnetic order, whereas NaOsO$_3$ is a $G$-type antiferromagnetic insulator with an orthorhombic $Pnma$ structure. By means of comparative first-principles DFT+$U$ calculations with the inclusion of the spin-orbit coupling, we ($i$) identify the origin of the different structural ($R3c$ vs. $Pnma$) properties using a symmetry-adapted soft mode analysis, ($ii$) provide evidence that all considered exchange-correlation functionals (LDA, PBE, PBEsol, SCAN, and HSE06) and the spin disordered polymorphous descriptions are unsatisfactory to accurately describe the electronic and magnetic properties of both systems simultaneously, and ($iii$) clarify that the distinct electronic (metallic vs. insulating) properties originates mainly from a cooperative steric and magnetic effect. Finally, we find that although at ambient pressure LiOsO$_3$ with a $Pnma$ symmetry and NaOsO$_3$ with a $R\bar{3}c$ symmetry are energetically unfavorable, they do not show soft phonons and therefore are dynamically stable. A pressure-induced structural phase transition from $R3c$ to $Pnma$ for LiOsO$_3$ is predicted, whereas for NaOsO$_3$ no symmetry change is discerned in the considered pressure range. △ Less

Submitted 14 February, 2020; originally announced February 2020.

Comments: 9 pages, 9 figures

Journal ref: Phys. Rev. Materials 4, 045001 (2020)

Accurate optical spectra through time-dependent density functional theory based on screening-dependent hybrid functionals

Authors: Alexey Tal, Peitao Liu, Georg Kresse, Alfredo Pasquarello

Abstract: We investigate optical absorption spectra obtained through time-dependent density functional theory (TD-DFT) based on nonempirical hybrid functionals that are designed to correctly reproduce the dielectric function. The comparison with state-of-the-art $GW$ calculations followed by the solution of the Bethe-Sapeter equation (BSE-$GW$) shows close agreement for both the transition energies and the… ▽ More We investigate optical absorption spectra obtained through time-dependent density functional theory (TD-DFT) based on nonempirical hybrid functionals that are designed to correctly reproduce the dielectric function. The comparison with state-of-the-art $GW$ calculations followed by the solution of the Bethe-Sapeter equation (BSE-$GW$) shows close agreement for both the transition energies and the main features of the spectra. We confront TD-DFT with BSE-$GW$ by focusing on the model dielectric function and the local exchange-correlation kernel. The present TD-DFT approach achieves the accuracy of BSE-$GW$ at a fraction of the computational cost. △ Less

Submitted 14 August, 2020; v1 submitted 6 February, 2020; originally announced February 2020.

Journal ref: Phys. Rev. Research 2, 032019 (2020)

Electron-Phonon Interactions Using the PAW Method and Wannier Functions

Authors: Manuel Engel, Martijn Marsman, Cesare Franchini, Georg Kresse

Abstract: We present an ab-initio density-functional-theory approach for calculating electron-phonon interactions within the projector augmented-wave method. The required electron-phonon matrix elements are defined as the second derivative of the one-electron energies in the PAW method. As the PAW method leads to a generalized eigenvalue problem, the resulting electron-phonon matrix elements lack some symme… ▽ More We present an ab-initio density-functional-theory approach for calculating electron-phonon interactions within the projector augmented-wave method. The required electron-phonon matrix elements are defined as the second derivative of the one-electron energies in the PAW method. As the PAW method leads to a generalized eigenvalue problem, the resulting electron-phonon matrix elements lack some symmetries that are usually present for simple eigenvalue problems and all-electron formulations. We discuss the relation between our definition of the electron-phonon matrix element and other formulations. To allow for efficient evaluation of physical properties, we introduce a Wannier-interpolation scheme, again adapted to generalized eigenvalue problems. To explore the method's numerical characteristics, the temperature-dependent band-gap renormalization of diamond is calculated and compared with previous publications. Furthermore, we apply the method to selected binary compounds and show that the obtained zero-point renormalizations agree well with other values found in literature and experiments. △ Less

Submitted 23 January, 2020; originally announced January 2020.

Journal ref: Phys. Rev. B 101, 184302 (2020)

Plane wave basis set correction methods for RPA correlation energies

Authors: Stefan Riemelmoser, Merzuk Kaltak, Georg Kresse

Abstract: Electronic correlation energies from the random-phase approximation converge slowly with respect to the plane wave basis set size. We study the conditions, under which a short-range local density functional can be used to account for the basis set incompleteness error. Furthermore, we propose a one-shot extrapolation scheme based on the Lindhard response function of the homogeneous electron gas. T… ▽ More Electronic correlation energies from the random-phase approximation converge slowly with respect to the plane wave basis set size. We study the conditions, under which a short-range local density functional can be used to account for the basis set incompleteness error. Furthermore, we propose a one-shot extrapolation scheme based on the Lindhard response function of the homogeneous electron gas. The different basis set correction methods are used to calculate equilibrium lattice constants for prototypical solids of different bonding types. △ Less

Submitted 16 July, 2023; v1 submitted 22 January, 2020; originally announced January 2020.

Comments: postprint, typographical errors corrected: lower frequency integration limit in Eqs. (13) and (B1) [correct is 0 instead of minus infinity]; corrected minus signs for Eqs. (B9) and (B10)

Journal ref: J. Chem. Phys. 152, 134103 (2020)

Optical response of an interacting polaron gas in strongly polar crystals

Authors: S. N. Klimin, J. Tempere, J. T. Devreese, C. Franchini, G. Kresse

Abstract: Optical conductivity of an interacting polaron gas is calculated within an extended random phase approximation which takes into account mixing of collective excitations of the electron gas with LO phonons. This mixing is important for the optical response of strongly polar crystals where the static dielectric constant is rather high: strontium titanate is the case. The present calculation sheds li… ▽ More Optical conductivity of an interacting polaron gas is calculated within an extended random phase approximation which takes into account mixing of collective excitations of the electron gas with LO phonons. This mixing is important for the optical response of strongly polar crystals where the static dielectric constant is rather high: strontium titanate is the case. The present calculation sheds light on unexplained features of experimentally observed optical conductivity spectra in $n$-doped SrTiO$_{3}$. These features appear to be due to dynamic screening of the electron-electron interaction by polar optical phonons and hence do not require additional mechanisms for the explanation. △ Less

Submitted 6 December, 2019; originally announced December 2019.

Comments: 8 pages, 1 figure

RPA natural orbitals and their application to post-Hartree-Fock electronic structure methods

Authors: Benjamin Ramberger, Zoran Sukurma, Tobias Schäfer, Georg Kresse

Abstract: We present a method to approximate post-Hartree-Fock correlation energies by using approximate natural orbitals obtained by the random phase approximation (RPA). We demonstrate the method by applying it to the helium atom, the hydrogen and fluorine molecule, and to diamond as an example of a periodic system. For these benchmark systems, we show that RPA natural orbitals converge the MP2 correlatio… ▽ More We present a method to approximate post-Hartree-Fock correlation energies by using approximate natural orbitals obtained by the random phase approximation (RPA). We demonstrate the method by applying it to the helium atom, the hydrogen and fluorine molecule, and to diamond as an example of a periodic system. For these benchmark systems, we show that RPA natural orbitals converge the MP2 correlation energy rapidly. Additionally, we calculated full configuration interaction energies for He and H$_2$, which are in excellent agreement with the literature and experimental values. We conclude that the proposed method may serve as a compromise to reach good approximations to correlation energies at moderate computational cost, and we expect the method to be especially useful for theoretical studies on surface chemistry by providing an efficient basis to correlated wave function based methods. △ Less

Submitted 18 March, 2021; v1 submitted 16 September, 2019; originally announced September 2019.

Comments: 10 pages, 18 figures, 16 Equations; v2: fixed typo in metadata abstract, changed bibliography style; v3: added reference (4), added Acknowledgement

Minimax Isometry Method: A compressive sensing approach for Matsubara summation in many-body perturbation theory

Authors: Merzuk Kaltak, Georg Kresse

Abstract: We present a compressive sensing approach for the long standing problem of Matsubara summation in many-body perturbation theory. By constructing low-dimensional, almost isometric subspaces of the Hilbert space we obtain optimum imaginary time and frequency grids that allow for extreme data compression of fermionic and bosonic functions in a broad temperature regime. The method is applied to the ra… ▽ More We present a compressive sensing approach for the long standing problem of Matsubara summation in many-body perturbation theory. By constructing low-dimensional, almost isometric subspaces of the Hilbert space we obtain optimum imaginary time and frequency grids that allow for extreme data compression of fermionic and bosonic functions in a broad temperature regime. The method is applied to the random phase and self-consistent $GW$ approximation of the grand potential. Integration and transformation errors are investigated for Si and SrVO$_3$. △ Less

Submitted 10 March, 2020; v1 submitted 4 September, 2019; originally announced September 2019.

Comments: 19 pages, 11 figures

On-the-fly machine learning force field generation: Application to melting points

Authors: Ryosuke Jinnouchi, Ferenc Karsai, Georg Kresse

Abstract: An efficient and robust on-the-fly machine learning force field method is developed and integrated into an electronic-structure code. This method realizes automatic generation of machine learning force fields on the basis of Bayesian inference during molecular dynamics simulations, where the first principles calculations are only executed, when new configurations out of already sampled datasets ap… ▽ More An efficient and robust on-the-fly machine learning force field method is developed and integrated into an electronic-structure code. This method realizes automatic generation of machine learning force fields on the basis of Bayesian inference during molecular dynamics simulations, where the first principles calculations are only executed, when new configurations out of already sampled datasets appear. The developed method is applied to the calculation of melting points of Al, Si, Ge, Sn and MgO. The applications indicate that more than 99 \% of the first principles calculations are bypassed during the force field generation. This allows the machine to quickly construct first principles datasets over wide phase spaces. Furthermore, with the help of the generated machine learning force fields, simulations are accelerated by a factor of thousand compared with first principles calculations. Accuracies of the melting points calculated by the force fields are examined by thermodynamic perturbation theory, and the examination indicates that the machine learning force fields can quantitatively reproduce the first principles melting points. △ Less

Submitted 5 May, 2019; v1 submitted 29 April, 2019; originally announced April 2019.

Comments: 15 pages, 7 figures and 2 tables

Journal ref: Phys. Rev. B 100, 014105 (2019)

Phase transitions of hybrid perovskites simulated by machine-learning force fields trained on-the-fly with Bayesian inference

Authors: Ryosuke Jinnouchi, Jonathan Lahnsteiner, Ferenc Karsai, Georg Kresse, Menno Bokdam

Abstract: Realistic finite temperature simulations of matter are a formidable challenge for first principles methods. Long simulation times and large length scales are required, demanding years of compute time. Here we present an on-the-fly machine learning scheme that generates force fields automatically during molecular dynamics simulations. This opens up the required time and length scales, while retaini… ▽ More Realistic finite temperature simulations of matter are a formidable challenge for first principles methods. Long simulation times and large length scales are required, demanding years of compute time. Here we present an on-the-fly machine learning scheme that generates force fields automatically during molecular dynamics simulations. This opens up the required time and length scales, while retaining the distinctive chemical precision of first principles methods and minimizing the need for human intervention. The method is widely applicable to multi-element complex systems. We demonstrate its predictive power on the entropy driven phase transitions of hybrid perovskites, which have never been accurately described in simulations. Using machine learned potentials, isothermal-isobaric simulations give direct insight into the underlying microscopic mechanisms. Finally, we relate the phase transition temperatures of different perovskites to the radii of the involved species, and we determine the order of the transitions in Landau theory. △ Less

Submitted 22 March, 2019; originally announced March 2019.

Journal ref: Phys. Rev. Lett. 122, 225701 (2019)

First-principle study of the melting temperature of MgO

Authors: Max Rang, Georg Kresse

Abstract: Using first-principles only, we calculate the melting point of MgO, also called periclase or magnesia. The random phase approximation (RPA) is used to include the exact exchange as well as local and non-local many-body correlation terms, in order to provide high accuracy. Using the free energy method, we obtain the melting temperature directly from the internal energies calculated with DFT. The fr… ▽ More Using first-principles only, we calculate the melting point of MgO, also called periclase or magnesia. The random phase approximation (RPA) is used to include the exact exchange as well as local and non-local many-body correlation terms, in order to provide high accuracy. Using the free energy method, we obtain the melting temperature directly from the internal energies calculated with DFT. The free energy differences between the ensembles generated by the molecular dynamics simulations are calulated with thermodynamic integration or thermodynamic perturbation theory. The predicted melting temperature is $T_m^\text{RPA} =3043\pm 86$ K, and the values obtained with the PBE and SCAN functionals are $T_m^\text{PBE} = 2747 \pm 59$ K and $T_m^\text{SCAN} = 3032\pm 53 $ K. △ Less

Submitted 4 April, 2019; v1 submitted 21 February, 2019; originally announced February 2019.

Comments: 5 pages, 4 figures

Journal ref: Phys. Rev. B 99, 184103 (2019)

Assessing model-dielectric-dependent hybrid functionals on the antiferromagnetic transition-metal monoxides MnO, FeO, CoO, and NiO

Authors: Peitao Liu, Cesare Franchini, Martijn Marsman, Georg Kresse

Abstract: Recently, two nonempirical hybrid functionals, dielectric-dependent range-separated hybrid functional based on the Coulomb-attenuating method (DD-RSH-CAM) and doubly screened hybrid functional (DSH), have been suggested by [Chen et al, Phys. Rev. Mater. 2, 073803 (2018)] and [Cui et al, J. Phys. Chem. Lett. 9, 2338 (2018)], respectively. These two hybrid functionals are both based on a common mode… ▽ More Recently, two nonempirical hybrid functionals, dielectric-dependent range-separated hybrid functional based on the Coulomb-attenuating method (DD-RSH-CAM) and doubly screened hybrid functional (DSH), have been suggested by [Chen et al, Phys. Rev. Mater. 2, 073803 (2018)] and [Cui et al, J. Phys. Chem. Lett. 9, 2338 (2018)], respectively. These two hybrid functionals are both based on a common model dielectric function approach, but differ in the way how to non-empirically obtain the range-separation parameter. By retaining the full short-range Fock exchange and a fraction of the long-range Fock exchange that equals the inverse of the dielectric constant, both DD-RSH-CAM and DSH turn out to perform very well in predicting the band gaps for a large variety of semiconductors and insulators. Here, we assess how these two hybrid functionals perform on challenging antiferromagnetic transition-metal monoxides MnO, FeO, CoO, and NiO by comparing them to other conventional hybrid functionals and the $GW$ method. We find that single-shot DD0-RSH-CAM and DSH0 improve the band gaps towards experiments as compared to conventional hybrid functionals. The magnetic moments are slightly increased, but the predicted dielectric constants are decreased. The valence band density of states (DOS) predicted by DD0-RSH-CAM and DSH0 are as satisfactory as HSE03 in comparison to experimental spectra, however, the conduction band DOS are shifted to higher energies by about 2 eV compared to HSE03. Self-consistent DD-RSH-CAM and DSH deteriorate the results with a significant overestimation of band gaps. △ Less

Submitted 20 August, 2019; v1 submitted 13 February, 2019; originally announced February 2019.

Comments: 11 pages, 5 figures

Journal ref: J. Phys.: Condens. Matter 32 (2019) 015502

Small Polarons in Transition Metal Oxides

Authors: Michele Reticcioli, Ulrike Diebold, Georg Kresse, Cesare Franchini

Abstract: The formation of polarons is a pervasive phenomenon in transition metal oxide compounds, with a strong impact on the physical properties and functionalities of the hosting materials. In its original formulation the polaron problem considers a single charge carrier in a polar crystal interacting with its surrounding lattice. Depending on the spatial extension of the polaron quasiparticle, originati… ▽ More The formation of polarons is a pervasive phenomenon in transition metal oxide compounds, with a strong impact on the physical properties and functionalities of the hosting materials. In its original formulation the polaron problem considers a single charge carrier in a polar crystal interacting with its surrounding lattice. Depending on the spatial extension of the polaron quasiparticle, originating from the coupling between the excess charge and the phonon field, one speaks of small or large polarons. This chapter discusses the modeling of small polarons in real materials, with a particular focus on the archetypal polaron material TiO2. After an introductory part, surveying the fundamental theoretical and experimental aspects of the physics of polarons, the chapter examines how to model small polarons using first principles schemes in order to predict, understand and interpret a variety of polaron properties in bulk phases and surfaces. Following the spirit of this handbook, different types of computational procedures and prescriptions are presented with specific instructions on the setup required to model polaron effects. △ Less

Submitted 15 February, 2019; v1 submitted 11 February, 2019; originally announced February 2019.

Comments: 36 pages, 12 figures

Journal ref: In: Andreoni W., Yip S. (eds) Handbook of Materials Modeling. Springer, Cham. (2019)

Superconductivity in SrTiO$_{3}$: dielectric function method for non-parabolic bands

Authors: S. N. Klimin, J. Tempere, J. T. Devreese, J. He, C. Franchini, G. Kresse

Abstract: The dielectric function method for superconductivity has been applied to SrTiO$_{3}$ accounting for the non-parabolic dispersion of charge carriers in the conduction band and for the dispersion of optical phonons based on density functional theory calculations. The obtained critical temperatures of the superconducting phase transition in SrTiO$_{3}$ are in agreement with experiments in the density… ▽ More The dielectric function method for superconductivity has been applied to SrTiO$_{3}$ accounting for the non-parabolic dispersion of charge carriers in the conduction band and for the dispersion of optical phonons based on density functional theory calculations. The obtained critical temperatures of the superconducting phase transition in SrTiO$_{3}$ are in agreement with experiments in the density range $n\sim5\times10^{18}$ to $5\times10^{20}$ cm$^{-3}$. The dielectric function method predicts also the sign of the anomalous isotope effect in strontium titanate, in line with recent observations. △ Less

Submitted 28 November, 2018; originally announced November 2018.

Comments: 5 pages, 2 figures, submitted to the Special Issue of the Journal of Superconductivity and Novel Magnetism

Journal ref: Journal of Superconductivity and Novel Magnetism 32, 2739-2744 (2019)

On the physisorption of water on graphene: Sub-chemical accuracy from many-body electronic structure methods

Authors: Jan Gerit Brandenburg, Andrea Zen, Martin Fitzner, Benjamin Ramberger, Georg Kresse, Theodoros Tsatsoulis, Andreas Grüneis, Angelos Michaelides, Dario Alfè

Abstract: Molecular adsorption on surfaces plays a central role in catalysis, corrosion, desalination, and many other processes of relevance to industry and the natural world. Few adsorption systems are more ubiquitous or of more widespread importance than those involving water and carbon, and for a molecular level understanding of such interfaces water monomer adsorption on graphene is a fundamental and re… ▽ More Molecular adsorption on surfaces plays a central role in catalysis, corrosion, desalination, and many other processes of relevance to industry and the natural world. Few adsorption systems are more ubiquitous or of more widespread importance than those involving water and carbon, and for a molecular level understanding of such interfaces water monomer adsorption on graphene is a fundamental and representative system. This system is particularly interesting as it calls for an accurate treatment of electron correlation effects, as well as posing a practical challenge to experiments. Here, we employ many-body electronic structure methodologies that can be rigorously converged and thus provide faithful references for the molecule-surface interaction. In particular, we use diffusion Monte-Carlo (DMC), coupled cluster (CCSD(T)), as well as the random phase approximation (RPA) to calculate the strength of the interaction between water and an extended graphene surface. We establish excellent, sub-chemical, agreement between the complementary high-level methodologies, and an adsorption energy estimate in the most stable configuration of approximately -100\,meV is obtained. We also find that the adsorption energy is rather insensitive to the orientation of the water molecule on the surface, despite different binding motifs involving qualitatively different interfacial charge reorganisation. In producing the first demonstrably accurate adsorption energies for water on graphene this work also resolves discrepancies amongst previously reported values for this widely studied system. It also paves the way for more accurate and reliable studies of liquid water at carbon interfaces with cheaper computational methods, such as density functional theory and classical potentials. △ Less

Submitted 19 November, 2018; originally announced November 2018.

Journal ref: J. Phys. Chem. Lett., 2019, 10, pp 358-368

Correlated excited states in the narrow gap band semiconductor FeSi and antiferromagnetic screening of local spin moments

Authors: Sergii Khmelevskyi, Georg Kresse, Peter Mohn

Abstract: The physical properties of the semiconductor FeSi with very narrow band gap, anomalous behavior of the magnetic susceptibility and metal-insulator transition at elevated temperatures attract gross interest due to the still controversial theoretical understanding of their origin. On one side the purely band like mechanism of the gap formation in FeSi at low temperature is well established, on other… ▽ More The physical properties of the semiconductor FeSi with very narrow band gap, anomalous behavior of the magnetic susceptibility and metal-insulator transition at elevated temperatures attract gross interest due to the still controversial theoretical understanding of their origin. On one side the purely band like mechanism of the gap formation in FeSi at low temperature is well established, on other side a number of experiments and their theoretical interpretation suggest a rich physics of strong correlations at finite temperature. In this work we use an ab-initio scheme based on the Random Phase Approximation and Local Spin Density Approximation (RPA@LSDA) to reveal the role of the electron correlation effects in FeSi extending it by applying a fixed spin moment constraint. In the parameter free framework we show that correlation effects essentially alter the one-electron LSDA results leading to the formation of an additional state with finite magnetic moment on Fe, whose energy is almost degenerate with the non-magnetic ground state. This explains the results of high field experiments, which found a first-order meta-magnetic phase transition into a metallic ferromagnetic state. Our results suggest a strongly correlated nature of the low-energy excitations in FeSi. From our super-cells calculations we reveal that these excitations are local and exhibit a Kondo-like behavior since a strong antiferromagnetic screening is present. △ Less

Submitted 6 September, 2018; originally announced September 2018.

Comments: accepted in Phys. Rev. B

Journal ref: Physical Review B 98, 125205 (2018)

Beyond the quasiparticle approximation: Fully self-consistent $GW$ calculations

Authors: Manuel Grumet, Peitao Liu, Merzuk Kaltak, Jiří Klimeš, Georg Kresse

Abstract: We present quasiparticle (QP) energies from fully self-consistent $GW$ (sc$GW$) calculations for a set of prototypical semiconductors and insulators within the framework of the projector-augmented wave methodology. To obtain converged results, both finite basis-set corrections and $k$-point corrections are included, and a simple procedure is suggested to deal with the singularity of the Coulomb ke… ▽ More We present quasiparticle (QP) energies from fully self-consistent $GW$ (sc$GW$) calculations for a set of prototypical semiconductors and insulators within the framework of the projector-augmented wave methodology. To obtain converged results, both finite basis-set corrections and $k$-point corrections are included, and a simple procedure is suggested to deal with the singularity of the Coulomb kernel in the long-wavelength limit, the so called head correction. It is shown that the inclusion of the head corrections in the sc$GW$ calculations is critical to obtain accurate QP energies with a reasonable $k$-point set. We first validate our implementation by presenting detailed results for the selected case of diamond, and then we discuss the converged QP energies, in particular the band gaps, for a set of gapped compounds and compare them to single-shot $G_0W_0$, QP self-consistent $GW$, and previously available sc$GW$ results as well as experimental results. △ Less

Submitted 29 October, 2018; v1 submitted 23 August, 2018; originally announced August 2018.

Comments: 9 pages, 5 figures

Journal ref: Phys. Rev. B 98, 155143 (2018)

Melting Si: beyond density functional theory

Authors: Florian Dorner, Zoran Sukurma, Christoph Dellago, Georg Kresse

Abstract: The melting point of silicon in the cubic diamond phase is calculated using the random phase approximation (RPA). The RPA includes exact exchange as well as an approximate treatment of local as well as non-local many body correlation effects of the electrons. We predict a melting temperature of about 1735 K and 1640 K without and with core polarization effects, respectively. Both values are within… ▽ More The melting point of silicon in the cubic diamond phase is calculated using the random phase approximation (RPA). The RPA includes exact exchange as well as an approximate treatment of local as well as non-local many body correlation effects of the electrons. We predict a melting temperature of about 1735 K and 1640 K without and with core polarization effects, respectively. Both values are within 3 % of the experimental melting temperature of 1687 K. In comparison, the commonly used gradient approximation to density functional theory predicts a melting point that is 200 K too low, and hybrid functionals overestimate the melting point by 150 K. We correlate the predicted melting point with the energy difference between cubic diamond and the beta-tin phase of silicon, establishing that this energy difference is an important benchmark for the development of approximate functionals. The current results establish that the RPA can be used to predict accurate finite temperature properties and underlines the excellent predictive properties of the RPA for condensed matter. △ Less

Submitted 24 September, 2018; v1 submitted 6 August, 2018; originally announced August 2018.

Comments: 5 pages, 3 figures

Journal ref: Phys. Rev. Lett. 121, 195701 (2018)

Ab initio phase diagram of PbSe crystals calculated with the Random Phase Approximation

Authors: Tobias Schäfer, Zhaochuan Fan, Michael Grünwald, Georg Kresse

Abstract: Understanding the phase behavior of semiconductor materials is important for applications in solid state physics and nanoscience. Accurate experimental data is often difficult to obtain due to strong kinetic effects. In this work, we calculate the temperature-pressure phase diagram of lead selenide (PbSe) using the random phase approximation (RPA), an accurate wavefunction based many-body techniqu… ▽ More Understanding the phase behavior of semiconductor materials is important for applications in solid state physics and nanoscience. Accurate experimental data is often difficult to obtain due to strong kinetic effects. In this work, we calculate the temperature-pressure phase diagram of lead selenide (PbSe) using the random phase approximation (RPA), an accurate wavefunction based many-body technique. We consider three crystalline phases, the low pressure B1 phase (NaCl-type), the intermediate B33 phase (CrB-type), and the high pressure B2 phase (CsCl-type). The electronic contributions to the free energy (at T=0K) are calculated in the Born-Oppenheimer approximation using the RPA, whereas phononic contributions are computed in the quasi-harmonic approximation using DFT and the PBEsol functional. At room temperature, we find transition pressures of 4.6 +/- 0.3 GPa for the B1-B33 transition and 18.7 +/- 0.3 GPa for the B33-B2 transition, in good agreement with experiments. In contrast to the interpretation of recent experiments, we observe a negative Clapeyron slope for both transitions. Gibbs free energy differences between competing structures have small gradients close to coexistence, consistent with pronounced hysteresis observed in experiments. The phase diagram presented in this work can serve as a reference for future studies of PbSe and should prove useful in the development of accurate and efficient force fields. △ Less

Submitted 30 July, 2018; originally announced July 2018.

Comments: 7 pages, 7 figures

Journal ref: Phys. Rev. B 98, 144103 (2018)

Relativistic $GW$+BSE study of the optical properties of Ruddlesden-Popper iridates

Authors: Peitao Liu, Bongjae Kim, Xing-Qiu Chen, D. D. Sarma, Georg Kresse, Cesare Franchini

Abstract: We study the optical properties of the Ruddlesden-Popper series of iridates Sr$_{n+1}$Ir$_n$O$_{3n+1}$ ($n$=1, 2 and $\infty$) by solving the Bethe-Salpeter equation (BSE), where the quasiparticle (QP) energies and screened interactions $W$ are obtained by the $GW$ approximation including spin-orbit coupling. The computed optical conductivity spectra show strong excitonic effects and reproduce ver… ▽ More We study the optical properties of the Ruddlesden-Popper series of iridates Sr$_{n+1}$Ir$_n$O$_{3n+1}$ ($n$=1, 2 and $\infty$) by solving the Bethe-Salpeter equation (BSE), where the quasiparticle (QP) energies and screened interactions $W$ are obtained by the $GW$ approximation including spin-orbit coupling. The computed optical conductivity spectra show strong excitonic effects and reproduce very well the experimentally observed double-peak structure, in particular for the spin-orbital Mott insulators Sr$_2$IrO$_4$ and Sr$_3$Ir$_2$O$_7$. However, $GW$ does not account well for the correlated metallic state of SrIrO$_3$ owing to a much too small band renormalization, and this affects the overall quality of the optical conductivity. Our analysis describes well the progressive redshift of the main optical peaks as a function of dimensionality ($n$), which is correlated with the gradual decrease of the electronic correlation (quantified by the constrained random phase approximation) towards the metallic $n=\infty$ limit. We have also assessed the quality of a computationally cheaper BSE approach that is based on a model dielectric function and conducted on top of DFT+$U$ one-electron energies. Unfortunately, this model BSE approach does not accurately reproduce the outcome of the full $GW$+BSE method and leads to larger deviations to the measured spectra. △ Less

Submitted 31 July, 2018; v1 submitted 19 June, 2018; originally announced June 2018.

Comments: 13 pages, 8 figures

Journal ref: Phys. Rev. Materials 2, 075003 (2018)

The Finite Temperature Structure of the MAPbI3 Perovskite: Comparing Density Functional Approximations and Force Fields to Experiment

Authors: Jonathan Lahnsteiner, Georg Kresse, Jurn Heinen, Menno Bokdam

Abstract: Determining the finite temperature structure of the hybrid perovskite MAPbI3 is a challenge for both experimental and theoretical methods. A very powerful computational method that can resolve the atomic structure is molecular dynamics (MD). The resulting structure depends on the density functional approximation (DFA) in the case of ab initio MD and the force field in classical MD. We compare the… ▽ More Determining the finite temperature structure of the hybrid perovskite MAPbI3 is a challenge for both experimental and theoretical methods. A very powerful computational method that can resolve the atomic structure is molecular dynamics (MD). The resulting structure depends on the density functional approximation (DFA) in the case of ab initio MD and the force field in classical MD. We compare the structure between 250K and 400K obtained with different DFAs and force fields in one consistent manner. The symmetry of the PbI3 framework is analyzed as well as the relative ordering of the neighboring organic molecules inside the framework. The distribution function of the molecules is used to map out an effective energy surface for the rotation of a single molecule. This surface is accurately modeled by a pair of cubic harmonics. Available experimental data in literature are discussed and compared to the structure obtained with the different methods. The spread in these data is still too large to uniquely determine the method that 'best' describes the perovskite, however promising candidates and outliers have been identified. △ Less

Submitted 23 April, 2018; originally announced April 2018.

Comments: 15 pages, 8 figures

Journal ref: Phys. Rev. Materials 2, 073604 (2018)

Charge Trapping at the Step Edges of TiO2 Anatase (101)

Authors: Martin Setvin, Xianfeng Hao, Benjamin Daniel, Jiri Pavelec, Zbynek Novotny, Gareth S. Parkinson, Michael Schmid, Georg Kresse, Cesare Franchini, Ulrike Diebold

Abstract: A combination of Photoemission, Atomic Force and Scanning Tunneling Microscopy/Spectroscopy measurements shows that excess electrons in TiO2 anatase (101) surface are trapped at step edges. Consequently, steps act as preferred adsorption sites for O2. In Density Functional Theory calculations electrons localize at clean step edges, this tendency is enhanced by O vacancies and hydroxylation. The re… ▽ More A combination of Photoemission, Atomic Force and Scanning Tunneling Microscopy/Spectroscopy measurements shows that excess electrons in TiO2 anatase (101) surface are trapped at step edges. Consequently, steps act as preferred adsorption sites for O2. In Density Functional Theory calculations electrons localize at clean step edges, this tendency is enhanced by O vacancies and hydroxylation. The results show the importance of defects for the wide-ranging applications of titania. △ Less

Submitted 20 April, 2018; originally announced April 2018.

Journal ref: Angew. Chem. Intl. Ed. 53, 4714 (2014)

Charge self-consistent many-body corrections using optimized projected localized orbitals

Authors: Malte Schüler, Oleg E. Peil, Gernot J. Kraberger, Ronald Pordzik, Martijn Marsman, Georg Kresse, Tim O. Wehling, Markus Aichhorn

Abstract: In order for methods combining ab initio density-functional theory and many-body techniques to become routinely used, a flexible, fast, and easy-to-use implementation is crucial. We present an implementation of a general charge self-consistent scheme based on projected localized orbitals in the projector augmented wave framework in the Vienna Ab Initio Simulation Package (VASP). We give a detailed… ▽ More In order for methods combining ab initio density-functional theory and many-body techniques to become routinely used, a flexible, fast, and easy-to-use implementation is crucial. We present an implementation of a general charge self-consistent scheme based on projected localized orbitals in the projector augmented wave framework in the Vienna Ab Initio Simulation Package (VASP). We give a detailed description on how the projectors are optimally chosen and how the total energy is calculated. We benchmark our implementation in combination with dynamical mean-field theory: first we study the charge-transfer insulator NiO using a Hartree-Fock approach to solve the many-body Hamiltonian. We address the advantages of the optimized against non-optimized projectors and furthermore find that charge self-consistency decreases the dependence of the spectral function - especially the gap - on the double counting. Second, using continuous-time quantum Monte Carlo we study a monolayer of SrVO$_3$, where strong orbital polarization occurs due to the reduced dimensionality. Using total-energy calculation for structure determination, we find that electronic correlations have a non-negligible influence on the position of the apical oxygens, and therefore on the thickness of the single SrVO$_3$ layer. △ Less

Submitted 5 November, 2018; v1 submitted 5 April, 2018; originally announced April 2018.

Comments: 11 pages, 6 figures

Journal ref: J. Phys.: Condens. Matter 30, 475901 (2018)

Laplace transformed MP2 for three dimensional periodic materials using stochastic orbitals in the plane wave basis and correlated sampling

Authors: Tobias Schäfer, Benjamin Ramberger, Georg Kresse

Abstract: We present an implementation and analysis of a stochastic high performance algorithm to calculate the correlation energy of three dimensional periodic systems in second-order Møller-Plesset perturbation theory (MP2). In particular we measure the scaling behavior of the sample variance and probe whether this stochastic approach is competitive if accuracies well below 1 meV per valence orbital are r… ▽ More We present an implementation and analysis of a stochastic high performance algorithm to calculate the correlation energy of three dimensional periodic systems in second-order Møller-Plesset perturbation theory (MP2). In particular we measure the scaling behavior of the sample variance and probe whether this stochastic approach is competitive if accuracies well below 1 meV per valence orbital are required, as it is necessary for calculations of adsorption, binding, or surface energies. The algorithm is based on the Laplace transformed MP2 (LTMP2) formulation in the plane wave basis. The time-dependent Hartree-Fock orbitals, appearing in the LTMP2 formulation, are stochastically rotated in the occupied and unoccupied Hilbert space. This avoids a full summation over all combinations of occupied and unoccupied orbitals, as inspired by the work of D. Neuhauser, E. Rabani, and R. Baer in J. Chem. Theory Comput. 9, 24 (2013). Additionally, correlated sampling is introduced, accelerating the statistical convergence significantly. △ Less

Submitted 9 February, 2018; v1 submitted 15 November, 2017; originally announced November 2017.

Comments: 11 pages, 6 figures

Journal ref: J. Chem. Phys. 148, 064103 (2018)

Converged GW quasiparticle energies for transition metal oxide perovskites

Authors: Zeynep Ergonenc, Bongjae Kim, Peitao Liu, Georg Kresse, Cesare Franchini

Abstract: The ab initio calculation of quasiparticle (QP) energies is a technically and computationally challenging problem. In condensed matter physics the most widely used approach to determine QP energies is the GW approximation. Although the GW method has been widely applied to many typical semiconductors and insulators, its application to more complex compounds such as transition metal oxide perovskite… ▽ More The ab initio calculation of quasiparticle (QP) energies is a technically and computationally challenging problem. In condensed matter physics the most widely used approach to determine QP energies is the GW approximation. Although the GW method has been widely applied to many typical semiconductors and insulators, its application to more complex compounds such as transition metal oxide perovskites has been comparatively rare, and its proper use is not well established from a technical point of view. In this work, we have applied the single-shot G0W0 method to a representative set of transition metal oxide perovskites including 3d (SrTiO3, LaScO3, SrMnO3, LaTiO3, LaVO3, LaCrO3, LaMnO3, and LaFeO3), 4d (SrZrO3, SrTcO3, and Ca2RuO4) and 5d (SrHfO3, KTaO3 and NaOsO3) compounds with different electronic configurations, magnetic orderings, structural characteristics and bandgaps ranging from 0.1 to 6.1 eV. We discuss the proper procedure to obtain well converged QP energies and accurate bandgaps within single-shot G0W0 by comparing the conventional approach based on an incremental variation of a specific set of parameters (number of bands, energy cutoff for the plane-wave expansion and number of k-points and the basis-set extrapolation scheme [Phys. Rev. B 90, 075125 (2014)]. In addition, we have inspected the difference between the adoption of norm-conserving and ultrasoft potentials in GW calculations. A minimal statistical analysis indicates that the correlation of the GW data with the DFT gap is more robust than the correlation with the experimental gaps; moreover we identify the static dielectric constant as alternative useful parameter for the approximation of GW gap in high-throughput automatic procedures. Finally, we compute the QP band structure and spectra within the random phase approximation and compare the results with available experimental data. △ Less

Submitted 26 January, 2018; v1 submitted 2 November, 2017; originally announced November 2017.

Comments: Physical Review Materials, accepted

Journal ref: Phys. Rev. Materials 2, 024601 (2018)