We develop a multi-reference perturbation theory for electronic structure calculations based on symmetries of the Hamiltonian. The reference Hamiltonian in the symmetry-based perturbation theory (SBPT) is chosen such that it possesses more symmetries than the original Hamiltonian, leading to a larger reduction of computational resources in terms of both the number of configurations in the configuration interaction expansion and the number of required qubits in quantum computing applications. We provide approximate, scalable solutions for the second-order correction, as well as an application to selected configuration interaction. We show that SBPT is an extension of other existing multi-reference perturbation theories and that it can give better results for some molecular systems in a robust way.
We develop a multi-reference perturbation theory for electronic structure calculations based on symmetries of the Hamiltonian. The reference Hamiltonian in the symmetry-based perturbation theory (SBPT) is chosen such that it possesses more symmetries than the original Hamiltonian, leading to a larger reduction of computational resources in terms of both the number of configurations in the configuration interaction expansion and the number of required qubits in quantum computing applications. We provide approximate, scalable solutions for the second-order correction, as well as an application to selected configuration interaction. We show that SBPT is an extension of other existing multi-reference perturbation theories and that it can give better results for some molecular systems in a robust way.
Quantum computation is one of the most promising new paradigms for the simulation of physical systems composed of electrons and atomic nuclei, with applications in chemistry, solid-state physics, materials science, and molecular biology. This requires a truncated representation of the electronic structure Hamiltonian using a finite number of orbitals. While it is, in principle, obvious how to improve on the representation by including more orbitals, this is usually unfeasible in practice (e.g., because of the limited number of qubits available) and severely compromises the accuracy of the obtained results. Here, we propose a quantum algorithm that improves on the representation of the physical problem by virtue of second-order perturbation theory. In particular, our quantum algorithm evaluates the second-order energy correction through a series of time-evolution steps under the unperturbed Hamiltonian. An important application is to go beyond the active-space approximation, allowing to include corrections of virtual orbitals, known as multireference perturbation theory. Here, we exploit the fact that the unperturbed Hamiltonian is diagonal for virtual orbitals and show that the number of qubits is independent of the number of virtual orbitals. This gives rise to more accurate energy estimates without increasing the number of qubits. Moreover, we demonstrate numerically for realistic chemical systems that the total runtime has highly favorable scaling in the number of virtual orbitals compared to previous studies. Numerical calculations confirm the necessity of the multireference perturbation theory energy corrections to reach accurate ground-state energy estimates. Our perturbation theory quantum algorithm can also be applied to symmetry-adapted perturbation theory. As such, we demonstrate that perturbation theory can help to reduce the quantum hardware requirements for quantum chemistry. Published by the American Physical Society 2024
The solution of the quantum many-body problem for medium-mass nuclei using realistic nuclear interactions poses a superbe challenge for nuclear structure research. Because an exact solution can only be provided for the lightest nuclei, one has to rely on approximate solutions when proceeding to heavier systems. Over the past years, tremendous progress has been made in the development and application of systematically improvable expansion methods and an accurate description of nuclear observables has become viable up to mass number A≈100. While closed-shell systems are consistently described via a plethora of different many-body methods, the extension to genuine open-shell systems still remains a major challenge and up to now there is no \textit{ab initio} many-body method which applies equally well to systems with even and odd mass numbers. The goal of this thesis is the development and implementation of innovative perturbative approaches with genuine open-shell capabilities. This requires the extension of well-known single-reference approaches to more general vacua. In this work we choose two complementary routes for the usage of generalized reference states. First, we derive a new \emph{ab initio} approach based on multi-configurational reference states that are conveniently derived from a prior no-core shell model calculation. Perturbative corrections are derived via second-order many-body perturbation theory, thus, merging configuration interaction and many-body perturbation theory. The generality of this ansatz enables for a treatment of medium-mass systems with arbitrary mass number, as well as the extension to low-lying excited states such that ground and excited states are treated on an equal footing. In a complementary approach, we use reference states that break a symmetry of the underlying Hamiltonian. In the simplest case this corresponds to the expansion around a particle-number-broken Hartree-Fock-Bogoliubov vacuum which is obtained from a mean-field calculation. Pairing correlations are already absorbed in the reference state. The mild scaling allows for investigations up to tin isotopic chains. All benchmarks were performed using state-of-the-art chiral two- plus three-nucleon interactions thus allowing for a universal description of nuclear observables in the medium-mass regime.
Abstract The accurate electronic structure calculation for strongly correlated chemical systems requires an adequate description for both static and dynamic electron correlation, and is a persistent challenge for quantum chemistry. In order to account for static and dynamic electron correlations accurately and efficiently, in this work we propose a new method by integrating the density matrix renormalization group (DMRG) method and multi-reference second-order Epstein–Nesbet perturbation theory (ENPT2) with a selected configuration interaction (SCI) approximation. Compared with previous DMRG-based dynamic correlation methods, the DMRG-ENPT2 method extends the range of applicability, allowing us to efficiently calculate systems with very large active space beyond 30 orbitals. We demonstrate this by performing calculations on H 2 S with an active space of (16e, 15o), hexacene with an active space of (26e, 26o) and 2D H 64 square lattice with an active space of (42e, 42o).
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).
This work describes resource‑efficient quantum algorithms for molecular physics that exploit geometric structure and correspondences between physical and qubit symmetries. It begins with a review of classical approaches to the molecular electronic‑structure problem and of the variational quantum eigensolver in quantum computing. We then discuss the symmetry‑adapted encodings we introduced to reduce qubit requirements by eliminating redundant degrees of freedom arising from point‑group and other Boolean symmetries in the quantum simulation of molecular systems. Building on these encodings, we describe a qubit‑based symmetry‑adapted complete active‑space method that permits the simulation of larger molecular systems in a resource‑efficient manner on near‑term quantum devices. We present the QuantumSymmetry software package and illustrate its application to representative molecular systems. Next, we address electron‑molecule scattering by introducing a quantum computing algorithm to solve the inner‑region problem within the R‑matrix framework. In particular, we propose a novel simultaneous diagonalisation algorithm with subspace optimisation that requires only a limited number of Hamiltonian evaluations and could be generalised to a wide range of problems. The final chapter derives the explicit form of qubit operators for particle‑number projection. Altogether, these developments advance practical quantum simulations of complex molecular phenomena by leveraging symmetry and geometry in quantum‑algorithm design.
Quantum-selected configuration interaction (QSCI) is a novel quantum-classical hybrid algorithm for quantum chemistry calculations. This method identifies electron configurations having large weights for the target state using quantum devices and allows CI calculations to be performed with the selected configurations on classical computers. In principle, the QSCI algorithm can take advantage of the ability to handle large configuration spaces while reducing the negative effects of noise on the calculated values. At present, QSCI calculations are limited by qubit noise during the input state preparation and measurement processes, restricting them to small active spaces. These limitations make it difficult to perform calculations with quantitative accuracy. The present study demonstrates a computational scheme based on multireference perturbation theory calculations on a classical computer, using the QSCI wave function as a reference. This method was applied to ground and excited state calculations for two typical aromatic molecules, naphthalene and tetracene. The incorporation of the perturbation treatment was found to provide improved accuracy. Extension of the reference space based on the QSCI-selected configurations as a means of further improvement was also investigated.