The GW formalism, a Green’s function many-body perturbation theory, is growing in popularity for the description of the electronic properties of condensed matter systems in solid-state physics, and more recently chemistry. Unfortunately, its application to complex systems of interest in nanosciences, chemistry, or even biology, is hampered by the large associated computing cost, in particular in the case of disordered systems, or systems immersed in an opened environment (a solvent, a molecular medium, an electrode, etc.) The goal of the present PhD thesis is to focus on the development of multiscale techniques, merging high-level many-body treatments of the subsystem of interest, with a simplified but fully ab initio description of the electrostatic and dielectric environment. Such approaches aim to go beyond classical parametrized models, mainly developed in the quantum chemistry community, which are based on a continuum (“polarizable continuum model”) or discrete (QM/MM) description of the environment.To reach such a goal, we adopt a divide-and-conquer fragmentation scheme for the environment, particularly suited to molecular systems. This leads to a block-diagonal non- interacting electron susceptibility, decreasing the algorithmic complexity from quartic to cubic. To reduce the prefactor associated with the inversion of the Dyson equation for the screened Coulomb potential W, we have further developed a compression algorithm for the susceptibility operator. The automatic computation of an extremely compact polarization basis set allows a large reduction of the size of the susceptibility blocks, associated to the fragments in the environment. Such a method enables us to compute the dielectric response of systems made of several hundred thousand atoms, with an excellent accuracy when it comes to reproduce the effect of the environment as a response to an excitation in the immersed subsystem. This approach is presented through the study of fullerene bulk, surface and subsurface crystals.While the GW formalism is dynamical, with a frequency-dependent screened Coulomb potential W, a first study is done adopting a static approximation (low-frequency limit) for the screening properties of the environment. Such an approach follows the traditional semi-empirical models of a polarizable environment. This PhD thesis assesses the validity of such an approximation, which assumes an instantaneous response (adiabatic limit) of the environment to an electronic excitation, thanks to an explicit comparison with a fully dynamical dielectric response of the environment. The study of a surface of fullerenes, as well as a water molecule inside a metallic carbon nanotube, show that a static description of the environment leads to errors on the polarization energy below 10%, provided that the “folding” of the environment is treated in a proper way.The fragment approach is also applied to covalent insulator crystals, and more particularly to hexagonal boron nitride (h-BN). We explain how to compute the energy levels of point defects in h-BN, in the true dilute limit, and we give the asymptotic scaling laws for the renormalization of these energy levels, from the monolayer to a (n)-layer system. This study highlights thus the possibility to apply the fragment approach to covalent insulator systems, a possibility hinging probably on the short range behavior of the susceptibility in these systems.All of these developments, extending ab initio many-body methods to increasingly complex systems, have been implemented in the massively parallel code beDeft, dedicated to the study of the electronic properties of large scale systems.
The GW formalism, a Green’s function many-body perturbation theory, is growing in popularity for the description of the electronic properties of condensed matter systems in solid-state physics, and more recently chemistry. Unfortunately, its application to complex systems of interest in nanosciences, chemistry, or even biology, is hampered by the large associated computing cost, in particular in the case of disordered systems, or systems immersed in an opened environment (a solvent, a molecular medium, an electrode, etc.) The goal of the present PhD thesis is to focus on the development of multiscale techniques, merging high-level many-body treatments of the subsystem of interest, with a simplified but fully ab initio description of the electrostatic and dielectric environment. Such approaches aim to go beyond classical parametrized models, mainly developed in the quantum chemistry community, which are based on a continuum (“polarizable continuum model”) or discrete (QM/MM) description of the environment.To reach such a goal, we adopt a divide-and-conquer fragmentation scheme for the environment, particularly suited to molecular systems. This leads to a block-diagonal non- interacting electron susceptibility, decreasing the algorithmic complexity from quartic to cubic. To reduce the prefactor associated with the inversion of the Dyson equation for the screened Coulomb potential W, we have further developed a compression algorithm for the susceptibility operator. The automatic computation of an extremely compact polarization basis set allows a large reduction of the size of the susceptibility blocks, associated to the fragments in the environment. Such a method enables us to compute the dielectric response of systems made of several hundred thousand atoms, with an excellent accuracy when it comes to reproduce the effect of the environment as a response to an excitation in the immersed subsystem. This approach is presented through the study of fullerene bulk, surface and subsurface crystals.While the GW formalism is dynamical, with a frequency-dependent screened Coulomb potential W, a first study is done adopting a static approximation (low-frequency limit) for the screening properties of the environment. Such an approach follows the traditional semi-empirical models of a polarizable environment. This PhD thesis assesses the validity of such an approximation, which assumes an instantaneous response (adiabatic limit) of the environment to an electronic excitation, thanks to an explicit comparison with a fully dynamical dielectric response of the environment. The study of a surface of fullerenes, as well as a water molecule inside a metallic carbon nanotube, show that a static description of the environment leads to errors on the polarization energy below 10%, provided that the “folding” of the environment is treated in a proper way.The fragment approach is also applied to covalent insulator crystals, and more particularly to hexagonal boron nitride (h-BN). We explain how to compute the energy levels of point defects in h-BN, in the true dilute limit, and we give the asymptotic scaling laws for the renormalization of these energy levels, from the monolayer to a (n)-layer system. This study highlights thus the possibility to apply the fragment approach to covalent insulator systems, a possibility hinging probably on the short range behavior of the susceptibility in these systems.All of these developments, extending ab initio many-body methods to increasingly complex systems, have been implemented in the massively parallel code beDeft, dedicated to the study of the electronic properties of large scale systems.
We present a many-body GW formalism for quantum subsystems embedded in discrete polarizable environments containing up to several hundred thousand atoms described at a fully ab initio random phase approximation level. Our approach is based on a fragment approximation in the construction of the Green's function and independent-electron susceptibilities. Further, the environing fragments susceptibility matrices are reduced to a minimal but accurate representation preserving low order polarizability tensors through a constrained minimization scheme. This approach dramatically reduces the cost associated with inverting the Dyson equation for the screened Coulomb potential W, while preserving the description of short to long-range screening effects. The efficiency and accuracy of the present scheme is exemplified in the paradigmatic cases of fullerene bulk, surface, subsurface, and slabs with varying number of layers.
Formalisme à N-corps GW environné dans une approche fragment : développements et applications à des systèmes complexes Le formalisme GW, dans le cadre des théories de perturbation à N-corps utilisant les fonctions de Green, gagne en popularité pour la description des propriétés électroniques des systèmes de la matière condensée en physique du solide, et plus récemment en chimie. Son application à des systèmes complexes d'intérêt en nanoscience, chimie, voire biologie, est freinée cependant par son coût numérique en particulier dans le cas de systèmes désordonnés, ou immergés dans un environnement ouvert (un solvant, un milieu moléculaire, une électrode, etc.) Le but de cette thèse est de développer des techniques multi-échelles, combinant des approches à N-corps de haut niveau pour le sous-système d'intérêt, avec une description simplifiée, mais tout de même totalement ab initio, d'un environnement électrostatique et diélectrique. Ces approches vont donc au-delà des modèles classiques paramétrés, développés en particulier dans la communauté chimie quantique, et basés sur une description continue (« polarizable continuum model ») ou discrète (QM/MM) de l’environnement.Pour atteindre cet objectif, nous adoptons une approche en fragments de l'environnement, particulièrement adaptée aux systèmes moléculaires. La susceptibilité électronique non-interagissante devient ainsi diagonale par blocs, permettant d'abaisser la complexité algorithmique de quartique à cubique. Pour réduire le pré-facteur associé à l’obtention du potentiel écranté W (équation de Dyson), nous avons développé un algorithme de compression de l’opérateur susceptibilité. L’obtention automatique d’une base de polarisation très compacte permet de réduire fortement la taille des blocs de susceptibilité associés aux fragments de l'environnement. Cette méthode permet de calculer la réponse diélectrique de systèmes contenant des centaines de milliers d'atomes avec une excellente précision. Cette approche est présentée via l'étude de cristaux de fullerènes en volume, en surface, et en sous-surface.Alors que le formalisme GW est dynamique par nature, avec ainsi un potentiel coulombien écranté W dépendant de la fréquence, une première étude est réalisée dans le cadre d'une approximation statique (limite basse fréquence) pour décrire l’écrantage par l’environnement. Une telle approche s'inscrit dans la continuité des modèles semi-empiriques traditionnels pour la description d’un milieu environnant polarisable. Cette thèse est donc l'occasion de mesurer la validité d'une telle approximation, qui suppose que l’environnement répond de façon instantanée à une excitation électronique, grâce à une comparaison explicite avec une description totalement dynamique de la réponse diélectrique de l'environnement. L'étude d'une surface de fullerènes, ainsi que d'une molécule d'eau dans un nanotube de carbone métallique, montrent qu'une description statique de l'environnement induit des erreurs sur l'énergie de polarisation inférieure à 10% sous condition que le « repliement » de l’environnement soit correctement effectué.L'approche fragment est également appliquée à des cristaux covalents isolants, et en particulier au nitrure de bore hexagonal (h-BN). Nous avons illustré en particulier comment calculer les niveaux d’énergie de défauts ponctuels dans h-BN, dans la vraie limite diluée, et donnons les lois d’échelle pour la renormalisation de ces niveaux de la monocouche vers un nombre (n) de couches. Cette étude démontre qu’à l’instar des systèmes moléculaires, la fragmentation de systèmes covalents isolants est possible, en lien sans doute avec le caractère très courte portée de la susceptibilité dans ces systèmes.Ces développements, permettant l’extension d’approches quantiques à N-corps à des systèmes de plus en plus complexes, ont été implémentés dans le code beDeft, un code massivement parallèle pour l’étude des propriétés électroniques de systèmes de grande taille.
We present a fragment-based many-body Green's function method suitable for treating large molecular systems in heterogeneous polarizable environments. The Green's function for a total system is approximated from fragment Green's functions and is expanded up to two-body terms. The screened Coulomb potential is approximated from the sum of intrafragment density-response functions, with interfragment polarization terms being neglected. The approximations for the Green's function and screened Coulomb potential lead to a many-body expansion of the self-energy. This expansion is essentially equivalent to the many-body expansion of the Fock matrix in the fragment molecular orbital method. To handle large molecular systems, the present implementation relies on the Coulomb hole plus screened exchange (COHSEX) approximation for the GW self-energy. The accuracy of the FMO-COHSEX method was demonstrated in comparison to conventional COHSEX results for organic molecular aggregates. We confirmed that the present fragmentation approximation can provide reasonably accurate results, and mean absolute errors for quasiparticle energies of less than 0.1 eV have been achieved for valence orbitals. We also assessed the accuracy of the COHSEX approximation to describe the effects of molecular aggregation of electronic states, by comparing them with the GW method. The COHSEX approximation has been shown to successfully describe the induced polarization and dispersion effects. As an illustrative application of the present method, we considered the electronic states of the pentacene thin film, which contains 1476 atoms. We investigated the impact of the induced polarization effect in the heterogeneous environment, highlighting the gap renormalization and the polarization-induced localization. This application shows that the present fragment-based method is useful for studying electronic structures of molecular aggregates in complex environments.
Electronic structure of materials from first principles is essential for one-to-one comparisons between theories and experiments. However, the many-body effects resulting from electron correlation and their interplay with other degrees of freedom pose a great challenge to design an ab initio approximate method with good accuracy and efficiency at the same time. In this thesis, we present our recent developments of ab initio many-body theory for solids based on the finite-temperature Green’s function formalism. We begin by reviewing the basic concepts and definitions of realistic electronic Hamiltonians (Ch. 2) as well as Green’s functions at finite temperature (Ch. 3). Overviews of self-consistent many-body perturbation theories and quantum embedding theories are then presented in Chs. 4 and 5 as complementary materials for the following chapters. In Ch. 6, we present algorithms and implementation details for the finite-temperature self-consistent GW (scGW) method in Gaussian-type orbitals for solids. Our implementation is based on the finite-temperature Green’s function formalism in which all equations are solved on the imaginary axis, without resorting to analytical continuation during the self-consistency. No quasiparticle approximation is employed and all matrix elements of self-energies are explicitly evaluated. We show agreement with other, differently formulated finite-temperature scGW implementations. By migrating computationally intensive calculations to GPUs, we obtain scalable results on large supercomputers with nearly optimal performance. Inclusion of relativistic effects in scGW is discussed in Ch. 7. Our formulation of relativistic scGW is based on the exact two-component formalism with one-electron approximation (X2C1e) and non-relativistic Coulomb interactions. Our theory allows us to study scalar relativistic effects, spin-orbit coupling, and their interplay with electron correlation without adjustable parameters. We examine the effect of the X2C1e by comparison to the established four-component formalism and reach excellent agreement. The simplicity of X2C1e enables construction of higher order theories, such as embedding theories, on top of perturbative calculations. In Ch. 8, we use self-energy embedding theory (SEET) to study the spectra of the prototypical undistorted cubic perovskites SrVO3 and SrMnO3. In the strongly correlated metallic SrVO3, we find that the usual attribution of the satellite peaks at -1.8 eV to Hund or Hubbard physics in the t2g orbitals is inconsistent with our calculations. In the strongly correlated insulator SrMnO3 we recover insulating behavior due to a feedback effect between the strongly correlated orbitals and the weakly correlated environment. Our calculation shows a systematic convergence of spectral features as the space of strongly correlated orbitals is enlarged, paving the way to a systematic parameter-free study of correlated perovskites. Lastly, in Ch. 9, we apply the Green’s function coupled cluster singles and doubles (GFCCSD) impurity solver to realistic impurity problems in strongly correlated solids within the framework of SEET. We describe the implementation details of our Green’s function coupled cluster (GFCC) solver, investigate its performance, and highlight potential advantages and problems on examples of impurities created for antiferromagnetic MnO and paramagnetic SrMnO3. GFCCSD provides satisfactory descriptions for weakly and moderately correlated impurities with sizes that are intractable by existing accurate impurity solvers such as exact diagonalization. However, when correlations become strong, the singles and doubles approximation used in GFCC could lead to instabilities in searching for the particle number present in impurity problems.
Conspectus Ab initio modeling of matter has become a pillar of chemical research: with ever-increasing computational power, simulations can be used to accurately predict, for example, chemical reaction rates, electronic and mechanical properties of materials, and dynamical properties of liquids. Many competing quantum mechanical methods have been developed over the years that vary in computational cost, accuracy, and scalability: density functional theory (DFT), the workhorse of solid-state electronic structure calculations, features a good compromise between accuracy and speed. However, approximate exchange-correlation functionals limit DFT's ability to treat certain phenomena or states of matter, such as charge-transfer processes or strongly correlated materials. Furthermore, conventional DFT is purely a ground-state theory: electronic excitations are beyond its scope. Excitations in molecules are routinely calculated using time-dependent DFT linear response; however applications to condensed matter are still limited. By contrast, many-electron wavefunction methods aim for a very accurate treatment of electronic exchange and correlation. Unfortunately, the associated computational cost renders treatment of more than a handful of heavy atoms challenging. On the other side of the accuracy spectrum, parametrized approaches like tight-binding can treat millions of atoms. In view of the different (dis-)advantages of each method, the simulation of complex systems seems to force a compromise: one is limited to the most accurate method that can still handle the problem size. For many interesting problems, however, compromise proves insufficient. A possible solution is to break up the system into manageable subsystems that may be treated by different computational methods. The interaction between subsystems may be handled by an embedding formalism. In this Account, we review embedded correlated wavefunction (CW) approaches and some applications. We first discuss our density functional embedding theory, which is formally exact. We show how to determine the embedding potential, which replaces the interaction between subsystems, at the DFT level. CW calculations are performed using a fixed embedding potential, that is, a non-self-consistent embedding scheme. We demonstrate this embedding theory for two challenging electron transfer phenomena: (1) initial oxidation of an aluminum surface and (2) hot-electron-mediated dissociation of hydrogen molecules on a gold surface. In both cases, the interaction between gas molecules and metal surfaces were treated by sophisticated CW techniques, with the remainder of the extended metal surface being treated by DFT. Our embedding approach overcomes the limitations of conventional Kohn-Sham DFT in describing charge transfer, multiconfigurational character, and excited states. From these embedding simulations, we gained important insights into fundamental processes that are crucial aspects of fuel cell catalysis (i.e., O2 reduction at metal surfaces) and plasmon-mediated photocatalysis by metal nanoparticles. Moreover, our findings agree very well with experimental observations, while offering new views into the chemistry. We finally discuss our recently formulated potential-functional embedding theory that provides a seamless, first-principles way to include back-action onto the environment from the embedded region.
Understanding emergent many-body phenomena in correlated materials remains one of the grandest challenges in solid-state quantum chemistry calculations. The individuality of electrons diminishes as correlation increases, which leads to the emergence of new features in the system. Seeking to address such problems accurately requires the aid of non-perturbative many-body techniques as well as the corresponding high-quality basis sets. In periodic systems, a given element may be present in different spatial arrangements displaying vastly different physical and chemical properties, and an elemental basis set independent of the physical properties of materials may lead to significant simulation inaccuracies. In fact, with the rapid progress of quantum chemistry methods in condensed-phase simulations, the need for a library of reliable Gaussian basis sets explicitly designed for periodic calculations has become urgent. In the meantime, the development of parameter-free, systematic and reliable quantum chemistry methods for simultaneously treating weak and strong correlations in periodic systems has never stopped. Their ability to treat realistic solid materials needs to be tested, and the obtained results can serve as benchmarks for empirical parameter-guided solid-state calculations. In this thesis, we present numerical studies of various solid-state systems, emphasizing the design and optimization of Gaussian basis sets and assessing the band structure properties of the periodic solid systems. The work begins by introducing the theoretical framework of the electronic structure theory and presenting the fundamentals of a fully self-consistent GW approximation and the self-energy embedding theory (SEET) in Chapters 2 and 3. It continues with a description of the Gaussian basis set optimization scheme we devised for solid-state quantum chemistry calculations. The scheme we present is designed to avoid a lack of material specificity within a given basis set by simultaneously minimizing the total energy of the system and optimizing the band energies when compared to the reference plane wave calculation while accounting for the overlap matrix condition number. We compare the quality of the Gaussian basis sets generated via our method against the existing basis sets. And the optimization scheme is tested to yield improved results. We finally present a quantitative study of electronic properties of realistic materials with the treatment of Green's function based weakly correlated self-consistent GW approximation and the strongly correlated self-energy embedding theory. In the transition metal oxide BiVO3, we found the inaccurate illustration of orbitals including t_{2g} in transition metal Vanadium by DFT. Our calculation shows a systematic trend of DFT failure in illustrating both insulating and metallic solutions, demonstrating the dangers of using DFT when the experimental data is scarce.