Pepijn Demol

In order to solve the A-body Schr\"odinger equation both accurately and efficiently for open-shell nuclei, a novel many-body method coined as Bogoliubov many-body perturbation theory (BMBPT) was recently formalized and applied at low orders. Based on the breaking of U(1) symmetry associated with particle-number conservation, this perturbation theory must operate under the constraint that the average number of particles is self-consistently adjusted at each perturbative order. The corresponding formalism is presently detailed with the goal to characterize the behavior of the associated Taylor series. BMBPT is, thus, investigated numerically up to high orders at the price of restricting oneself to a small, i.e. schematic, portion of Fock space. While low-order results only differ by 2 - 3 % from those obtained via a configuration interaction (CI) diagonalization, the series is shown to eventually diverge. The application of a novel resummation method coined as eigenvector continuation further increase the accuracy when built from low-order BMBPT corrections and quickly converges towards the CI result when applied at higher orders. Furthermore, the numerically-costly self-consistent particle number adjustment procedure is shown to be safely bypassed via the use of a computationally cheap a posteriori correction method. Eventually, the present work validates the fact that low order BMBPT calculations based on an a posteriori (average) particle number correction deliver controlled results and demonstrates that they can be optimally complemented by the eigenvector continuation method to provide results with sub-percent accuracy. This approach is, thus, planned to become a workhorse for realistic ab initio calculations of open-shell nuclei in the near future.

Quantum many-body theory has witnessed tremendous progress in various fields, ranging from atomic and solid-state physics to quantum chemistry and nuclear structure. Due to the inherent computational burden linked to the ab initio treatment of microscopic fermionic systems, it is desirable to obtain accurate results through low-order perturbation theory. In atomic nuclei, however, effects such as strong short-range repulsion between nucleons can spoil the convergence of the expansion and make the reliability of perturbation theory unclear. Mathematicians have devised an extensive machinery to overcome the problem of divergent expansions by making use of so-called resummation methods. In large-scale many-body applications, such schemes are often of limited use since no a priori analytical knowledge of the expansion is available. We present here eigenvector continuation as an alternative resummation tool that is both efficient and reliable because it is based on robust and simple mathematical principles.