Showing 1–50 of 77 results for author: Quaglioni, S

Ab initio calculation of the $^3$He$(αあるふぁ,γがんま)^7$Be astrophysical S factor with chiral two- and three-nucleon forces

Authors: M. C. Atkinson, K. Kravvaris, S. Quaglioni, P. Navrátil

Abstract: The $^3$He$(αあるふぁ,γがんま)^7$Be radiative capture reaction plays a key role in the creation of elements in stars as well as in the production of solar neutrinos, the observation of which is one of the main tools to study the properties of our sun. Since accurate experimental measurements of this fusion cross section at solar energies are difficult due to the strong Coulomb repulsion between the reactants, t… ▽ More The $^3$He$(αあるふぁ,γがんま)^7$Be radiative capture reaction plays a key role in the creation of elements in stars as well as in the production of solar neutrinos, the observation of which is one of the main tools to study the properties of our sun. Since accurate experimental measurements of this fusion cross section at solar energies are difficult due to the strong Coulomb repulsion between the reactants, the onus falls on theory to provide a robust means for extrapolating from the region where experimental data is available down to the desired astrophysical regime. We present the first microscopic calculations of $^3$He$(αあるふぁ,γがんま)^7$Be with explicit inclusion of three-nucleon forces. Our prediction of the astrophysical $S$ factor qualitatively agrees with experimental data. We further incorporate experimental bound-state and scattering information in our calculation to arrive at a more quantitative description. This process reveals that our current model lacks sufficient repulsion in the $1/2^+$ channel of our model space to simultaneously reproduce elastic-scattering data. This deficit suggests that $^3$He$(αあるふぁ,γがんま)^7$Be probes aspects of the nuclear force that are not currently well-constrained. △ Less

Submitted 14 September, 2024; originally announced September 2024.

Comments: 7 pages, 4 figures, 2 tables

Evaluation of phase shifts for non-relativistic elastic scattering using quantum computers

Authors: Francesco Turro, Kyle A. Wendt, Sofia Quaglioni, Francesco Pederiva, Alessandro Roggero

Abstract: Simulations of scattering processes are essential in understanding the physics of our universe. Computing relevant scattering quantities from ab initio methods is extremely difficult on classical devices because of the substantial computational resources needed. This work reports the development of an algorithm that makes it possible to obtain phase shifts for generic non-relativistic elastic scat… ▽ More Simulations of scattering processes are essential in understanding the physics of our universe. Computing relevant scattering quantities from ab initio methods is extremely difficult on classical devices because of the substantial computational resources needed. This work reports the development of an algorithm that makes it possible to obtain phase shifts for generic non-relativistic elastic scattering processes on a quantum computer. This algorithm is based on extracting phase shifts from the direct implementation of the real-time evolution. The algorithm is improved by a variational procedure, making it more accurate and resistant to the quantum noise. The reliability of the algorithm is first demonstrated by means of classical numerical simulations for different potentials, and later tested on existing quantum hardware, specifically on IBM quantum processors. △ Less

Submitted 4 July, 2024; originally announced July 2024.

Comments: 17 pages, 15 figures, 4 tables

Report number: LLNL-JRNL-865762-DRAFT; IQuS@UW-21-080

Tunable-fidelity wave functions for the \textit{ab initio} description of scattering and reactions

Authors: Konstantinos Kravvaris, Sofia Quaglioni, Petr Navratil

Abstract: The no-core shell model (NCSM) is an \textit{ab initio} method that solves the nuclear many-body problem by expanding the many-particle wave function into a (typically) harmonic oscillator basis and minimizing the energy to obtain the expansion coefficients. Extensions of the NCSM, such as its coupling with microscopic-cluster basis states, further allow for an \textit{ab initio} treatment of ligh… ▽ More The no-core shell model (NCSM) is an \textit{ab initio} method that solves the nuclear many-body problem by expanding the many-particle wave function into a (typically) harmonic oscillator basis and minimizing the energy to obtain the expansion coefficients. Extensions of the NCSM, such as its coupling with microscopic-cluster basis states, further allow for an \textit{ab initio} treatment of light-ion nuclear reactions of interest for both astrophysics and nuclear technology applications. A downside of the method is the exponential scaling of the basis size with increasing number of nucleons and excitation quanta, which limits its applicability to mass $A\lesssim 16$ nuclei, except for variants where the basis is further down-selected via some truncation scheme. We consider a basis selection method for the NCSM that captures the essential degrees of freedom of the nuclear wave function leading to a favorable complexity scaling for calculations and enabling \textit{ab initio} reaction calculations in $sd$-shell nuclei. The particle configurations within the NCSM basis are ordered based on their contribution to the first moment of the Hamiltonian matrix that results from the projection onto the many-body basis. The truncation scheme then consists in retaining only the lowest-first-moment configurations, which typically contain only few many-body basis states (Slater determinants). We present calculations for $^7$Li and $n+^{12}$C scattering using nucleon-nucleon interactions derived from chiral effective field theory and softened using the similarity renormalization group method. The obtained energy levels invariably demonstrate exponential convergence with the size of the basis, and we find improved convergence in scattering calculations. To demonstrate the possibilities enabled by the approach, we also present a first calculation for the scattering of neutrons from $^{24}$Mg. △ Less

Submitted 4 February, 2024; originally announced February 2024.

Comments: 13 pages, 8 figures

Ab initio investigation of the $^7$Li($p,e^+e^-$)$^8$Be process and the X17 boson

Authors: P. Gysbers, P. Navratil, K. Kravvaris, G. Hupin, S. Quaglioni

Abstract: Observations of anomalies in the electron-positron angular correlations in high-energy decays in $^4$He, $^8$Be, and $^{12}$C have been reported recently by the ATOMKI collaboration. These could be explained by the creation and subsequent decay of a new boson with a mass of ${\sim}17$ MeV. Theoretical understanding of pair creation in the proton capture reactions used in these experiments is impor… ▽ More Observations of anomalies in the electron-positron angular correlations in high-energy decays in $^4$He, $^8$Be, and $^{12}$C have been reported recently by the ATOMKI collaboration. These could be explained by the creation and subsequent decay of a new boson with a mass of ${\sim}17$ MeV. Theoretical understanding of pair creation in the proton capture reactions used in these experiments is important for the interpretation of the anomalies. We apply the ab initio No-Core Shell Model with Continuum (NCSMC) to the proton capture on $^7$Li. The NCSMC describes both bound and unbound states in light nuclei in a unified way with chiral two- and three-nucleon interactions as the only input. We investigate the structure of $^8$Be, the $p+^7$Li elastic scattering, the $^7$Li($p,γがんま$)$^8$Be cross section and the internal pair creation $^7$Li($p,e^+ e^-$)$^8$Be. We discuss the impact of a proper treatment of the initial scattering state on the electron-positron angular correlation spectrum and compare our results to available ATOMKI data sets. Finally, we calculate $^7$Li($p,X$)$^8$Be cross sections for several proposed models of the hypothetical X17 particle. △ Less

Submitted 12 July, 2024; v1 submitted 25 August, 2023; originally announced August 2023.

Comments: 32 pages, 13 figures. Matches more closely the published version

Journal ref: Physical Review C 110, 015503 (2024)

Impact of the $^6$Li asymptotic normalization constant onto $αあるふぁ$-induced reactions of astrophysical interest

Authors: Chloë Hebborn, Melina L. Avila, Konstantinos Kravvaris, Gregory Potel, Sofia Quaglioni

Abstract: Indirect methods have become the predominant approach in experimental nuclear astrophysics for studying several low-energy nuclear reactions occurring in stars, as direct measurements of many of these relevant reactions are rendered infeasible due to their low reaction probability. Such indirect methods, however, require theoretical input that in turn can have significant poorly-quantified uncerta… ▽ More Indirect methods have become the predominant approach in experimental nuclear astrophysics for studying several low-energy nuclear reactions occurring in stars, as direct measurements of many of these relevant reactions are rendered infeasible due to their low reaction probability. Such indirect methods, however, require theoretical input that in turn can have significant poorly-quantified uncertainties, which can then be propagated to the reaction rates and have a large effect on our quantitative understanding of stellar evolution and nucleosynthesis processes. We present two such examples involving $αあるふぁ$-induced reactions, $^{13}$C($αあるふぁ,n)^{16}$O and $^{12}$C$(αあるふぁ,γがんま)^{16}$O, for which the low-energy cross sections have been constrained with $(^6$Li$,d)$ transfer data. In this Letter, we discuss how a first-principle calculation of $^6$Li leads to a 21% reduction of the $^{12}$C$(αあるふぁ,γがんま)^{16}$O cross sections with respect to a previous estimation. This calculation further resolves the discrepancy between recent measurements of the $^{13}$C$(αあるふぁ,n)^{16}$O reaction and points to the need for improved theoretical formulations of nuclear reactions. △ Less

Submitted 3 June, 2024; v1 submitted 11 July, 2023; originally announced July 2023.

Comments: 7 pages (including references) and 3 figures

Report number: LLNL-JRNL-851130

Journal ref: Phys. Rev. C 109, L061601 (2024)

Quantum Information Science and Technology for Nuclear Physics. Input into U.S. Long-Range Planning, 2023

Authors: Douglas Beck, Joseph Carlson, Zohreh Davoudi, Joseph Formaggio, Sofia Quaglioni, Martin Savage, Joao Barata, Tanmoy Bhattacharya, Michael Bishof, Ian Cloet, Andrea Delgado, Michael DeMarco, Caleb Fink, Adrien Florio, Marianne Francois, Dorota Grabowska, Shannon Hoogerheide, Mengyao Huang, Kazuki Ikeda, Marc Illa, Kyungseon Joo, Dmitri Kharzeev, Karol Kowalski, Wai Kin Lai, Kyle Leach , et al. (76 additional authors not shown)

Abstract: In preparation for the 2023 NSAC Long Range Plan (LRP), members of the Nuclear Science community gathered to discuss the current state of, and plans for further leveraging opportunities in, QIST in NP research at the Quantum Information Science for U.S. Nuclear Physics Long Range Planning workshop, held in Santa Fe, New Mexico on January 31 - February 1, 2023. The workshop included 45 in-person pa… ▽ More In preparation for the 2023 NSAC Long Range Plan (LRP), members of the Nuclear Science community gathered to discuss the current state of, and plans for further leveraging opportunities in, QIST in NP research at the Quantum Information Science for U.S. Nuclear Physics Long Range Planning workshop, held in Santa Fe, New Mexico on January 31 - February 1, 2023. The workshop included 45 in-person participants and 53 remote attendees. The outcome of the workshop identified strategic plans and requirements for the next 5-10 years to advance quantum sensing and quantum simulations within NP, and to develop a diverse quantum-ready workforce. The plans include resolutions endorsed by the participants to address the compelling scientific opportunities at the intersections of NP and QIST. These endorsements are aligned with similar affirmations by the LRP Computational Nuclear Physics and AI/ML Workshop, the Nuclear Structure, Reactions, and Astrophysics LRP Town Hall, and the Fundamental Symmetries, Neutrons, and Neutrinos LRP Town Hall communities. △ Less

Submitted 28 February, 2023; originally announced March 2023.

Comments: A white paper for the 2023 nuclear physics long-range planning activity, emerging from the workshop "Quantum Information Science for U.S. Nuclear Physics Long Range Planning'', held in Santa Fe, New Mexico on January 31 - February 1, 2023. 26 pages with 7 figures

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.

Comments: 12 pages, 10 figures

Report number: IQuS@UW-21-043; LLNL-JRNL-853700

Ab initio investigations of A=8 nuclei: $αあるふぁ{-}αあるふぁ$ scattering, deformation in $^8$He, radiative capture of protons on $^7$Be and $^7$Li and the X17 boson

Authors: P. Navratil, K. Kravvaris, P. Gysbers, C. Hebborn, G. Hupin, S. Quaglioni

Abstract: We apply the No-Core Shell Model with Continuum (NCSMC) that is capable of describing both bound and unbound states in light nuclei in a unified way with chiral two- and three-nucleon interactions as the only input. The NCSMC can predict structure and dynamics of light nuclei and, by comparing to available experimental data, test the quality of chiral nuclear forces. We discuss applications of NCS… ▽ More We apply the No-Core Shell Model with Continuum (NCSMC) that is capable of describing both bound and unbound states in light nuclei in a unified way with chiral two- and three-nucleon interactions as the only input. The NCSMC can predict structure and dynamics of light nuclei and, by comparing to available experimental data, test the quality of chiral nuclear forces. We discuss applications of NCSMC to the $αあるふぁ{-}αあるふぁ$ scattering and the structure of $^8$Be, the p+$^7$Be and p+$^7$Li radiative capture and the production of the hypothetical X17 boson claimed in ATOMKI experiments. The $^7$Be(p,$γがんま$)$^8$B reaction plays a role in Solar nucleosynthesis and Solar neutrino physics and has been subject of numerous experimental investigations. We also highlight our investigation of the neutron rich exotic $^8$He that has been recently studied experimentally at TRIUMF with an unexpected deformation reported. △ Less

Submitted 30 November, 2022; originally announced December 2022.

Comments: 6 pages, 4 figures, contribution to the proceedings of the International Nuclear Physics Conference 2022 (INPC 2022)

Quantifying uncertainties due to irreducible three-body forces in deuteron-nucleus reactions

Authors: Linda Hlophe, Konstantinos Kravvaris, Sofia Quaglioni

Abstract: \noindent{\bf Background:} Deuteron-induced nuclear reactions are an essential tool for probing the structure of nuclei as well as astrophysical information such as $(n,γがんま)$ cross sections. The deuteron-nucleus system is typically described within a Faddeev three-body model consisting of a neutron ($n$), a proton ($p$), and the target nucleus ($A$) interacting through pairwise phenomenological pote… ▽ More \noindent{\bf Background:} Deuteron-induced nuclear reactions are an essential tool for probing the structure of nuclei as well as astrophysical information such as $(n,γがんま)$ cross sections. The deuteron-nucleus system is typically described within a Faddeev three-body model consisting of a neutron ($n$), a proton ($p$), and the target nucleus ($A$) interacting through pairwise phenomenological potentials. While Faddeev techniques enable the exact description of the three-body dynamics, their predictive power is limited in part by the omission of irreducible neutron-proton-nucleus three-body force ($n$-$p$-$A$ 3BF). {\bf Results:} By comparing the Faddeev and NCSM/RGM results, we show that the irreducible $n$-$p$-$αあるふぁ$ 3BF has a non-negligible effect on bound state and scattering observables alike. Specifically, the Faddeev approach %are yields a $^6$Li ground state that is approximately $600$~keV shallower than the one obtained with the NCSM/RGM. Additionally, the Faddeev calculations for $d$+$αあるふぁ$ scattering yield a $3^+$ resonance that is located approximately $400$~keV higher in energy compared to the NCSM/RGM result. The shape of the $d$+$αあるふぁ$ angular distributions computed using the two approaches also differ, owing to the discrepancy in the predictions of the $3^+$ resonance energy. △ Less

Submitted 22 August, 2022; originally announced August 2022.

Ab Initio Nuclear Reaction Theory with Applications to Astrophysics

Authors: Petr Navratil, Sofia Quaglioni

Abstract: We present an introduction to ab initio nuclear theory with a focus on nuclear reactions. After a high-level overview of ab initio approaches in nuclear physics, we give a more detailed description of the no-core shell model technique equivalent to a large extent to configuration-interaction methods applied in quantum chemistry. We then introduce the no-core shell model with continuum approach tha… ▽ More We present an introduction to ab initio nuclear theory with a focus on nuclear reactions. After a high-level overview of ab initio approaches in nuclear physics, we give a more detailed description of the no-core shell model technique equivalent to a large extent to configuration-interaction methods applied in quantum chemistry. We then introduce the no-core shell model with continuum approach that provides a quantum many-body description of nuclear reactions. After a brief review of nuclear reactions important for astrophysics, we present examples of results of ab initio calculations of radiative capture and transfer reactions in light nuclei. △ Less

Submitted 3 April, 2022; originally announced April 2022.

Comments: Contribution to the "Handbook of Nuclear Physics", Springer, 2022, edited by I. Tanihata, H. Toki, and T. Kajino

Ab initio prediction of the $^4{\rm He}(d,γがんま)\,^6\rm Li$ big bang radiative capture

Authors: Chloë Hebborn, Guillaume Hupin, Konstantinos Kravvaris, Sofia Quaglioni, Petr Navrátil, Peter Gysbers

Abstract: The rate at which helium ($^4$He) and deuterium ($d$) fuse together to produce lithium-6 ($^6$Li) and a $γがんま$ ray, $^4$He$(d,γがんま)^6$Li, is a critical puzzle piece in resolving the roughly three orders of magnitude discrepancy between big bang predictions and astronomical observations for the primordial abundance of $^6$Li. The accurate determination of this radiative capture rate requires the quantita… ▽ More The rate at which helium ($^4$He) and deuterium ($d$) fuse together to produce lithium-6 ($^6$Li) and a $γがんま$ ray, $^4$He$(d,γがんま)^6$Li, is a critical puzzle piece in resolving the roughly three orders of magnitude discrepancy between big bang predictions and astronomical observations for the primordial abundance of $^6$Li. The accurate determination of this radiative capture rate requires the quantitative and predictive description of the fusion probability across the big bang energy window ($30$ keV $\lesssim E\lesssim 400$ keV), where measurements are hindered by low counting rates. We present first-principles (or, ab initio) predictions of the $^4$He$(d,γがんま)^6$Li astrophysical S-factor using validated nucleon-nucleon and three-nucleon interactions derived within the framework of chiral effective field theory. By employing the ab initio no-core shell model with continuum to describe $^4{\rm He}$-$d$ scattering dynamics and bound $^6\rm Li$ product on an equal footing, we accurately and consistently determine the contributions of the main electromagnetic transitions driving the radiative capture process. Our results reveal an enhancement of the capture probability below 100 keV owing to previously neglected magnetic dipole (M1) transitions and reduce by an average factor of 7 the uncertainty of the thermonuclear capture rate between $0.002$ and $2$ GK. △ Less

Submitted 4 July, 2022; v1 submitted 29 March, 2022; originally announced March 2022.

Comments: 5 pages, 2 figures and 6 pages supplemental material

Report number: LLNL-JRNL-832850

Ab initio calculation of the $βべーた$-decay from $^{11}$Be to a p${+}^{10}$Be resonance

Authors: M. C. Atkinson, P. Navrátil, G. Hupin, K. Kravvaris, S. Quaglioni

Abstract: The exotic $βべーた$-delayed proton emission is calculated in $^{11}$Be from first principles using chiral two- and three-nucleon forces. To investigate the unexpectedly-large branching ratio measured in [PRL 123, 082501 (2019)] we calculate the proposed $(1/2^+,1/2)$ proton resonance in $^{11}$B using the no-core shell model with continuum. This calculation helps to address whether this enhancement is… ▽ More The exotic $βべーた$-delayed proton emission is calculated in $^{11}$Be from first principles using chiral two- and three-nucleon forces. To investigate the unexpectedly-large branching ratio measured in [PRL 123, 082501 (2019)] we calculate the proposed $(1/2^+,1/2)$ proton resonance in $^{11}$B using the no-core shell model with continuum. This calculation helps to address whether this enhancement is caused by unknown dark decay modes or an unobserved proton resonance. We report a branching ratio of $b_p = (1.3\pm0.5)\times10^{-6}$, suggesting that its unexpectedly-large value is caused by an unobserved proton resonance in $^{11}$B. △ Less

Submitted 26 March, 2022; originally announced March 2022.

Comments: 12 pages, 5 figures

Journal ref: Physical Review C 105, 054316 (2022)

Ab initio informed evaluation of the radiative capture of protons on $^7$Be

Authors: Konstantinos Kravvaris, Petr Navrátil, Sofia Quaglioni, Chloë Hebborn, Guillaume Hupin

Abstract: The radiative capture of protons by $^7$Be, which is the source of $^8$B that $βべーた$-decays emitting the majority of solar neutrinos measured on earth, has not yet been measured at astrophysically relevant energies. The recommended value for its zero-energy S-factor, $S_{17}$(0) = 20.8$\pm$(0.7)exp$\pm$(1.4)theory eV$\cdot$b, relies on theoretical extrapolations from higher-energy measurements, a pro… ▽ More The radiative capture of protons by $^7$Be, which is the source of $^8$B that $βべーた$-decays emitting the majority of solar neutrinos measured on earth, has not yet been measured at astrophysically relevant energies. The recommended value for its zero-energy S-factor, $S_{17}$(0) = 20.8$\pm$(0.7)exp$\pm$(1.4)theory eV$\cdot$b, relies on theoretical extrapolations from higher-energy measurements, a process that leads to significant uncertainty. We performed a set of first-principle (or, ab initio) calculations of the $^7$Be($p$, $γがんま$)$^8$B reaction to provide an independent prediction of the low-energy S-factor with quantified uncertainties. We demonstrate underlying features in the predicted S-factor allowing the combination of theoretical calculations and measurements to produce an evaluated S-factor of $S_{17}$(0) = 19.8$\pm$0.3 eV$\cdot$b. We expect the calculations and uncertainty quantification process described here to set a new standard for the evaluation of light-ion astrophysical reactions. △ Less

Submitted 23 February, 2022; originally announced February 2022.

Comments: 7 pages, 3 figures

Report number: LLNL-JRNL-832005

Journal ref: Physics Letters B Volume 845 (2023)

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.

Comments: 12 pages, 10 figures

Report number: LLNL-JRNL-829224-DRAFT

Efficacy of the symmetry-adapted basis for ab initio nucleon-nucleus interactions for light- and intermediate-mass nuclei

Authors: Alexis Mercenne, Kristina D. Launey, Tomas Dytrych, Jutta E. Escher, Sofia Quaglioni, Grigor H. Sargsyan, Jerry P. Draayer

Abstract: We study the efficacy of a new ab initio framework that combines the symmetry-adapted (SA) no-core shell-model approach with the resonating group method (RGM) for unified descriptions of nuclear structure and reactions. We obtain ab initio neutron-nucleus interactions for $^4$He, $^{16}$O, and $^{20}$Ne targets, starting with realistic nucleon-nucleon potentials. We discuss the effect of increasin… ▽ More We study the efficacy of a new ab initio framework that combines the symmetry-adapted (SA) no-core shell-model approach with the resonating group method (RGM) for unified descriptions of nuclear structure and reactions. We obtain ab initio neutron-nucleus interactions for $^4$He, $^{16}$O, and $^{20}$Ne targets, starting with realistic nucleon-nucleon potentials. We discuss the effect of increasing model space sizes and symmetry-based selections on the SA-RGM norm and direct potential kernels, as well as on phase shifts, which are the input to calculations of cross sections. We demonstrate the efficacy of the SA basis and its scalability with particle numbers and model space dimensions, with a view toward ab initio descriptions of nucleon scattering and capture reactions up through the medium-mass region. △ Less

Submitted 2 November, 2021; originally announced November 2021.

Comments: 10 pages, 6 figures

Proton inelastic scattering reveals deformation in $^8$He

Authors: M. Holl, R. Kanungo, Z. H. Sun, G. Hagen, J. A. Lay, A. M. Moro, P. Navrátil, T. Papenbrock, M. Alcorta, D. Connolly, B. Davids, A. Diaz Varela, M. Gennari, G. Hackman, J. Henderson, S. Ishimoto, A. I. Kilic, R. Krücken, A. Lennarz, J. Liang, J. Measures, W. Mittig, O. Paetkau, A. Psaltis, S. Quaglioni , et al. (5 additional authors not shown)

Abstract: A measurement of proton inelastic scattering of $^8$He at $8.25A$~MeV at TRIUMF shows a resonance at 3.54(6)~MeV with a width of 0.89(11)~MeV. The energy of the state is in good agreement with coupled cluster and no-core shell model with continuum calculations, with the latter successfully describing the measured resonance width as well. Its differential cross section analyzed with phenomenologica… ▽ More A measurement of proton inelastic scattering of $^8$He at $8.25A$~MeV at TRIUMF shows a resonance at 3.54(6)~MeV with a width of 0.89(11)~MeV. The energy of the state is in good agreement with coupled cluster and no-core shell model with continuum calculations, with the latter successfully describing the measured resonance width as well. Its differential cross section analyzed with phenomenological collective excitation form factor and microscopic coupled reaction channels framework consistently reveals a large deformation parameter $βべーた_2$ = 0.40(3), consistent with no-core shell model predictions of a large neutron deformation. This deformed double-closed shell at the neutron drip-line opens a new paradigm. △ Less

Submitted 4 October, 2021; originally announced October 2021.

Comments: 9 pages, 5 figures. Accepted for publication in Physics Letters B

Journal ref: Phys. Lett. B 822, 136710 (2021)

Using spin alignment of inelastically-excited fast beams to make spin assignments: the spectroscopy of 13O as a test case

Authors: R. J. Charity, T. B. Webb, J. M. Elson, D. E. M. Hoff, C. D. Pruitt, L. G. Sobotka, P. Navratil, G. Hupin, K. Kravvaris, S. Quaglioni, K. W. Brown, G. Cerizza, J. Estee, W. G. Lynch, J. Manfredi, P. Morfouace, C. Santamaria, S. Sweany, M. B. Tsang, T. Tsang, K. Zhu, S. A. Kuvin, D. McNeel, J. Smith, A. H. Wousmaa , et al. (1 additional authors not shown)

Abstract: Excited states in 13O were investigated using inelastic scattering of an E/A=69.5-MeV 13O beam off of a 9Be target. The excited states were identified in the invariant-mass spectra of the decay products. Both single proton and sequential two-proton decays of the excited states were examined. For a number of the excited states, the protons were emitted with strong anisotropy where emissions transve… ▽ More Excited states in 13O were investigated using inelastic scattering of an E/A=69.5-MeV 13O beam off of a 9Be target. The excited states were identified in the invariant-mass spectra of the decay products. Both single proton and sequential two-proton decays of the excited states were examined. For a number of the excited states, the protons were emitted with strong anisotropy where emissions transverse to the beam axis are favored. The measured proton-decay angular distributions were compared to predictions from distorted-wave born-approximation (DWBA) calculations of the spin alignment which was shown to be largely independent of the excitation mechanism. The deduced $^{13}$O level scheme is compared to ab initio no-core shell model with continuum (NCSMC) predictions. The lowest-energy excited states decay isotropically consistent with predictions of strong proton 1s1/2 structure. Above these states in the level scheme, we observed a number of higher-spin states not predicted within the model. Possibly these are associated with rotational bands built on deformed cluster configurations predicted by antisymmetrized molecular dynamics (AMD) calculations. The spin alignment mechanism is shown to be useful for making spin assignments and may have widespread use. △ Less

Submitted 7 July, 2021; originally announced July 2021.

Journal ref: Phys. Rev. C 104, 024325 (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.

Comments: 12 pages, 8 figures, 3 tables

Report number: LLNL-JRNL-827423-DRAFT

Nuclear Resonances, Scattering and Reactions from First Principles: Progress and Prospects

Authors: Sofia Quaglioni, Petr Navratil

Abstract: We present a brief overview of recent developments in ab initio calculations of nuclear scattering and reactions with a focus on applications of the no-core shell model with continuum method. We present a brief overview of recent developments in ab initio calculations of nuclear scattering and reactions with a focus on applications of the no-core shell model with continuum method. △ Less

Submitted 19 January, 2021; originally announced January 2021.

Comments: 9 pages, 7 figures (one more figure than in the published version)

Journal ref: Nuclear Physics News, 30:4, 12-16 (2020)

Ab initio framework for nuclear scattering and reactions induced by light projectiles

Authors: Konstantinos Kravvaris, Sofia Quaglioni, Guillaume Hupin, Petr Navratil

Abstract: A quantitative and predictive microscopic theoretical framework that can describe reactions induced by $αあるふぁ$ particles ($^4$He nuclei) and heavier projectiles is currently lacking. Such a framework would contribute to reducing uncertainty in the modeling of stellar evolution and nucleosynthesis and provide the basis for achieving a comprehensive understanding of the phenomenon of nuclear clustering… ▽ More A quantitative and predictive microscopic theoretical framework that can describe reactions induced by $αあるふぁ$ particles ($^4$He nuclei) and heavier projectiles is currently lacking. Such a framework would contribute to reducing uncertainty in the modeling of stellar evolution and nucleosynthesis and provide the basis for achieving a comprehensive understanding of the phenomenon of nuclear clustering (the organization of protons and neutrons into distinct substructures within a nucleus). We have developed an efficient and general configuration-interaction framework for the description of low-energy reactions and clustering in light nuclei. The new formalism takes full advantage of powerful second-quantization techniques, enabling the description of $αあるふぁ$-$αあるふぁ$ scattering and an exploration of clustering in the exotic $^{12}$Be nucleus. We find that the $^4$He($αあるふぁ$, $αあるふぁ$)$^4$He differential cross section computed with non-locally regulated chiral interactions is in good agreement with experimental data. Our results for $^{12}$Be indicate the presence of strongly mixed helium-cluster states consistent with a molecular-like picture surviving far above the $^6$He+$^6$He threshold, and reveal the strong influence of neutron decay in both the $^{12}$Be spectrum and in the $^6$He($^6$He,$αあるふぁ$)$^8$He cross section. We expect that this approach will enable the description of helium burning cross sections and provide insight on how three-nucleon forces influence the emergence of clustering in nuclei. △ Less

Submitted 30 November, 2020; originally announced December 2020.

Report number: LLNL-JRNL-816124

Microscopic investigation of the $^8$Li($n, γがんま$)$^9$Li reaction

Authors: Callum McCracken, Petr Navratil, Anna McCoy, Sofia Quaglioni, Guillaume Hupin

Abstract: The $^8$Li($n,γがんま$)$^9$Li reaction plays an important role in several astrophysics scenarios. It cannot be measured directly and indirect experiments have so far provided only cross section limits. Theoretical predictions differ by an order of magnitude. In this work we study the properties of $^9$Li bound states and low-lying resonances and calculate the $^8$Li($n,γがんま$)$^9$Li cross section within the… ▽ More The $^8$Li($n,γがんま$)$^9$Li reaction plays an important role in several astrophysics scenarios. It cannot be measured directly and indirect experiments have so far provided only cross section limits. Theoretical predictions differ by an order of magnitude. In this work we study the properties of $^9$Li bound states and low-lying resonances and calculate the $^8$Li($n,γがんま$)$^9$Li cross section within the no-core shell model with continuum (NCSMC) with chiral nucleon-nucleon and three-nucleon interactions as the only input. The NCSMC is an ab initio method applicable to light nuclei that provides a unified description of bound and scattering states well suited to calculate low-energy nuclear scattering and reactions. Our calculations reproduce the experimentally known bound states as well as the lowest $5/2^-$ resonance of $^9$Li. We predict a $3/2^-$ spin-parity assignment for the resonance observed at 5.38 MeV. In addition to the a very narrow $7/2^-$ resonance corresponding presumably to the experimental 6.43 MeV state, we find several other broad low-lying resonances. Our calculated $^8$Li($n,γがんま$)$^9$Li cross section is within the limits derived from the 1998 National Superconducting Cyclotron Laboratory Coulomb-dissociation experiment [Phys. Rev. C {\bf 57}, 959 (1998)]. However, it is higher than cross sections obtained in recent phenomenological studies. It is dominated by a direct E1 capture to the ground state with a resonant contribution at $\sim0.2$ MeV due to E2/M1 radiation enhanced by the $5/2^-$ resonance. △ Less

Submitted 10 February, 2021; v1 submitted 1 September, 2020; originally announced September 2020.

Comments: 12 pages, 9 figures, version accepted for publication

Report number: LLNL-JRNL-814114-DRAFT

Journal ref: Phys. Rev. C 103, 035801 (2021)

Coupled-cluster calculations of neutrinoless double-beta decay in $^{48}$Ca

Authors: S. J. Novario, P. Gysbers, J. Engel, G. Hagen, G. R. Jansen, T. D. Morris, P. Navrátil, T. Papenbrock, S. Quaglioni

Abstract: We use coupled-cluster theory and nuclear interactions from chiral effective field theory to compute the nuclear matrix element for the neutrinoless double-beta decay of $^{48}$Ca. Benchmarks with the no-core shell model in several light nuclei inform us about the accuracy of our approach. For $^{48}$Ca we find a relatively small matrix element. We also compute the nuclear matrix element for the t… ▽ More We use coupled-cluster theory and nuclear interactions from chiral effective field theory to compute the nuclear matrix element for the neutrinoless double-beta decay of $^{48}$Ca. Benchmarks with the no-core shell model in several light nuclei inform us about the accuracy of our approach. For $^{48}$Ca we find a relatively small matrix element. We also compute the nuclear matrix element for the two-neutrino double-beta decay of $^{48}$Ca with a quenching factor deduced from two-body currents in recent ab-initio calculation of the Ikeda sum-rule in $^{48}$Ca [Gysbers et al., Nature Physics 15, 428-431 (2019)]. △ Less

Submitted 12 May, 2021; v1 submitted 21 August, 2020; originally announced August 2020.

Comments: 14 pages, 13 figures; Version accepted for publication, Supplemental material also updated

Journal ref: Phys. Rev. Lett. 126, 182502 (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.

Comments: 12 pages, 11 figures

Report number: LLNL-JRNL-807174-DRAFT

Journal ref: Phys. Rev. C 102, 024616 (2020)

From bound states to the continuum

Authors: Calvin W. Johnson, Kristina D. Launey, Naftali Auerbach, Sonia Bacca, Bruce R. Barrett, Carl Brune, Mark A. Caprio, Pierre Descouvemont, W. H. Dickhoff, Charlotte Elster, Patrick J. Fasano, Kevin Fossez, Heiko Hergert, Morten Hjorth-Jensen, Linda Hlophe, Baishan Hu, Rodolfo M. Id Betan, Andrea Idini, Sebastian König, Konstantinos Kravvaris, Dean Lee, Jin Lei, Alexis Mercenne, Rodrigo Navarro Perez, Witold Nazarewicz , et al. (13 additional authors not shown)

Abstract: This white paper reports on the discussions of the 2018 Facility for Rare Isotope Beams Theory Alliance (FRIB-TA) topical program "From bound states to the continuum: Connecting bound state calculations with scattering and reaction theory". One of the biggest and most important frontiers in nuclear theory today is to construct better and stronger bridges between bound state calculations and calcul… ▽ More This white paper reports on the discussions of the 2018 Facility for Rare Isotope Beams Theory Alliance (FRIB-TA) topical program "From bound states to the continuum: Connecting bound state calculations with scattering and reaction theory". One of the biggest and most important frontiers in nuclear theory today is to construct better and stronger bridges between bound state calculations and calculations in the continuum, especially scattering and reaction theory, as well as teasing out the influence of the continuum on states near threshold. This is particularly challenging as many-body structure calculations typically use a bound state basis, while reaction calculations more commonly utilize few-body continuum approaches. The many-body bound state and few-body continuum methods use different language and emphasize different properties. To build better foundations for these bridges, we present an overview of several bound state and continuum methods and, where possible, point to current and possible future connections. △ Less

Submitted 16 November, 2020; v1 submitted 1 December, 2019; originally announced December 2019.

Comments: 48 pages, 23 figures; paper from the FRIB Theory Alliance Workshop on "From Bound States to the Continuum", June 11-22, 2018

Journal ref: J. Phys. G: Nucl. Part. Phys. 47 123001 (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.

Report number: LLNL-JRNL-787600

Journal ref: Phys. Rev. A 101, 062307 (2020)

$^7$Be and $^7$Li nuclei within the no-core shell model with continuum

Authors: Matteo Vorabbi, Petr Navrátil, Sofia Quaglioni, Guillaume Hupin

Abstract: The production of $^7$Be and $^7$Li nuclei plays an important role in primordial nucleosynthesis, nuclear astrophysics, and fusion energy generation. The $^3\mathrm{He}(αあるふぁ, γがんま) ^7\mathrm{Be}$ and $^3\mathrm{H}(αあるふぁ, γがんま) ^7\mathrm{Li}$ radiative-capture processes are important to determine the $^7$Li abundance in the early universe and to predict the correct fraction of pp-chain branches resulting in… ▽ More The production of $^7$Be and $^7$Li nuclei plays an important role in primordial nucleosynthesis, nuclear astrophysics, and fusion energy generation. The $^3\mathrm{He}(αあるふぁ, γがんま) ^7\mathrm{Be}$ and $^3\mathrm{H}(αあるふぁ, γがんま) ^7\mathrm{Li}$ radiative-capture processes are important to determine the $^7$Li abundance in the early universe and to predict the correct fraction of pp-chain branches resulting in $^7$Be versus $^8$B neutrinos. In this work we study the properties of $^7$Be and $^7$Li within the no-core shell model with continuum (NCSMC) method, using chiral nucleon-nucleon interactions as the only input, and analyze all the binary mass partitions involved in the formation of these systems. The NCSMC is an ab initio method applicable to light nuclei that provides a unified description of bound and scattering states and thus is well suited to investigate systems with many resonances and pronounced clustering like $^7$Be and $^7$Li. Our calculations reproduce all the experimentally known states of the two systems and provide predictions for several new resonances of both parities. Some of these new possible resonances are built on the ground states of $^6$Li and $^6$He, and thus represent a robust prediction. We do not find any resonance in the p${+}^6$Li mass partition near the threshold. On the other hand, in the p${+}^6$He mass partition of $^7$Li we observe an $S$-wave resonance near the threshold producing a very pronounced peak in the calculated S factor of the $^6\mathrm{He} (\mathrm{p},γがんま) ^7\mathrm{Li}$ radiative-capture reaction, which could be relevant for astrophysics and its implications should be investigated. △ Less

Submitted 21 June, 2019; originally announced June 2019.

Report number: LLNL-JRNL-777984-DRAFT

Journal ref: Phys. Rev. C 100, 024304 (2019)

Application of an ab-initio S-matrix to data analysis of transfer to the continuum reactions populating 11Be

Authors: A. Bonaccorso, F. Cappuzzello, D. Carbone, M. Cavallaro, G. Hupin, P. Navratil, S. Quaglioni

Abstract: Recently, the bound and continuum spectrum of 11Be has been calculated within the ab-initio no-core shell model with continuum (NCSMC) method successfully reproducing the parity inversion in the ground state. The continuum spectrum obtained is in agreement with known experimental levels. The S-matrix contained in the NCSMC continuum wave functions of the n+10Be system is used in this work for the… ▽ More Recently, the bound and continuum spectrum of 11Be has been calculated within the ab-initio no-core shell model with continuum (NCSMC) method successfully reproducing the parity inversion in the ground state. The continuum spectrum obtained is in agreement with known experimental levels. The S-matrix contained in the NCSMC continuum wave functions of the n+10Be system is used in this work for the first time in a Transfer-to-the-Continuum (TC) reaction calculation. The TC approach is applied to study the excitation energy spectrum of 11Be measured in the 9Be(18O,16O)11Be reaction at 84 MeV. Previously known levels are confirmed and theoretical and experimental evidence for a 9/2+ state at Ex=5.8 MeV is given, whose configuration is thought to be 10Be(2+)+n(d5/2). △ Less

Submitted 1 September, 2019; v1 submitted 30 April, 2019; originally announced April 2019.

Comments: 7 pages, 2 figures

Journal ref: Phys. Rev. C 100, 024617 (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.

Comments: 20 pages, 18 figures

Journal ref: Nature Physics 15, 428 (2019)

Ab initio predictions for polarized DT thermonuclear fusion

Authors: Guillaume Hupin, Sofia Quaglioni, Petr Navrátil

Abstract: The fusion of deuterium (D) with tritium (T) is the most promising of the reactions that could power the thermonuclear reactors of the future. Already favored for its low activation energy and high yield, it may lead to even more efficient energy generation if obtained in a polarized state, i.e. with the spin of the reactants aligned. While the DT fusion rate has been measured extensively, very li… ▽ More The fusion of deuterium (D) with tritium (T) is the most promising of the reactions that could power the thermonuclear reactors of the future. Already favored for its low activation energy and high yield, it may lead to even more efficient energy generation if obtained in a polarized state, i.e. with the spin of the reactants aligned. While the DT fusion rate has been measured extensively, very little is known of the effects of polarization. Meanwhile, arriving at a fundamental understanding of the fusion process in terms of the laws of quantum mechanics and the underlying theory of the strong force has been a daunting challenge. We use nuclear forces derived from chiral effective field theory and apply the ab initio reaction method known as no-core shell model with continuum to predict, for the first time from first principles, the enhancement factor of the polarized DT fusion rate and anisotropy of the emitted neutron and αあるふぁ particle. △ Less

Submitted 30 January, 2019; v1 submitted 30 March, 2018; originally announced March 2018.

Report number: LLNL-JRNL-745269-DRAFT

Journal ref: Nature Communications volume 10, Article number: 351 (2019)

Structure of the exotic $^9$He nucleus from the no-core shell model with continuum

Authors: Matteo Vorabbi, Angelo Calci, Petr Navratil, Michael K. G Kruse, Sofia Quaglioni, Guillaume Hupin

Abstract: The exotic $^9$He nucleus, which presents one of the most extreme neutron-to-proton ratios, belongs to the $N=7$ isotonic chain famous for the phenomenon of ground-state parity inversion with decreasing number of protons. Consequently, it would be expected to have an unnatural (positive) parity ground state similar to $^{11}$Be and $^{10}$Li. Despite many experimental and theoretical investigation… ▽ More The exotic $^9$He nucleus, which presents one of the most extreme neutron-to-proton ratios, belongs to the $N=7$ isotonic chain famous for the phenomenon of ground-state parity inversion with decreasing number of protons. Consequently, it would be expected to have an unnatural (positive) parity ground state similar to $^{11}$Be and $^{10}$Li. Despite many experimental and theoretical investigations, its structure remains uncertain. Apart from the fact that it is unbound, other properties including the spin and parity of its ground state and the very existence of additional low-lying resonances are still a matter of debate. In this work we study the properties of $^9$He by analyzing the $n+^8$He continuum in the context of the ab initio no-core shell model with continuum (NCSMC) formalism with chiral interactions as the only input. The NCSMC is a state-of-the-art approach for the ab initio description of light nuclei. With its capability to predict properties of bound states, resonances, and scattering states in a unified framework, the method is particularly well suited for the study of unbound nuclei such as $^9$He. Our analysis produces an unbound $^9$He nucleus. Two resonant states are found at the energies of ${\sim}1$ and ${\sim}3.5$ MeV, respectively, above the $n+^8$He breakup threshold. The first state has a spin-parity assignment of $J^πぱい = {1/2}^-$ and can be associated with the ground state of $^9$He, while the second, broader state has a spin-parity of ${3/2}^-$. No resonance is found in the ${1/2}^+$ channel, only a very weak attraction. We find that the $^9$He ground-state resonance has a negative parity and thus breaks the parity-inversion mechanism found in the $^{11}$Be and $^{10}$Li nuclei of the same $N=7$ isotonic chain. △ Less

Submitted 15 December, 2017; originally announced December 2017.

Report number: LLNL-JRNL-743601-DRAFT

Journal ref: Phys. Rev. C 97, 034314 (2018)

Three-cluster dynamics within the ab initio no-core shell model with continuum: How many-body correlations and $αあるふぁ$-clustering shape $^6$He

Authors: Sofia Quaglioni, Carolina Romero-Redondo, Petr Navratil, Guillaume Hupin

Abstract: We realize the treatment of bound and continuum nuclear systems in the proximity of a three-body breakup threshold within the ab initio framework of the no-core shell model with continuum. Many-body eigenstates obtained from the diagonalization of the Hamiltonian within the harmonic-oscillator expansion of the no-core shell model are coupled with continuous microscopic three-cluster states to corr… ▽ More We realize the treatment of bound and continuum nuclear systems in the proximity of a three-body breakup threshold within the ab initio framework of the no-core shell model with continuum. Many-body eigenstates obtained from the diagonalization of the Hamiltonian within the harmonic-oscillator expansion of the no-core shell model are coupled with continuous microscopic three-cluster states to correctly describe the nuclear wave function both in the interior and asymptotic regions. We discuss the formalism in detail and give algebraic expressions for the case of core+$n$+$n$ systems. Using similarity-renormalization-group evolved nucleon-nucleon interactions, we analyze the role of $^4$He+$n$+$n$ clustering and many-body correlations in the ground and low-lying continuum states of the Borromean $^6$He nucleus, and study the dependence of the energy spectrum on the resolution scale of the interaction. We show that $^6$He small binding energy and extended radii compatible with experiment can be obtained simultaneously, without recurring to extrapolations. We also find that a significant portion of the ground-state energy and the narrow width of the first $2^+$ resonance stem from many-body correlations that can be interpreted as core-excitation effects. △ Less

Submitted 19 October, 2017; originally announced October 2017.

Comments: 25 pages, 14 figures

Report number: LLNL-JRNL-739751-DRAFT

Journal ref: Phys. Rev. C 97, 034332 (2018)

The nuclear force imprints revealed on the elastic scattering of protons with $^{10}$C

Authors: A. Kumar, R. Kanungo, A. Calci, P. Navratil, A. Sanetullaev, M. Alcorta, V. Bildstein, G. Christian, B. Davids, J. Dohet-Eraly, J. Fallis, A. T. Gallant, G. Hackman, B. Hadinia, G. Hupin, S. Ishimoto, R. Krücken, A. T. Laffoley, J. Lighthall, D. Miller, S. Quaglioni, J. S. Randhawa, E. T. Rand, A. Rojas, R. Roth , et al. (4 additional authors not shown)

Abstract: How does nature hold together protons and neutrons to form the wide variety of complex nuclei in the universe? Describing many-nucleon systems from the fundamental theory of quantum chromodynamics has been the greatest challenge in answering this question. The chiral effective field theory description of the nuclear force now makes this possible but requires certain parameters that are not uniquel… ▽ More How does nature hold together protons and neutrons to form the wide variety of complex nuclei in the universe? Describing many-nucleon systems from the fundamental theory of quantum chromodynamics has been the greatest challenge in answering this question. The chiral effective field theory description of the nuclear force now makes this possible but requires certain parameters that are not uniquely determined. Defining the nuclear force needs identification of observables sensitive to the different parametrizations. From a measurement of proton elastic scattering on $^{10}$C at TRIUMF and ab initio nuclear reaction calculations we show that the shape and magnitude of the measured differential cross section is strongly sensitive to the nuclear force prescription. △ Less

Submitted 15 May, 2017; originally announced May 2017.

Comments: 7 pages, 6 figures

Report number: LLNL-JRNL-708097

${}^{7}$Li($d$,$p$)${}^{8}$Li transfer reaction in the NCSM/RGM approach

Authors: F. Raimondi, G. Hupin, P. Navrátil, S. Quaglioni

Abstract: Recently, we applied an $ab$ $initio$ method, the no-core shell model combined with the resonating group method, to the transfer reactions with light p-shell nuclei as targets and deuteron as the projectile. In particular, we studied the elastic scattering of deuterium on $^7$Li and the ${}^{7}$Li($d$,$p$)${}^{8}$Li transfer reaction starting from a realistic two-nucleon interaction. In this contr… ▽ More Recently, we applied an $ab$ $initio$ method, the no-core shell model combined with the resonating group method, to the transfer reactions with light p-shell nuclei as targets and deuteron as the projectile. In particular, we studied the elastic scattering of deuterium on $^7$Li and the ${}^{7}$Li($d$,$p$)${}^{8}$Li transfer reaction starting from a realistic two-nucleon interaction. In this contribution, we review of our main results on the ${}^{7}$Li($d$,$p$)${}^{8}$Li transfer reaction, and we extend the study of the relevant reaction channels, by showing the dominant resonant phase shifts of the scattering matrix. We assess also the impact of the polarization effects of the deuteron below the breakup on the positive-parity resonant states in the reaction. For this purpose, we perform an analysis of the convergence trend of the phase and eigenphase shifts, with respect to the number of deuteron pseudostates included in the model space. △ Less

Submitted 15 February, 2017; originally announced February 2017.

Comments: Contribution to TNPI2016

Can Ab Initio Theory Explain the Phenomenon of Parity Inversion in ${}^{11}$Be?

Authors: Angelo Calci, Petr Navrátil, Robert Roth, Jérémy Dohet-Eraly, Sofia Quaglioni, Guillaume Hupin

Abstract: The weakly bound exotic ${}^{11}$Be nucleus, famous for its ground-state parity inversion and distinct n+ ${}^{10}$Be halo structure, is investigated from first principles using chiral two- and three-nucleon forces. An explicit treatment of continuum effects is found to be indispensable. We study the sensitivity of the ${}^{11}$Be spectrum to the details of the three-nucleon force and demonstrate… ▽ More The weakly bound exotic ${}^{11}$Be nucleus, famous for its ground-state parity inversion and distinct n+ ${}^{10}$Be halo structure, is investigated from first principles using chiral two- and three-nucleon forces. An explicit treatment of continuum effects is found to be indispensable. We study the sensitivity of the ${}^{11}$Be spectrum to the details of the three-nucleon force and demonstrate that only certain chiral interactions are capable of reproducing the parity inversion. With such interactions, the extremely large E1 transition between the bound states is reproduced. We compare our photodisintegration calculations to conflicting experimental data and predict a distinct dip around the $3/2^-_1$ resonance energy. Finally, we predict low-lying $3/2^+$ and $9/2^+$ resonances that are not or not sufficiently measured in experiments. △ Less

Submitted 30 November, 2016; v1 submitted 10 August, 2016; originally announced August 2016.

Comments: 6 pages, 5 figures, plus Supplemental Material

Report number: LLNL-JRNL-833846

Journal ref: Phys. Rev. Lett. 117, 242501 (2016)

How many-body correlations and $αあるふぁ$-clustering shape $^6$He

Authors: Carolina Romero-Redondo, Sofia Quaglioni, Petr Navratil, Guillaume Hupin

Abstract: The Borromean $^6$He nucleus is an exotic system characterized by two `halo' neutrons orbiting around a compact $^4$He (or $αあるふぁ$) core, in which the binary subsystems are unbound. The simultaneous reproduction of its small binding energy and extended matter and point-proton radii has been a challenge for {\em ab initio} theoretical calculations based on traditional bound-state methods. Using soft nu… ▽ More The Borromean $^6$He nucleus is an exotic system characterized by two `halo' neutrons orbiting around a compact $^4$He (or $αあるふぁ$) core, in which the binary subsystems are unbound. The simultaneous reproduction of its small binding energy and extended matter and point-proton radii has been a challenge for {\em ab initio} theoretical calculations based on traditional bound-state methods. Using soft nucleon-nucleon interactions based on chiral effective field theory potentials, we show that supplementing the model space with $^4$He+$n$+$n$ cluster degrees of freedom largely solves this issue. We analyze the role played by the $αあるふぁ$-clustering and many-body correlations, and study the dependence of the energy spectrum on the resolution scale of the interaction. △ Less

Submitted 31 May, 2016; originally announced June 2016.

Comments: 5 pages, 4 figures

Report number: LLNL-JRNL-693078-DRAFT

Journal ref: Phys. Rev. Lett. 117, 222501 (2016)

Deuteron-induced nucleon transfer reactions within an ab initio framework: First application to p-shell nuclei

Authors: Francesco Raimondi, Guillaume Hupin, Petr Navrátil, Sofia Quaglioni

Abstract: Background: Low-energy transfer reactions in which a proton is stripped from a deuteron projectile and dropped into a target play a crucial role in the formation of nuclei in both primordial and stellar nucleosynthesis, as well as in the study of exotic nuclei using radioactive beam facilities and inverse kinematics. Ab initio approaches have been successfully applied to describe the $^3$H… ▽ More Background: Low-energy transfer reactions in which a proton is stripped from a deuteron projectile and dropped into a target play a crucial role in the formation of nuclei in both primordial and stellar nucleosynthesis, as well as in the study of exotic nuclei using radioactive beam facilities and inverse kinematics. Ab initio approaches have been successfully applied to describe the $^3$H$(d,n)^4$He and $^3$He$(d,p)^4$He fusion processes. Purpose: An ab initio treatment of transfer reactions would also be desirable for heavier targets. In this work, we extend the ab initio description of $(d,p)$ reactions to processes with light $p$-shell nuclei. As a first application, we study the elastic scattering of deuterium on $^7$Li and the ${}^{7}$Li($d$,$p$)${}^{8}$Li transfer reaction based on a two-body Hamiltonian. Methods: We use the no-core shell model to compute the wave functions of the nuclei involved in the reaction, and describe the dynamics between targets and projectiles with the help of microscopic-cluster states in the spirit of the resonating group method. Results: The shape of the excitation functions for deuteron impinging on ${}^{7}$Li are qualitatively reproduced up to the deuteron breakup energy. The interplay between $d$-$^7$Li and $p$-$^8$Li particle-decay channels determines some features of the ${}^{9}$Be spectrum above the $d$+${}^{7}$Li threshold. Our prediction for the parity of the 17.298 MeV resonance is at odds with the experimental assignment Conclusions: Deuteron stripping reactions with $p$-shell targets can now be computed ab initio, but calculations are very demanding. A quantitative description of the ${}^{7}$Li($d$,$p$)${}^{8}$Li reaction will require further work to include the effect of three-nucleon forces and additional decay channels, and improve the convergence rate of our calculations. △ Less

Submitted 13 February, 2016; originally announced February 2016.

Comments: 16 pages, 11 figures, submitted to Phys. Rev. C

Report number: LLNL-JRNL-682741-DRAF

Unified ab initio approaches to nuclear structure and reactions

Authors: Petr Navratil, Sofia Quaglioni, Guillaume Hupin, Carolina Romero-Redondo, Angelo Calci

Abstract: The description of nuclei starting from the constituent nucleons and the realistic interactions among them has been a long-standing goal in nuclear physics. In addition to the complex nature of the nuclear forces, with two-, three- and possibly higher many-nucleon components, one faces the quantum-mechanical many-nucleon problem governed by an interplay between bound and continuum states. In recen… ▽ More The description of nuclei starting from the constituent nucleons and the realistic interactions among them has been a long-standing goal in nuclear physics. In addition to the complex nature of the nuclear forces, with two-, three- and possibly higher many-nucleon components, one faces the quantum-mechanical many-nucleon problem governed by an interplay between bound and continuum states. In recent years, significant progress has been made in ab initio nuclear structure and reaction calculations based on input from QCD-employing Hamiltonians constructed within chiral effective field theory. After a brief overview of the field, we focus on ab initio many-body approaches - built upon the No-Core Shell Model - that are capable of simultaneously describing both bound and scattering nuclear states, and present results for resonances in light nuclei, reactions important for astrophysics and fusion research. In particular, we review recent calculations of resonances in the $^6$He halo nucleus, of five- and six-nucleon scattering, and an investigation of the role of chiral three-nucleon interactions in the structure of $^9$Be. Further, we discuss applications to the $^7$Be$(p,γがんま)^8$B radiative capture. Finally, we highlight our efforts to describe transfer reactions including the $^3$H$(d,n)^4$He fusion. △ Less

Submitted 21 April, 2016; v1 submitted 14 January, 2016; originally announced January 2016.

Comments: Contribution to the Special Physica Scripta Edition - 40 year anniversary - Nobel Prize '75, 71 pages, 29 figures

Report number: LLNL-JRNL-680282

Journal ref: Physica Scripta 91 (2016) 053002

${^3{\rm He}}(αあるふぁ,γがんま){^7{\rm Be}}$ and ${^3{\rm H}}(αあるふぁ,γがんま){^7{\rm Li}}$ astrophysical $S$ factors from the no-core shell model with continuum

Authors: Jérémy Dohet-Eraly, Petr Navrátil, Sofia Quaglioni, Wataru Horiuchi, Guillaume Hupin, Francesco Raimondi

Abstract: The ${^3{\rm He}}(αあるふぁ,γがんま){^7{\rm Be}}$ and ${^3{\rm H}}(αあるふぁ,γがんま){^7{\rm Li}}$ astrophysical $S$ factors are calculated within the no-core shell model with continuum using a renormalized chiral nucleon-nucleon interaction. The ${^3{\rm He}}(αあるふぁ,γがんま){^7{\rm Be}}$ astrophysical $S$ factors agree reasonably well with the experimental data while the ${^3{\rm H}}(αあるふぁ,γがんま){^7{\rm Li}}$ ones are overestimated. The seven… ▽ More The ${^3{\rm He}}(αあるふぁ,γがんま){^7{\rm Be}}$ and ${^3{\rm H}}(αあるふぁ,γがんま){^7{\rm Li}}$ astrophysical $S$ factors are calculated within the no-core shell model with continuum using a renormalized chiral nucleon-nucleon interaction. The ${^3{\rm He}}(αあるふぁ,γがんま){^7{\rm Be}}$ astrophysical $S$ factors agree reasonably well with the experimental data while the ${^3{\rm H}}(αあるふぁ,γがんま){^7{\rm Li}}$ ones are overestimated. The seven-nucleon bound and resonance states and the $αあるふぁ+{^3{\rm He}}/{^3{\rm H}}$ elastic scattering are also studied and compared with experiment. The low-lying resonance properties are rather well reproduced by our approach. At low energies, the $s$-wave phase shift, which is non-resonant, is overestimated. △ Less

Submitted 13 February, 2016; v1 submitted 26 October, 2015; originally announced October 2015.

Comments: 8 pages, submitted to Phys. Lett. B

Report number: LLNL-JRNL-678557

Ab initio calculations of reactions with light nuclei

Authors: S. Quaglioni, G. Hupin, A. Calci, P. Navratil, R. Roth

Abstract: An {\em ab initio} (i.e., from first principles) theoretical framework capable of providing a unified description of the structure and low-energy reaction properties of light nuclei is desirable to further our understanding of the fundamental interactions among nucleons, and provide accurate predictions of crucial reaction rates for nuclear astrophysics, fusion-energy research, and other applicati… ▽ More An {\em ab initio} (i.e., from first principles) theoretical framework capable of providing a unified description of the structure and low-energy reaction properties of light nuclei is desirable to further our understanding of the fundamental interactions among nucleons, and provide accurate predictions of crucial reaction rates for nuclear astrophysics, fusion-energy research, and other applications. In this contribution we review {\em ab initio} calculations for nucleon and deuterium scattering on light nuclei starting from chiral two- and three-body Hamiltonians, obtained within the framework of the {\em ab initio} no-core shell model with continuum. This is a unified approach to nuclear bound and scattering states, in which square-integrable energy eigenstates of the $A$-nucleon system are coupled to $(A-a)+a$ target-plus-projectile wave functions in the spirit of the resonating group method to obtain an efficient description of the many-body nuclear dynamics both at short and medium distances and at long ranges. △ Less

Submitted 30 September, 2015; originally announced September 2015.

Comments: 9 pages, 5 figures, proceedings of the 21st International Conference on Few-Body Problems in Physics

Report number: LLNL-PROC-676805

Towards an ab initio description of the light-nuclei radiative captures

Authors: Jérémy Dohet-Eraly, Petr Navrátil, Sofia Quaglioni, Wataru Horiuchi, Guillaume Hupin

Abstract: The ${^3{\rm He}}(αあるふぁ,γがんま){^7{\rm Be}}$ and ${^3{\rm H}}(αあるふぁ,γがんま){^7{\rm Li}}$ astrophysical $S$ factors are evaluated at low collision energies (less than 2.5 MeV in the centre-of-mass frame) within the no-core shell model with continuum approach using a renormalized chiral nucleon-nucleon interaction. The ${^3{\rm He}}(αあるふぁ,γがんま){^7{\rm Be}}$ and ${^3{\rm H}}(αあるふぁ,γがんま){^7{\rm Li}}$ astrophysical $S$ factors are evaluated at low collision energies (less than 2.5 MeV in the centre-of-mass frame) within the no-core shell model with continuum approach using a renormalized chiral nucleon-nucleon interaction. △ Less

Submitted 23 October, 2015; v1 submitted 3 September, 2015; originally announced September 2015.

Comments: 4 pages, submitted as a proceeding of the 21st International Conference on Few-Body Problems in Physics

Report number: LLNL-PROC-676838

Advances in the ab initio description of nuclear three-cluster systems

Authors: Carolina Romero-Redondo, Sofia Quaglioni, Petr Navrátil, Guillaume Hupin

Abstract: We introduce the extension of the ab initio no-core shell model with continuum to describe three-body cluster systems. We present results for the ground state of 6He and show improvements with respect to the description obtained within the no-core shell model and the no-core shell model/resonating group methods. We introduce the extension of the ab initio no-core shell model with continuum to describe three-body cluster systems. We present results for the ground state of 6He and show improvements with respect to the description obtained within the no-core shell model and the no-core shell model/resonating group methods. △ Less

Submitted 2 September, 2015; originally announced September 2015.

Comments: Proceedings of the 21st International Conference on Few-Body Problems in Physics. May 18-22, 2015. Chicago, Illinois, USA

Report number: LLNL-PROC-676587

Microscopic Study of $αあるふぁ+N$ Bremsstrahlung from Effective and Realistic Multi-nucleon Interactions

Authors: Jérémy Dohet-Eraly, Sofia Quaglioni, Petr Navrátil, Guillaume Hupin

Abstract: Based on an effective nucleon-nucleon interaction, a microscopic cluster model of the nucleus-nucleus bremsstrahlung, including implicitly a part of the effects of meson-exchange currents via an extension of the Siegert theorem is applied to the $αあるふぁ+p$ and $αあるふぁ+n$ systems. The contributions of the $E1$ and $E2$ transitions to the bremsstrahlung cross sections are evaluated and their relative importan… ▽ More Based on an effective nucleon-nucleon interaction, a microscopic cluster model of the nucleus-nucleus bremsstrahlung, including implicitly a part of the effects of meson-exchange currents via an extension of the Siegert theorem is applied to the $αあるふぁ+p$ and $αあるふぁ+n$ systems. The contributions of the $E1$ and $E2$ transitions to the bremsstrahlung cross sections are evaluated and their relative importance for the mirror systems $αあるふぁ+p$ and $αあるふぁ+n$ is compared. Another approach based on realistic two- and three-nucleon interactions and the No-Core Shell Model/Resonating-Group Method is also investigated. Some preliminary results for the $αあるふぁ+p$ bremsstrahlung are displayed. △ Less

Submitted 12 January, 2015; originally announced January 2015.

Comments: Contribution to the proceedings of the Conference on "Advances in Radioactive Isotope Science" (ARIS2014). 4 pages, 2 figures

Report number: LLNL-PROC-663865

Unified description of $^6$Li structure and deuterium-$^4$He dynamics with chiral two- and three-nucleon forces

Authors: Guillaume Hupin, Sofia Quaglioni, Petr Navrátil

Abstract: Prototype for the study of weakly bound projectiles colliding on stable targets, the scattering of deuterium ($d$) on $^4$He ($αあるふぁ$) is an important milestone in the search for a fundamental understanding of low-energy reactions. At the same time, it is also important for its role in the Big-bang nucleosynthesis of $^6$Li and applications in the characterization of deuterium impurities in materials.… ▽ More Prototype for the study of weakly bound projectiles colliding on stable targets, the scattering of deuterium ($d$) on $^4$He ($αあるふぁ$) is an important milestone in the search for a fundamental understanding of low-energy reactions. At the same time, it is also important for its role in the Big-bang nucleosynthesis of $^6$Li and applications in the characterization of deuterium impurities in materials. We present the first unified {\em ab initio} study of the $^6$Li ground state and $d$-$^4$He elastic scattering using two- and three-nucleon forces derived within the framework of chiral effective field theory. The six-nucleon bound-state and scattering observables are calculated by means of the no-core shell model with continuum. %and are compared to available experimental data. We analyze the influence of the dynamic polarization of the deuterium and of the chiral three-nucleon force, and examine the role of the continuum degrees of freedom in shaping the low-lying spectrum of $^6$Li. We find that the adopted Hamiltonian correctly predicts the binding energy of $^6$Li, yielding an asymptotic $D$- to $S$-state ratio of the $^6$Li wave function in $d+αあるふぁ$ configuration of $-0.027$ in agreement with the value determined from a phase shift analysis of $^6$Li+$^4$He elastic scattering, but overestimates the excitation energy of the first $3^+$ state by $350$ keV. The bulk of the computed differential cross section is in good agreement with data. △ Less

Submitted 12 December, 2014; originally announced December 2014.

Comments: 5 pages, 5 figures

Report number: LLNL-JRNL-665027

Journal ref: Phys. Rev. Lett. 114, 212502 (2015)

Continuum and Three-Nucleon Force Effects on 9Be Energy Levels

Authors: Joachim Langhammer, Petr Navratil, Sofia Quaglioni, Guillaume Hupin, Angelo Calci, Robert Roth

Abstract: We extend the recently proposed ab initio no-core shell model with continuum to include three-nucleon (3N) interactions beyond the few-body domain. The extended approach allows for the assessment of effects of continuum degrees of freedom as well as of the 3N force in ab initio calculations of structure and reaction observables of p- and lower-sd-shell nuclei. As first application we concentrate o… ▽ More We extend the recently proposed ab initio no-core shell model with continuum to include three-nucleon (3N) interactions beyond the few-body domain. The extended approach allows for the assessment of effects of continuum degrees of freedom as well as of the 3N force in ab initio calculations of structure and reaction observables of p- and lower-sd-shell nuclei. As first application we concentrate on energy levels of the 9Be system for which all excited states lie above the n-8Be threshold. For all energy levels, the inclusion of the continuum significantly improves the agreement with experiment, which was an issue in standard no-core shell model calculations. Furthermore, we find the proper treatment of the continuum indispensable for reliable statements about the quality of the adopted 3N interaction from chiral effective field theory. In particular, we find the 1/2+ resonance energy, which is of astrophysical interest, in good agreement with experiment. △ Less

Submitted 10 November, 2014; originally announced November 2014.

Comments: 6 pages, 5 figures

Report number: LLNL-JRNL-663595

A Predictive Theory for Elastic Scattering and Recoil of Protons from $^4$He

Authors: Guillaume Hupin, Sofia Quaglioni, Petr Navrátil

Abstract: Low-energy cross sections for elastic scattering and recoil of protons from $^4$He nuclei (also known as $αあるふぁ$ particles) are calculated directly by solving the Schrödinger equation for five nucleons interacting through accurate two- and three-nucleon forces derived within the framework of chiral effective field theory. Precise knowledge of these processes at various proton backscattering/recoil ang… ▽ More Low-energy cross sections for elastic scattering and recoil of protons from $^4$He nuclei (also known as $αあるふぁ$ particles) are calculated directly by solving the Schrödinger equation for five nucleons interacting through accurate two- and three-nucleon forces derived within the framework of chiral effective field theory. Precise knowledge of these processes at various proton backscattering/recoil angles and energies is needed for the ion-beam analysis of numerous materials, from the surface layers of solids, to thin films, to fusion-reactor materials. Indeed, the same elastic scattering process, in two different kinematic configurations, can be used to probe concentrations and depth profiles of either hydrogen or helium. We compare our results to available experimental data and show that direct calculations with modern nuclear potentials can help to resolve remaining inconsistencies among different data sets and can be used to predict these cross sections when measurements are not available. △ Less

Submitted 2 September, 2014; originally announced September 2014.

Comments: 5 pages, 5 figures

Report number: LLNL-JRNL-659467

Journal ref: Phys. Rev. C 90, 061601 (2014)

$^4{\rm He}$+$n$+$n$ continuum within an ab initio framework

Authors: Carolina Romero-Redondo, Sofia Quaglioni, Petr Navrátil, Guillaume Hupin

Abstract: The low-lying continuum spectrum of the $^6{\rm He}$ nucleus is investigated for the first time within an ab initio framework that encompasses the $^4{\rm He}$+$n$+$n$ three-cluster dynamics characterizing its lowest decay channel. This is achieved through an extension of the no-core-shell model combined with the resonating-group method, in which energy-independent non-local interactions among thr… ▽ More The low-lying continuum spectrum of the $^6{\rm He}$ nucleus is investigated for the first time within an ab initio framework that encompasses the $^4{\rm He}$+$n$+$n$ three-cluster dynamics characterizing its lowest decay channel. This is achieved through an extension of the no-core-shell model combined with the resonating-group method, in which energy-independent non-local interactions among three nuclear fragments can be calculated microscopically starting from realistic nucleon-nucleon interactions and consistent ab initio many-body wave functions of the clusters. The three-cluster Schrödinger equation is solved with three-body scattering boundary conditions by means of the hyperspherical-harmonic method on a Lagrange mesh. Using a soft similarity-renormalization-group evolved chiral nucleon-nucleon potential, we find the known $J^πぱい= 2^+$ resonance as well as a result consistent with a new low-lying second $2^+$ resonance recently observed at GANIL at $2.6$ MeV above the $^6$He ground state. We also find resonances in the $2^-$, $1^+$ and $0^-$ channels, while no low-lying resonances are present in the $0^+$ and $1^-$ channels. △ Less

Submitted 7 April, 2014; originally announced April 2014.

Comments: 5 pages 4 figures

Report number: LLNL-JRNL-652701

Journal ref: Phys. Rev. Lett. 113, 032503 (2014)

Operator evolution for ab initio nuclear theory

Authors: Micah D Schuster, Sofia Quaglioni, Calvin W. Johnson, Eric D. Jurgenson, Petr Navratil

Abstract: The past two decades have seen a revolution in ab initio calculations of nuclear properties. One key element has been the development of a rigorous effective interaction theory, applying unitary transformations to soften the nuclear Hamiltonian and hence accelerate the convergence as a function of the model space size. For consistency, however, one ought to apply the same transformation to other o… ▽ More The past two decades have seen a revolution in ab initio calculations of nuclear properties. One key element has been the development of a rigorous effective interaction theory, applying unitary transformations to soften the nuclear Hamiltonian and hence accelerate the convergence as a function of the model space size. For consistency, however, one ought to apply the same transformation to other operators when calculating transitions and mean values from the eigenstates of the renormalized Hamiltonian. Working in a translationally-invariant harmonic oscillator basis for the two- and three-nucleon systems, we evolve the Hamiltonian, square-radius and total dipole strength operators by the similarity renormalization group (SRG). The inclusion of up to three-body matrix elements in the 4He nucleus all but completely restores the invariance of the expectation values under the transformation. We also consider a Gaussian operator with adjustable range and find at short ranges an increased contribution from such induced three-body terms. △ Less

Submitted 27 February, 2014; originally announced February 2014.

Comments: 4 pages, 3 figures

Report number: LLNL-JRNL-650041

Journal ref: Phys. Rev. C 90, 011301(R) (2014)

Progress on Light-Ion Fusion Reactions with Three-Nucleon Forces

Authors: Guillaume Hupin, Sofia Quaglioni, Joachim Langhammer, Petr Navrátil, Angelo Calci, Robert Roth

Abstract: The description of structural and dynamical properties of nuclei starting from the fundamental interaction between nucleons has been a long-standing goal in nuclear physics. The ab initio No-Core Shell Model combined with the Resonating-Group Method (NCSM/RGM) is capable of addressing both structural and reaction properties of light-nuclei. While promising results have already been achieved starti… ▽ More The description of structural and dynamical properties of nuclei starting from the fundamental interaction between nucleons has been a long-standing goal in nuclear physics. The ab initio No-Core Shell Model combined with the Resonating-Group Method (NCSM/RGM) is capable of addressing both structural and reaction properties of light-nuclei. While promising results have already been achieved starting from a two-body Hamiltonian, a truly realistic prediction of nuclear observables requires the treatment of the three-nucleon interaction. Using similarity-renormalization-group evolved two- and three-nucleon interactions, we will present recent applications to n-4He scattering process when accounting for the chiral two- plus three-nucleon interaction versus the chiral two-nucleon interaction. We compare our results to phase shifts obtained from R-matrix analysis of data up to 16 MeV neutron energy, below the d-3H threshold. △ Less

Submitted 1 January, 2014; originally announced January 2014.

Comments: 4 pages, 2 figures, Proceedings of the 22nd European Conference on Few-Body Problems in Physics, Krakow, Poland, 9-13 September 2013

Report number: LLNL-PROC-645788

Ab initio NCSM/RGM for three-body cluster systems and application to $^4$He+n+n

Authors: Carolina Romero-Redondo, Petr Navratil, Sofia Quaglioni, Guillaume Hupin

Abstract: We introduce an extension of the ab initio no-core shell model/resonating group method (NCSM/RGM) in order to describe three-body cluster states. We present results for the $^6$He ground state within a $^4$He+n+n cluster basis as well as first results for the phase shifts of different channels of the $^4$He+n+n system which provide information about low-lying resonances of this nucleus. We introduce an extension of the ab initio no-core shell model/resonating group method (NCSM/RGM) in order to describe three-body cluster states. We present results for the $^6$He ground state within a $^4$He+n+n cluster basis as well as first results for the phase shifts of different channels of the $^4$He+n+n system which provide information about low-lying resonances of this nucleus. △ Less

Submitted 18 November, 2013; originally announced November 2013.

Comments: 4 pages, 3 figures, Proceedings of the 22nd European Conference on Few-Body Problems in Physics, Krakow, Poland, 9 - 13 September 2013

Report number: LLNL-PROC-652666

Ab initio many-body calculations of nucleon-4He scattering with three-nucleon forces

Authors: Guillaume Hupin, Joachim Langhammer, Petr Navrátil, Sofia Quaglioni, Angelo Calci, Robert Roth

Abstract: We extend the ab initio no-core shell model/resonating-group method to include three-nucleon (3N) interactions for the description of nucleon-nucleus collisions. We outline the formalism, give algebraic expressions for the 3N-force integration kernels, and discuss computational aspects of two alternative implementations. The extended theoretical framework is then applied to nucleon-4He scattering… ▽ More We extend the ab initio no-core shell model/resonating-group method to include three-nucleon (3N) interactions for the description of nucleon-nucleus collisions. We outline the formalism, give algebraic expressions for the 3N-force integration kernels, and discuss computational aspects of two alternative implementations. The extended theoretical framework is then applied to nucleon-4He scattering using similarity-renormalization-group (SRG) evolved nucleon-nucleon plus three-nucleon potentials derived from chiral effective field theory. We analyze the convergence properties of the calculated phase shifts and explore their dependence upon the SRG evolution parameter. We include up to six excited states of the 4He target and find significant effects from the inclusion of the chiral 3N force, e.g., it enhances the spin-orbit splitting between the 3/2- and 1/2- resonances and leads to an improved agreement with the phase shifts obtained from an accurate R-matrix analysis of the five-nucleon experimental data. We find remarkably good agreement with measured differential cross sections at various energies, while analyzing powers manifest larger deviations from experiment for certain energies and angles. △ Less

Submitted 12 August, 2013; originally announced August 2013.

Comments: 18 pages, 12 figures

Report number: LLNL-JRNL-641544