Source author record

A. Kievsky

A. Kievsky 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

61works
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

61 published item(s)

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

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

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

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

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.

preprint2016arXiv

Matching universal behavior with potential models

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

preprint2015arXiv

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

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

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

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

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

preprint2015arXiv

Weakly bound states with spin-isospin symmetry

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

preprint2014arXiv

$N$-boson spectrum from a Discrete Scale Invariance

We present the analysis of the $N$-boson spectrum computed using a soft two-body potential the strength of which has been varied in order to cover an extended range of positive and negative values of the two-body scattering length $a$ close to the unitary limit. The spectrum shows a tree structure of two states, one shallow and one deep, attached to the ground-state of the system with one less particle. It is governed by an unique universal function, $Δ(ξ)$, already known in the case of three bosons. In the three-particle system the angle $ξ$, determined by the ratio of the two- and three-body binding energies $E_3/E_2=\tan^2ξ$, characterizes the Discrete Scale Invariance of the system. Extending the definition of the angle to the $N$-body system as $E_N/E_2=\tan^2ξ$, we study the $N$-boson spectrum in terms of this variable. The analysis of the results, obtained for up to $N=16$ bosons, allows us to extract a general formula for the energy levels of the system close to the unitary limit. Interestingly, a linear dependence of the universal function as a function of $N$ is observed at fixed values of $a$. We show that the finite-range nature of the calculations results in the range corrections that generate a shift of the linear relation between the scattering length $a$ and a particular form of the universal function. We also comment on the limits of applicability of the universal relations.

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

Efimov spectrum in bosonic systems with increasing number of particles

It is well-known that three-boson systems show the Efimov effect when the two-body scattering length $a$ is large with respect to the range of the two-body interaction. This effect is a manifestation of a discrete scaling invariance (DSI). In this work we study DSI in the $N$-body system by analysing the spectrum of $N$ identical bosons obtained with a pairwise gaussian interaction close to the unitary limit. We consider different universal ratios such as $E_N^0/E_3^0$ and $E_N^1/E_N^0$, with $E_N^i$ being the energy of the ground ($i=0$) and first-excited ($i=1$) state of the system, for $N\le16$. We discuss the extension of the Efimov radial law, derived by Efimov for $N=3$, to general $N$.

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

Universality and scaling in the $N$-body sector of Efimov physics

Universal behaviour has been found inside the window of Efimov physics for systems with $N=4,5,6$ particles. Efimov physics refers to the emergence of a number of three-body states in systems of identical bosons interacting {\it via} a short-range interaction becoming infinite at the verge of binding two particles. These Efimov states display a discrete scale invariance symmetry, with the scaling factor independent of the microscopic interaction. Their energies in the limit of zero-range interaction can be parametrized, as a function of the scattering length, by a universal function. We have found, using a particular form of finite-range scaling, that the same universal function can be used to parametrize the energies of $N\le6$ systems inside the Efimov-physics window. Moreover, we show that the same finite-scale analysis reconciles experimental measurements of three-body binding energies with the universal theory.

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

Recombination rates from potential models close to the unitary limit

We investigate universal behavior in the recombination rate of three bosons close to threshold. Using the He-He system as a reference, we solve the three-body Schrödinger equation above the dimer threshold for different potentials having large values of the two-body scattering length $a$. To this aim we use the hyperspherical adiabatic expansion and we extract the $S$-matrix through the integral relations recently derived. The results are compared to the universal form, $α\approx 67.1\sin^2[s_0\ln(κ_*a)+γ]$, for different values of $a$ and selected values of the three-body parameter $κ_*$. A good agreement with the universal formula is obtained after introducing a particular type of finite-range corrections, which have been recently proposed by two of the authors in Ref.[1]. Furthermore, we analyze the validity of the above formula in the description of a very different system: neutron-neutron-proton recombination. Our analysis confirms the universal character of the process in systems of very different scales having a large two-body scattering length.

preprint2013arXiv

Universal nature and finite-range corrections in elastic atom-dimer scattering below the dimer breakup threshold

We investigate universal behavior in elastic atom-dimer scattering below the dimer breakup threshold calculating the atom-dimer effective-range function $ak\cotδ$. Using the He-He system as a reference, we solve the Schrödinger equation for a family of potentials having different values of the two-body scattering length $a$ and we compare our results to the universal zero-range form deduced by Efimov, $ak\cotδ=c_1(ka)+c_2(ka)\cot[s_0\ln(κ_*a)+ϕ(ka)]$, for selected values of the three-body parameter $κ_*$. Using the parametrization of the universal functions $c_1,c_2,ϕ$ given in the literature, a good agreement with the universal formula is obtained after introducing a particular type of finite-range corrections. Furthermore, we show that the same parametrization describes a very different system: nucleon-deuteron scattering below the deuteron breakup threshold. Our analysis confirms the universal character of the process, and relates the pole energy in the effective-range function of nucleon-deuteron scattering to the three-body parameter $κ_*$.

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 helium trimer with soft-core potentials

The helium trimer is studied using two- and three-body soft-core potentials. Realistic helium-helium potentials present an extremely strong short-range repulsion and support a single, very shallow, bound state. The description of systems with more than two helium atoms is difficult due to the very large cancellation between kinetic and potential energy. We analyze the possibility of describing the three helium system in the ultracold regime using a gaussian representation of a widely used realistic potential, the LM2M2 interaction. However, in order to describe correctly the trimer ground state a three-body force has to be added to the gaussian interaction. With this potential model the two bound states of the trimer and the low energy scattering helium-dimer phase shifts obtained with the LM2M2 potential are well reproduced.

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

Analysis of three-nucleon forces effects in the $A=3$ system

Using modern nucleon-nucleon interactions in the description of the $A=3,4$ nuclear systems the $χ^2$ per datum results to be much bigger than one. In particular it is not possible to reproduce the three- and four-nucleon binding energies and the $n-d$ scattering length simultaneously. This is one manifestation of the necessity of including a three-nucleon force in the nuclear Hamiltonian. In this paper we perform an analysis of some, widely used, three-nucleon force models. We analyze their capability to describe the aforementioned quantities and, to improve their description, we propose modifications in the parametrization of the models. The effects of these new parametrization are studied in some polarization observables at low energies.

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

Recent Developments in Few-Nucleon Scattering

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. Several models of the three-nucleon force exist and are applied in the description of light nuclei. However, as it is discussed here, a simultaneous description of the three- and four-body binding energies and the $n-d$ doublet scattering length seems to be problematic. Accordingly, a comparative study of some of these models is performed. In a different analysis, we study applications of the Kohn Variational Principle, formulated in terms of integral relations, to describe $N-d$ scattering processes.

preprint2010arXiv

Relativistic description of 3He(e,e'p)2H

The Relativistic Distorted-Wave Impulse Approximation is used to describe the $^3$He($e,e^\prime p$)$^2$H process. We describe the $^3$He nucleus within the adiabatic hyperspherical expansion method with realistic nucleon-nucleon interactions. The overlap between the $^3$He and the deuteron wave functions can be accurately computed from a three-body calculation. The nucleons are described by solutions of the Dirac equation with scalar and vector (S-V) potentials. The wave function of the outgoing proton is obtained by solving the Dirac equation with a S-V optical potential fitted to elastic proton scattering data on the residual nucleus. Within this theoretical framework, we compute the cross section of the reaction and other observables like the transverse-longitudinal asymmetry, and compare them with the available experimental data measured at JLab.

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

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.

preprint2004arXiv

Variational Estimates using a Discrete Variable Representation

The advantage of using a Discrete Variable Representation (DVR) is that the Hamiltonian of two interacting particles can be constructed in a very simple form. However the DVR Hamiltonian is approximate and, as a consequence, the results cannot be considered as variational ones. We will show that the variational character of the results can be restored by performing a reduced number of integrals. In practice, for a variational description of the lowest n bound states only n(n+1)/2 integrals are necessary whereas D(D+1)/2 integrals are enough for the scattering states (D is the dimension of the S matrix). Applications of the method to the study of dimers of He, Ne and Ar, for both bound and scattering states, are presented.

preprint2003arXiv

Microscopic study of the He2-SF6 trimers

The He2-SF6 trimers, in their different He isotopic combinations, are studied both in the framework of the correlated Jastrow approach and of the Correlated Hyperspherical Harmonics expansion method. The energetics and structure of the He-SF6 dimers are analyzed, and the existence of a characteristic rotational band in the excitation spectrum is discussed, as well as the isotopic differences. The binding energies and the spatial properties of the trimers, in their ground and lowest lying excited states, obtained by the Jastrow ansatz are in excellent agreement with the results of the converged CHH expansion. The introduction of the He-He correlation makes all trimers bound by largely suppressing the short range He-He repulsion. The structural properties of the trimers are qualitatively explained in terms of the shape of the interactions, Pauli principle and masses of the constituents.

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.