Search

In this work, we present the elastic electron scattering longitudinal form factors of $A\le 10$ nuclei computed in a variational Monte Carlo approach. We employ the Norfolk family of local chiral interactions and a consistent electromagnetic charge operator. Our calculations are compared both to data and past theoretical evaluations. This work represents, to our knowledge, the first exact many-body calculation of longitudinal form factors in the $7\,\leq\,A\,\leq\,10$ mass range. Finally, we identify $^9$Be and $^{10}$B as candidate targets for renewed experimental interest, as they exhibit the potential to provide more stringent constraints on the theoretical models., Comment: 12 pages, 6 figures

We present Quantum Monte Carlo calculations of magnetic moments, form factors, and densities of $A\le 10$ nuclei within a chiral effective field theory approach. We use the Norfolk two- and three-body chiral potentials and their consistent electromagnetic one- and two-nucleon current operators. We find that two-body contributions to the magnetic moment can be large (up to $\sim33\%$ in $A=9$ systems). We study the model dependence of these observables and place particular emphasis on investigating their sensitivity to using different cutoffs to regulate the many-nucleon operators. Calculations of elastic magnetic form factors for $A\leq 10$ nuclei show excellent agreement with the data out to momentum transfers $q\approx 3$ fm$^{-1}$., Comment: 25 pages, 14 figures, 6 tables

We present Quantum Monte Carlo calculations of magnetic form factors in $A=6-10$ nuclei, based on Norfolk two- and three-nucleon interactions, and associated one- and two-body electromagnetic currents. Agreement with the available experimental data for $^6$Li, $^7$Li, $^9$Be and $^{10}$B up to values of momentum transfer $q\sim 3$ fm$^{-1}$ is achieved when two-nucleon currents are accounted for. We present a set of predictions for the magnetic form factors of $^7$Be, $^8$Li, $^9$Li, and $^9$C. In these systems, two-body currents account for $\sim40-60\%$ of the total magnetic strength. Measurements in any of these radioactive systems would provide valuable insights on the nuclear magnetic structure emerging from the underlying many-nucleon dynamics. A particularly interesting case is that of $^7$Be, as it would enable investigations of the magnetic structure of mirror nuclei., Comment: 5 pages, 3 figures, 1 table

The homogeneous Bethe-Salpeter equation (hBSE), describing a bound system in a genuinely relativistic quantum-field theory framework, was solved for the first time by using a D-Wave quantum annealer. After applying standard techniques of discretization, the hBSE, in ladder approximation, can be formally transformed in a generalized eigenvalue problem (GEVP), with two square matrices: one symmetric and the other non symmetric. The latter matrix poses the challenge of obtaining a suitable formal approach for investigating the non symmetric GEVP by means of a quantum annealer, i.e to recast it as a quadratic unconstrained binary optimization problem. A broad numerical analysis of the proposed algorithms, applied to matrices of dimension up to 64, was carried out by using both the proprietary simulated-anneaing package and the D-Wave Advantage 4.1 system. The numerical results very nicely compare with those obtained with standard classical algorithms, and also show interesting scalability features., Comment: 20 pages, 11 figures

We compute ground-state and dynamical properties of $^4$He and $^{16}$O nuclei using as input high-resolution, phenomenological nucleon-nucleon and three-nucleon forces that are local in coordinate space. The nuclear Schr\"odinger equation for both nuclei is accurately solved employing the auxiliary-field diffusion Monte Carlo approach. For the $^4$He nucleus, detailed benchmarks are carried out with the hyperspherical harmonics method. In addition to presenting results for the binding energies and radii, we also analyze the momentum distributions of these nuclei and their Euclidean response function corresponding to the isoscalar density transition. The latter quantity is particularly relevant for lepton-nucleus scattering experiments, as it paves the way to quantum Monte Carlo calculations of electroweak response functions of $^{16}$O., Comment: 14 pages, 9 figures

In pursuing the essential elements of nuclear binding, we compute ground-state properties of atomic nuclei with up to $A=20$ nucleons, using as input a leading order pionless effective field theory Hamiltonian. A variational Monte Carlo method based on a new, highly-expressive, neural-network quantum state ansatz is employed to solve the many-body Schr\"odinger equation in a systematically improvable fashion. In addition to binding energies and charge radii, we accurately evaluate the magnetic moments of these nuclei, as they reveal the self-emergence of the shell structure, which is not a priori encoded in the neural-network ansatz. To this aim, we introduce a novel computational protocol based on adding an external magnetic field to the nuclear Hamiltonian, which allows the neural network to learn the preferred polarization of the nucleus within the given magnetic field., Comment: 7 pages, 3 figures

We compute the muon capture on deuteron in the doublet hyperfine state for a variety of nuclear interactions and consistent nuclear currents. Our analysis includes a detailed examination of the theoretical uncertainties coming from different sources: the single-nucleon axial form factor, the truncation of the interaction and current chiral expansion, and the model dependence. Moreover, we study the impact of the use of different power counting scheme for the electroweak currents on the truncation error. To estimate the truncation error of the chiral expansion of interactions and currents we use the most modern techniques based on Bayesian analysis. This method enables us to give a clear statistical interpretation of the computed theoretical uncertainties. Finally, we provide the differential capture rate as function of the kinetic energy of the outgoing neutron which may be measured in future experiments. Our recommended theoretical value for the total doublet capture rate is $\Gamma_{\rm th}=395\pm 10$ s$^{-1}$ ($68\%$ confidence level). We calculated also the capture rate in the quartet hyperfine state, which turns out to be in the range $[13.3-13.8]$ s$^{-1}$ depending on the adopted nuclear interaction., Comment: 14 pages, 10 figures, Supplemental Materials available under request. Accepted version on PRC

It is a fascinating phenomenon in nuclear physics that states with a pronounced few-body structure can emerge from the complex dynamics of many nucleons. Such halo or cluster states often appear near the boundaries of nuclear stability. As such, they are an important part of the experimental program beginning at the Facility for Rare Isotope Beams (FRIB). A concerted effort of theory and experiment is necessary both to analyze experiments involving effective few-body states, as well as to constrain and refine theories of the nuclear force in light of new data from these experiments. As a contribution to exactly this effort, this paper compiles a collection of ``perspectives'' that emerged out of the Topical Program ``Few-body cluster structures in exotic nuclei and their role in FRIB experiments'' that was held at FRIB in August 2022 and brought together theorists and experimentalists working on this topic., Comment: This version should match the published one

The muon capture reaction $\mu^- + d \rightarrow n + n + \nu_{\mu}$ in the doublet hyperfine state is studied using nuclear potentials and consistent currents derived in chiral effective field theory, which are local and expressed in coordinate-space (the so-called Norfolk models). Only the largest contribution due to the $^1S_0$ $nn$ scattering state is considered. Particular attention is given to the estimate of the theoretical uncertainty, for which four sources have been identified: (i) the model dependence, (ii) the chiral order convergence for the weak nuclear current, (iii) the uncertainty in the single-nucleon axial form factor, and (iv) the numerical technique adopted to solve the bound and scattering $A=2$ systems. This last source of uncertainty has turned out essentially negligible. The $^1S_0$ doublet muon capture rate $\Gamma^D(^1S_0)$ has been found to be $\Gamma^D(^1S_0)=255.8(0.6)(4.4)(2.9)$ s$^{-1}$, where the three errors come from the first three sources of uncertainty. The value for $\Gamma^D(^1S_0)$ obtained within this local chiral framework is compared with previous calculations and found in very good agreement., Comment: 33 pages, 3 figures

The five low-energy constants (LECs) in the electromagnetic current derived in chiral effective field theory ($\chi$EFT) up to one loop are determined by a simultaneous fit to the $A\,$=$\,2$--3 nuclei magnetic moments and to the deuteron magnetic form factor and threshold electrodisintegration at backward angles over a wide range of momentum transfers. The resulting parametrization then yields predictions for the $^3$He/$^3$H magnetic form factors in excellent accord with the experimental values for momentum transfers ranging up to $\approx 0.8$ GeV/c, beyond the expected regime of validity of the $\chi$EFT approach. The calculations are based on last-generation two-nucleon interactions including high orders in the chiral expansion and derived by Entem, Macheleidt, and Nosyk [Phys.\ Rev.\ C {\bf 96}, 024004 (2017)] and by Piarulli {\it et al.} [Phys.\ Rev.\ C {\bf 94}, 054007 (2016)], using different $\chi$EFT formulations. In the $A\,$=$\,3$ calculations, (chiral) three-nucleon interactions are also accounted for. The model dependence resulting from these different formulations of the interactions is found to be mild for momentum transfer below $\approx0.8$ GeV/c. An analysis of the convergence of the chiral expansion is also provided., Comment: 12 pages, 8 figures

The ground-breaking works of Weinberg have opened the way to calculations of atomic nuclei that are based on systematically improvable Hamiltonians. Solving the associated many-body Schr\"odinger equation involves non-trivial difficulties, due to the non-perturbative nature and strong spin-isospin dependence of nuclear interactions. Artificial neural networks have proven to be able to compactly represent the wave functions of nuclei with up to $A=4$ nucleons. In this work, we extend this approach to $^6$Li and $^6$He nuclei, using as input a leading-order pionless effective field theory Hamiltonian. We successfully benchmark their binding energies, point-nucleon densities, and radii with the highly accurate hyperspherical harmonics method., Comment: 15 pages, 3 figures, invited contribution to the special issue in Few-Body Systems "Celebrating 30 years of Steven Weinberg's papers on Nuclear Forces from Chiral Lagrangians"

We report on a study of the Gamow-Teller matrix element contributing to ${}^6$He $\beta$-decay with similarity renormalization group (SRG) versions of momentum- and configuration-space two-nucleon interactions. These interactions are derived from two different formulations of chiral effective field theory ($\chi$EFT) -- without and with the explicit inclusion of $\Delta$-isobars. We consider evolution parameters $\Lambda_{\rm SRG}$ in the range between 1.2 and 2.0 fm$^{-1}$ and, for the $\Delta$-less case, also the unevolved (bare) interaction. The axial current contains one- and two-body terms, consistently derived at tree level (no loops) in the two distinct $\chi$EFT formulations we have adopted here. The ${}^6$He and ${}^6$Li ground-state wave functions are obtained from hyperspherical-harmonics (HH) solutions of the nuclear many-body problem. In $A\,$=$\,6$ systems, the HH method is limited at present to treat only two-body interactions and non-SRG evolved currents. Our results exhibit a significant dependence on $\Lambda_{\text{SRG}}$ of the contributions associated with two-body currents, suggesting that a consistent SRG-evolution of these is needed in order to obtain reliable estimates. We also show that the contributions from one-pion-exchange currents depend strongly on the model (chiral) interactions and on the momentum- or configuration-space cutoffs used to regularize them. These results might prove helpful in clarifying the origin of the sign difference recently found in No-Core-Shell-Model and Quantum Monte Carlo calculations of the ${}^6$He Gamow-Teller matrix element., Comment: 16 pages, 4 figures, accepted for publication on Phys. Rev. C

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 $\chi^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., Comment: 23 pages, 6 figures

Standard Big Bang Nucleosynthesis (BBN) predicts the abundances of the light elements in the early universe. Even if the overall agreement with the experimental data is good, still some discrepancies exist on the relic abundances of ${}^7$Li and ${}^6$Li. In order to exclude or confirm these scenarios, the BBN model needs precise input parameters, in particular the cross-sections of the BBN nuclear reaction network. However, the suppression of the cross-sections due to the Coulomb barrier makes the measurement very difficult and so affected by large systematic errors. Therefore, reliable theoretical calculations result fundamental in order to reduce the uncertainties. In this work we present a theoretical study of two nuclear reactions connected to ${}^6$Li abundance and recently the $\alpha$+d$\rightarrow$ ${}^6$Li + $\gamma$ and the p+${}^6$Li$\rightarrow$${}^7$Be+$\gamma$ radiative captures. For the first reaction we use a so-called ab-initio approach in which we solve the full six-body problem by using realistic nuclear potentials to describe the nucleon interactions. In particular we concentrate on the calculation and characterization of the final state of the reaction, the ${}^6$Li ground state, focusing on the electromagnetic static structure and the quantities relevant from the astrophysical point of view such as the asymptotic normalization coefficient. For doing this we use the Hyperspherical Harmonic approach developed by the Pisa group providing for the first time the possibility of using this approach beyond A = 4 nuclear systems. The second reaction is instead studied by using a two-body cluster approach where the proton and ${}^6$Li are considered as structureless particles. The angular distribution of the emitted photon obtained in this work were used by the LUNA Collaboration to determine the efficiency of the detector used in the measurement of the reaction., Comment: PhD Thesis

Electric dipole moments of nuclei, diamagnetic atoms, and certain molecules are induced by CP-violating nuclear forces. Naive dimensional analysis predicts these forces to be dominated by long-range one-pion-exchange processes, with short-range forces entering only at next-to-next-to-leading order in the chiral expansion. Based on renormalization arguments we argue that a consistent picture of CP-violating nuclear forces requires a leading-order short-distance operator contributing to ${}^1S_0$-${}^3P_0$ transitions, due to the attractive and singular nature of the strong tensor force in the ${}^3P_0$ channel. The short-distance operator leads to $\mathcal O(1)$ corrections to static and oscillating, relevant for axion searches, electric dipole moments. We discuss strategies how the finite part of the associated low-energy constant can be determined in the case of CP violation from the QCD theta term by the connection to charge-symmetry violation in nuclear systems.

We have studied the solution of the six-nucleon bound state problem using the hyperspherical harmonic (HH) approach. For this study we have considered only two-body nuclear forces. In particular we have used a chiral nuclear potential evolved with the similarity renormalization group unitary transformation. A restricted basis has been selected by performing a careful analysis of the convergence of different HH classes. Finally, the binding energy and other properties of ${}^6$Li ground state are calculated and compared with the results obtained by other techniques. Then, we present a calculation of matrix elements relevant for direct dark matter search involving ${}^6$Li. The results obtained demonstrate the feasibility of using the HH method to perform calculation beyond $A = 4$., Comment: 21 pages, 9 figures

Parity-violating and time-reversal conserving (PVTC) and parity-violating and time-reversal-violating (PVTV) forces in nuclei form only a tiny component of the total interaction between nucleons. The study of these tiny forces can nevertheless be of extreme interest because they allow to obtain information on fundamental symmetries using nuclear systems. The PVTC interaction derives from the weak interaction between the quarks inside nucleons and nuclei and the study of PVTC effects opens a window on the quark-quark weak interaction. The PVTV interaction is sensitive to more exotic interactions at the fundamental level, in particular to strong CP violation in the Standard Model Lagrangian, or even to exotic phenomena predicted in various beyond-the-Standard-Model scenarios. The presence of these interactions can be revealed either by studying various asymmetries in polarized scattering of nuclear systems, or by measuring the presence of non-vanishing permanent electric dipole moments of nucleons, nuclei and diamagnetic atoms and molecules. In this contribution, we review the derivation of the nuclear PVTC and PVTV interactions within various frameworks. We focus in particular on the application of chiral effective field theory, which allows for a more strict connection with the fundamental interactions at the quark level. We investigate PVTC and PVTV effects induced by these potential on several few-nucleon observables, such as the longitudinal asymmetry in proton-proton scattering and radiative neutron-proton capture, and the electric dipole momentsof the deuteron and the trinucleon system., Comment: 68 pages, contribution to Frontiers in Physics

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., Comment: 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"

Time Reversal Violation (TRV) interactions between quarks which appear in Standard Model (SM) and beyond-SM theories can induce TRV components in the nucleon-nucleon potential. The effects of these components can be studied by measuring the electric dipole moment (EDM) of light nuclei. In this work we present a complete derivation of the TRV nucleon-nucleon and three-nucleon potential up to next-to-next-to leading order (N2LO) in a chiral effective field theory (chiEFT) framework. The TRV interaction is then used to evaluate the EDM of 2H, 3H and 3he focusing in particular on the effects of the TRV three-body force and on the calculation of the theoretical errors. In case of a measurement of the EDM of these nuclei, the result of present work would permit to determine the values of the low energy constants and to identify the source of TRV., Comment: 19 pages, 4 figures

We present a new calculation of the ${}^6{\rm Li}(p,\gamma){}^7{\rm Be}$ radiative capture astrophysical $S$-factor in a cluster model framework. We consider several intercluster potentials, adjusted to reproduce the ${}^7{\rm Be}$ bound state properties and the $p-{}^6{\rm Li}$ elastic scattering phase shifts. Using these potentials, we calculate the astrophysical $S$-factor, obtaining a good agreement with available data, and the photon angular distribution. Finally, we discuss the consequences of a hypothetical resonance-like structure on the $S$-factor., Comment: 27 pages, 9 figures