Manipulating expressions in many-body perturbation theory becomes unwieldily with increasing order of the perturbation theory. Here I derive a set of theorems for efficient simplification of such expressions. The derived rules are specifically designed for implementing with symbolic algebra tools. As an illustration, we count the numbers of Brueckner-Goldstone diagrams in the first several orders of many-body perturbation theory for matrix elements between two states of a mono-valent system.
Manipulating expressions in many-body perturbation theory becomes unwieldily with increasing order of the perturbation theory. Here I derive a set of theorems for efficient simplification of such expressions. The derived rules are specifically designed for implementing with symbolic algebra tools. As an illustration, we count the numbers of Brueckner-Goldstone diagrams in the first several orders of many-body perturbation theory for matrix elements between two states of a mono-valent system.
In the context of many-body quantum mechanics, one often encounters the need to evaluate expectation values of operators or their matrix elements. For instance, given a Hamiltonian https://www.w3.org/1998/Math/MathML" display="inline"> H ^ , we typically want to evaluate its expectation value, which is the energy of the system, https://www.w3.org/1998/Math/MathML" display="inline"> E = ⟨ H ^ ⟩ . In second quantization, the matrix element of an operator can always be expressed as a vacuum expectation value. This generally yields an expression which is hard to evaluate. Fortunately, a powerful tool to simplify these calculations exists: Wick's theorem. This theorem provides a systematic way to break down expressions involving products of creation and annihilation operators into sums of simpler terms.
Computational many-body theories in quantum chemistry, condensed matter, and nuclear physics aim to provide sufficiently accurate description of and insights into the motion of many interacting particles. Due to their intrinsic complexity, the development of such theories generally involves very complex, tedious, and error-prone symbolic manipulations. Here a complete solution to automate the symbolics in many-body theory development is attempted. General data structures based on an existing computer algebra system are designed to specifically address the symbolic problems for which there is currently no satisfactory handling. Based on the data structures, algorithms are given to accomplish common symbolic manipulations and simplifications. Noncommutative algebraic systems, tensors with symmetry, and symbolic summations can all enjoy deep simplifications efficient enough for theories of very complex form. After the symbolic derivation, novel algorithms for automatic optimization of tensor contractions and their sums are devised, which can be used together with automatic code generation tools. In this way, the burden of symbolic tasks in theory development can be vastly reduced, with the potential to spare scientists more time and energy for the actual art and science of many-body theories.
Scope The book focuses primarily on many-body (or better, many-electron) methods for electron correlation. These include Rayleigh–Schrödinger perturbation theory (RSPT) , particularly in its diagrammatic representation (referred to as many-body perturbation theory , or MBPT ), and coupled-cluster (CC) theory; their relationship to configuration interaction (CI) is included. Further extensions address properties other than the energy, and also excited states and multireference CC and MBPT methods. The many-body algebraic and diagrammatic methods used in electronic structure theory have their origin in quantum field theory and in the study of nuclear matter and nuclear structure. The second-quantization formalism was first introduced in a treatment of quantized fields by Dirac (1927) and was extended to fermion systems by Jordan and Klein (1927) and by Jordan and Wigner (1928). This formalism is particularly useful in field theory, in scattering problems and in the study of infinite systems because it easily handles problems involving infinite, indefinite or variable numbers of particles. The diagrammatic approach was introduced into field theory by Feynman (1949a,b) and applied to many-body systems by Hugenholtz (1957) and by Goldstone (1957). Many-body perturbation theory and its linked-diagram formalism were first introduced by Brueckner and Levinson (1955) and by Brueckner (1955), and were formalized by Goldstone (1957). Other important contributions to the methodology, first in field theory and then in the theory of nuclear structure, are due to Dyson (1949a,b), Wick (1950), Hubbard (1957, 1958a,b) and Frantz and Mills (1960). Applications to the electronic structure of atoms and molecules began with the work of Kelly (1963, 1964a,b, 1968), and molecular applications using finite analytical basis sets appeared in the work of Bartlett and Silver (1974a, b).
Matrix elements between nonorthogonal Slater determinants represent an essential component of many emerging electronic structure methods. However, evaluating nonorthogonal matrix elements is conceptually and computationally harder than their orthogonal counterparts. While several different approaches have been developed, these are predominantly derived from the first-quantized generalized Slater-Condon rules and usually require biorthogonal occupied orbitals to be computed for each matrix element. For coupling terms between nonorthogonal excited configurations, a second-quantized approach such as the nonorthogonal Wick's theorem is more desirable, but this fails when the two reference determinants have a zero many-body overlap. In this contribution, we derive an entirely generalized extension to the nonorthogonal Wick's theorem that is applicable to all pairs of determinants with nonorthogonal orbitals. Our approach creates a universal methodology for evaluating any nonorthogonal matrix element and allows Wick's theorem and the generalized Slater-Condon rules to be unified for the first time. Furthermore, we present a simple well-defined protocol for deriving arbitrary coupling terms between nonorthogonal excited configurations. In the case of overlap and one-body operators, this protocol recovers efficient formulas with reduced scaling, promising significant computational acceleration for methods that rely on such terms.
Background There are two stages in the study of perturbation theory and related techniques (although they are mixed intimately in most derivations in the literature). The first is the formal development , carried out in terms of the total Hamiltonian and total wave function (and total zero-order wave function), without attempt to express anything in terms of one- and two-body quantities (components of Ĥ, orbitals, integrals over orbitals etc.). We can make a considerable amount of progress in this way before considering the detailed form of Ĥ. The second is the many-body development , where all expressions are obtained in terms of orbitals (one-electron states) and oneand two-electron integrals. We shall try to keep these separate for a while and begin with a consideration of formal perturbation theory. Another aspect of the study of many-body techniques is the large variety of approaches, notations and derivations that have been used. Each different approach has contributed to the lore and the language of many-body theory, and each tends to illuminate some aspects better than the other approaches. If we want to be able to read the literature in this field, we should be familiar with several alternative formulations. Therefore, we shall occasionally derive some results in more than one way and, in particular, we shall derive the basic perturbation-theory equations and their many-body representations in several complementary ways. Classical derivation of Rayleigh–Schrödinger perturbation theory The perturbation Ansatz We begin with a classical textbook derivation of formal Rayleigh–Schrödinger perturbation theory (RSPT).