Search

Electron-phonon coupling (EPC) is key for understanding many properties of materials such as superconductivity and electric resistivity. Although first principles density-functional-theory (DFT) based EPC calculations are used widely, their efficacy is limited by the accuracy and efficiency of the underlying exchange-correlation functionals. These limitations become exacerbated in complex $d$- and $f$-electron materials, where beyond-DFT approaches and empirical corrections, such as the Hubbard $U$, are commonly invoked. Here, using the examples of CoO and NiO, we show how the efficient r2scan density functional correctly captures strong EPC effects in transition-metal oxides without requiring the introduction of empirical parameters. We also demonstrate the ability of r2scan to accurately model phonon-mediated superconducting properties of the main group compounds (e.g., MgB$_2$), with improved electronic bands and phonon dispersions over those of traditional density functionals. Our study provides a pathway for extending the scope of accurate first principles modeling of electron-phonon interactions to encompass complex $d$-electron materials., Comment: 9 pages, 3 figures, 1 table

The GW approximation represents the state-of-the-art ab-initio method for computing excited-state properties. Its execution requires control over a larger number of (often interdependent) parameters, and therefore its application in high-throughput studies is hindered by the intricate and time-consuming convergence process across a multi-dimensional parameter space. To address these challenges, here we develop a fully-automated open-source workflow for G$_0$W$_0$ calculations within the AiiDA-VASP plugin architecture. The workflow is based on an efficient estimation of the errors on the quasi-particle (QP) energies due to basis-set truncation and the pseudo-potential norm violation, which allows a reduction of the dimensionality of the parameter space and avoids the need for multi-dimensional convergence searches. Protocol validation is conducted through a systematic comparison against established experimental and state-of-the-art GW data. To demonstrate the effectiveness of the approach, we construct a database of QP energies for a diverse dataset of over 320 bulk structures. The openly accessible workflow and resulting dataset can serve as a valuable resource and reference for conducting accurate data-driven research., Comment: 11 pages, 6 figures, 2 tables

We investigate the relationship between the K-edge fine structure of isolated single-wall carbon nanotubes (SWCNTs) and the Van Hove singularities (VHSs) in the conduction band density of states. To this end, we model X-ray absorption spectra of SWCNTs using the final-state approximation and the Bethe-Salpeter equation (BSE) method. Both methods can reproduce the experimental fine structure, where the BSE results improve on peak positions and amplitude rations compared to the final-state approximation. When the fine structure in the modeled spectra is related to the VHSs significant differences are found. We suggest that these differences arise due to modifications of the core exciton wavefunctions induced by the confinement along the circumference. Additionally, we analyze the character of core excitons in SWCNTs and find that the first bright excitons are Frenkel excitons while higher-lying excitons are charge resonance states. Finally, we suggest that the qualitative picture based on VHSs in the density of states holds when there is a large energy gap between successive VHSs., Comment: 12 pages, 8 figures

Constructing a self-consistent first-principles framework that accurately predicts the properties of electron transfer reactions through finite-temperature molecular dynamics simulations is a dream of theoretical electrochemists and physical chemists. Yet, predicting even the absolute standard hydrogen electrode potential, the most fundamental reference for electrode potentials, proves to be extremely challenging. Here, we show that a hybrid functional incorporating 25 % exact exchange enables quantitative predictions when statistically accurate phase-space sampling is achieved via thermodynamic integrations and thermodynamic perturbation theory calculations, utilizing machine-learned force fields and $\Delta$-machine learning models. The application to seven redox couples, including molecules and transition metal ions, demonstrates that the hybrid functional can predict redox potentials across a wide range of potentials with an average error of 80 mV., Comment: 10 pages, 4 figures

The magnetite/water interface is commonly found in nature and plays a crucial role in various technological applications. However, our understanding of its structural and dynamical properties at the molecular scale remains still limited. In this study, we develop an efficient Behler-Parrinello neural network potential (NNP) for the magnetite/water system, paying particular attention to the accurate generation of reference data with density functional theory. Using this NNP, we performed extensive molecular dynamics simulations of the magnetite (001) surface across a wide range of water coverages, from the single molecule to bulk water. Our simulations revealed several new ground states of low coverage water on the Subsurface Cation Vacancy (SCV) model and yielded a density profile of water at the surface that exhibits marked layering. By calculating mean square displacements, we obtained quantitative information on the diffusion of water molecules on the SCV for different coverages, revealing significant anisotropy. Additionally, our simulations provided qualitative insights into the dissociation mechanisms of water molecules at the surface.

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%, worse than LODE but still very good. For real materials, we perform tests for liquid NaCl, rock salt NaCl, and solid zirconia. For NaCl, the present descriptors improve on short-range density descriptors, reducing errors by a factor of two to three and coming close to message-passing networks. However, for solid zirconia, no improvements are observed with the present approach, while message-passing networks reduce the error by almost a factor of two to three. Possible shortcomings of the present model are briefly discussed.

Macroscopic properties of materials stem from fundamental atomic-scale details, yet for insulators, resolving surface structures remains a challenge. The basal (0001) plane of ${\alpha}$-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{\deg}$ 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}$.

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., Comment: 9 pages, 4 figures

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., Comment: 13 pages, 6 figures

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.

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 $\beta$- to $\lambda$-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 $\beta$-$\lambda$ 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 $\beta$-$\lambda$ transition, but also presents useful strategies and methods for tackling other complex structural phase transitions., Comment: 26 pages,23 figures (including Supporting Information)

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 $\Delta$-machine learning to determine the redox potentials for computationally expensive hybrid functionals. The combination of these approaches allows one to obtain statistically accurate results.

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.

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)., Comment: Main text: 23 pages, 4 figures. Supplementary: 68 pages. Nature Review Physics 2023

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., Comment: 13 pages, 7 figures

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.

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., Comment: 10 pages, 6 figures. To be published in Advanced Quantum Technologies. Editor: Luigi Ranalli

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., Comment: 16 pages, 2 figures

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 $\Delta$-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., Comment: 6 pages, 2 figures

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.

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\"ohlich 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., Comment: 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

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.

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 $\Delta$-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.

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.

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$.

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\"uter equation. We investigate the accuracy of the QPA on Kohn-Sham band gaps and dielectric constants, and comment on the issue of self-consistency., Comment: postprint

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 {\AA} and 0.13 {\AA}, respectively. The frequency of the highest optical mode at the $\Gamma$ 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 $\Gamma$ 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 {\AA} 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., Comment: 7 pages, 5 figures

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., Comment: 6 pages, 2 figures

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 $\alpha$-$\beta$ 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., Comment: 10 pages, 7 figures

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.

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., Comment: 9 pages, 9 figures

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.

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.

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., Comment: 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)

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., Comment: 10 pages, 18 figures, 16 Equations; v2: fixed typo in metadata abstract, changed bibliography style; v3: added reference (4), added Acknowledgement

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$., Comment: 19 pages, 11 figures

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., Comment: 15 pages, 7 figures and 2 tables

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.

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., Comment: 5 pages, 4 figures

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., Comment: 11 pages, 5 figures

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., Comment: 36 pages, 12 figures

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.

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., Comment: accepted in Phys. Rev. B

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., Comment: 9 pages, 5 figures

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., Comment: 5 pages, 3 figures

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., Comment: 7 pages, 7 figures

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., Comment: 13 pages, 8 figures

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., Comment: 15 pages, 8 figures

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.

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., Comment: 11 pages, 6 figures