Showing 1–50 of 144 results for author: Kievsky, A

Emulation of Proton-Deuteron Scattering via the Reduced Basis Method and Active Learning: Detailed Description

Authors: Alex Gnech, Xilin Zhang, Christian Drischler, R. J. Furnstahl, Alessandro Grassi, Alejandro Kievsky, Laura E. Marcucci, Michele Viviani

Abstract: Nucleon-deuteron ($Nd$) scattering can be used to constrain three-nucleon forces in chiral effective field theory ($χ$EFT). However, high-fidelity calculations, such as the Hyperspherical Harmonic (HH) method, are computationally expensive, making it difficult or even prohibitive to explore the vast parameter space of $χ$EFT\xspace. To address this challenge, specifically for proton-deuteron (… ▽ More Nucleon-deuteron ($Nd$) scattering can be used to constrain three-nucleon forces in chiral effective field theory ($χ$EFT). However, high-fidelity calculations, such as the Hyperspherical Harmonic (HH) method, are computationally expensive, making it difficult or even prohibitive to explore the vast parameter space of $χ$EFT\xspace. To address this challenge, specifically for proton-deuteron ($pd$) scattering below the deuteron breakup threshold, we developed model-driven emulators based on the Reduced Basis Method (RBM) and active learning techniques, as presented in \href{https://arxiv.org/abs/2511.01844}{arXiv:2511.01844}. The method exploits the similarities between solutions at different parameter points to significantly reduce computational costs. In this companion paper, we provide a comprehensive description of our HH-based high-fidelity calculations and implementation of both variational-method-based and Galerkin-projection-based scattering emulators. We demonstrate the effectiveness of active learning in the form of greedy algorithms for selecting optimal training points in the parameter space, and the high accuracy and speed of the emulators, for two different nucleon forces and two scattering channels (${1/2}^+$ and ${1/2}^-$). For example, in a two-dimensional parameter space, the relative emulation errors can be reduced to $10^{-7}$ with fewer than 10 training points. Our work paves the way for the efficient calibration of $χ$EFT\xspace nucleon interactions using Bayesian statistics, and the methodology can be applied to other nuclear scattering processes (including neutron-deuteron scattering), as well as other finite quantum systems. △ Less

Submitted 13 November, 2025; originally announced November 2025.

Comments: 24 page, 16 figures, 1 table

Accurate and Efficient Emulation of Proton-Deuteron Scattering via the Reduced Basis Method and Active Learning

Authors: Alex Gnech, Xilin Zhang, Christian Drischler, R. J. Furnstahl, Alessandro Grassi, Alejandro Kievsky, Laura E. Marcucci, Michele Viviani

Abstract: We introduce highly accurate and efficient emulators for proton-deuteron scattering below the deuteron breakup threshold. We explore two different reduced-basis method strategies: one based on the Kohn variational principle and another on Galerkin projections of the underlying system of linear equations. We use the adaptive greedy algorithm previously developed for two-body scattering for optimal… ▽ More We introduce highly accurate and efficient emulators for proton-deuteron scattering below the deuteron breakup threshold. We explore two different reduced-basis method strategies: one based on the Kohn variational principle and another on Galerkin projections of the underlying system of linear equations. We use the adaptive greedy algorithm previously developed for two-body scattering for optimal selection of high-fidelity training points in the input parameter space. We demonstrate that these emulators reproduce ab initio hyperspherical harmonics calculations of $R$-matrix elements with remarkable precision, achieving relative errors as low as $10^{-7}$ with a small number of training points, even in regions of strong nonlinear parameter dependence. They also dramatically accelerate the exploration of the scattering predictions in the parameter space, a capability highly desired for calibrating (chiral) three-nucleon forces against scattering measurements. Our formalism can be further generalized to handle nucleon-deuteron scattering above the breakup threshold. These emulator developments will provide valuable tools to accelerate uncertainty quantification and rigorous parameter inference in the study of nuclear forces. △ Less

Submitted 3 November, 2025; originally announced November 2025.

Comments: 7 pages, 2 figures

The excited state of the $α$-particle: a benchmark study

Authors: Pierre-Yves Duerinck, Arnoldas Deltuva, Jérémy Dohet-Eraly, Mario Gattobigio, Alejandro Kievsky, Rimantas Lazauskas, Darius Likandrovas, Michele Viviani

Abstract: A benchmark study is performed for the excited state of $^4$He. When the Coulomb interaction is switched off, the $^4$He nucleus exhibits a bound excited state in the vicinity of $p-{}^3$H threshold. As the Coulomb interaction is gradually introduced, the excited state crosses the threshold and eventually becomes a resonant state. Using three numerical methods, we track the evolution of this exc… ▽ More A benchmark study is performed for the excited state of $^4$He. When the Coulomb interaction is switched off, the $^4$He nucleus exhibits a bound excited state in the vicinity of $p-{}^3$H threshold. As the Coulomb interaction is gradually introduced, the excited state crosses the threshold and eventually becomes a resonant state. Using three numerical methods, we track the evolution of this excited state and determine the resonance energy and width. Comparisons of the theoretical predictions reveal a significant discrepancy with commonly used $R$-matrix values based on the analysis of the experimental data. We explain the origin for this discrepancy. Additionally, the two-level energy spectrum of $^4$He in the absence of the Coulomb force exhibits characteristics linked to Efimov physics, suggesting a reduced sensitivity to interaction details. △ Less

Submitted 11 August, 2025; originally announced August 2025.

Convergence of the $ppp$ correlation function within the hyperspherical adiabatic basis

Authors: E. Garrido, A. Kievsky, R. Del Grande, L. Serksnyte, M. Viviani, L. E. Marcucci

Abstract: The computation of the three-particle correlation function involving three hadrons started just recently after the first publications of ALICE measurements. Key elements to be considered are the correct description of the asymptotics, antisymmetrization issues and, in most cases, the treatment of the Coulomb interaction. In the case of the $ppp$ correlation function, a first analysis was done wher… ▽ More The computation of the three-particle correlation function involving three hadrons started just recently after the first publications of ALICE measurements. Key elements to be considered are the correct description of the asymptotics, antisymmetrization issues and, in most cases, the treatment of the Coulomb interaction. In the case of the $ppp$ correlation function, a first analysis was done where the hyperspherical adiabatic method was used to determine the $ppp$ wave function at different energies. Although the asymptotic behavior, antisymmetrization issues and the treatment of the Coulomb interaction were discussed in detail, the convergence properties of the adiabatic basis were studied at low energies around the formation of the correlation peak determined mainly by the $J^π=1/2^-$ and $3/2^-$ three-body states. Since many and very precise data have been taken or are planned to be measured at energies beyond the peak, we present an analysis of the convergence characteristics of the basis as the energy of the process increases. We show that in order to describe correctly the correlation tail it is necessary to consider three-body states up to $J^π=21/2^-$ whereas higher states can be considered as free. Once those states are incorporated solving the associate dynamical equations, the agreement with the experimental data is found to be excellent. △ Less

Submitted 28 May, 2025; originally announced May 2025.

Comments: Submitted for publication, six pages, three figures

Nucleon-nucleon correlation functions from different interactions in comparison

Authors: Matthias Göbel, Alejandro Kievsky

Abstract: Correlation functions as they can be observed in heavy-ion collisions using the femtoscopy technique are a powerful tool to study the interaction among different baryons or mesons. Specifically, the multi-nucleon correlation functions have been under intense experimental and theoretical investigation in the recent years. Due to the interest of using this observable as an input in the construction… ▽ More Correlation functions as they can be observed in heavy-ion collisions using the femtoscopy technique are a powerful tool to study the interaction among different baryons or mesons. Specifically, the multi-nucleon correlation functions have been under intense experimental and theoretical investigation in the recent years. Due to the interest of using this observable as an input in the construction of potentials between hadrons we revisit the nucleon-nucleon correlation function and calculate it using different nuclear interactions at high precision. Since the nucleon-nucleon potential is determined to reproduce the two-nucleon scattering data, we would like to critically evaluate the amount of this information captured by the correlation function. We study the dependence of the correlations on the nuclear force giving detailed insights into the calculations, in particular the convergence behavior in the partial waves. The coupling between the different partial-wave channels is taken into account and the relevance of this effect is quantified. To make contact with precedent studies the results based on the Argonne V18 interaction are presented. Then we consider also the Norfolk NV2-IIa and NV2-IIb chiral EFT interactions. The analysis of the differences between the correlations of the various interactions shows that for momenta between 0 and 500 MeV there are variations of up to 5.9 % for the $nn$ system, of up to 1.8 % for the $np$ system, and of 1.4 % for the $pp$ system. △ Less

Submitted 19 May, 2025; originally announced May 2025.

Comments: Paper + supplementary material: 9+3 pages, 10+5 figures

ToMCCA-3: A realistic 3-body coalescence model

Authors: Maximilian Mahlein, Bhawani Singh, Michele Viviani, Francesca Bellini, Laura Fabbietti, Alejandro Kievsky, Laura Elisa Marcucci

Abstract: The formation of light nuclei in high-energy collisions provides valuable insights into the underlying dynamics of the strong interaction and the structure of the particle-emitting source. Understanding this process is crucial not only for nuclear physics but also for astrophysical studies, where the production of rare antinuclei could serve as a probe for new physics. This work presents a three-b… ▽ More The formation of light nuclei in high-energy collisions provides valuable insights into the underlying dynamics of the strong interaction and the structure of the particle-emitting source. Understanding this process is crucial not only for nuclear physics but also for astrophysical studies, where the production of rare antinuclei could serve as a probe for new physics. This work presents a three-body coalescence model based on the Wigner function formalism, offering a refined description of light-nucleus production. By incorporating realistic two- and three-body nuclear interaction potentials constrained by modern scattering and femtoscopic correlation data, our approach improves on traditional coalescence models. The framework is validated using event generators applied to proton-proton collisions at $\sqrt{s}=13$ TeV to predict the momentum spectra of light (anti) nuclear nuclei with mass number $A=3$, which are then compared with the experimental data from ALICE. Our results demonstrate the sensitivity of light nucleus yields to the choice of nuclear wave functions, emphasizing the importance of an accurate description of the coalescence process. This model lays the foundation for the extension of coalescence studies of $A=3$ light nuclei to a wider range of collision systems and energies. △ Less

Submitted 3 April, 2025; originally announced April 2025.

Comments: 23 pages, 4 Figures

The 3-$α$ and 4-$α$ particle systems within short-range Effective Field Theory

Authors: E. Filandri, M. Viviani, L. Girlanda, A. Kievsky, L. E. Marcucci

Abstract: ${}^{12}{\rm C}$ and ${}^{16}{\rm O}$ nuclei represent essential elements for life on Earth. The study of their formation plays a key role in understanding heavy element nucleosynthesis and stellar evolution. In this paper we present the study of ${}^{12}{\rm C}$ and ${}^{16}{\rm O}$ nuclei as systems composed of $α… ▽ More ${}^{12}{\rm C}$ and ${}^{16}{\rm O}$ nuclei represent essential elements for life on Earth. The study of their formation plays a key role in understanding heavy element nucleosynthesis and stellar evolution. In this paper we present the study of ${}^{12}{\rm C}$ and ${}^{16}{\rm O}$ nuclei as systems composed of $α$-particle clusters using the short-range effective field theory approach. The fundamental and excited states of the studied nuclei are calculated within an ab-initio approach, using the Hyperspherical Harmonics method. Thanks to the two-body potential and fine-tuning of the three-body force, we have found the ${}^{12}{\rm C}$ system nicely reproduced by theory. However, for the ${}^{16}{\rm O}$ case, it is necessary to include a 4-body force in order to achieve agreement with the experimental data. △ Less

Submitted 24 January, 2025; originally announced January 2025.

The X17 boson and the d(p,e+ e-)3He and d(n,e+ e-)3H processes: a theoretical analysis

Authors: M. Viviani, E. Filandri, L. Girlanda, C. Gustavino, A. Kievsky, L. E. Marcucci

Abstract: The present work deals with the e+ e- pair production in the d(p,e+e-)3He and d(n,e+ e-)3H processes, in order to evidentiate possible effects due to the exchange of a hypothetical low-mass boson, the so-called X17. These processes are studied for energies of the incident beams in the range 18-30 MeV, in order to have a sufficient energy to produce such a boson, whose mass is estimated to be aroun… ▽ More The present work deals with the e+ e- pair production in the d(p,e+e-)3He and d(n,e+ e-)3H processes, in order to evidentiate possible effects due to the exchange of a hypothetical low-mass boson, the so-called X17. These processes are studied for energies of the incident beams in the range 18-30 MeV, in order to have a sufficient energy to produce such a boson, whose mass is estimated to be around 17 MeV. We first analyze them as a purely electromagnetic processes, in the context of a state-of-the-art approach to nuclear strong-interaction dynamics and nuclear electromagnetic currents, derived from chiral effective field theory chiEFT. Next, we examine how the exchange of a hypothetical low-mass boson would impact the cross sections for such processes. We consider several possibilities, that this boson is either a scalar, pseudoscalar, vector, or axial particle. The main aim of the study is to exploit the specular structure of the 3He and 3H nuclei to investigate the isospin dependency of the X17-nucleon interaction, as the alleged "proto-phobicity". △ Less

Submitted 15 April, 2025; v1 submitted 29 August, 2024; originally announced August 2024.

Comments: 19 pages, 15 figures, new version with more details, added two appendices. arXiv admin note: text overlap with arXiv:2104.07808

Journal ref: Phys. Rev. C 111, 034002 (2025)

The $pΛ$ and $ppΛ$ correlation functions

Authors: E. Garrido, A. Kievsky, M. Gattobigio, M. Viviani, L. E. Marcucci, R. Del Grande, L. Fabbietti, D. Melnichenko

Abstract: In this work we present the study of $pΛ$ and $ppΛ$ scattering processes using femtoscopic correlation functions. This observable has been recently used to access the low-energy interaction of hadrons emitted in the final state of high-energy collisions, delivering unprecedented precision information of the interaction among strange hadrons. The formalism for particle pairs is well established and… ▽ More In this work we present the study of $pΛ$ and $ppΛ$ scattering processes using femtoscopic correlation functions. This observable has been recently used to access the low-energy interaction of hadrons emitted in the final state of high-energy collisions, delivering unprecedented precision information of the interaction among strange hadrons. The formalism for particle pairs is well established and it relates the measured correlation functions with the scattering wave function and the emission source. In the present work we analyze the $NNΛ$ scattering in free space and relate the corresponding wave function to the $ppΛ$ correlation measurement performed by the ALICE collaboration. The three-body problem is solved using the hyperspherical adiabatic basis. Regarding the $pΛ$ and $ppΛ$ interactions, different models are used and their impact on the correlation function is studied. The three body force considered in this work is anchored to describe the binding energy of the hypertriton and to give a good description of the two four-body hypernuclei. As a main result we have observed a huge, low-energy peak in the $ppΛ$ correlation function, mainly produced by the $J^π=1/2^+$ three-body state. The study of this peak from an experimental as well as a theoretical point of view will provide important constraints to the two- and three-body interactions. △ Less

Submitted 3 August, 2024; originally announced August 2024.

Comments: Submitted for publication

Study of the alpha-particle monopole transition form factor

Authors: M. Viviani, A. Kievsky, L. E. Marcucci, L. Girlanda

Abstract: The 4He monopole form factor is studied by computing the transition matrix element of the electromagnetic charge operator between the 4He ground-state and the p+3H and n+3He scattering states. The nuclear wave functions are calculated using the hyperspherical harmonic method, by starting from Hamiltonians including two- and three-body forces derived in chiral effective field theory. The electromag… ▽ More The 4He monopole form factor is studied by computing the transition matrix element of the electromagnetic charge operator between the 4He ground-state and the p+3H and n+3He scattering states. The nuclear wave functions are calculated using the hyperspherical harmonic method, by starting from Hamiltonians including two- and three-body forces derived in chiral effective field theory. The electromagnetic charge operator retains, beyond the leading order (impulse approximation) term, also higher order contributions, as relativistic corrections and meson-exchange currents. The results for the monopole form factor are in fairly agreement with recent MAMI data. Comparison with other theoretical calculations are also provided. △ Less

Submitted 2 June, 2024; originally announced June 2024.

Comments: 22 pages, 6 figures

The 3He(\vec n,p)3H parity-conserving asymmetry

Authors: M. Viviani, S. Baeßler, L. Barrón-Palos, N. Birge, J. D. Bowman, J. Calarco, V. Cianciolo, C. E. Coppola, C. B. Crawford, G. Dodson, N. Fomin, I. Garishvili, M. T. Gericke, L. Girlanda, G. L. Greene, G. M. Hale, J. Hamblen, C. Hayes, E. B. Iverson, M. L. Kabir, A. Kievsky, L. E. Marcucci, M. McCrea, E. Plemons, A. Ramírez-Morales , et al. (6 additional authors not shown)

Abstract: Recently, the n$^3$He collaboration reported a measurement of the parity-violating (PV) proton directional asymmetry $A_{\mathrm {PV}} = (1.55\pm 0.97~\mathrm {(st\ at)} \pm 0.24~\mathrm {(sys)})\times 10^{-8}$ in the capture reaction of ${}^3$He$(\vec {n},{\mathrm p}){}^3$H at meV incident neutron energies. The result increased the limited inventory of precisely measured and calculable PV observa… ▽ More Recently, the n$^3$He collaboration reported a measurement of the parity-violating (PV) proton directional asymmetry $A_{\mathrm {PV}} = (1.55\pm 0.97~\mathrm {(st\ at)} \pm 0.24~\mathrm {(sys)})\times 10^{-8}$ in the capture reaction of ${}^3$He$(\vec {n},{\mathrm p}){}^3$H at meV incident neutron energies. The result increased the limited inventory of precisely measured and calculable PV observables in few-body systems required to further understand the structure of hadronic weak interaction. In this letter, we report the experimental and theoretical investigation of a parity conserving (PC) asymmetry $A_{\mathrm {PC}}$ in the same reaction (the first ever measured PC observable at meV neutron energies). As a result of S- and P-wave mixing in the reaction, the $A_{\mathrm {PC}}$ is inversely proportional to the neutron wavelength $λ$. The experimental value is $(λ\times A_{\mathrm {PC}})\equivβ= (-1.97 \pm 0.28~\mathrm{(stat)}\pm 0.12~\mathrm{(sys)}) \times 10^{-6}$ Amstrongs. We present results for a theoretical analysis of this reaction by solving the four-body scattering problem within the hyperspherical harmonic method. We find that in the ${}^3$He$(\vec {n},{\mathrm p}){}^3$H reaction, $A_{\mathrm {PC}}$ depends critically on the energy and width of the close $0^-$ resonant state of ${}^4$He, resulting in a large sensitivity to the spin-orbit components of the nucleon-nucleon force and even to the three-nucleon force. The analysis of the accurately measured $A_{\mathrm {PC}}$ and $A_{\mathrm {PV}}$ using the same few-body theoretical models gives essential information needed to interpret the PV asymmetry in the ${}^3$He$(\vec {n}, {\mathrm p}){}^3$H reaction. △ Less

Submitted 16 May, 2024; originally announced May 2024.

Comments: 8 pages, 3 figures

Momentum dependent nucleon-nucleon contact interactions and their effect on p-d scattering observables

Authors: E. Filandri, L. Girlanda, A. Kievsky, L. E. Marcucci, M. Viviani

Abstract: Starting from a complete set of relativistic nucleon-nucleon contact operators up to order $O(p^4)$ of the expansion in the soft (relative or nucleon) momentum $p$, we show that non-relativistic expansions of relativistic operators involve twenty-six independent combinations, two starting at $O(p^0)$, seven at order $O(p^2)$ and seventeen at order $O(p^4)$. This demonstrates the existence of two l… ▽ More Starting from a complete set of relativistic nucleon-nucleon contact operators up to order $O(p^4)$ of the expansion in the soft (relative or nucleon) momentum $p$, we show that non-relativistic expansions of relativistic operators involve twenty-six independent combinations, two starting at $O(p^0)$, seven at order $O(p^2)$ and seventeen at order $O(p^4)$. This demonstrates the existence of two low-energy free constants that parameterize interactions dependent on the total momentum of the pair of nucleons $P$. The latter, through the use of a unitary transformation, can be removed in the two-nucleon fourth-order contact interaction of the Chiral Effective Field Theory, generating a three-nucleon interaction at the same order. Within a hybrid approach in which this interaction is considered together with the phenomenological potential AV18, we show that the LECs involved can be used to fit very accurate data on the polarization observables of the low-energy $p-d$ scattering, in particular the $A_y$ asymmetry. △ Less

Submitted 17 April, 2024; originally announced April 2024.

Comments: arXiv admin note: text overlap with arXiv:2302.03468

The $ppp$ correlation function with a screened Coulomb potential

Authors: A. Kievsky, E. Garrido, M. Viviani, M. Gattobigio

Abstract: The correlation function is a useful tool to study the interaction between hadrons. The theoretical description of this observable requires the knowledge of the scattering wave function, whose asymptotic part is distorted when two or more particles are charged. For a system of three (or more) particles, with more than two particles asymptotically free and at least two of them charged, the asymptot… ▽ More The correlation function is a useful tool to study the interaction between hadrons. The theoretical description of this observable requires the knowledge of the scattering wave function, whose asymptotic part is distorted when two or more particles are charged. For a system of three (or more) particles, with more than two particles asymptotically free and at least two of them charged, the asymptotic part of the wave function is not known in a closed form. In the present study we introduce a screened Coulomb potential and analyze the impact of the screening radius on the correlation function. As we will show, when a sufficiently large screening radius is used, the correlation function results almost unchanged if compared to the case in which the unscreened Coulomb potential is used. This fact allows the use of free asymptotic matching conditions in the solution of the scattering equation simplifying noticeably the calculation of the correlation function. As an illustration we discuss the $pp$ and $ppp$ correlation functions. △ Less

Submitted 9 February, 2024; originally announced February 2024.

Comments: 14 pages, 4 figures, belonging to the collection Advances on Few-Body Problems in Physics, Few-Body Systems

The $nnn$ and $ppp$ correlation functions

Authors: A. Kievsky, E. Garrido, M. Viviani, L. E. Marcucci, L. Serksnyte, R. Del Grande

Abstract: Scattering experiments with three free nucleons in the ingoing channel are extremely challenging in terrestrial laboratories. Recently, the ALICE Collaboration has successfully measured the scattering of three protons indirectly, by using the femtoscopy method in high-energy proton-proton collisions at the Large Hadron Collider. In order to establish a connection with current and future measuremen… ▽ More Scattering experiments with three free nucleons in the ingoing channel are extremely challenging in terrestrial laboratories. Recently, the ALICE Collaboration has successfully measured the scattering of three protons indirectly, by using the femtoscopy method in high-energy proton-proton collisions at the Large Hadron Collider. In order to establish a connection with current and future measurements of femtoscopic three-particle correlation functions, we analyse the scenarios involving $nnn$ and $ppp$ systems using the hyperspherical adiabatic basis. The correlation function is a convolution of the source function and the corresponding scattering wave function. The finite size of the source allows for the use of the free scattering wave function in most of the adiabatic channels except the lowest ones. The scattering wave function has been computed using two different potential models: $(i)$ a spin-dependent Gaussian potential with parameters fixed to reproduce the scattering length and effective range and $(ii)$ the Argonne $v_{18}$ nucleon-nucleon interaction. Moreover, in the case of three protons, the Coulomb interaction has been considered in its hypercentral form. The results presented here have to be considered as a first step in the description of three-particle correlation functions using the hyperspherical adiabatic basis, opening the door to the investigation of other systems, such as the $ppΛ$ system. For completeness, the comparison with the measurement by the ALICE Collaboration is shown assuming different values of the source radius. △ Less

Submitted 3 July, 2025; v1 submitted 16 October, 2023; originally announced October 2023.

Comments: 18 pages, 12 figures

Role of three-body dynamics in nucleon-deuteron correlation functions

Authors: M. Viviani, S. König, A. Kievsky, L. E. Marcucci, B. Singh, O. Vázquez Doce

Abstract: Correlation functions of hadrons can be accessed in high-energy collisions of atomic nuclei, revealing information about the underlying interaction. This work complements experimental efforts to study nucleon-deuteron $Nd$ -- with $N=p$ (proton) or $N=n$ (neutron) -- correlations with theory evaluations using different techniques. The correlation functions $C_{nd}$ and $C_{pd}$ are calculated base… ▽ More Correlation functions of hadrons can be accessed in high-energy collisions of atomic nuclei, revealing information about the underlying interaction. This work complements experimental efforts to study nucleon-deuteron $Nd$ -- with $N=p$ (proton) or $N=n$ (neutron) -- correlations with theory evaluations using different techniques. The correlation functions $C_{nd}$ and $C_{pd}$ are calculated based on scattering wave function, extending previous benchmarks for the $Nd$ scattering matrix to this new observable. We use hyperspherical harmonics and Faddeev techniques with one of the widely used nucleon-nucleon ($NN$) interactions, the Argonne $v_{18}$ potential. Moreover, in the low-energy region we perform additional calculations in the framework of pionless effective field theory. The $pd$ correlation function is computed in the large-energy region to make contact with a recent measurement by the ALICE Collaboration. We show that the scattering wave function has the proper dynamical input to describe an initial rise and subsequent oscillations of $C_{pd}$ as a function of the energy. Effects on the observables using different $NN$ and three-nucleon potentials are evaluated with the conclusion that variations of around $2\%$ are observed. Although these effects are small, future measurements can go beyond this accuracy allowing for new detailed studies of strong interaction in light nuclear systems. The present study supports the current efforts devoted to the measurement of correlation functions in systems dominated by the strong interactions, such as $pd$, $ppp$, $Λd$ and $ppΛ$. △ Less

Submitted 4 June, 2023; originally announced June 2023.

Comments: 23 pages, 8 figures

The Fate of excited state of $^4\text{He}$

Authors: Mario Gattobigio, Alejandro Kievsky

Abstract: We investigate the properties of the excited state of $^4\mathrm{He}$, $^4\mathrm{He}^*$, within the framework of Efimov physics and its connection to the unitary point of the nuclear interaction. We explore two different approaches to track the trajectory of $^4\mathrm{He}^*$ as it crosses the $^3\mathrm{H}$+p threshold and potentially becomes a resonant state. The first approach involves an anal… ▽ More We investigate the properties of the excited state of $^4\mathrm{He}$, $^4\mathrm{He}^*$, within the framework of Efimov physics and its connection to the unitary point of the nuclear interaction. We explore two different approaches to track the trajectory of $^4\mathrm{He}^*$ as it crosses the $^3\mathrm{H}$+p threshold and potentially becomes a resonant state. The first approach involves an analytical continuation of the energy with respect to the Coulomb coupling, while the second approach introduces an artificial four-body force that it is gradually released. By utilizing Padé approximants and extrapolation techniques, we estimate the energy and width of the resonance. Our results suggest a central energy value of $E_R=0.060(3)$ MeV and a width of $Γ/2=0.036(6)$ MeV using the Coulomb analysis, and $E_R=0.068(1)$ MeV and $Γ/2=0.007(5)$ MeV with the four-body force analysis. Interestingly, these results are consistent with calculations based on {\it ab-initio} nuclear interactions but differ from the accepted values of the $0^+$ resonance energy and width. This highlights the challenges in accurately determining the properties of resonant states in light nuclei and calls for further investigations and refinements in theoretical approaches. △ Less

Submitted 26 May, 2023; originally announced May 2023.

Effect of the N3LO three-nucleon contact interaction on p-d scattering observables

Authors: L. Girlanda, E. Filandri, A. Kievsky, L. E. Marcucci, M. Viviani

Abstract: A unitary transformation allows to remove redundant terms in the two-nucleon (2N) contact interaction at the fourth order (N3LO) in the low-energy expansion of Chiral Effective Field Theory. In so doing a three-nucleon (3N) interaction is generated. We express its short-range component in terms of five combinations of low-energy constants (LECs) parametrizing the N3LO 2N contact Lagrangian. Within… ▽ More A unitary transformation allows to remove redundant terms in the two-nucleon (2N) contact interaction at the fourth order (N3LO) in the low-energy expansion of Chiral Effective Field Theory. In so doing a three-nucleon (3N) interaction is generated. We express its short-range component in terms of five combinations of low-energy constants (LECs) parametrizing the N3LO 2N contact Lagrangian. Within a hybrid approach, in which this interaction is considered in conjunction with the phenomenological AV18 2N potential, we show that the involved LECs can be used to fit very accurate data on polarization observables of low-energy $p-d$ scattering, in particular the $A_y$ asymmetry. The resulting interaction is of the right order of magnitude for a N3LO contribution. △ Less

Submitted 7 February, 2023; originally announced February 2023.

Comments: 9 pages, 2 figures

Theoretical study of the d(d,p)3H and d(d,n)3He processes at low energies

Authors: M. Viviani, L. Girlanda, A. Kievsky, D. Logoteta, L. E. Marcucci

Abstract: We present a theoretical study of the processes d(d,p)3H and d(d,n)3He at energies of interest for energy production and for big-bang nucleosynthesis. We accurately solve the four body scattering problem using the ab-initio hyperspherical harmonic method, starting from nuclear Hamiltonians which include modern two- and three-nucleon interactions, derived in chiral effective field theory. We report… ▽ More We present a theoretical study of the processes d(d,p)3H and d(d,n)3He at energies of interest for energy production and for big-bang nucleosynthesis. We accurately solve the four body scattering problem using the ab-initio hyperspherical harmonic method, starting from nuclear Hamiltonians which include modern two- and three-nucleon interactions, derived in chiral effective field theory. We report results for the astrophysical factor, the quintet suppression factor, and various single and double polarized observables. An estimate of the "theoretical uncertainty" for all these quantities is provided by varying the cutoff parameter used to regularize the chiral interactions at high momentum. △ Less

Submitted 4 July, 2022; originally announced July 2022.

Comments: 5 pages, 4 figures

Subleading contributions to $N$-boson systems inside the universal window

Authors: Paolo Recchia, Alejandro Kievsky, Luca Girlanda, Mario Gattobigio

Abstract: We study bosonic systems in the regime in which the two-body system has a shallow bound state or, equivalently, a large value of the two-body scattering length. Using the effective field theory framework as a guide, we construct a series of potential terms which have decreasing importance in the description of the binding energy of the systems. The leading order potential terms consist of a two-bo… ▽ More We study bosonic systems in the regime in which the two-body system has a shallow bound state or, equivalently, a large value of the two-body scattering length. Using the effective field theory framework as a guide, we construct a series of potential terms which have decreasing importance in the description of the binding energy of the systems. The leading order potential terms consist of a two-body term, usually attractive, plus a three-body term, usually repulsive; this last term is required to prevent the collapse of systems with more than two particles. At this order, the parametrization of the two-body potential is done to obtain a correct description of the scattering length, which governs the dynamics in this regime, whereas the three-body term fixes a three-body datum. We investigate the role of the cut-off in the leading order description and we extend the exploration beyond the leading order by including the next-to-leading order terms in both, the two- and three-body potentials. We use the requirement of the stability of the N-body system, whose energy is variationally estimated, to introduce the three-body forces. The potential parametrization, as a function of the cut-off, is fixed to describe the energy of 4 He clusters up to seven particles within the expected accuracy. Finally, we also explore the possibility to describe at the same time the atom-dimer scattering length. △ Less

Submitted 13 June, 2022; originally announced June 2022.

Nuclear Forces for Precision Nuclear Physics -- a collection of perspectives

Authors: Ingo Tews, Zohreh Davoudi, Andreas Ekström, Jason D. Holt, Kevin Becker, Raúl Briceño, David J. Dean, William Detmold, Christian Drischler, Thomas Duguet, Evgeny Epelbaum, Ashot Gasparyan, Jambul Gegelia, Jeremy R. Green, Harald W. Grießhammer, Andrew D. Hanlon, Matthias Heinz, Heiko Hergert, Martin Hoferichter, Marc Illa, David Kekejian, Alejandro Kievsky, Sebastian König, Hermann Krebs, Kristina D. Launey , et al. (20 additional authors not shown)

Abstract: This is a collection of perspective pieces contributed by the participants of the Institute of Nuclear Theory's Program on Nuclear Physics for Precision Nuclear Physics which was held virtually from April 19 to May 7, 2021. The collection represents the reflections of a vibrant and engaged community of researchers on the status of theoretical research in low-energy nuclear physics, the challenges… ▽ More This is a collection of perspective pieces contributed by the participants of the Institute of Nuclear Theory's Program on Nuclear Physics for Precision Nuclear Physics which was held virtually from April 19 to May 7, 2021. The collection represents the reflections of a vibrant and engaged community of researchers on the status of theoretical research in low-energy nuclear physics, the challenges ahead, and new ideas and strategies to make progress in nuclear structure and reaction physics, effective field theory, lattice QCD, quantum information, and quantum computing. The contributed pieces solely reflect the perspectives of the respective authors and do not represent the viewpoints of the Institute for Nuclear theory or the organizers of the program. △ Less

Submitted 2 February, 2022; originally announced February 2022.

Comments: Perspective pieces of the virtual INT program 21-1b "Nuclear Forces for Precision Nuclear Physics", 107 pages

Report number: INT-PUB-22-002, LA-UR-22-20419, UMD-PP-022-02, FERMILAB-PUB-22-090-T

Journal ref: Few-Body Systems 63, 67 (2022)

A Many-Body Density Energy Functional

Authors: A. Kievsky, G. Orlandini, M. Gattobigio

Abstract: The Hohenberg-Kohn theorem and the Kohn-Sham equations, which are at the basis of the Density Functional Theory, are reformulated in terms of a particular many-body density, which is translational invariant and therefore is relevant for self-bound systems. In a similar way that there is a unique relation between the one-body density and the external potential that gives rise to it, we demonstrate… ▽ More The Hohenberg-Kohn theorem and the Kohn-Sham equations, which are at the basis of the Density Functional Theory, are reformulated in terms of a particular many-body density, which is translational invariant and therefore is relevant for self-bound systems. In a similar way that there is a unique relation between the one-body density and the external potential that gives rise to it, we demonstrate that there is a unique relation between that particular many-body density and a definite many-body potential. The energy is then a functional of this density and its minimization leads to the ground-state energy of the system. As a proof of principle, the analogous of the Kohn-Sham equation is solved in the specific case of $^4$He atomic clusters, to put in evidence the advantages of this new formulation in terms of physical insights. △ Less

Submitted 24 June, 2021; originally announced June 2021.

Comments: 7 pages and 2 figures

Journal ref: Phys. Rev. A 104, 030801 (2021)

The X17 boson and the $^3$H$(p,e^+ e^-)^4$He and $^3$He$(n,e^+ e^-)^4$He processes: a theoretical analysis

Authors: M. Viviani, E. Filandri, L. Girlanda, C. Gustavino, A. Kievsky, L. E. Marcucci, R. Schiavilla

Abstract: The present work deals with $e^+$-$e^-$ pair production in the four-nucleon system. We first analyze the process as a purely electromagnetic one in the context of a state-of-the-art approach to nuclear strong-interaction dynamics and nuclear electromagnetic currents, derived from chiral effective field theory ($χ$EFT). Next, we examine how the exchange of a hypothetical low-mass boson would impact… ▽ More The present work deals with $e^+$-$e^-$ pair production in the four-nucleon system. We first analyze the process as a purely electromagnetic one in the context of a state-of-the-art approach to nuclear strong-interaction dynamics and nuclear electromagnetic currents, derived from chiral effective field theory ($χ$EFT). Next, we examine how the exchange of a hypothetical low-mass boson would impact the cross section for such a process. We consider several possibilities, that this boson is either a scalar, pseudoscalar, vector, or axial particle. The ab initio calculations use exact hyperspherical-harmonics methods to describe the bound state and low-energy spectrum of the $A\,$=$\,4$ continuum, and fully account for initial state interaction effects in the $3+1$ clusters. While electromagnetic interactions are treated to high orders in the chiral expansion, the interactions of the hypothetical boson with nucleons are modeled in leading-order $χ$EFT (albeit, in some instances, selected subleading contributions are also accounted for). We also provide an overview of possible future experiments probing pair production in the $A\,$=$\,4$ system at a number of candidate facilities.facilities. △ Less

Submitted 18 June, 2021; v1 submitted 15 April, 2021; originally announced April 2021.

Comments: 32 pages, 17 figures. Corrected version, having fixed a problem with a phase in the multipolar expansions. Only Figs. 4 and 10 change significantly, the rest only marginally. The conclusions are unchanged. We include also a number of improvements: 1) better wave functions, 2) the isoscalar pseudoscalar case, and 3) the comparison with the new ATOMKI data (see new Fig. 17)

Efimov Physics and Connections to Nuclear Physics

Authors: A. Kievsky, L. Girlanda, M. Gattobigio, M. Viviani

Abstract: Physical systems characterized by a shallow two-body bound or virtual state are governed at large distances by a continuous-scale invariance, which is broken to a discrete one when three or more particles come into play. This symmetry induces a universal behavior for different systems, independent of the details of the underlying interaction, rooted in the smallness of the ratio $\ell/a_B \ll 1$,… ▽ More Physical systems characterized by a shallow two-body bound or virtual state are governed at large distances by a continuous-scale invariance, which is broken to a discrete one when three or more particles come into play. This symmetry induces a universal behavior for different systems, independent of the details of the underlying interaction, rooted in the smallness of the ratio $\ell/a_B \ll 1$, where the length $a_B$ is associated to the binding energy of the two-body system $E_2=\hbar^2/m a_B^2$ and $\ell$ is the natural length given by the interaction range. Efimov physics refers to this universal behavior, which is often hidden by the on-set of system-specific non-universal effects. In this work we identify universal properties by providing an explicit link of physical systems to their unitary limit, in which $a_B\rightarrow\infty$, and show that nuclear systems belong to this class of universality. △ Less

Submitted 26 February, 2021; originally announced February 2021.

Comments: 27 pages, 8 figures

Two- and three-nucleon contact interactions and ground-state energies of light- and medium-mass nuclei

Authors: R. Schiavilla, L. Girlanda, A. Gnech, A. Kievsky, A. Lovato, L. E. Marcucci, M. Piarulli, M. Viviani

Abstract: Classes of two-nucleon ($2N$) contact interactions are developed in configuration space at leading order (LO), next-to-leading order (NLO), and next-to-next-to-next-to-leading order (N3LO) by fitting the experimental singlet $np$ scattering length and deuteron binding energy at LO, and $np$ and $pp$ scattering data in the laboratory-energy range 0--15 MeV at NLO and 0--25 MeV at N3LO. These intera… ▽ More Classes of two-nucleon ($2N$) contact interactions are developed in configuration space at leading order (LO), next-to-leading order (NLO), and next-to-next-to-next-to-leading order (N3LO) by fitting the experimental singlet $np$ scattering length and deuteron binding energy at LO, and $np$ and $pp$ scattering data in the laboratory-energy range 0--15 MeV at NLO and 0--25 MeV at N3LO. These interactions are regularized by including two Gaussian cutoffs, one for $T\,$=$\,0$ and the other for $T\,$=$\,1$ channels. The cutoffs are taken to vary in the ranges $R_0\,$=$(1.5$--2.3) fm and $R_1\,$=$(1.5$--3.0) fm. The 780 (1,100) data points up to 15 (25) MeV energy, primarily differential cross sections, are fitted by the NLO (N3LO) models with a $χ^2$/datum about 1.7 or less (well below 1.5), when harder cutoff values are adopted. As a first application, we report results for the binding energies of nuclei with mass numbers $A\,$=$\,3$--6 and 16 obtained with selected LO and NLO $2N$ models both by themselves as well as in combination with a LO three-nucleon ($3N$) contact interaction. The latter is characterized by a single low-energy constant that is fixed to reproduce the experimental $^3$H binding energy. The inclusion of the $3N$ interaction largely removes the sensitivity to cutoff variations in the few-nucleon systems and leads to predictions for the $^3$He and $^4$He binding energies that cluster around 7.8 MeV and 30 MeV, respectively. However, in $^{16}$O this cutoff sensitivity remains rather strong. Finally, predictions at LO only are also reported for medium-mass nuclei with $A\,$=$\,40$, 48, and 90. △ Less

Submitted 3 February, 2021; originally announced February 2021.

Comments: 23 pages, 6 figures

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

A Comprehensive Study of the Three- and Four-Neutron Systems at Low Energies

Authors: Michael D. Higgins, Chris H. Greene, Alejandro Kievsky, Michele Viviani

Abstract: This work presents further analysis of the three- and four-neutron systems in the low energy regime using adiabatic hyperspherical methods. In our previous Phys. Rev. Lett. article (Phys. Rev. Lett. 125, 052501 (2020)), the low-energy behavior of these neutron systems was treated in the adiabatic approximation, neglecting the off-diagonal non-adiabatic couplings. A thorough analysis of the density… ▽ More This work presents further analysis of the three- and four-neutron systems in the low energy regime using adiabatic hyperspherical methods. In our previous Phys. Rev. Lett. article (Phys. Rev. Lett. 125, 052501 (2020)), the low-energy behavior of these neutron systems was treated in the adiabatic approximation, neglecting the off-diagonal non-adiabatic couplings. A thorough analysis of the density of states through a multi-channel treatment of the three-and four-neutron scattering near the scattering continuum threshold is performed, showing no evidence of a 4n resonance at low energy. A detailed analysis of the long-range behavior of the lowest few adiabatic hyperspherical potentials shows there is an attractive $ρ^{-3}$ universal behavior which dominates in the low-energy regime of the multi-channel scattering. This long-range behavior leads to a divergent behavior of the density of state for $E\rightarrow0$ that could account for the low-energy signal observed in the 2016 experiment by Kisamori et al. (Phys. Rev. Lett. 116, 052501 (2016)). △ Less

Submitted 23 November, 2020; originally announced November 2020.

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

Gaussian characterization of the unitary window for $N=3$: bound, scattering and virtual states

Authors: A. Deltuva. A. Kievsky, M. Gattobigio, M. Viviani

Abstract: The three-body system inside the unitary window is studied for three equal bosons and three equal fermions having $1/2$ spin-isospin symmetry. We perform a gaussian characterization of the window using a gaussian potential to define trajectories for low-energy quantities as binding energies and phase shifts. On top of this trajectories experimental values are placed or, when not available, quantit… ▽ More The three-body system inside the unitary window is studied for three equal bosons and three equal fermions having $1/2$ spin-isospin symmetry. We perform a gaussian characterization of the window using a gaussian potential to define trajectories for low-energy quantities as binding energies and phase shifts. On top of this trajectories experimental values are placed or, when not available, quantities calculated using realistic potentials that are known to reproduce experimental values. The intention is to show that the gaussian characterization of the window, thought as a contact interaction plus range corrections, captures the main low-energy properties of real systems as for example three helium atoms or three nucleons. The mapping of real systems on the gaussian trajectories is taken as indication of universal behavior. The trajectories continuously link the physical points to the unitary limit allowing for the explanation of strong correlations between observables appearing in real systems and which are known to exist in that limit. In the present study we focus on low-energy bound, scattering and virtual states. △ Less

Submitted 13 November, 2020; originally announced November 2020.

Comments: 10 pages, 7 figures

Unitary ambiguity of NN contact interactions and the 3N force

Authors: L. Girlanda, A. Kievsky, L. E. Marcucci, M. Viviani

Abstract: We identify a redundancy between two- and three-nucleon contact interactions at the fourth and fifth order of the chiral expansion respectively. In particular we show that tensor-type and spin-orbit three-nucleon contact interactions effectively account for that part of the two-nucleon interaction which depends on the total center-of-mass momentum and is unconstrained by relativity. This might giv… ▽ More We identify a redundancy between two- and three-nucleon contact interactions at the fourth and fifth order of the chiral expansion respectively. In particular we show that tensor-type and spin-orbit three-nucleon contact interactions effectively account for that part of the two-nucleon interaction which depends on the total center-of-mass momentum and is unconstrained by relativity. This might give the chiral effective field theory enough flexibility to successfully address $A=3$ scattering observables already at N3LO. △ Less

Submitted 10 November, 2020; v1 submitted 8 July, 2020; originally announced July 2020.

Comments: 9 pages, discussion extended to OPE 3NF

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

From few to many bosons inside the unitary window: a transition between universal to non-universal behavior

Authors: A. Kievsky, A. Polls, B. Juliá-Díaz, N. K. Timofeyuk, M. Gattobigio

Abstract: Universal behaviour in few-bosons systems close to the unitary limit, where two bosons become unbound, has been intensively investigated in recent years both experimentally and theoretically. In this particular region, called the unitary window, details of the inter-particle interactions are not important and observables, such as binding energies, can be characterized by a few parameters. With an… ▽ More Universal behaviour in few-bosons systems close to the unitary limit, where two bosons become unbound, has been intensively investigated in recent years both experimentally and theoretically. In this particular region, called the unitary window, details of the inter-particle interactions are not important and observables, such as binding energies, can be characterized by a few parameters. With an increasing number of particles the short-range repulsion, present in all atomic, molecular or nuclear interactions, gradually induces deviations from the universal behaviour. In the present letter we discuss for the first time a simple way of incorporating non-universal behaviour through one specific parameter which controls the smooth transition of the system from universal to non-universal regime. Using a system of $N$ helium atoms as an example we calculate their ground state energies as trajectories within the unitary window and also show that the control parameters can be used to determine the energy per particle in homogeneous systems when $N \rightarrow \infty$. △ Less

Submitted 17 June, 2020; originally announced June 2020.

Comments: 6 pages, 3 figures

Journal ref: Phys. Rev. A 102, 063320 (2020)

Non-resonant Density of States Enhancement at Low Energies for Three or Four Neutrons

Authors: Michael D. Higgins, Chris H. Greene, Alejandro Kievsky, Michele Viviani

Abstract: The low energy systems of three or four neutrons are treated within the adiabatic hyperspherical framework, yielding an understanding of the low energy quantum states in terms of an adiabatic potential energy curve. The dominant low energy potential curve for each system, computed here using widely accepted nucleon-nucleon interactions with and without the inclusion of a three-nucleon force, shows… ▽ More The low energy systems of three or four neutrons are treated within the adiabatic hyperspherical framework, yielding an understanding of the low energy quantum states in terms of an adiabatic potential energy curve. The dominant low energy potential curve for each system, computed here using widely accepted nucleon-nucleon interactions with and without the inclusion of a three-nucleon force, shows no sign of a low energy resonance. However, both systems exhibit a low energy enhancement of the density of states, or of the Wigner-Smith time-delay, which derives from long-range universal physics analogous to the Efimov effect. That enhancement could be relevant to understanding the low energy excess of correlated 4-neutron ejection events observed experimentally in a nuclear reaction by Kisamori et al. △ Less

Submitted 26 June, 2020; v1 submitted 10 May, 2020; originally announced May 2020.

Journal ref: Phys. Rev. Lett. 125, 052501 (2020)

Study of $n+{}^3$H, $p+{}^3$He, $p+{}^3$H, and $n+{}^3$He scattering with the HH method

Authors: M. Viviani, L. Girlanda, A. Kievsky, L. E. Marcucci

Abstract: The $n+{}^3$H, $p+{}^3$He, $p+{}^3$H, and $n+{}^3$He elastic and charge exchange reactions at low energies are studied by means of the hyperspherical harmonic method. The considered nuclear Hamiltonians include modern two- and three-nucleon interactions, in particular results are reported in case of chiral two-nucleon potentials, with and without the inclusion of chiral three-nucleon (3N) interact… ▽ More The $n+{}^3$H, $p+{}^3$He, $p+{}^3$H, and $n+{}^3$He elastic and charge exchange reactions at low energies are studied by means of the hyperspherical harmonic method. The considered nuclear Hamiltonians include modern two- and three-nucleon interactions, in particular results are reported in case of chiral two-nucleon potentials, with and without the inclusion of chiral three-nucleon (3N) interactions. A detailed study of the convergence and numerical stability of the method is presented. We have found that the effect the 3N force is in general tiny except for $p+{}^3$H scattering below the opening of the $n+{}^3$He channel. In such a case, the effect of 3N forces is appreciable and a clear dependence on the cutoff used to regularize the high-momentum tail of the interactions is observed. Such a dependence is related to the presence of the poorly known sharp $0^+$ resonance, considered to be the first excited state of ${}^4$He. △ Less

Submitted 31 March, 2020; originally announced March 2020.

Comments: 33 pages, 15 figures

The Hyperspherical Harmonics method: a tool for testing and improving nuclear interaction models

Authors: L. E. Marcucci, J. Dohet-Eraly, L. Girlanda, A. Gnech, A. Kievsky, M. Viviani

Abstract: The Hyperspherical Harmonics (HH) method is one of the most accurate techniques to solve the quantum mechanical problem for nuclear systems with $A\le 4$. In particular, by applying the Rayleigh-Ritz or Kohn variational principle, both bound and scattering states can be addressed, using either local or non-local interactions. Thanks to this versatility, the method can be used to test the two- and… ▽ More The Hyperspherical Harmonics (HH) method is one of the most accurate techniques to solve the quantum mechanical problem for nuclear systems with $A\le 4$. In particular, by applying the Rayleigh-Ritz or Kohn variational principle, both bound and scattering states can be addressed, using either local or non-local interactions. Thanks to this versatility, the method can be used to test the two- and three-nucleon components of the nuclear interaction. In the present review we introduce the formalism of the HH method, both for bound and scattering states. In particular, we describe the implementation of the method to study the $A=3$ and $4$ scattering problem. Second, we present a selected choice of results of the last decade, most representative of the latest achievements. Finally, we conclude with a discussion of what we believe will be the most significant developments within the HH method for the next five-to-ten years. △ Less

Submitted 20 December, 2019; originally announced December 2019.

Comments: 25 pages, 9 figures, 4 tables, submitted to Frontiers in Physics (Research Topic) "The Long-Lasting Quest for Nuclear Interactions: The Past, the Present and the Future"

Embedding nuclear physics inside the unitary window

Authors: Mario Gattobigio, Alejandro Kievsky, Michele Viviani

Abstract: The large values of the singlet and triplet scattering lengths locate the two-nucleon system close to the unitary limit, the limit in which these two values diverge. As a consequence, the system shows a continuous scale invariance which strongly constrains the values of the observables, a well-known fact already noticed a long time ago. The three-nucleon system shows a discrete scale invariance th… ▽ More The large values of the singlet and triplet scattering lengths locate the two-nucleon system close to the unitary limit, the limit in which these two values diverge. As a consequence, the system shows a continuous scale invariance which strongly constrains the values of the observables, a well-known fact already noticed a long time ago. The three-nucleon system shows a discrete scale invariance that can be observed by correlations of the triton binding energy with other observables as the doublet nucleon-deuteron scattering length or the alpha-particle binding energy. The low-energy dynamics of these systems is universal; it does not depend on the details of the particular way in which the nucleons interact. Instead, it depends on a few control parameters, the large values of the scattering lengths and the triton binding energy. Using a potential model with variable strength set to give values to the control parameters, we study the spectrum of $A=2,3,4,6$ nuclei in the region between the unitary limit and their physical values. In particular, we analyze how the binding energies emerge from the unitary limit forming the observed levels. △ Less

Submitted 21 March, 2019; originally announced March 2019.

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

More on the universal equation for Efimov states

Authors: M. Gattobigio, M. Göbel, H. -W. Hammer, A. Kievsky

Abstract: Efimov states are a sequence of shallow three-body bound states that arise when the two-body scattering length is much larger than the range of the interaction. The binding energies of these states are described as a function of the scattering length and one three-body parameter by a transcendental equation involving a universal function of one angular variable. We provide an accurate and convenie… ▽ More Efimov states are a sequence of shallow three-body bound states that arise when the two-body scattering length is much larger than the range of the interaction. The binding energies of these states are described as a function of the scattering length and one three-body parameter by a transcendental equation involving a universal function of one angular variable. We provide an accurate and convenient parametrization of this function. Moreover, we discuss the effective treatment of range corrections in the universal equation and compare with a strictly perturbative scheme. △ Less

Submitted 13 March, 2019; originally announced March 2019.

Comments: 13 pages, 5 figures

Journal ref: Few-Body Syst (2019) 60: 40

Short-range three-nucleon interaction from A=3 data and its hierarchical structure

Authors: L. Girlanda, A. Kievsky, M. Viviani, L. E. Marcucci

Abstract: We construct accurate models of three-nucleon (3N) interaction by fitting, in a hybrid phenomenological approach, the low-energy constants parametrizing the subleading 3N contact operators to the triton binding energy, n-d scattering lengths, cross section and polarization observables of p-d scattering at 2 MeV center-of-mass energy. These models lead to a satisfactory description of polarized p-d… ▽ More We construct accurate models of three-nucleon (3N) interaction by fitting, in a hybrid phenomenological approach, the low-energy constants parametrizing the subleading 3N contact operators to the triton binding energy, n-d scattering lengths, cross section and polarization observables of p-d scattering at 2 MeV center-of-mass energy. These models lead to a satisfactory description of polarized p-d scattering data in the whole energy range below the deuteron breakup threshold. In particular, the long-standing $A_y$ puzzle seems to be solved thanks to the new terms considered in the 3N force. Two types of hierarchies among the subleading contact operators are also derived, based on the large-$N_c$ counting and on a recently proposed relativistic counting. We test these hierarchies against the same experimental data and show that they are respected at a reasonable level. △ Less

Submitted 23 November, 2018; originally announced November 2018.

Comments: 35 pages, 9 figures

Journal ref: Phys. Rev. C 99, 054003 (2019)

Local chiral interactions and magnetic structure of few-nucleon systems

Authors: R. Schiavilla, A. Baroni, S. Pastore, M. Piarulli, L. Girlanda, A. Kievsky, A. Lovato, L. E. Marcucci, Steven C. Pieper, M. Viviani, R. B. Wiringa

Abstract: The magnetic form factors of $^2$H, $^3$H, and $^3$He, deuteron photodisintegration cross sections at low energies, and deuteron threshold electrodisintegration cross sections at backward angles in a wide range of momentum transfers, are calculated with the chiral two-nucleon (and three-nucleon) interactions including $Δ$ intermediate states that have recently been constructed in configuration spa… ▽ More The magnetic form factors of $^2$H, $^3$H, and $^3$He, deuteron photodisintegration cross sections at low energies, and deuteron threshold electrodisintegration cross sections at backward angles in a wide range of momentum transfers, are calculated with the chiral two-nucleon (and three-nucleon) interactions including $Δ$ intermediate states that have recently been constructed in configuration space. The $A\,$=$\,$3 wave functions are obtained from hyperspherical-harmonics solutions of the Schrödinger equation. The electromagnetic current includes one- and two-body terms, the latter induced by one- and two-pion exchange (OPE and TPE, respectively) mechanisms and contact interactions. The contributions associated with $Δ$ intermediate states are only retained at the OPE level, and are neglected in TPE loop (tree-level) corrections to two-body (three-body) current operators. Expressions for these currents are derived and regularized in configuration space for consistency with the interactions. The low-energy constants that enter the contact few-nucleon systems. The predicted form factors and deuteron electrodisintegration cross section are in excellent agreement with experiment for momentum transfers up to 2--3 fm$^{-1}$. However, the experimental values for the deuteron photodisintegration cross section are consistently underestimated by theory, unless use is made of the Siegert form of the electric dipole transition operator. A complete analysis of the results is provided, including the clarification of the origin of the aforementioned discrepancy. △ Less

Submitted 26 September, 2018; originally announced September 2018.

Comments: 24 pages, 13 figures

Journal ref: Phys. Rev. C 99, 034005 (2019)

Local chiral interactions, the tritium Gamow-Teller matrix element, and the three-nucleon contact term

Authors: A. Baroni, R. Schiavilla, L. E. Marcucci, L. Girlanda, A. Kievsky, A. Lovato, S. Pastore, M. Piarulli, Steven C. Pieper, M. Viviani, R. B. Wiringa

Abstract: The Gamow-Teller (GT) matrix element contributing to tritium $β$ decay is calculated with trinucleon wave functions obtained from hyperspherical-harmonics solutions of the Schrödinger equation with the chiral two- and three-nucleon interactions including $Δ$ intermediate states that have recently been constructed in configuration space. Predictions up to N3LO in the chiral expansion of the axial c… ▽ More The Gamow-Teller (GT) matrix element contributing to tritium $β$ decay is calculated with trinucleon wave functions obtained from hyperspherical-harmonics solutions of the Schrödinger equation with the chiral two- and three-nucleon interactions including $Δ$ intermediate states that have recently been constructed in configuration space. Predictions up to N3LO in the chiral expansion of the axial current (with $Δ$'s) overestimate the empirical value by 1--4 \%. By exploiting the relation between the low-energy constant (LEC) in the contact three-nucleon interaction and two-body axial current, we provide new determinations of the LECs $c_D$ and $c_E$ that characterize this interaction by fitting the trinucleon binding energy and tritium GT matrix element. Some of the implications that the resulting models of three-nucleon interactions have on the spectra of light nuclei and the equation of state of neutron matter are briefly discussed. We also provide a partial analysis, which ignores $Δ$'s, of the contributions due to loop corrections in the axial current at N4LO. Finally, explicit expressions for the axial current up to N4LO have been derived in configuration space, which other researchers in the field may find useful. △ Less

Submitted 26 June, 2018; originally announced June 2018.

Comments: 16 pages, 3 figures

Journal ref: Phys. Rev. C 98, 044003 (2018)

Correlations imposed by the unitary limit between few-nucleon systems, nuclear matter and neutron stars

Authors: A. Kievsky, M. Viviani, D. Logoteta, I. Bombaci, L. Girlanda

Abstract: The large values of the singlet and triplet two-nucleon scattering lengths locate the nuclear system close to the unitary limit. This particular position strongly constrains the low-energy observables in the three-nucleon system as depending on one parameter, the triton binding energy, and introduces correlations in the low energy sector of light nuclei. Here we analyze the propagation of these co… ▽ More The large values of the singlet and triplet two-nucleon scattering lengths locate the nuclear system close to the unitary limit. This particular position strongly constrains the low-energy observables in the three-nucleon system as depending on one parameter, the triton binding energy, and introduces correlations in the low energy sector of light nuclei. Here we analyze the propagation of these correlations to infinite nuclear matter showing that its saturation properties, the equation of state of $β$-stable nuclear matter and several properties of neutron stars, as their maximum mass, are well determined solely by a few number of low-energy quantities of the two- and three-nucleon systems. In this way we make a direct link between the universal behavior observed in the low-energy region of few-nucleon systems and fundamental properties of nuclear matter and neutron stars. △ Less

Submitted 31 July, 2018; v1 submitted 7 June, 2018; originally announced June 2018.

Comments: 5 pages, 3 figures

Journal ref: Phys. Rev. Lett. 121, 072701 (2018)

Saturation properties of helium drops from a Leading Order description

Authors: A. Kievsky, A. Polls, B. Juliá-Díaz, N. K. Timofeyuk

Abstract: Saturation properties are directly linked to the short-range scale of the two-body interaction of the particles. The case of helium is particular, from one hand the two-body potential has a strong repulsion at short distances. On the other hand, the extremely weak binding of the helium dimer locates this system very close to the unitary limit allowing for a description based on an effective theory… ▽ More Saturation properties are directly linked to the short-range scale of the two-body interaction of the particles. The case of helium is particular, from one hand the two-body potential has a strong repulsion at short distances. On the other hand, the extremely weak binding of the helium dimer locates this system very close to the unitary limit allowing for a description based on an effective theory. At leading order of this theory a two- and a three-body term appear, each one characterized by a low energy constant. In a potential model this description corresponds to a soft potential model with a two-body term purely attractive plus a three-body term purely repulsive constructed to describe the dimer and trimer binding energies. Here we analyse the capability of this model to describe the saturation properties making a direct link between the low energy scale and the short-range correlations. We will show that the energy per particle, $E_N/N$, can be obtained with reasonable accuracy at leading order extending the validity of this approximation, characterizing universal behavior in few-boson systems close to the unitary limit, to the many-body system. △ Less

Submitted 19 July, 2017; v1 submitted 18 July, 2017; originally announced July 2017.

Comments: 5 pages, 3 figures

Journal ref: Phys. Rev. A 96, 040501 (2017)

Light-nuclei spectra from chiral dynamics

Authors: M. Piarulli, A. Baroni, L. Girlanda, A. Kievsky, A. Lovato, Ewing Lusk, L. E. Marcucci, Steven C. Pieper, R. Schiavilla, M. Viviani, R. B. Wiringa

Abstract: A major goal of nuclear theory is to explain the spectra and stability of nuclei in terms of effective many-body interactions amongst the nucleus' constituents-the nucleons, i.e., protons and neutrons. Such an approach, referred to below as the basic model of nuclear theory, is formulated in terms of point-like nucleons, which emerge as effective degrees of freedom, at sufficiently low energy, as… ▽ More A major goal of nuclear theory is to explain the spectra and stability of nuclei in terms of effective many-body interactions amongst the nucleus' constituents-the nucleons, i.e., protons and neutrons. Such an approach, referred to below as the basic model of nuclear theory, is formulated in terms of point-like nucleons, which emerge as effective degrees of freedom, at sufficiently low energy, as a result of a decimation process, starting from the fundamental quarks and gluons, described by Quantum Chromodynamics (QCD). A systematic way to account for the constraints imposed by the symmetries of QCD, in particular chiral symmetry, is provided by chiral effective field theory, in the framework of a low-energy expansion. Here we show, in quantum Monte Carlo calculations accurate to $\leq\!2\%$ of the binding energy, that two- and three-body chiral interactions fitted {\sl only} to bound- and scattering-state observables in, respectively, the two- and three-nucleon sectors, lead to predictions for the energy levels and level ordering of nuclei in the mass range $A\,$=$\,$4-12 in very satisfactory agreement with experimental data. Our findings provide strong support for the fundamental assumptions of the basic model, and pave the way to its systematic application to the electroweak structure and response of these systems as well as to more complex nuclei. △ Less

Submitted 10 July, 2017; originally announced July 2017.

Comments: 21 pages, 3 figures

Journal ref: Phys. Rev. Lett. 120, 052503 (2018)

Universal behavior of few-boson systems using potential models

Authors: A. Kievsky, M. Viviani, R. Álvarez-Rodrí guez, M. Gattobigio, A. Deltuva

Abstract: The universal behavior of a three-boson system close to the unitary limit is encoded in a simple dependence of many observables in terms of few parameters. For example the product of the three-body parameter $κ_*$ and the two-body scattering length $a$, $κ_* a$ depends on the angle $ξ$ defined by $E_3/E_2=\tan^2ξ$. A similar dependence is observed in the ratio $a_{AD}/a$ with $a_{AD}$ the boson-di… ▽ More The universal behavior of a three-boson system close to the unitary limit is encoded in a simple dependence of many observables in terms of few parameters. For example the product of the three-body parameter $κ_*$ and the two-body scattering length $a$, $κ_* a$ depends on the angle $ξ$ defined by $E_3/E_2=\tan^2ξ$. A similar dependence is observed in the ratio $a_{AD}/a$ with $a_{AD}$ the boson-dimer scattering length. We use a two-parameter potential to determine this simple behavior and, as an application, to compute $a_{AD}$ for the case of three $^4$He atoms. △ Less

Submitted 21 February, 2017; originally announced February 2017.

Comments: presented at the 23rd European Conference on Few-Body Problems in Physics

Journal ref: Few-Body Syst (2017) 58: 66

Implications of Efimov physics for the description of three and four nucleons in chiral effective field theory

Authors: A. Kievsky, M. Viviani, M. Gattobigio, L. Girlanda

Abstract: In chiral effective field theory the leading order (LO) nucleon-nucleon potential includes two contact terms, in the two spin channels $S=0,1$, and the one-pion-exchange potential. When the pion degrees of freedom are integrated out, as in the pionless effective field theory, the LO potential includes two contact terms only. In the three-nucleon system, the pionless theory includes a three-nucleon… ▽ More In chiral effective field theory the leading order (LO) nucleon-nucleon potential includes two contact terms, in the two spin channels $S=0,1$, and the one-pion-exchange potential. When the pion degrees of freedom are integrated out, as in the pionless effective field theory, the LO potential includes two contact terms only. In the three-nucleon system, the pionless theory includes a three-nucleon contact term interaction at LO whereas the chiral effective theory does not. Accordingly arbitrary differences could be observed in the LO description of three- and four-nucleon binding energies. We analyze the two theories at LO and conclude that a three-nucleon contact term is necessary at this order in both theories. In turn this implies that subleading three-nucleon contact terms should be promoted to lower orders. Furthermore this analysis shows that one single low energy constant might be sufficient to explain the large values of the singlet and triplet scattering lengths. △ Less

Submitted 23 January, 2017; v1 submitted 31 October, 2016; originally announced October 2016.

Comments: 5 pages, 3 figures

Benchmark calculation of p-3H and n-3He scattering

Authors: M. Viviani, A. Deltuva, R. Lazauskas, A. C. Fonseca, A. Kievsky, L. E. Marcucci

Abstract: p-3H and n-3He scattering in the energy range above the n-3He but below the d-d thresholds is studied by solving the 4-nucleon problem with a realistic nucleon-nucleon interaction. Three different methods -- Alt, Grassberger and Sandhas, Hyperspherical Harmonics, and Faddeev-Yakubovsky -- have been employed and their results for both elastic and charge-exchange processes are compared. We observe a… ▽ More p-3H and n-3He scattering in the energy range above the n-3He but below the d-d thresholds is studied by solving the 4-nucleon problem with a realistic nucleon-nucleon interaction. Three different methods -- Alt, Grassberger and Sandhas, Hyperspherical Harmonics, and Faddeev-Yakubovsky -- have been employed and their results for both elastic and charge-exchange processes are compared. We observe a good agreement between the three different methods, thus the obtained results may serve as a benchmark. A comparison with the available experimental data is also reported and discussed. △ Less

Submitted 28 October, 2016; originally announced October 2016.

Comments: 13 pages, 6 figures. arXiv admin note: text overlap with arXiv:1109.3625

Journal ref: Phys. Rev. C 95, 034003 (2017)

Nuclear matter properties from local chiral interactions with $Δ$ isobar intermediate states

Authors: Domenico Logoteta, Ignazio Bombaci, Alejandro Kievsky

Abstract: Using two-nucleon and three-nucleon interactions derived in the framework of chiral perturbation theory (ChPT) with and without the explicit $Δ$ isobar contributions, we calculate the energy per particle of symmetric nuclear matter and pure neutron matter in the framework of the microscopic Brueckner-Hartree-Fock approach. In particular, we present for the first time nuclear matter calculations us… ▽ More Using two-nucleon and three-nucleon interactions derived in the framework of chiral perturbation theory (ChPT) with and without the explicit $Δ$ isobar contributions, we calculate the energy per particle of symmetric nuclear matter and pure neutron matter in the framework of the microscopic Brueckner-Hartree-Fock approach. In particular, we present for the first time nuclear matter calculations using the new fully local in coordinate-space two-nucleon interaction at the next-to-next-to-next-to-leading-order (N3LO) of ChPT with $Δ$ isobar intermediate states (N3LO$Δ$) recently developed by Piarulli et al. [arXiv:1606:06335]. We find that using this N3LO$Δ$ potential, supplemented with a local N2LO three-nucleon interaction with explicit $Δ$ isobar degrees of freedom, it is possible to obtain a satisfactory saturation point of symmetric nuclear matter. For this combination of two- and three-nucleon interactions we also calculate the nuclear symmetry energy and we compare our results with the empirical constraints on this quantity obtained using the excitation energies to isobaric analog states in nuclei and using experimental data on the neutron skin thickness of heavy nuclei, finding a very good agreement with these empirical constraints in all the considered nucleonic density range. In addition, we find that the explicit inclusion of $Δ$ isobars diminishes the strength of the three-nucleon interactions needed the get a good saturation point of symmetric nuclear matter. We also compare the results of our calculations with those obtained by other research groups using chiral nuclear interactions with different many-body methods, finding in many cases a very satisfactory agreement. △ Less

Submitted 2 September, 2016; originally announced September 2016.

Local chiral potentials and the structure of light nuclei

Authors: Maria Piarulli, Luca Girlanda, Rocco Schiavilla, Alejandro Kievsky, Alessandro Lovato, Laura E. Marcucci, Steven C. Pieper, Michele Viviani, Robert B. Wiringa

Abstract: We present fully local versions of the minimally non-local nucleon-nucleon potentials constructed in a previous paper [M.\ Piarulli {\it et al.}, Phys.\ Rev.\ C {\bf 91}, 024003 (2015)], and use them in hypersperical-harmonics and quantum Monte Carlo calculations of ground and excited states of $^3$H, $^3$He, $^4$He, $^6$He, and $^6$Li nuclei. The long-range part of these local potentials includes… ▽ More We present fully local versions of the minimally non-local nucleon-nucleon potentials constructed in a previous paper [M.\ Piarulli {\it et al.}, Phys.\ Rev.\ C {\bf 91}, 024003 (2015)], and use them in hypersperical-harmonics and quantum Monte Carlo calculations of ground and excited states of $^3$H, $^3$He, $^4$He, $^6$He, and $^6$Li nuclei. The long-range part of these local potentials includes one- and two-pion exchange contributions without and with $Δ$-isobars in the intermediate states up to order $Q^3$ ($Q$ denotes generically the low momentum scale) in the chiral expansion, while the short-range part consists of contact interactions up to order $Q^4$. The low-energy constants multiplying these contact interactions are fitted to the 2013 Granada database in two different ranges of laboratory energies, either 0--125 MeV or 0--200 MeV, and to the deuteron binding energy and $nn$ singlet scattering length. Fits to these data are performed for three models characterized by long- and short-range cutoffs, $R_{\rm L}$ and $R_{\rm S}$ respectively, ranging from $(R_{\rm L},R_{\rm S})=(1.2,0.8)$ fm down to $(0.8,0.6)$ fm. The long-range (short-range) cutoff regularizes the one- and two-pion exchange (contact) part of the potential. △ Less

Submitted 20 June, 2016; originally announced June 2016.

Comments: 29 pages, 3 figures

Journal ref: Phys. Rev. C 94, 054007 (2016)

Tritium $β$-decay in chiral effective field theory

Authors: A. Baroni, L. Girlanda, A. Kievsky, L. E. Marcucci, R. Schiavilla, M. Viviani

Abstract: We evaluate the Fermi and Gamow-Teller (GT) matrix elements in tritium $β$-decay by including in the charge-changing weak current the corrections up to one loop recently derived in nuclear chiral effective field theory ($χ$ EFT). The trinucleon wave functions are obtained from hyperspherical-harmonics solutions of the Schrodinger equation with two- and three-nucleon potentials corresponding to ei… ▽ More We evaluate the Fermi and Gamow-Teller (GT) matrix elements in tritium $β$-decay by including in the charge-changing weak current the corrections up to one loop recently derived in nuclear chiral effective field theory ($χ$ EFT). The trinucleon wave functions are obtained from hyperspherical-harmonics solutions of the Schrodinger equation with two- and three-nucleon potentials corresponding to either $χ$ EFT (the N3LO/N2LO combination) or meson-exchange phenomenology (the AV18/UIX combination). We find that contributions due to loop corrections in the axial current are, in relative terms, as large as (and in some cases, dominate) those from one-pion exchange, which nominally occur at lower order in the power counting. We also provide values for the low-energy constants multiplying the contact axial current and three-nucleon potential, required to reproduce the experimental GT matrix element and trinucleon binding energies in the N3LO/N2LO and AV18/UIX calculations. △ Less

Submitted 4 January, 2017; v1 submitted 5 May, 2016; originally announced May 2016.

Comments: 19 pages,6 figures, corrections to Text as suggested by Referee added; Erratum: 4 pages, 3 figures, corrections to Eq.(20), Tables I, II, III, Figures 4, 5, conclusions unchanged

Matching universal behavior with potential models

Authors: R. Álvarez-Rodríguez, A. Deltuva. M. Gattobigio, A. Kievsky

Abstract: Two-, three-, and four-boson systems are studied close to the unitary limit using potential models constructed to reproduce the minimal information given by the two-body scattering length $a$ and the two-body binding energy or virtual state energy $E_2$. The particular path used to reach the unitary limit is given by varying the potential strength. In this way the energy spectrum in the three- and… ▽ More Two-, three-, and four-boson systems are studied close to the unitary limit using potential models constructed to reproduce the minimal information given by the two-body scattering length $a$ and the two-body binding energy or virtual state energy $E_2$. The particular path used to reach the unitary limit is given by varying the potential strength. In this way the energy spectrum in the three- and four-boson systems is computed. The lowest energy states show finite-range effects absorbed in the construction of level functions that can be used to study real systems. Higher energy levels are free from finite-range effects, therefore the corresponding level functions tend to the zero-range universal function. Using this property a zero-range equation for the four-boson system is proposed and the four-boson universal function is computed. △ Less

Submitted 18 March, 2016; originally announced March 2016.

Comments: 8 pages, 4 figures

Journal ref: Phys. Rev. A 93, 062701 (2016)

Efimov physics with $1/2$ spin-isospin fermions

Authors: A. Kievsky, M. Gattobigio

Abstract: The structure of few-fermion systems having $1/2$ spin-isospin symmetry is studied using potential models. The strength and range of the two-body potentials are fixed to describe low energy observables in the angular momentum $L=0$ state and spin $S=0,1$ channels of the two-body system. Successively the strength of the potentials are varied in order to explore energy regions in which the two-body… ▽ More The structure of few-fermion systems having $1/2$ spin-isospin symmetry is studied using potential models. The strength and range of the two-body potentials are fixed to describe low energy observables in the angular momentum $L=0$ state and spin $S=0,1$ channels of the two-body system. Successively the strength of the potentials are varied in order to explore energy regions in which the two-body scattering lengths are close to the unitary limit. This study is motivated by the fact that in the nuclear system the singlet and triplet scattering lengths are both large with respect to the range of the interaction. Accordingly we expect evidence of universal behavior in the three- and four-nucleon systems that can be observed from the study of correlations between observables. In particular we concentrate in the behavior of the first excited state of the three-nucleon system as the system moves away from the unitary limit. We also analyze the dependence on the range of the three-body force of some low-energy observables in the three- and four-nucleon systems. △ Less

Submitted 30 November, 2015; originally announced November 2015.

Comments: 12 pages, submitted to the special issue weakly bound exotic nuclei on Few-Body Systems

Implication of the proton-deuteron radiative capture for Big Bang Nucleosynthesis

Authors: L. E. Marcucci, G. Mangano, A. Kievsky, M. Viviani

Abstract: The astrophysical $S$-factor for the radiative capture $d(p,γ)^3$He in the energy-range of interest for Big Bang Nucleosynthesis (BBN) is calculated using an {\it ab-initio} approach. The nuclear Hamiltonian retains both two- and three-nucleon interactions - the Argonne $v_{18}$ and the Urbana IX, respectively. Both one- and many-body contributions to the nuclear current operator are included. The… ▽ More The astrophysical $S$-factor for the radiative capture $d(p,γ)^3$He in the energy-range of interest for Big Bang Nucleosynthesis (BBN) is calculated using an {\it ab-initio} approach. The nuclear Hamiltonian retains both two- and three-nucleon interactions - the Argonne $v_{18}$ and the Urbana IX, respectively. Both one- and many-body contributions to the nuclear current operator are included. The former retain for the first time, besides the $1/m$ leading order contribution ($m$ is the nucleon mass), also the next-to-leading order term, proportional to $1/m^3$. The many-body currents are constructed in order to satisfy the current conservation relation with the adopted Hamiltonian model. The hyperspherical harmonics technique is applied to solve the $A=3$ bound and scattering states. A particular attention is used in this second case in order to obtain, in the energy range of BBN, an uncertainty on the astrophysical $S$-factor of the order or below $\sim$1 %. Then, in this energy range, the $S$-factor is found to be $\sim$10 % larger than the currently adopted values.Part of this increase (1-3 %) is due to the $1/m^3$ one-body operator, while the remaining is due to the new more accurate scattering wave functions. We have studied the implication of this new determination for the $d(p,γ)^3$He $S$-factor on deuterium primordial abundance. We find that the predicted theoretical value for $^2$H/H is in excellent agreement with its experimental determination, using the most recent determination of baryon density of Planck experiment, and with a standard number of relativistic degrees of freedom $N_{\rm eff}=3.046$ during primordial nucleosynthesis. △ Less

Submitted 27 October, 2015; originally announced October 2015.

Comments: 5 pages, 2 figures, submitted to Phys. Rev. Lett

Journal ref: Phys. Rev. Lett. 116, 102501 (2016)

Weakly bound states with spin-isospin symmetry

Authors: A. Kievsky, M. Gattobigio

Abstract: We discuss weakly bound states of a few-fermion system having spin-isospin symmetry. This corresponds to the nuclear physics case in which the singlet, $a_0$, and triplet, $a_1$, $n-p$ scattering lengths are large with respect to the range of the nuclear interaction. The ratio of the two is about $a_0/a_1\approx-4.31$. This value defines a plane in which $a_0$ and $a_1$ can be varied up to the uni… ▽ More We discuss weakly bound states of a few-fermion system having spin-isospin symmetry. This corresponds to the nuclear physics case in which the singlet, $a_0$, and triplet, $a_1$, $n-p$ scattering lengths are large with respect to the range of the nuclear interaction. The ratio of the two is about $a_0/a_1\approx-4.31$. This value defines a plane in which $a_0$ and $a_1$ can be varied up to the unitary limit, $1/a_0=0$ and $1/a_1=0$, maintaining its ratio fixed. Using a spin dependant potential model we estimate the three-nucleon binding energy along that plane. This analysis can be considered an extension of the Efimov plot for three bosons to the case of three $1/2$-spin-isospin fermions. △ Less

Submitted 7 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

Universal range corrections to the Efimov trimer for a class of paths to the unitary limit

Authors: A. Kievsky, M. Gattobigio

Abstract: Using potential models we analyze range corrections to the universal law dictated by the Efimov theory of three bosons. In the case of finite-range interactions we have observed that, at first order, it is necessary to supplement the theory with one finite-range parameter, $Γ_n^3$, for each specific $n$-level [Kievsky and Gattobigio, Phys. Rev. A {\bf 87}, 052719 (2013)]. The value of $Γ_n^3$ depe… ▽ More Using potential models we analyze range corrections to the universal law dictated by the Efimov theory of three bosons. In the case of finite-range interactions we have observed that, at first order, it is necessary to supplement the theory with one finite-range parameter, $Γ_n^3$, for each specific $n$-level [Kievsky and Gattobigio, Phys. Rev. A {\bf 87}, 052719 (2013)]. The value of $Γ_n^3$ depends on the way the potentials is changed to tune the scattering length toward the unitary limit. In this work we analyze a particular path in which the length $r_B=a-a_B$, measuring the difference between the two-body scattering length $a$ and the energy scattering length $a_B$, results almost constant. Analyzing systems with very different scales, as atomic or nuclear systems, we observe that the finite-range parameter remains almost constant along the path with a numerical value of $Γ_0^3\approx 0.87$ for the ground state level. This observation suggests the possibility of constructing a single universal function that incorporate finite-range effects for this class of paths. The result is used to estimate the three-body parameter $κ_*$ in the case of real atomic systems brought to the unitary limit thought a broad Feshbach resonances. Furthermore, we show that the finite-range parameter can be put in relation with the two-body contact $C_2$ at the unitary limit. △ Less

Submitted 2 December, 2015; v1 submitted 13 July, 2015; originally announced July 2015.

Comments: accepted Physical Review A

Journal ref: Phys. Rev. A 92, 062715 (2015)