Semilocal Density Functional Theory Research Articles

Eliminating Delocalization Error to Improve Heterogeneous Catalysis Predictions with Molecular DFT + U

Approximate semilocal density functional theory (DFT) is known to underestimate surface formation energies yet paradoxically overbind adsorbates on catalytic transition-metal oxide surfaces due to delocalization error. The low-cost DFT + U approach only improves surface formation energies for early transition-metal oxides or adsorption energies for late transition-metal oxides. In this work, we demonstrate that this inefficacy arises due to the conventional usage of metal-centered atomic orbitals as projectors within DFT + U. We analyze electron density rearrangement during surface formation and O atom adsorption on rutile transition-metal oxides to highlight that a standard DFT + U correction fails to tune properties when the corresponding density rearrangement is highly delocalized across both metal and oxygen sites. To improve both surface properties simultaneously while retaining the simplicity of a single-site DFT + U correction, we systematically construct multi-atom-centered molecular-orbital-like projectors for DFT + U. We demonstrate this molecular DFT + U approach for tuning adsorption energies and surface formation energies of minimal two-dimensional models of representative early (i.e., TiO2) and late (i.e., PtO2) transition-metal oxides. Molecular DFT + U simultaneously corrects adsorption energies and surface formation energies of multilayer models of rutile TiO2(110) and PtO2(110) to resolve the paradoxical description of surface stability and surface reactivity of semilocal DFT.

Elastic stability of Ga2O3 : Addressing the β to α phase transition from first principles

The elastic and structural properties of $\ensuremath{\beta}\text{\ensuremath{-}}{\mathrm{Ga}}_{2}{\mathrm{O}}_{3}$ and $\ensuremath{\alpha}\text{\ensuremath{-}}{\mathrm{Ga}}_{2}{\mathrm{O}}_{3}$ are investigated from first principles. The full elastic tensors and elastic moduli of both phases at $0\phantom{\rule{0.16em}{0ex}}\mathrm{K}$ are computed in the framework of semilocal density-functional theory. We determine mechanical instabilities of $\ensuremath{\beta}\text{\ensuremath{-}}{\mathrm{Ga}}_{2}{\mathrm{O}}_{3}$ by evaluating the full stiffness tensor under load for a range of hydrostatic pressure values. While a phase transition from the $\ensuremath{\beta}$ to $\ensuremath{\alpha}$ phase is found to be energetically favored at $2.6\phantom{\rule{0.16em}{0ex}}\mathrm{G}\mathrm{Pa}$, we show that the $\ensuremath{\beta}$ phase is mechanically unstable only for much higher pressures $(>30\phantom{\rule{0.16em}{0ex}}\mathrm{G}\mathrm{Pa})$, which agrees well with experimental results. Our employed approach is based on the Born stability criterion, is independent of crystal symmetry, and thus can be readily applied to different materials.

Beyond-dipole van der Waals contributions within the many-body dispersion framework

By introducing a suitable range-separation of the Coulomb coupling in analogy to Ambrosetti et al (2014 J. Chem. Phys. 140 18A508), here we extend the many-body dispersion approach to include beyond-dipole van der Waals (vdW) interactions at a full many-body level, in combination with semi-local density functional theory. A reciprocal-space implementation is further introduced in order to efficiently treat periodic systems. Consistent reliability is found from molecular dimers to large supramolecular complexes and two-dimensional systems. The large weight of both many-body effects and multipolar terms illustrates how a correct description of vdW forces in large-scale systems requires full account of both contributions, beyond standard pairwise dipolar approaches.

Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory-Based Ab Initio Molecular Dynamics II: Extensions to the Isobaric-Isoenthalpic and Isobaric-Isothermal Ensembles.

In the previous paper of this series [Ko, H.-Y. et al. J. Chem. Theory Comput. 2020, 16, 3757-3785], we presented a theoretical and algorithmic framework based on a localized representation of the occupied space that exploits the inherent sparsity in the real-space evaluation of the exact exchange (EXX) interaction in finite-gap systems. This was accompanied by a detailed description of exx, a massively parallel hybrid message-passing interface MPI/OpenMP implementation of this approach in Quantum ESPRESSO (QE) that enables linear scaling hybrid density functional theory (DFT)-based ab initio molecular dynamics (AIMD) in the microcanonical/canonical (NVE/NVT) ensembles of condensed-phase systems containing 500-1000 atoms (in fixed orthorhombic cells) with a wall time cost comparable to semi-local DFT. In this work, we extend the current capabilities of exx to enable hybrid DFT-based AIMD simulations of large-scale condensed-phase systems with general and fluctuating cells in the isobaric-isoenthalpic/isobaric-isothermal (NpH/NpT) ensembles. The theoretical extensions to this approach include an analytical derivation of the EXX contribution to the stress tensor for systems in general simulation cells with a computational complexity that scales linearly with system size. The corresponding algorithmic extensions to exx include optimized routines that (i) handle both static and fluctuating simulation cells with non-orthogonal lattice symmetries, (ii) solve Poisson's equation in general/non-orthogonal cells via an automated selection of the auxiliary grid directions in the Natan-Kronik representation of the discrete Laplacian operator, and (iii) evaluate the EXX contribution to the stress tensor. Using this approach, we perform a case study on a variety of condensed-phase systems (including liquid water, a benzene molecular crystal polymorph, and semi-conducting crystalline silicon) and demonstrate that the EXX contributions to the energy and stress tensor simultaneously converge with an appropriate choice of exx parameters. This is followed by a critical assessment of the computational performance of the extended exx module across several different high-performance computing architectures via case studies on (i) the computational complexity due to lattice symmetry during NpT simulations of three different ice polymorphs (i.e., ice Ih, II, and III) and (ii) the strong/weak parallel scaling during large-scale NpT simulations of liquid water. We demonstrate that the robust and highly scalable implementation of this approach in the extended exx module is capable of evaluating the EXX contribution to the stress tensor with negligible cost (<1%) as well as all other EXX-related quantities needed during NpT simulations of liquid water (with a very tight 150 Ry planewave cutoff) in ≈5.2 s ((H2O)128) and ≈6.8 s ((H2O)256) per AIMD step. As such, the extended exx module presented in this work brings us another step closer to routinely performing hybrid DFT-based AIMD simulations of sufficient duration for large-scale condensed-phase systems across a wide range of thermodynamic conditions.

Assessing the G0W0Γ0(1) Approach: Beyond G0W0 with Hedin's Full Second-Order Self-Energy Contribution.

We present and benchmark a self-energy approach for quasiparticle energy calculations that goes beyond Hedin's GW approximation by adding the full second-order self-energy (FSOS-W) contribution. The FSOS-W diagram involves two screened Coulomb interaction (W) lines, and adding the FSOS-W to the GW self-energy can be interpreted as first-order vertex correction to GW (GWΓ(1)). Our FSOS-W implementation is based on the resolution-of-identity technique and exhibits better than O(N5) scaling with system size for small- to medium-sized molecules. We then present one-shot GWΓ(1) (G0W0Γ0(1)) benchmarks for the GW100 test set and a set of 24 acceptor molecules. For semilocal or hybrid density functional theory starting points, G0W0Γ0(1) systematically outperforms G0W0 for the first vertical ionization potentials and electron affinities of both test sets. Finally, we demonstrate that a static FSOS-W self-energy significantly underestimates the quasiparticle energies.

Machine-learning semilocal density functional theory for many-body lattice models at zero and finite temperature

We introduce a machine-learning density-functional-theory formalism for the spinless Hubbard model in one dimension at both zero and finite temperature. In the zero-temperature case this establishes a one-to-one relation between the site occupation and the total energy, which is then minimised at the ground-state occupation. In contrast, at finite temperature the same relation is defined between the Helmholtz free energy and the equilibrium site occupation. Most importantly, both functionals are semi-local, so that they are independent from the size of the system under investigation and can be constructed over exact data for small systems. These 'exact' functionals are numerically defined by neural networks. We also define additional neural networks for finite-temperature thermodynamical quantities, such as the entropy and heat capacity. These can be either a functional of the ground-state site occupation or of the finite-temperature equilibrium site occupation. In the first case their equilibrium value does not correspond to an extremal point of the functional, while it does in the second case. Our work gives us access to finite-temperature properties of many-body systems in the thermodynamic limit.

Efficient Band Structure Calculation of Two-Dimensional Materials from Semilocal Density Functionals.

The experimental and theoretical realization of two-dimensional (2D) materials is of utmost importance in semiconducting applications. Computational modeling of these systems with satisfactory accuracy and computational efficiency is only feasible with semilocal density functional theory methods. In the search for the most useful method in predicting the band gap of 2D materials, we assess the accuracy of recently developed semilocal exchange–correlation (XC) energy functionals and potentials. Though the explicit forms of exchange–correlation (XC) potentials are very effective against XC energy functionals for the band gap of bulk solids, their performance needs to be investigated for 2D materials. In particular, the LMBJ [J. Chem. Theory Comput.2020, 16, 265432097004] and GLLB-SC [Phys. Rev. B82, 2010, 115106] potentials are considered for their dominance in bulk band gap calculation. The performance of recently developed meta generalized gradient approximations, like TASK [Phys. Rev. Res.1, 2019, 033082] and MGGAC [Phys. Rev. B. 100, 2019, 155140], is also assessed. We find that the LMBJ potential constructed for 2D materials is not as successful as its parent functional, i.e., MBJ [Phys. Rev. Lett.102, 2009, 22640119658882] in bulk solids. Due to a contribution from the derivative discontinuity, the band gaps obtained with GLLB-SC are in a certain number of cases, albeit not systematically, larger than those obtained with other methods, which leads to better agreement with the quasi-particle band gap obtained from the GW method. The band gaps obtained with TASK and MGGAC can also be quite accurate.

Revisiting Understanding of Electrochemical CO2 Reduction on Cu(111): Competing Proton-Coupled Electron Transfer Reaction Mechanisms Revealed by Embedded Correlated Wavefunction Theory.

Copper (Cu) electrodes, as the most efficacious of CO2 reduction reaction (CO2RR) electrocatalysts, serve as prototypes for determining and validating reaction mechanisms associated with electrochemical CO2 reduction to hydrocarbons. As in situ electrochemical mechanism determination by experiments is still out of reach, such mechanistic analysis typically is conducted using density functional theory (DFT). The semilocal exchange-correlation (XC) approximations most often used to model such catalysis unfortunately engender a basic error: predicting the wrong adsorption site for CO (a key CO2RR intermediate) on the most ubiquitous facet of Cu, namely, Cu(111). This longstanding inconsistency casts lingering doubt on previous DFT predictions of the attendant CO2RR kinetics. Here, we apply embedded correlated wavefunction (ECW) theory, which corrects XC functional error, to study the CO2RR on Cu(111) via both surface hydride (*H) transfer and proton-coupled electron transfer (PCET). We predict that adsorbed CO (*CO) reduces almost equally to two intermediates, namely, hydroxymethylidyne (*COH) and formyl (*CHO) at -0.9 V vs the RHE. In contrast, semilocal DFT approximations predict a strong preference for *COH. With increasing applied potential, the dominance of *COH (formed via potential-independent surface *H transfer) diminishes, switching to the competitive formation of both *CHO and *COH (both formed via potential-dependent PCET). Our results also demonstrate the importance of including explicitly modeled solvent molecules in predicting electron-transfer barriers and reveal the pitfalls of overreliance on simple surface *H transfer models of reduction reactions.

Photoelectron spectra of early 3d-transition metal dioxide molecular anions from GW calculations.

Photoelectron spectra of early 3d-transition metal dioxide anions, ScO2 -, TiO2 -, VO2 -, CrO2 -, and MnO2 -, are calculated using semilocal and hybrid density functional theory (DFT) and many-body perturbation theory within the GW approximation using one-shot perturbative and eigenvalue self-consistent formalisms. Different levels of theory are compared with each other and with available photoelectron spectra. We show that one-shot GW with a PBE0 starting point (G0W0@PBE0) consistently provides very good agreement for all experimentally measured binding energies (within 0.1 eV-0.2 eV or less). We attribute this to the success of PBE0 in mitigating self-interaction error and providing good quasiparticle wave functions, which renders a first-order perturbative GW correction effective. One-shot GW calculations with a Perdew-Burke-Ernzerhof (PBE) starting point do poorly in predicting electron removal energies by underbinding orbitals with typical errors near 1.5 eV. A higher exact exchange amount of 50% in the DFT starting point of one-shot GW does not provide very good agreement with experiment by overbinding orbitals with typical errors near 0.5 eV. While not as accurate as G0W0@PBE0, the G-only eigenvalue self-consistent GW scheme with W fixed to the PBE level provides a reasonably predictive level of theory (typical errors near 0.3 eV) to describe photoelectron spectra of these 3d-transition metal dioxide anions. Adding eigenvalue self-consistency also in W, on the other hand, worsens the agreement with experiment overall. Our findings on the performance of various GW methods are discussed in the context of our previous studies on other transition metal oxide molecular systems.

R2SCAN-3c: A "Swiss army knife" composite electronic-structure method.

The recently proposed r2SCAN meta-generalized-gradient approximation (mGGA) of Furness and co-workers is used to construct an efficient composite electronic-structure method termed r2SCAN-3c. To this end, the unaltered r2SCAN functional is combined with a tailor-made triple-ζ Gaussian atomic orbital basis set as well as with refitted D4 and geometrical counter-poise corrections for London-dispersion and basis set superposition error. The performance of the new method is evaluated for the GMTKN55 database covering large parts of chemical space with about 1500 data points, as well as additional benchmarks for non-covalent interactions, organometallic reactions, and lattice energies of organic molecules and ices, as well as for the adsorption on polar salt and non-polar coinage-metal surfaces. These comprehensive tests reveal a spectacular performance and robustness of r2SCAN-3c: It by far surpasses its predecessor B97-3c at only twice the cost and provides one of the best results of all semi-local density-functional theory (DFT)/QZ methods ever tested for the GMTKN55 database at one-tenth of the cost. Specifically, for reaction and conformational energies as well as non-covalent interactions, it outperforms prominent hybrid-DFT/QZ approaches at two to three orders of magnitude lower cost. Perhaps, the most relevant remaining issue of r2SCAN-3c is self-interaction error (SIE), owing to its mGGA nature. However, SIE is slightly reduced compared to other (m)GGAs, as is demonstrated in two examples. After all, this remarkably efficient and robust method is chosen as our new group default, replacing previous composite DFT and partially even expensive high-level methods in most standard applications for systems with up to several hundreds of atoms.

Random-Phase Approximation in Many-Body Noncovalent Systems: Methane in a Dodecahedral Water Cage.

The many-body expansion (MBE) of energies of molecular clusters or solids offers a way to detect and analyze errors of theoretical methods that could go unnoticed if only the total energy of the system was considered. In this regard, the interaction between the methane molecule and its enclosing dodecahedral water cage, CH4···(H2O)20, is a stringent test for approximate methods, including density functional theory (DFT) approximations. Hybrid and semilocal DFT approximations behave erratically for this system, with three- and four-body nonadditive terms having neither the correct sign nor magnitude. Here, we analyze to what extent these qualitative errors in different MBE contributions are conveyed to post-Kohn-Sham random-phase approximation (RPA), which uses approximate Kohn-Sham orbitals as its input. The results reveal a correlation between the quality of the DFT input states and the RPA results. Moreover, the renormalized singles energy (RSE) corrections play a crucial role in all orders of the many-body expansion. For dimers, RSE corrects the RPA underbinding for every tested Kohn-Sham model: generalized-gradient approximation (GGA), meta-GGA, (meta-)GGA hybrids, as well as the optimized effective potential at the correlated level. Remarkably, the inclusion of singles in RPA can also correct the wrong signs of three- and four-body nonadditive energies as well as mitigate the excessive higher-order contributions to the many-body expansion. The RPA errors are dominated by the contributions of compact clusters. As a workable method for large systems, we propose to replace those compact contributions with CCSD(T) energies and to sum up the remaining many-body contributions up to infinity with supermolecular or periodic RPA. As a demonstration of this approach, we show that for RPA(PBE0)+RSE it suffices to apply CCSD(T) to dimers and 30 compact, hydrogen-bonded trimers to get the methane-water cage interaction energy to within 1.6% of the reference value.

New insights into the 1D carbon chain through the RPA.

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.

Critical assessment of Co-Cu phase diagram from first-principles calculations

Using first-principles alloy theory, we perform a systematic study of the Co-Cu phase diagram. Calculations are carried out for ferromagnetic and paramagnetic ${\mathrm{Co}}_{1\ensuremath{-}x}{\mathrm{Cu}}_{x}$ solid solutions with face-centered-cubic (fcc) crystal structure. We find that the equilibrium volumes and magnetic states are crucial for a quantitative description of the thermodynamics of the Co-Cu system at temperatures up to 1400 K. In particular, the paramagnetic state of Cu-rich alloys with persisting local magnetic moments is shown to be responsible for the solubility of a small amount of Co in fcc Cu whereas the excess entropy in the ferromagnetic Co-rich region critically depends on the adopted lattice parameters. None of the common local or semilocal density functional theory approximations have the necessary accuracy for the lattice parameters when compared to the experimental data. The predicted ab initio Co-Cu phase diagram is in good agreement with the measurements and CALPHAD data, making it possible to gain a deep insight into the various contributions to the Gibbs free energy. The present study provides an atomic-level description of the thermodynamic quantities controlling the limited mutual solubility of Co and Cu and highlights the importance of high-temperature magnetism.

Termination Effects in Aluminosilicate and Aluminogermanate Imogolite Nanotubes: A Density Functional Theory Study

We investigate termination effects in aluminosilicate (AlSi) and aluminogermanate (AlGe) imogolite nanotubes (NTs) by means of semi-local and range-corrected hybrid Density Functional Theory (DFT) simulations. Following screening and identification of the smallest finite model capable of accommodating full relaxation of the NT terminations around an otherwise geometrically and electrostatically unperturbed core region, we quantify and discuss the effects of physical truncation on the structure, relative energy, electrostatics and electronic properties of differently terminated, finite-size models of the NTs. In addition to composition-dependent changes in the valence (VB) and conduction band (CB) edges and resultant band gap (BG), the DFT simulations uncover longitudinal band bending and separation in the finite AlSi and AlGe models. Depending on the given termination of the NTs, such longitudinal effects manifest in conjunction with the radial band separation typical of fully periodic AlSi and AlGe NTs. The strong composition dependence of the longitudinal and radial band bending in AlSi and AlGe NTs suggests different mechanisms for the generation, relaxation and separation of photo-generated holes in AlSi and AlGe NTs, inviting further research in the untapped potential of imogolite compositional and structural flexibility for photo-catalytic applications.

Photophysics and Electronic Structure of Lateral Graphene/MoS2 and Metal/MoS2 Junctions.

Integration of semiconducting transition metal dichalcogenides (TMDs) into functional optoelectronic circuitries requires an understanding of the charge transfer across the interface between the TMD and the contacting material. Here, we use spatially resolved photocurrent microscopy to demonstrate electronic uniformity at the epitaxial graphene/molybdenum disulfide (EG/MoS2) interface. A 10× larger photocurrent is extracted at the EG/MoS2 interface when compared to the metal (Ti/Au)/MoS2 interface. This is supported by semi-local density functional theory (DFT), which predicts the Schottky barrier at the EG/MoS2 interface to be ∼2× lower than that at Ti/MoS2. We provide a direct visualization of a 2D material Schottky barrier through combination of angle-resolved photoemission spectroscopy with spatial resolution selected to be ∼300 nm (nano-ARPES) and DFT calculations. A bending of ∼500 meV over a length scale of ∼2-3 μm in the valence band maximum of MoS2 is observed via nano-ARPES. We explicate a correlation between experimental demonstration and theoretical predictions of barriers at graphene/TMD interfaces. Spatially resolved photocurrent mapping allows for directly visualizing the uniformity of built-in electric fields at heterostructure interfaces, providing a guide for microscopic engineering of charge transport across heterointerfaces. This simple probe-based technique also speaks directly to the 2D synthesis community to elucidate electronic uniformity at domain boundaries alongside morphological uniformity over large areas.

P‐Type Transparent Quadruple Perovskite Halide Conductors: Fact or Fiction?

AbstractHigh‐performance p‐type transparent conductors (TCs) can be used for a variety of optoelectronic applications, particularly transparent electronics with complementary metal–oxide–semiconductor circuits. Recently, several quadruple perovskite halides, such as Cs4CdSb2Cl12 and Cs4CdBi2Cl12, have been theoretically predicted to be excellent p‐type TCs that surpass all the conventional p‐type TCs and hence attracted a lot of attention. In this work, experiments and hybrid density functional theory (DFT) calculations are combined to evaluate the p‐type dopability in these quadruple perovskite halides. The experimental results reveal that these materials are optically transparent but electrically insulating, which contradicts the previous prediction. The hybrid DFT calculations show that the difficulty in forming p‐type doping in these compounds is because of the easy formation of the compensating n‐type intrinsic defects in the p‐type region, which is resulted from too deep valence band maximum levels. It is further demonstrated that the bandgap underestimation associated with the semilocal DFT and the unintentional chemical boundary overestimation should be the origins of the p‐type conductivity predicted previously. The results highlight the importance of reasonable bandgap description and chemical boundary determination in predicting defect thermodynamics.

Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based Ab Initio Molecular Dynamics. 1. Theory, Algorithm, and Performance.

By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable description of the underlying electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we describe a linear-scaling approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions (MLWFs)) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based ab initio molecular dynamics (AIMD) simulations of large-scale condensed-phase systems of interest at the hybrid DFT level. We focus our theoretical discussion on the integration of this approach into the framework of Car-Parrinello AIMD, and highlight the central role played by the MLWF-product potential (i.e., the solution of Poisson's equation for each corresponding MLWF-product density) in the evaluation of the EXX energy and wave function forces. We then provide a comprehensive description of the exx algorithm implemented in the open-source Quantum ESPRESSO program, which employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize the high-performance computing (HPC) resources available on current- and next-generation supercomputer architectures. This is followed by a critical assessment of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach when AIMD simulations of liquid water are performed in the canonical (NVT) ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT-based AIMD simulations of condensed-phase systems containing 500-1000 atoms (e.g., (H2O)256) with a wall time cost that is comparable to that of semilocal DFT. In doing so, exx takes us one step closer to routinely performing AIMD simulations of complex and large-scale condensed-phase systems for sufficiently long time scales at the hybrid DFT level of theory.

Charge Density and Redox Potential of LiNiO2 Using Ab Initio Diffusion Quantum Monte Carlo

Electronic structure of layered LiNiO2 has been controversial despite numerous theoretical and experimental reports regarding its nature. We investigate the charge densities, lithium intercalation potentials and Li diffusion barrier energies of LixNiO2 (0.0 < x < 1.0) system using a truly ab-initio method, diffusion quantum Monte Carlo (DMC). We compare the charge densities from DMC and density functional theory (DFT) and show that local and semi-local DFT functionals yield spin polarization densities with incorrect sign on the oxygen atoms. SCAN functional and Hubbard-U correction improves the polarization density around Ni and O atoms, resulting in smaller deviations from the DMC densities. DMC accurately captures the p-d hybridization between the Ni-O atoms, yielding accurate lithium intercalation voltages, polarization densities and reaction barriers.

Electron Confinement and Magnetism of (LaTiO3)1/(SrTiO3)5 Heterostructure: A Diffusion Quantum Monte Carlo Study.

We have applied the diffusion quantum Monte Carlo (DMC) method to study the electron confinement and magnetic structure in the (LaTiO3)1/(SrTiO3)5 heterostructure. The DMC results were compared with various density functional theory (DFT) methods, including local density approximation (LDA), generalized gradient approximation (GGA), LDA+U, and GGA+U, as well as the recently proposed strongly constrained appropriately normed (SCAN) and van der Waals-Bayesian error estimation functional (vdW-BEEF). We found that many-body correlations are crucial to accurately describe the localization of the two-dimensional (2D) electron gas around the lanthanum planes. DMC predicts 20% more electron density within the first interfacial titanium layer in (LaTiO3)1/(SrTiO3)5 than LDA+U, suggesting that the degree of confinement of the 2D electron gas in the interfacial region is underestimated with semilocal DFT approximations. DMC yields the ferromagnetic (FM) state as the ground state of (LaTiO3)1/(SrTiO3)5 and the antiferromagnetic (AFM) and nonmagnetic (NM) states that are higher in energy by 37(15) and 238(15) meV per lanthanum atom at the interface, respectively. Most DFT methods yield the FM and NM states within less than 25 meV in energy and could not find the AFM state.

Impact of Approximate DFT Density Delocalization Error on Potential Energy Surfaces in Transition Metal Chemistry.

For approximate density functional theory (DFT) to be useful in catalytic applications of transition metal complexes, modeling strategies must simultaneously address electronic, geometric, and energetic properties of the relevant species. We show that for representative transition metal triatomics (MO2, where M = Cr, Mn, Fe, Co, or Ni) and related diatomics the incorporation of Hartree-Fock (HF) exchange in most cases improves the properties of the Born-Oppenheimer potential energy surface (PES) with respect to accurate experimental or CCSD(T) references. We rationalize this observation by noting reduced delocalization obtained with hybrid functionals (20-40% HF exchange), as evidenced by reduced hybridization of non-bonding orbitals and increases in metal partial charges. Although we show that the optimal exchange fraction is both property and system specific, incorporating HF exchange synergistically improves properties of density, structure, and energetics within a single PES characterized by moderately covalent bonding. The same improvement is not observed in the ordering of MO2 spin states, as good agreement of semi-local DFT spin state ordering is worsened by over-stabilization of higher spin states when HF exchange is added. More work is needed to understand minimal functional forms capable of improving multiple properties with respect to semi-local DFT descriptions of transition metal chemistry.