Tensor-decomposition techniques for ab initio nuclear structure calculations: From chiral nuclear potentials to ground-state energies

A. Tichai, R. Schutski, G. E. Scuseria, and T. Duguet
Phys. Rev. C 99, 034320 – Published 18 March 2019

Abstract

Background: The computational resources needed to generate the ab initio solution of the nuclear many-body problem for increasing mass number and/or accuracy necessitates innovative developments to improve upon (i) the storage of many-body operators and (ii) the scaling of many-body methods used to evaluate nuclear observables. The storing and efficient handling of many-body operators with high particle ranks is currently one of the major bottlenecks limiting the applicability range of ab initio studies with respect to mass number and accuracy. Recently, the application of tensor decomposition techniques to many-body tensors has proven highly beneficial to reduce the computational cost of ab initio calculations in quantum chemistry and solid-state physics.

Purpose: The impact of applying state-of-the-art tensor factorization techniques to modern nuclear Hamiltonians derived from chiral effective field theory is investigated. Subsequently, the error induced by the tensor decomposition of the input Hamiltonian on ground-state energies of closed-shell nuclei calculated via second-order many-body perturbation theory is benchmarked.

Methods: The first proof-of-principles application of tensor-decomposition techniques to the nuclear Hamiltonian is performed. Two different tensor formats are investigated by systematically benchmarking the approximation error on matrix elements stored in various bases of interest. The analysis is achieved while including normal-ordered three-nucleon interactions that are nowadays used as input to the most advanced ab initio calculations in medium-mass nuclei. With the aid of the factorized Hamiltonian, the second-order perturbative correction to ground-state energies is decomposed and the scaling properties of the underlying tensor network are discussed.

Results: The employed tensor formats are found to lead to an efficient data compression of two-body matrix elements of the nuclear Hamiltonian. In particular, the sophisticated tensor hypercontraction scheme yields low tensor ranks with respect to both harmonic-oscillator and Hartree-Fock single-particle bases. It is found that the tensor rank depends on the two-body total angular momentum J for which one performs the decomposition, which is itself directly related to the sparsity of the corresponding tensor. Furthermore, including normal-ordered two-body contributions originating from three-body interactions does not compromise the efficient data compression. Ultimately, the use of factorized matrix elements authorizes controlled approximations of the exact second-order ground-state energy corrections. In particular, a small enough error is obtained from low-rank factorizations in He4,O16, and Ca40.

Conclusions: It is presently demonstrated that tensor-decomposition techniques can be efficiently applied to systematically approximate the nuclear many-body Hamiltonian in terms of lower-rank tensors. Beyond the input Hamiltonian, tensor-decomposition techniques can be envisioned to scale down the cost of state-of-the-art nonperturbative many-body methods in order to extend ab initio studies to (i) higher precisions, (ii) larger masses, and (iii) nuclei of doubly open-shell character.

  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
6 More
  • Received 19 October 2018

DOI:https://doi.org/10.1103/PhysRevC.99.034320

©2019 American Physical Society

Physics Subject Headings (PhySH)

Nuclear Physics

Authors & Affiliations

A. Tichai1,*, R. Schutski2,3,†, G. E. Scuseria2,4,‡, and T. Duguet5,6,§

  • 1ESNT, CEA-Saclay, DRF, IRFU, Département de Physique Nucléaire, Université Paris-Saclay, F-91191 Gif-sur-Yvette, France
  • 2Department of Chemistry, Rice University, Houston, Texas 77005-1892, USA
  • 3CDISE, Skolkovo Institute of Science and Technology, 3 Nobel St., Moscow, Russia
  • 4Department of Physics and Astronomy, Rice University, Houston, Texas 77005-1892, USA
  • 5IRFU, CEA, Université Paris-Saclay, 91191 Gif-sur-Yvette, France
  • 6KU Leuven, Instituut voor Kern- en Stralingsfysica, 3001 Leuven, Belgium

  • *alexander.tichai@cea.fr
  • r.schutski@skoltech.ru
  • guscus@rice.edu
  • §thomas.duguet@cea.fr

Article Text (Subscription Required)

Click to Expand

References (Subscription Required)

Click to Expand
Issue

Vol. 99, Iss. 3 — March 2019

Reuse & Permissions
Access Options
CHORUS

Article Available via CHORUS

Download Accepted Manuscript
Author publication services for translation and copyediting assistance advertisement

Authorization Required


×
×

Images

×

Sign up to receive regular email alerts from Physical Review C

Log In

Cancel
×

Search


Article Lookup

Paste a citation or DOI

Enter a citation
×