Hyperspherical Harmonics Expansion Research Articles

Few-body model approach to the lowest bound S-state of non-symmetric exotic atoms and ions

Lowest bound S-state energy of Coulombic three-body systems (N Z+ μ − e −) consisting of a positively charged nucleus of charge number Z (N Z+), a negatively charged muon (μ −) and an electron (e −), is investigated in the framework of few-body (i.e., two- and three-body) cluster model approach. For the three-body cluster model, we adopted the hyperspherical harmonics expansion (HHE) method. An approximated two-body model calculation is also performed for all the three-body systems considered here. A Yukawa-type screened Coulomb potential with an arbitrary screening parameter (λ) is chosen for the two-body subsystems of the three-body system. In the resulting Schrödinger equation (SE), the three-body relative wave function is expanded in the complete set of hyperspherical harmonics (HH). The use of the orthonormality of HH in the SE leads to a set of coupled differential equations (CDEs) which are solved numerically for a manageable basis size to get the energy (E). The pattern of convergence in energy relative to increasing basis size is also investigated. Results are compared with some of those found in the literature.

Hyperspherical Cluster Model for Bosons: Application to Sub-threshold Halo States in Helium Drops

To describe long-range behaviour of one particle removed from a few- or a many-body system, a hyperspherical cluster model has been developed. It has been applied to the ground and first excited states of helium drops with five, six, eight and ten atoms interacting via a two-body soft gaussian potential. Convergence of the hyperspherical cluster harmonics expansion is studied for binding energies, root-mean-squared radii and overlaps of the wave functions of two helium drops differing by one atom. It was shown that with increasing model space the functional form of such overlaps at large distances converges to the correct asymptotic behaviour. The asymptotic normalization coefficients that quantify the overlaps’ amplitudes in this region are calculated. It was also shown that in the first excited state one helium atom stays far apart from the rest forming a two-body molecule, or a halo. The probability of finding the halo atom in the classically-forbidden region of space depends on the definition of the latter and on the valence atom binding energy. The total norm of the overlap integrals, the spectroscopic factor, represents the number of partitions of a many-body state into a chosen state of the system with one particle removed. The spectroscopic factors have been calculated and their sum rules are discussed giving a further insight into the structure of helium drops.

Study of resonances in exotic neutron-rich 14Be by the use of isospectral potential

In this work, Hyperspherical Harmonic Expansion (HHE) formalism aided by supersymmetric quantum mechanics is used to study bound and resonant states of 14Be in the three-body (12Be+ n + n) cluster model, and the analysis of the resonant states close to the binding threshold utilizes supersymmetric quantum mechanics (SSQM). GPT nucleon-nucleon potential together with SBB core-nucleon potential is chosen for the solution of three-body Schrödinger equation to get the lowest bound state energy and wave function. In the next stage, energy and wave function of the bound state is used to derive an isospectral potential that exhibit deep well following a strong barrier facilitating confinement (or trapping) of particle inside it at some positive energy (E>0). The trapping probability when plotted against energy shows a prominent peak at the energy of resonance. WKB approximation is used to determine the width of resonance width. Calculated resonance energy and width of resonance are compared with those measured experimentally for the promising neutron halo candidate, the 14Be.

Onset of the bound states in strange Λ, Ξ-hypernuclei and hyperspherical three-body model calculation for the corresponding double hyperon hypernuclei

In this work, we have explored ground state properties of Be hypernucleus in the light of recent observation of the experiment J-PARC E07 and the onset of stability for Be hypernucleus in the three-body framework using Hyperspherical Harmonics Expansion (HHE) method. An effective core-Y (Y= Λ, Ξ) Woods-Saxon potential with adjustable depth parameter is used to reproduce the observed ground state energy of the two-body subsystem by solving the two-body Schrödinger equation. The core-Y potential so obtained is clubbed with the YY potential model of the Nijmegen group to obtain the observables of the double-hyperon hypernuclei. Computed YY separation energy and YY bond energy are compared with the data available in the literature.

Three-body Faddeev calculations for and hypernuclei* *M.R.H. was Supported by the National Science Foundation under Grant No. NSF-PHY-2000029 with Central State University

We study the ground-state properties of the double hyperon for and nuclei in a three-body model . We solve two coupled Faddeev equations corresponding to the three-body configurations and in configuration space with the hyperspherical harmonics expansion method by employing the most recent hyperon-hyperon interactions obtained from lattice QCD simulations. Our numerical analysis for , using three lattice interaction models, leads to a ground state binding energy in the MeV domain and the separations and in the domains of fm and fm, respectively. The binding energy of the double- hypernucleus leads to MeV and consequently to smaller separations fm and fm. In addition to geometrical properties, we study the structure of ground-state wave functions and show that the main contributions are from the wave channels. Our results are consistent with the existing theoretical and experimental data.

Hyperspherical three-body model calculation for the bound and resonant states of the neutron dripline nuclei 42,44Mg using isospectral potential

In this paper, an effective theoretical scheme is presented to investigate ground and resonant levels of weakly bound nuclei near the neutron-dripline. The hyperspherical harmonics expansion method is employed for the bound state while supersymmetric quantum mechanics (SSQM) is used for the resonant state. The three-body Schrödinger equation is solved for the lowest bound state with chosen set of two-body potentials to obtain the energy and normalized wavefunction. The three-body effective potential has a very shallow well following a low-wide barrier giving a broad resonance due to far extended wavefunction in the asymptotic region. Numerical computation of such broad resonance is very challenging due to large computational errors. Hence, a one-parameter family of isospectral potential is constructed with the ground state energy and wavefunction, by the use of the algebra of SSQM. This enhanced potential facilitates a more accurate computation of the resonance energy and width. Effectiveness of the scheme has been checked for the J=π0+ bound and resonant states of neutron-rich nuclei 42,44Mg in the coreNN three-body cluster model where core=40,42Mg. In the present calculation, two-neutron separation energies obtained with very large number of partial waves for the J=π0+ bound states of neutron-rich nuclei 42,44Mg are found to be 0.41433 MeV and 0.44975 MeV respectively. Results of the calculation have been compared with those found in the literature.

Implementation of Local Chiral Interactions in the Hyperspherical Harmonics Formalism

With the goal of using chiral interactions at various orders to explore the 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 3H, 3He, and 4He 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.

Hyperspherical Harmonics Expansion on Lagrange Meshes for Bosonic Systems in One Dimension

A one-dimensional system of bosons interacting with contact and single-Gaussian forces is studied with an expansion in hyperspherical harmonics. The hyperradial potentials are calculated using the link between the hyperspherical harmonics and the single-particle harmonic-oscillator basis while the coupled hyperradial equations are solved with the Lagrange-mesh method. Extensions of this method are proposed to achieve good convergence with small numbers of mesh points for any truncation of hypermomentum. The convergence with hypermomentum strongly depends on the range of the two-body forces: it is very good for large ranges but deteriorates as the range decreases, being the worst for the contact interaction. In all cases, the lowest-order energy is within 4.5% of the exact solution and shows the correct cubic asymptotic behaviour at large boson numbers. Details of the convergence studies are presented for 3, 5, 20 and 100 bosons. A special treatment for three bosons was found to be necessary. For single-Gaussian interactions, the convergence rate improves with increasing boson number, similar to what happens in the case of three-dimensional systems of bosons.

Convergence of the hyperspherical-harmonics expansion with increasing number of particles for bosonic systems. II. Inclusion of the three-body force

This paper presents a numerical convergence study of a hyperspherical-harmonics expansion for binding energies of a system of $4\ensuremath{\le}N\ensuremath{\le}728$ helium atoms using a phenomenological soft attractive two-body He-He potential and a repulsive three-body force aimed at compensating for the absence of the two-body repulsive core. Earlier calculations with such a potential have shown an improved convergence when $N$ increases from four to six. The present study reveals that the improved convergence occurs only for a limited range of $N$ determined by the range of the three-body repulsion. For a soft repulsive three-body force, the convergence is fast for $N\ensuremath{\le}20$, while for a short-range three-body repulsion it deteriorates at $N\ensuremath{\ge}10$. The reasons for this deterioration are discussed. The range of the three-body force also determines the binding energy behavior with $N$, and it is also responsible for binding the excited states. The long-range force binds all first excited ${0}^{+}$ states but strongly underbinds the systems of $N$ helium atoms at large $N$. The short-range force does not bind the first ${0}^{+}$ states for $A\ensuremath{\le}7$ but gives better predictions of binding energies as compared to the calculations of other authors though overestimating them. Some options to improve both the description of the binding energies and the convergence of the hyperspherical-harmonics expansion using phenomenological forces are discussed. It is pointed out that a fast convergence is very much needed for the reliable predictions of states with nonzero angular momentum, examples of which are also given.

Solution of Faddeev integral equations in configuration space using the hyperspherical harmonics expansion method

In this paper, a method has been developed to solve three-particle Faddeev equations in the configuration space making use of a series expansion in hyperspherical harmonics. The following parameters of the bound state of triton and helium-3 nuclei have been calculated: the binding energies, the weights of symmetric and mixed-symmetry components of the wave function, the magnetic moments, and the charge radii.

Application of hyperspherical harmonics expansion method to the low-lying bound S-states of exotic two-muon three-body systems

In this paper, energies of the low-lying bound S-states (L = 0) of exotic three-body systems, consisting a nuclear core of charge +Ze (Z being atomic number of the core) and two negatively charged valence muons, have been calculated by hyperspherical harmonics expansion method (HHEM). The three-body Schrödinger equation is solved assuming purely Coulomb interaction among the binary pairs of the three-body systems X Z+μ-μ- for Z = 1 to 54. Convergence pattern of the energies have been checked with respect to the increasing number of partial waves Λmax. For available computer facilities, calculations are feasible up to Λmax = 28 partial waves, however, calculation for still higher partial waves have been achieved through an appropriate extrapolation scheme. The dependence of bound state energies has been checked against increasing nuclear charge Z and finally, the calculated energies have been compared with the ones of the literature.

Low-Lying S-States of Two-Electron Systems

The energies of the low-lying bound S-states of some two-electron systems (treating them as three-body systems) like negatively charged hydrogen, neutral helium, positively charged-lithium, beryllium, carbon, oxygen, neon, argon and negatively charged muonium and exotic positronium ions have been calculated employing hyperspherical harmonics expansion method. The matrix elements of two-body interactions involve Raynal–Revai coefficients which are particularly essential for the numerical solution of three-body Schrődinger equation when the two-body potentials are other from Coulomb or harmonic. The technique has been applied for to two-electron ions 1H− (Z = 1) to 40Ar16+ (Z = 18), negatively charged-muonium Mu− and exotic positronium ion Ps−(e+e−e−) considering purely Coulomb interaction. The available computer facility restricted reliable calculations up to 28 partial waves (i.e. Km = 28) and energies for higher Km have been obtained by applying an extrapolation scheme suggested by Schneider.

Four-nucleon Scattering

Recent developments in the theoretical description of the four-nucleon scattering are discussed. Accurate results for neutron-3H and proton-3He scattering with realistic interactions below inelastic threshold are obtained using the variational hyperspherical harmonics expansion, and exact scattering equations in coordinate or momentum space. The momentum-space framework is successfully applied also to coupled proton-3H, neutron-3He and deuteron-deuteron reactions below three-cluster threshold as well as to neutron-3H and proton-3He scattering above the four-nucleon breakup threshold.

Convergence of the hyperspherical harmonic expansion for crystallographic texture

Advances in instrumentation allow a material texture to be measured as a collection of spatially resolved crystallite orientations rather than as a collection of pole figures. However, the hyperspherical harmonic expansion of a collection of spatially resolved crystallite orientations is subject to significant truncation error, resulting in ringing artifacts (spurious oscillations around sharp transitions) and false peaks in the orientation distribution function. This article finds that the ringing artifacts and the accompanying regions of negative probability density may be mitigated or removed entirely by modifying the coefficients of the hyperspherical harmonic expansion by a simple multiplicative factor. An addition theorem for the hyperspherical harmonics is derived as an intermediate result.

Nuclear Polarization Corrections to theμHe+4Lamb 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 $\delta_{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.

On a generalization of the generating function for Gegenbauer polynomials

A generalization of the generating function for Gegenbauer polynomials is introduced whose coefficients are given in terms of associated Legendre functions of the second kind. We discuss how our expansion represents a generalization of several previously derived formulae such as Heine's formula and Heine's reciprocal square-root identity. We also show how this expansion can be used to compute hyperspherical harmonic expansions for power-law fundamental solutions of the polyharmonic equation.

Nuclear electric polarizability of6He

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.

Convergence of the hyperspherical-harmonics expansion with increasing number of particles for bosonic systems

A formal proof is given that for $N$-boson systems in which two-body interaction potentials are described by a single Gaussian the ratio of higher-order to the lowest-order hyper-radial potentials decreases at $N\ensuremath{\rightarrow}\ensuremath{\infty}$ as ${N}^{\ensuremath{-}1/2}$ or faster. As a result, for such potentials, the convergence of binding energies for ground and several lowest excited states, obtained in expansion of the $N$-body wave function over the hyperspherical-harmonics basis, improves with increasing number of bosons. For a phenomenological three-body repulsive potential, introduced to account for the missing hard core, the ratio of higher-order to lowest-order hyper-radial potential, corresponding to this three-body potential, also decreases as ${N}^{\ensuremath{-}1/2}$ or faster when $N\ensuremath{\rightarrow}\ensuremath{\infty}$. Although adding the three-body contributions leads to increased influence from the total nondiagonal couplings around the node of the lowest hyper-radial potential, the arguments are given that this should not dramatically deteriorate the convergence if the range of repulsion is properly chosen. This means that the hyperspherical-harmonics expansion with soft two-body and repulsive three-body effective forces may become an attractive tool for studying the spectra of many-body systems. It is suggested that fine tuning of the three-body repulsion to reproduce the binding energies over a large region of $N$ is possible. In particular, it has been shown that an $N$-independent choice for the three-body repulsion exists for which the ground-state binding energies of $N\ensuremath{\le}112$ atoms of helium, obtained in the lowest-order approximation of the hyperspherical-harmonics expansion, are close to the prediction of the Green's function Monte Carlo method with a hard-core He-He potential.

Matter and charge radius of6He in the hyperspherical-harmonics approach

We present ab-initio calculations of the binding energy and radii of the two-neutron halo nucleus 6He using two-body low-momentum interactions based on chiral effective field theory potentials. Calculations are performed via a hyperspherical harmonics expansion where the convergence is sped up introducing an effective interaction for non-local potentials. The latter is essential to reach a satisfactory convergence of the extended matter radius and of the point-proton radius. The dependence of the results on the resolution scale is studied. A correlation is found between the radii and the two-neutron separation energy. The importance of three-nucleon forces is pointed out comparing our results and previous calculations to experiment.

The Quadrupole Response of Borromean Bosonic Trimers

The inelastic response of a Borromean trimer composed of three identical bosons is explored. To this end we use the quadrupole approximation as the excitation mechanism of the bosonic system from its ground state. We utilize the hyperspherical-harmonics expansion to solve the Schroedinger equation and the Lorentz integral transform method to calculate the reaction. It is found that the magnitude of the response function and corresponding sum rules increase as a power law when approaching the three-body threshold. It is also found that this increase is governed by unnatural exponents.