Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2105.03935 · Physics Letters B

Low-rank matrix decompositions for ab initio nuclear structure

A. Tichai, P. Arthuis, K. Hebeler, M. Heinz, J. Hoppe, A. Schwenk

Published Sep 6, 2021• 11 citations
arXivHTMLPDF

Abstract

The extension of ab initio quantum many-body theory to higher accuracy and larger systems is intrinsically limited by the handling of large data objects in form of wave-function expansions and/or many-body operators. In this work we present matrix factorization techniques as a systematically improvable and robust tool to significantly reduce the computational cost in many-body applications at the price of introducing a moderate decomposition error. We demonstrate the power of this approach for the nuclear two-body systems, for many-body perturbation theory calculations of symmetric nuclear matter, and for non-perturbative in-medium similarity renormalization group simulations of finite nuclei. Establishing low-rank expansions of chiral nuclear interactions offers possibilities to reformulate many-body methods in ways that take advantage of tensor factorization strategies.

DOI
10.1016/j.physletb.2021.136623

Related papers

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

A. Tichai, Roman Schutski, Gustavo E. Scuseria +1 more

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 $^{4}\mathrm{He},^{16}\mathrm{O}$, and $^{40}\mathrm{Ca}$.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.

arXiv:1810.08419•Physical Review C•Mar 18, 2019•22 citations

Dimensionality reduction through tensor factorization: Application to ab initio nuclear physics calculations

Mikaël Frosini, T. Duguet, Pierre Tamagno +1 more

The construction of predictive models of atomic nuclei from first principles is a challenging (yet necessary) task towards the systematic generation of theoretical predictions (and associated uncertainties) to support nuclear data evaluation. The consistent description of the rich phenomenology of nuclear systems indeed requires the introduction of reductionist approaches that construct nuclei directly from interacting nucleons by solving the associated quantum many-body problem. In this context, so-called ab initio methods offer a promising route by deriving controlled (and systematically improvable) approximations both to the inter-nucleon interaction and to the solutions of the many-body problem. From a technical point of view, approximately solving the many-body Schrödinger equation in heavy open-shell systems typically requires the construction and contraction of large mode-4 (mode-6) tensors that need to be stored repeatedly. Recently, a new dimensionality reduction method based on randomized singular value decomposition has been introduced to reduce the numerical cost of many-body perturbation theory. This work applies this lightweight formalism to the study of the Germanium isotopic chain, where standard approaches would be too expansive to run. Inclusion of triaxiality is found to improve the overall agreement with experimental data on differential quantities.

arXiv:2407.17971•EPJ Web of Conferences•Jan 1, 2024•1 citation

Tensor factorization in ab initio many-body calculations: Triaxially-deformed (B) MBPT calculations in large bases

Mikaël Frosini, T. Duguet, Pierre Tamagno

Whether for fundamental studies or nuclear data evaluations, first-principle calculations of atomic nuclei constitute the path forward. Today, performing \textit{ab initio} calculations (a) of heavy nuclei, (b) of doubly open-shell nuclei or (c) with a sub-percent accuracy is at the forefront of nuclear structure theory. While combining any two of these features constitutes a major challenge, addressing the three at the same time is currently impossible. From a numerical standpoint, these challenges relate to the necessity to handle (i) very large single bases and (ii) mode-6, \textit{i.e.} three-body, tensors (iii) that must be stored repeatedly. Performing second-order many-body perturbation theory(ies) calculations based on triaxially deformed and superfluid reference states of doubly open-shell nuclei up to mass $A=72$, the present work achieves a significant step forward by addressing challenge (i). To do so, the memory and computational cost associated with the handling of large tensors is scaled down via the use of tensor factorization techniques. The presently used factorization format is based on a randomized singular value decomposition that does not require the computation and storage of the very large initial tensor. The procedure delivers an inexpensive and controllable approximation to the original problem, as presently illustrated for calculations that could not be performed without tensor factorization. With the presently developed technology at hand, one can envision to perform calculations of yet heavier doubly open-shell nuclei at sub-percent accuracy in a foreseeable future.

arXiv:2404.08532•Apr 12, 2024•0 citations

Pre-processing the nuclear many-body problem: Importance truncation versus tensor factorization techniques

A. Tichai, Ripoche, Julien, T. Duguet

The solution of the nuclear A-body problem encounters severe limitations from the size of many-body operators. These limitations are typically related to both the (iterative) storing of the associated tensors and to the computational time related to their multiple contractions in the calculation of various quantities of interest. However, not all the degrees of freedom encapsulated into these tensors equally contribute to the description of many-body observables. Identifying systematic and dominating patterns, a relevant objective is to achieve an \emph{a priori} reduction to the most relevant degrees of freedom via a pre-processing of the A-body problem. The present paper is dedicated to the analysis of two different paradigms to do so. The factorization of tensors in terms of lower-rank ones, whose know-how has been recently transferred to the realm of nuclear structure, is compared to a reduction of the tensors' index size based on an importance truncation. While the objective is to eventually utilize these pre-processing tools in the context of non-perturbative many-body methods, benchmark calculations are presently performed within the frame of perturbation theory. More specifically, we employ the recently introduced Bogoliubov many-body perturbation theory that is systematically applicable to open-shell nuclei displaying strong correlations. This extended perturbation theory serves as a jumpstart for non-perturbative Bogoliubov coupled cluster and Gorkov self-consistent Green's function theories. Results obtained in "small" model spaces are equally encouraging for tensor factorization and importance truncation techniques. While the former requires significant numerical developments to be applied in large model spaces, the latter is presently applied in this context and demonstrates great potential to enable high-accuracy calculations at a much reduced computational cost.

arXiv:1902.09043•Feb 25, 2019•0 citations

Natural orbitals for many-body expansion methods

J. Hoppe, A. Tichai, M. Heinz +2 more

The nuclear many-body problem for medium-mass systems is commonly addressed using wave-function expansion methods that build upon a second-quantized representation of many-body operators with respect to a chosen computational basis. While various options for the computational basis are available, perturbatively constructed natural orbitals recently have been shown to lead to significant improvement in many-body applications yielding faster model-space convergence and lower sensitivity to basis set parameters in large-scale no-core shell model diagonalizations. This work provides a detailed comparison of single-particle basis sets and a systematic benchmark of natural orbitals in nonperturbative many-body calculations using the in-medium similarity renormalization group approach. As a key outcome we find that the construction of natural orbitals in a large single-particle basis enables for performing the many-body calculation in a reduced space of much lower dimension, thus offering significant computational savings in practice that help extend the reach of ab initio methods towards heavier masses and higher accuracy.

arXiv:2009.04701•Physical Review C•Jan 28, 2021•35 citations

Singular value decomposition and similarity renormalization group evolution of nuclear interactions

Bin Zhu, R. Wirth, H. Hergert

One of the main challenges for ab initio nuclear many-body theory in the coming decade is the growth of computational and storage costs as calculations are extended to increasingly heavy, exotic, and structurally complex nuclei. Here, we investigate the factorization of nuclear interactions as a means to address this issue. We perform singular value decompositions (SVD) of current nucleon-nucleon interactions in partial wave representation and study the dependence of the singular value spectrum on interaction characteristics like regularization schemes and resolution scales. Next, we develop and implement the similarity renormalization group (SRG) evolution of the interaction in terms of the relevant singular vectors, and demonstrate that this SVD-SRG approach accurately preserves two-nucleon observables. We find that low-resolution interactions naturally allow the truncation of the SVD at low rank, and that a small number of relevant components is sufficient to capture the nuclear interaction and perform an accurate SRG evolution, while the Coulomb interaction requires special consideration. The rank is uniform across all partial waves and almost independent of the basis choice in the tested cases. This suggests an interpretation of the relevant singular components as mere representations of a small set of abstract operators that can describe the interaction and its SRG flow. Following the traditional workflow for nuclear interactions, we discuss how the transformation between the center-of-mass and laboratory frames creates redundant copies of the partial wave components when implemented in matrix representation, and we discuss strategies for mitigation. Finally, we test the low-rank approximation to the SRG-evolved interactions in many-body calculations using the in-medium SRG. By including nuclear radii in our analysis, we verify that the implementation of the SRG using the singular vectors of the interaction does not spoil the evolution of other observables.

arXiv:2106.01302•Physical Review C•Oct 11, 2021•9 citations
The extension of ab initio quantum many-body theory to higher accuracy and larger systems is intrinsically limited by the handling of large data objects in form of wave-function expansions and/or many-body operators. In this work we present matrix factorization techniques as a systematically improvable and robust tool to significantly reduce the computational cost in many-body applications at the price of introducing a moderate decomposition error. We demonstrate the power of this approach for the nuclear two-body systems, for many-body perturbation theory calculations of symmetric nuclear matter, and for non-perturbative in-medium similarity renormalization group simulations of finite nuclei. Establishing low-rank expansions of chiral nuclear interactions offers possibilities to reformulate many-body methods in ways that take advantage of tensor factorization strategies.
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 JJJ 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 4He,16O^{4}\mathrm{He},^{16}\mathrm{O}4He,16O, and 40Ca^{40}\mathrm{Ca}40Ca.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.
The construction of predictive models of atomic nuclei from first principles is a challenging (yet necessary) task towards the systematic generation of theoretical predictions (and associated uncertainties) to support nuclear data evaluation. The consistent description of the rich phenomenology of nuclear systems indeed requires the introduction of reductionist approaches that construct nuclei directly from interacting nucleons by solving the associated quantum many-body problem. In this context, so-called ab initio methods offer a promising route by deriving controlled (and systematically improvable) approximations both to the inter-nucleon interaction and to the solutions of the many-body problem. From a technical point of view, approximately solving the many-body Schrödinger equation in heavy open-shell systems typically requires the construction and contraction of large mode-4 (mode-6) tensors that need to be stored repeatedly. Recently, a new dimensionality reduction method based on randomized singular value decomposition has been introduced to reduce the numerical cost of many-body perturbation theory. This work applies this lightweight formalism to the study of the Germanium isotopic chain, where standard approaches would be too expansive to run. Inclusion of triaxiality is found to improve the overall agreement with experimental data on differential quantities.
Whether for fundamental studies or nuclear data evaluations, first-principle calculations of atomic nuclei constitute the path forward. Today, performing \textit{ab initio} calculations (a) of heavy nuclei, (b) of doubly open-shell nuclei or (c) with a sub-percent accuracy is at the forefront of nuclear structure theory. While combining any two of these features constitutes a major challenge, addressing the three at the same time is currently impossible. From a numerical standpoint, these challenges relate to the necessity to handle (i) very large single bases and (ii) mode-6, \textit{i.e.} three-body, tensors (iii) that must be stored repeatedly. Performing second-order many-body perturbation theory(ies) calculations based on triaxially deformed and superfluid reference states of doubly open-shell nuclei up to mass A=72A=72A=72, the present work achieves a significant step forward by addressing challenge (i). To do so, the memory and computational cost associated with the handling of large tensors is scaled down via the use of tensor factorization techniques. The presently used factorization format is based on a randomized singular value decomposition that does not require the computation and storage of the very large initial tensor. The procedure delivers an inexpensive and controllable approximation to the original problem, as presently illustrated for calculations that could not be performed without tensor factorization. With the presently developed technology at hand, one can envision to perform calculations of yet heavier doubly open-shell nuclei at sub-percent accuracy in a foreseeable future.
The solution of the nuclear A-body problem encounters severe limitations from the size of many-body operators. These limitations are typically related to both the (iterative) storing of the associated tensors and to the computational time related to their multiple contractions in the calculation of various quantities of interest. However, not all the degrees of freedom encapsulated into these tensors equally contribute to the description of many-body observables. Identifying systematic and dominating patterns, a relevant objective is to achieve an \emph{a priori} reduction to the most relevant degrees of freedom via a pre-processing of the A-body problem. The present paper is dedicated to the analysis of two different paradigms to do so. The factorization of tensors in terms of lower-rank ones, whose know-how has been recently transferred to the realm of nuclear structure, is compared to a reduction of the tensors' index size based on an importance truncation. While the objective is to eventually utilize these pre-processing tools in the context of non-perturbative many-body methods, benchmark calculations are presently performed within the frame of perturbation theory. More specifically, we employ the recently introduced Bogoliubov many-body perturbation theory that is systematically applicable to open-shell nuclei displaying strong correlations. This extended perturbation theory serves as a jumpstart for non-perturbative Bogoliubov coupled cluster and Gorkov self-consistent Green's function theories. Results obtained in "small" model spaces are equally encouraging for tensor factorization and importance truncation techniques. While the former requires significant numerical developments to be applied in large model spaces, the latter is presently applied in this context and demonstrates great potential to enable high-accuracy calculations at a much reduced computational cost.
The nuclear many-body problem for medium-mass systems is commonly addressed using wave-function expansion methods that build upon a second-quantized representation of many-body operators with respect to a chosen computational basis. While various options for the computational basis are available, perturbatively constructed natural orbitals recently have been shown to lead to significant improvement in many-body applications yielding faster model-space convergence and lower sensitivity to basis set parameters in large-scale no-core shell model diagonalizations. This work provides a detailed comparison of single-particle basis sets and a systematic benchmark of natural orbitals in nonperturbative many-body calculations using the in-medium similarity renormalization group approach. As a key outcome we find that the construction of natural orbitals in a large single-particle basis enables for performing the many-body calculation in a reduced space of much lower dimension, thus offering significant computational savings in practice that help extend the reach of ab initio methods towards heavier masses and higher accuracy.
One of the main challenges for ab initio nuclear many-body theory in the coming decade is the growth of computational and storage costs as calculations are extended to increasingly heavy, exotic, and structurally complex nuclei. Here, we investigate the factorization of nuclear interactions as a means to address this issue. We perform singular value decompositions (SVD) of current nucleon-nucleon interactions in partial wave representation and study the dependence of the singular value spectrum on interaction characteristics like regularization schemes and resolution scales. Next, we develop and implement the similarity renormalization group (SRG) evolution of the interaction in terms of the relevant singular vectors, and demonstrate that this SVD-SRG approach accurately preserves two-nucleon observables. We find that low-resolution interactions naturally allow the truncation of the SVD at low rank, and that a small number of relevant components is sufficient to capture the nuclear interaction and perform an accurate SRG evolution, while the Coulomb interaction requires special consideration. The rank is uniform across all partial waves and almost independent of the basis choice in the tested cases. This suggests an interpretation of the relevant singular components as mere representations of a small set of abstract operators that can describe the interaction and its SRG flow. Following the traditional workflow for nuclear interactions, we discuss how the transformation between the center-of-mass and laboratory frames creates redundant copies of the partial wave components when implemented in matrix representation, and we discuss strategies for mitigation. Finally, we test the low-rank approximation to the SRG-evolved interactions in many-body calculations using the in-medium SRG. By including nuclear radii in our analysis, we verify that the implementation of the SRG using the singular vectors of the interaction does not spoil the evolution of other observables.