https://doi.org/10.1140/epja/s10050-023-00999-5
Regular Article –Theoretical Physics
Modewise Johnson–Lindenstrauss embeddings for nuclear many-body theory
1
Department of Computational Mathematics, Science and Engineering, Michigan State University, 48824, East Lansing, MI, USA
2
Facility for Rare Isotope Beams, Michigan State University, 48824, East Lansing, MI, USA
3
Department of Mathematics, Michigan State University, 48824, East Lansing, MI, USA
4
Department of Physics and Astronomy, Michigan State University, 48824, East Lansing, MI, USA
Received:
7
November
2022
Accepted:
6
April
2023
Published online:
2
May
2023
In the present work, we initiate a program that explores modewise Johnson–Lindenstrauss embeddings (JLEs) as a tool to reduce the computational cost and memory requirements of (nuclear) many-body methods. These embeddings are randomized projections of high-dimensional data tensors onto low-dimensional subspaces that preserve structural features like norms and inner products. An appealing feature of randomized embedding techniques is that they allow for the oblivious and incremental compression of large tensors, e.g., the nuclear Hamiltonian or wave functions amplitudes, into significantly smaller random sketches that still allow for the accurate calculation of ground-state energies and other observables. In particular, the oblivious nature of randomized JLE techniques makes it possible to compress a tensor without knowing in advance exactly what observables one might want to approximate at a later time. This opens the door for the use of tensors that are much too large to store in memory, e.g., untruncated three-nucleon forces in current approaches, or complete two- plus three-nucleon Hamiltonians in large, symmetry-unrestricted bases. Such compressed Hamiltonians can be stored and used later on with relative ease. As a first step, we perform a detailed analysis of a JLE’s impact on the second-order Many-Body Perturbation Theory (MBPT) corrections for nuclear ground-state observables like the energy and the radius, noting that these will be the dominant corrections in a well-behaved perturbative expansion, and highly important implicit contributions even in nonperturbative approaches. Numerical experiments for a wide range of closed-shell nuclei, model spaces and state-of-the-art nuclear interactions demonstrate the validity and potential of the proposed approach: We can compress nuclear Hamiltonians hundred- to thousandfold while only incurring mean relative errors of 1% or less in ground-state observables. Importantly, we show that JLEs capture the relevant physical information contained in the highly structured Hamiltonian tensor despite their random characteristics. In addition to the significant storage savings, the achieved compressions imply multiple order-of-magnitude reductions in computational effort when the compressed Hamiltonians are used in higher-order MBPT or nonperturbative many-body methods.
Copyright comment Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
© The Author(s), under exclusive licence to Società Italiana di Fisica and Springer-Verlag GmbH Germany, part of Springer Nature 2023. Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
