A realistic quantum many-body system, characterized by a generic microscopic Hamiltonian, is accessible only through approximation methods. The mean field theories, as the simplest practices of approximation methods, commonly serve as a powerful tool, but unfortunately often violate the symmetry of the Hamiltonian. The conventional BCS theory, as an excellent mean field approach, violates the particle number conservation and completely erases quantumness characterized by concurrence and quantum discord between different modes. We restore the symmetry by using the projected BCS theory and the exact numerical solution and find that the lost quantumness is synchronously reestablished. We show that while entanglement remains unchanged with the particle numbers, quantum discord behaves as an extensive quantity with respect to the system size. Surprisingly, discord is hardly dependent on the interaction strengths. The new feature of discord offers promising applications in modern quantum technologies.
A realistic quantum many-body system, characterized by a generic microscopic Hamiltonian, is accessible only through approximation methods. The mean field theories, as the simplest practices of approximation methods, commonly serve as a powerful tool, but unfortunately often violate the symmetry of the Hamiltonian. The conventional BCS theory, as an excellent mean field approach, violates the particle number conservation and completely erases quantumness characterized by concurrence and quantum discord between different modes. We restore the symmetry by using the projected BCS theory and the exact numerical solution and find that the lost quantumness is synchronously reestablished. We show that while entanglement remains unchanged with the particle numbers, quantum discord behaves as an extensive quantity with respect to the system size. Surprisingly, discord is hardly dependent on the interaction strengths. The new feature of discord offers promising applications in modern quantum technologies.
The mean-field approximation based on effective interactions or density functionals plays a pivotal role in the description of finite quantum many-body systems that are too large to be treated by ab initio methods. Some examples are strongly interacting medium and heavy mass atomic nuclei and mesoscopic condensed matter systems. In this approach, the linear Schrodinger equation for the exact many-body wave function is mapped onto a non-linear density-dependent one-body potential problem. This approximation, not only provides computationally very simple solutions even for systems with many particles, but due to the non-linearity, it also allows for obtaining solutions that break essential symmetries of the system, often connected with phase transitions. In this way, additional correlations are subsumed in the system. However, the mean-field approach suffers from the drawback that the corresponding wave functions do not have sharp quantum numbers and, therefore, many results cannot be compared directly with experimental data. In this article, we discuss general group-theory techniques to restore the broken symmetries, and provide detailed expressions on the restoration of translational, rotational, spin, isospin, parity and gauge symmetries, where the latter corresponds to the restoration of the particle number. In order to avoid the numerical complexity of exact projection techniques, various approximation methods available in the literature are examined. Applications of the projection methods are presented for simple nuclear models, realistic calculations in relatively small configuration spaces, nuclear energy density functional theory, as well as in other mesoscopic systems. We also discuss applications of projection techniques to quantum statistics in order to treat the averaging over restricted ensembles with fixed quantum numbers. Further, unresolved problems in the application of the symmetry restoration methods to the energy density functional theories are highlighted in the present work.
Abstract The mean-field approximation based on effective interactions or density functionals plays a pivotal role in the description of finite quantum many-body systems that are too large to be treated by ab initio methods. Some examples are strongly interacting medium and heavy mass atomic nuclei and mesoscopic condensed matter systems. In this approach, the linear Schrödinger equation for the exact many-body wave function is mapped onto a non-linear one-body potential problem. This approximation, not only provides computationally very simple solutions even for systems with many particles, but due to the non-linearity, it also allows for obtaining solutions that break essential symmetries of the system, often connected with phase transitions. In this way, additional correlations are subsumed in the system. However, the mean-field approach suffers from the drawback that the corresponding wave functions do not have sharp quantum numbers and, therefore, many results cannot be compared directly with experimental data. In this article, we discuss general group-theory techniques to restore the broken symmetries, and provide detailed expressions on the restoration of translational, rotational, spin, isospin, parity and gauge symmetries, where the latter corresponds to the restoration of the particle number. In order to avoid the numerical complexity of exact projection techniques, various approximation methods available in the literature are examined. Applications of the projection methods are presented for simple nuclear models, realistic calculations in relatively small configuration spaces, nuclear energy density functional (EDF) theory, as well as in other mesoscopic systems. We also discuss applications of projection techniques to quantum statistics in order to treat the averaging over restricted ensembles with fixed quantum numbers. Further, unresolved problems in the application of the symmetry restoration methods to the EDF theories are highlighted in the present work.
The dynamics of interacting many-particle systems on quantum mechanical scale is a broad subject of research in modern physics. The understanding of the quantum correlations, or entanglement, between the particles in an isolated many-body quantum system may however be challenging, as the amount of interactions may be large. To be able to describe the system in a simpler way, mean field theory is often applied; the interactions of all individual particles are substituted by an averaged effect. First and second order mean field approximations are applied to the equations of motion of the amplitudes <ân> for each mode of the many-body system. In first order mean field, under the assumption that each mode has the same magnitude of the (constant) amplitude, the Kuramoto equation follows for the phases of the complex amplitudes. To study this first order approximation and to refine the mean field solution, the equations of motion have been extended to second order mean field. After that, numerical integration is used to solve the system of coupled differential equations for two modes. The solutions seem unstable, which is verified by defining quantum fluctuations beyond mean field. These are shown not to be negligible for a two mode system.It was discovered that the mean field equations of motion show divergent fluctuations for a system of two modes, which opens a whole set of questions on the validity of the first order mean field approximation and therefore the use of the Kuramoto model for quantum many-body dynamics. The extension to the second order mean field solution might be promising as this does include correlations between the operators.
We present a pedagogical showcase of BRST quantization for restoring symmetries in nuclear many-body systems as a reformulation and potential alternative to conventional projection methods. The formalism is illustrated using translational invariance for a simple system of two interacting masses in one dimension, but with an eye toward generalizing to more particles, higher dimensions, and other symmetries. We explain the considerations underlying particular choices within the BRST construction, and develop both Hamiltonian and path integral formulations to provide guidance for the variety of many-body and effective field theory contexts where gauge fixing for symmetry restoration might be useful. For the demonstration system we show how to diagonalize within the extended BRST phase space, how variation after projection is recovered for product reference states, how the corresponding gauge-fixed functional integral is constructed, and how collective zero modes are isolated and controlled. Throughout, we keep in mind extensions of BRST symmetry to various approximation schemes as a guide for consistent symmetry restoration.
Background: Ab initio many-body methods have been developed over the past ten years to address mid-mass nuclei. In their best current level of implementation, their accuracy is of the order of a few percent error on the ground-state correlation energy. Recently implemented variants of these methods are operating a breakthrough in the description of medium-mass open-shell nuclei at a polynomial computational cost while putting state-of-the-art models of internucleon interactions to the test.Purpose: As progress in the design of internucleon interactions is made, and as questions one wishes to answer are refined in connection with increasingly available experimental data, further efforts must be made to tailor many-body methods that can reach an even higher precision for an even larger number of observable quantum states or nuclei. The objective of the present work is to contribute to such a quest by designing and testing a new many-body scheme.Methods: We formulate a truncated configuration-interaction method that consists of diagonalizing the Hamiltonian in a highly truncated subspace of the total N-body Hilbert space. The reduced Hilbert space is generated via the particle-number projected BCS state along with projected seniority-zero two- and four-quasiparticle excitations. Furthermore, the extent by which the underlying BCS state breaks U(1) symmetry is optimized in the presence of the projected two- and four-quasiparticle excitations. This constitutes an extension of the so-called restricted variation after projection method in use within the frame of multireference energy density functional calculations. The quality of the newly designed method is tested against exact solutions of the so-called attractive pairing Hamiltonian problem.Results: By construction, the method reproduces exact results for N=2 and N=4. For N=(8,16,20), the error in the ground-state correlation energy is less than (0.006%, 0.1%, 0.15%) across the entire range of internucleon coupling defining the pairing Hamiltonian and driving the normal-to-superfluid quantum phase transition. The presently proposed method offers the advantage of automatic access to the low-lying spectroscopy, which it does with high accuracy.Conclusions: The numerical cost of the newly designed variational method is polynomial (N6) in system size. This method achieves unprecedented accuracy for the ground-state correlation energy, effective pairing gap, and one-body entropy as well as for the excitation energy of low-lying states of the attractive pairing Hamiltonian. This constitutes a sufficiently strong motivation to envision its application to realistic nuclear Hamiltonians in view of providing a complementary, accurate, and versatile ab initio description of mid-mass open-shell nuclei in the future.
This thesis explores the application of the Symmetry-Breaking/Symmetry-Restoration methodology on quantum computers to better approximate a Hamiltonian's ground state energy within a variational framework in many-body physics. This involves intentionally breaking and restoring the symmetries of the wave function ansatz at different stages of the variational search for the ground state. The Variational Quantum Eigensolver (VQE) is utilized for the variational component together with an ansatz inspired by the Bardeen-Cooper-Schrieffer (BCS) theory. The applications were demonstrated using the pairing and Hubbard Hamiltonians. Two approaches were identified with the VQE method: varying the symmetry-breaking ansatz parameters before or after symmetry restoration, termed Quantum Projection After Variation and Quantum Variation After Projection, respectively. The main contribution of this thesis was the development of a variety of symmetry restoration techniques based on the principles of the Quantum Phase Estimation algorithm, the notion of a Quantum "Oracle," and the Classical Shadow formalism. In the final part, hybrid quantum-classical techniques were introduced to extract an approximation of the low-lying spectrum of a Hamiltonian. Assuming accurate Hamiltonian moment extraction from their generating function with a quantum computer, two methods were presented for spectral analysis: the t-expansion method and the Krylov method, which provides, in particular, information about the evolution of the survival probability. Furthermore, the Quantum Krylov method was introduced, offering similar insights without the need to estimate Hamiltonian moments, a task that can be difficult on near-term quantum computers.