Wave-packet continuum discretisation for nucleon-nucleon scattering predictions

Abstract

In this paper we analyse the efficiency, precision, and accuracy of computing elastic nucleon-nucleon (NN) scattering amplitudes with the wave-packet continuum discretisation method (WPCD). This method provides approximate scattering solutions at multiple scattering energies simultaneously. We therefore utilise a graphics processing unit (GPU) to explore the benefits of this inherent parallelism. From a theoretical perspective, the WPCD method promises a speedup compared to a standard matrix-inversion method. We use the chiral NNLO${}_{\mathrm{opt}}$ interaction to demonstrate that WPCD enables efficient computation of NN scattering amplitudes provided one can tolerate an averaged method error of $1-5$ mb in the total cross section at scattering energies $0-350$ MeV in the laboratory frame of reference. Considering only scattering energies $\sim 40-350$ MeV, we find a smaller method error of $\lesssim 1-2$ mb. By increasing the number of wave-packets we can further reduce the overall method error. However, the parallel leverage of the WPCD method will be offset by the increased size of the resulting discretisation mesh. In practice, a GPU-implementation is mainly advantageous for matrices that fit in the fast on-chip shared memory. We find that WPCD is a promising method for computationally efficient, statistical analyses of nuclear interactions from effective field theory, where we can utilise Bayesian inference methods to incorporate relevant uncertainties.

arXiv comment33 pages, 9 figures, v2 corresponds to published version

References in corpus (21)

• Modern Theory of Nuclear Forces
• Chiral effective field theory and nuclear forces
• Accurate Charge-Dependent Nucleon-Nucleon Potential at Fourth Order of Chiral Perturbation Theory
• Effective Field Theory for Few-Nucleon Systems
• Nuclear effective field theory: status and perspectives
• An optimized chiral nucleon-nucleon interaction at next-to-next-to-leading order
• Semilocal momentum-space regularized chiral two-nucleon potentials up to fifth order
• Minimally non-local nucleon-nucleon potentials with chiral two-pion exchange including $\Delta$'s
• Bayesian truncation errors in chiral effective field theory: nucleon-nucleon observables
• Uncertainty analysis and order-by-order optimization of chiral nuclear interactions
• Eigenvector continuation with subspace learning
• Coarse grained Potential analysis of neutron-proton and proton-proton scattering below pion production threshold
• Bayesian parameter estimation for effective field theories
• Eigenvector Continuation as an Efficient and Accurate Emulator for Uncertainty Quantification
• Dissecting the NVidia Turing T4 GPU via Microbenchmarking
• Efficient emulators for scattering using eigenvector continuation
• Bound state techniques to solve the multiparticle scattering problem
• Fast & accurate emulation of two-body scattering observables without wave functions
• Fast emulation of quantum three-body scattering
• Discrete Wave-Packet Representation in Nuclear Matter Calculations
• New general approach in few-body scattering calculations: Solving discretized Faddeev equations on a graphics processing unit

Cited by in corpus (3)

• Nuclear Forces for Precision Nuclear Physics -- a collection of perspectives
• Posterior predictive distributions of neutron-deuteron cross sections
• Active learning emulators for nuclear two-body scattering in momentum space