Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userADG: Automated generation and evaluation of many-body diagrams II. Particle-number projected Bogoliubov many-body perturbation theory
108
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We describe the second version (v2.0.0) of the code ADG that automatically (1) generates all valid off-diagonal Bogoliubov many-body perturbation theory diagrams at play in particle-number projected Bogoliubov many-body perturbation theory (PNP-BMBPT) and (2) evaluates their algebraic expression to be implemented for numerical applications. This is achieved at any perturbative order $p$ for a Hamiltonian containing both two-body (four-legs) and three-body (six-legs) interactions (vertices). All valid off-diagonal BMBPT diagrams of order $p$ are systematically generated from the set of diagonal, i.e., unprojected, BMBPT diagrams. The production of the latter were described at length in https://doi.org/10.1016/j.cpc.2018.11.023 dealing with the first version of ADG. The automated evaluation of off-diagonal BMBPT diagrams relies both on the application of algebraic Feynmans rules and on the identification of a powerful diagrammatic rule providing the result of the remaining $p$-tuple time integral. The new diagrammatic rule generalizes the one already identified in https://doi.org/10.1016/j.cpc.2018.11.023 to evaluate diagonal BMBPT diagrams independently of their perturbative order and topology. The code ADG is written in Python3 and uses the graph manipulation package NetworkX. The code is kept flexible enough to be further expanded throughout the years to tackle the diagrammatics at play in various many-body formalisms that already exist or are yet to be formulated.
rate research
Read More
The goal of the present paper is twofold. First, a novel expansion many-body method applicable to superfluid open-shell nuclei, the so-called Bogoliubov in-medium similarity renormalization group (BIMSRG) theory, is formulated. This generalization of standard single-reference IMSRG theory for closed-shell systems parallels the recent extensions of coupled cluster, self-consistent Greens function or many-body perturbation theory. Within the realm of IMSRG theories, BIMSRG provides an interesting alternative to the already existing multi-reference IMSRG (MR-IMSRG) method applicable to open-shell nuclei. The algebraic equations for low-order approximations, i.e., BIMSRG(1) and BIMSRG(2), can be derived manually without much difficulty. However, such a methodology becomes already impractical and error prone for the derivation of the BIMSRG(3) equations, which are eventually needed to reach high accuracy. Based on a diagrammatic formulation of BIMSRG theory, the second objective of the present paper is thus to describe the third version (v3.0.0) of the ADG code that automatically (1) generates all valid BIMSRG(n) diagrams and (2) evaluates their algebraic expressions in a matter of seconds. This is achieved in such a way that equations can easily be retrieved for both the flow equation and the Magnus expansion formulations of BIMSRG. Expanding on this work, the first future objective is to numerically implement BIMSRG(2) (eventually BIMSRG(3)) equations and perform ab initio calculations of mid-mass open-shell nuclei.
While coupled cluster theory accurately models weakly correlated quantum systems, it often fails in the presence of strong correlations where the standard mean-field picture is qualitatively incorrect. In many cases, these failures can be largely ameliorated by permitting the mean-field reference to break physical symmetries. Symmetry-broken coupled cluster, e.g. Bogoliubov coupled cluster, theory can indeed provide reasonably accurate energetic predictions, but the broken symmetry can compromise the quality of the resulting wave function and predictions of observables other than the energy. Merging symmetry projection and coupled cluster theory is therefore an appealing way to describe strongly correlated systems. Independently, two different but related formalisms have been recently proposed to achieve this goal. The two formalisms are contrasted in this manuscript, with results tested on the Richardson pairing Hamiltonian. Both formalisms are based on the disentangled cluster representation of the symmetry-rotated coupled cluster wavefunction. However, they differ in the way that the disentangled clusters are solved. One approach sets up angle-dependent coupled cluster equations, while the other involves first-order ordinary differential equations. The latter approach yields energies and occupation probabilities significantly better than those of number-projected BCS and BCS coupled cluster and, when the disentangled clusters are truncated at low excitation levels, has a computational cost not too much larger than that of BCS coupled cluster. The high quality of results presented in this manuscript indicates that symmetry-projected coupled cluster is a promising method that can accurately describe both weakly and strongly correlated finite many-fermion systems.
We present a new method which uses Feynman-like diagrams to calculate the statistical quantities of embedded many-body random matrix problems. The method provides a promising alternative to existing techniques and offers many important simplifications. We use it here to find the fourth, sixth and eighth moments of the level density for k fermions or bosons interacting through a random hermitian potential in the limit where the number of possible single-particle states is taken to infinity. All share the same transition, starting immediately after 2k = m, from moments arising from a semi-circular level density to gaussian moments. The results also reveal a striking feature; the domain of the 2nth moment is naturally divided into n subdomains specified by the points 2k = m, 3k = m, ..., nk = m.
The order-by-order renormalization of the self-consistent mean-field potential in many-body perturbation theory for normal Fermi systems is investigated in detail. Building on previous work mainly by Balian and de Dominicis, as a key result we derive a thermodynamic perturbation series that manifests the consistency of the adiabatic zero-temperature formalism with perturbative statistical mechanics---for both isotropic and anisotropic systems---and satisfies at each order and for all temperatures the thermodynamic relations associated with Fermiliquid theory. These properties are proved to all orders.
We investigate the order-by-order convergence behavior of many-body perturbation theory (MBPT) as a simple and efficient tool to approximate the ground-state energy of closed-shell nuclei. To address the convergence properties directly, we explore perturbative corrections up to 30th order and highlight the role of the partitioning for convergence. The use of a simple Hartree-Fock solution to construct the unperturbed basis leads to a convergent MBPT series for soft interactions, in contrast to, e.g., a harmonic oscillator basis. For larger model spaces and heavier nuclei, where a direct high-order MBPT calculation in not feasible, we perform third-order calculation and compare to advanced ab initio coupled-cluster calculations for the same interactions and model spaces. We demonstrate that third-order MBPT provides ground-state energies for nuclei up into tin isotopic chain that are in excellent agreement with the best available coupled-cluster results at a fraction of the computational cost.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
