Source author record

M. Viviani

M. Viviani 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

62works
15topics
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

62 published item(s)

preprint2026arXiv

$A=2,3,4$ nuclear contact coefficients in the generalized contact formalism

We present a theoretical calculation for the $A = 2, 3$ and 4 nuclear contact coefficients within the generalized contact formalism, using both local and non-local chiral potentials. The Hyperspherical Harmonics method is employed to calculate the nuclear wave functions, from which we derive two-body momentum distributions and density functions to extract the contact coefficients. We have extracted the contact coefficients from two-body momentum distributions or from density functions, for a given nucleus and potential, and we have found that the generalized contact formalism predictions are verified in the triplet spin channel for local and non-local potentials. On the other hand, some significant tensions exist for the singlet channels, especially when studied with non-local potentials. We have also analyzed the model-independence of the ratios between the contact coefficients, and we have found to be quite satisfied. This study extends previous works based on local interaction models only.

preprint2022arXiv

Heisenberg's uncertainty principle in the PTOLEMY project: a theory update

We discuss the consequences of the quantum uncertainty on the spectrum of the electron emitted by the $β$-processes of a tritium atom bound to a graphene sheet. We analyze quantitatively the issue recently raised in [Cheipesh et al., Phys. Rev. D 104, 116004 (2021)], and discuss the relevant time scales and the degrees of freedom that can contribute to the intrinsic spread in the electron energy. We perform careful calculations of the potential between tritium and graphene with different coverages and geometries. With this at hand, we propose possible avenues to mitigate the effect of the quantum uncertainty.

preprint2022arXiv

Implementation and Optimization of the PTOLEMY Transverse Drift Electromagnetic Filter

The PTOLEMY transverse drift filter is a new concept to enable precision analysis of the energy spectrum of electrons near the tritium beta-decay endpoint. This paper details the implementation and optimization methods for successful operation of the filter. We present the first demonstrator that produces the required magnetic field properties with an iron return-flux magnet. Two methods for the setting of filter electrode voltages are detailed. The challenges of low-energy electron transport in cases of low field are discussed, such as the growth of the cyclotron radius with decreasing magnetic field, which puts a ceiling on filter performance relative to fixed filter dimensions. Additionally, low pitch angle trajectories are dominated by motion parallel to the magnetic field lines and introduce non-adiabatic conditions and curvature drift. To minimize these effects and maximize electron acceptance into the filter, we present a three-potential-well design to simultaneously drain the parallel and transverse kinetic energies throughout the length of the filter. These optimizations are shown, in simulation, to achieve low-energy electron transport from a 1 T iron core (or 3 T superconducting) starting field with initial kinetic energy of 18.6 keV drained to <10 eV (<1 eV) in about 80 cm. This result for low field operation paves the way for the first demonstrator of the PTOLEMY spectrometer for measurement of electrons near the tritium endpoint to be constructed at the Gran Sasso National Laboratary (LNGS) in Italy.

preprint2021arXiv

Hunting down the cause of solar magnetism

To understand solar and stellar dynamos combining local and global numerical modelling with long-term observations is a challenging task: even with state of the art computational methods and resources, the stellar parameter regime remains unattainable. Our goal is to relax some approximations, in order to simulate more realistic systems, and try to connect the results with theoretical predictions and state-of-the-art observations. We present here the first test-field measurements from our higher-resolution runs with improved heat conduction description. They indicate significant changes in the profiles of the most crucial inductive effect related to solar and stellar dynamo mechanisms. Higher resolution runs, currently undertaken, will bring us into an even more turbulent regime, in which we will be able to study, for the first time, the interaction of small- and large-scale dynamos in a quantitative way.

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

Parity- and time-reversal-violating nuclear forces

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.

preprint2020arXiv

Physically motivated heat conduction treatment in simulations of solar-like stars: effects on dynamo transitions

Context. Results from global magnetoconvection simulations of solar-like stars are at odds with observations in many respects: They show a surplus of energy in the kinetic power spectrum at large scales, anti-solar differential rotation profiles, with accelerated poles and a slow equator, for the solar rotation rate, and a transition from axi- to non-axisymmetric dynamos at a much lower rotation rate than what is observed. Even though the simulations reproduce the observed active longitudes in fast rotators, their motion in the rotational frame (the so-called azimuthal dynamo wave, ADW) is retrograde, in contrast to the prevalent prograde motion in observations. Aims. We study the effect of a more realistic treatment of heat conductivity in alleviating the discrepancies between observations and simulations. Methods. We use physically-motivated heat conduction, by applying Kramers opacity law, on a semi-global spherical setup describing convective envelopes of solar-like stars, instead of a prescribed heat conduction profile from mixing-length arguments. Results. We find that some aspects of the results now better correspond to observations: The axi- to non-axisymmetric transition point is shifted towards higher rotation rates. We also find a change in the propagation direction of ADWs so that also prograde waves are now found. The transition from anti-solar to solar-like rotation profile, however, is also shifted towards higher rotation rates, leaving the models into an even more unrealistic regime. Conclusions. Although a Kramers-based heat conduction does not help in reproducing the solar rotation profile, it does help in the faster rotation regime, where the dynamo solutions now match better with observations.

preprint2020arXiv

Shapes of stellar activity cycles

Magnetic activity cycles are an important phenomenon in both the Sun and other stars. The shape of the solar cycle is commonly characterised by a fast rise and a slower decline, but not much attention has been paid to the shape of cycles in other stars. Our aim is to study whether the asymmetric shape of the solar cycle is common in other stars as well, and compare the cycle asymmetry to other stellar parameters. We also study the differences in the shape of the solar cycle, depending on what activity indicator is used. The observations are also compared to simulated activity cycles. We use the chromospheric Ca II H&K data from the Mount Wilson Observatory HK Project. From this data set we identify 47 individual cycles from 18 stars. We use the statistical skewness of a cycle as a measure of its asymmetry, and compare this to other stellar parameters. A similar analysis has been done to magnetic cycles extracted from direct numerical magnetohydrodynamic simulations of solar-type convection zones. The shape of the solar cycle (fast rise and slower decline) is common in other stars as well, although the Sun has particularly asymmetric cycles. Cycle-to-cycle variations are large, but the average shape of a cycle is still fairly well represented by a sinusoid. We find only slight correlations between the cycle asymmetry and other stellar parameters. There are large differences in the shape of the solar cycle, depending on what activity indicator is used. In the simulated cycles, there is a difference in the symmetry of global simulations covering the full longitudinal range, hence capable of exciting non-axisymmetric dynamo modes, versus wedge simulations covering a partial extent in longitude, where only axisymmetric modes are possible. The former produce preferentially positive skewness, while the latter a negative one.

preprint2020arXiv

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

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.

preprint2020arXiv

Subleading contributions to the three-nucleon contact interaction

We obtain a minimal form of the two-derivative three-nucleon contact Lagrangian, by imposing all constraints deriving from discrete symmetries, Fierz identities and Poincare' covariance. The resulting interaction, depending on 13 unknown low-energy constants, leads to a three-nucleon potential which we give in a local form in configuration space. We also consider the leading (no-derivative) four-nucleon interaction and show that there exists only one independent operator.

preprint2020arXiv

Unitary ambiguity of NN contact interactions and the 3N force

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.

preprint2019arXiv

Time Reversal Violation in Light Nuclei

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.

preprint2016arXiv

Advances in solving the two-fermion homogeneous Bethe-Salpeter equation in Minkowski space

Actual solutions of the Bethe-Salpeter equation for a two-fermion bound system are becoming available directly in Minkowski space, by virtue of a novel technique, based on the so-called Nakanishi integral representation of the Bethe-Salpeter amplitude and improved by expressing the relevant momenta through light-front components, i.e. $k^\pm=k^0 \pm k^3$. We solve a crucial problem that widens the applicability of the method to real situations by providing an analytically exact treatment of the singularities plaguing the two-fermion problem in Minkowski space, irrespective of the complexity of the irreducible Bethe-Salpeter kernel. This paves the way for feasible numerical investigations of relativistic composite systems, with any spin degrees of freedom. We present a thorough comparison with existing numerical results, evaluated in both Minkowski and Euclidean space, fully corroborating our analytical treatment, as well as fresh light-front amplitudes illustrating the potentiality of non perturbative calculations performed directly in Minkowski space.

preprint2016arXiv

Bethe-Salpeter bound-state structure in Minkowski space

The quantitative investigation of the scalar Bethe-Salpeter equation in Minkowski space, within the ladder-approximation framework, is extended to include the excited states. This study has been carried out for an interacting system composed by two massive bosons exchanging a massive scalar, by adopting (i) the Nakanishi integral representation of the Bethe-Salpeter amplitude, and (ii) the formally exact projection onto the null plane. Our analysis, on one hand, confirms the reliability of the method already applied to the ground state and, on the other one, extends the investigation from the valence distribution in momentum space to the corresponding quantity in the impact-parameter space, pointing out some relevant features, like (i) the equivalence between Minkowski and Euclidean transverse-momentum amplitudes, and (ii) the leading exponential fall-off of the valence wave function in the impact-parameter space.

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

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

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

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.

preprint2015arXiv

Solving the inhomogeneous Bethe-Salpeter Equation in Minkowski space: the zero-energy limit

For the first time, the inhomogeneous Bethe-Salpeter Equation for an interacting system, composed by two massive scalars exchanging a massive scalar, is numerically investigated in ladder approximation, directly in Minkowski space, by using an approach based on the Nakanishi integral representation. In this paper, the limiting case of zero-energy states is considered, extending the approach successfully applied to bound states. The numerical values of scattering lengths, are calculated for several values of the Yukawa coupling constant, by using two different integral equations that stem within the Nakanishi framework. Those low-energy observables are compared with (i) the analogous quantities recently obtained in literature, within a totally different framework and (ii) the non relativistic evaluations, for illustrating the relevance of a non perturbative, genuine field theoretical treatment in Minkowski space, even in the low-energy regime. Moreover, dynamical functions, like the Nakanishi weight functions and the distorted part of the zero-energy Light-front wave functions are also presented. Interestingly, a highly non trivial issue related to the abrupt change in the width of the support of the Nakanishi weight function, when the zero-energy limit is approached, is elucidated, ensuring a sound basis to the forthcoming evaluation of phase-shifts.

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

Breakup of three particles within the adiabatic expansion method

General expressions for the breakup cross sections in the lab frame for $1+2$ reactions are given in terms of the hyperspherical adiabatic basis. The three-body wave function is expanded in this basis and the corresponding hyperradial functions are obtained by solving a set of second order differential equations. The ${\cal S}$-matrix is computed by using two recently derived integral relations. Even though the method is shown to be well suited to describe $1+2$ processes, there are nevertheless particular configurations in the breakup channel (for example those in which two particles move away close to each other in a relative zero-energy state) that need a huge number of basis states. This pathology manifests itself in the extremely slow convergence of the breakup amplitude in terms of the hyperspherical harmonic basis used to construct the adiabatic channels. To overcome this difficulty the breakup amplitude is extracted from an integral relation as well. For the sake of illustration, we consider neutron-deuteron scattering. The results are compared to the available benchmark calculations.

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.

preprint2013arXiv

Effect of three-nucleon interaction in p-3He elastic scattering

We present a detailed study of the effect of different three-nucleon interaction models in p-3He elastic scattering at low energies. In particular, two models have been considered: one derived from effective field theory at next-to-next-to-leading order and one derived from a more phenomenological point of view -- the so-called Illinois model. The four-nucleon scattering observables are calculated using the Kohn variational principle and the hyperspherical harmonics technique and the results are compared with available experimental data. We have found that the inclusion of either one of the other force model improves the agreement with the experimental data, in particular for the proton vector analyzing power.

preprint2013arXiv

JLab Measurement of the $^4$He Charge Form Factor at Large Momentum Transfers

The charge form factor of $^$4He has been extracted in the range 29 fm$^{-2}$ $\le Q^2 \le 77$ fm$^{-2}$ from elastic electron scattering, detecting $^4$He nuclei and electrons in coincidence with the High Resolution Spectrometers of the Hall A Facility of Jefferson Lab. The results are in qualitative agreement with realistic meson-nucleon theoretical calculations. The data have uncovered a second diffraction minimum, which was predicted in the $Q^2$ range of this experiment, and rule out conclusively long-standing predictions of dimensional scaling of high-energy amplitudes using quark counting.

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

Effect of three nucleon forces in p-3He scattering

The effect of the inclusion of different models of three nucleon (3N) forces in p-3He elastic scattering at low energies is studied. Two models have been considered: one derived from effective field theory at next-to-next-to-leading order and one derived from a more phenomenological point of view -- the so-called Illinois model. The four nucleon scattering observables are calculated using the Kohn variational principle and the hyperspherical harmonic technique and the results are compared with available experimental data. We have found that with the inclusion of both 3N force models the agreement with the experimental data is improved, in particular for the proton vector analyzing power A_y.

preprint2012arXiv

Efimov Physics in small bosonic clusters

We study small clusters of bosons, $A=2,3,4,5,6$, characterized by a resonant interaction. Firstly, we use a soft-gaussian interaction that reproduces the values of the dimer binding energy and the atom-atom scattering length obtained with LM2M2 potential, a widely used $^4$He-$^4$He interaction. We change the intensity of the potential to explore the clusters' spectra in different regions with large positive and large negative values of the two-body scattering length and we report the clusters' energies on Efimov plot, which makes the scale invariance explicit. Secondly, we repeat our calculation adding a repulsive three-body force to reproduce the trimer binding energy. In all the region explored, we have found that these systems present two states, one deep and one shallow close to the $A-1$ threshold, and scale invariance has been investigated for these states. The calculations are performed by means of Hyperspherical Harmonics basis set.

preprint2012arXiv

Energy spectra of small bosonic clusters having a large two-body scattering length

In this work we investigate small clusters of bosons using the hyperspherical harmonic basis. We consider systems with $A=2,3,4,5,6$ particles interacting through a soft inter-particle potential. In order to make contact with a real system, we use an attractive gaussian potential that reproduces the values of the dimer binding energy and the atom-atom scattering length obtained with one of the most widely used $^4$He-$^4$He interactions, the LM2M2 potential. The intensity of the potential is varied in order to explore the clusters' spectra in different regions with large positive and large negative values of the two-body scattering length. In addition, we include a repulsive three-body force to reproduce the trimer binding energy. With this model, consisting in the sum of a two- and three-body potential, we have calculated the spectrum of the four, five and six particle systems. In all the region explored, we have found that these systems present two bound states, one deep and one shallow close to the $A-1$ threshold. Some universal relations between the energy levels are extracted; in particular, we have estimated the universal ratios between thresholds of the three-, four-, and five-particle continuum using the two-body gaussian

preprint2012arXiv

Integral relations and the adiabatic expansion method for 1+2 reactions above the breakup threshold: Helium trimers with soft-core potentials

The integral relations formalism introduced in \cite{bar09,rom11}, and designed to describe 1+$N$ reactions, is extended here to collision energies above the threshold for the target breakup. These two relations are completely general, and in this work they are used together with the adiabatic expansion method for the description of 1+2 reactions. The neutron-deuteron breakup, for which benchmark calculations are available, is taken as a test of the method. The s-wave collision between the $^4$He atom and $^4$He$_2$ dimer above the breakup threshold and the possibility of using soft-core two-body potentials plus a short-range three-body force will be investigated. Comparisons to previous calculations for the three-body recombination and collision dissociation rates will be shown.

preprint2012arXiv

Six-bodies calculations using the Hyperspherical Harmonics method

In this work we show results for light nuclear systems and small clusters of helium atoms using the hyperspherical harmonics basis. We use the basis without previous symmetrization or antisymmetrization of the state. After the diagonalization of the Hamiltonian matrix, the eigenvectors have well defined symmetry under particle permutation and the identification of the physical states is possible. We show results for systems composed up to six particles. As an example of a fermionic system, we consider a nucleon system interacting through the Volkov potential, used many times in the literature. For the case of bosons, we consider helium atoms interacting through a potential model which does not present a strong repulsion at short distances. We have used an attractive gaussian potential to reproduce the values of the dimer binding energy, the atom-atom scattering length, and the effective range obtained with one of the most widely used He-He interaction, the LM2M2 potential. In addition, we include a repulsive hypercentral three-body force to reproduce the trimer binding energy.

preprint2012arXiv

Three-nucleon contact interaction at the next-to-leading order

A satisfactory description of bound and scattering states of the three-nucleon (3N) system is still lacking, contrary to what happens in the two-nucleon case. In the framework of chiral effective theory, it is possible that a realistic 3N interaction will require the inclusion of subleading contact terms, which are unconstrained by chiral symmetry. We construct the minimal 3N contact Lagrangian imposing all constraints from the discrete symmetries, Fierz identities and Poincaré covariance, and show that it consists of 10 independent operators. By integrating out the pions from the effective theory we derive the pion-induced contributions to the corresponding low-energy constants (LECs), and find no clear signal of convergence of the pionful theory between N2LO and N3LO.

preprint2011arXiv

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

The n-3H and p-3He elastic phase-shifts below the trinucleon disintegration thresholds are calculated by solving the 4-nucleon problem with three different realistic nucleon-nucleon interactions (the I-N3LO model by Entem and Machleidt, the Argonne v18 potential model, and a low-k model derived from the CD-Bonn potential). Three different methods -- Alt, Grassberger and Sandhas, Hyperspherical Harmonics, and Faddeev-Yakubovsky -- have been used and their respective results are compared. For both n-3H and p-3He we observe a rather good agreement between the three different theoretical methods. We also compare the theoretical predictions with the available experimental data, confirming the large underprediction of the p-3He analyzing power.

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

General integral relations for the description of scattering states using the hyperspherical adiabatic basis

In this work we investigate 1+2 reactions within the framework of the hyperspherical adiabatic expansion method. To this aim two integral relations, derived from the Kohn variational principle, are used. A detailed derivation of these relations is shown. The expressions derived are general, not restricted to relative $s$ partial waves, and with applicability in multichannel reactions. The convergence of the ${\cal K}$-matrix in terms of the adiabatic potentials is investigated. Together with a simple model case used as a test for the method, we show results for the collision of a $^4$He atom on a \dimer dimer (only the elastic channel open), and for collisions involving a $^6$Li and two $^4$He atoms (two channels open).

preprint2011arXiv

Spectra of helium clusters with up to six atoms using soft core potentials

In this work we investigate small clusters of helium atoms using the hyperspherical harmonic basis. We consider systems with $A=2,3,4,5,6$ atoms with an inter-particle potential which does not present a strong repulsion at short distances. We use an attractive gaussian potential that reproduces the values of the dimer binding energy, the atom-atom scattering length, and the effective range obtained with one of the widely used He-He interactions, the LM2M2 potential. In systems with more than two atoms we consider a repulsive three-body force that, by construction, reproduces the trimer binding energy of the LM2M2 potential. With this model, consisting in the sum of a two- and three-body potential, we have calculated the spectrum of clusters formed by four, five and six helium atoms. We have found that these systems present two bound states, one deep and one shallow close to the threshold fixed by the energy of the $(A-1)$-atom system. Universal relations between the energies of the excited state of the $A$-atom system and the ground state energy of the $(A-1)$-atom system are extracted as well as the ratio between the ground state of the $A$-atom system and the ground state energy of the trimer.

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.

preprint2011arXiv

Theoretical description of three- and four-nucleon scattering states using bound-state-like wave functions

Bound-state-like wave functions are used to determine the scattering matrix corresponding to low energy $N-d$ and $p-^3$He collisions. To this end, the coupled channel form of the integral relations derived from the Kohn variational principle is used. The construction of degenerate bound-state-like wave functions belonging to the continuum spectrum of the Hamiltonian is discussed. Examples are shown using realistic nucleon-nucleon forces.

preprint2010arXiv

Comparative study of three-nucleon force models in $A=3,4$ systems

Using modern nucleon-nucleon interactions in the description of the $A=3,4$ nuclei, it is not possible to reproduce both the three- and four-nucleon binding energies simultaneously. This is one manifestation of the necessity of including a three-nucleon force in the nuclear Hamiltonian. In this paper we will perform a comparative study of some, widely used, three-nucleon force models. We will analyze their capability to describe the aforementioned binding energies as well as the $n-d$ doublet scattering length. A correct description of these quantities can be considered a stringent requirement for a nuclear Hamiltonian containing two- and three-nucleon interaction terms. As we will show, this requirement is not fulfilled by several of the models available in the literature. To satisfy it, we propose modifications in the parametrization of the three-nucleon forces and we study their effects on few selected $N-d$ low energy scattering observables.

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

Non-symmetrized hyperspherical harmonic basis for $A$-bodies

The use of the hyperspherical harmonic (HH) basis in the description of bound states in an $A$-body system composed by identical particles is normally preceded by a symmetrization procedure in which the statistic of the system is taken into account. This preliminary step is not strictly necessary; the direct use of the HH basis is possible, even if the basis has not a well defined behavior under particle permutations. In fact, after the diagonalization of the Hamiltonian matrix, the eigenvectors reflect the symmetries present in it. They have well defined symmetry under particle permutation and the identification of the physical states is possible, as it will be shown in specific cases. The problem related to the large degeneration of the basis is circumvented by constructing the Hamiltonian matrix as a sum of products of sparse matrices. This particular representation of the Hamiltonian is well suited for a numerical iterative diagonalization, where only the action of the matrix on a vector is needed. As an example we compute bound states for systems with $A=3-6$ particles interacting through a short-range central interaction. We also consider the case in which the potential is restricted to act in relative s-waves with and without the inclusion of the Coulomb potential. This very simple model predicts results in qualitative good agreement with the experimental data and it represents a first step in a project dedicated to the use of the HH basis to describe bound and low energy scattering states in light nuclei.

preprint2010arXiv

Proton-3He elastic scattering at low energies and the "A_y Puzzle"

The Kohn variational principle and the hyperspherical harmonic technique are applied to study p-3He elastic scattering at low energies. Preliminary results obtained using several interaction models are reported. The calculations are compared to a recent phase shift analysis performed at the Triangle University Nuclear Laboratory and to the available experimental data. Using a three-nucleon interaction derived from chiral perturbation theory at N2LO, we have found a noticeable reduction of the discrepancy observed for the A_y observable.

preprint2010arXiv

Relativistic and spin effects in elastic backward p-d scattering

The elastic backward proton-deuteron scattering is analyzed including relativistic effects in the deuteron and the mechanism of this reaction which includes the graphs corresponding to the emission, rescattering and absorption of the virtual pion by a deuteron nucleon in addition to the one-nucleon exchange graph. It allows one to obtain a rather satisfactory description of all the experimental data on the differential cross section, tensor analyzing power of the deuteron and transfer polarization in this reaction.

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

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.

preprint2010arXiv

Variational description of continuum states in terms of integral relations

Two integral relations derived from the Kohn Variational Principle (KVP) are used for describing scattering states. In usual applications the KVP requires the explicit form of the asymptotic behavior of the scattering wave function. This is not the case when the integral relations are applied since, due to their short range nature, the only condition for the scattering wave function $Ψ$ is that it be the solution of $(H-E)Ψ=0$ in the internal region. Several examples are analyzed for the computation of phase-shifts from bound state type wave functions or, in the case of the scattering of charged particles, it is possible to obtain phase-shifts using free asymptotic conditions. As a final example we discuss the use of the integral relations in the case of the Hyperspherical Adiabatic method.

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.

preprint2009arXiv

Integral relations for three-body continuum states with the adiabatic expansion

Application of the Hyperspherical Adiabatic expansion to describe three-body scattering states suffers the problem of a very slow convergence. Contrary to what happens for bound states, a huge number of hyperradial equations has to be solved, and even if done, the extraction of the scattering amplitude is problematic. In this paper we show how to obtain accurate scattering phase shifts using the Hyperspherical Adiabatic expansion. To this aim two integral relations, derived from the Kohn Variational Principle, are used. The convergence of this procedure is as fast as for bound states.

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.

preprint2006arXiv

Testing nuclear forces by polarization transfer coefficients in d(\vec p, \vec p)d and d(\vec p,\vec d)p reactions at E^{lab}_p = 22.7 MeV

The proton to proton polarization transfer coefficients K_x^{x'}, K_y^{y'}, K_z^{x'} and the proton to deuteron polarization transfer coefficients K_x^{x'}, K_y^{y'}, K_z^{x'}, K_x^{y'z'}, K_y^{z'z'}, K_z^{y'z'}, K_y^{x'z'} and K_y^{x'x'-y'y'} have been measured in d(\vec p, \vec p)d and d(\vec p, \vec d)p reactions at E^{lab}_p = 22.7 MeV, respectively. The data have been compared to predictions of modern nuclear forces obtained by solving the three-nucleon Faddeev equations in momentum space. Realistic (semi) phenomenological nucleon-nucleon potentials combined with model three-nucleon forces and modern chiral nuclear forces have been used. The AV18, CD Bonn, Nijm I and II nucleon-nucleon interactions have been applied alone or combined with the Tucson-Melbourne 99 three-nucleon force, adjusted separately for each potential to reproduce the triton binding energy. For the AV18 potential also the Urbana IX three-nucleon force have been used. In addition chiral NN potentials in the next-to-leading-order and chiral two- and three-nucleon forces in the next-to-next-to-leading-order have been applied. Only when three-nucleon forces are included a satisfactory description of all data results. For the chiral approach the restriction to the forces in the next-to-leading order is insufficient. Only when going over to the next-to-next-to-leading order one gets a satisfactory description of the data, similar to the one obtained with the (semi) phenomenological forces.

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

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.

preprint2005arXiv

Variational Calculation on A=3 and 4 Nuclei with Non-Local Potentials

The application of the hyperspherical harmonic approach to the case of non-local two-body potentials is described. Given the properties of the hyperspherical harmonic functions, there are no difficulties in considering the approach in both coordinate and momentum space. The binding energies and other ground state properties of A=3 and 4 nuclei are calculated using the CD Bonn 2000 and N3LO two-body potentials. The results are shown to be in excellent agreement with corresponding ones obtained by other accurate techniques.

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

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.

preprint2003arXiv

Recent Developments in Few-Nucleon Systems

N-d elastic scattering is studied at different energies using one of the modern NN interactions, the Argonne v_{18} which explicitly includes the magnetic moment interaction between two nucleons. This interaction, which has been often neglected in the description of the few-nucleon continuum, produces sizable modifications in some elastic observables. Its effects, as well as those produced by the Coulomb potential, are analyzed as a function of energy. The magnetic moment interaction produces appreciable effects in $p-d$ scattering at low energies butthey are very small above 10 MeV. Above 65 MeV Coulomb effects can be observed only in specific observables as for example $T_{21}$.

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

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

Benchmark Test Calculation of a Four-Nucleon Bound State

In the past, several efficient methods have been developed to solve the Schroedinger equation for four-nucleon bound states accurately. These are the Faddeev-Yakubovsky, the coupled-rearrangement-channel Gaussian-basis variational, the stochastic variational, the hyperspherical variational, the Green's function Monte Carlo, the no-core shell model and the effective interaction hyperspherical harmonic methods. In this article we compare the energy eigenvalue results and some wave function properties using the realistic AV8' NN interaction. The results of all schemes agree very well showing the high accuracy of our present ability to calculate the four-nucleon bound state.

preprint2000arXiv

Phase Shifts and Mixing Parameters for Elastic Proton-Deuteron Scattering

The eigenphase shifts and mixing parameters for elastic p-d scattering below the deuteron-breakup threshold are calculated for the AV14 two-body potential using two different numerical methods. The excellent agreement confirms that it is possible to perform accurate numerical studies of p-d scattering as well as n-d scattering for low energies. The numerical results can be considered as benchmarks for future calculations.

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.