Recent mathematical advances in coupled cluster theory

This article presents an in-depth educational overview of the latest mathematical developments in coupled cluster (CC) theory, beginning with Schneider's seminal work from 2009 that introduced the first local analysis of CC theory. We offer a tutorial review of second quantization and the CC ansatz, laying the groundwork for understanding the mathematical basis of the theory. This is followed by a detailed exploration of the most recent mathematical advancements in CC theory.Our review starts with an in-depth look at the local analysis pioneered by Schneider which has since been applied to analyze various CC methods. We then move on to discuss the graph-based framework for CC methods developed by Csirik and Laestadius. This framework provides a comprehensive platform for comparing different CC methods, including multireference approaches. Next, we delve into the latest numerical analysis results analyzing the single reference CC method developed by Hassan, Maday, and Wang. This very general approach is based on the invertibility of the CC function's Fr\'echet derivative. We conclude the article with a discussion on the recent incorporation of algebraic geometry into CC theory, highlighting how this novel and fundamentally different mathematical perspective has furthered our understanding and provides exciting pathways to new computational approaches

One-Dimensional Lieb-Oxford Bounds

We investigate and prove Lieb-Oxford bounds in one dimension by studying convex potentials that approximate the ill-defined Coulomb potential. A Lieb-Oxford inequality establishes a bound of the indirect interaction energy for electrons in terms of the one-body particle density $\rho_\psi$ of a wave function $\psi$. Our results include modified soft Coulomb potential and regularized Coulomb potential. For these potentials, we establish Lieb-Oxford-type bounds utilizing logarithmic expressions of the particle density. Furthermore, a previous conjectured form $I_\mathrm{xc}(\psi)\geq - C_1 \int_{\mathbb R} \rho_\psi(x)^{2} \mathrm{d}x$ is discussed for different convex potentials

Homotopy continuation methods for coupled-cluster theory in quantum chemistry

Homotopy methods have proven to be a powerful tool for understanding the multitude of solutions provided by the coupled-cluster polynomial equations. This endeavor has been pioneered by quantum chemists that have undertaken both elaborate numerical as well as mathematical investigations. Recently, from the perspective of applied mathematics, new interest in these approaches has emerged using both topological degree theory and algebraically oriented tools. This article provides an overview of describing the latter development

The Coupled-Cluster Formalism - a Mathematical Perspective

The Coupled-Cluster theory is one of the most successful high precision methods used to solve the stationary Schr\"odinger equation. In this article, we address the mathematical foundation of this theory with focus on the advances made in the past decade. Rather than solely relying on spectral gap assumptions (non-degeneracy of the ground state), we highlight the importance of coercivity assumptions - G\aa rding type inequalities - for the local uniqueness of the Coupled-Cluster solution. Based on local strong monotonicity, different sufficient conditions for a local unique solution are suggested. One of the criteria assumes the relative smallness of the total cluster amplitudes (after possibly removing the single amplitudes) compared to the G\aa rding constants. In the extended Coupled-Cluster theory the Lagrange multipliers are wave function parameters and, by means of the bivariational principle, we here derive a connection between the exact cluster amplitudes and the Lagrange multipliers. This relation might prove useful when determining the quality of a Coupled-Cluster solution. Furthermore, the use of an Aubin-Nitsche duality type method in different Coupled-Cluster approaches is discussed and contrasted with the bivariational principle

Coupled cluster theory: Towards an algebraic geometry formulation

Coupled cluster theory produced arguably the most widely used high-accuracy computational quantum chemistry methods. Despite the approach's overall great computational success, its mathematical understanding is so far limited to results within the realm of functional analysis. The coupled cluster amplitudes, which are the targeted objects in coupled cluster theory, correspond to solutions to the coupled cluster equations, which is a system of polynomial equations of at most degree four. The high dimensionality of the electronic Schr\"odinger equation and the non-linearity of the coupled cluster ansatz have so far stalled a formal analysis of this polynomial system. In this article, we present algebraic investigations that shed light on the coupled cluster equations and the root structure of this ansatz. This is of importance for the a posteriori evaluation of coupled cluster calculations. To that end, we investigate the root structure by means of Newton polytopes. We derive a general v-description, which is subsequently turned into an h-description for explicit examples. This perspective reveals an apparent connection between Pauli's exclusion principle and the geometrical structure of the Newton polytopes. We also propose an alternative characterization of the coupled cluster equations projected onto singles and doubles as cubic polynomials on an algebraic variety with certain sparsity patterns. Moreover, we provide numerical simulations of two computationally tractable systems, namely, the two electrons in four spin-orbitals system and the three electrons in six spin-orbitals system. These simulations provide novel insight into the root structure of the coupled cluster solutions when the coupled cluster ansatz is truncated

Analysis of The Tailored Coupled-Cluster Method in Quantum Chemistry

In quantum chemistry, one of the most important challenges is the static correlation problem when solving the electronic Schr\"odinger equation for molecules in the Born--Oppenheimer approximation. In this article, we analyze the tailored coupled-cluster method (TCC), one particular and promising method for treating molecular electronic-structure problems with static correlation. The TCC method combines the single-reference coupled-cluster (CC) approach with an approximate reference calculation in a subspace [complete active space (CAS)] of the considered Hilbert space that covers the static correlation. A one-particle spectral gap assumption is introduced, separating the CAS from the remaining Hilbert space. This replaces the nonexisting or nearly nonexisting gap between the highest occupied molecular orbital and the lowest unoccupied molecular orbital usually encountered in standard single-reference quantum chemistry. The analysis covers, in particular, CC methods tailored by tensor-network states (TNS-TCC methods). The problem is formulated in a nonlinear functional analysis framework, and, under certain conditions such as the aforementioned gap, local uniqueness and existence are proved using Zarantonello's lemma. From the Aubin--Nitsche-duality method, a quadratic error bound valid for TNS-TCC methods is derived, e.g., for linear-tensor-network TCC schemes using the density matrix renormalization group method

The SS-diagnostic -- an a posteriori error assessment for single-reference coupled-cluster methods

We propose a novel a posteriori error assessment for the single-reference coupled-cluster (SRCC) method called the $S$-diagnostic. We provide a derivation of the $S$-diagnostic that is rooted in the mathematical analysis of different SRCC variants. We numerically scrutinized the $S$-diagnostic, testing its performance for (1) geometry optimizations, (2) electronic correlation simulations of systems with varying numerical difficulty, and (3) the square-planar copper complexes [CuCl$_4$]$^{2-}$, [Cu(NH$_3$)$_4$]$^{2+}$, and [Cu(H$_2$O)$_4$]$^{2+}$. Throughout the numerical investigations, the $S$-diagnostic is compared to other SRCC diagnostic procedures, that is, the $T_1$, $D_1$, and $D_2$ diagnostics as well as different indices of multi-determinantal and multi-reference character in coupled-cluster theory. Our numerical investigations show that the $S$-diagnostic outperforms the $T_1$, $D_1$, and $D_2$ diagnostics and is comparable to the indices of multi-determinantal and multi-reference character in coupled-cluster theory in their individual fields of applicability. The experiments investigating the performance of the $S$-diagnostic for geometry optimizations using SRCC reveal that the $S$-diagnostic correlates well with different error measures at a high level of statistical relevance. The experiments investigating the performance of the $S$-diagnostic for electronic correlation simulations show that the $S$-diagnostic correctly predicts strong multi-reference regimes. The $S$-diagnostic moreover correctly detects the successful SRCC computations for [CuCl$_4$]$^{2-}$, [Cu(NH$_3$)$_4$]$^{2+}$, and [Cu(H$_2$O)$_4$]$^{2+}$, which have been known to be misdiagnosed by $T_1$ and $D_1$ diagnostics in the past. This shows that the $S$-diagnostic is a promising candidate for an a posteriori diagnostic for SRCC calculations

Interacting models for twisted bilayer graphene: a quantum chemistry approach

The nature of correlated states in twisted bilayer graphene (TBG) at the magic angle has received intense attention in recent years. We present a numerical study of an interacting Bistritzer-MacDonald (IBM) model of TBG using a suite of methods in quantum chemistry, including Hartree-Fock, coupled cluster singles, doubles (CCSD), and perturbative triples (CCSD(T)), as well as a quantum chemistry formulation of the density matrix renormalization group method (DMRG). Our treatment of TBG is agnostic to gauge choices, and hence we present a new gauge-invariant formulation to detect the spontaneous symmetry breaking in interacting models. To benchmark our approach, we focus on a simplified spinless, valleyless IBM model. At integer filling ($\nu=0$), all numerical methods agree in terms of energy and $C_{2z} \mathcal{T}$ symmetry breaking. Additionally, as part of our benchmarking, we explore the impact of different schemes for removing ``double-counting'' in the IBM model. Our results at integer filling suggest that cross-validation of different IBM models may be needed for future studies of the TBG system. After benchmarking our approach at integer filling, we perform the first systematic study of the IBM model near integer filling (for $|\nu|< 0.2$). In this regime, we find that the ground state can be in a metallic and $C_{2z} \mathcal{T}$ symmetry breaking phase. The ground state appears to have low entropy, and therefore can be relatively well approximated by a single Slater determinant. Furthermore, we observe many low entropy states with energies very close to the ground state energy in the near integer filling regime