In this work, the ground-state energy of the lithium atom is systematically investigated using both time-independent perturbation theory and the variational method to provide a comprehensive pedagogical analysis of many-body atomic systems. The unperturbed Hamiltonian is initially constructed by neglecting electron-electron interactions, treating the system as three independent hydrogen-like electrons to yield a zeroth-order energy baseline of -275.51 eV. The antisymmetric fermionic nature of the exact wave function is rigorously enforced through the Slater determinant formalism. First-order perturbation theory is applied to evaluate static inter-electronic repulsion using exact Coulomb and exchange integrals, refining the energy state to -192.01 eV. To account for dynamical electronic correlation, second-order perturbation theory is computed numerically for virtual single-electron s-orbital transitions, leading to a total perturbative energy of -196.36 eV. A brief discussion of two-electron excitations is also included to encapsulate further physical realism within the framework. Furthermore, a non-orthogonal two-parameter variational approach is employed to model the shell-specific shielding effect. By optimizing the effective nuclear charges, the variational method establishes a superior upper bound energy of -201.187 eV. The results of both methods are comprehensively contrasted against each other and the reference baseline to provide critical insights into the nature of electron correlation and screening in multi-electron atoms.
In this work, the ground-state energy of the lithium atom is systematically investigated using both time-independent perturbation theory and the variational method to provide a comprehensive pedagogical analysis of many-body atomic systems. The unperturbed Hamiltonian is initially constructed by neglecting electron-electron interactions, treating the system as three independent hydrogen-like electrons to yield a zeroth-order energy baseline of -275.51 eV. The antisymmetric fermionic nature of the exact wave function is rigorously enforced through the Slater determinant formalism. First-order perturbation theory is applied to evaluate static inter-electronic repulsion using exact Coulomb and exchange integrals, refining the energy state to -192.01 eV. To account for dynamical electronic correlation, second-order perturbation theory is computed numerically for virtual single-electron s-orbital transitions, leading to a total perturbative energy of -196.36 eV. A brief discussion of two-electron excitations is also included to encapsulate further physical realism within the framework. Furthermore, a non-orthogonal two-parameter variational approach is employed to model the shell-specific shielding effect. By optimizing the effective nuclear charges, the variational method establishes a superior upper bound energy of -201.187 eV. The results of both methods are comprehensively contrasted against each other and the reference baseline to provide critical insights into the nature of electron correlation and screening in multi-electron atoms.
The techniques of many-body Green's function theory have been applied to the calculation of properties of the lithium atom. Using basis set expansion techniques and a second-order expansion of the self-energy, Dyson's equation was solved, and the first-order reduced density matrix and the first ionization potential for the lithium ground state were computed. The natural spin orbitals of s symmetry were in excellent argeement with those computed from wavefunctions. The ionization potential was within 0.01 eV of the exact experimental result.
In chapters 13 and 14 we focused on perturbation theories for time-independent and time-dependent perturbations. For time-independent perturbations we obtained first-order and second-order corrections to energy eigenvalues and eigenfunctions. For time-dependent perturbations we obtained the probability for a system to make a transition from an ith initial state to fth final state. We have also discussed how to study the effect of adiabatic and sudden perturbations. In the present chapter we consider the variational method. The variational method can be used to construct an approximate ground state eigenfunction. In some cases the method can be used for finding approximately the energy of first few discrete states of a bound system. The perturbation theory uses a set of unperturbed states and the perturbation Hamiltonian H(1) to construct the eigenfunctions and eigenvalues through highly formalized approximations. In contrast, the variational method is based on a guess for the state-a trial state which is used to determine the corresponding energy. Any well behaved trial state introduced will yield an appropriate energy, although some trial states will be better than other, but whatever the guess, we can always make systematic improvements. We can construct a reasonably accurate eigenfunction often by a physical insight and ingenuity.
In this lecture we present many-body perturbation theory as a method to determine quasiparticle excitations in solids, especially electronic band structures, accurately from first principles.The main ingredient is the electronic self-energy that, in principle, contains all many-body exchange and correlation effects beyond the Hartree potential.As its exact mathematical expression is unknown, approximations must be used in practical calculations.The approximation is obtained using a systematic algebraic approach on the basis of Green function techniques.It constitutes an expansion of the self-energy up to linear order in the screened Coulomb potential, which describes the interaction between the quasiparticles and includes dynamic screening through the creation of exchange-correlation holes around the bare particles.The implementation of the approximation relies on a perturbative treatment starting from density functional theory.Besides a detailed mathematical discussion we focus on the underlying physical concepts and show some illustrative applications.
In order to properly understand the utility of many-body perturbation theory as applied to finite systems, we use the configuration interaction approach to compute the electronic self-energy. The validity of the commonly used GW approximation from many-body perturbation theory is tested by comparing the self-energy explicitly computed according to the diagrammatic expansions to that from the inversion of the Dyson equation. It is constructed as a functional of the interacting Green function G and screened Coulomb interaction W. W is explicitly computed through the diagonalization of the many-body Hamiltonian and, thus, takes polarization effects into account beyond the random phase approximation. The Na9+ cluster and the dissociation of the C2 molecule are discussed as examples. We find that the GW approximation yields accurate results for weakly correlated systems with a predominantly single-determinant ground state such as Na9+ clusters. The C2 molecule is a pathological system with multideterminantal ground state where the GW approach ceases to be valid.
This dissertation concerns the quantum many-body problem, which is the problem of predicting the properties of systems of several quantum particles from the first principles of quantum mechanics. Included under this umbrella are various problems of fundamental importance in quantum chemistry, condensed matter physics, and materials science. Of particular interest is the electronic structure problem, the problem of determining the state of the electrons in a system with fixed atomic nuclei. Since direct numerical solution of the many-body Schrödinger equation is intractable even for systems of moderate size, a diverse array of approximate methods has been developed. The broad goals of this dissertation are to improve the mathematical understanding of certain widely-used approximations, as well as to propose new methods. Roughly speaking, we consider three (overlapping) categories of methods: Green's function methods, embedding methods, and variational methods.One can understand Green's function methods in terms of many-body perturbation theory, which computes series expansions of physical quantities about a non-interacting reference system. These expansions can be expressed graphically in terms of Feynman diagrams, which can in turn be reorganized, in some cases, into an expansion in terms of so-called bold diagrams. Green's function methods can be specified by choosing a subset of bold diagrams to approximate the sum. At the same time, such methods can be understood in terms of an object known as the Luttinger-Ward (LW) functional, which admits a representation in terms of the bold diagrams. Many aspects of these constructions are purely formal, and indeed the existence of the fermionic LW functional as a single-valued functional has recently been called into question. To contribute to the understanding of these issues, we provide rigorous proofs of the combinatorial construction and analytic interpretation of the bold diagrams in the simplified setting of a classical field theory. In this setting we also provide a rigorous non-perturbative construction of the LW functional via convex duality and prove several key properties, including continuity up to the boundary of its domain and asymptotics in the limit of large interaction.Quantum embedding methods, meanwhile, view a large system as being composed of smaller fragments that are treated with high accuracy and embedded in the larger system in a mutually consistent way. Inspired by a connection between the boundary analysis of the LW functional and embedding, we perform similar analysis for the 1-RDM theory for fermionic systems, which is also developed via convex duality, illustrating a relation to fermionic embedding methods such as the density matrix embedding theory (DMET).Another embedding method of note is the dynamical mean-field theory (DMFT), which is at the same time a Green's function method that can be understood in terms of the LW functional. DMFT relies on the solution of impurity problems, which specify the embedding of an interacting system into a non-interacting bath. Underlying DMFT is a result about the sparsity pattern of the self-energy matrix for impurity problems, which to our knowledge has not been proved in the literature. We provide a rigorous proof of this result in various classical and quantum settings. We go on to investigate the fermionic DMFT in depth, identifying the key mathematical structures that appear in the algorithmic loop for solving it and using these to prove the well-posedness of this loop, in a certain sense.Finally, we introduce a suite of new approaches to the quantum many-body problem that provide variational lower bounds to the ground-state energy. These methods, which combine the themes of convexity and embedding, are based on novel convex relaxations of the variational principles for the ground-state energies of many-body systems. To begin, we recover a second-quantized version of the formalism of strictly correlated electrons (SCE), which yields an exact expression for the exchange-correlation functional in Kohn-Sham density functional theory in the limit of infinite Coulomb repulsion in terms of the solution of a multi-marginal optimal transport problem. We introduce a semidefinite relaxation method for approximately solving this problem and obtaining a lower bound for the ground-state energy. The ideas underlying this relaxation are generalized considerably, outside the context of SCE, to yield much tighter lower bounds, which we validate numerically for both quantum spin systems and fermionic systems. We also describe how these relaxation methods can be interpreted as embedding methods via convex duality.
In this thesis we develop the variational two-electron reduced-density-matrix method for extended systems. Extended systems are represented in two ways: i) lattice models describing the dominant valence electronic structure with periodic boundaries to account for their extended nature and ii) a crystalline-orbital basis built from atomic orbitals using the generalization of molecular orbital theory to polymers. The first part of this thesis (Ch. 3-4) examines the performance of the variational 2-RDM method on lattice systems with tunable electron correlation. The first of these systems is the classic Hubbard model with linear and ladder lattice topologies. Because electron correlation functions, such as charge- and spin-ordering, are linear functions of the 2-RDM, the difference in electronic structure between one- and quasi-one-dimensional systems is accurately characterized. The second model contains only two-body interactions and is unique among typical spin models in that it does not have a mean-field reference wave function. The ground state wave functions from all Hamiltonians in the model have the same 1-electron reduced density matrix; consequently, one-electron theories are largely inapplicable. The superconducting eta-pairing ground states make the model a unique tool for demonstrating the necessary N-representability in highly correlated environments. The second part of this thesis (Ch. 5-6) develops a formalism for modeling materials by solving the full Schrodinger equation. Crystalline-orbital Hartree-Fock provides a set of orbitals and integral tensors for the variational 2-RDM method. We demonstrate that time-reversal symmetry, which is implicitly included in position space electronic structure calculations, must be explicitly included as an N-representability constraint on the 2-RDM when using a momentum space basis. The necessity of these equality constraints is demonstrated by the accurate recovery of the binding energy of two polymers and the symmetry of their natural orbital occupations.