-
16O electroweak response functions from first principles
Authors:
Bijaya Acharya,
Joanna E. Sobczyk,
Sonia Bacca,
Gaute Hagen,
Weiguang Jiang
Abstract:
We present calculations of various electroweak response functions for the 16O nucleus obtained using coupled-cluster theory in conjunction with the Lorentz integral transform method. We employ nuclear forces derived at next-to-leading order and next-to-next-to-leading order in chiral effective field theory and perform a Bayesian analysis to assess uncertainties. Our results are in good agreement w…
▽ More
We present calculations of various electroweak response functions for the 16O nucleus obtained using coupled-cluster theory in conjunction with the Lorentz integral transform method. We employ nuclear forces derived at next-to-leading order and next-to-next-to-leading order in chiral effective field theory and perform a Bayesian analysis to assess uncertainties. Our results are in good agreement with available electron-scattering data at q~ 326 MeV/c. Additionally, we provide several predictions for the weak response functions in the quasi-elastic peak region at q= 300 and 400 MeV/c, which are critical for long-baseline neutrino experiments.
△ Less
Submitted 8 October, 2024;
originally announced October 2024.
-
Electric dipole polarizability of $^{58}$Ni
Authors:
I. Brandherm,
F. Bonaiti,
P. von Neumann-Cosel,
S. Bacca,
G. Colò,
G. R. Jansen,
Z. Z. Li,
H. Matsubara,
Y. F. Niu,
P. -G. Reinhard,
A. Richter,
X. Roca-Maza,
A. Tamii
Abstract:
The electric dipole strength distribution in $^{58}$Ni between 6 and 20 MeV has been determined from proton inelastic scattering experiments at very forward angles at RCNP, Osaka. The experimental data are rather well reproduced by quasiparticle random-phase approximation calculations including vibration coupling, despite a mild dependence on the adopted Skyrme interaction. They allow an estimate…
▽ More
The electric dipole strength distribution in $^{58}$Ni between 6 and 20 MeV has been determined from proton inelastic scattering experiments at very forward angles at RCNP, Osaka. The experimental data are rather well reproduced by quasiparticle random-phase approximation calculations including vibration coupling, despite a mild dependence on the adopted Skyrme interaction. They allow an estimate of the experimentally inaccessible high-energy contribution above 20 MeV, leading to an electric dipole polarizability $α_\mathrm{D}(^{58}{\rm Ni}) = 3.48(31)$ fm$^3$. This serves as a test case for recent extensions of coupled-cluster calculations with chiral effective field theory interactions to nuclei with two nucleons on top of a closed-shell system.
△ Less
Submitted 1 October, 2024;
originally announced October 2024.
-
Measurement of the $\mathrm{{}^{12}C}(e,e')$ cross sections at $Q^2=0.8\,\mathrm{GeV}^2/c^2$
Authors:
M. Mihovilovič,
L. Doria,
P. Achenbach,
A. M. Ankowski,
S. Bacca,
D. Bosnar,
A. Denig,
M. O. Distler,
A. Esser,
I. Friščić,
C. Giusti,
M. Hoek,
S. Kegel,
M. Littich,
G. D. Megias,
H. Merkel,
U. Muller,
J. Pochodzalla,
B. S. Schlimme,
M. Schoth,
C. Sfienti,
S. Širca,
J. E. Sobczyk,
Y. Stottinger,
M. Thiel
Abstract:
We present the findings of a study based on a new inelastic electron-scattering experiment on the ${}^{12}\mathrm{C}$ nucleus focusing on the kinematic region of $Q^2=0.8\,\mathrm{GeV}^2/{c}^2$. The measured cross section is sensitive to the transverse response function and provides a stringent test of theoretical models, as well as of the theoretical assumptions made in Monte-Carlo event-generato…
▽ More
We present the findings of a study based on a new inelastic electron-scattering experiment on the ${}^{12}\mathrm{C}$ nucleus focusing on the kinematic region of $Q^2=0.8\,\mathrm{GeV}^2/{c}^2$. The measured cross section is sensitive to the transverse response function and provides a stringent test of theoretical models, as well as of the theoretical assumptions made in Monte-Carlo event-generator codes developed for the interpretation of neutrino-nucleus experiments, such as DUNE and HyperK. We find that modern generators such as GENIE and GiBUU reproduce our new experimental data within 10$\%$.
△ Less
Submitted 23 June, 2024;
originally announced June 2024.
-
Electromagnetic observables of open-shell nuclei from coupled-cluster theory
Authors:
Francesca Bonaiti,
Sonia Bacca,
Gaute Hagen,
Gustav R. Jansen
Abstract:
We develop a new method to describe electromagnetic observables of open-shell nuclei with two nucleons outside a closed shell. This approach combines the equation-of-motion coupled-cluster method for such systems and the Lorentz integral transform technique, expanding the applicability of coupled-cluster theory for these properties beyond closed-shell nuclei. To validate this new approach, we comp…
▽ More
We develop a new method to describe electromagnetic observables of open-shell nuclei with two nucleons outside a closed shell. This approach combines the equation-of-motion coupled-cluster method for such systems and the Lorentz integral transform technique, expanding the applicability of coupled-cluster theory for these properties beyond closed-shell nuclei. To validate this new approach, we compute the non-energy-weighted dipole sum rule and the dipole polarizability of $^{16,24}$O in both the closed-shell and the new equation-of-motion coupled-cluster frameworks, finding agreement within error bars. We then analyze the evolution of the dipole polarizability along the oxygen and calcium isotopic chains. Our predictions agree well with available experimental data and other available theoretical calculations for the closed-shell $^{16,22}$O and the open-shell $^{18}$O. In the calcium isotopes, we observe that our dipole polarizability predictions for open-shell nuclei are lower than those of closed-shell nuclei. Our predictions for $^{24}$O and $^{54,56}$Ca will motivate future experimental studies at the dripline.
△ Less
Submitted 11 October, 2024; v1 submitted 9 May, 2024;
originally announced May 2024.
-
Low-energy dipole strength in 8He
Authors:
Francesca Bonaiti,
Sonia Bacca
Abstract:
In this work, we present new ab initio coupled-cluster calculations of dipole-excited state properties of 8He based on the chiral effective field theory interaction 1.8/2.0 (EM). We focus on the dipole polarizability, and compare the results to our previous study [Phys. Rev. C 105, 034313 (2022)] and subsequent theoretical work. With the aim of connecting the presence of low-lying dipole strength…
▽ More
In this work, we present new ab initio coupled-cluster calculations of dipole-excited state properties of 8He based on the chiral effective field theory interaction 1.8/2.0 (EM). We focus on the dipole polarizability, and compare the results to our previous study [Phys. Rev. C 105, 034313 (2022)] and subsequent theoretical work. With the aim of connecting the presence of low-lying dipole strength to structure properties of 8He, we compute the point-neutron radius, finding excellent agreement with available experimental data, and investigate its correlation with the dipole polarizability.
△ Less
Submitted 20 February, 2024;
originally announced February 2024.
-
Revisiting the helium isotope-shift puzzle with improved uncertainties from nuclear structure corrections
Authors:
Simone Salvatore Li Muli,
Thomas R. Richardson,
Sonia Bacca
Abstract:
Measurements of the difference between the squared charge radii of the helion ($^3$He nucleus) and the $α$-particle ($^4$He nucleus) have been characterized by longstanding tensions, recently spotlighted in the 3.6 $σ$ discrepancy of the extractions from ordinary atoms versus those from muonic atoms. Here, we present a novel analysis of uncertainties in nuclear structure corrections that must be s…
▽ More
Measurements of the difference between the squared charge radii of the helion ($^3$He nucleus) and the $α$-particle ($^4$He nucleus) have been characterized by longstanding tensions, recently spotlighted in the 3.6 $σ$ discrepancy of the extractions from ordinary atoms versus those from muonic atoms. Here, we present a novel analysis of uncertainties in nuclear structure corrections that must be supplied by theory to enable the extraction of the difference in radii from spectroscopic experiments. We use modern Bayesian inference techniques to quantify uncertainties stemming from the truncation of the chiral effective field theory expansion of the nuclear force for both muonic and ordinary atoms. With the new nuclear structure input, the helium isotope-shift puzzle cannot be explained, rather it is reinforced to a 4 $σ$ discrepancy.
△ Less
Submitted 24 January, 2024;
originally announced January 2024.
-
Response functions and giant monopole resonances for light to medium-mass nuclei from the \textit{ab initio} symmetry-adapted no-core shell model
Authors:
M. Burrows,
R. B. Baker,
S. Bacca,
K. D. Launey,
T. Dytrych,
D. Langr
Abstract:
Using the \textit{ab initio} symmetry-adapted no-core shell model, we compute sum rules and response functions for light to medium-mass nuclei, starting from interactions that are derived in the chiral effective field theory. We investigate electromagnetic transitions of monopole, dipole and quadrupole nature for symmetric nuclei such as $^4$He, $^{16}$O, $^{20}$Ne and $^{40}$Ca. Furthermore, we s…
▽ More
Using the \textit{ab initio} symmetry-adapted no-core shell model, we compute sum rules and response functions for light to medium-mass nuclei, starting from interactions that are derived in the chiral effective field theory. We investigate electromagnetic transitions of monopole, dipole and quadrupole nature for symmetric nuclei such as $^4$He, $^{16}$O, $^{20}$Ne and $^{40}$Ca. Furthermore, we study giant monopole resonance, which can provide information on the incompressibility of symmetric nuclear matter.
△ Less
Submitted 15 December, 2023;
originally announced December 2023.
-
Impact of two-body currents on magnetic dipole moments of nuclei
Authors:
T. Miyagi,
X. Cao,
R. Seutin,
S. Bacca,
R. F. Garcia Ruiz,
K. Hebeler,
J. D. Holt,
A. Schwenk
Abstract:
We investigate the effects of two-body currents on magnetic dipole moments of medium-mass and heavy nuclei using the valence-space in-medium similarity renormalization group with chiral effective field theory interactions and currents. Focusing on near doubly magic nuclei from oxygen to bismuth, we have found that the leading two-body currents globally improve the agreement with experimental magne…
▽ More
We investigate the effects of two-body currents on magnetic dipole moments of medium-mass and heavy nuclei using the valence-space in-medium similarity renormalization group with chiral effective field theory interactions and currents. Focusing on near doubly magic nuclei from oxygen to bismuth, we have found that the leading two-body currents globally improve the agreement with experimental magnetic moments. Moreover, our results show the importance of multi-shell effects for $^{41}$Ca, which suggest that the $Z=N=20$ gap in $^{40}$Ca is not as robust as in $^{48}$Ca. The increasing contribution of two-body currents in heavier systems is explained by the operator structure of the center-of-mass dependent Sachs term.
△ Less
Submitted 10 June, 2024; v1 submitted 24 November, 2023;
originally announced November 2023.
-
Magnetic dipole transition in $^{48}$Ca
Authors:
B. Acharya,
B. S. Hu,
S. Bacca,
G. Hagen,
P. Navrátil,
T. Papenbrock
Abstract:
The magnetic dipole transition strength $B(M1)$ of $^{48}$Ca is dominated by a single resonant state at an excitation energy of 10.23 MeV. Experiments disagree about $B(M1)$ and this impacts our understanding of spin flips in nuclei. We performed ab initio computations based on chiral effective field theory and found that $B(M1:0^+\rightarrow1^+)$ lies in the range from $7.0$ to $10.2~μ_N^2$. This…
▽ More
The magnetic dipole transition strength $B(M1)$ of $^{48}$Ca is dominated by a single resonant state at an excitation energy of 10.23 MeV. Experiments disagree about $B(M1)$ and this impacts our understanding of spin flips in nuclei. We performed ab initio computations based on chiral effective field theory and found that $B(M1:0^+\rightarrow1^+)$ lies in the range from $7.0$ to $10.2~μ_N^2$. This is consistent with a $(γ,n)$ experiment but larger than results from $(e,e^\prime)$ and $(p,p')$ scattering. Two-body currents yield no quenching of the $B(M1)$ strength and continuum effects reduce it by about 10%. For a validation of our approach, we computed magnetic moments in $^{47,49}$Ca and performed benchmark calculations in light nuclei.
△ Less
Submitted 7 June, 2024; v1 submitted 19 November, 2023;
originally announced November 2023.
-
$^{40}$Ca transverse response function from coupled-cluster theory
Authors:
Joanna Sobczyk,
Bijaya Acharya,
Sonia Bacca,
Gaute Hagen
Abstract:
We present calculations of the $^{40}$Ca transverse response function obtained from coupled-cluster theory used in conjunction with the Lorentz integral transform method. We employ nuclear forces derived at next-to-next-to leading order in chiral effective field theory with and without $Δ$ degrees of freedom. We first benchmark this approach on the $^4$He nucleus and compare both the transverse su…
▽ More
We present calculations of the $^{40}$Ca transverse response function obtained from coupled-cluster theory used in conjunction with the Lorentz integral transform method. We employ nuclear forces derived at next-to-next-to leading order in chiral effective field theory with and without $Δ$ degrees of freedom. We first benchmark this approach on the $^4$He nucleus and compare both the transverse sum rule and the response function to earlier calculations based on different methods. As expected from the power counting of the chiral expansion of electromagnetic currents and from previous studies, our results retaining only one-body term underestimate the experimental data for $^4$He by about $20\%$. However, when the method is applied to $^{40}$Ca at the same order of the expansion, response functions do not lack strength and agree well with the world electron scattering data. We discuss various sources of theoretical uncertainties and comment on the comparison of our results with the available experiments.
△ Less
Submitted 4 October, 2023;
originally announced October 2023.
-
$^{16}$O spectral function from coupled-cluster theory: applications to lepton-nucleus scattering
Authors:
Joanna E. Sobczyk,
Sonia Bacca
Abstract:
We calculate the $^{16}$O spectral function by combining coupled-cluster theory with a Gaussian integral transform and by expanding the integral kernel in terms of Chebyshev polynomials to allow for a quantification of the theoretical uncertainties. We perform an analysis of the spectral function and employ it to predict lepton-nucleus scattering. Our results well describe the $^{16}$O electron sc…
▽ More
We calculate the $^{16}$O spectral function by combining coupled-cluster theory with a Gaussian integral transform and by expanding the integral kernel in terms of Chebyshev polynomials to allow for a quantification of the theoretical uncertainties. We perform an analysis of the spectral function and employ it to predict lepton-nucleus scattering. Our results well describe the $^{16}$O electron scattering data in the quasi-elastic peak for momentum transfers $|\mathbf{q}|\gtrapprox500$ MeV and electron energies up to 1.2 GeV, extending therefore the so-called first principles approach to lepton-nucleus cross sections well into the relativistic regime. To prove the applicability of this method to neutrino-nucleus cross sections, we implement our $^{16}$O spectral functions in the NuWro Monte Carlo event generator and provide a comparison with recently published T2K neutrino data.
△ Less
Submitted 1 September, 2023;
originally announced September 2023.
-
Magnetic dipole operator from chiral effective field theory for many-body expansion methods
Authors:
R. Seutin,
O. J. Hernandez,
T. Miyagi,
S. Bacca,
K. Hebeler,
S. König,
A. Schwenk
Abstract:
Many-body approaches for atomic nuclei generally rely on a basis expansion of the nuclear states, interactions, and current operators. In this work, we derive the representation of the magnetic dipole operator in plane-wave and harmonic-oscillator basis states, as needed for Faddeev calculations of few-body systems or many-body calculations within, e.g., the no-core shell model, the in-medium reno…
▽ More
Many-body approaches for atomic nuclei generally rely on a basis expansion of the nuclear states, interactions, and current operators. In this work, we derive the representation of the magnetic dipole operator in plane-wave and harmonic-oscillator basis states, as needed for Faddeev calculations of few-body systems or many-body calculations within, e.g., the no-core shell model, the in-medium renormalization group, coupled-cluster theory, or the nuclear shell model. We focus in particular on the next-to-leading-order two-body contributions derived from chiral effective field theory. We provide detailed benchmarks and also comparisons with quantum Monte Carlo results for three-body systems. The derived operator matrix elements represent the basic input for studying magnetic properties of atomic nuclei based on chiral effective field theory.
△ Less
Submitted 4 December, 2023; v1 submitted 31 July, 2023;
originally announced August 2023.
-
Electric dipole polarizability of $^{40}$Ca
Authors:
R. W. Fearick,
P. von Neumann-Cosel,
S. Bacca,
J. Birkhan,
F. Bonaiti,
I. Brandherm,
G. Hagen,
H. Matsubara,
W. Nazarewicz,
N. Pietralla,
V. Yu. Ponomarev,
P. -G. Reinhard,
X. Roca-Maza,
A. Richter,
A. Schwenk,
J. Simonis,
A. Tamii
Abstract:
The electric dipole strength distribution in $^{40}$Ca between 5 and 25 MeV has been determined at RCNP, Osaka, from proton inelastic scattering experiments at very forward angles. Combined with total photoabsorption data at higher excitation energy, this enables an extraction of the electric dipole polarizability $α_\mathrm{D}$($^{40}$Ca) = 1.92(17) fm$^3$. Together with the measured $α_{\rm D}$…
▽ More
The electric dipole strength distribution in $^{40}$Ca between 5 and 25 MeV has been determined at RCNP, Osaka, from proton inelastic scattering experiments at very forward angles. Combined with total photoabsorption data at higher excitation energy, this enables an extraction of the electric dipole polarizability $α_\mathrm{D}$($^{40}$Ca) = 1.92(17) fm$^3$. Together with the measured $α_{\rm D}$ in $^{48}$Ca, it provides a stringent test of modern theoretical approaches, including coupled cluster calculations with chiral effective field theory interactions and state-of-the art energy density functionals. The emerging picture is that for this medium-mass region dipole polarizabilities are well described theoretically, with important constraints for the neutron skin in $^{48}$Ca and related equation of state quantities.
△ Less
Submitted 18 April, 2023; v1 submitted 15 February, 2023;
originally announced February 2023.
-
Comprehensive theory of the Lamb shift in light muonic atoms
Authors:
K. Pachucki,
V. Lensky,
F. Hagelstein,
S. S. Li Muli,
S. Bacca,
R. Pohl
Abstract:
A comprehensive theory of the Lamb shift in light muonic atoms such as $μ$H, $μ$D, $μ^3$He$^+$, and $μ^4$He$^+$ is presented, with all quantum electrodynamic corrections included at the precision level constrained by the uncertainty of nuclear structure effects. This analysis can be used in the global adjustment of fundamental constants and in the determination of nuclear charge radii. Further imp…
▽ More
A comprehensive theory of the Lamb shift in light muonic atoms such as $μ$H, $μ$D, $μ^3$He$^+$, and $μ^4$He$^+$ is presented, with all quantum electrodynamic corrections included at the precision level constrained by the uncertainty of nuclear structure effects. This analysis can be used in the global adjustment of fundamental constants and in the determination of nuclear charge radii. Further improvements in the understanding of electromagnetic interactions of light nuclei will allow for a promising test of fundamental interactions by comparison with ``normal" atomic spectroscopy, in particular, with H-D and $^3$He-$^4$He isotope shifts.
△ Less
Submitted 21 February, 2024; v1 submitted 28 December, 2022;
originally announced December 2022.
-
Muonic-Atom Spectroscopy and Impact on Nuclear Structure and Precision QED Theory
Authors:
Aldo Antognini,
Sonia Bacca,
Andreas Fleischmann,
Loredana Gastaldo,
Franziska Hagelstein,
Paul Indelicato,
Andreas Knecht,
Vadim Lensky,
Ben Ohayon,
Vladimir Pascalutsa,
Nancy Paul,
Randolf Pohl,
Frederik Wauters
Abstract:
Recent progress in laser and x-ray spectroscopy of muonic atoms offers promising long-term possibilities at the intersection of atomic, nuclear and particle physics. In muonic hydrogen, laser spectroscopy measurements will determine the ground-state hyperfine splitting (HFS) and additionally improve the Lamb shift by a factor of 5. Precision spectroscopy with cryogenic microcalorimeters has the po…
▽ More
Recent progress in laser and x-ray spectroscopy of muonic atoms offers promising long-term possibilities at the intersection of atomic, nuclear and particle physics. In muonic hydrogen, laser spectroscopy measurements will determine the ground-state hyperfine splitting (HFS) and additionally improve the Lamb shift by a factor of 5. Precision spectroscopy with cryogenic microcalorimeters has the potential to significantly improve the charge radii of the light nuclei in the $Z=3-8$ range. Complementary progress in precision should be achieved on the theory of nucleon- and nuclear-structure effects. The impact of this muonic-atom spectroscopy program will be amplified by the upcoming results from H and He$^+$ spectroscopy, simple molecules such as HD$^+$ and Penning trap measurements. In this broader context, one can test ab-initio nuclear theories, bound-state QED for two- or three-body systems, and determine fundamental constants, such as the Rydberg ($R_\infty$) and the fine-structure ($α$) constants.
△ Less
Submitted 30 October, 2022;
originally announced October 2022.
-
$β^{+}$ Gamow-Teller strengths from unstable $^{14}$O via the $(d,{}^2\text{He})$ reaction in inverse kinematics
Authors:
S. Giraud,
J. C. Zamora,
R. Zegers,
D. Bazin,
Y. Ayyad,
S. Bacca,
S. Beceiro-Novo,
B. A. Brown,
A. Carls,
J. Chen,
M. Cortesi,
M. DeNudt,
G. Hagen,
C. Hultquist,
C. Maher,
W. Mittig,
F. Ndayisabye,
S. Noji,
S. J. Novario,
J. Pereira,
Z. Rahman,
J. Schmitt,
M. Serikow,
L. J. Sun,
J. Surbrook
, et al. (2 additional authors not shown)
Abstract:
For the first time, the $(d,{}^2\text{He})$ reaction was successfully used in inverse kinematics to extract the Gamow-Teller transition strength in the $β^{+}$ direction from an unstable nucleus. The nucleus studied was $^{14}$O, and the Gamow-Teller transition strength to $^{14}$N was extracted up to an excitation energy of 22 MeV. The measurement of the $(d,{}^2\text{He})$ reaction in inverse ki…
▽ More
For the first time, the $(d,{}^2\text{He})$ reaction was successfully used in inverse kinematics to extract the Gamow-Teller transition strength in the $β^{+}$ direction from an unstable nucleus. The nucleus studied was $^{14}$O, and the Gamow-Teller transition strength to $^{14}$N was extracted up to an excitation energy of 22 MeV. The measurement of the $(d,{}^2\text{He})$ reaction in inverse kinematics was made possible by the combination of an active target time projection chamber and a magnetic spectrometer. The data were used to test shell-model and state-of-the-art coupled cluster calculations. Shell-model calculations reproduce the measured Gamow-Teller strength distribution up to about 15 MeV reasonably well, after the application of a phenomenological quenching factor. Coupled-cluster calculation reproduces the full strength distribution well without such quenching, owing to the large model space, the inclusion of strong correlations, and the coupling of the weak interaction to two nucleons through two-body currents. This indicates that such calculations provide a very promising path for answering long-standing questions about the observed quenching of Gamow-Teller strengths in nuclei.
△ Less
Submitted 27 October, 2022;
originally announced October 2022.
-
Uncertainty quantification in electromagnetic observables of nuclei
Authors:
Bijaya Acharya,
Sonia Bacca,
Francesca Bonaiti,
Simone Salvatore Li Muli,
Joanna E. Sobczyk
Abstract:
We present strategies to quantify theoretical uncertainties in modern ab-initio calculations of electromagnetic observables in light and medium-mass nuclei. We discuss how uncertainties build up from various sources, such as the approximations introduced by the few- or many-body solver and the truncation of the chiral effective field theory expansion. We review the recent progress encompassing a b…
▽ More
We present strategies to quantify theoretical uncertainties in modern ab-initio calculations of electromagnetic observables in light and medium-mass nuclei. We discuss how uncertainties build up from various sources, such as the approximations introduced by the few- or many-body solver and the truncation of the chiral effective field theory expansion. We review the recent progress encompassing a broad range of electromagnetic observables in stable and unstable nuclei.
△ Less
Submitted 10 October, 2022;
originally announced October 2022.
-
Spectral function for $^4$He using the Chebyshev expansion in coupled-cluster theory
Authors:
J. E. Sobczyk,
S. Bacca,
G. Hagen,
T. Papenbrock
Abstract:
We compute spectral function for $^4$He by combining coupled-cluster theory with an expansion of integral transforms into Chebyshev polynomials. Our method allows to estimate the uncertainty of spectral reconstruction. The properties of the Chebyshev polynomials make the procedure numerically stable and considerably lower in memory usage than the typically employed Lanczos algorithm. We benchmark…
▽ More
We compute spectral function for $^4$He by combining coupled-cluster theory with an expansion of integral transforms into Chebyshev polynomials. Our method allows to estimate the uncertainty of spectral reconstruction. The properties of the Chebyshev polynomials make the procedure numerically stable and considerably lower in memory usage than the typically employed Lanczos algorithm. We benchmark our predictions with other calculations in the literature and with electron scattering data in the quasi-elastic peak. The spectral function formalism allows one to extend ab-initio lepton-nucleus cross sections into the relativistic regime. This makes it a promising tool for modeling this process at higher energy transfers. The results we present open the door for studies of heavier nuclei, important for the neutrino oscillation programs.
△ Less
Submitted 7 May, 2022;
originally announced May 2022.
-
Bayesian analysis of nuclear polarizability corrections to the Lamb shift of muonic H-atoms and He-ions
Authors:
Simone Salvatore Li Muli,
Bijaya Acharya,
Oscar Javier Hernandez,
Sonia Bacca
Abstract:
The extraction of nuclear charge radii from spectroscopy experiments in muonic atoms is currently limited by the large uncertainties associated with the theoretical evaluation of the nuclear polarizability effects. To facilitate calculations, these polarizability corrections are conventionally expressed as an expansion in a dimensionless parameter $η$, which has been argued in previous literature…
▽ More
The extraction of nuclear charge radii from spectroscopy experiments in muonic atoms is currently limited by the large uncertainties associated with the theoretical evaluation of the nuclear polarizability effects. To facilitate calculations, these polarizability corrections are conventionally expressed as an expansion in a dimensionless parameter $η$, which has been argued in previous literature to hold an approximate value of 0.33 in light-nuclear systems. In this work, we check this claim by doing a Bayesian analysis of the nuclear-polarizability corrections to the Lamb shift in $μ^2$H and $μ^3$H atoms and in $μ^3$He$^+$ and $μ^4$He$^+$ ions at various orders in the $η$-expansion. Our analysis supports the claim that $η\ll 1$ in these systems and finds truncation uncertainties that are similar to the previous estimate, the only exception being the truncation uncertainties in the $μ^3$He$^+$ ion, which are found to be larger.
△ Less
Submitted 7 October, 2022; v1 submitted 21 March, 2022;
originally announced March 2022.
-
Theoretical tools for neutrino scattering: interplay between lattice QCD, EFTs, nuclear physics, phenomenology, and neutrino event generators
Authors:
L. Alvarez Ruso,
A. M. Ankowski,
S. Bacca,
A. B. Balantekin,
J. Carlson,
S. Gardiner,
R. Gonzalez-Jimenez,
R. Gupta,
T. J. Hobbs,
M. Hoferichter,
J. Isaacson,
N. Jachowicz,
W. I. Jay,
T. Katori,
F. Kling,
A. S. Kronfeld,
S. W. Li,
H. -W. Lin,
K. -F. Liu,
A. Lovato,
K. Mahn,
J. Menendez,
A. S. Meyer,
J. Morfin,
S. Pastore
, et al. (36 additional authors not shown)
Abstract:
Maximizing the discovery potential of increasingly precise neutrino experiments will require an improved theoretical understanding of neutrino-nucleus cross sections over a wide range of energies. Low-energy interactions are needed to reconstruct the energies of astrophysical neutrinos from supernovae bursts and search for new physics using increasingly precise measurement of coherent elastic neut…
▽ More
Maximizing the discovery potential of increasingly precise neutrino experiments will require an improved theoretical understanding of neutrino-nucleus cross sections over a wide range of energies. Low-energy interactions are needed to reconstruct the energies of astrophysical neutrinos from supernovae bursts and search for new physics using increasingly precise measurement of coherent elastic neutrino scattering. Higher-energy interactions involve a variety of reaction mechanisms including quasi-elastic scattering, resonance production, and deep inelastic scattering that must all be included to reliably predict cross sections for energies relevant to DUNE and other accelerator neutrino experiments. This white paper discusses the theoretical status, challenges, required resources, and path forward for achieving precise predictions of neutrino-nucleus scattering and emphasizes the need for a coordinated theoretical effort involved lattice QCD, nuclear effective theories, phenomenological models of the transition region, and event generators.
△ Less
Submitted 20 April, 2022; v1 submitted 16 March, 2022;
originally announced March 2022.
-
Electron Scattering and Neutrino Physics
Authors:
A. M. Ankowski,
A. Ashkenazi,
S. Bacca,
J. L. Barrow,
M. Betancourt,
A. Bodek,
M. E. Christy,
L. Doria. S. Dytman,
A. Friedland,
O. Hen,
C. J. Horowitz,
N. Jachowicz,
W. Ketchum,
T. Lux,
K. Mahn,
C. Mariani,
J. Newby,
V. Pandey,
A. Papadopoulou,
E. Radicioni,
F. Sánchez,
C. Sfienti,
J. M. Udías,
L. Weinstein,
L. Alvarez-Ruso
, et al. (28 additional authors not shown)
Abstract:
A thorough understanding of neutrino-nucleus scattering physics is crucial for the successful execution of the entire US neutrino physics program. Neutrino-nucleus interaction constitutes one of the biggest systematic uncertainties in neutrino experiments - both at intermediate energies affecting long-baseline Deep Underground Neutrino Experiment (DUNE), as well as at low energies affecting cohere…
▽ More
A thorough understanding of neutrino-nucleus scattering physics is crucial for the successful execution of the entire US neutrino physics program. Neutrino-nucleus interaction constitutes one of the biggest systematic uncertainties in neutrino experiments - both at intermediate energies affecting long-baseline Deep Underground Neutrino Experiment (DUNE), as well as at low energies affecting coherent scattering neutrino program - and could well be the difference between achieving or missing discovery level precision. To this end, electron-nucleus scattering experiments provide vital information to test, assess and validate different nuclear models and event generators intended to be used in neutrino experiments. In this white paper, we highlight connections between electron- and neutrino-nucleus scattering physics at energies ranging from 10s of MeV to a few GeV, review the status of ongoing and planned electron scattering experiments, identify gaps, and layout a path forward that benefits the neutrino community. We also highlight the systemic challenges with respect to the divide between the nuclear and high-energy physics communities and funding that presents additional hurdle in mobilizing these connections to the benefit of neutrino programs.
△ Less
Submitted 10 May, 2023; v1 submitted 14 March, 2022;
originally announced March 2022.
-
Measurement of the α-particle monopole transition form factor challenges theory: a low-energy puzzle for nuclear forces?
Authors:
S. Kegel,
P. Achenbach,
S. Bacca,
N. Barnea,
J. Bericic,
D. Bosnar,
L. Correa,
M. O. Distler,
A. Esser,
H. Fonvieille,
I. Friscic,
M. Heilig,
P. Herrmann,
M. Hoek,
P. Klag,
T. Kolar,
W. Leidemann,
H. Merkel,
M. Mihovilovic,
J. Mueller,
U. Mueller,
G. Orlandini,
J. Pochodzalla,
B. S. Schlimme,
M. Schoth
, et al. (8 additional authors not shown)
Abstract:
We perform a systematic study of the $α$-particle excitation from its ground state $0_1^+$ to the $0_2^+$ resonance. The so-called monopole transition form factor is investigated via an electron scattering experiment in a broad $Q^2$-range (from $0.5$ to $5.0$ fm$^{-2}$). The precision of the new data dramatically superseeds that of older sets of data, each covering only a portion of the $Q^2$-ran…
▽ More
We perform a systematic study of the $α$-particle excitation from its ground state $0_1^+$ to the $0_2^+$ resonance. The so-called monopole transition form factor is investigated via an electron scattering experiment in a broad $Q^2$-range (from $0.5$ to $5.0$ fm$^{-2}$). The precision of the new data dramatically superseeds that of older sets of data, each covering only a portion of the $Q^2$-range. The new data allow the determination of two coefficients in a low-momentum expansion leading to a new puzzle. By confronting experiment to state-of-the-art theoretical calculations we observe that modern nuclear forces, including those derived within chiral effective field theory which are well tested on a variety of observables, fail to reproduce the excitation of the $α$-particle.
△ Less
Submitted 20 December, 2021;
originally announced December 2021.
-
Ab-initio coupled-cluster calculations of ground and dipole excited states in 8He
Authors:
Francesca Bonaiti,
Sonia Bacca,
Gaute Hagen
Abstract:
We perform coupled-cluster calculations of ground- and dipole excited-state properties of the 8He halo nucleus with nucleon-nucleon and three-nucleon interactions from chiral effective field theory, both with and without explicit delta degrees of freedom. By increasing the precision in our coupled-cluster calculations via the inclusion of leading order three-particle three-hole excitations in the…
▽ More
We perform coupled-cluster calculations of ground- and dipole excited-state properties of the 8He halo nucleus with nucleon-nucleon and three-nucleon interactions from chiral effective field theory, both with and without explicit delta degrees of freedom. By increasing the precision in our coupled-cluster calculations via the inclusion of leading order three-particle three-hole excitations in the cluster operator, we obtain a ground-state energy and a charge radius that are consistent with experiment, albeit with a slight under-binding. We also investigate the excited states induced by the electric dipole operator and present a discussion on the Thomas-Reiche-Kuhn and cluster sum rules. Finally, we compute the electric dipole polarizability, providing a theoretical benchmark for future experimental determinations that will study this exotic nucleus.
△ Less
Submitted 14 March, 2022; v1 submitted 15 December, 2021;
originally announced December 2021.
-
Nuclear Charge Radii of the Nickel Isotopes $^{58-68,70}$Ni
Authors:
S. Malbrunot-Ettenauer,
S. Kaufmann,
S. Bacca,
C. Barbieri,
J. Billowes,
M. L. Bissell,
K. Blaum,
B. Cheal,
T. Duguet,
R. F. Garcia Ruiz,
W. Gins,
C. Gorges,
G. Hagen,
H. Heylen,
J. D. Holt,
G. R. Jansen,
A. Kanellakopoulos,
M. Kortelainen,
T. Miyagi,
P. Navrátil,
W. Nazarewicz,
R. Neugart,
G. Neyens,
W. Nörtershäuser,
S. J. Novario
, et al. (16 additional authors not shown)
Abstract:
Collinear laser spectroscopy is performed on the nickel isotopes $^{58-68,70}$Ni, using a time-resolved photon counting system. From the measured isotope shifts, nuclear charge radii $R_c$ are extracted and compared to theoretical results. Three ab initio approaches all employ, among others, the chiral interaction NNLO$_{\rm sat}$, which allows an assessment of their accuracy. We find agreement wi…
▽ More
Collinear laser spectroscopy is performed on the nickel isotopes $^{58-68,70}$Ni, using a time-resolved photon counting system. From the measured isotope shifts, nuclear charge radii $R_c$ are extracted and compared to theoretical results. Three ab initio approaches all employ, among others, the chiral interaction NNLO$_{\rm sat}$, which allows an assessment of their accuracy. We find agreement with experiment in differential radii $δ\left\langle r_\mathrm{c}^2 \right\rangle$ for all employed ab initio methods and interactions, while the absolute radii are consistent with data only for NNLO$_{\rm sat}$. Within nuclear density functional theory, the Skyrme functional SV-min matches experiment more closely than the Fayans functional Fy($Δr$,HFB).
△ Less
Submitted 6 December, 2021;
originally announced December 2021.
-
Removing center of mass effects in response function and sum rule calculations based on the harmonic oscillator basis
Authors:
K. Akdogan,
D. Layh,
I. Weinberger,
J. Simonis,
N. Barnea,
S. Bacca
Abstract:
Response functions are at the heart of any comparison of theory with experiment in studies of the nuclear dynamics with electroweak probes. Calculations performed in the laboratory frame often suffer from center of mass contaminations that need to be removed. By confining the system in a harmonic oscillator, we derive a set of analytical formulas to subtract the center of mass effects from calcula…
▽ More
Response functions are at the heart of any comparison of theory with experiment in studies of the nuclear dynamics with electroweak probes. Calculations performed in the laboratory frame often suffer from center of mass contaminations that need to be removed. By confining the system in a harmonic oscillator, we derive a set of analytical formulas to subtract the center of mass effects from calculations of response functions and associated sum rules. After a general analytical derivation, we first deal specifically with the longitudinal response function appearing in electron scattering and provide expressions for the center of mass correcting functions. Next, we present a proof of principle study for the case of the electric dipole sum rules in a two-body problem with a numerical implementation of our formalism. These steps pave the way to applying the proposed method to heavier nuclei in the future.
△ Less
Submitted 29 October, 2021;
originally announced October 2021.
-
Gaussian process error modeling for chiral effective-field-theory calculations of $np\leftrightarrow dγ$ at low energies
Authors:
Bijaya Acharya,
Sonia Bacca
Abstract:
We calculate the energy-dependent cross section of the $np\leftrightarrow dγ$ process in chiral effective field theory and apply state-of-the-art tools for quantification of theory uncertainty. We focus on the low-energy regime, where the magnetic dipole and the electric dipole transitions cross over, including the range relevant for big-bang nucleosynthesis. Working with the leading one- and two-…
▽ More
We calculate the energy-dependent cross section of the $np\leftrightarrow dγ$ process in chiral effective field theory and apply state-of-the-art tools for quantification of theory uncertainty. We focus on the low-energy regime, where the magnetic dipole and the electric dipole transitions cross over, including the range relevant for big-bang nucleosynthesis. Working with the leading one- and two-body electromagnetic currents, we study the order-by-order convergence of this observable in the chiral expansion of the nuclear potential. We find that the Gaussian process error model describes the observed convergence very well, allowing us to present Bayesian credible intervals for the truncation error with correlations between the cross sections at different energies taken into account. We obtain a 1$σ$ estimate of about 0.2\% for the uncertainty from the truncation of the nuclear potential. This is an important step towards calculations with statistically interpretable uncertainties for astrophysical reactions involving light nuclei.
△ Less
Submitted 23 March, 2022; v1 submitted 28 September, 2021;
originally announced September 2021.
-
Ab initio computation of the longitudinal response function in $^{40}$Ca
Authors:
J. E. Sobczyk,
B. Acharya,
S. Bacca,
G. Hagen
Abstract:
We present a consistent \emph{ab initio} computation of the longitudinal response function $R_L$ in $^{40}$Ca using the coupled-cluster and Lorentz integral transform methods starting from chiral nucleon-nucleon and three-nucleon interactions. We validate our approach by comparing our results for $R_L$ in $^4$He and the Coulomb sum rule in $^{40}$Ca against experimental data and other calculations…
▽ More
We present a consistent \emph{ab initio} computation of the longitudinal response function $R_L$ in $^{40}$Ca using the coupled-cluster and Lorentz integral transform methods starting from chiral nucleon-nucleon and three-nucleon interactions. We validate our approach by comparing our results for $R_L$ in $^4$He and the Coulomb sum rule in $^{40}$Ca against experimental data and other calculations. For $R_L$ in $^{40}$Ca we obtain a very good agreement with experiment in the quasi-elastic peak up to intermediate momentum transfers, and we find that final state interactions are essential for an accurate description of the data. This work presents a milestone towards \emph{ab initio} computations of neutrino-nucleus cross sections relevant for experimental long-baseline neutrino programs.
△ Less
Submitted 12 March, 2021; v1 submitted 11 March, 2021;
originally announced March 2021.
-
Implementation of local chiral interactions in the hyperspherical harmonics formalism
Authors:
Simone Salvatore Li Muli,
Sonia Bacca,
Nir Barnea
Abstract:
With the goal of using chiral interactions at various orders to explore properties of the few-body nuclear systems, we write the recently developed local chiral interactions as spherical irreducible tensors and implement them in the hyperspherical harmonics expansion method. We devote particular attention to three-body forces at next-to-next-to leading order, which play an important role in reprod…
▽ More
With the goal of using chiral interactions at various orders to explore properties of the few-body nuclear systems, we write the recently developed local chiral interactions as spherical irreducible tensors and implement them in the hyperspherical harmonics expansion method. We devote particular attention to three-body forces at next-to-next-to leading order, which play an important role in reproducing experimental data. We check our implementation by benchmarking the ground-state properties of $^3$H, $^3$He and $^4$He against the available Monte Carlo calculations. We then confirm their order-by-order truncation error estimates and further investigate uncertainties in the charge radii obtained by using the precise muonic atom data for single-nucleon radii. Having local chiral Hamiltonians at various orders implemented in our hyperspherical harmonics suites of codes opens up the possibility to test such interactions on other light-nuclei properties, such as electromagnetic reactions.
△ Less
Submitted 22 February, 2021;
originally announced February 2021.
-
Dispersive evaluation of the Lamb shift in muonic deuterium from chiral effective field theory
Authors:
Bijaya Acharya,
Vadim Lensky,
Sonia Bacca,
Mikhail Gorchtein,
Marc Vanderhaeghen
Abstract:
We merge the dispersive relation approach and the ab initio method to compute nuclear structure corrections to the Lamb shift in muonic deuterium. We calculate the deuteron response functions and corresponding uncertainties up to next-to-next-to-next-to-leading order in chiral effective field theory and compare our results to selected electromagnetic data to test the validity of the theory. We the…
▽ More
We merge the dispersive relation approach and the ab initio method to compute nuclear structure corrections to the Lamb shift in muonic deuterium. We calculate the deuteron response functions and corresponding uncertainties up to next-to-next-to-next-to-leading order in chiral effective field theory and compare our results to selected electromagnetic data to test the validity of the theory. We then feed response functions calculated over a wide range of kinematics to the dispersion-theory formalism and show that an improved accuracy is obtained compared to that with the use of available experimental data in the dispersive analysis. This opens up the possibility of applying this hybrid method to other light muonic atoms and supplementing experimental data with ab initio theory for kinematics where data are scarce or difficult to measure with the goal of reducing uncertainties in estimates of nuclear structure effects in atomic spectroscopy.
△ Less
Submitted 1 February, 2021; v1 submitted 21 October, 2020;
originally announced October 2020.
-
Coulomb sum rule for $^4$He and $^{16}$O from coupled-cluster theory
Authors:
J. E. Sobczyk,
B. Acharya,
S. Bacca,
G. Hagen
Abstract:
We demonstrate the capability of coupled-cluster theory to compute the Coulomb sum rule for the $^4$He and $^{16}$O nuclei using interactions from chiral effective field theory. We perform several checks, including a few-body benchmark for $^4$He. We provide an analysis of the center-of-mass contaminations, which we are able to safely remove. We then compare with other theoretical results and expe…
▽ More
We demonstrate the capability of coupled-cluster theory to compute the Coulomb sum rule for the $^4$He and $^{16}$O nuclei using interactions from chiral effective field theory. We perform several checks, including a few-body benchmark for $^4$He. We provide an analysis of the center-of-mass contaminations, which we are able to safely remove. We then compare with other theoretical results and experimental data available in the literature, obtaining a fair agreement. This is a first and necessary step towards initiating a program for computing neutrino-nucleus interactions from first principles and supporting the experimental long-baseline neutrino program with a state-of-the-art theory that can reach medium-mass nuclei.
△ Less
Submitted 3 September, 2020;
originally announced September 2020.
-
Role of chiral two-body currents in $^6$Li magnetic properties in light of a new precision measurement with the relative self-absorption technique
Authors:
U. Friman-Gayer,
C. Romig,
T. Hüther,
K. Albe,
S. Bacca,
T. Beck,
M. Berger,
J. Birkhan,
K. Hebeler,
O. J. Hernandez,
J. Isaak,
S. König,
N. Pietralla,
P. C. Ries,
J. Rohrer,
R. Roth,
D. Savran,
M. Scheck,
A. Schwenk,
R. Seutin,
V. Werner
Abstract:
A direct measurement of the decay width of the excited $0^+_1$ state of $^6$Li using the relative self-absorption technique is reported. Our value of $Γ_{γ, 0^+_1 \to 1^+_1} = 8.17(14)_\mathrm{stat.}(11)_\mathrm{syst.} \mathrm{eV}$ provides sufficiently low experimental uncertainties to test modern theories of nuclear forces. The corresponding transition rate is compared to the results of ab initi…
▽ More
A direct measurement of the decay width of the excited $0^+_1$ state of $^6$Li using the relative self-absorption technique is reported. Our value of $Γ_{γ, 0^+_1 \to 1^+_1} = 8.17(14)_\mathrm{stat.}(11)_\mathrm{syst.} \mathrm{eV}$ provides sufficiently low experimental uncertainties to test modern theories of nuclear forces. The corresponding transition rate is compared to the results of ab initio calculations based on chiral effective field theory that take into account contributions to the magnetic dipole operator beyond leading order. This enables a precision test of the impact of two-body currents that enter at next-to-leading order.
△ Less
Submitted 20 January, 2021; v1 submitted 15 May, 2020;
originally announced May 2020.
-
Charge radius of the short-lived $^{68}$Ni and correlation with the dipole polarizability
Authors:
S. Kaufmann,
J. Simonis,
S. Bacca,
J. Billowes,
M. L. Bissell,
K. Blaum,
B. Cheal,
R. F. Garcia Ruiz,
W. Gins,
C. Gorges,
G. Hagen,
H. Heylen,
A. Kanellakopoulos,
S. Malbrunot-Ettenauer,
M. Miorelli,
R. Neugart,
G. Neyens,
W. Nörtershäuser,
R. Sánchez,
S. Sailer,
A. Schwenk,
T. Ratajczyk,
L. V. Rodríguez,
L. Wehner,
C. Wraith
, et al. (4 additional authors not shown)
Abstract:
We present the first laser spectroscopic measurement of the neutron-rich nucleus $^{68}$Ni at the \mbox{$N=40$} subshell closure and extract its nuclear charge radius. Since this is the only short-lived isotope for which the dipole polarizability $α_{\rm D}$ has been measured, the combination of these observables provides a benchmark for nuclear structure theory. We compare them to novel coupled-c…
▽ More
We present the first laser spectroscopic measurement of the neutron-rich nucleus $^{68}$Ni at the \mbox{$N=40$} subshell closure and extract its nuclear charge radius. Since this is the only short-lived isotope for which the dipole polarizability $α_{\rm D}$ has been measured, the combination of these observables provides a benchmark for nuclear structure theory. We compare them to novel coupled-cluster calculations based on different chiral two- and three-nucleon interactions, for which a strong correlation between the charge radius and dipole polarizability is observed, similar to the stable nucleus $^{48}$Ca. Three-particle--three-hole correlations in coupled-cluster theory substantially improve the description of the experimental data, which allows to constrain the neutron radius and neutron skin of $^{68}$Ni.
△ Less
Submitted 13 March, 2020;
originally announced March 2020.
-
Benchmark calculations of electromagnetic sum rules with a symmetry-adapted basis and hyperspherical harmonics
Authors:
R. B. Baker,
K. D. Launey,
S. Bacca,
N. Nevo Dinur,
T. Dytrych
Abstract:
We demonstrate the ability to calculate electromagnetic sum rules with the \textit{ab initio} symmetry-adapted no-core shell model. By implementing the Lanczos algorithm, we compute non-energy weighted, energy weighted, and inverse energy weighted sum rules for electric monopole, dipole, and quadrupole transitions in $^4$He using realistic interactions. We benchmark the results with the hyperspher…
▽ More
We demonstrate the ability to calculate electromagnetic sum rules with the \textit{ab initio} symmetry-adapted no-core shell model. By implementing the Lanczos algorithm, we compute non-energy weighted, energy weighted, and inverse energy weighted sum rules for electric monopole, dipole, and quadrupole transitions in $^4$He using realistic interactions. We benchmark the results with the hyperspherical harmonics method and show agreement within $2σ$, where the uncertainties are estimated from the use of the many-body technique. We investigate the dependence of the results on three different interactions, including chiral potentials, and we report on the $^4$He electric dipole polarizability calculated in the SA-NCSM that reproduces the experimental data and earlier theoretical outcomes. We also detail a novel use of the Lawson procedure to remove the spurious center-of-mass contribution to the sum rules that arises from using laboratory-frame coordinates. We further show that this same technique can be applied in the Lorentz integral transform method, with a view toward studies of electromagnetic reactions for light through medium-mass nuclei.
△ Less
Submitted 12 March, 2020;
originally announced March 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.
-
Neutrino-deuteron scattering: Uncertainty quantification and new $L_{1,A}$ constraints
Authors:
Bijaya Acharya,
Sonia Bacca
Abstract:
We study neutral- and charged-current (anti)neutrino-induced dissociation of the deuteron at energies from threshold up to 150 MeV by employing potentials, as well as one- and two-body currents, derived in chiral effective field theory ($χ$EFT). We provide uncertainty estimates from $χ$EFT truncations of the electroweak current, dependences on the $χ$EFT cutoff and variations in the pool of fit da…
▽ More
We study neutral- and charged-current (anti)neutrino-induced dissociation of the deuteron at energies from threshold up to 150 MeV by employing potentials, as well as one- and two-body currents, derived in chiral effective field theory ($χ$EFT). We provide uncertainty estimates from $χ$EFT truncations of the electroweak current, dependences on the $χ$EFT cutoff and variations in the pool of fit data used to fix the low-energy constants of $χ$EFT. At 100 MeV of incident (anti)neutrino energy, these uncertainties amount to about 2-3\% and are smaller than the sensitivity of the cross sections to the single-nucleon axial form factor, which amounts to 5\% if one varies the range of the nucleon axial radius within the bands determined by recent lattice quantum chromodynamics evaluations and phenomenological extractions. We conclude that a precise determination of the nucleon axial form factor is required for a high-precision calculation of the neutrino-deuteron cross sections at energies higher than 100 MeV. By matching our low-energy $χ$EFT results to those of pionless effective field theory ($\cancelπ$EFT), we provide new constraints for the counterterm $L_{1,A}$ that parameterizes the strength of the axial two-body current in $\cancelπ$EFT. We obtain a value of $4.9^{+1.9}_{-1.5}\mathrm{fm}^3$ at renormalization scale set to pion mass, which is compatible with, albeit narrower than, previous experimental determinations, and comparable to a recent lattice quantum chromodynamics calculation.
△ Less
Submitted 29 January, 2020; v1 submitted 28 November, 2019;
originally announced November 2019.
-
Muonic Lithium atoms: nuclear structure corrections to the Lamb shift
Authors:
S. Li Muli,
A. Poggialini,
S. Bacca
Abstract:
In view of the future plans to measure the Lamb shift in muonic Lithium atoms we address the microscopic theory of the $μ$-$^6$Li$^{2+}$ and $μ$-$^7$Li$^{2+}$ systems. The goal of the CREMA collaboration is to measure the Lamb shift to extract the charge radius with high precision and compare it to electron scattering data or atomic spectroscopy to see if interesting puzzles, such as the proton an…
▽ More
In view of the future plans to measure the Lamb shift in muonic Lithium atoms we address the microscopic theory of the $μ$-$^6$Li$^{2+}$ and $μ$-$^7$Li$^{2+}$ systems. The goal of the CREMA collaboration is to measure the Lamb shift to extract the charge radius with high precision and compare it to electron scattering data or atomic spectroscopy to see if interesting puzzles, such as the proton and deuteron radius puzzles, arise. For this experiment to be successful, theoretical information on the nuclear structure corrections to the Lamb shift is needed. For $μ$-$^6$Li$^{2+}$ and $μ$-$^7$Li$^{2+}$ there exist only estimates of nuclear structure corrections based on experimental data that suffer from very large uncertainties. We present the first steps towards an ab initio computation of these quantities using few-body techniques.
△ Less
Submitted 15 December, 2019; v1 submitted 31 October, 2019;
originally announced October 2019.
-
Probing uncertainties of nuclear structure corrections in light muonic atoms
Authors:
Oscar Javier Hernandez,
Chen Ji,
Sonia Bacca,
Nir Barnea
Abstract:
Recent calculations of nuclear structure corrections to the Lamb shift in light muonic atoms are based on an expansion in a parameter eta, where only terms up to second order are retained. The parameter eta can be shown to be proportional to the square root of the muon/proton mass ratio, so that it is small and the expansion is expected to converge. However, practical implementations show that the…
▽ More
Recent calculations of nuclear structure corrections to the Lamb shift in light muonic atoms are based on an expansion in a parameter eta, where only terms up to second order are retained. The parameter eta can be shown to be proportional to the square root of the muon/proton mass ratio, so that it is small and the expansion is expected to converge. However, practical implementations show that the eta convergence may be slower than expected. In this work we probe the uncertainties due to this expansion using a different formalism, which is based on a multipole expansion of the longitudinal and transverse response functions and was first introduced by Leidemann and Rosenfelder. We refer to this alternative expansion as the eta-less formalism. We generalize this formalism to account for the cancellation of elastic terms such as the third Zemach moment (or Friar moment) and embed it in a computationally efficient framework. We implement and test this approach in the case of muonic deuterium. The comparison of results in the point nucleon limit for both methods achieve sub-percent agreement. When nucleon form factors are introduced we find a 4% and 2% difference in the third Zemach moment and nuclear polarizability, respectively, compared to the eta-less expansion, indicating that the nucleon form factor approximations should be improved. However, we find that the sum of these terms removes this dependence and the uncertainty due to the eta-expansion and the related second-order approximation in the nucleon form factors amounts only to 0.2% and thus is fully justified in muonic deuterium. This computationally efficient framework paves the way to further studies in light muonic systems with more than two nucleons, where controlling and reducing uncertainties in nuclear structure corrections is key to the experimental efforts of the CREMA collaboration.
△ Less
Submitted 12 September, 2019;
originally announced September 2019.
-
Coherent elastic neutrino-nucleus scattering on 40Ar from first principles
Authors:
C. G. Payne,
S. Bacca,
G. Hagen,
W. Jiang,
T. Papenbrock
Abstract:
Coherent elastic neutrino scattering on the 40Ar nucleus is computed with coupled-cluster theory based on nuclear Hamiltonians inspired by effective field theories of quantum chromodynamics. Our approach is validated by calculating the charge form factor and comparing it to data from electron scattering. We make predictions for the weak form factor, the neutron radius, and the neutron skin, and es…
▽ More
Coherent elastic neutrino scattering on the 40Ar nucleus is computed with coupled-cluster theory based on nuclear Hamiltonians inspired by effective field theories of quantum chromodynamics. Our approach is validated by calculating the charge form factor and comparing it to data from electron scattering. We make predictions for the weak form factor, the neutron radius, and the neutron skin, and estimate systematic uncertainties. The neutron-skin thickness of 40Ar40 is consistent with results from density functional theory. Precision measurements from coherent elastic neutrino-nucleus scattering could potentially be used to extract these observables and help to constrain nuclear models.
△ Less
Submitted 26 August, 2019;
originally announced August 2019.
-
First principles electromagnetic responses in medium-mass nuclei
Authors:
Johannes Simonis,
Sonia Bacca,
Gaute Hagen
Abstract:
We review the recent progress made in the computation of electromagnetic response functions in light and medium-mass nuclei using coupled-cluster theory. We show how a many-body formulation of the Lorentz integral transform method allows to calculate the photoabsorption cross sections of $^{16, 22}$O and $^{40}$Ca. Then, we discuss electromagnetic sum rules, with particular emphasis on the electri…
▽ More
We review the recent progress made in the computation of electromagnetic response functions in light and medium-mass nuclei using coupled-cluster theory. We show how a many-body formulation of the Lorentz integral transform method allows to calculate the photoabsorption cross sections of $^{16, 22}$O and $^{40}$Ca. Then, we discuss electromagnetic sum rules, with particular emphasis on the electric dipole polarizability, $α_D$. By including triples corrections in coupled-cluster theory, we revisit $^{48}$Ca, for which, beside the electric dipole polarizability, we had previously investigated the neutron and proton radii, as well as the size of the neutron-skin thickness [1]. We show that correlations among these observables still hold, albeit a better agreement with experiment is obtained for $α_D$ and the prediction of a small neutron-skin thickness is further corroborated.
△ Less
Submitted 6 May, 2019;
originally announced May 2019.
-
A statistical analysis of the nuclear structure uncertainties in $μ$D
Authors:
Oscar J. Hernandez,
Sonia Bacca,
Nir Barnea,
Nir Nevo-Dinur,
Andreas Ekström,
Chen Ji
Abstract:
The charge radius of the deuteron (D), was recently determined to three times the precision compared with previous measurements using the measured Lamb shift in muonic deuterium (muD). However, the muD value is 5.6 $σ$ smaller than the world averaged CODATA-2014 value [1]. To shed light on this discrepancy we analyze the uncertainties of the nuclear structure calculations of the Lamb shift in muD…
▽ More
The charge radius of the deuteron (D), was recently determined to three times the precision compared with previous measurements using the measured Lamb shift in muonic deuterium (muD). However, the muD value is 5.6 $σ$ smaller than the world averaged CODATA-2014 value [1]. To shed light on this discrepancy we analyze the uncertainties of the nuclear structure calculations of the Lamb shift in muD and conclude that nuclear theory uncertainty is not likely to be the source of the discrepancy.
△ Less
Submitted 6 March, 2019;
originally announced March 2019.
-
Zemach moments and radii of 2,3H and 3,4He
Authors:
N. Nevo Dinur,
O. J. Hernandez,
S. Bacca,
N. Barnea,
C. Ji,
S. Pastore,
M. Piarulli,
R. B. Wiringa
Abstract:
We present benchmark calculations of Zemach moments and radii of 2,3H and 3,4He using various few-body methods. Zemach moments are required to interpret muonic atom data measured by the CREMA collaboration at the Paul Scherrer Institute. Conversely, radii extracted from spectroscopic measurements can be compared with ab initio computations, posing stringent constraints on the nuclear model. For a…
▽ More
We present benchmark calculations of Zemach moments and radii of 2,3H and 3,4He using various few-body methods. Zemach moments are required to interpret muonic atom data measured by the CREMA collaboration at the Paul Scherrer Institute. Conversely, radii extracted from spectroscopic measurements can be compared with ab initio computations, posing stringent constraints on the nuclear model. For a given few-body method, different numerical procedures can be applied to compute these quantities. A detailed analysis of the numerical uncertainties entering the total theoretical error is presented. Uncertainties from the few-body method and the calculational procedure are found to be smaller than the dependencies on the dynamical modeling and the single nucleon inputs, which are found to be <= 2%. When relativistic corrections and two-body currents are accounted for, the calculated moments and radii are in very good agreement with the available experimental data.
△ Less
Submitted 26 December, 2018;
originally announced December 2018.
-
Electromagnetic Excitations and Responses in Nuclei from First Principles
Authors:
R. B. Baker,
K. D. Launey,
N. Nevo Dinur,
S. Bacca,
J. P. Draayer,
T. Dytrych
Abstract:
We discuss the role of clustering on monopole, dipole, and quadrupole excitations in nuclei in the framework of the ab initio symmetry-adapted no-core shell model (SA-NCSM). The SA-NCSM starts from nucleon-nucleon potentials and, by exploring symmetries known to dominate the nuclear dynamics, can reach nuclei up through the calcium region by accommodating ultra-large model spaces critical to descr…
▽ More
We discuss the role of clustering on monopole, dipole, and quadrupole excitations in nuclei in the framework of the ab initio symmetry-adapted no-core shell model (SA-NCSM). The SA-NCSM starts from nucleon-nucleon potentials and, by exploring symmetries known to dominate the nuclear dynamics, can reach nuclei up through the calcium region by accommodating ultra-large model spaces critical to descriptions of clustering and collectivity. The results are based on calculations of electromagnetic sum rules and discretized responses using the Lanczos algorithm, that can be used to determine response functions, and for 4He are benchmarked against exact solutions of the hyperspherical harmonics method. In particular, we focus on He, Be, and O isotopes, including giant resonances and monopole sum rules.
△ Less
Submitted 19 December, 2018;
originally announced December 2018.
-
Ab initio calculation of nuclear structure corrections in muonic atoms
Authors:
Chen Ji,
Sonia Bacca,
Nir Barnea,
Oscar Javier Hernandez,
Nir Nevo-Dinur
Abstract:
The measurement of the Lamb shift in muonic hydrogen and the subsequent emergence of the proton-radius puzzle have motivated an experimental campaign devoted to measuring the Lamb shift in other light muonic atoms, such as muonic deuterium and helium. For these systems it has been shown that two-photon exchange nuclear structure corrections are the largest source of uncertainty and consequently th…
▽ More
The measurement of the Lamb shift in muonic hydrogen and the subsequent emergence of the proton-radius puzzle have motivated an experimental campaign devoted to measuring the Lamb shift in other light muonic atoms, such as muonic deuterium and helium. For these systems it has been shown that two-photon exchange nuclear structure corrections are the largest source of uncertainty and consequently the bottleneck for exploiting the experimental precision to extract the nuclear charge radius. Utilizing techniques and methods developed to study electromagnetic reactions in light nuclei, recent calculations of nuclear structure corrections to the muonic Lamb shift have reached unprecedented precision, reducing the uncertainty with respect to previous estimates by a factor of 5 in certain cases. These results will be useful for shedding light on the nature of the proton-radius puzzle and other open questions pertaining to it. Here, we review and update calculations for muonic deuterium and tritium atoms, and for muonic helium-3 and helium-4 ions. We present a thorough derivation of the formalism and discuss the results in relation to other approaches where available. We also describe how to assess theoretical uncertainties, for which the language of chiral effective field theory furnishes a systematic approach that could be further exploited in the future.
△ Less
Submitted 8 June, 2018;
originally announced June 2018.
-
Computing the dipole polarizability of 48Ca with increased precision
Authors:
Mirko Miorelli,
Sonia Bacca,
Gaute Hagen,
Thomas Papenbrock
Abstract:
We compute the electric dipole polarizability of 48Ca with an increased precision by including more correlations than in previous studies. Employing the coupled-cluster method we go beyond singles and doubles excitations and include leading-order three-particle-three-hole (3p-3h) excitations for the ground state, excited states, and the similarity transformed operator. We study electromagnetic sum…
▽ More
We compute the electric dipole polarizability of 48Ca with an increased precision by including more correlations than in previous studies. Employing the coupled-cluster method we go beyond singles and doubles excitations and include leading-order three-particle-three-hole (3p-3h) excitations for the ground state, excited states, and the similarity transformed operator. We study electromagnetic sum rules, such as the bremsstrahlung sum rule m_0 and the polarizability sum rule alpha_D using interactions from chiral effective field theory. To gauge the quality of our coupled-cluster approximations we perform several benchmarks with the effective interaction hyperspherical harmonics approach in 4He and with self consistent Green's function in 16O. We compute the dipole polarizability of 48Ca employing the chiral interaction N2LOsat [Ekstroem et al., Phys. Rev. C 91, 051301 (2015)] and the 1.8/2.0 (EM) [Hebeler et al., Phys. Rev. C 83, 031301 (2011)]. We find that the effect of 3p-3h excitations in the ground state is small for 1.8/2.0 (EM) but non-negligible for N2LOsat. The addition of these new correlations allows us to improve the precision of our 48Ca calculations and reconcile the recently reported discrepancy between coupled-cluster results based on these interactions and the experimentally determined alpha_D from proton inelastic scattering in 48Ca [Birkhan et al., Phys. Rev. Lett. 118, 252501 (2017)]. For the computation of electromagnetic and polarizability sum rules, the inclusion of leading-order 3p-3h excitations in the ground state is important, while less so for the excited states and the similarity-transformed dipole operator.
△ Less
Submitted 29 June, 2018; v1 submitted 5 April, 2018;
originally announced April 2018.
-
Recent developments in nuclear structure theory: an outlook on the muonic atom program
Authors:
Oscar Javier Hernandez,
Sonia Bacca,
Kyle Andrew Wendt
Abstract:
The discovery of the proton-radius puzzle and the subsequent deuteron-radius puzzle is fueling an on-going debate on possible explanations for the difference in the observed radii obtained from muonic atoms and from electron-nucleus systems. Atomic nuclei have a complex internal structure that must be taken into account when analyzing experimental spectroscopic results. Ab initio nuclear structure…
▽ More
The discovery of the proton-radius puzzle and the subsequent deuteron-radius puzzle is fueling an on-going debate on possible explanations for the difference in the observed radii obtained from muonic atoms and from electron-nucleus systems. Atomic nuclei have a complex internal structure that must be taken into account when analyzing experimental spectroscopic results. Ab initio nuclear structure theory provided the so far most precise estimates of important corrections to the Lamb shift in muonic atoms and is well poised to also investigate nuclear structure corrections to the hyperfine splitting in muonic atoms. Independently on whether the puzzle is due to beyond-the-standard-model physics or not, nuclear structure corrections are a necessary theoretical input to any experimental extraction of electric and magnetic radii from precise muonic atom measurements.
Here, we review the status of the calculations performed by the TRIUMF-Hebrew University group, focusing on the deuteron, and discuss preliminary results on magnetic sum rules calculated with two-body currents at next-to-leading order. Two-body currents will be an important ingredient in future calculations of nuclear structure corrections to the hyperfine splitting in muonic atoms.
△ Less
Submitted 14 December, 2017;
originally announced December 2017.
-
The deuteron-radius puzzle is alive: A new analysis of nuclear structure uncertainties
Authors:
Oscar Javier Hernandez,
Andreas Ekström,
Nir Nevo Dinur,
Chen Ji,
Sonia Bacca,
Nir Barnea
Abstract:
To shed light on the deuteron radius puzzle we analyze the theoretical uncertainties of the nuclear structure corrections to the Lamb shift in muonic deuterium. We find that the discrepancy between the calculated two-photon exchange correction and the corresponding experimentally inferred value by Pohl et al. [1] remain. The present result is consistent with our previous estimate, although the dis…
▽ More
To shed light on the deuteron radius puzzle we analyze the theoretical uncertainties of the nuclear structure corrections to the Lamb shift in muonic deuterium. We find that the discrepancy between the calculated two-photon exchange correction and the corresponding experimentally inferred value by Pohl et al. [1] remain. The present result is consistent with our previous estimate, although the discrepancy is reduced from 2.6 $σ$ to 2 $σ$. The error analysis includes statistic as well as systematic uncertainties stemming from the use of nucleon-nucleon interactions derived from chiral effective field theory at various orders. We therefore conclude that nuclear theory uncertainty is more likely not the source of the discrepancy.
△ Less
Submitted 11 June, 2018; v1 submitted 3 November, 2017;
originally announced November 2017.
-
Electromagnetic reactions from coupled-cluster theory
Authors:
Sonia Bacca,
Mirko Miorelli,
Gaute Hagen
Abstract:
We review recent results for electromagnetic reactions and related sum rules in light and medium-mass nuclei obtained from coupled-cluster theory. In particular, we highlight our recent computations of the photodisintegration cross section of 40Ca and of the electric dipole polarizability for oxygen and calcium isotopes. We also provide new results for the Coulomb sum rule for 4He and 16O. For 4He…
▽ More
We review recent results for electromagnetic reactions and related sum rules in light and medium-mass nuclei obtained from coupled-cluster theory. In particular, we highlight our recent computations of the photodisintegration cross section of 40Ca and of the electric dipole polarizability for oxygen and calcium isotopes. We also provide new results for the Coulomb sum rule for 4He and 16O. For 4He we perform a thorough comparison of coupled-cluster theory with exact hyperspherical harmonics.
△ Less
Submitted 26 October, 2017;
originally announced October 2017.
-
Electric dipole polarizability of $^{48}$Ca and implications for the neutron skin
Authors:
J. Birkhan,
M. Miorelli,
S. Bacca,
S. Bassauer,
C. A. Bertulani,
G. Hagen,
H. Matsubara,
P. von Neumann-Cosel,
T. Papenbrock,
N. Pietralla,
V. Yu. Ponomarev,
A. Richter,
A. Schwenk,
A. Tamii
Abstract:
The electric dipole strength distribution in Ca-48 between 5 and 25 MeV has been determined at RCNP, Osaka, from proton inelastic scattering experiments at forward angles. Combined with photoabsorption data at higher excitation energy, this enables for the first time the extraction of the electric dipole polarizability alpha_D(Ca-48) = 2.07(22) fm^3. Remarkably, the dipole response of Ca-48 is fou…
▽ More
The electric dipole strength distribution in Ca-48 between 5 and 25 MeV has been determined at RCNP, Osaka, from proton inelastic scattering experiments at forward angles. Combined with photoabsorption data at higher excitation energy, this enables for the first time the extraction of the electric dipole polarizability alpha_D(Ca-48) = 2.07(22) fm^3. Remarkably, the dipole response of Ca-48 is found to be very similar to that of Ca-40, consistent with a small neutron skin in Ca-48. The experimental results are in good agreement with ab initio calculations based on chiral effective field theory interactions and with state-of-the-art density-functional calculations, implying a neutron skin in Ca-48 of 0.14 - 0.20 fm.
△ Less
Submitted 4 May, 2017; v1 submitted 21 November, 2016;
originally announced November 2016.
-
Update on nuclear structure effects in light muonic atoms
Authors:
Oscar Javier Hernandez,
Nir Nevo Dinur,
Chen Ji,
Sonia Bacca,
Nir Barnea
Abstract:
We present calculations of the nuclear structure corrections to the Lamb shift in light muonic atoms, using state-of-the-art nuclear potentials. We outline updated results on finite nucleon size contributions.
We present calculations of the nuclear structure corrections to the Lamb shift in light muonic atoms, using state-of-the-art nuclear potentials. We outline updated results on finite nucleon size contributions.
△ Less
Submitted 21 April, 2016;
originally announced April 2016.
-
Structure models: from shell model to ab initio methods
Authors:
Sonia Bacca
Abstract:
A brief review of models to describe nuclear structure and reactions properties is presented, starting from the historical shell model picture and encompassing modern ab initio approaches. A selection of recent theoretical results on observables for exotic light and medium-mass nuclei is shown. Emphasis is given to the comparison with experiment and to what can be learned about three-body forces a…
▽ More
A brief review of models to describe nuclear structure and reactions properties is presented, starting from the historical shell model picture and encompassing modern ab initio approaches. A selection of recent theoretical results on observables for exotic light and medium-mass nuclei is shown. Emphasis is given to the comparison with experiment and to what can be learned about three-body forces and continuum properties.
△ Less
Submitted 19 April, 2016;
originally announced April 2016.