Source author record

Nir Barnea

Nir Barnea 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

35works
12topics
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

35 published item(s)

preprint2022arXiv

Spectrum of light nuclei in a finite volume

Lattice quantum chromodynamics calculations of multi-baryon systems with physical quark masses would start a new age of ab initio predictions in nuclear physics. Performed on a finite grid, such calculations demand extrapolation of their finite volume numerical results to free-space physical quantities. Such extraction of the physical information can be carried out fitting effective field theories (EFTs) directly to the finite-volume results or utilizing the Lüscher free-space formula or its generalizations for extrapolating the lattice data to infinite volume. To understand better the effect of periodic boundary conditions on the binding energy of few nucleon systems we explore here light nuclei with physical masses in a finite box and in free space. The stochastic variational method is used to solve the few-body systems. Substantial optimizations of the method are introduced to enable efficient calculations in a periodic box. With the optimized code, we perform accurate calculations of light nuclei $A \le 4$ within leading order pionless EFT. Using Lüscher formula for the two-body system, and its generalization for 3- and 4-body systems, we examine the box effect and explore possible limitations of these formulas for the considered nuclear systems.

preprint2021arXiv

Implementation of local chiral interactions in the hyperspherical harmonics formalism

With the goal of using chiral interactions at various orders to explore properties of the few-body nuclear systems, we write the recently developed local chiral interactions as spherical irreducible tensors and implement them in the hyperspherical harmonics expansion method. We devote particular attention to three-body forces at next-to-next-to leading order, which play an important role in reproducing experimental data. We check our implementation by benchmarking the ground-state properties of $^3$H, $^3$He and $^4$He against the available Monte Carlo calculations. We then confirm their order-by-order truncation error estimates and further investigate uncertainties in the charge radii obtained by using the precise muonic atom data for single-nucleon radii. Having local chiral Hamiltonians at various orders implemented in our hyperspherical harmonics suites of codes opens up the possibility to test such interactions on other light-nuclei properties, such as electromagnetic reactions.

preprint2020arXiv

The collinear helium atom and two-electron ions

Collinear configurations of the helium-like atomic systems, relevant, e.g., for the quasifree mechanism of the double photoionization of helium, are studied, parameterized by the single scalar parameter $-1\leq λ\leq1$ ("collinear parameter") where $λ=0$ corresponds to the electron-nucleus ($\textbf{e-n}$) coalescence and $λ=1$ corresponds to the electron-electron ($\textbf{e-e}$) coalescence. In general, $λ>0$ corresponds to the \textbf{n-e-e} configuration, and $λ<0$ to the \textbf{e-n-e} configuration. Simple mathematical representations of the expectation values of the Dirac delta function relevant for the collinear configurations are derived and calculated from fully three-body dynamics without approximation for the two-electron atomic wave functions with nuclear charge $1\leq Z\leq5$. Simple formulas for calculating the expectation values of the kinetic and potential energy operators in collinear configurations are derived. Unusual physical properties of the \textbf{n-e-e} collinear configurations found for certain ranges of $λ$ are presented. The first few angular Fock coefficients for collinear configurations are derived as functions of $λ$. Highly accurate model wave functions describing the ground states of the two-electron atoms with collinear arrangement of the particles are constructed. All results are illustrated by tables and figures.

preprint2020arXiv

Universal Short Range Correlations in Bosonic Helium Clusters

Short-range correlations in bosonic Helium clusters, composed of $^4$He atoms, are studied utilizing the generalized contact formalism. The emergence of universal $n$-body short range correlations is formulated and demonstrated numerically via Monte Carlo simulations. The values of the $n$-particle contacts are evaluated for $n\le5$. In the thermodynamic limit, the two-body contact is extracted from available experimental measurements of the static structure factor of liquid $^4$He at high momenta, and found in a good agreement with the value extracted from our calculations.

preprint2019arXiv

Removing the Wigner bound in non-perturbative effective field theory

The Wigner bound, setting an upper limit on the scattering effective range, is examined at different orders of contact effective field theory. Using cutoff regulator we show that the bound loosens when higher orders of the theory are considered. For a sharp and a Gaussian regulators, we conjecture an analytic formula for the dependence of the Wigner bound on the theory's order. It follows that the bound vanishes in the limit of infinite order. Using a concrete numerical example we demonstrate that the above surmise still holds after renormalization at finite cutoff. Studying the 3-body system with this example, we have found that limiting the permissible range of cutoffs by the Wigner bound, we avoid the Thomas collapse, and don't need to promote the 3-body force to leading order.

preprint2016arXiv

Short range correlations - The important role of few-body dynamics in many-body systems

For many-body systems with short range interaction a series of relations were derived connecting many properties of the system to the dynamics of a closely packed few-body subsystems. Some of these relations were experimentally verified in ultra cold atomic gases. Here we shall review the implications of these developments on our understanding of nuclear one and two-body momentum distributions, and on the electron scattering Coulomb sum rule.

preprint2016arXiv

Subcomponents $ψ_{4,1}^{(2d)}$ and $ψ_{3,0}^{(2c)}$ of the angular Fock coefficients: Explicit analytic representation

The Fock expansion [1] describes the $S$-state wave function of the two-electron atomic system in the vicinity of the triple coalescence point. The present work constitutes the additional appendix to our paper [2] devoted to refinement and further development of calculation of the angular Fock coefficients of the Fock expansion. We derive the explicit analytic expressions for the complicated subcomponents $ψ_{4,1}^{(2d)}$ and $ψ_{3,0}^{(2c)}$ of the angular Fock coefficients. An unusual method of using the Wolfram Mathematica was applied.

preprint2016arXiv

The Generalized Nuclear Contact and its Application to the Photoabsorption Cross-Section

Using the zero-range model, it was demonstrated recently that Levinger's quasi-deuteron model can be utilized to extract the nuclear neutron-proton contact. Going beyond the zero-range approximation and considering the full nuclear contact formalism, we rederive here the quasi-deuteron model for the nuclear photoabsorption cross-section and utilize it to establish relations and constrains for the general contact matrix. We also define and demonstrate the importance of the diagonalized nuclear contacts, which can be also relevant to further applications of the nuclear contacts.

preprint2015arXiv

Angular Fock coefficients. Fixing the errors, and further development

The angular coefficients $ψ_{k,p}(α,θ)$ of the Fock expansion characterizing the S-state wave function of the two-electron atomic system, are calculated in hyperspherical angular coordinates $α$ and $θ$. To solve the problem the Fock recurrence relations separated into the independent individual equations associated with definite power $j$ of the nucleus charge $Z$, are applied. The "pure" $j$-components of the angular Fock coefficients, orthogonal to of the hyperspherical harmonics $Y_{kl}$, are found for even values of $k$. To this end, the specific coupling equation is proposed and applied. Effective techniques for solving the individual equations with simplest nonseparable and separable right-hand sides are proposed. Some mistakes/misprints made earlier in representations of $ψ_{2,0}$, were noted and corrected. All $j$-components of $ψ_{4,1}$ and the majority of components and subcomponents of $ψ_{3,0}$ are calculated and presented for the first time. All calculations were carried out with the help of the Wolfram \emph{Mathematica}.

preprint2015arXiv

Correlations in light nuclei and their relation to fine tuning and uncertainty quantifications of many body forces in low-energy nuclear physics

The large nucleon-nucleon scattering length, and the isospin approximate symmetry, are low energy properties of quantum chromodynamics (QCD). These entail correlations in the binding energies of light nuclei, e.g., the A=3 iso-multiplet, and Tjon's correlation between the binding energy of three and four body nuclei. Using a new representation of these, we establish that they translate into a correlation between different short-range contributions to three body forces in chiral effective field theory of low-energy nuclear physics. We demonstrate that these correlations should be taken into account in order to avoid fine-tuning in the calibration of three body forces. We relate this to the role of correlations in uncertainty quantification of non-renormalizable effective field theories of the nuclear regime. In addition, we show that correlations can be useful in assessing the importance of forces induced by renormalization group (RG) transformations. We give numerical evidence that such RG transformations can be represented effectively by adding a constant to the pure three nucleon contact low energy constant $c_E$.

preprint2015arXiv

Electric dipole polarizability: from few- to many-body systems

We review the Lorentz integral transform coupled-cluster method for the calculation of the electric dipole polarizability. We benchmark our results with exact hyperspherical harmonics calculations for 4He and then we move to a heavier nucleus studying 16O. We observe that the implemented chiral nucleon-nucleon interaction at next-to-next-to-next-to-leading order underestimates the electric dipole polarizability.

preprint2015arXiv

Generalized nuclear contacts and the nucleon's momentum distributions

The general nuclear contact matrices are defined, taking into consideration all partial waves and finite-range interactions, extending Tan's work for the zero range model. The properties of these matrices are discussed and the relations between the contacts and the one-nucleon and two-nucleon momentum distributions are derived. Using these relations, a new asymptotic connection between the one-nucleon and two-nucleon momentum distributions, describing the two-body short-range correlations in nuclei, is obtained. Using available numerical data, we extract few connections between the different contacts and verify their relations to the momentum distributions. The numerical data also allows us to identify the main nucleon momentum range affected by two-body short-range correlations. Utilizing these relations and the numerical data, we also verify a previous independent prediction connecting between the Levinger constant and the contacts. This work provides an important indication for the relevance of the contact formalism to nuclear systems, and should open the path for revealing more useful relations between the contacts and interesting quantities of nuclei and nuclear matter.

preprint2015arXiv

Photoassociation of Universal Efimov Trimers

In view of recent experiments in ultracold atomic systems, the photoassociation of Efimov trimers, composed of three identical bosons, is studied utilizing the multipole expansion. We study both the normal hierarchy case, where one-body current is dominant, and the strong hierarchy case, relevant for photoassociation in ultracold atoms, where two-body current is dominant. For identical particles in the normal hierarchy case, the leading contribution comes from the r^2 s-mode operator and from the quadrupole d-mode operator. The s-mode reaction is found to be dominant at low temperature, while as the temperature increases the d-mode becomes as significant. For the strong hierarchy case, the leading contribution comes from a 2-body s-wave δ operator. In both cases log periodic oscillations are found in the cross section. For large but finite scattering length the amplitude of the oscillations becomes larger in comparison to infinite scattering length case. We apply our theory to photoassociation of 7Li ultracold atoms and show a good fit to the available experimental results.

preprint2015arXiv

Spectra and Scattering of Light Lattice Nuclei from Effective Field Theory

An effective field theory is used to describe light nuclei, calculated from quantum chromodynamics on a lattice at unphysically large pion masses. The theory is calibrated at leading order to two available data sets on two- and three-body nuclei for two pion masses. At those pion masses we predict the quartet and doublet neutron-deuteron scattering lengths, and the alpha-particle binding energy. For $m_π=510~$MeV we obtain, respectively, $^4a_{\rm nD}=2.3\pm 1.3~$fm, $^2a_{\rm nD}=2.2\pm 2.1~$fm, and $B_α^{}=35\pm 22~$MeV, while for $m_π=805~$MeV $^4a_{\rm nD}=1.6\pm 1.3~$fm, $^2a_{\rm nD}=0.62\pm 1.0~$fm, and $B_α^{}=94\pm 45~$MeV are found. Phillips- and Tjon-like correlations to the triton binding energy are established. Higher-order effects on the respective correlation bands are found insensitive to the pion mass. As a benchmark, we present results for the physical pion mass, using experimental two-body scattering lengths and the triton binding energy as input. Hints of subtle changes in the structure of the triton and alpha particle are discussed.

preprint2015arXiv

Understanding the proton radius puzzle: Nuclear structure effects in light muonic atoms

We present calculations of nuclear structure effects to the Lamb shift in light muonic atoms. We adopt a modern ab-initio approach by combining state-of-the-art nuclear potentials with the hyperspherical harmonics method. Our calculations are instrumental to the determination of nuclear charge radii in the Lamb shift measurements, which will shed light on the proton radius puzzle.

preprint2014arXiv

An efficient method for evaluating energy-dependent sum rules

Energy-dependent sum rules are useful tools in many fields of physics. In nuclear physics, they typically involve an integration of the response function over the nuclear spectrum with a weight function composed of integer powers of the energy. More complicated weight functions are also encountered, e.g., in nuclear polarization corrections of atomic spectra. Using the Lorentz integral transform method and the Lanczos algorithm, we derive a computationally efficient technique for evaluating such sum rules that avoids the explicit calculation of both the continuum states and the response function itself. Our numerical results for electric dipole sum rules of the Helium-4 nucleus with various energy-dependent weights show rapid convergence with respect to the number of Lanczos steps. This demonstrates the usefulness of the method in a variety of electroweak reactions.

preprint2014arXiv

Analytic matrix elements for the two-electron atomic basis with logarithmic terms

The two-electron problem for the helium-like atom/ions in $S$-state is considered. The basis containing the integer powers of $\ln r$, where $r$ is a radial variable of the Fock expansion, is studied. In this basis, the analytic expressions for the matrix elements of the corresponding Hamiltonian are presented. These expressions include only special functions presented by the built-in $Mathematica$ codes, what enables very fast and accurate computation of the matrix elements. The decisive contribution of the correct logarithmic terms to the behavior of the two-electron wave function in the vicinity of the triple-coalescence point is reaffirmed.

preprint2014arXiv

Experimental evaluation of the nuclear neutron-proton contact

The nuclear neutron-proton contact is introduced, generalizing Tan's work, and evaluated from medium energy nuclear photodisintegration experiments. To this end we reformulate the quasi-deuteron model of nuclear photodisintegration and establish the bridge between the Levinger constant and the contact. Using experimental evaluations of Levinger's constant we extract the value of the neutron-proton contact in finite nuclei and in symmetric nuclear matter. Assuming isospin symmetry we propose to evaluate the neutron-neutron contact through measurement of photonuclear spin correlated neutron-proton pairs.

preprint2014arXiv

Improved estimates of the nuclear structure corrections in $μ$D

We calculate the nuclear structure corrections to the Lamb shift in muonic deuterium by using state-of-the-art nucleon-nucleon potentials derived from chiral effective field theory. Our calculations complement previous theoretical work obtained from phenomenological potentials and the zero range approximation. The study of the chiral convergence order-by-order and the dependence on cutoff variations allows us to improve the estimates on the nuclear structure corrections and the theoretical uncertainty coming from nuclear potentials. This will enter the determination of the nuclear radius from ongoing muonic deuterium experiments at PSI.

preprint2014arXiv

In-medium antikaon and eta-meson interactions and bound states

The role played by subthreshold meson-baryon dynamics is demonstrated in kaonic-atom, Kbar-nuclear and eta-nuclear bound-state calculations within in-medium models of Kbar-N and eta-N interactions. New analyses of kaonic atom data reveal appreciable multi-nucleon contributions. Calculations of eta-nuclear bound states show, in particular, that the eta-N scattering length is not a useful indicator of whether or not eta mesons bind in nuclei nor of the widths anticipated for such states.

preprint2014arXiv

Is the first excited state of the $α$-particle a breathing mode?

The isoscalar monopole excitation of 4He is studied within a few-body ab initio approach. We consider the transition density to the low-lying and narrow 0+ resonance, as well as various sum rules and the strength energy distribution itself at different momentum transfers q. Realistic nuclear forces of chiral and phenomenological nature are employed. Various indications for a collective breathing mode are found: i) the specific shape of the transition density, ii) the high degree of exhaustion of the non-energy-weighted sum rule at low q and iii) the complete dominance of the resonance peak in the excitation spectrum. For the incompressibility K of the alpha-particle values between 20 and 30 MeV are found.

preprint2014arXiv

Photo Reactions with Universal Trimers

Considering one-body and two-body currents, we study the photoassociation and photodissociation of universal bosonic trimers. Analyzing the relative importance of these currents we identify two physical scenarios (i) Normal hierarchy, where naive power counting holds and the one-body current dominates, and (ii) Strong hierarchy, where the one-body current is suppressed. For both scenarios we observe that at the high frequency tail, the response function exhibits log periodic oscillations in transition to or from any continuum state regardless of the reaction partial wave channel. In contrast, near threshold log periodic oscillations appear only in the leading $s$-wave components. These oscillations are the fingerprints of universal Efimov physics. We discuss the relevance of this effect to contemporary experiments in ultracold atoms.

preprint2013arXiv

First principles description of the giant dipole resonance in 16O

We present an ab-initio calculation of the giant dipole resonance in 16O based on a nucleon-nucleon (NN) interaction from chiral effective field theory that reproduces NN scattering data with high accuracy. By merging the Lorentz integral transform and the coupled-cluster methods, we extend the previous theoretical limits for break-up observables in light nuclei with mass numbers (A<=7), and address the collective giant dipole resonance of 16O. We successfully benchmark the new approach against virtually exact results from the hyper-spherical harmonics method in 4He. Our results for 16O reproduce the position and the total strength (bremsstrahlung sum rule) of the dipole response very well. When compared to the cross section from photo-absorption experiments the theoretical curve exhibits a smeared form of the peak. The tail region between 40 and 100 MeV is reproduced within uncertainties.

preprint2013arXiv

Log-Periodic Oscillations in the Photo Response of Efimov Trimers

The photoassociation of Efimov trimer, composed of three identical bosons, is studied utilizing the multipole expansion. For identical particles the leading contribution comes from the r^2 s-mode operator and from the quadrupole d-mode operator. Log-periodic oscillations are found in the photoassociation response function, both near the energy threshold for the leading s-wave reaction, and in the high frequency tail for all partial waves.

preprint2013arXiv

Nuclear polarization corrections in the mu-4He+ Lamb shift

Stimulated by the proton radius conundrum, measurements of the Lamb shift in various light muonic atoms are planned at PSI. The aim is to extract the rms charge radius with high precision, limited by the uncertainty in the nuclear polarization corrections. We present an ab-initio calculation of the nuclear polarization for mu-4He+ leading to an energy correction in the 2S-2P transitions of $δ_{pol}=-2.47$ meV $\pm$ 6%. We use two different state-of-the-art nuclear Hamiltonians and utilize the Lorentz integral transform with hyperspherical harmonics expansion as few-body methods. We take into account the leading multipole contributions, plus Coulomb, relativistic and finite-nucleon-size corrections. Our main source of uncertainty is the nuclear Hamiltonian, which currently limits the attainable accuracy. Our predictions considerably reduce the uncertainty with respect to previous estimates and should timely serve the mu-4He+ experiment planned for 2013.

preprint2013arXiv

The nuclear electric polarizability of 6He

We present an estimate of the nuclear electric polarizability of the 6He halo nucleus based on six-body microscopic calculations. Wave functions are obtained from semi-realistic two-body interactions using the hyperspherical harmonics expansion method. The polarizability is calculated as a sum rule of the dipole response function using the Lanczos algorithm and also by integrating the photo-absorption cross section calculated via the Lorentz integral transform method. We obtain alpha_E=1.00(14) fm^3, which is much smaller than the published value 1.99(40) fm^3 extracted from experimental data. This points towards a potential disagreement between microscopic theories and experimental observations.

preprint2013arXiv

Theoretical study of the ${}^{4}\rm{He}(γ,p)^3\rm{H}$ and ${}^{4}\rm{He}(γ,n)^3\rm{He}$ reactions

{\it Ab initio} calculation of the total cross section for the reactions $^{4}\rm{He}(γ,p)^3\rm{H}$ and $^{4}\rm{He}(γ,n)^3\rm{He}$ is presented, using state-of-the-art nuclear forces. The Lorentz integral transform (LIT) method is applied, which allows exact treatment of the final state interaction (FSI). The dynamic equations are solved using the effective interaction hyperspherical harmonics (EIHH) method. In this calculation of the cross sections the three-nucleon force is fully taken into account, except in the source term of the LIT equation for the FSI transition matrix element.

preprint2012arXiv

Multipole Analysis of Radio-Frequency Reactions in Ultra-Cold Atoms

Using the multipole expansion we analyze photo induced reactions in an ultra-cold atomic gas composed of identical neutral bosons. While the Frank-Condon factor dominates the photo induced spin-flip reactions, we have found that for frozen-spin process where the atomic spins are conserved the reaction rate is governed by the monopole $r^2$ and the quadrupole terms. Consequently, the dependence of the frozen-spin reaction rate on the photon wave number $k$ acquires an extra $k^4$ factor in comparison to the spin-flip process. Comparing the relative strength of the $r^2$ and quadrupole modes in dimer photoassociation we predict that the mutual importance of these two modes changes with temperature and scattering length.

preprint2012arXiv

Nonsymmetrized hyperspherical harmonics with realistic NN potentials

The Schroedinger equation is solved for an A-nucleon system using an expansion of the wave function in nonsymmetrized hyperspherical harmonics. Our approach is both an extension and a modification of the formalism developed by Gattobigio et al.. The extension consists in the inclusion of spin and isospin degrees of freedom such that a calculation with more realistic NN potential models becomes possible, whereas the modification allows a much simpler determination of the fermionic ground state. The approach is applied to four- and six-body nuclei (4He, 6Li) with various NN potential models. It is shown that the results for ground-state energy and radius agree well with those from the literature.

preprint2011arXiv

Gapless Superfluid in $SU(2n+1)$ Fermions Systems

We investigate the generalized Hubbard model of $(2n+1)$ Fermion species interacting via a symmetric contact attraction potential. We prove that the ground state of such system is a gapless superfluid, where a full Fermi surface coexists with a superfluid. Moreover, doing so we prove the existence of a free mode in a strongly interacting system, regardless of the potential strength. This proof holds at the mean field level. A Grassmannian Gaussian integration technique is used to deal with the problem. Our predictions may be relevant to future high-spin cold atoms experiments.

preprint2011arXiv

Transition states and the critical parameters of central potentials

Transition states or quantum states of zero energy appear at the boundary between the discrete part of the spectrum of negative energies and the continuum part of positive energy states. As such, transition states can be regarded as a limiting case of a bound state with vanishing binding energy, emerging for a particular set of critical potential parameters. In this work we study the properties of these critical parameters for short range central potentials. To this end we develop two exact methods and also utilize the first and second order WKB approximations. Using these methods we have calculated the critical parameters for several widely used central potentials. The general analytic expressions for the asymptotic representations of the critical parameters were derived for cases where either the orbital quantum number $l$ or the number $n$ of bound states approaches infinity. The above mathematical models enable us to answer the following physical (quantum mechanical) questions: i) what is the number of bound states for a given central potential and given orbital quantum number $l$; ii) what is the maximum value of $l$ which can provide a bound state for the given central potential; iii) what is the order of energy levels for the given form of the central potential. It is revealed that the ordering of energy levels depends on the potential singularity at the origin.

preprint2007arXiv

Inclusive electron scattering off 4He

Inclusive electron scattering off 4He is investigated for low and medium energy and momentum transfers. The final state interaction, given by the simple semirealistic Malfliet-Tjon potential, is treated rigorously applying the Lorentz Integral Transform (LIT) method. Besides the nonrelativistic one-body current a consistent meson exchange current is constructed and implemented. Results are presented for both longitudinal and transverse response functions at various momentum transfers. Good agreement with experimental data is found for the longitudinal response function, while some strength is missing in the transverse response function on the low-energy side of the quasi-elastic peak.

preprint2004arXiv

Incorporation of three-nucleon force in the effective interaction hyperspherical harmonic approach

It is shown how a bare three-nucleon force is incorporated into the formalism of the effective interaction approach for hyperspherical harmonics. As a practical example we calculate the ground state properties of 3H and 3He using the Argonne V18 nucleon-nucleon potential and the Urbana IX three-nucleon force. A very good convergence of binding energies and matter radii is obtained. We also find a very good agreement of our results compared to other high precision calculations.