Source author record

R. Schiavilla

R. Schiavilla appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

54works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

54 published item(s)

preprint2022arXiv

Partial muon capture rates in $A=3$ and $A=6$ nuclei with chiral effective field theory

Searches for neutrinoless-double beta decay rates are crucial in addressing questions within fundamental symmetries and neutrino physics. The rates of these decays depend not only on unknown parameters associated with neutrinos, but also on nuclear properties. In order to reliably extract information about the neutrino, one needs an accurate treatment of the complex many-body dynamics of the nucleus. Neutrinoless-double beta decays take place at momentum transfers on the order of 100 MeV/$c$ and require both nuclear electroweak vector and axial current matrix elements. Muon capture, a process in the same momentum transfer regime, has readily available experimental data to validate these currents. In this work, we present results of {\it ab initio} calculations of partial muon capture rates for $^3$He and $^6$Li nuclei using variational and Green's Function Monte Carlo computational methods. We estimate the impact of the three-nucleon interactions, the cutoffs used to regularize two-nucleon ($2N$) interactions, and the energy range of $2N$ scattering data used to fit these interactions.

preprint2021arXiv

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

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.

preprint2020arXiv

Chiral Effective Field Theory Calculations of Weak Transitions in Light Nuclei

We report Quantum Monte Carlo calculations of weak transitions in $A\leq 10$ nuclei, based on the Norfolk two- and three-nucleon chiral interactions, and associated one- and two-body axial currents. We find that the contribution from two-body currents is at the $2$ - $3\%$ level, with the exception of matrix elements entering the rates of $^8$Li, $^8$B, and $^8$He beta decays. These matrix elements are suppressed in impulse approximation based on the (leading order) Gamow Teller transition operator alone; two-body currents provide a $20$ - $30\%$ correction, which is, however, insufficient to bring theory in agreement with experimental data. For the other transitions, the agreement with the data is satisfactory, and the results exhibit a negligible to mild model dependence when different combinations of Norfolk interactions are utilized to construct the nuclear wave functions. We report a complete study of two-body weak transition densities which reveals the expected universal behavior of two-body currents at short distances throughout the range of $A\,$=$\,3$ to $A\,$=$\,10$ systems considered here.

preprint2016arXiv

Electromagnetic response of 12C: a first-principles calculation

The longitudinal and transverse electromagnetic response functions of $^{12}$C are computed in a "first-principles" Green's function Monte Carlo calculation, based on realistic two- and three-nucleon interactions and associated one- and two-body currents. We find excellent agreement between theory and experiment and, in particular, no evidence for the quenching of measured versus calculated longitudinal response. This is further corroborated by a re-analysis of the Coulomb sum rule, in which the contributions from the low-lying $J^π\,$=$\, 2^+$, $0^+_2$ (Hoyle), and $4^+$ states in $^{12}$C are accounted for explicitly in evaluating the total inelastic strength.

preprint2016arXiv

Nuclear Axial Currents in Chiral Effective Field Theory

Two-nucleon axial charge and current operators are derived in chiral effective field theory up to one loop. The derivation is based on time-ordered perturbation theory, and accounts for cancellations between the contributions of irreducible diagrams and the contributions due to non-static corrections from energy denominators of reducible diagrams. Ultraviolet divergencies associated with the loop corrections are isolated in dimensional regularization. The resulting axial current is finite and conserved in the chiral limit, while the axial charge requires renormalization. A complete set of contact terms for the axial charge up to the relevant order in the power counting is constructed.

preprint2015arXiv

Ab initio calculation of the electromagnetic and neutral-weak response functions of 4He and 12C

Precise measurement of neutrino oscillations, and hence the determination of their masses demands a quantitative understanding of neutrino-nucleus interactions. To this aim, two-body meson-exchange currents have to be accounted for along within realistic models of nuclear dynamics. We summarize our progresses towards the construction of a consistent framework, based on quantum Monte Carlo methods and on the spectral function approach, that can be exploited to accurately describe neutrino interactions with atomic nuclei over the broad kinematical region covered by neutrino experiments.

preprint2015arXiv

Electromagnetic and neutral-weak response functions of 4He and 12C

Ab initio calculations of the quasi-elastic electromagnetic and neutral-weak response functions of 4He and 12C are carried out for the first time. They are based on a realistic approach to nuclear dynamics, in which the strong interactions are described by two- and three-nucleon potentials and the electroweak interactions with external fields include one- and two-body terms. The Green's function Monte Carlo method is used to calculate directly the Laplace transforms of the response functions, and maximum-entropy techniques are employed to invert the resulting imaginary-time correlation functions with associated statistical errors. The theoretical results, confirmed by experiment in the electromagnetic case, show that two-body currents generate excess transverse strength from threshold to the quasi-elastic to the dip region and beyond. These findings challenge the conventional picture of quasi-elastic inclusive scattering as being largely dominated by single-nucleon knockout processes.

preprint2015arXiv

Electromagnetic Structure of Few-Nucleon Ground States

Experimental form factors of the hydrogen and helium isotopes, extracted from an up-to-date global analysis of cross sections and polarization observables measured in elastic electron scattering from these systems, are compared to predictions obtained in three different theoretical approaches: the first is based on realistic interactions and currents, including relativistic corrections (labeled as the conventional approach); the second relies on a chiral effective field theory description of the strong and electromagnetic interactions in nuclei (labeled $χ$EFT); the third utilizes a fully relativistic treatment of nuclear dynamics as implemented in the covariant spectator theory (labeled CST). For momentum transfers below $Q \lesssim 5$ fm$^{-1}$ there is satisfactory agreement between experimental data and theoretical results in all three approaches. However, at $Q \gtrsim 5$ fm$^{-1}$, particularly in the case of the deuteron, a relativistic treatment of the dynamics, as is done in the CST, is necessary. The experimental data on the deuteron $A$ structure function extend to $Q \simeq 12$ fm$^{-1}$, and the close agreement between these data and the CST results suggests that, even in this extreme kinematical regime, there is no evidence for new effects coming from quark and gluon degrees of freedom at short distances.

preprint2015arXiv

Measurement of the Target-Normal Single-Spin Asymmetry in Quasi-Elastic Scattering from the Reaction $^3$He$^\uparrow(e,e^\prime)$

We report the first measurement of the target single-spin asymmetry, $A_y$, in quasi-elastic scattering from the inclusive reaction $^3$He$^{\uparrow}(e,e^\prime)$ on a $^3$He gas target polarized normal to the lepton scattering plane. Assuming time-reversal invariance, this asymmetry is strictly zero for one-photon exchange. A non-zero $A_y$ can arise from the interference between the one- and two-photon exchange processes which is sensitive to the details of the sub-structure of the nucleon. An experiment recently completed at Jefferson Lab yielded asymmetries with high statistical precision at $Q^{2}=$ 0.13, 0.46 and 0.97 GeV$^{2}$. These measurements demonstrate, for the first time, that the $^3$He asymmetry is clearly non-zero and negative with a statistical significance of (8-10)$σ$. Using measured proton-to-$^{3}$He cross-section ratios and the effective polarization approximation, neutron asymmetries of $-$(1-3)% were obtained. The neutron asymmetry at high $Q^2$ is related to moments of the Generalized Parton Distributions (GPDs). Our measured neutron asymmetry at $Q^2=0.97$ GeV$^2$ agrees well with a prediction based on two-photon exchange using a GPD model and thus provides a new, independent constraint on these distributions.

preprint2015arXiv

Minimally non-local nucleon-nucleon potentials with chiral two-pion exchange including $Δ$'s

We construct a coordinate-space chiral potential, including $Δ$-isobar intermediate states in its two-pion-exchange component. The contact interactions entering at next-to-leading and next-to-next-to-next-to-leading orders ($Q^2$ and $Q^4$, respectively, $Q$ denoting generically the low momentum scale) are rearranged by Fierz transformations to yield terms at most quadratic in the relative momentum operator of the two nucleons. The low-energy constants multiplying these contact interactions are fitted to the 2013 Granada database, consisting of 2309 $pp$ and 2982 $np$ data (including, respectively, 148 and 218 normalizations) in the laboratory-energy range 0--300 MeV. For the total 5291 $pp$ and $np$ data in this range, we obtain a $χ^2$/datum of roughly 1.3 for a set of 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.

preprint2015arXiv

Neutrino-nucleus interactions and the short-range structure of nuclei

Improvements in theoretical modeling of Short Range structures and phenomena, and comparisons with data, will require sustained collaboration between nuclear theorists and neutrino experimentalists. The extensive history of studying this area of nuclear physics in electron- and hadron-scattering experiments, coupled with the transformative capabilities of LArTPCs to identify neutrinos, will provide a ripe opportunity for new discoveries that will further our understanding of the nucleus.

preprint2015arXiv

Quantum Monte Carlo methods for nuclear physics

Quantum Monte Carlo methods have proved very valuable to study the structure and reactions of light nuclei and nucleonic matter starting from realistic nuclear interactions and currents. These ab-initio calculations reproduce many low-lying states, moments and transitions in light nuclei, and simultaneously predict many properties of light nuclei and neutron matter over a rather wide range of energy and momenta. We review the nuclear interactions and currents, and describe the continuum Quantum Monte Carlo methods used in nuclear physics. These methods are similar to those used in condensed matter and electronic structure but naturally include spin-isospin, tensor, spin-orbit, and three-body interactions. We present a variety of results including the low-lying spectra of light nuclei, nuclear form factors, and transition matrix elements. We also describe low-energy scattering techniques, studies of the electroweak response of nuclei relevant in electron and neutrino scattering, and the properties of dense nucleonic matter as found in neutron stars. A coherent picture of nuclear structure and dynamics emerges based upon rather simple but realistic interactions and currents.

preprint2014arXiv

A chiral effective field theory study of hadronic parity violation in few-nucleon systems

We reconsider the derivation of the nucleon-nucleon parity-violating (PV) potential within a chiral effective field theory framework. We construct the potential up to next-to-next-to-leading order by including one-pion-exchange, two-pion-exchange, contact, and 1/M (M being the nucleon mass) terms, and use dimensional regularization to renormalize the pion-loop corrections. A detailed analysis of the number of independent low-energy constants (LEC's) entering the potential is carried out. We find that it depends on six LEC's: the pion-nucleon PV coupling constant $h^1_π$ and five parameters multiplying contact interactions. We investigate PV effects induced by this potential on several few-nucleon observables, including the $\vec{p}$-$p$ longitudinal asymmetry, the neutron spin rotation in $\vec{n}$-$p$ and $\vec{n}$-$d$ scattering, and the longitudinal asymmetry in the $^3$He$(\vec{n},p)^3$H charge-exchange reaction. An estimate for the range of values of the various LEC's is provided by using available experimental data.

preprint2014arXiv

Measurement of double-polarization asymmetries in the quasi-elastic $^3\vec{\mathrm{He}}(\vec{\mathrm{e}},\mathrm{e}'\mathrm{d})$ process

We present a precise measurement of double-polarization asymmetries in the $^3\vec{\mathrm{He}}(\vec{\mathrm{e}},\mathrm{e}'\mathrm{d})$ reaction. This particular process is a uniquely sensitive probe of hadron dynamics in $^3\mathrm{He}$ and the structure of the underlying electromagnetic currents. The measurements have been performed in and around quasi-elastic kinematics at $Q^2 = 0.25\,(\mathrm{GeV}/c)^2$ for missing momenta up to $270\,\mathrm{MeV}/c$. The asymmetries are in fair agreement with the state-of-the-art calculations in terms of their functional dependencies on $p_\mathrm{m}$ and $ω$, but are systematically offset. Beyond the region of the quasi-elastic peak, the discrepancies become even more pronounced. Thus, our measurements have been able to reveal deficiencies in the most sophisticated calculations of the three-body nuclear system, and indicate that further refinement in the treatment of their two- and/or three-body dynamics is required.

preprint2014arXiv

Neutral weak current two-body contributions in inclusive scattering from $^{12}$C

An ab initio calculation of the sum rules of the neutral weak response functions in $^{12}$C is reported, based on a realistic Hamiltonian, including two- and three-nucleon potentials, and on realistic currents, consisting of one- and two-body terms. We find that the sum rules of the response functions associated with the longitudinal and transverse components of the (space-like) neutral current are largest and that a significant portion ($\simeq 30$ %) of the calculated strength is due to two-body terms. This fact may have implications for the MiniBooNE and other neutrino quasi-elastic scattering data on nuclei.

preprint2014arXiv

Quantum Monte Carlo calculations of electromagnetic transitions in 8Be with meson-exchange currents derived from chiral effective field theory

We report quantum Monte Carlo calculations of electromagnetic transitions in 8Be. The realistic Argonne v18 two-nucleon and Illinois-7 three-nucleon potentials are used to generate the ground state and nine excited states, with energies that are in excellent agreement with experiment. A dozen M1 and eight E2 transition matrix elements between these states are then evaluated. The E2 matrix elements are computed only in impulse approximation, with those transitions from broad resonant states requiring special treatment. The M1 matrix elements include two-body meson-exchange currents derived from chiral effective field theory, which typically contribute 20--30% of the total expectation value. Many of the transitions are between isospin-mixed states; the calculations are performed for isospin-pure states and then combined with empirical mixing coefficients to compare to experiment. Alternate mixings are also explored. In general, we find that transitions between states that have the same dominant spatial symmetry are in reasonable agreement with experiment, but those transitions between different spatial symmetries are often underpredicted.

preprint2014arXiv

The $^3$He$(e,e^\prime p)^2$H and $^4$He$(e,e^\prime p)^3$H reactions at high momentum transfer

We present updated calculations for observables in the processes $^3$He($e,e^\prime p$)$^2$H, $^4$He($e,e^\prime p$)$^3$H, and $^4$He($\vec{e},e^\prime \vec{p}\,$)$^3$H. This update entails the implementation of improved nucleon-nucleon ($NN$) amplitudes to describe final state interactions (FSI) within a Glauber approximation and includes full spin-isospin dependence in the profile operator. In addition, an optical potential, which has also been updated since previous work, is utilized to treat FSI for the $^4$He($e,e^\prime p$)$^3$H and $^4$He($\vec{e},e^\prime \vec{p}\,$)$^3$H reactions. The calculations are compared with experimental data and show good agreement between theory and experiment. Comparisons are made between the various approximations in the Glauber treatment, including model dependence due to the $NN$ scattering amplitudes, rescattering contributions, and spin dependence. We also analyze the validity of the Glauber approximation at the kinematics the data is available, by comparing to the results obtained with the optical potential.

preprint2013arXiv

Charge form factor and sum rules of electromagnetic response functions in Carbon-12

An "ab initio" calculation of the Carbon-12 elastic form factor, and sum rules of longitudinal and transverse response functions measured in inclusive (e,e') scattering, is reported, based on realistic nuclear potentials and electromagnetic currents. The longitudinal elastic form factor and sum rule are found to be in satisfactory agreement with available experimental data. A direct comparison between theory and experiment is difficult for the transverse sum rule. However, it is shown that the calculated one has large contributions from two-body currents, indicating that these mechanisms lead to a significant enhancement of the quasi-elastic transverse response. This fact may have implications for the anomaly observed in recent neutrino quasi-elastic charge-changing scattering data off Carbon-12.

preprint2013arXiv

GFMC calculations of electromagnetic moments and M1 transitions in $A\leq 9$ nuclei

We present recent Green's function Monte Carlo calculations of magnetic moments and M1 transitions in $A \leq 9$ nuclei, which include corrections arising from two-body meson-exchange electromagnetic currents. Two-body effects provide significant corrections to the calculated observables, bringing them in excellent agreement with the experimental data. In particular, we find that two-body corrections are especially large in the A=9, $T=3/2$ systems, in which they account for up to $\sim 20%$ ($\sim 40 %$) of the total predicted value for the $^9$Li ($^9$C) magnetic moment.

preprint2013arXiv

Nucleon and nucleon-pair momentum distributions in A <= 12 nuclei

We report variational Monte Carlo calculations of single-nucleon momentum distributions for A <= 12 nuclei and nucleon-pair and nucleon-cluster momentum distributions for A <= 8. The wave functions have been generated for a Hamiltonian containing the Argonne v18 two-nucleon and Urbana X three-nucleon potentials. The single-nucleon and nucleon-pair momentum distributions exhibit universal features attributable to the one-pion-exchange tensor interaction. The single-nucleon distributions are broken down into proton and neutron components and spin-up and spin-down components where appropriate. The nucleon-pair momentum distributions are given separately for pp and pn pairs. The nucleon-cluster momentum distributions include dp in 3He, tp and dd in 4He, alpha-d in 6Li, alpha-t in 7Li, and alpha-alpha in 8Be. Detailed tables are provided on-line for download.

preprint2013arXiv

Quantum Monte Carlo calculations of electromagnetic moments and transitions in A <= 9 nuclei with meson-exchange currents derived from chiral effective field theory

Quantum Monte Carlo calculations of electromagnetic moments and transitions are reported for A <= 9 nuclei. The realistic Argonne v18 two-nucleon and Illinois-7 three-nucleon potentials are used to generate the nuclear wave functions. Contributions of two-body meson-exchange current (MEC) operators are included for magnetic moments and M1 transitions. The MEC operators have been derived in both a standard nuclear physics approach and a chiral effective field theory formulation with pions and nucleons including up to one-loop corrections. The two-body MEC contributions provide significant corrections and lead to very good agreement with experiment. Their effect is particularly pronounced in the A=9, T=3/2 systems, in which they provide up to 20% (40 %) of the total predicted value for the 9Li (9C) magnetic moment.

preprint2013arXiv

The proton-proton weak capture in chiral effective field theory

The astrophysical S-factor for proton-proton weak capture is calculated in chiral effective field theory over the center-of-mass relative-energy range 0--100 keV. The chiral two-nucleon potential derived up to next-to-next-to-next-to leading order is augmented by the full electromagnetic interaction including, beyond Coulomb, two-photon and vacuum-polarization corrections. The low-energy constants (LEC's) entering the weak current operators are fixed so as to reproduce the A=3 binding energies and magnetic moments, and the Gamow-Teller matrix element in tritium beta decay. Contributions from S and P partial waves in the incoming two-proton channel are retained. The S-factor at zero energy is found to be S(0)=(4.030 +/- 0.006) x 10^{-23} MeV fm^2, with a P-wave contribution of 0.020 x 10^{-23} MeV fm^2. The theoretical uncertainty is due to the fitting procedure of the LEC's and to the cutoff dependence.

preprint2012arXiv

Inclusive neutrino scattering off deuteron from threshold to GeV energies

Background: Neutrino-nucleus quasi-elastic scattering is crucial to interpret the neutrino oscillation results in long baseline neutrino experiments. There are rather large uncertainties in the cross section, due to insufficient knowledge on the role of two-body weak currents. Purpose: Determine the role of two-body weak currents in neutrino-deuteron quasi-elastic scattering up to GeV energies. Methods: Calculate cross sections for inclusive neutrino scattering off deuteron induced by neutral and charge-changing weak currents, from threshold up to GeV energies, using the Argonne $v_{18}$ potential and consistent nuclear electroweak currents with one- and two-body terms. Results: Two-body contributions are found to be small, and increase the cross sections obtained with one-body currents by less than 10% over the whole range of energies. Total cross sections obtained by describing the final two-nucleon states with plane waves differ negligibly, for neutrino energies $\gtrsim 500$ MeV, from those in which interaction effects in these states are fully accounted for. The sensitivity of the calculated cross sections to different models for the two-nucleon potential and/or two-body terms in the weak current is found to be weak. Comparing cross sections to those obtained in a naive model in which the deuteron is taken to consist of a free proton and neutron at rest, nuclear structure effects are illustrated to be non-negligible. Conclusion: Contributions of two-body currents in neutrino-deuteron quasi-elastic scattering up to GeV are found to be smaller than 10%. Finally, it should be stressed that the results reported in this work do not include pion production channels.

preprint2012arXiv

Quantum Monte Carlo calculations of magnetic moments and M1 transitions in $A \le 9$ nuclei

We present Quantum Monte Carlo calculations of magnetic moments and M1 transitions in $A\le 9$ nuclei which take into account contributions of two-body electromagnetic currents. The Hamiltonian utilized to generate the nuclear wave functions includes the realistic Argonne-{\it v}$_{18}$ two-nucleon and the Illinois-7 three-nucleon interactions. The nuclear two-body electromagnetic currents are derived from a pionful chiral effective field theory including up to one-loop corrections. These currents involve unknown Low Energy Constants which have been fixed so as to reproduce a number of experimental data for the two- and three-nucleon systems, such as $np$ phase shifts and deuteron, triton, and $^3$He magnetic moments. This preliminary study shows that two-body contributions provide significant corrections which are crucial to bring the theory in agreement with the experimental data in both magnetic moments and M1 transitions.

preprint2011arXiv

Chiral effective field theory predictions for muon capture on deuteron and 3He

The muon-capture reactions 2H(μ^-,ν_μ)nn and 3He(μ^-,ν_μ)3H are studied with nuclear strong-interaction potentials and charge-changing weak currents, derived in chiral effective field theory. The low-energy constants (LEC's) c_D and c_E, present in the three-nucleon potential and (c_D) axial-vector current, are constrained to reproduce the A=3 binding energies and the triton Gamow-Teller matrix element. The vector weak current is related to the isovector component of the electromagnetic current via the conserved-vector-current constraint, and the two LEC's entering the contact terms in the latter are constrained to reproduce the A=3 magnetic moments. The muon capture rates on deuteron and 3He are predicted to be 399(3) sec^{-1} and 1494 (21) sec^{-1}, respectively, where the spread accounts for the cutoff sensitivity as well as uncertainties in the LEC's and electroweak radiative corrections. By comparing the calculated and precisely measured rates on 3He, a value for the induced pseudoscalar form factor is obtained in good agreement with the chiral perturbation theory prediction.

preprint2011arXiv

The two-nucleon electromagnetic charge operator in chiral effective field theory ($χ$EFT) up to one loop

The electromagnetic charge operator in a two-nucleon system is derived in chiral effective field theory ($χ$EFT) up to order $e\, Q$ (or N4LO), where $Q$ denotes the low-momentum scale and $e$ is the electric charge. The specific form of the N3LO and N4LO corrections from, respectively, one-pion-exchange and two-pion-exchange depends on the off-the-energy-shell prescriptions adopted for the non-static terms in the corresponding potentials. We show that different prescriptions lead to unitarily equivalent potentials and accompanying charge operators. Thus, provided a consistent set is adopted, predictions for physical observables will remain unaffected by the non-uniqueness associated with these off-the-energy-shell effects.

preprint2010arXiv

Muon capture on deuteron and 3He

The muon capture reactions 2H(μ^-,ν_μ)nn and 3He(μ^-,ν_μ)3H are studied with conventional or chiral realistic potentials and consistent weak currents. The initial and final A=2 and 3 nuclear wave functions are obtained from the Argonne v18 or chiral N3LO two-nucleon potential, in combination with, respectively, the Urbana IX or chiral N2LO three-nucleon potential in the case of A=3. The weak current consists of polar- and axial-vector components. The former are related to the isovector piece of the electromagnetic current via the conserved-vector-current hypothesis. These and the axial currents are derived either in a meson-exchange or in a chiral effective field theory (chiEFT) framework. There is one parameter (either the N-to-Δaxial coupling constant in the meson-exchange model, or the strength of a contact term in the chiEFT model) which is fixed by reproducing the Gamow-Teller matrix element in tritium beta-decay. The model dependence relative to the adopted interactions and currents (and cutoff sensitivity in the chiEFT currents) is weak, resulting in total rates of 392.0 +/- 2.3 Hz for A=2, and 1484 +/- 13 Hz for A=3, where the spread accounts for this model dependence.

preprint2010arXiv

Relativity constraints on the two-nucleon contact interaction

We construct the most general, relativistically invariant, contact Lagrangian at order Q^2 in the power counting, Q denoting the low momentum scale. A complete, but non-minimal, set of (contact) interaction terms is identified, which upon non-relativistic reduction generate 2 leading independent operator combinations of order Q^0 and 7 sub-leading ones of order Q^2 - a result derived previously in the heavy-baryon formulation of effective field theories (EFT's). We show that Poincare covariance of the theory requires that additional terms with fixed coefficients be included, in order to describe the two-nucleon potential in reference frames other than the center-of-mass frame. These terms will contribute in systems with mass number A>2, and their impact on EFT calculations of binding energies and scattering observables in these systems should be studied.

preprint2010arXiv

Solar fusion cross sections II: the pp chain and CNO cycles

We summarize and critically evaluate the available data on nuclear fusion cross sections important to energy generation in the Sun and other hydrogen-burning stars and to solar neutrino production. Recommended values and uncertainties are provided for key cross sections, and a recommended spectrum is given for 8B solar neutrinos. We also discuss opportunities for further increasing the precision of key rates, including new facilities, new experimental techniques, and improvements in theory. This review, which summarizes the conclusions of a workshop held at the Institute for Nuclear Theory, Seattle, in January 2009, is intended as a 10-year update and supplement to Reviews of Modern Physics 70 (1998) 1265.

preprint2010arXiv

The parity-violating asymmetry in the 3He(n,p)3H reaction

The longitudinal asymmetry induced by parity-violating (PV) components in the nucleon-nucleon potential is studied in the charge-exchange reaction 3He(n,p)3H at vanishing incident neutron energies. An expression for the PV observable is derived in terms of T-matrix elements for transitions from the {2S+1}L_J=1S_0 and 3S_1 states in the incoming n-3He channel to states with J=0 and 1 in the outgoing p-3H channel. The T-matrix elements involving PV transitions are obtained in first-order perturbation theory in the hadronic weak-interaction potential, while those connecting states of the same parity are derived from solutions of the strong-interaction Hamiltonian with the hyperspherical-harmonics method. The coupled-channel nature of the scattering problem is fully accounted for. Results are obtained corresponding to realistic or chiral two- and three-nucleon strong-interaction potentials in combination with either the DDH or pionless EFT model for the weak-interaction potential. The asymmetries, predicted with PV pion and vector-meson coupling constants corresponding (essentially) to the DDH "best values" set, range from -9.44 to -2.48 in units of 10^{-8}, depending on the input strong-interaction Hamiltonian. This large model dependence is a consequence of cancellations between long-range (pion) and short-range (vector-meson) contributions, and is of course sensitive to the assumed values for the PV coupling constants.

preprint2010arXiv

Thermal neutron captures on $d$ and $^3$He

We report on a study of the $nd$ and $n\,^3$He radiative captures at thermal neutron energies, using wave functions obtained from either chiral or conventional two- and three-nucleon realistic potentials with the hyperspherical harmonics method, and electromagnetic currents derived in chiral effective field theory up to one loop. The predicted $nd$ and $n\,^3$He cross sections are in good agreement with data, but exhibit a significant dependence on the input Hamiltonian. A comparison is also made between these and new results for the $nd$ and $n\,^3$He cross sections obtained in the conventional framework for both potentials and currents.

preprint2009arXiv

Electromagnetic processes in a $χ$EFT framework

Recently, we have derived a two--nucleon potential and consistent nuclear electromagnetic currents in chiral effective field theory with pions and nucleons as explicit degrees of freedom. The calculation of the currents has been carried out to include N$^3$LO corrections, consisting of two--pion exchange and contact contributions. The latter involve unknown low-energy constants (LECs), some of which have been fixed by fitting the $np$ S- and P-wave phase shifts up to 100 MeV lab energies. The remaining LECs entering the current operator are determined so as to reproduce the experimental deuteron and trinucleon magnetic moments, as well as the $np$ cross section. This electromagnetic current operator is utilized to study the $nd$ and $n^3$He radiative captures at thermal neutron energies. Here we discuss our results stressing on the important role played by the LECs in reproducing the experimental data.

preprint2009arXiv

Electromagnetic Structure and Reactions of Few-Nucleon Systems in $χ$EFT

We summarize our recent work dealing with the construction of the nucleon-nucleon potential and associated electromagnetic currents up to one loop in chiral effective field theory ($χ$EFT). The magnetic dipole operators derived from these currents are then used in hybrid calculations of static properties and low-energy radiative capture processes in few-body nuclei. A preliminary set of results are presented for the magnetic moments of the deuteron and trinucleons and thermal neutron captures on $p$, $d$, and $^3$He.

preprint2008arXiv

Electromagnetic Processes in $χ$EFT

Nuclear electromagnetic currents derived in a chiral-effective-field-theory framework including explicit nucleons, $Δ$ isobars, and pions up to N$^2$LO, {\it i.e.} ignoring loop corrections, are used in a study of neutron radiative captures on protons and deuterons at thermal energies, and of $A$=2 and 3 nuclei magnetic moments. With the strengths of the $Δ$-excitation currents determined to reproduce the $n$-$p$ cross section and isovector combination of the trinucleon magnetic moments, we find that the cross section and photon circular polarization parameter, measured respectively in $n$-$d$ and $\vec{n}$-$d$ processes, are significantly underpredicted by theory.

preprint2008arXiv

Neutron spin rotation in \vec{n}-d scattering

The neutron spin rotation induced by parity-violating (PV) components in the nucleon-nucleon potential is studied in \vec{n}-d scattering at zero energy. Results are obtained corresponding to the Argonne v_{18} two-nucleon and Urbana-IX three-nucleon strong-interaction potentials in combination with either the DDH or pionless EFT model for the weak-interaction potential. We find that this observable is dominated by the contribution of the long-range part of the PV potential associated with pion exchange. Thus its measurement could provide a further constraint, complementary to that coming from measurements of the photon asymmetry in \vec{n}-p radiative capture, on the strength of this component of the hadronic weak interaction.

preprint2008arXiv

Quantum Monte Carlo calculations of magnetic moments and M1 transitions in A<=7 nuclei including meson-exchange currents

Green's function Monte Carlo calculations of magnetic moments and M1 transitions including two-body meson-exchange current (MEC) contributions are reported for A<=7 nuclei. The realistic Argonne v18 two-nucleon and Illinois-2 three-nucleon potentials are used to generate the nuclear wave functions. The two-body meson-exchange operators are constructed to satisfy the continuity equation with the Argonne v18 potential. The MEC contributions increase the A=3,7 isovector magnetic moments by 16% and the A=6,7 M1 transition rates by 17--34%, bringing them into very good agreement with the experimental data.

preprint2005arXiv

Electromagnetic structure of A=2 and 3 nuclei and the nuclear current operator

Different models for conserved two- and three-body electromagnetic currents are constructed from two- and three-nucleon interactions, using either meson-exchange mechanisms or minimal substitution in the momentum dependence of these interactions. The connection between these two different schemes is elucidated. A number of low-energy electronuclear observables, including (i) $np$ radiative capture at thermal neutron energies and deuteron photodisintegration at low energies, (ii) $nd$ and $pd$ radiative capture reactions, and (iii) isoscalar and isovector magnetic form factors of $^3$H and $^3$He, are calculated in order to make a comparative study of these models for the current operator. The realistic Argonne $v_{18}$ two-nucleon and Urbana IX or Tucson-Melbourne three-nucleon interactions are taken as a case study. For $A$=3 processes, the bound and continuum wave functions, both below and above deuteron breakup threshold, are obtained with the correlated hyperspherical-harmonics method. Three-body currents give small but significant contributions to some of the polarization observables in the $^2$H($p,γ$)$^3$He process and the $^2$H($n,γ$)$^3$H cross section at thermal neutron energies. It is shown that the use of a current which did not exactly satisfy current conservation with the two- and three-nucleon interactions in the Hamiltonian was responsible for some of the discrepancies reported in previous studies between the experimental and theoretical polarization observables in $pd$ radiative capture.

preprint2005arXiv

Induced Polarization in the $^2$H($γ,\vec n$)$^1$H Reaction at Low Energy

The induced polarization, $P^\prime_y$, of the neutron in the deuteron photo-disintegration from threshold up to 30 MeV is calculated using a variety of different, latest-generation potentials--Argonne $v_{18}$, Bonn 2000, and Nijmegen I-- and a realistic model for the nuclear electromagnetic current operator, including one- and two-body terms. The model dependence of the theoretical predictions is found to be very small. These predictions are systematically larger in magnitude than the measured $P^\prime_y$ values, and corroborate the conclusions of an earlier, and much older, study. There is considerable scatter in the available experimental data. New and more accurate measurements of the induced polarization in the $^2$H($γ,\vec n$)$^1$H reaction are needed in order to establish unequivocally whether there is a discrepancy between theory and experiment.

preprint2005arXiv

Two-Body Electrodisintegration of $^3$He at High Momentum Transfer

The $^3$He$(e,e^\prime p)$$d$ reaction is studied using an accurate three-nucleon bound state wave function, a model for the electromagnetic current operator including one- and two-body terms, and the Glauber approximation for the treatment of final state interactions. In contrast to earlier studies, the profile operator in the Glauber expansion is derived from a nucleon-nucleon scattering amplitude, which retains its full spin and isospin dependence and is consistent with phase-shift analyses of two-nucleon scattering data. The amplitude is boosted from the center-of-mass frame, where parameterizations for it are available, to the frame where rescattering occurs. Exact Monte Carlo methods are used to evaluate the relevant matrix elements of the electromagnetic current operator. The predicted cross section is found to be in quantitative agreement with the experimental data for values of the missing momentum $p_{\rm m}$ in the range (0--700) MeV/c, but underestimates the data at $p_{\rm m} \simeq 1$ GeV/c by about a factor of two. However, the longitudinal-transverse asymmetry, measured up to $p_{\rm m} \simeq$ 600 MeV/c, is well reproduced by theory. A critical comparison is carried out between the results obtained in the present work and those of earlier studies.

preprint2004arXiv

Electrodisintegration of $^3$He below and above deuteron breakup threshold

Recent advances in the study of electrodisintegration of 3He are presented and discussed. The pair-correlated hyperspherical harmonics method is used to calculate the initial and final state wave functions, with a realistic Hamiltonian consisting of the Argonne v18 two-nucleon and Urbana IX three-nucleon interactions. The model for the nuclear current and charge operators retains one- and many-body contributions. Particular attention is made in the construction of the two-body current operators arising from the momentum-dependent part of the two-nucleon interaction. Three-body current operators are also included so that the full current operator is strictly conserved. The present model for the nuclear current operator is tested comparing theoretical predictions and experimental data of pd radiative capture cross section and spin observables.

preprint2004arXiv

Parity-Violating Interaction Effects in the np System

We investigate parity-violating observables in the np system, including the longitudinal asymmetry and neutron-spin rotation in np elastic scattering, the photon asymmetry in np radiative capture, and the asymmetries in deuteron photo-disintegration d(gamma,n)p in the threshold region and electro-disintegration d(e,e`)np in quasi-elastic kinematics. To have an estimate of the model dependence for the various predictions, a number of different, latest-generation strong-interaction potentials--Argonne v18, Bonn 2000, and Nijmegen I--are used in combination with a weak-interaction potential consisting of pi-, rho-, and omega-meson exchanges--the model known as DDH. The complete bound and scattering problems in the presence of parity-conserving, including electromagnetic, and parity-violating potentials is solved in both configuration and momentum space. The issue of electromagnetic current conservation is examined carefully. We find large cancellations between the asymmetries induced by the parity-violating interactions and those arising from the associated pion-exchange currents. In the np capture, the model dependence is nevertheless quite small, because of constraints arising through the Siegert evaluation of the relevant E1 matrix elements. In quasi-elastic electron scattering these processes are found to be insignificant compared to the asymmetry produced by gamma-Z interference on individual nucleons.

preprint2004arXiv

Polarization transfer in $^4$He$(\vec{e},e^\prime\vec{p})$$^3$H: is the ratio $G_{Ep}/G_{Mp}$ modified in medium ?

Polarization observables in the $^4$He$(\vec{e},e^\prime\vec{p})$$^3$H reaction are calculated using accurate three- and four-nucleon bound-state wave functions, a realistic model for the nuclear electromagnetic current operator, and a treatment of final-state-interactions with an optical potential. In contrast to earlier studies, no significant discrepancies are found between theory and experiment both for the ratio of transverse to longitudinal polarization transfers and for the induced polarization, when free-nucleon electromagnetic form factors are used in the current operator. The present results challenge the current interpretation of the experimental data in terms of medium-modified form factors.

preprint2002arXiv

Elastic e-d Scattering Data and the Deuteron Wave Function

What range of momentum components in the deuteron wave function are available e d elastic scattering data sensitive to ? This question is addressed within the context of a model calculation of the deuteron form factors, based on realistic interactions and currents. It is shown that the data on the $A(q)$, $B(q)$, and $T_{20}(q)$ observables at $q \leq 6$ fm$^{-1}$ essentially probe momentum components up to $\simeq 4 m_π$.

preprint2002arXiv

Parameter-free effective field theory calculation for the solar proton-fusion and hep processes

Spurred by the recent complete determination of the weak currents in two-nucleon systems up to ${\cal O}(Q^3)$ in heavy-baryon chiral perturbation theory, we carry out a parameter-free calculation of the threshold $S$-factors for the solar $pp$ (proton-fusion) and $hep$ processes in an effective field theory that {\it combines} the merits of the standard nuclear physics method and systematic chiral expansion. The power of the EFT adopted here is that one can correlate in a unified formalism the weak-current matrix elements of two-, three- and four-nucleon systems. Using the tritium $β$-decay rate as an input to fix the only unknown parameter in the theory, we can evaluate the threshold $S$ factors with drastically improved precision; the results are $S_{pp}(0) = 3.94\times(1 \pm 0.004) \times 10^{-25} {MeV-b}$ and $S_{hep}(0) = (8.6\pm 1.3)\times 10^{-20} {keV-b}$. The dependence of the calculated $S$-factors on the momentum cutoff parameter $Λ$ has been examined for a physically reasonable range of $Λ$. This dependence is found to be extremely small for the $pp$ process, and to be within acceptable levels for the $hep$ process, substantiating the consistency of our calculational scheme.

preprint2002arXiv

Parity Violating Interactions and Currents in the Deuteron

We investigate parity-violating asymmetries in polarized n p radiative capture, and deuteron electro-disintegration in quasi-elastic kinematics, using the DDH model for the parity-violating nucleon-nucleon interaction. We find dramatic cancellations between the asymmetries induced by the parity-violating interaction and those arising from the associated parity-violating pion-exchange currents. In np capture, the model-dependence of the result is nevertheless quite small because of constraints arising through the Siegert evaluation of the relevant E1 matrix element. In quasi-elastic electron scattering these processes are found to be insignificant comared to the asymmetry produced by the gamma-Z interference on individual nucleons. These two experiments, then, provide clean probes of different aspects of weak-interaction physics associated with parity violation in the np system.

preprint2002arXiv

Theoretical Study of the ^3He(mu^-,nu_mu)^3H Capture

The ^3He(mu^-,nu_mu)^3H weak capture is studied using correlated-hyperspherical-harmonics wave functions, obtained from realistic Hamiltonians consisting of the Argonne $v_{14}$ or Argonne $v_{18}$ two-nucleon, and Tucson-Melbourne or Urbana-IX three-nucleon interactions. The nuclear weak charge and current operators have vector and axial-vector components that include one- and two-body contributions. The strength of the leading two-body operator in the axial-vector current is adjusted to reproduce the Gamow-Teller matrix element in tritium $β$-decay. The calculated total capture rate is in excellent agreement with the most recent experimental determination $1496\pm 4$ sec$^{-1}$, when the PCAC value is adopted for the induced pseudo-scalar coupling constant $g_{PS}$. The predictions for the capture rate and angular correlation parameters $A_v$, $A_t$, and $A_Δ$ are found to be only very weakly dependent on the model input Hamiltonian. The variation of these observables with $g_{PS}$ and the theoretical uncertainties deriving from the model-dependent procedure used to constrain the axial current are investigated.

preprint2001arXiv

Longitudinal and Transverse Quasi-Elastic Response Functions of Light Nuclei

The $^3$He and $^4$He longitudinal and transverse response functions are determined from an analysis of the world data on quasi-elastic inclusive electron scattering. The corresponding Euclidean response functions are derived and compared to those calculated with Green's function Monte Carlo methods, using realistic interactions and currents. Large contributions associated with two-body currents are found, particularly in the $^4$He transverse response, in agreement with data. The contributions of two-body charge and current operators in the $^3$He, $^4$He, and $^6$Li response functions are also studied via sum-rule techniques. A semi-quantitative explanation for the observed systematics in the excess of transverse quasi-elastic strength, as function of mass number and momentum transfer, is provided. Finally, a number of model studies with simplified interactions, currents, and wave functions is carried out to elucidate the role played, in the full calculation, by tensor interactions and correlations.

preprint2001arXiv

Parity-Violating Interaction Effects I: the Longitudinal Asymmetry in pp Elastic Scattering

The proton-proton parity-violating longitudinal asymmetry is calculated in the lab-energy range 0--350 MeV, using a number of different, latest-generation strong-interaction potentials--Argonne V18, Bonn-2000, and Nijmegen-I--in combination with a weak-interaction potential consisting of rho- and omega-meson exchanges--the model known as DDH. The complete scattering problem in the presence of parity-conserving, including Coulomb, and parity-violating potentials is solved in both configuration- and momentum-space. The predicted parity-violating asymmetries are found to be only weakly dependent upon the input strong-interaction potential adopted in the calculation. Values for the rho- and omega-meson weak coupling constants $h^{pp}_ρ$ and $h^{pp}_ω$ are determined by reproducing the measured asymmetries at 13.6 MeV, 45 MeV, and 221 MeV.

preprint2000arXiv

Electro-weak Reactions for Astrophysics

The status of "ab initio" microscopic calculations of the 2H(p,gamma)3He and 3He(p,e+nu_e)4He reactions is reviewed. The methods used to generate accurate nuclear ground- and scattering-state wave functions, and to construct realistic electro-weak transition operators are described. The uncertainties in the theoretical predictions, particularly those relevant to the p-3He weak capture, are discussed. For the dp radiative capture, the theoretical results are compared with the TUNL data in the energy range 0--100 keV.

preprint2000arXiv

Parity-Violating Electron-Deuteron Scattering

The longitudinal asymmetry due to $Z^0$ exchange is calculated in quasi-elastic electron-deuteron scattering at momentum transfers $|Q^2| \simeq 0.1$ GeV$^2$ relevant for the SAMPLE experiment. The deuteron and $pn$ scattering-state wave functions are obtained from solutions of a Schrödinger equation with the Argonne $v_{18}$ potential. Electromagnetic and weak neutral one- and two-nucleon currents are included in the calculation. The two-nucleon currents of pion range are shown to be identical to those derived in Chiral Perturbation Theory. The results indicate that two-body contributions to the asymmetry are small ($\simeq$ 0.2%) around the quasi-elastic peak, but become relatively more significant ($\simeq$ 3%) in the high-energy wing of the quasi-elastic peak.

preprint1999arXiv

Photo- and Electro-Disintegration of 3He at Threshold and pd Radiative Capture

The present work reports results for: pd radiative capture observables measured at center-of-mass (c.m.) energies in the range 0--100 keV and at 2 MeV by the TUNL and Wisconsin groups, respectively; contributions to the Gerasimov-Drell-Hearn (GDH) integral in 3He from the two- up to the three-body breakup thresholds, compared to experimental determinations by the TUNL group in this threshold region; longitudinal, transverse, and interference response functions measured in inclusive polarized electron scattering off polarized 3He at excitation energies below the threshold for breakup into ppn, compared to unpolarized longitudinal and transverse data from the Saskatoon group. The calculations are based on a realistic Hamiltonian with two- and three-nucleon interactions and a realistic current operator, including one- and two-body components. The theoretical predictions obtained by including only one-body currents are in violent disagreement with data. These differences between theory and experiment are, to a large extent, removed when two-body currents are taken into account, although some rather large discrepancies remain in the c.m. energy range 0--100 keV, particularly for the pd differential cross section and tensor analyzing power at small angles, and contributions to the GDH integral. A rather detailed analysis indicates that these discrepancies have, in large part, a common origin, and can be traced back to an excess strength obtained in the theoretical calculation of the E1 reduced matrix element associated with the pd channel having L,S,J=1,1/2,3/2. It is suggested that this lack of E1 strength observed experimentally might have implications for the nuclear interaction at very low energies. Finally, the validity of the long-wavelength approximation for electric dipole transitions is discussed.

preprint1998arXiv

Electromagnetic Structure of the Trinucleons

The electromagnetic form factors of the trinucleons 3H and 3He are calculated with wave functions obtained with the Argonne AV18 two-nucleon and Urbana IX three-nucleon interactions. Full account is taken of the two-body currents required by current conservation with the AV18 interaction as well as those associated with N-Delta transition currents and the currents of Delta resonance components in the wave functions. Explicit three-nucleon current operators associated with the two-pion exchange three-nucleon interaction arising from irreducible S-wave pion-nucleon scattering is constructed and shown to have very little effect on the calculated magnetic form factors. The calculated magnetic form factor of 3H, and charge form factors of both 3H and 3He are in satisfactory agreement with the experimental data. However, the position of the zero in the magnetic form factor of 3He is slightly underpredicted.