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

ACM/IEEE Conference on Supercomputing (SC)

Quantum mechanics---Ab initio nuclear structure determination

David J. Dean

Published Jan 1, 2006• 0 citations
Publisher

Abstract

Solving the nuclear many-body problem engenders significant theoretical and computational challenges. I will briefly discuss the science that goes into and comes out of the nuclear problem, and then turn to the numerical implementations of coupled-cluster theory for the problem. Specific to the nuclear problem is the presence of real three-body forces which are computationally represented by six-dimensional tensors. Both memory and communications become issues as one begins to scale up the underlying coupled-cluster algorithms to thousands of processors. Another interesting feature of the nuclear problem is weak binding that occurs in very unstable nuclei. When the single-particle continuum is present, the nuclear Hamiltonian becomes complex. I also will describe our attempts to solve the complex coupled-cluster problem in the presence of complex interactions. Efficient parallelism requires spreading the interaction matrices globally and the coupled-cluster amplitudes in a more local fashion.

DOI
10.1145/1188455.1188520

Related papers

Computational aspects of nuclear coupled-cluster theory

David J. Dean, G. Hagen, Morten Hjorth-Jensen +1 more

Coupled-cluster (CC) theory represents an important theoretical tool that we use to solve the quantum many-body problem. CC theory also lends itself to computation in a parallel computing environment. In this paper, we present selected results from ab initio studies of stable and weakly bound nuclei utilizing computational techniques that we employ to solve CC theory. We also outline several perspectives for future research directions in this area.

Computational Science & Discovery•Nov 20, 2008•6 citations

1Computational aspects of nuclear coupled-cluster theory

G. Hagen, Hai Ah Nam

We discuss computational aspects of the spherical coupled-cluster method specific to the nuclear many-body problem. Using chiral nucleon-nucleon interaction at next-to-next-to-next-to leading order (N3LO) with cutoff Λ = 500MeV, we present coupled-cluster results for the ground state of 40Ca. Scaling and performance studies are presented together with challenges we meet with when extending the coupled-cluster effort to nuclei mass hundred and beyond. The low-energy nuclear many-body problem is a challenging undertaking, how-ever, in the last decade there has been significant progress in first principle cal-culations of nuclei.1), 2), 3), 4), 5), 6), 7), 8), 9) With the advance of chiral effective field theory10), 11), 12), 13), 14), 15) which allows for a systematic and consistent derivation of

•Jan 1, 2012•3 citations

Computational Aspects of Nuclear Coupled-Cluster Theory

Gaute Hagen, Hai Ah Nam

We discuss computational aspects of the spherical coupled-cluster method specific to the nuclear many-body problem. Using chiral nucleon-nucleon interaction at next-to-next-to-next-to leading order (N3LO) with cutoff Λ= 500 MeV, we present coupled-cluster results for the ground state of 40Ca. Scaling and performance studies are presented together with challenges we meet with when extending the coupled-cluster effort to nuclei mass hundred and beyond.

arXiv:1203.3765•Progress of Theoretical Physics Supplement•Jan 1, 2012•6 citations

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

Many-Body Perturbation Theory for Ab Initio Nuclear Structure

A. Tichai

The solution of the quantum many-body problem for medium-mass nuclei using realistic nuclear interactions poses a superbe challenge for nuclear structure research. Because an exact solution can only be provided for the lightest nuclei, one has to rely on approximate solutions when proceeding to heavier systems. Over the past years, tremendous progress has been made in the development and application of systematically improvable expansion methods and an accurate description of nuclear observables has become viable up to mass number $A \approx 100$. While closed-shell systems are consistently described via a plethora of different many-body methods, the extension to genuine open-shell systems still remains a major challenge and up to now there is no \textit{ab initio} many-body method which applies equally well to systems with even and odd mass numbers. The goal of this thesis is the development and implementation of innovative perturbative approaches with genuine open-shell capabilities. This requires the extension of well-known single-reference approaches to more general vacua. In this work we choose two complementary routes for the usage of generalized reference states. First, we derive a new \emph{ab initio} approach based on multi-configurational reference states that are conveniently derived from a prior no-core shell model calculation. Perturbative corrections are derived via second-order many-body perturbation theory, thus, merging configuration interaction and many-body perturbation theory. The generality of this ansatz enables for a treatment of medium-mass systems with arbitrary mass number, as well as the extension to low-lying excited states such that ground and excited states are treated on an equal footing. In a complementary approach, we use reference states that break a symmetry of the underlying Hamiltonian. In the simplest case this corresponds to the expansion around a particle-number-broken Hartree-Fock-Bogoliubov vacuum which is obtained from a mean-field calculation. Pairing correlations are already absorbed in the reference state. The mild scaling allows for investigations up to tin isotopic chains. All benchmarks were performed using state-of-the-art chiral two- plus three-nucleon interactions thus allowing for a universal description of nuclear observables in the medium-mass regime.

•Jan 1, 2017•1 citation

Wavefunction matching for solving quantum many-body problems

Serdar Elhatisari, Lukas Bovermann, Yuan-Zhuo Ma +15 more

Abstract Ab initio calculations have an essential role in our fundamental understanding of quantum many-body systems across many subfields, from strongly correlated fermions 1–3 to quantum chemistry 4–6 and from atomic and molecular systems 7–9 to nuclear physics 10–14 . One of the primary challenges is to perform accurate calculations for systems where the interactions may be complicated and difficult for the chosen computational method to handle. Here we address the problem by introducing an approach called wavefunction matching. Wavefunction matching transforms the interaction between particles so that the wavefunctions up to some finite range match that of an easily computable interaction. This allows for calculations of systems that would otherwise be impossible owing to problems such as Monte Carlo sign cancellations. We apply the method to lattice Monte Carlo simulations 15,16 of light nuclei, medium-mass nuclei, neutron matter and nuclear matter. We use high-fidelity chiral effective field theory interactions 17,18 and find good agreement with empirical data. These results are accompanied by insights on the nuclear interactions that may help to resolve long-standing challenges in accurately reproducing nuclear binding energies, charge radii and nuclear-matter saturation in ab initio calculations 19,20 .

Nature•May 15, 2024•62 citations
Solving the nuclear many-body problem engenders significant theoretical and computational challenges. I will briefly discuss the science that goes into and comes out of the nuclear problem, and then turn to the numerical implementations of coupled-cluster theory for the problem. Specific to the nuclear problem is the presence of real three-body forces which are computationally represented by six-dimensional tensors. Both memory and communications become issues as one begins to scale up the underlying coupled-cluster algorithms to thousands of processors. Another interesting feature of the nuclear problem is weak binding that occurs in very unstable nuclei. When the single-particle continuum is present, the nuclear Hamiltonian becomes complex. I also will describe our attempts to solve the complex coupled-cluster problem in the presence of complex interactions. Efficient parallelism requires spreading the interaction matrices globally and the coupled-cluster amplitudes in a more local fashion.
We discuss computational aspects of the spherical coupled-cluster method specific to the nuclear many-body problem. Using chiral nucleon-nucleon interaction at next-to-next-to-next-to leading order (N3LO) with cutoff Λ = 500MeV, we present coupled-cluster results for the ground state of 40Ca. Scaling and performance studies are presented together with challenges we meet with when extending the coupled-cluster effort to nuclei mass hundred and beyond. The low-energy nuclear many-body problem is a challenging undertaking, how-ever, in the last decade there has been significant progress in first principle cal-culations of nuclei.1), 2), 3), 4), 5), 6), 7), 8), 9) With the advance of chiral effective field theory10), 11), 12), 13), 14), 15) which allows for a systematic and consistent derivation of
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 solution of the quantum many-body problem for medium-mass nuclei using realistic nuclear interactions poses a superbe challenge for nuclear structure research. Because an exact solution can only be provided for the lightest nuclei, one has to rely on approximate solutions when proceeding to heavier systems. Over the past years, tremendous progress has been made in the development and application of systematically improvable expansion methods and an accurate description of nuclear observables has become viable up to mass number A≈100A \approx 100A≈100. While closed-shell systems are consistently described via a plethora of different many-body methods, the extension to genuine open-shell systems still remains a major challenge and up to now there is no \textit{ab initio} many-body method which applies equally well to systems with even and odd mass numbers. The goal of this thesis is the development and implementation of innovative perturbative approaches with genuine open-shell capabilities. This requires the extension of well-known single-reference approaches to more general vacua. In this work we choose two complementary routes for the usage of generalized reference states. First, we derive a new \emph{ab initio} approach based on multi-configurational reference states that are conveniently derived from a prior no-core shell model calculation. Perturbative corrections are derived via second-order many-body perturbation theory, thus, merging configuration interaction and many-body perturbation theory. The generality of this ansatz enables for a treatment of medium-mass systems with arbitrary mass number, as well as the extension to low-lying excited states such that ground and excited states are treated on an equal footing. In a complementary approach, we use reference states that break a symmetry of the underlying Hamiltonian. In the simplest case this corresponds to the expansion around a particle-number-broken Hartree-Fock-Bogoliubov vacuum which is obtained from a mean-field calculation. Pairing correlations are already absorbed in the reference state. The mild scaling allows for investigations up to tin isotopic chains. All benchmarks were performed using state-of-the-art chiral two- plus three-nucleon interactions thus allowing for a universal description of nuclear observables in the medium-mass regime.
Abstract Ab initio calculations have an essential role in our fundamental understanding of quantum many-body systems across many subfields, from strongly correlated fermions 1–3 to quantum chemistry 4–6 and from atomic and molecular systems 7–9 to nuclear physics 10–14 . One of the primary challenges is to perform accurate calculations for systems where the interactions may be complicated and difficult for the chosen computational method to handle. Here we address the problem by introducing an approach called wavefunction matching. Wavefunction matching transforms the interaction between particles so that the wavefunctions up to some finite range match that of an easily computable interaction. This allows for calculations of systems that would otherwise be impossible owing to problems such as Monte Carlo sign cancellations. We apply the method to lattice Monte Carlo simulations 15,16 of light nuclei, medium-mass nuclei, neutron matter and nuclear matter. We use high-fidelity chiral effective field theory interactions 17,18 and find good agreement with empirical data. These results are accompanied by insights on the nuclear interactions that may help to resolve long-standing challenges in accurately reproducing nuclear binding energies, charge radii and nuclear-matter saturation in ab initio calculations 19,20 .