Ab initio Bogoliubov many-body perturbation theory: closed-form constraint on the average particle number

TL;DR

This paper introduces a general formalism for Bogoliubov many-body perturbation theory that allows constraining the average particle number at any order, enabling more accurate ab initio nuclear structure calculations.

Contribution

The paper presents a new algorithm to impose particle number constraints in BMBPT at arbitrary order, overcoming previous limitations to second order.

Findings

Successfully implemented BMBPT(3) calculations for calcium isotopes.

Demonstrated the polynomial equation approach for particle number constraint.

Extended the applicability of BMBPT to more accurate nuclear calculations.

Abstract

Bogoliubov many-body perturbation theory (BMBPT) relying on the breaking of U(1) global gauge symmetry has been recently formulated and applied to extend the applicability of standard perturbation theory to ab initio calculations of atomic nuclei away from shell closures. So far, practical applications have been limited to second-order calculations due to the lack of a generic algorithm to constrain the average particle number of the symmetry-broken state. This limitation is presently lifted and a general BMBPT formalism is presented that allows to constrain the particle-number expectation value at arbitrary order P. The constraint can be incorporated in closed form by solving a polynomial equation of degree P-1. The numerical procedure is illustrated through BMBPT(3) calculations of calcium isotopes using a nuclear Hamiltonian derived within chiral effective field theory.