Electric Multipole Polarizabilities of Quantum Bound Systems in the Transition Matrix Formalism

A new general formalism for determining the electric multipole polarizabilities of quantum (atomic and nuclear) bound systems based on the use of the transition matrix in momentum space has been developed. As distinct from the conventional approach with the application of the spectral expansion of the total Green’s function, our approach does not require preliminary determination of the entire unperturbated spectrum; instead, it makes possible to calculate the polarizability of a few-body bound complex directly based on solving integral equations for the wave function of the ground bound state and the transition matrix at negative energy, both of them being real functions of momenta. A formula for the multipole polarizabilities of a two-body bound complex formed by a central interaction potential has been derived and studied. To test, the developed t-matrix formalism has been applied to the calculation of the dipole, quadrupole and octupole polarizabilities of the hydrogen atom.


Introduction
Study of the few-body quantum systems persists a perspective line of the advancement of the modern physics.In atomic physics, the Efimov effect for the three-body system, the existence of which was predicted in 1970 [1], has been experimentally verified using ultracold atoms and tuning the atom-atom scattering length by the magnetic field near a Feshbach resonance [2] (see, for example, the reviews [3][4][5]).In nuclear physics, extensive investigations of few-body nuclei advantageously undertaken over the past fifty years make it possible to gain new important information about the nuclear force.
The investigation of the behaviour of the few-body nuclei in the external electromagnetic field permits to obtain additional data on their properties, specifically, on the electric polarizabilities and magnetic susceptibilities as important fundamental quantities of complex systems.In spite of the fact that the study of the electric polarizabilities of the few-body nuclei has attracted considerable interest from both the experimental [6][7][8][9][10][11] and theoretical [11][12][13][14][15][16][17][18][19][20][21][22][23][24] points of view and definite progress in the field has been reached, a great deal needs to be done in this area.
Up to now, there exists a discrepancy between the direct experimental result for the electric dipole polarizability of the nucleus 3 He (determined by measuring deviations from the Rutherford scattering law of the low-energy elastic 3 He scattering by the Coulomb field of the heavy nucleus 208 Pb [8]) and the result deduced from the data for the total 3 He photoabsorption cross section using the sum rule 2 [9].

 
The conventional formula for the electric dipole polarizability of the -particle bound complex has the form where 1 11 D  r E is the dipole moment operator,  and   are the energy and the wave function of the v-excited state, 1 and 0 r ρ are the radius vectors of the particle 1 relative to the center of mass of the bound complex and the center of mass of the complex relative to the charged particle 0 creating the electric field, and the summation is taken over all possible excited discrete bound and continuum states.
Although the Formula (1) using the spectral expansion can be practically applied in the case of the two-body complex   2 N  , the direct determination of V. F. KHARCHENKO 100 introducing an additional function F , where F being the first-order correction to the unperturbated state satisfies an inhomogeneous differential equation.Dalgarno-Lewis method has been effectively applied in [12] and [22] to calculate the electric polarizabilities of hydrogen and helium isotopes.Though the result for 1 E obtained therewith supports the value deduced using the sum rule 2   3


  [9], the subject of the discrepancy of the experimental data for 3 He under discussion may not be considered as conclusively established since indispensable consistent calculations of the electric dipole polarizabilities of the three-body nuclei on the basis of the rigorous mathematical Faddeev's [26] formalism have not be performed yet.
In this paper we formulate the direct t-matrix approach to determination of the electric polarizabilities of a bound system that relies on the solution of the integral equations for both the bound-state wave function and corresponding components of the partial transition matrix of the system.In Section 2, following the Watson-Feshbach method, we express the polarization potential describing interaction between a charge particle and a bound complex consisting of particles in terms of a "truncated" Green's operator of the system.In this way we obtain a general expression for the electric multipole polarizability of the system.In Section 3 the t-matrix approach to determination of the electric polarizabilities is formulated and simplifications of the general formula assuming conservation of the space parity and the total orbital moment of the system are considered.Section 4 is devoted to the application of the elaborated formalism to the two-body bound systems with the central interaction.It is shown that the electric N 2  -pole polarizability of the two-particle -state bound complex contains information both on derivatives (of the order S  and lower) of the wave function in momentum space and on the partial component of the transition matrix that corresponds to the orbital state with l   .Section 5 contains the ap- plication of the t-matrix formalism to determination of the electric multipole polarizabilities of the hydrogen atom and discussion.Section 6 is devoted to conclusion.

Polarization Potential
The formula for the polarization potential describing the interaction between a charged particle 0 and a bound complex consisting of particles follows immediately when treating the body problem, namely, the low-energy scattering of the complex by the Coulomb field of the charged particle 0 with the kinetic energy of the relative motion of the particle 0 and the complex being well below the breakup threshold energy of the complex.Experimentally, to determine the polarizability of the complex in the direct way a heavy nucleus (the particle 0) is used as a source of the intense electric field.For simplicity sake assume that the -body complex contains only one charged particle 1, the other particles The Hamiltonian of whole system of particle has the form (3) where 0 is the operator of the kinetic energy of the relative motion of the particle 0 and the center of mass of the complex, H is the total Hamiltonian of the body complex, is the kinetic energy operator of the relative motion of the particles inside the complex, is the total interaction potential of the  , is the operator of the Cou-complex, lomb interaction between the charge particles 0 and 1 (with the charges 0 and 1 e ), ij is the operator of the interaction between the particles of the complex and .Using the known Watson-Feshbach projection technique [2,27,28] with the projection operators 0 0 and  is the wave function of the ground bound state of the complex of particles N   1, 2, , N  B with the binding energy 0 (normalized to one), the operator of the effective interaction may be written as where the averaging is taken over the variables of the relative motion of the particles of the complex.The operator

 
  satisfies the integral equation in which the potential energy operator is the "external" Coulomb interaction potential between the field source (the particle 0) and the charged particle 1 of the complex, 01 , and the propagator is a "truncated" Green's operator containing the "internal" interaction potential of the complex , H where is the free Hamiltonian of the complete system (of  particles), 0 , is the total energy of the system, Introducing the transition operator which satisfies the integral Lippmann-Schwinger equation   int   can be written in the form the operator where the free propagator is given by     .
The polarization potential that corresponds to the second order of the perturbation expansion (in powers of the Coulomb interaction 01 ) of the operator in the expression for the effective potential (4) is given by where According to the uncertainty principle, in the case of asymptotically large (in comparison with the complex size) distances between the particle 0 and the complex, 0  , the momentum variables of relative motion among particles inside the complex considerably exceed the momentum variable of relative motion of the particle 0 and the center of mass of the complex.In such a case, in the expression for the "truncated" Green's function int in Equation ( 9) we can neglect with the variable quantity that corresponds to the operator Then the polarization potential takes the known local form where each of coefficients in the sum (11), E  , that characterizes the strength of the individual term of the polarization potential with the asymptotic 0 is the electric polarizability of the -particle complex of the multipolarity Here M  is the multipole moment operator of the charged constituent particle 1, 1 is the radius vector of the particle 1 relative to the center of mass of the bound complex; the radius vector 0 r ρ describes the position of the center of mass of the bound complex relative to the charged particle 0 creating the electric field,  is Legendre polinomial (the unit vectors are marked with a hat, ); the quantity 0 is the "truncated" Green's operator of the complex at the bound state Green's function of the complex.Note that Formula (1) readily follows from the general Formula ( 12) after applying the spectral expansion of the total Green's operator in the complete set of the eigenfunctions of the total -particle Hamiltonian H .

T-Matrix Method of Determination of the Electric Multipole Polarizabilities of a Bound Complex
In this Section we formulate a method for the direct calculation of the polarizability of a few-particle bound system starting immediately from the definition ( 12) without recourse to expansion in terms of a set of excited discrete and continuum wave functions.
The general Formula ( 12) for the electric multipole polarizability of a bound complex of particles N   at the energy of the bound state 0 .The wave function  satisfies the Schrödinger equation or the equivalent homogeneous integral equation We shall restrict our consideration to the case that the complex is in the ground bound state.
The "truncated" propagator , as distinct from the total one , does not contain the pole singularity at   .In an explicit form, we find the expression for in terms of the free propagator and the transition operator ,


where the operator In view of Equation ( 17) and the equality the expression ( 14) takes the form Separating out from the total transition operator the singular pole part that corresponds to the ground bound state of the system at the energy   , we write where the vertex function 0 is expressed through the wave function of the ground bound state of the system, satisfying the homogeneous integral equation and denotes the smooth (non-singular) part of the transition operator.
Taking into account Expression (21) we write Formula (20) for the operator in the form Cancelling the pole terms in Equation ( 24) with the use of the identity we may write Expression (24) in the form At the point of the negative energy of the bound state of the -particle complex, 0 , the operator 0 in the general formula for the electric dipole polarizability (12) becomes With the use of Expression (27) the formula for the polarizability (12) takes the form

         
A noticeable simplification of the general formula for the electric multipole polarizability (28) takes place, if the interaction potential V is invariant relative to the space reflection (for example, for the systems with the Coulomb or nuclear interactions).In this case the wave function of the bound complex is characterized by a definite parity.The conservation of the parity leads to nullification of the matrix elements   0 0 0 0 0 0 0 0 , and , which are present in the second and third summands of Formula ( 28), at odd values of  .
Formula ( 28) is also simplified in the case of the invariance of the interaction relative to rotations that leads to conservation of the angular momentum.In such a situation the matrix elements (29) are proportional to the Clebsh-Gordan coefficient  satisfying the triangle condition . Specifically, they vanish after integrating in angular variables for all   , if the total orbital moment of the complex is equal to zero, 0 L  .
In the case when the ground bound state of the complex is characterized with the total orbital momentum 0 L  , what is realized for the simplest atomic and few-hadron systems, the formula for the electric multipole polarizability (28) is simplified to the expression As is evident from (28) or (30), to determine the electric dipole polarizability of the bound complex, it is necessary to know not only its wave function, but the smooth part of the transition matrix at negative energy of the bound state   E B   N 0 as well.The corresponding transition matrix can be determined by solving the Lippmann-Schwinger integral equation for the two-particle system, the Faddeev integral equations [26] for the threeparticle system or the Faddeev-Yakubovsky equations [29] for more complex systems.

Electric Multipole Polarizabilities of the Two-Particle Bound Complex
In the case of two-particle bound complexes, the derived formula for the electric multipole polarizability (30) is simplified.Considered here is a stable bound complex consisting of a charged particle 1 and a neutral particle 2. The interaction potential between the particles is taken to be central, the complex is in S-wave ground bound state.We denote quantities describing the two-particle system by small letters as distinct to notations by capital letters in the general case of quantities for the -particle complex used in the foregoing Sections:  are the binding energy and the wave function of the two-particle bound complex, is the free Green's operator,   is the reduced mass of the particles 1 and 2. Further consideration we perform in the momentum space.
Starting from (30) and taking into consideration that the smooth part of the two-particle transition matrix 0 has form of the sum of the smooth part of its -wave partial component (with ), 0 0 , and the sum of all higher partial orbital components (with ), 0 , which are non-singular at the energy of the bound state, we write the formula for the electric multipole polarizability of the two-particle system in the form The wave function of the complex in the S-wave bound state satisfies the homogeneous integral equation and each of the partial components of two transition matrix l satisfies the inhomogeneous integral Lippmann-Schwinger equations where l is the partial of the interaction potential, is the energy of the relative motion of the particles.
The solution of the equation for the partial transition matrix (34) that corresponds to the orbital moment of the ground state has the form of the sum of the pole and smooth operators where the vertex function satisfies the homogeneous integral equations that follows from the Equation (33), Using the operator 0 0 , which is determined by the inhomogeneous equation we write the smooth part of the transition operator, At the point 0    , the inhomogeneous Equation (37) becomes homogeneous one, Since the kernels of the equations for the function 0  (36) and the operator 0 0 (37) coincide, the solution of the operator Equation (39) may be written as According to (38), the operator 0 0 , which is contained in the formula for polarizability (32), is related to the operator , we obtain the expression for the operator , where In case of the central interaction, the part of the transition matrix 0  in (32) may be written in the form of expansion in the set of the spherical functions of the angular variable momenta (with ), where the notations and 12 12 are used, we perform the integration with respect to angular variables.
Notice that the contribution of the smooth part of the -wave components of the transition matrix,   t b   l 0 0 , which is of the separable form (43), proves to be equal to zero as a result of the conservation of the space parity and the orbital moment of the relative motion (similar to zero contributions in the general case from other factorable terms of the "truncated" Green's operator (27), see comments below Formula ( 28)).Non vanishing contribution in (32) where  is the partial component of the interaction potential.
The final result for the electric multipole polarizability of the two-particle complex with a central interaction between particles following from Equation ( 32) is written in the form According to (47), the functions .
Formula (49) derived in the case of the central interaction between the constituents of the two-body complex demonstrates that the electric  -pole polarizability of the complex contains information not only on the derivatives (of the order  and lower) of its wave function, but on the partial component of the transition matrix,   It is advantageous to write Formula (49) in the form more convenient for the practical use.In order to do this, instead of the two-particle t-matrix, which is the function of two momentum variables and , we introduce in (49) the function of one variable , as well.
that satisfies the inhomogeneous integral equation with the kernel of the Lippmann-Schwinger Equation (48) in which the free term is determined by the formula In such a case, Formula (49) takes the form

Approach and Discussion
Advantages of the t-matrix approach are manifested in calculation of polarizabilities of a quantum system, even if it is two-particle (the deuteron nucleus or the hydrogen atom).In the special case that the interaction between particles is purely S-wave (for example, the separable interaction potential [30], in Formula (49) only the first term persists.The corresponding reduced formula for the dipole polarizability of the two-particle bound complex has been earlier derived in the framework of the threebody formalism of the effective interaction of a charged particle and a complex [13,14,31].Applications of the t-matrix approach to calculation of the electric dipole polarizabilities of the deuteron, the triton and the lambda hypertriton as two cluster systems with purely -wave interaction have been carried out in the preceding our papers [23,24].
For the deuteron with noncentral interaction between the proton and the neutron (in the state with the total angular momentum 1) the polarization in the electric field is anisotropic.Using the separable tensor potential [32] the longitudinal and transverse (relatively to the direction of the electric field) components of the dipole polarizability of the deuteron have been calculated in the work [15].The results of the further calculations of the components of the deuteron electric dipole polarizability obtained in the framework of the chiral effective field theory [18] are in agreement with the results of [15].
In the case of the hydrogen atom H (assuming that the proton mass is infinitely great when compared to the electron mass), the exact values of the electric multipole polarizabilities E are known [25].This is why it is possible to test directly the validity of the general formula for the polarizability of the two-particle system (49)both of the term with free propagation in intermediate state and of the terms with the multiple scattering in the higher orbital intermediate states (with for the dipole, quadrupole, octupole and higher multipole polarizabilities, Inserting the analytical expression for the partial components of the Coulomb transition matrix, obtained from the representation for the three-dimensional t-matrix derived in [33] with the use of the O(4) rotation symmetry in four-dimensional Fock space [34], into Formula (49) and separating out the Born term from the Coulomb -matrix, we write as a sum of three terms, where  We show in the The elaborated approach to determination of the electric multipole polarizabilities that relies on the two-particle t-matrix differs radically from the traditional approach based on the spectral expansion of the Green's function of the system that considers contributions from an infinite number of excited bound and continuum states.As evidenced by the results of the spectral-expansion calculations by Castillejo et al. [35], 65.8% of the magnitude of the electric dipole polarizability of the hydrogen atom comes from the -wave excited bound states.Inclusion of all the excited bound states accounts for 81.4%.The rest 18.6% is provided taking into consideration the continuum states.P N

Conclusions
In conclusion, a new approach to the determination of the electric multipole polarizabilities of quantum bound systems (both atomic and nuclear) that is based on the use of the transition matrix has been developed.The most important advantage of the t-matrix formalism is its application allows to avoid cumbersome calculations of individual contributions from infinitely large number of discrete and continuum excited states.Instead, the determination of one or several (depending on the interaction mode) partial components of the transition matrix is now required.It is also essential that the transition matrix in the main Formula (30) depends on the negative energy of the bound state being therefore a real function of momenta.
To test the validity of the proposed method, the calculation of the electric dipole, quadrupole and octupole polarizabilities of the hydrogen atom has been performed.The obtained results precisely reproduce the known analytically derived values of the polarizabilities.In nuclear physics, in the special case of purely S-wave separable interaction potential, the general Formula (49) is simplified and reduced to one obtained earlier for the deuteron.
The developed method can be immediately extended to more complicated interactions between constituents, specifically, to the case of tensor that gives rise to anisotropic polarization properties of the system.Finally, it is quite important that the proposed approach is suitable for more complicated three-and -body systems described by the Faddeev and Faddeev-Yakubovsky integral equations, since the electric polarizabilities of the few-body nuclei are important characteristics containing additional independent information about the fundamental nuclear force.In the first place, we plan to apply the t-matrix approach to study the deformation properties of the 3 H and 3 He nuclei in the electric field.

0
can be deduced by performing the differentiation of the Equation (37) with respect to and the inverted transition from to t and  , the point  , which appears in view of the relation (40) and the normalization condition

)
Substituting the expressions (43) and (45) for   for the electric polarizability(32) and taking account of the action of the multi- makes only the partial component with the Lippmann-Schwinger integral Equation (34),

 1 
49) are expressed in terms of the derivatives of the wave func-tion of the ground bound state with respect to the relative momentum variable .In the cases of the dipole 

Table 1
the values of the components