Modewise Johnson-Lindenstrauss Embeddings for Nuclear Many-Body Theory

Abstract

In this work, we explore modewise Johnson-Lindenstrauss embeddings (JLEs) as a tool to reduce the computational cost and memory requirements of nuclear many-body methods. JLEs are randomized projections of high-dimensional data tensors onto low-dimensional subspaces that preserve key structural features. Such embeddings allow for the oblivious and incremental compression of large tensors, e.g., the nuclear Hamiltonian, into significantly smaller random sketches that still allow for the accurate calculation of ground-state energies and other observables. Their oblivious character 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., 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 analyze the JLE's impact on the second-order Many-Body Perturbation Theory (MBPT) corrections for nuclear ground-state observables. 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 thousand-fold 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.

23 pages, 14 figures

References in corpus (19)

• Improved nuclear matter calculations from chiral low-momentum interactions
• Sketching as a Tool for Numerical Linear Algebra
• A Guided Tour of Ab Initio Nuclear Many-Body Theory
• Ab initio predictions link the neutron skin of ${}^{208}$Pb to nuclear forces
• Coupled-cluster theory for three-body Hamiltonians
• Converged ab initio calculations of heavy nuclei
• Charge radii of exotic neon and magnesium isotopes
• In-Medium Similarity Renormalization Group for Closed and Open-Shell Nuclei
• Angular-momentum projection in coupled-cluster theory: structure of ${}^{34}$Mg
• Multi-reference many-body perturbation theory for nuclei III -- Ab initio calculations at second order in PGCM-PT
• Rank-reduced coupled-cluster III. Tensor hypercontraction of the doubles amplitudes
• Implementation of the full CCSDT electronic structure model with tensor decompositions
• Natural orbitals for many-body expansion methods
• Normal ordering of three-nucleon interactions for ab initio calculations of heavy nuclei
• Spectra of Open-Shell Nuclei with Padé-Resummed Degenerate Perturbation Theory
• Gorkov algebraic diagrammatic construction formalism at third order
• Low-rank matrix decompositions for ab initio nuclear structure
• Least-square approach for singular value decompositions of scattering problems
• Symmetry restoration methods

Cited by in corpus (2)

• Randomized Low-Rank Decompositions of Nuclear Three-Body Interactions
• In-Medium Similarity Renormalization Group at Finite Temperature