-
Illuminating Systematic Trends in Nuclear Data with Generative Machine Learning Models
Authors:
Jordan M. R. Fox,
Kyle A. Wendt
Abstract:
We introduce a novel method for studying systematic trends in nuclear reaction data using generative adversarial networks. Libraries of nuclear cross section evaluations exhibit intricate systematic trends across the nuclear landscape, and predictive models capable of reproducing and analyzing these trends are valuable for many applications. We have developed a predictive model using deep generati…
▽ More
We introduce a novel method for studying systematic trends in nuclear reaction data using generative adversarial networks. Libraries of nuclear cross section evaluations exhibit intricate systematic trends across the nuclear landscape, and predictive models capable of reproducing and analyzing these trends are valuable for many applications. We have developed a predictive model using deep generative adversarial networks to learn trends from the inelastic neutron scattering channel of TENDL for even-even nuclei. The system predicts cross sections based on adding/subtracting particles to/from the target nucleus. It can thus help identify cross sections that break from expected trends and predict beyond the limit of current experiments. Our model can produce good predictions for cross section curves for many nuclides, and it is most robust near the line of stability. We also create an ensemble of predictions to leverage different correlations and estimate model uncertainty. This research marks an important first step in computer generation of nuclear cross-section libraries.
△ Less
Submitted 29 April, 2024; v1 submitted 24 March, 2024;
originally announced March 2024.
-
Direct pulse-level compilation of arbitrary quantum logic gates on superconducting qutrits
Authors:
Yujin Cho,
Kristin M. Beck,
Alessandro R. Castelli,
Kyle A. Wendt,
Bram Evert,
Matthew J. Reagor,
Jonathan L DuBois
Abstract:
Advanced simulations and calculations on quantum computers require high-fidelity implementations of quantum operations. The universal gateset approach builds complex unitaries from a small set of primitive gates, often resulting in a long gate sequence which is typically a leading factor in the total accumulated error. Compiling a complex unitary for processors with higher-dimensional logical elem…
▽ More
Advanced simulations and calculations on quantum computers require high-fidelity implementations of quantum operations. The universal gateset approach builds complex unitaries from a small set of primitive gates, often resulting in a long gate sequence which is typically a leading factor in the total accumulated error. Compiling a complex unitary for processors with higher-dimensional logical elements, such as qutrits, exacerbates the accumulated error per unitary, since an even longer gate sequence is required. Optimal control methods promise time and resource efficient compact gate sequences and, therefore, higher fidelity. These methods generate pulses that can directly implement any complex unitary on a quantum device. In this work, we demonstrate any arbitrary qubit and qutrit gate can be realized with high-fidelity, which can significantly reduce the length of a gate sequence. We generate and test pulses for a large set of randomly selected arbitrary unitaries on several quantum processing units (QPUs): the LLNL Quantum Device and Integration Testbed (QuDIT) standard QPU and three of Rigetti QPUs: Ankaa-2, Ankaa-9Q-1, and Aspen-M-3. On the QuDIT platform's standard QPU, the average fidelity of random qutrit gates is 97.9+-0.5% measured with conventional QPT and 98.8+-0.6% from QPT with gate folding. Rigetti's Ankaa-2 achieves random qubit gates with an average fidelity of 98.4+-0.5% (conventional QPT) and 99.7+-0.1% (QPT with gate folding). On Ankaa-9Q-1 and Aspen-M-3, the average fidelities with conventional qubit QPT measurements were higher than 99%. We show that optimal control gates are robust to drift for at least three hours and that the same calibration parameters can be used for all implemented gates. Our work promises the calibration overheads for optimal control gates can be made small enough to enable efficient quantum circuits based on this technique.
△ Less
Submitted 28 June, 2024; v1 submitted 7 March, 2023;
originally announced March 2023.
-
A quantum-classical co-processing protocol towards simulating nuclear reactions on contemporary quantum hardware
Authors:
Francesco Turro,
Trevor Chistolini,
Akel Hashim,
Yosep Kim,
William Livingston,
Kyle. A. Wendt,
Jonathan L Dubois,
Francesco Pederiva,
Sofia Quaglioni,
David I. Santiago,
Irfan Siddiqi
Abstract:
Quantum computers hold great promise for arriving at exact simulations of nuclear dynamical processes (e.g., scattering and reactions) that are paramount to the study of nuclear matter at the limit of stability and to explaining the formation of chemical elements in stars. However, quantum simulations of the unitary (real) time dynamics of fermionic many-body systems require a currently prohibitiv…
▽ More
Quantum computers hold great promise for arriving at exact simulations of nuclear dynamical processes (e.g., scattering and reactions) that are paramount to the study of nuclear matter at the limit of stability and to explaining the formation of chemical elements in stars. However, quantum simulations of the unitary (real) time dynamics of fermionic many-body systems require a currently prohibitive number of reliable and long-lived qubits. We propose a co-processing algorithm for the simulation of real-time dynamics in which the time evolution of the spatial coordinates is carried out on a classical processor, while the evolution of the spin degrees of freedom is carried out on a quantum processor. This hybrid algorithm is demonstrated by a quantum simulation of the scattering of two neutrons performed at the Lawrence Berkeley National Laboratory's Advanced Quantum Testbed. We show that, after implementation of error mitigation strategies to improve the accuracy of the algorithm in addition to the use of either circuit compression techniques or tomography as methods to elucidate the onset of decoherence, this initial demonstration validates the principle of the proposed co-processing scheme. We anticipate that a generalization of this present scheme will open the way for (real-time) path integral simulations of nuclear scattering.
△ Less
Submitted 26 September, 2023; v1 submitted 13 February, 2023;
originally announced February 2023.
-
Towards Precise and Accurate Calculations of Neutrinoless Double-Beta Decay: Project Scoping Workshop Report
Authors:
V. Cirigliano,
Z. Davoudi,
J. Engel,
R. J. Furnstahl,
G. Hagen,
U. Heinz,
H. Hergert,
M. Horoi,
C. W. Johnson,
A. Lovato,
E. Mereghetti,
W. Nazarewicz,
A. Nicholson,
T. Papenbrock,
S. Pastore,
M. Plumlee,
D. R. Phillips,
P. E. Shanahan,
S. R. Stroberg,
F. Viens,
A. Walker-Loud,
K. A. Wendt,
S. M. Wild
Abstract:
We present the results of a National Science Foundation (NSF) Project Scoping Workshop, the purpose of which was to assess the current status of calculations for the nuclear matrix elements governing neutrinoless double-beta decay and determine if more work on them is required. After reviewing important recent progress in the application of effective field theory, lattice quantum chromodynamics, a…
▽ More
We present the results of a National Science Foundation (NSF) Project Scoping Workshop, the purpose of which was to assess the current status of calculations for the nuclear matrix elements governing neutrinoless double-beta decay and determine if more work on them is required. After reviewing important recent progress in the application of effective field theory, lattice quantum chromodynamics, and ab initio nuclear-structure theory to double-beta decay, we discuss the state of the art in nuclear-physics uncertainty quantification and then construct a road map for work in all these areas to fully complement the increasingly sensitive experiments in operation and under development. The road map contains specific projects in theoretical and computational physics as well as an uncertainty-quantification plan that employs Bayesian Model Mixing and an analysis of correlations between double-beta-decay rates and other observables. The goal of this program is a set of accurate and precise matrix elements, in all nuclei of interest to experimentalists, delivered together with carefully assessed uncertainties. Such calculations will allow crisp conclusions from the observation or non-observation of neutrinoless double-beta decay, no matter what new physics is at play.
△ Less
Submitted 3 July, 2022;
originally announced July 2022.
-
Quantum state preparation by adiabatic evolution with customized gates
Authors:
E. A. Coello Perez,
J. Bonitati,
D. Lee,
S. Quaglioni,
K. A. Wendt
Abstract:
Quantum state preparation by adiabatic evolution is currently rendered ineffective by the long implementation times of the underlying quantum circuits, comparable to the decoherence time of present and near-term quantum devices. These implementation times can be significantly reduced by realizing the evolution with a minimal number of customized gates. Employing a realistic model of a two-qubit pr…
▽ More
Quantum state preparation by adiabatic evolution is currently rendered ineffective by the long implementation times of the underlying quantum circuits, comparable to the decoherence time of present and near-term quantum devices. These implementation times can be significantly reduced by realizing the evolution with a minimal number of customized gates. Employing a realistic model of a two-qubit processor, we carried out classical device-level simulations of the adiabatic evolution of a two-spin system implemented with customized two-qubit gates. These device-level simulations were compared with (experimental) ones solving the same problem on IBMQ systems. When used to emulate the IBMQ quantum circuit, our device-level simulations reached state fidelities ranging from 65% to 85%, similar to the actual performance of a diverse set of IBMQ devices.When we reduced the implementation time by using a minimal number of customized gates, however, the loss of fidelity was reduced by at least a factor of four, allowing us to accurately extract the energy of the target state. This improvement is enough to render adiabatic evolution useful for quantum state preparation for small systems or as a preconditioner for other state preparation methods.
△ Less
Submitted 23 November, 2021;
originally announced November 2021.
-
Imaginary Time Propagation on a Quantum Chip
Authors:
Francesco Turro,
Alessandro Roggero,
Valentina Amitrano,
Piero Luchi,
Kyle A. Wendt,
Jonathan L DuBois,
Sofia Quaglioni,
Francesco Pederiva
Abstract:
Evolution in imaginary time is a prominent technique for finding the ground state of quantum many-body systems, and the heart of a number of numerical methods that have been used with great success in quantum chemistry, condensed matter and nuclear physics. We propose an algorithm to implement imaginary time propagation on a quantum computer. Our algorithm is devised in the context of an efficient…
▽ More
Evolution in imaginary time is a prominent technique for finding the ground state of quantum many-body systems, and the heart of a number of numerical methods that have been used with great success in quantum chemistry, condensed matter and nuclear physics. We propose an algorithm to implement imaginary time propagation on a quantum computer. Our algorithm is devised in the context of an efficient encoding into an optimized gate, drawing on the underlying characteristics of the quantum device, of a unitary operation in an extended Hilbert space. However, we proved that for simple problems it can be successfully applied to standard digital quantum machines. This work paves the way for porting quantum many-body methods based on imaginary-time propagation to near-term quantum devices, enabling the future quantum simulation of the ground states of a broad class of microscopic systems.
△ Less
Submitted 12 November, 2021; v1 submitted 24 February, 2021;
originally announced February 2021.
-
Quantifying uncertainties in neutron-alpha scattering with chiral nucleon-nucleon and three-nucleon forces
Authors:
Konstantinos Kravvaris,
Kevin R. Quinlan,
Sofia Quaglioni,
Kyle A. Wendt,
Petr Navratil
Abstract:
Background: Modern ab initio theory combined with high-quality nucleon-nucleon (NN) and three-nucleon (3N) interactions from chiral effective field theory (EFT) can provide a predictive description of low-energy light-nuclei reactions relevant for astrophysics and fusion-energy applications. However, the high cost of computations has so far impeded a complete analysis of the uncertainty budget of…
▽ More
Background: Modern ab initio theory combined with high-quality nucleon-nucleon (NN) and three-nucleon (3N) interactions from chiral effective field theory (EFT) can provide a predictive description of low-energy light-nuclei reactions relevant for astrophysics and fusion-energy applications. However, the high cost of computations has so far impeded a complete analysis of the uncertainty budget of such calculations. Purpose: Starting from NN potentials up to fifth order (N4LO) combined with leading-order 3N forces, we study how the order-by-order convergence of the chiral expansion and confidence intervals for the 3N contact and contact-plus-one-pion-exchange low-energy constants (cE and cD) contribute to the overall uncertainty budget of many-body calculations of neutron-He elastic scattering. Methods: We compute structure and reaction observables for three-, four- and five-nucleon systems within the ab initio frameworks of the no-core shell model an no-core shell model with continuum. Using a small set of design runs, we construct a Gaussian process model (GPM) that acts as a statistical emulator for the theory. With this, we gain insight into how uncertainties in the 3N low-energy constants propagate throughout the calculation and determine the Bayesian posterior distribution of these parameters with Markov-Chain Monte-Carlo.
△ Less
Submitted 17 April, 2020;
originally announced April 2020.
-
Optimal Control for the Quantum Simulation of Nuclear Dynamics
Authors:
Eric T. Holland,
Kyle A. Wendt,
Konstantinos Kravvaris,
Xian Wu,
W. Erich Ormand,
Jonathan L DuBois,
Sofia Quaglioni,
Francesco Pederiva
Abstract:
We propose a method for enacting the unitary time propagation of two interacting neutrons at leading order of chiral effective field theory by efficiently encoding the nuclear dynamics into a single multi-level quantum device. The emulated output of the quantum simulation shows that, by applying a single gate that draws on the underlying characteristics of the device, it is possible to observe mul…
▽ More
We propose a method for enacting the unitary time propagation of two interacting neutrons at leading order of chiral effective field theory by efficiently encoding the nuclear dynamics into a single multi-level quantum device. The emulated output of the quantum simulation shows that, by applying a single gate that draws on the underlying characteristics of the device, it is possible to observe multiple cycles of the nucleons' dynamics before the onset of decoherence. Owing to the signal's longevity, we can then extract spectroscopic properties of the simulated nuclear system. This allows us to validate the encoding of the nuclear Hamiltonian and the robustness of the simulation in the presence of quantum-hardware noise by comparing the extracted spectroscopic information to exact calculations. This work paves the way for transformative calculations of dynamical properties of nuclei on near-term quantum devices.
△ Less
Submitted 22 August, 2019;
originally announced August 2019.
-
Discrepancy between experimental and theoretical $β$-decay rates resolved from first principles
Authors:
P. Gysbers,
G. Hagen,
J. D. Holt,
G. R. Jansen,
T. D. Morris,
P. Navratil,
T. Papenbrock,
S. Quaglioni,
A. Schwenk,
S. R. Stroberg,
K. A. Wendt
Abstract:
$β$-decay, a process that changes a neutron into a proton (and vice versa), is the dominant decay mode of atomic nuclei. This decay offers a unique window to physics beyond the standard model, and is at the heart of microphysical processes in stellar explosions and the synthesis of the elements in the Universe. For 50 years, a central puzzle has been that observed $β…
▽ More
$β$-decay, a process that changes a neutron into a proton (and vice versa), is the dominant decay mode of atomic nuclei. This decay offers a unique window to physics beyond the standard model, and is at the heart of microphysical processes in stellar explosions and the synthesis of the elements in the Universe. For 50 years, a central puzzle has been that observed $β$-decay rates are systematically smaller than theoretical predictions. This was attributed to an apparent quenching of the fundamental coupling constant $g_A \simeq $ 1.27 in the nucleus by a factor of about 0.75 compared to the $β$-decay of a free neutron. The origin of this quenching is controversial and has so far eluded a first-principles theoretical understanding. Here we address this puzzle and show that this quenching arises to a large extent from the coupling of the weak force to two nucleons as well as from strong correlations in the nucleus. We present state-of-the-art computations of $β$-decays from light to heavy nuclei. Our results are consistent with experimental data, including the pioneering measurement for $^{100}$Sn. These theoretical advances are enabled by systematic effective field theories of the strong and weak interactions combined with powerful quantum many-body techniques. This work paves the way for systematic theoretical predictions for fundamental physics problems. These include the synthesis of heavy elements in neutron star mergers and the search for neutrino-less double-$β$-decay, where an analogous quenching puzzle is a major source of uncertainty in extracting the neutrino mass scale.
△ Less
Submitted 28 February, 2019;
originally announced March 2019.
-
Recent developments in nuclear structure theory: an outlook on the muonic atom program
Authors:
Oscar Javier Hernandez,
Sonia Bacca,
Kyle Andrew Wendt
Abstract:
The discovery of the proton-radius puzzle and the subsequent deuteron-radius puzzle is fueling an on-going debate on possible explanations for the difference in the observed radii obtained from muonic atoms and from electron-nucleus systems. Atomic nuclei have a complex internal structure that must be taken into account when analyzing experimental spectroscopic results. Ab initio nuclear structure…
▽ More
The discovery of the proton-radius puzzle and the subsequent deuteron-radius puzzle is fueling an on-going debate on possible explanations for the difference in the observed radii obtained from muonic atoms and from electron-nucleus systems. Atomic nuclei have a complex internal structure that must be taken into account when analyzing experimental spectroscopic results. Ab initio nuclear structure theory provided the so far most precise estimates of important corrections to the Lamb shift in muonic atoms and is well poised to also investigate nuclear structure corrections to the hyperfine splitting in muonic atoms. Independently on whether the puzzle is due to beyond-the-standard-model physics or not, nuclear structure corrections are a necessary theoretical input to any experimental extraction of electric and magnetic radii from precise muonic atom measurements.
Here, we review the status of the calculations performed by the TRIUMF-Hebrew University group, focusing on the deuteron, and discuss preliminary results on magnetic sum rules calculated with two-body currents at next-to-leading order. Two-body currents will be an important ingredient in future calculations of nuclear structure corrections to the hyperfine splitting in muonic atoms.
△ Less
Submitted 14 December, 2017;
originally announced December 2017.
-
Unexpectedly large charge radii of neutron-rich calcium isotopes
Authors:
R. F. Garcia Ruiz,
M. L. Bissell,
K. Blaum,
A. Ekstrom,
N. Frommgen,
G. Hagen,
M. Hammen,
K. Hebeler,
J. D. Holt,
G. R. Jansen,
M. Kowalska,
K. Kreim,
W. Nazarewicz,
R. Neugart,
G. Neyens,
W. Nortershauser,
T. Papenbrock,
J. Papuga,
A. Schwenk,
J. Simonis,
K. A. Wendt,
D. T. Yordanov
Abstract:
Despite being a complex many-body system, the atomic nucleus exhibits simple structures for certain "magic" numbers of protons and neutrons. The calcium chain in particular is both unique and puzzling: evidence of doubly-magic features are known in 40,48Ca, and recently suggested in two radioactive isotopes, 52,54Ca. Although many properties of experimentally known Ca isotopes have been successful…
▽ More
Despite being a complex many-body system, the atomic nucleus exhibits simple structures for certain "magic" numbers of protons and neutrons. The calcium chain in particular is both unique and puzzling: evidence of doubly-magic features are known in 40,48Ca, and recently suggested in two radioactive isotopes, 52,54Ca. Although many properties of experimentally known Ca isotopes have been successfully described by nuclear theory, it is still a challenge to predict their charge radii evolution. Here we present the first measurements of the charge radii of 49,51,52Ca, obtained from laser spectroscopy experiments at ISOLDE, CERN. The experimental results are complemented by state-of-the-art theoretical calculations. The large and unexpected increase of the size of the neutron-rich calcium isotopes beyond N = 28 challenges the doubly-magic nature of 52Ca and opens new intriguing questions on the evolution of nuclear sizes away from stability, which are of importance for our understanding of neutron-rich atomic nuclei.
△ Less
Submitted 25 February, 2016;
originally announced February 2016.
-
Effective field theory in the harmonic oscillator basis
Authors:
S. Binder,
A. Ekström,
G. Hagen,
T. Papenbrock,
K. A. Wendt
Abstract:
We develop interactions from chiral effective field theory (EFT) that are tailored to the harmonic oscillator basis. As a consequence, ultraviolet convergence with respect to the model space is implemented by construction and infrared convergence can be achieved by enlarging the model space for the kinetic energy. In oscillator EFT, matrix elements of EFTs formulated for continuous momenta are eva…
▽ More
We develop interactions from chiral effective field theory (EFT) that are tailored to the harmonic oscillator basis. As a consequence, ultraviolet convergence with respect to the model space is implemented by construction and infrared convergence can be achieved by enlarging the model space for the kinetic energy. In oscillator EFT, matrix elements of EFTs formulated for continuous momenta are evaluated at the discrete momenta that stem from the diagonalization of the kinetic energy in the finite oscillator space. By fitting to realistic phase shifts and deuteron data we construct an effective interaction from chiral EFT at next-to-leading order. Many-body coupled-cluster calculations of nuclei up to 132Sn exhibit a fast convergence of ground-state energies and radii in feasible model spaces.
△ Less
Submitted 29 April, 2016; v1 submitted 11 December, 2015;
originally announced December 2015.
-
Charge, neutron, and weak size of the atomic nucleus
Authors:
G. Hagen,
A. Ekström,
C. Forssén,
G. R. Jansen,
W. Nazarewicz,
T. Papenbrock,
K. A. Wendt,
S. Bacca,
N. Barnea,
B. Carlsson,
C. Drischler,
K. Hebeler,
M. Hjorth-Jensen,
M. Miorelli,
G. Orlandini,
A. Schwenk,
J. Simonis
Abstract:
What is the size of the atomic nucleus? This deceivably simple question is difficult to answer. While the electric charge distributions in atomic nuclei were measured accurately already half a century ago, our knowledge of the distribution of neutrons is still deficient. In addition to constraining the size of atomic nuclei, the neutron distribution also impacts the number of nuclei that can exist…
▽ More
What is the size of the atomic nucleus? This deceivably simple question is difficult to answer. While the electric charge distributions in atomic nuclei were measured accurately already half a century ago, our knowledge of the distribution of neutrons is still deficient. In addition to constraining the size of atomic nuclei, the neutron distribution also impacts the number of nuclei that can exist and the size of neutron stars. We present an ab initio calculation of the neutron distribution of the neutron-rich nucleus $^{48}$Ca. We show that the neutron skin (difference between radii of neutron and proton distributions) is significantly smaller than previously thought. We also make predictions for the electric dipole polarizability and the weak form factor; both quantities are currently targeted by precision measurements. Based on ab initio results for $^{48}$Ca, we provide a constraint on the size of a neutron star.
△ Less
Submitted 23 September, 2015;
originally announced September 2015.
-
Uncertainty analysis and order-by-order optimization of chiral nuclear interactions
Authors:
B. D. Carlsson,
A. Ekström,
C. Forssén,
D. Fahlin Strömberg,
G. R. Jansen,
O. Lilja,
M. Lindby,
B. A. Mattsson,
K. A. Wendt
Abstract:
Chiral effective field theory (chiEFT) provides a systematic approach to describe low-energy nuclear forces. Moreover, chiEFT is able to provide well-founded estimates of statistical and systematic uncertainties -- although this unique advantage has not yet been fully exploited. We fill this gap by performing an optimization and statistical analysis of all the low-energy constants (LECs) up to nex…
▽ More
Chiral effective field theory (chiEFT) provides a systematic approach to describe low-energy nuclear forces. Moreover, chiEFT is able to provide well-founded estimates of statistical and systematic uncertainties -- although this unique advantage has not yet been fully exploited. We fill this gap by performing an optimization and statistical analysis of all the low-energy constants (LECs) up to next-to-next-to-leading order. Our optimization protocol corresponds to a simultaneous fit to scattering and bound-state observables in the pion-nucleon, nucleon-nucleon, and few-nucleon sectors, thereby utilizing the full model capabilities of chiEFT. We study the effect on other observables by demonstrating error-propagation methods that can easily be adopted by future works. We employ mathematical optimization and implement automatic differentiation to attain efficient and machine-precise first- and second-order derivatives of the objective function with respect to the LECs. We use power-counting arguments to estimate the systematic uncertainty that is inherent to chiEFT and we construct chiral interactions at different orders with quantified uncertainties. Statistical error propagation is compared with Monte Carlo sampling showing that statistical errors are in general small compared to systematic ones. In conclusion, we find that a simultaneous fit to different sets of data is critical to (i) identify the optimal set of LECs, (ii) capture all relevant correlations, (iii) reduce the statistical uncertainty, and (iv) attain order-by-order convergence in chiEFT. Furthermore, certain systematic uncertainties in the few-nucleon sector are shown to get substantially magnified in the many-body sector; in particlar when varying the cutoff in the chiral potentials. The methodology and results presented in this Paper open a new frontier for uncertainty quantification in ab initio nuclear theory.
△ Less
Submitted 27 January, 2016; v1 submitted 8 June, 2015;
originally announced June 2015.
-
Infrared length scale and extrapolations for the no-core shell model
Authors:
K. A. Wendt,
C. Forssén,
T. Papenbrock,
D. Sääf
Abstract:
We precisely determine the infrared (IR) length scale of the no-core shell model (NCSM). In the NCSM, the $A$-body Hilbert space is truncated by the total energy, and the IR length can be determined by equating the intrinsic kinetic energy of $A$ nucleons in the NCSM space to that of $A$ nucleons in a $3(A-1)$-dimensional hyper-radial well with a Dirichlet boundary condition for the hyper radius.…
▽ More
We precisely determine the infrared (IR) length scale of the no-core shell model (NCSM). In the NCSM, the $A$-body Hilbert space is truncated by the total energy, and the IR length can be determined by equating the intrinsic kinetic energy of $A$ nucleons in the NCSM space to that of $A$ nucleons in a $3(A-1)$-dimensional hyper-radial well with a Dirichlet boundary condition for the hyper radius. We demonstrate that this procedure indeed yields a very precise IR length by performing large-scale NCSM calculations for $^{6}$Li. We apply our result and perform accurate IR extrapolations for bound states of $^{4}$He, $^{6}$He, $^{6}$Li, $^{7}$Li. We also attempt to extrapolate NCSM results for $^{10}$B and $^{16}$O with bare interactions from chiral effective field theory over tens of MeV.
△ Less
Submitted 19 June, 2015; v1 submitted 24 March, 2015;
originally announced March 2015.
-
Accurate nuclear radii and binding energies from a chiral interaction
Authors:
A. Ekström,
G. R. Jansen,
K. A. Wendt,
G. Hagen,
T. Papenbrock,
B. D. Carlsson,
C. Forssén,
M. Hjorth-Jensen,
P. Navrátil,
W. Nazarewicz
Abstract:
With the goal of developing predictive ab-initio capability for light and medium-mass nuclei, two-nucleon and three-nucleon forces from chiral effective field theory are optimized simultaneously to low-energy nucleon-nucleon scattering data, as well as binding energies and radii of few-nucleon systems and selected isotopes of carbon and oxygen. Coupled-cluster calculations based on this interactio…
▽ More
With the goal of developing predictive ab-initio capability for light and medium-mass nuclei, two-nucleon and three-nucleon forces from chiral effective field theory are optimized simultaneously to low-energy nucleon-nucleon scattering data, as well as binding energies and radii of few-nucleon systems and selected isotopes of carbon and oxygen. Coupled-cluster calculations based on this interaction, named NNLO$_{\rm sat}$, yield accurate binding energies and radii of nuclei up to $^{40}$Ca, and are consistent with the empirical saturation point of symmetric nuclear matter. In addition, the low-lying collective $J^π=3^-$ states in $^{16}$O and $^{40}$Ca are described accurately, while spectra for selected $p$- and $sd$-shell nuclei are in reasonable agreement with experiment.
△ Less
Submitted 1 May, 2015; v1 submitted 16 February, 2015;
originally announced February 2015.
-
Uncertainty Quantification of the Pion-Nucleon Low-Energy Coupling Constants up to Fourth Order in Chiral Perturbation Theory
Authors:
K. A. Wendt,
B. D. Carlsson,
A. Ekström
Abstract:
We extract the statistical uncertainties for the pion-nucleon ($πN$) low energy constants (LECs) up to fourth order $\mathcal{O}(Q^4)$ in the chiral expansion of the nuclear effective Lagrangian. The LECs are optimized with respect to experimental scattering data. For comparison, we also present an uncertainty quantification that is based solely on \pin{} scattering phase shifts. Statistical error…
▽ More
We extract the statistical uncertainties for the pion-nucleon ($πN$) low energy constants (LECs) up to fourth order $\mathcal{O}(Q^4)$ in the chiral expansion of the nuclear effective Lagrangian. The LECs are optimized with respect to experimental scattering data. For comparison, we also present an uncertainty quantification that is based solely on \pin{} scattering phase shifts. Statistical errors on the LECs are critical in order to estimate the subsequent uncertainties in \textit{ab initio} modeling of light and medium mass nuclei which exploit chiral effective field theory. As an example of the this, we present the first complete predictions with uncertainty quantification of peripheral phase shifts of elastic proton-neutron scattering.
△ Less
Submitted 2 October, 2014;
originally announced October 2014.
-
Infrared extrapolations for atomic nuclei
Authors:
R. J. Furnstahl,
G. Hagen,
T. Papenbrock,
K. A. Wendt
Abstract:
Harmonic oscillator model-space truncations introduce systematic errors to the calculation of binding energies and other observables. We identify the relevant infrared scaling variable and give values for this nucleus-dependent quantity. We consider isotopes of oxygen computed with the coupled-cluster method from chiral nucleon-nucleon interactions at next-to-next-to-leading order and show that th…
▽ More
Harmonic oscillator model-space truncations introduce systematic errors to the calculation of binding energies and other observables. We identify the relevant infrared scaling variable and give values for this nucleus-dependent quantity. We consider isotopes of oxygen computed with the coupled-cluster method from chiral nucleon-nucleon interactions at next-to-next-to-leading order and show that the infrared component of the error is sufficiently understood to permit controlled extrapolations. By employing oscillator spaces with relatively large frequencies, well above the energy minimum, the ultraviolet corrections can be suppressed while infrared extrapolations over tens of MeVs are accurate for ground-state energies. However, robust uncertainty quantification for extrapolated quantities that fully accounts for systematic errors is not yet developed.
△ Less
Submitted 1 August, 2014;
originally announced August 2014.
-
Statistical uncertainties of a chiral interaction at next-to-next-to leading order
Authors:
A. Ekström,
B. D. Carlsson,
K. A. Wendt,
C. Forssén,
M. Hjorth-Jensen,
R. Machleidt,
S. M. Wild
Abstract:
We have quantified the statistical uncertainties of the low-energy coupling-constants (LECs) of an optimized nucleon-nucleon (NN) interaction from chiral effective field theory ($χ$EFT) at next-to-next-to-leading order (NNLO). In addition, we have propagated the impact of the uncertainties of the LECs to two-nucleon scattering phase shifts, effective range parameters, and deuteron observables.
We have quantified the statistical uncertainties of the low-energy coupling-constants (LECs) of an optimized nucleon-nucleon (NN) interaction from chiral effective field theory ($χ$EFT) at next-to-next-to-leading order (NNLO). In addition, we have propagated the impact of the uncertainties of the LECs to two-nucleon scattering phase shifts, effective range parameters, and deuteron observables.
△ Less
Submitted 30 June, 2014; v1 submitted 26 June, 2014;
originally announced June 2014.
-
Effects of three-nucleon forces and two-body currents on Gamow-Teller strengths
Authors:
A. Ekström,
G. R. Jansen,
K. A. Wendt,
G. Hagen,
T. Papenbrock,
S. Bacca,
B. Carlsson,
D. Gazit
Abstract:
We optimize chiral interactions at next-to-next-to leading order to observables in two- and three-nucleon systems, and compute Gamow-Teller transitions in carbon-14, oxygen-22 and oxygen-24 using consistent two-body currents. We compute spectra of the daughter nuclei nitrogen-14, fluorine-22 and fluorine-24 via an isospin-breaking coupled-cluster technique, with several predictions. The two-body c…
▽ More
We optimize chiral interactions at next-to-next-to leading order to observables in two- and three-nucleon systems, and compute Gamow-Teller transitions in carbon-14, oxygen-22 and oxygen-24 using consistent two-body currents. We compute spectra of the daughter nuclei nitrogen-14, fluorine-22 and fluorine-24 via an isospin-breaking coupled-cluster technique, with several predictions. The two-body currents reduce the Ikeda sum rule, corresponding to a quenching factor q^2 ~ 0.84-0.92 of the axial-vector coupling. The half life of carbon-14 depends on the energy of the first excited 1+ state, the three-nucleon force, and the two-body current.
△ Less
Submitted 18 June, 2014;
originally announced June 2014.
-
Coupled-cluster calculations of nucleonic matter
Authors:
G. Hagen,
T. Papenbrock,
A. Ekström,
K. A. Wendt,
G. Baardsen,
S. Gandolfi,
M. Hjorth-Jensen,
C. J. Horowitz
Abstract:
Background: The equation of state (EoS) of nucleonic matter is central for the understanding of bulk nuclear properties, the physics of neutron star crusts, and the energy release in supernova explosions.
Purpose: This work presents coupled-cluster calculations of infinite nucleonic matter using modern interactions from chiral effective field theory (EFT). It assesses the role of correlations be…
▽ More
Background: The equation of state (EoS) of nucleonic matter is central for the understanding of bulk nuclear properties, the physics of neutron star crusts, and the energy release in supernova explosions.
Purpose: This work presents coupled-cluster calculations of infinite nucleonic matter using modern interactions from chiral effective field theory (EFT). It assesses the role of correlations beyond particle-particle and hole-hole ladders, and the role of three-nucleon-forces (3NFs) in nuclear matter calculations with chiral interactions.
Methods: This work employs the optimized nucleon-nucleon NN potential NNLOopt at next-to-next-to leading-order, and presents coupled-cluster computations of the EoS for symmetric nuclear matter and neutron matter. The coupled-cluster method employs up to selected triples clusters and the single-particle space consists of a momentum-space lattice. We compare our results with benchmark calculations and control finite-size effects and shell oscillations via twist-averaged boundary conditions.
Results: We provide several benchmarks to validate the formalism and show that our results exhibit a good convergence toward the thermodynamic limit. While neutron matter with interactions from chiral EFT is perturbative, symmetric nuclear matter requires nonperturbative approaches. Correlations beyond the standard particle-particle ladder approximation yield non-negligible contributions. The saturation point of symmetric nuclear matter is sensitive to the employed 3NFs and the employed regularization scheme. 3NFs with nonlocal cutoffs exhibit a considerably improved convergence than their local cousins. We are unable to find values for the parameters of the short-range part of the local 3NF that simultaneously yield acceptable values for the saturation point in symmetric nuclear matter and the binding energies of light nuclei.
△ Less
Submitted 12 November, 2013;
originally announced November 2013.
-
Similarity Renormalization Group Evolution of Three-Nucleon Forces in a Hyperspherical Momentum Representation
Authors:
Kyle A. Wendt
Abstract:
A new framework for computing the Similarity Renormalization Group (SRG) evolution of three-nucleon forces (3NF) in momentum representation is presented. The use of antisymmetric three-particle hyperspherical momentum states ensures unitary evolutions within certain basis truncations, much like antisymmetric harmonic oscillator SRG evolutions. Additionally, in each partial wave the…
▽ More
A new framework for computing the Similarity Renormalization Group (SRG) evolution of three-nucleon forces (3NF) in momentum representation is presented. The use of antisymmetric three-particle hyperspherical momentum states ensures unitary evolutions within certain basis truncations, much like antisymmetric harmonic oscillator SRG evolutions. Additionally, in each partial wave the $T_{\mbox{rel}}$-SRG regulator is exactly represented, similar to recent 3NF momentum representation evolutions. Unitary equivalence is demonstrated for the triton using several chiral two- plus three-nucleon interactions. This method allows for a clean visualization of the evolution of the three-nucleon forces, which manifests the SRG decoupling pattern and low-momentum universality.
△ Less
Submitted 4 April, 2013;
originally announced April 2013.
-
Local Projections of Low-Momentum Potentials
Authors:
K. A. Wendt,
R. J. Furnstahl,
S. Ramanan
Abstract:
Nuclear interactions evolved via renormalization group methods to lower resolution become increasingly non-local (off-diagonal in coordinate space) as they are softened. This inhibits both the development of intuition about the interactions and their use with some methods for solving the quantum many-body problem. By applying "local projections", a softened interaction can be reduced to a local ef…
▽ More
Nuclear interactions evolved via renormalization group methods to lower resolution become increasingly non-local (off-diagonal in coordinate space) as they are softened. This inhibits both the development of intuition about the interactions and their use with some methods for solving the quantum many-body problem. By applying "local projections", a softened interaction can be reduced to a local effective interaction plus a non-local residual interaction. At the two-body level, a local projection after similarity renormalization group (SRG) evolution manifests the elimination of short-range repulsive cores and the flow toward universal low-momentum interactions. The SRG residual interaction is found to be relatively weak at low energy, which motivates a perturbative treatment.
△ Less
Submitted 4 April, 2012; v1 submitted 27 March, 2012;
originally announced March 2012.
-
Decoupling of Spurious Deep Bound States with the Similarity Renormalization Group
Authors:
K. A. Wendt,
R. J. Furnstahl,
R. J. Perry
Abstract:
The Similarity Renormalization Group (SRG) is a continuous series of unitary transformations that can be implemented as a flow equation. When the relative kinetic energy ($\Trel$) is used in the SRG generator, nuclear structure calculations have shown greatly improved convergence with basis size because of the decoupling of high-energy and low-energy physics. However this generator can sometimes b…
▽ More
The Similarity Renormalization Group (SRG) is a continuous series of unitary transformations that can be implemented as a flow equation. When the relative kinetic energy ($\Trel$) is used in the SRG generator, nuclear structure calculations have shown greatly improved convergence with basis size because of the decoupling of high-energy and low-energy physics. However this generator can sometimes be problematic. A test case is provided by a study of initial interactions from chiral effective field theories with large cutoffs, which can lead to spurious deep bound states. We would like the SRG to decouple these from the physical shallow bound states. However, with $\Trel$ the high- and low-energy bound states are not decoupled in the usual sense. Replacing $\Trel$ by the momentum-space diagonal of the Hamiltonian ($H_d$) in the SRG generator does produce decoupling, such that the shallow states are in the low-momentum region and the deep bound states are at higher momentum. The flow toward universal low-momentum interactions is also restored.
△ Less
Submitted 13 January, 2011;
originally announced January 2011.
-
Toward large-scale Hybrid Monte Carlo simulations of the Hubbard model on graphics processing units
Authors:
Kyle A. Wendt,
Joaquín E. Drut,
Timo A. Lähde
Abstract:
The performance of the Hybrid Monte Carlo algorithm is determined by the speed of sparse matrix-vector multiplication within the context of preconditioned conjugate gradient iteration. We study these operations as implemented for the fermion matrix of the Hubbard model in d+1 space-time dimensions, and report a performance comparison between a 2.66 GHz Intel Xeon E5430 CPU and an NVIDIA Tesla C106…
▽ More
The performance of the Hybrid Monte Carlo algorithm is determined by the speed of sparse matrix-vector multiplication within the context of preconditioned conjugate gradient iteration. We study these operations as implemented for the fermion matrix of the Hubbard model in d+1 space-time dimensions, and report a performance comparison between a 2.66 GHz Intel Xeon E5430 CPU and an NVIDIA Tesla C1060 GPU using double-precision arithmetic. We find speedup factors ranging between 30-350 for d = 1, and in excess of 40 for d = 3. We argue that such speedups are of considerable impact for large-scale simulational studies of quantum many-body systems.
△ Less
Submitted 20 July, 2010;
originally announced July 2010.
