Beyond time-dependent exact-exchange: the need for long-range correlation

In the description of the interaction between electrons beyond the classical Hartree picture, bare exchange often yields a leading contribution. Here we discuss its effect on optical spectra of solids, comparing three different frameworks: time-dependent Hartree-Fock, a recently introduced combined density-functional and Green's functions approach applied to the bare exchange self-energy, and time-dependent exact-exchange within time-dependent density-functional theory (TD-EXX). We show that these three approximations give rise to identical excitonic effects in solids; these effects are drastically overestimated for semiconductors. They are partially compensated by the usual overestimation of the quasiparticle band gap within Hartree-Fock. The physics that lacks in these approaches can be formulated as screening. We show that the introduction of screening in TD-EXX indeed leads to a formulation that is equivalent to previously proposed functionals derived from Many-Body Perturbation Theory. It can be simulated by reducing the long-range part of the Coulomb interaction: this produces absorption spectra of semiconductors in good agreement with experiment.Comment: 12 pages, 3 figures, 1 tabl

Effects of quantum confinement on excited state properties of SrTiO3_3 from ab initio many-body perturbation theory

The Ruddlesden-Popper (RP) homologous series Sr$_{n+1}$Ti$_{n}$O$_{3n+1}$ provides a useful template for the study and control of the effects of dimensionality and quantum confinement on the excited state properties of the complex oxide SrTiO$_3$. We use ab initio many-body perturbation theory within the $GW$ approximation and the Bethe-Salpeter equation approach to calculate quasiparticle energies and absorption spectrum of Sr$_{n+1}$Ti$_{n}$O$_{3n+1}$ for $n=1-5$ and $\infty$. Our computed direct and indirect optical gaps are in excellent agreement with spectroscopic measurements. The calculated optical spectra reproduce the main experimental features and reveal excitonic structure near the gap edge. We find that electron-hole interactions are important across the series, leading to significant exciton binding energies that increase for small $n$ and reach a value of 330~meV for $n=1$, a trend attributed to increased quantum confinement. We find that the lowest-energy singlet exciton of Sr$_2$TiO$_4$ ($n=1$) localizes in the 2D plane defined by the TiO$_2$ layer, and explain the origin of its localization

A systematic benchmark of the ab initio Bethe-Salpeter equation approach for low-lying optical excitations of small organic molecules

The predictive power of the ab initio Bethe-Salpeter equation (BSE) approach, rigorously based on many-body Green's function theory but incorporating information from density functional theory, has already been demonstrated for the optical gaps and spectra of solid-state systems. Interest in photoactive hybrid organic/inorganic systems has recently increased, and so has the use of the BSE for computing neutral excitations of organic molecules. However, no systematic benchmarks of the BSE for neutral electronic excitations of organic molecules exist. Here, we study the performance of the BSE for the 28 small molecules in Thiel's widely-used time-dependent density functional theory benchmark set [M. Schreiber et al. J. Chem. Phys. 128, 134110 (2008)]. We observe that the BSE produces results that depend critically on the mean-field starting point employed in the perturbative approach. We find that this starting point dependence is mainly introduced through the quasiparticle energies obtained at the intermediate GW step, and that with a judicious choice of starting mean-field, singlet excitation energies obtained from BSE are in excellent quantitative agreement with higher-level wavefunction methods. The quality of the triplet excitations is slightly less satisfactory

Point defect modeling in materials: coupling ab initio and elasticity approaches

International audienceModeling point defects at an atomic scale requires careful treatment of the long-range atomic relaxations. This elastic field can strongly affect point defect properties calculated in atomistic simulations because of the finite size of the system under study. This is an important restriction for ab initio methods which are limited to a few hundred atoms. We propose an original approach coupling ab initio calculations and linear elasticity theory to obtain the properties of an isolated point defect for reduced supercell sizes. The reliability and benefit of our approach are demonstrated for three problematic cases: the self-interstitial in zirconium, clusters of self-interstitials in iron, and the neutral vacancy in silicon

GWGW density matrix to estimate self-consistent GWGW total energy in solids

The $GW$ approximation is a well-established method for calculating ionization potentials and electron affinities in solids and molecules. For numerous years, obtaining self-consistent $GW$ total energies in solids has been a challenging objective that is not accomplished yet. However, it was shown recently that the linearized $GW$ density matrix permits a reliable prediction of the self-consistent $GW$ total energy for molecules [F. Bruneval et. al. J. Chem. Theory Comput. 17, 2126 (2021)] for which self-consistent $GW$ energies are available. Here we implement, test, and benchmark the linearized $GW$ density matrix for several solids. We focus on the total energy, lattice constant, and bulk modulus obtained from the $GW$ density matrix and compare our findings to more traditional results obtained within the random phase approximation (RPA). We conclude on the improved stability of the total energy obtained from the linearized $GW$ density matrix with respect to the mean-field starting point. We bring compelling clues that the RPA and the $GW$ density matrix total energies are certainly close to the self-consistent $GW$ total energy in solids if we use hybrid functionals with enriched exchange as a starting point.Comment: 16 pages, 8 figure

Effects of polarization on the band-structure of delafossite transparent conductive oxides

We use hybrid functionals and restricted self-consistent GW, state-of-the-art theoretical approaches for quasiparticle band structures, to study the electronic states of delafossite Cu(Al,In)O$_2$, the first p-type and bipolar transparent conductive oxides. We show that self-consistent GW gives remarkably wider band gaps than all the other approaches used so far. Accounting for polaronic effects in the GW scheme we recover a very nice agreement with experiments. Furthermore, the modifications with respect to the Kohn-Sham bands are strongly k-dependent, which makes questionable the common practice of using a scissor operator. Finally, our results support the view that the low energy structures found in optical experiments, and initially attributed to an indirect transition, are due to intrinsic defects in the samples.Comment: publishe

Vascular Smooth Muscle Cell Plasticity and Autophagy in Dissecting Aortic Aneurysms.

Objective- Recent studies suggested the occurrence of phenotypic switching of vascular smooth muscle cells (VSMCs) during the development of aortic aneurysm (AA). However, lineage-tracing studies are still lacking, and the behavior of VSMCs during the formation of dissecting AA is poorly understood. Approach and Results- We used multicolor lineage tracing of VSMCs to track their fate after injury in murine models of Ang II (angiotensin II)-induced dissecting AA. We also addressed the direct impact of autophagy on the response of VSMCs to AA dissection. Finally, we studied the relevance of these processes to human AAs. Here, we show that a subset of medial VSMCs undergoes clonal expansion and that VSMC outgrowths are observed in the adventitia and borders of the false channel during Ang II-induced development of dissecting AA. The clonally expanded VSMCs undergo phenotypic switching with downregulation of VSMC differentiation markers and upregulation of phagocytic markers, indicative of functional changes. In particular, autophagy and endoplasmic reticulum stress responses are activated in the injured VSMCs. Loss of autophagy in VSMCs through deletion of autophagy protein 5 gene ( Atg5) increases the susceptibility of VSMCs to death, enhances endoplasmic reticulum stress activation, and promotes IRE (inositol-requiring enzyme) 1α-dependent VSMC inflammation. These alterations culminate in increased severity of aortic disease and higher incidence of fatal AA dissection in mice with VSMC-restricted deletion of Atg5. We also report increased expression of autophagy and endoplasmic reticulum stress markers in VSMCs of human dissecting AAs. Conclusions- VSMCs undergo clonal expansion and phenotypic switching in Ang II-induced dissecting AAs in mice. We also identify a critical role for autophagy in regulating VSMC death and endoplasmic reticulum stress-dependent inflammation with important consequences for aortic wall homeostasis and repair

A fast and reliable method for the calculation of band structure of solids with hybrid functionals

A simple approximation within the framework of the hybrid methods for the calculation of the electronic structure of solids is presented. By considering only the diagonal elements of the perturbation operator (Hartree-Fock exchange minus semilocal exchange) calculated in the basis of the semilocal orbitals, the computational time is drastically reduced, while keeping the accuracy of the results obtained with hybrid functionals when applied without any approximations

ABINIT: Overview and focus on selected capabilities

Paper published as part of the special topic on Electronic Structure SoftwareABINIT is probably the first electronic-structure package to have been released under an open-source license about 20 years ago. It implements density functional theory, density-functional perturbation theory (DFPT), many-body perturbation theory (GW approximation and Bethe–Salpeter equation), and more specific or advanced formalisms, such as dynamical mean-field theory (DMFT) and the “temperaturedependent effective potential” approach for anharmonic effects. Relying on planewaves for the representation of wavefunctions, density, and other space-dependent quantities, with pseudopotentials or projector-augmented waves (PAWs), it is well suited for the study of periodic materials, although nanostructures and molecules can be treated with the supercell technique. The present article starts with a brief description of the project, a summary of the theories upon which ABINIT relies, and a list of the associated capabilities. It then focuses on selected capabilities that might not be present in the majority of electronic structure packages either among planewave codes or, in general, treatment of strongly correlated materials using DMFT; materials under finite electric fields; properties at nuclei (electric field gradient, Mössbauer shifts, and orbital magnetization); positron annihilation; Raman intensities and electro-optic effect; and DFPT calculations of response to strain perturbation (elastic constants and piezoelectricity), spatial dispersion (flexoelectricity), electronic mobility, temperature dependence of the gap, and spin-magnetic-field perturbation. The ABINIT DFPT implementation is very general, including systems with van der Waals interaction or with noncollinear magnetism. Community projects are also described: generation of pseudopotential and PAW datasets, high-throughput calculations (databases of phonon band structure, second-harmonic generation, and GW computations of bandgaps), and the library LIBPAW. ABINIT has strong links with many other software projects that are briefly mentioned.This work (A.H.R.) was supported by the DMREF-NSF Grant No. 1434897, National Science Foundation OAC-1740111, and U.S. Department of Energy DE-SC0016176 and DE-SC0019491 projects. N.A.P. and M.J.V. gratefully acknowledge funding from the Belgian Fonds National de la Recherche Scientifique (FNRS) under Grant No. PDR T.1077.15-1/7. M.J.V. also acknowledges a sabbatical “OUT” grant at ICN2 Barcelona as well as ULiège and the Communauté Française de Belgique (Grant No. ARC AIMED G.A. 15/19-09). X.G. and M.J.V. acknowledge funding from the FNRS under Grant No. T.0103.19-ALPS. X.G. and G.-M. R. acknowledge support from the Communauté française de Belgique through the SURFASCOPE Project (No. ARC 19/24-057). X.G. acknowledges the hospitality of the CEA DAM-DIF during the year 2017. G.H. acknowledges support from the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Materials Sciences and Engineering Division under Contract No. DE-AC02-05-CH11231 (Materials Project Program No. KC23MP). The Belgian authors acknowledge computational resources from supercomputing facilities of the University of Liège, the Consortium des Equipements de Calcul Intensif (Grant No. FRS-FNRS G.A. 2.5020.11), and Zenobe/CENAERO funded by the Walloon Region under Grant No. G.A. 1117545. M.C. and O.G. acknowledge support from the Fonds de Recherche du Québec Nature et Technologie (FRQ-NT), Canada, and the Natural Sciences and Engineering Research Council of Canada (NSERC) under Grant No. RGPIN-2016-06666. The implementation of the libpaw library (M.T., T.R., and D.C.) was supported by the ANR NEWCASTLE project (Grant No. ANR-2010-COSI-005-01) of the French National Research Agency. M.R. and M.S. acknowledge funding from Ministerio de Economia, Industria y Competitividad (MINECO-Spain) (Grants Nos. MAT2016-77100-C2-2-P and SEV-2015-0496) and Generalitat de Catalunya (Grant No. 2017 SGR1506). This work has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 Research and Innovation program (Grant Agreement No. 724529). P.G. acknowledges support from FNRS Belgium through PDR (Grant No. HiT4FiT), ULiège and the Communauté française de Belgique through the ARC project AIMED, the EU and FNRS through M.ERA.NET project SIOX, and the European Funds for Regional Developments (FEDER) and the Walloon Region in the framework of the operational program “Wallonie-2020.EU” through the project Multifunctional thin films/LoCoTED. The Flatiron Institute is a division of the Simons Foundation. A large part of the data presented in this paper is available directly from the Abinit Web page www.abinit.org. Any other data not appearing in this web page can be provided by the corresponding author upon reasonable request.Peer reviewe