studies of electron - positron pair production in strong external fields are an important tool to reveal the structure of the qed vacuum .
dynamical pair creation by two coulomb fields in relativistic heavy ion collisions has been analysed in great detail , both theoretically and experimentally ( see , e.g. , @xcite for a topical review ) . in recent years , due to the sustained progressing of laser technology , pair creation in strong laser fields is encountering a growing interest . in several theoretical papers pair production in the standing wave formed by two counterpropagating laser beams
has been considered @xcite . in most cases a laser field strength of order or above the critical value @xmath0 v / cm is required in order to observe the process in experiment @xcite , which still is by four orders of magnitude larger than the highest laser field strengths achievable today @xcite .
more promising in view of its experimental realization is pair production in the combination of a laser and a coulomb field . here
one can exploit the possibility to let a nucleus or a highly charged ion collide at large value of the relativistic lorentz factor @xmath1 with the laser beam .
then , in the nuclear rest frame , the laser s field strength @xmath2 and frequency @xmath3 are enhanced by a factor @xmath4 .
observable pair creation rates should result for @xmath5 or @xmath6 in this frame ( here , @xmath7 , @xmath8 and @xmath9 are planck s constant , light velocity and electron mass , respectively ) . with the help of the upcoming novel laser @xcite and accelerator @xcite facilities
these conditions can be met .
the large hadron collider , that is presently under construction at cern ( geneva , switzerland ) , will accelerate protons to an energy of 7 tev ( @xmath10 ) @xcite .
when such protons are brought into collision with the strongest laser beams available today ( @xmath11 at @xmath12 ev ) , then the laser field strength in the proton s rest frame approaches or even exceeds the critical value . on the other hand , the x - ray free - electron laser ( xfel ) facilities
presently being developed at desy ( hamburg , germany ) and slac ( stanford , usa ) are proposed to provide spatially coherent and highly brilliant beams of synchrotron radiation with single - photon energies of up to @xmath13 kev at maximum field strengths of @xmath14 @xcite . when combined with a moderately relativistic ion beam
, the laser photon energy in the projectile frame can reach the electron s rest energy .
we should mention that a few years ago nonlinear electron - positron pair creation in the collision of an ultrarelativistic electron beam ( @xmath15 ) and an intense optical laser pulse was observed at slac @xcite .
in this situation , there are two channels for pair production : a direct channel , via laser - photon absorption in the coulomb field of the electron , and an indirect two - step channel ( so - called breit - wheeler process ) , where first a high - energy photon is generated by compton backscattering which afterwards creates the pair via absorption of laser photons @xcite . in the experiment
the contribution of the indirect channel was shown to be dominant . in the present paper , however , the indirect production mechanism will not be addressed .
it is of importance only in the case of a light projectile beam .
in fact , when the incoming electrons were replaced by a beam of heavy particles ( e.g. , protons ) the indirect production channel would be strongly suppressed and direct nonlinear pair production ( by a virtual photon from the nuclear coulomb field ) could be investigated .
the latter process is the subject of this paper .
the first calculation of pair creation by a coulomb and a strong laser field is due to yakovlev @xcite .
he treated the case of circular laser polarization and derived asymptotic formulae for the production rate in certain parameter regimes .
particularly , the total cross sections in the perturbative regime ( @xmath16 , with the intensity parameter @xmath17 given in eq.([xi ] ) below ) as well as in the overcritical - field limit ( @xmath18 and @xmath19 ) were obtained .
twenty years later , mittleman @xcite considered the process in a linearly polarized laser wave .
his calculation of the leading term for the total cross section in the perturbative multiphoton regime ( @xmath16 and @xmath20 ) led him to the conclusion that , for the laser intensities and frequencies available at that time , this cross section was the smallest one on record . despite of this pessimistic result , another twenty years later the great headway in laser technology that was following the invention of chirped - pulse amplification @xcite triggered newly revived interest in the process .
dietz and prbsting @xcite investigated the influence of a highly charged ion of charge @xmath21 within a nonperturbative numerical approach .
however , to make their coupled - channel calculations feasible they treated the laser field in dipole approximation . while in the earlier work the nucleus was always assumed at rest , mller , voitkiv , and grn explored in detail the improved opportunities arising from laser - nucleus collisions . in the articles @xcite they performed numerical calculations for tunneling pair production by an ultrarelativistic nucleus colliding with a superintense near - infrared laser beam and for multiphoton pair creation by a moderately relativistic ion beam in combination with an intense x - ray laser wave .
moreover , they considered the variant of bound - free pair production where the electron is created in a bound state of the projectile nucleus @xcite .
free pair creation in the collision of a relativistic nucleus with an ultrastrong x - ray laser of circular polarization was also treated by avetissian , avetissian , mkrtchian , and sedrakian in ref.@xcite , where analytical formulae in the high - intensity limit and some numerical results on the energy spectrum of the created particles in the nonperturbative multiphoton regime ( @xmath22 and @xmath23 ) are given .
most of the treatments mentioned above fully account for the influence of the external laser field by using the volkov solutions to the dirac equation as basis states in a perturbative calculation of first order with respect to the nuclear coulomb field . within this framework , the fully differential pair creation rate can be expressed as a fourier series over the number of absorbed laser photons , the coefficients of which are essentially given by ordinary or generalized bessel functions for circular or linear laser polarization , respectively . because of the large photon numbers involved , taking this sum is one of the main obstacles in this kind of approach , in particular for a linearly polarized laser wave . as a consequence , it is difficult to get useable expressions for the total production rate .     pictorial equation in terms of feynman diagrams describing @xmath24 pair production in combined laser and coulomb fields . in our approach
, the coulomb field is treated within the first order of perturbation theory , while the effect of the laser wave is nonperturbatively taken into account to all orders .
this is expressed by the feynman graph on the left - hand side , where the double lines represent the exact lepton wave - functions in the laser field ( volkov states ) , while the dashed line stands for their interaction with the coulomb field .
expanding the volkov states with respect to the lepton - laser coupling results in a perturbation series , some typical low - order terms of which are shown on the right - hand side .
the wavy lines symbolize the laser photons and the arrows indicate whether the respective laser photon is emitted or absorbed during the process . the first diagram on the right - hand side ,
e.g. , describes the leading order of pair production by the net - absorption of one laser photon ( with @xmath25 ) . for later reference
, we note that the interference between the first and the fourth diagram leads to a @xmath26-correction to the probability for one - photon pair creation [ cf .
eq.([delta ] ) below].,width=302 ]    in the present paper a different approach based on the optical theorem is used to calculate the total rate for pair creation in combined laser and coulomb fields .
this approach employs the explicit form of the polarization operator of a photon in a laser field found by baier , milstein , and strakhovenko @xcite by means of an operator technique .
an alternative form of this polarization operator was derived independently by becker and mitter @xcite by means of another method .
both results are in agreement with each other .
similar to refs.@xcite the interaction of the leptons with the laser field is taken into account to all orders while the effect of the coulomb field is treated in first order ( cf .
1 ) . the main advantage of our method is the possibility to derive , independent of the laser s polarization state , in an analytical manner compact formulae for the total production rate that only involve low - dimensional integrals and no additional summations .
based on these representations one can rather easily find all different kinds of asymptotics . in this way we will confirm and extend previous results for a circularly polarized laser wave and give for the first time corresponding expressions for elliptic laser polarization . as a main result , we show that there are essential differences between the general case of an elliptically polarized laser field and the special case of circular polarization . furthermore , because of its interest for near - future experiments , particular emphasis is placed on the nonlinear process of pair creation by the simultaneous absorption of two photons from the laser wave in a situation , where the energy of a single photon is not sufficient to produce a pair and where the contribution from higher photon orders is neglibly small .
we note that , in general , contributions stemming from different net - numbers of absorbed laser photons can be distinguished via the momentum spectra of the produced particles .
in addition , the strong - field regimes of pair production shall be examined .
the paper is organized as follows . in the next section
we describe how the polarization operator is related to the process of pair creation , and how it can be used to calculate the total pair production rate in the superposition of a coulomb and a laser field .
general formulae are given in terms of two - fold integrals for both circular and elliptic laser polarization .
in the following section we apply these expressions to various parameter regions and determine the corresponding asymptotic behaviour of the pair production rate .
special attention is paid to the nonlinear process of pair production by two - photon absorption in the limit of low laser intensity , for which in particular the behaviour close to the energetic threshold is investigated . moreover , the pair - creation cross - section in the high - frequency limit ( @xmath27 ) is analysed , where we calculate the next - to - leading order correction to one - photon pair production at low intensity and analyse the limiting case of high laser intensity .
finally , we consider the overcritical - field regime as well as the quasiclassical limit of high intensity and low frequency , where the rate shows an exponential tunneling behaviour . we finish with a conclusion where the significance of our results in view of possible experimental investigations is discussed .
the polarization operator @xmath28 describes the propagation of a photon of four - momentum @xmath29 in a background field ( e.g. , the qed vacuum ) , including self - energy corrections .
it is related to the corresponding exact photon propagator @xmath30 via dyson s equation @xcite @xmath31 where @xmath32 denotes the free photon propagator . in refs.@xcite
the ( properly renormalized ) polarization operator was evaluated in the combined fields of the qed vacuum and a classical electromagnetic plane wave . the interaction with the vacuum
was taken into account to first order in the fine - structure constant @xmath33 , with @xmath34 being the electron charge , while the interaction with the plane - wave field was included to all orders . in the following
, the shape of the electromagnetic wave will be taken as @xmath35 with @xmath36 and the dimensionless intensity parameters @xmath37 here and henceforth we use relativistic units with @xmath38 and write @xmath39 for the product of two four - vectors .
the respective result on the polarization operator of a photon in a constant electric field instead of a plane electromagnetic wave can be found in ref.@xcite .
apart from its theoretical significance , the knowledge of the polarization operator allows for a number of applications .
for example , the total probability for @xmath24 pair production by a photon of momentum @xmath40 in an external field ( e.g. , a plane laser wave ) is related to the imaginary part of the corresponding polarization operator via @xcite @xmath41 where @xmath42 denotes the photon s polarization four - vector . in fig.2
we give an intuitive , pictorial explanation of this formula .
lowest order ( in @xmath43 ) feynman graph for the polarization operator of a photon in an external laser field . in our case
the photon is a virtual one stemming from a coulomb field ( outer dashed lines ) .
the diagram can be viewed as describing the elastic forward scattering of the photon through an intermediate laser - dressed @xmath24 state ( double line ) . by performing ,
as indicated , a cut along the central dashed line , it is seen that the graph represents the product of the amplitude for the production of a laser - dressed @xmath24 pair by the virtual photon and its complex conjugate , including a sum over all electron states .
application of the optical theorem ( i.e. , the unitarity of the scattering matrix ) therefore relates the imaginary part of the polarization operator with the total probability for pair creation , as expressed by eqs.([ppph ] ) and ( [ w]).,width=94 ]      the calculation in ref.@xcite was performed for an arbitrary photon momentum @xmath29 , including the case @xmath44 .
therefore , the polarization operator found there not only applies to the combination `` photon + laser wave '' .
instead , the photon may be replaced by any additional external field @xmath45 , e.g. a nuclear coulomb field .
since a coulomb field can transfer momentum but not energy , here one has @xmath46 . to get , at a given value of @xmath47 , the differential probability ( per unit time ) for @xmath24 pair creation one has to replace the photon wave - function @xmath48 in eq.([ppph ] ) by the fourier transform of the coulomb field @xmath49 . here
, @xmath50 denotes the nuclear charge number .
the total production rate then is found by integration over all possible momenta : @xmath51 note that in eq.([w ] ) the additional factor of 1/2 compensates for the double - counting of the contributions from the momenta @xmath47 and @xmath52 .
for the case of a constant homogeneous electromagnetic field the relation ( [ w ] ) was used in ref.@xcite , where the probability for pair production was found within the quasiclassical approximation for the polarization operator .    according to the results of ref.@xcite , the required component of the polarization operator in eq.([w ] )
can be written as @xmath53 \right\ } \nonumber\\ & & \!\!\!\!\ !
\times{1\over\cos^2\theta } \left [ \sin^2\theta\ , ( d_3 \cos^2\phi + d_4 \sin^2\phi ) - d_5 \right],\end{aligned}\ ] ] where we introduced the angles @xmath54 and @xmath55 , and write @xmath56 . the coefficients @xmath57 ( @xmath58 ) in eq.([pi00 ] ) are given by eq.(2.32 ) in ref.@xcite ; they read @xmath59 \nonumber\\ & & + \xi_1 ^ 2\sin^2\rho{1+v^2\over 1-v^2}j_0(\zeta ) \nonumber\\ & & + { 1\over 4}\left[{q^2\over m^2}-i{|\lambda|(1-v^2)\over\rho}\right][j_0(\zeta)-{\rm e}^{iy } ] \nonumber\\ d_4 & = & d_1(\xi_1 \leftrightarrow \xi_2)\nonumber\\ d_5 & = & { q^2\over 4m^2}(1-v^2)[j_0(\zeta)-{\rm e}^{iy}]\end{aligned}\ ] ] with the abbreviations @xmath60    for simplicity , we first consider the case of a circularly polarized laser wave ( @xmath61 ) .
then the coefficients in eq.([pi00 ] ) can be simplified to read @xcite
@xmath62(1-{\rm e}^{iy})\nonumber\\ d_5 & = & { q^2\over 4m^2}(1-v^2)(1-{\rm e}^{iy})\end{aligned}\ ] ] with @xmath63 employing spherical coordinates such that @xmath64 , the integration over @xmath65 in eq.([w ] ) is trivial and the total rate becomes @xmath66\right ] \nonumber\\ & & \times \bigg\ { ( 1-t^2)\xi^2\sin^2\rho { 1+v^2\over 1-v^2 } { \rm e}^{-iy }   + \bigg [ q^2(1-v^2 ) \nonumber\\ & &   + { 1-t^2\over 2}\left [ 1-q^2(1+v^2 ) \right ] \bigg](1-{\rm e}^{-iy } ) \bigg\},\end{aligned}\ ] ] where we have introduced the dimensionless quantities @xmath67 , @xmath68 , and @xmath69 .
going over to the new variables @xmath70 , @xmath71 considerably simplifies the @xmath72-dependence such that also this integral can be taken : @xmath73   \bigg\ { \xi^2\sin^2\rho \nonumber\\ & & \!\!\!\ !
\times\left[{2\over 3}(2+t^2){\sqrt{1-t^2}\over t^2}- \ln\left({1+\sqrt{1-t^2}\over 1-\sqrt{1-t^2}}\right)\right ] { \rm e}^{-iy}\nonumber\\ & & \!\!\!\ !
3}{(1-t^2)^{3/2}\over t^2}+{q^2\over 2 } \ln\left({1+\sqrt{1-t^2}\over 1-\sqrt{1-t^2}}\right)\right](1-{\rm e}^{-iy } ) \bigg\ } \nonumber\\ \\end{aligned}\ ] ] where now @xmath74 .
alternatively , the @xmath75-integration can be expressed by neumann functions @xmath76 : @xmath77 \nonumber\\ & & + { 2\over 3}\left [ y_1(\beta)-{1\over z}y_1(\beta z)\right]{(1-t^2)^{3/2}\over t^2 } \nonumber\\ & & + \big [ y_1(\beta)- z y_1(\beta z)\big ] \ln\left({1+\sqrt{1-t^2}\over 1-\sqrt{1-t^2}}\right ) \bigg\}\end{aligned}\ ] ] with @xmath78 and @xmath79^{1/2}$ ] . equation ( [ circ2 ] ) is a very compact , general expression for the total pair production rate in a circularly polarized laser field . for applications , both eqs.([circ1 ] ) and ( [ circ2 ] ) will prove to be useful .    in the general case of an elliptically polarized laser wave the integrations over @xmath65 and @xmath72 can be taken in the same manner as before yielding @xmath80 { \rm e}^{-iy } \nonumber\\ & & \!\!\!\ !
\times\bigg\ { \left(\xi_1 ^ 2+\xi_2 ^ 2\right)j_0(\zeta )              \bigg [ { 2\over 3}a_1(1-t^2)^{3/2 } \nonumber\\ & & \!\!\!\ !
+ \sin^2\rho \bigg({2\over 3}(2+t^2)\sqrt{1-t^2 }               - t^2 \ln\left({1+\sqrt{1-t^2}\over 1-\sqrt{1-t^2}}\right)\bigg ) \bigg ] \nonumber\\ & & \!\!\!\ !
+ i\left(\xi_1 ^ 2-\xi_2 ^ 2\right)j_1(\zeta ) \bigg [ { 2\over 3}a_1(1-t^2)^{3/2 } \nonumber\\   & & \!\!\!\ !
-\sin^2\rho \left(2\sqrt{1-t^2 }               - t^2 \ln\left({1+\sqrt{1-t^2}\over 1-\sqrt{1-t^2}}\right)\right ) \bigg]\nonumber\\ & & \!\!\!\ !
+ \left[j_0(\zeta)-{\rm e}^{iy } \right ]   \bigg[{2\over 3}(1-t^2)^{3/2}\left(q^2-i{qt\over a\rho}\right ) \nonumber\\ & & \!\!\!\ !
-q^2t^2\ln\left({1+\sqrt{1-t^2}\over 1-\sqrt{1-t^2}}\right)\ \bigg ] \bigg\}.\end{aligned}\ ] ] the integral over @xmath75 can again be taken analytically : @xmath81 \nonumber\\ & & \!\!\!\ !
\times \bigg [ { 2\over 3}a_1(1-t^2)^{3/2 } \nonumber\\ & & \!\!\!\ !
+ \sin^2\rho \bigg({2\over 3}(2+t^2)\sqrt{1-t^2 }               - t^2 \ln\left({1+\sqrt{1-t^2}\over 1-\sqrt{1-t^2}}\right)\bigg ) \bigg ] \nonumber\\ & & \!\!\!\ ! - { \xi_1 ^ 2-\xi_2 ^ 2\over \sqrt{g_0 ^ 2-g_1 ^ 2 } }              \left [ \varphi_- j_1'(\varphi_-)y_1(\varphi_+ )                    - \varphi_+ j_1(\varphi_-)y_1'(\varphi_+ ) \right ] \nonumber\\ & & \!\!\!\ !
\times { \rm sign}(g_1 ) \bigg [ { 2\over 3}a_1(1-t^2)^{3/2 } \nonumber\\   & & \!\!\!\ ! -\sin^2\rho \left(2\sqrt{1-t^2}-t^2
\ln\left({1+\sqrt{1-t^2}\over 1-\sqrt{1-t^2}}\right)\right ) \bigg]\nonumber\\ & & \!\!\!\
! - \left [ \varphi_- j_1(\varphi_-)y_0(\varphi_+ )                      + \varphi_+ j_0(\varphi_-)y_1(\varphi_+ ) - \beta y_1(\beta)\right ] \nonumber\\ & & \!\!\!\ !
\times\bigg[{2\over 3}(1-t^2)^{3/2 } - t^2 \ln\left({1+\sqrt{1-t^2}\over 1-\sqrt{1-t^2}}\right)\bigg ] \nonumber\\ & & \!\!\!\ !
-{4\over 3 } \left [ j_0(\varphi_-)y_0(\varphi_+ ) - y_0(\beta ) \right ] ( 1-t^2)^{3/2 } \bigg\}\end{aligned}\ ] ] with @xmath82 , @xmath83 , and @xmath84 .
the total production rate can be converted into a cross section according to the relation @xmath85 with @xmath86 being the photon flux .
let us consider the case of small intensity ( @xmath16 ) and frequency @xmath87 .
the leading term in the @xmath26-expansion corresponds to pair creation by a single photon in a coulomb field . under these circumstances
the total production probability is independent of the laser s polarization state and it is convenient to derive the asymptotic behaviour by using the formula ( [ circ1 ] ) for circular laser polarization . assuming @xmath16 , we can expand the exponentials to lowest order with respect to @xmath88 .
afterwards we perform the integration over @xmath89 with the result @xmath90 \nonumber\\ & & \times \left(1-{\gamma\over 2}\right ) \bigg\},\end{aligned}\ ] ] where @xmath91 and @xmath92 is the step function , which restricts the integration to the kinematically allowed region @xmath93 and @xmath94 with @xmath95 . integrating over @xmath75
gives : @xmath96 \nonumber\\ & & - 2\left ( 1+{a^2\over t^2}\right)\sqrt{{t^2\over a^2}-1 } \ln\left({1+\sqrt{1-t^2}\over 1-\sqrt{1-t^2}}\right ) \nonumber\\ & & + { 2a\over t}\ln\left({1+\sqrt{1-t^2}\over 1-\sqrt{1-t^2}}\right ) \ln\left({{t\over a}+\sqrt{{t^2\over a^2}-1}}\right ) \bigg\}. \nonumber\\ \ \end{aligned}\ ] ] the corresponding expression for the cross section [ see eq.([flux ] ) ] agrees with the well - known result of bethe and heitler @xcite .
accordingly , in the ultrarelativistic limit ( @xmath97 ) we find @xmath98\end{aligned}\ ] ] with the classical electron radius @xmath99 , and in the nonrelativistic limit ( @xmath100 ) @xmath101      it is interesting to consider the cross section for two - photon pair production at @xmath16 , which represents a nonlinear strong - field process of lowest possible order in the number of absorbed photons .
we assume that @xmath102 , which excludes the possibilty of one - photon pair creation .
going one step further in the @xmath88-expansion of eq.([circ1 ] ) and performing the @xmath89-integration , we obtain for the circular polarization case @xmath103 + { 1\over 2}\left(1-{\gamma\over 4}\right)^2{a^2\over q^2t^2 } \nonumber\\ & & \!\!\!\!\!\!\!\ !
\times \left[{2\over 3}{(1-t^2)^{3/2}\over t^2}+q^2 \ln\!\left({1+\sqrt{1-t^2}\over 1-\sqrt{1-t^2}}\right)\right]\bigg\}.\end{aligned}\ ] ] proceeding as before , we get @xmath104 \nonumber\\ & & + { ( 4t^2-a^2)^{5/2}\over 45at^4}{(1-t^2)^{3/2}\over t^2 } + { a\over 6t^4 } \nonumber\\   & & \times\ln\left({1+\sqrt{1-t^2}\over 1-\sqrt{1-t^2}}\right ) \bigg [ ( 2t^2+a^2)\sqrt{4t^2-a^2 } \nonumber\\ & &   -3a^2t\,\ln\!\left({2t\over a}+\sqrt{{4t^2\over a^2}-1}\right)\bigg]\bigg\}.\end{aligned}\ ] ] after conversion into a cross section , eq.([2ph ] ) reads @xmath105 with @xmath106^{3/2}\over ( xt)^2 } \nonumber\\ & & \times \left[{2\over 3}(2+t^2){\sqrt{1-t^2}\over t^2}- \ln\left({1+\sqrt{1-t^2}\over 1-\sqrt{1-t^2}}\right)\right ]   \nonumber\\ & & + { 4\over 15}{[(xt)^2 - 1]^{5/2}\over ( xt)^4}{(1-t^2)^{3/2}\over t^2 } + { 1\over ( xt)^4 } \nonumber\\ & & \times\ln\left({1+\sqrt{1-t^2}\over 1-\sqrt{1-t^2}}\right ) \big [ \left(2+(xt)^2\right)\sqrt{(xt)^2 - 1 }   \nonumber\\ & & -3xt\ln\left ( xt+\sqrt{(xt)^2 - 1}\right ) \big ] \bigg\}.\end{aligned}\ ] ] in the limit @xmath107 , eq.([f ] ) becomes @xmath108 hence , the threshold behaviour of the cross section for two - photon pair production in a circularly polarized laser field reads @xmath109 we note that , based on the results of their numerical calculations , the authors of ref.@xcite found that the rate for two - photon pair creation within the whole frequency range @xmath102 approximately scales as @xmath110
. however , eq.([2phlimit ] ) now shows that the exact scaling in the threshold region ( @xmath111 and @xmath112 ) is given by @xmath113 .
further , in ref.@xcite a total production rate of @xmath114 sec@xmath115 ( in the nuclear rest frame ) was reported for @xmath116 , @xmath117 kev , and @xmath118 . applying eq.([2ph ] ) to these parameters , perfectly reproduces this result .    in the general case of elliptic polarization
, the corresponding asymptotic cross section reads @xmath119 where @xmath120 with the function @xmath121 from eq.([f ] ) and @xmath122 \nonumber\\ & & -\left [ ( xt)^2 - 4 \right]{\sqrt{(xt)^2 - 1}\over xt } \nonumber\\ & & \times \left ( 2\sqrt{1-t^2}-t^2\ln\left({1+\sqrt{1-t^2}\over 1-\sqrt{1-t^2}}\right ) \right ) \nonumber\\ & & + { 1\over 2(xt)^3 }         \big [ \left ( ( xt)^2 - 16\right)\sqrt{(xt)^2 - 1 } \nonumber\\ & & + 12xt\ln\left(xt+\sqrt{(xt)^2 - 1}\right)\big ] \nonumber\\   & & \times \bigg[{2\over 3}(1-t^2)^{3/2 } - t^2\ln\!\left({1+\sqrt{1-t^2}\over 1-\sqrt{1-t^2}}\right)\bigg]\end{aligned}\ ] ] and the ellipticity parameter @xmath123 ^ 2 $ ] .
this time the threshold behaviour ( @xmath124 ) is found to be @xmath125 hence , in the general case , the cross section for two - photon pair creation close to threshold raises like the second power of the excess energy @xmath126 . only in the special case of circular polarization ( @xmath127 )
the threshold cross section scales with the fourth power as given in eq.([2phlimit ] ) .
we notice , however , that our first - order treatment of the coulomb field requires the parameter @xmath128 to be small , where @xmath129 denote the velocities of the created particles . for pair production very close to the threshold ,
coulomb corrections are expected to be important .
a summary of our results on two - photon pair production for linear and circular polarization including various approximations is displayed in fig.3 .
the function @xmath130 in the intervall @xmath131 for @xmath132 ( linear polarization , solid line ) and @xmath127 ( circular polarization , short - dashed line ) .
the long - dashed and dotted lines show the respective threshold laws for @xmath133 according to eqs.([flim ] ) and ( [ fmu]).,width=226 ]    ( 0,0)(0,0 ) ( -23,-5)@xmath134 ( -153,3)@xmath135 ( -195,50 )    when applying our general formula ( [ sigmagen ] ) to the same parameters as above but for a linearly polarized laser wave ( @xmath132 , @xmath136 , @xmath137 , @xmath117 kev , and @xmath116 ) , we get a total production rate of 290 sec@xmath115 , which is in rather good agreement with the corresponding value of 240 sec@xmath115 found in ref .
the difference is most likely due to small inaccuracies in the 5-dimensional integration along with the numerical evaluation of the generalized bessel functions in ref .
as is seen from fig.3 , the pair production rate for a linearly polarized laser wave is noticeably larger than the corresponding rate for circular laser polarization within the whole frequency range @xmath138 .      in this subsection
we analyse the high - energy behaviour of strong - field pair creation , i.e. , we shall assume @xmath97 but allow for arbitrary value of @xmath17 . let us first consider the circular polarization case and start from eq.([circ2 ] ) .
we exchange the order of integration and first integrate over @xmath139 .
since @xmath97 , we can introduce a splitting parameter @xmath140 with @xmath141 , and divide the @xmath139-integration into ranges from @xmath142 to @xmath140 and from @xmath140 to 1 .
in the first region we have @xmath143 .
substituting @xmath144 we get @xmath145 \bigg\}.\end{aligned}\ ] ] with the help of the identity @xmath146 this is readily integrated to give @xmath147 \nonumber\\ & & \!\!\!\!\!\!\!\!\!\!\!\!\!\!\
! \times \left ( \ln{a\rho\over\epsilon } + c
\right ) + \left ( 4\sin^2\rho - { 1\over\xi^2 } \right)\ln z \bigg\}\end{aligned}\ ] ] where the relation @xmath148 $ ] was used , which is valid for @xmath149 ; @xmath150 denotes euler s constant . in the second range
we have @xmath151 and we can exploit the small - argument behaviour of the neumann function : @xmath152 . then the integration is easily performed : @xmath153 \nonumber\\ & & - { 19\over 3}\sin^2\rho - { 4\over 3}\left(1-{\sin^2\rho\over\rho^2}\right ) \bigg\}.\end{aligned}\ ] ] putting both parts together , @xmath140 drops out as it should and we arrive at @xmath154 \left ( \ln{2\over
a\rho z } - c \right ) \nonumber\\ & & + { z^2\over\xi^2}\ln z - { 19\over 3}\sin^2\rho - { 4\over 3}\left(1-{\sin^2\rho\over\rho^2}\right ) \bigg\}.\end{aligned}\ ] ] applying an integration by parts to one of the terms , this can be rewritten as @xmath155 \left ( \ln{2\over a\rho z } - c \right ) \nonumber\\ & & + \left ( 1-{\sin 2\rho\over 2\rho}\right ) - { 19\over 3}\sin^2\rho - { 7\over 3}\left(1-{\sin^2\rho\over\rho^2}\right ) \bigg\}.\nonumber\\ \\end{aligned}\ ] ] at @xmath16 , the cross section resulting from eq.([rho ] ) in the leading order in an expansion with respect to @xmath26 coincides with eq.([bh ] ) .
the next - to - leading order term for the cross section reads @xmath156.\end{aligned}\ ] ] the correction @xmath157 is negative .
this is due to the fact that it not only contains the process of two - photon pair production but also includes the @xmath26-corrections to one - photon pair creation ( cf .
the latter is negative and has a larger absolute value than the contribution from two - photon pair creation .
now we consider the case @xmath18 , which corresponds to a supercritical laser field strength @xmath19 .
note that , for the very high frequencies ( @xmath97 ) assumed , such extremely large field strengths are far beyond present and near - future technical capabilities .
in this situation , the main contribution to the integral over @xmath89 in eq.([rho ] ) comes from the region @xmath158 .
performing the according expansion of the integrand and taking the integral , we obtain @xmath159.\end{aligned}\ ] ] we notice that the leading logarithm in eq.([sigmahe ] ) has already been derived by yakovlev @xcite .
further , eq.([sigmahe ] ) agrees with the corresponding asymptotic limits found in ref.@xcite for pair creation by a coulomb field and a constant crossed field and in ref.@xcite for pair creation in a coulomb field and a constant homogeneous field .
it is worth stressing that the @xmath17-dependence in eq.([sigmahe ] ) can not simply be expressed in terms of the effective ( laser - dressed ) electron mass @xmath160 @xcite .
for the general case of an elliptically polarized laser wave the expression corresponding to eq.([whe ] ) reads @xmath161 \nonumber\\ & & \times \bigg[{\xi_1 ^ 2+\xi_2 ^ 2\over g_2}(a_1 + 2\sin^2\rho ) - 2\ln{\chi\over 2 } \nonumber\\ & & + { 1\over a_0}\left({g_0\over g_2}-1\right)(a_1 - 3\sin^2\rho)\bigg ] \nonumber\\ & & + { \xi_1 ^ 2+\xi_2 ^ 2\over g_2 }     \bigg[a_1\left(\ln{g_0+g_2\over 2g_2 } - { 4\over 3}\right ) \nonumber\\ & & + 2\sin^2\rho\left ( \ln{g_0+g_2\over 2g_2 } - { 19\over 12 } \right)\bigg ] \nonumber\\ & & -\ln{\chi\over 2}\left(\ln{\chi\over 2}-{5\over 3}\right ) -{1\over
2}l\left({g_2-g_0\over g_2+g_0}\right ) \nonumber\\ & & + { 1\over a_0 g_2}\bigg [ a_1\left(g_0\ln{g_0+g_2\over 2g_2}-{4\over 3}(g_0-g_2)\right ) \nonumber\\ & & -3\sin^2\rho\left(g_0\ln{g_0+g_2\over 2g_2}-{3\over 2}(g_0-g_2)\right)\bigg]\bigg\}\end{aligned}\ ] ] with @xmath162 , @xmath163 , and spencer s function @xmath164 in the limit @xmath165 , we again can expand the integrand for small values of @xmath89 and , within logarithmic accuracy , arrive at the following cross section : @xmath166 here , @xmath167 denotes a hypergeometric function , which smoothly depends on the ellipticity parameter @xmath168 and monotonously decreases from unity for @xmath127 ( circular polarization ) to @xmath169 for @xmath132 ( linear polarization ) .
now we consider the limit @xmath18 .
note that today s most powerful lasers achieve @xmath170 . in this limit
, the main contribution to the integral over @xmath89 in eqs.([circ2 ] ) and ( [ gen2 ] ) comes again from the region @xmath158 .
hence , for the case of circular polarization , we can write @xmath171 \nonumber\\ & & + { 2\over 3}\left [ y_1(\beta)-{1\over u}y_1(\beta u)\right]{(1-t^2)^{3/2}\over t^2 } \nonumber\\ & & + \big [ y_1(\beta)- u y_1(\beta u)\big ] \ln\left({1+\sqrt{1-t^2}\over 1-\sqrt{1-t^2}}\right ) \bigg\},\end{aligned}\ ] ] where now @xmath172 with @xmath173 and @xmath174 . concerning the laser field parameters , the total rate in eq.([circquasi ] ) only depends on the dimensionless ratio @xmath175 , which is the so - called quasiclassical parameter . in the limit @xmath176
we arrive at the production rate @xmath177,\end{aligned}\ ] ] which corresponds to the cross section ( [ sigmahe ] ) obtained in sec.iii.c . in the opposite limit @xmath178 , which corresponds to the existing high - power lasers in the optical or infrared frequency range , the integral in eq.([circquasi ] ) can be evaluated using laplace s method .
the resulting pair creation rate @xmath179 exhibits a tunneling behaviour .
equation([sigmaquasi ] ) is in agreement with corresponding results on pair production in a coulomb field and a constant crossed field @xcite or a coulomb field and a constant homegenous field @xcite .    in a similar way we can derive from eq.([gen2 ] )
the @xmath178 limit for an elliptically polarized laser wave . to this end , we assume @xmath180 and @xmath181 , where @xmath182 for @xmath183 .
the latter condition means , that the laser field is considerably different from a circularly polarized wave
. then we find @xmath184 in comparison with the result for circular polarization in eq.([sigmaquasi ] ) , the production rate in eq.([ratebigeta ] ) is suppressed by an additional factor of @xmath185 .
this suppression has the same nature as in the case of @xmath24 pair production by a real photon in a laser field @xcite .
the reason for the suppression is that , in general , the modulus of the field strength is truely oscillating , while in a circularly polarized wave it has a constant value .
we note that eq.([ratebigeta ] ) can also be derived by averaging eq.([sigmaquasi ] ) over one oscillation cycle .
we have investigated in detail the total probability for electron - positron pair creation by a nuclear coulomb field and an intense laser field . employing the optical theorem in connection with the known polarization operator of a photon in an electromagnetic plane wave , we have derived compact , general expressions for the total pair production rate and found explicit formulae for various intensity and frequency regimes of interest .
we have demonstrated that significant differences exist between the general case of an elliptically polarized laser wave and the special case of circular polarization [ see , e.g. , eqs.([fmu ] ) and ( [ ratebigeta ] ) versus eqs.([flim ] ) and ( [ sigmaquasi ] ) ] .
particular emphasize was placed on the nonlinear process of two - photon pair creation .
let us estimate the feasibility of @xmath24 pair production using an xfel beam @xcite and a beam of relativistic protons .
we will compare the following possible regimes of interaction : 1 ) the one - photon perturbative regime , where @xmath16 and @xmath186 [ see eq.([bh ] ) ] , 2 ) the two - photon perturbative regime , where @xmath16 and @xmath102 [ see eq.([sigmagen ] ) ] , 3 ) the tunneling regime , where @xmath18 and @xmath187 [ see eq.([sigmaquasi ] ) ] , and 4 ) the overcritical - field regime , where @xmath18 and @xmath188 [ see eq.([rateloweta ] ) ] .
we have to take into account that the mentioned formulae are valid in the rest frame of the proton , which moves with a large lorentz - factor @xmath1 . in the rest frame of the proton , the photon frequency is increased with respect to its value in the lab frame : @xmath189 , while the xfel pulse duration is decreased : @xmath190 . here ,
@xmath191 and @xmath192 denote the laser frequency and pulse duration in the lab frame , respectively .
note that , nowadays , the proton beam at desy can be accelerated up to 920 gev , which corresponds to @xmath193 .    for the following estimates
, we will throughout assume an xfel pulse duration of @xmath194 fs .
the xfel beam collides with a bunch of protons moving at @xmath195 , containing 10@xmath196 particles , and having a beam radius of 30 @xmath168 m . in the one - photon perturbative regime , assuming an xfel photon energy of @xmath197 kev , an intensity parameter of @xmath198 , and a beam radius of @xmath199 @xmath168 m , we find that @xmath200 pairs are produced per collision . for a given laser intensity ,
the number of created particles is independent of the laser s polarization state . in the second regime of two - photon pair creation , employing a linearly polarized xfel beam is more favorable .
supposing the collision parameters @xmath201 kev , @xmath202 , and @xmath199 @xmath168 m , the process proceeds far above threshold and one pair can be produced per shot . in the overcritical - field regime ,
circular polarization of the xfel beam leads to higher production yields . in this regime ,
a large @xmath17-value is necessary , which can be achieved by focusing the xfel beam @xcite . for @xmath203 , @xmath197 kev and @xmath204 nm
, one pair can be produced per collision .
here one has to take into account that the xfel pulse will interact only with a very small fraction of the proton bunch . for the assumed parameters ,
the latter has a length of 10 cm and contains a particle density of 10@xmath205 @xmath206 .
thus , we see that , in principle , experimental studies of the pair production process are possible by using an xfel and a relativistic proton beam . employing a non - focused xfel beam ( cf . the so - called design and available parameters in table 1 of ref.@xcite ) and an ultra - relativistic proton beam ,
the perturbative regimes of pair creation can be realized , including nonlinear pair production via two - photon absorption .
the realization of the strong - field regimes with @xmath18 will require focusing of the xfel beam ( cf . the goal parameters in table 1 of ref.@xcite ) . in this case
the interaction volume , i.e. the region of beam overlap , is extremely small which substantially decreases the total production probability . as a result
, the pair creation seems unrealistic in the tunneling regime but , at least in principle , feasible in the overcritical - field regime .    in summary , with the novel sources of spatially coherent x - rays along with the relativistic proton beams of desy , the results on pair production presented in this paper could be tested .
a.i.m . thanks the max - planck - institute for nuclear physics for its kind hospitality during an extended guest researcher appointment when this work was done . u.d.j . acknowledges support from the deutsche forschungsgemeinschaft ( heisenberg program ) .
see , e.g. , e. brezin and c. itzykson , phys .
d * 2 * , 1191 ( 1970 ) ; v. s. popov , zh .
. fiz . * 61 * , 1334 ( 1971 ) [ sov .
jetp * 34 * , 709 ( 1972 ) ] ; yad . fiz . * 19 * , 1140 ( 1974 ) [ sov
* 19 * , 584 ( 1974 ) ] ; pisma zh .
. fiz . * 74 * , 151 ( 2001 ) [ jetp lett .
* 74 * , 133 ( 2001 ) ] ; v. m. mostepanenko and v. m. frolov , yad . fiz . * 19 * , 885 ( 1974 ) [ sov .
* 19 * , 451 ( 1974 ) ] ; a. di piazza , phys . rev .
d * 70 * , 053013 ( 2004 ) .
s. adler , ann . of phys . * 67 * , 599 ( 1971 ) ; v. o. papanyan and v. i. ritus , zh .
. fiz . * 61 * , 2231 ( 1971 ) [ sov .
jetp * 34 * , 1195 ( 1972 ) ] ; a. shabad , nuovo cim . lett . * 3 * , 457 ( 1972 ) ; v. n. baier , v. m. katkov , and v. m. strakhovenko , zh .
. fiz . * 68 * , 405 ( 1975 ) [ sov .
jetp * 41 * , 198 ( 1975 ) ] .


particles which locally strongly interact with their environment polarize this latter and finish up in self - trapped states . on a local level , the environment surrounding the particles can be described in terms of a coherent superposition ( in form of glauber states ) of the elementary bosonic excitations of the bare environment . as a result , new composite entities form , given by the charge carriers and their surrounding clouds of bosonic excitations of the environment .
it is at present a matter of dispute whether in the normal state such composite entities can exist as itinerant bloch - like states or are purely diffusive .
the retarded interaction between the charge carriers and the bosonic excitations of the environment favors the picture of diffusive motion caused by a dephasing of the dynamics of the two constituents .
the question we want to address here is whether this process of dephasing can be blocked , leading to a delocalization when such self - trapped composite entities condense into a _ macroscopic _ coherent quantum state .
we investigate such a possibility in terms of a model of very general form , given by bosonic charge carriers strongly coupled to local bosonic excitations of the environment .
the bosonic charge carriers can , amongst others , represent electron pairs induced by strong electron - phonon or electron - spin fluctuation interaction , electron pairs in double valence fluctuation systems ( stabilized by their ligand environment ) , frenkel excitons ( stabilized by their interaction with the surrounding lattice deformations ) and , in a more remote sense , he@xmath0 atoms in a bath of he@xmath1 with whose excitations they interact .
such bosonic entities can in principle condense into a superfluid state .
our interest here is to show how in such a _ macroscopic _ coherent quantum state the internal structure of the _ local _ coherent quantum state of the bosonic excitations of the environment is modified .
we shall investigate this problem on the basis of a model of itinerant electrons , coupled to some purely local vibrational modes .
for moderately strong and strong coupling , it results in self - trapped composite entities , given by on - site electron pairs , surrounded by local lattice deformations which are described by clouds of phonons in local coherent quantum states . in the intermediate coupling limit these localized entities can be considered as two - particle resonant states of the underlying system of itinerant electrons .
they form out of pairs of such uncorrelated itinerant electrons and ultimately decay into them .
a model hamiltonian which captures such physics is given by the boson - fermion model ( bfm ) @xmath2 \nonumber \\ & & -   \hbar \omega_0 \alpha \sum_i \left ( \rho_i^z + \frac{1}{2 } \right ) ( a_i+a_i^{+ } ) + \hbar \omega_0 \sum_i \left(a^{+}_i a_i + \frac{1}{2}\right ) \ , .
\nonumber \\\end{aligned}\ ] ]    here @xmath3 denote the creation and annihilation operators of hard - core bosons ( characterizing the self - trapped localized electron - pairs ) on some effective sites @xmath4 , which should be understood as of being made up of small clusters of atoms rather than single atomic sites .
such hard - core bosons have spin-@xmath5 commutation relations @xmath6_-=2 \rho_i^z$ ] and @xmath7_+= \delta_{ij}$ ] .
@xmath8 denote the fermionic operators referring to annihilation ( creation ) of itinerant electrons with spin @xmath9 , @xmath10 being the number operator for such fermions . @xmath11
denote annihilation ( creation ) operators of the excitations of the environment , related to the local lattice displacements @xmath12 where @xmath13 is some atomic mass characterizing the effective sites , and @xmath14 is the frequency of the local deformations of the lattice .
the bare hopping integral for the itinerant electrons is given by @xmath15 , the bare electron half bandwidth by @xmath16 ( which will serve as energy unit ) , @xmath17 denoting the coordination number .
@xmath18 denotes the coupling of the hard - core bosons to the deformations of the surrounding medium and @xmath19 the exchange coupling between the bosons and the pairs of itinerant fermions .
the bare boson energy level is denoted by @xmath20 .
the chemical potential @xmath21 is common to fermions and bosons such as to guarantee overall charge conservation .
the model defined by the hamiltonian ( 1 ) was introduced shortly after the proposition of the concept of _ bipolaronic superconductivity _
@xcite ( appropriate to the extreme strong coupling anti - adiabatic regime ) in an attempt to extend this concept to the regime of intermediate electron - lattice deformation coupling .
this model is based on a conjecture , rather than on a derivation from an underlying basic electron - phonon hamiltonian , such as for instance the holstein model@xcite .
we shall in the following develop the ideas which had initially led us to this conjecture .
let us start from a simple holstein model with purely local einstein modes for the lattice vibrations .
the corresponding hamiltonian is given by @xmath22 the meaning of the various operators and coupling constants is the same as in the hamiltonian ( 1 ) , except for the value of @xmath23 which is equal to @xmath24 .
the electron - phonon interaction with coupling strength @xmath23 is taken to be local , in view of the interaction with primarily local phonon modes which favours small polaron formation .
@xmath25 denotes some effective coulomb repulsion - not to be confused with a hubbard - type on - site interaction - on the effective sites , which , we stress again , ought to be considered as being composed of small molecular clusters such as a diatomic molecular unit or similar .
electrons described by this hamiltonian can gain energy either by becoming itinerant and remaining essentially uncoupled to the lattice or , alternatively , getting self - trapped . in the first case
the energy gain is described by the intrinsic band dispersion @xmath26 @xmath27 denoting lattice vectors linking nearest - neighbor sites . for self - trapped polaronic states ,
the local part of the hamiltonian ( [ h ] ) can be transformed by a shift in the phonon variables @xmath28 upon which it acquires the form @xmath29 with the polaron binding energy @xmath30 denoting the gain in energy by the corresponding polaron level shift .
the effective on - site interaction given by @xmath31 controls the possible formation of bipolaronic states @xmath32 when @xmath33 , with the single polaron states lying above such bipolaron states . for @xmath34
we have the inverse situation with the bipolaronic level lying above that of the single polarons .
what we are interested here is the first case .
let us begin by considering the situation of sufficiently strong interaction such that the bipolaronic level lies below the bottom of the unrenormalized electron band , i.e. , for @xmath35 , as illustrated in fig.1
. in that case bipolarons can acquire itinerancy by decaying into virtual states of pairs of itinerant electrons and reassembling subsequently on some neighbouring site . within the usual lang - firsov approximation , neglecting phonon creation as well as annihilation processes during the charge transfer from one site to the other , the corresponding effective bipolaron hopping integral is then given by @xmath36@xcite . within
this approximation one obtains a system of itinerant bipolarons on a lattice which , at low temperatures , can condense into a superfluid state , giving rise to what has been termed _
bipolaronic superconductivity_. this theoretical proposition has ever since remained a matter of dispute , questioning if such a normal - state itinerant bose liquid can be achieved in real materials .
the main arguments against it is that the stability of local bipolarons requires a relatively high value of @xmath37 for realistic values of @xmath38 , @xmath25 and @xmath14 and hence leads to exceedingly small values of the bipolaron bandwidth .
a further hampering factor for the realization of such bipolaronic itinerant band states is related to the retardation effects between the motion of the charge carriers and the lattice deformations which surround them .
this gives rise to dephasing between the two , as shown by exact diagonalization approaches@xcite and dmrg studies@xcite , indicating that incoherent processes in the optical conductivity are dominant .
all these results suggest that coherent bipolaron motion is most unlikely to occur in real systems and it was for that reason that the boson - fermion model was initially proposed@xcite .
the aim was to capture the physics outside the regime of the extreme strong coupling antiadiabatic limit and to describe the intermediate coupling case , where the bipolaron level lies inside the band of the itinerant electrons ( see fig.2 ) and where the superconducting state is controlled by the phase correlations of these resonant bosonic bipolaron states .
quite generally , this model can be considered as a paradigm for cross - over situations . for the specific case
we are interested here , it describes a cross - over between a bcs - type superconductor of cooper pairs in the weak coupling regime and the hypothetical , yet to be experimentally verified , _ bipolaronic superconductivity _ in the strong coupling regime . in a totally different physical situation ,
namely that of a strongly correlated electron system such as described by the hubbard model , it has been recently shown that the intermediate coupling case can be mapped onto such a boson - fermion scenario@xcite .
such an effective boson - fermion hamiltonian was also derived recently for the exchange interaction between spinon singlets of the rvb electron - pairs and holons@xcite .
the question to what extent such a bfm can emerge from a quite general class of interacting fermion systems , has been addressed from the point of view of a bosonization procedure of an intrinsically fermionic system@xcite .
the presently much discussed molecular bose - einstein condensate , involving entangled atoms in squeezed states has also recently been analyzed@xcite in the framework of such a bfm .
a derivation of the bfm for moderately strongly coupled electron - phonon systems , starting from the holstein model , is so far not available . in order to provide the physically intuitive picture which initially had led us to propose this model ,
let us consider the case where the local bipolaron states are stable relative to single polarons .
this eventuality can only be decided on the basis of direct experimental evidence for such local bipolaronic entities .
the classical cases which initially had led to consider such a situation are systems like amorphous semiconductors@xcite , ti@xmath39o@xmath40@xcite , na@xmath41v@xmath42o@xmath43@xcite , wo@xmath44@xcite and double valence fluctuation systems@xcite like pbte(tl ) , just to name a few of them , and all of which represent such bipolarons as _ exclusively _ localized self - trapped states .    when these systems are doped , the electrons responsible for pairing are at the same time responsible for the - albeit diffusive - transport .
it is hence tempting to hypothesize that such bipolaronic states can occur as resonant states inside a band of itinerant electrons . in more complex systems , such as ternary and quaternary compounds , different electrons come into play .
some of them are prone to form bipolarons and others will remain itinerant , but being coupled together via some hybridization term . in systems like that ,
such a boson - fermion scenario imposes itself more directly .
it is the eventuality of bipolaronic resonant states rather than localized states , which was at the origin of our proposal of the boson - fermion scenario@xcite .
the level scheme for the local bipolaronic states inside the band of the bare itinerant electrons for this situation is depicted in fig.2 .
filling up such a system gradually with electrons , it is evident from fig.2 that , as long as the fermi level remains below the level of the local bipolaronic states , the charge carriers will be essentially bare itinerant electrons . but upon further increasing the concentration of charge carriers , the itinerant electrons will eventually overlap in energy with the intrinsically localized bipolaronic level .
the major supposition leading to this bfm is now to assume that this overlap can be simulated by some effective on - site hybridization describing a charge exchange between the bipolarons and pairs of itinerant electrons in the form of @xmath45 \ ; , \ ] ] with the itinerant electrons and the localized bipolarons , eq.(5 ) , treated as commuting fields .
such a supposition is supported by the fact that the spectrum of electrons , moderately strongly coupled to local modes , is known to be given by the bare spectrum for all @xmath46-vectors , except where their energy crosses the level of localized selftrapped electrons . in this regime of energies
the electron spectrum is strongly modified , as shown by exact diagonalization studies@xcite as well as variational methods@xcite , with the spectrum becoming rather flat and approaching the local level of the self - trapped polaronic states .
these are features which are characteristic of systems where itinerant charge carriers are hybridized with localized charge carriers , and this was the center of our initial assumption of an effective hybridization between localized bipolarons and pairs of itinerant electrons , as given by eq.(6 ) .    on a local level the dynamics of a polaronic system , involving two - particle resonant states , is controlled by the atomic limit of such a system , which we thus propose to describe by the following hamiltonian    @xmath47 + \sum_i \hbar \omega_0 \left(a^{+}_ia_i + \frac{1}{2}\right ) \\ & \equiv & ( u-2\mu ) \sum_i
\left(\rho_i^z+\frac{1}{2}\right)+ v\sum_i \left[\rho^+_ic_{i\downarrow}c_{i\uparrow } + \rho_i^- c^+_{i\uparrow}c^+_{i\downarrow}\right]- \hbar \omega_0 \alpha \sum_i \left(\rho_i^z + \frac{1}{2}\right ) \left(a_i+a_i^{+}\right)+\sum_i \hbar
\omega_0 \left(a^{+}_ia_i + \frac{1}{2}\right ) , \nonumber\end{aligned}\ ] ]    with a bare bipolaron level given by @xmath48 .    the equivalent formulation of this hamiltonian [ the second part of eq.(7 ) ]
is obtained after employing an inverse lang - firsov - type transformation @xmath49 with @xmath50 . adding to this
local hamiltonian the term describing the kinetic energy of the itinerant electrons , leads immediately to the bfm hamiltonian given by eq.(1 ) .
the basic physics contained in this bfm model originates from a competition of the local electron pairing , acting as a prerequisite for superconductivity , with the delocalizing effect of the itinerant electrons which favors non - local cooper pairing as the temperature is lowered and the superconducting phase is approached . the change - over from local pairing to cooper pairing
is not only manifest in the electronic properties of such systems but is expected to be accompanied by corresponding lattice deformations , the study of which is the issue of the present paper .    we shall examine this feature for the two following extreme limits :    \i ) the normal - state phase , characterized by the existence of a pseudogap in the density of states ( dos ) of the itinerant electrons , which will be treated here in an approximate fashion within the atomic limit of the bfm ,    \ii ) the superconducting phase , which will be treated on a mean - field level , taking into account the phase locking of the bound electron pairs .    throughout this study
we choose as a representative example a set of parameters presenting a situation close to the fully symmetric case .
this implies an effective bipolaron level ( corresponding to the energy level for the bound pairs ) given by @xmath51 , resulting in a final state bipolarons level lying close to the center of the band of the itinerant electrons , after the boson - phonon coupling has been switched on . we moreover want to consider the case @xmath52 , giving a concentration of electrons @xmath53 and a concentration of bound pairs @xmath54 and choose @xmath55 , @xmath56 and @xmath57 .
let us start from the purely local physics , described by the atomic limit ( @xmath58 ) of the bfm corresponding to the hamiltonian @xmath59 , eq.(7 ) .
this limit has been studied by us previously@xcite in connection with the pseudogap phenomenon and the incoherent background in the electron spectral function , which , in such a scenario , arise from phonon shake - off processes . for the choice of parameters adopted here ,
the physics controlling the change - over from the normal into the superconducting phase can essentially be represented by the one- , two- and three - particle eigenstates of @xmath59 @xmath60 with @xmath61 denoting the @xmath4-particle @xmath62-th eigenstate .
the opening of the pseudogap at a certain temperature @xmath63 shows up in the rapid spectral weight increase , with decreasing temperature below @xmath63 , of the scattering amplitude involving transitions from the lowest energy two - particle bonding state @xmath64 to the lowest energy single - particle states @xmath65 and from the lowest energy three - particle states @xmath66 to @xmath64 .
these contributions to the dos of the itinerant electrons are manifest by the wings ( see fig.3 ) for energies below and above that of the single particle non - bonding contributions , which have spectral weight around the chemical potential ( equal to @xmath67 ) and which arise from transitions not involving any of the two - particle states .
decreasing the temperature below @xmath63 , the spectral weight at @xmath67 decreases rapidly and transfers to the wings of the dos , thus resulting in the opening of a pseudogap . taking into account the dispersion of the itinerant electrons in an approximate fashion , such as within a coherent potential approximation ( cpa ) approach@xcite ,
the main features seen in fig.3 are maintained , simply leading to a smearing of the dos .    ) , as the temperature decreases below @xmath63.,width=302 ]    the two - particle eigenstates of @xmath59 , eq.(9 ) , are obtained by expanding the local oscillator states in terms of the eigenstates @xmath68 of the undeformed harmonic oscillator at a given site @xmath4 . @xmath69 and @xmath70 denote respectively the vacuum states of the itinerant electrons and of the bound electron - pairs on such a site .
the coefficients @xmath71 and @xmath72 are determined by exact diagonalization of this atomic limit problem within a restricted hilbert space of those excited oscillator states .
it is generally sufficient to take into account up to @xmath73 such states in order to get converging results .
the qualitative changes in the lattice deformations expected to occur when going from the pseudogap phase into the superconducting phase are essentially described by the lowest energy one- , two- and three - particle states , which contain the information on the relative weight of the bound electron - pairs ( @xmath74 ) and of the induced pairing in the subsystem of the itinerant electrons ( @xmath75 ) .
the first ones are intrinsically coupled to the local lattice deformations while the second ones acquire such local lattice deformations due to the charge exchange mechanism acting between the two . as we shall see below , the efficiency of this transfer of polaronic features is noticeably different in the normal and in the superconducting phase .
it shows up in the composition @xmath76 of the coherent phonon states describing the local phonon clouds surrounding the charge carriers on a given site @xmath4 , which in the normal state ( atomic limit ) is given by @xmath77 @xmath78 denoting the partition function and @xmath79 the eigenvalues of @xmath59 associated with the eigenstates @xmath80 . upon decreasing the temperature
, @xmath81 shows a redistribution of weight from high to small values of the phonon number @xmath82 [ see fig.4(a ) ] .     in the atomic limit and in the superconducting state for different temperatures @xmath83.,width=302 ]    in the temperature regime of interest here ,
i.e. , from slightly above @xmath63 all the way down to @xmath84 , the contributions to @xmath76 come essentially from the lowest energy one- , two- and three - particle states , eqs.(8 - 10 ) . as illustrated in fig.5 ,
these various contributions have noticeably different features as far as the phonon distributions are concerned .
the one - particle states are totally decoupled from the lattice and hence have a contribution only for @xmath85 , while the three - particle contribution are maximally coupled to the lattice showing a broad peak in this distribution which corresponds to that of a shifted oscillator ground state .
the two - particle bonding state @xmath64 has features of a displaced oscillator which is particularly noticeable at low temperatures where it is the main contribution , together with a much reduced contribution from the two - particle antibonding state @xmath86 . that latter contribution to @xmath76 , considered without thermal weighting factor , clearly shows strong contributions for low values of @xmath82 , which is indicative of weak correlations between the electronic degrees of freedom and those of the lattice .     in the normal state ( full squares ) into contributions coming from the one- , two- and three - particle states ( @xmath87 , @xmath88 and @xmath89 , respectively ) . @xmath90 and @xmath91 represent the contributions coming respectively from the lowest energy two - particle state ( bonding ) and from the first excited two - particle state ( antibonding).,width=264 ]    the onset of the local coherence between bound pairs and of the induced pairing in the itinerant subsystem can be tracked in the pair distribution function ( pdf ) @xmath92 , given in the normal state ( atomic limit ) by @xmath93 which describes the probability of finding a local deformation of size @xmath94 and which can be studied by exafs and pulsed neutron scattering experiments@xcite .
@xmath95 is obtained by expressing the eigenstates of @xmath59 in terms of a real - space representation involving the harmonic oscillator wavefunctions . in fig.6(a )
we illustrate the temperature evolution of the pdf . its decomposition into the contributions arising from    , measured by the dimensionless parameter @xmath96 , in the atomic limit and in the superconducting phase for different @xmath83.,width=264 ]    , @xmath97 and @xmath98 , respectively ) .
@xmath99 and @xmath100 represent the contributions coming respectively from the two lowest bonding and antibonding two - particle eigenstates of @xmath59.,width=302 ]    the one- , two- and three- particle states given in eqs.(8 - 10 ) are illustrated in fig.7 .
the one- and three - particle states contribute to the pdf in form of peaks centered around @xmath101 and around some finite value of @xmath94 , respectively .
this expresses the fact that the one - particle states ( 8) are completely decoupled from the lattice , while the three - particle states ( 10 ) contain a bound electron - pair and thus give a maximal shift in the pdf .
the two - particle bonding and anti - bonding states , eq.(9 ) with @xmath102 , respectively , give contributions to the pdf with two peaks lying very close together , which is a result of the strong correlation between the bound electron pairs and the induced pairs in the itinerant electron system , characterizing the normal - state phase . at high temperatures ,
the pdf is thus given by a very broadened smeared out double - peak structure .
but as the temperature is lowered below @xmath63 , local coherence between the bound electron pairs and the pairs of the itinerant electrons sets in . as a result , the opening of the pseudogap in the local density of states is accompanied by a change - over in the pdf from this very smeared out double - peak structure to a rather well defined single - peak structure .
the latter is essentially due to the lowest energy bonding two - particle eigenstate contributing to a pdf with a maximum at some finite value of @xmath94 .
this is a signature of the strong correlations , which gradually build up as we go below @xmath63 , between the bound pairs and the induced pairing in the itinerant electron subsystem leading to roughly equal weight coming from the coefficients @xmath75 and @xmath74 .
let us now put such single units , as described by @xmath59 , together in an infinite lattice , and let them interact with each other via the hopping term of the itinerant electrons .
when the system undergoes a transition into a superfluid state of the charge carriers , the on - site correlations between bound pairs and pairs of itinerant electrons will be weakened , since the @xmath103-space pairing of the itinerant electrons enters into competition with the strong on - site pairing which controls the physics in the normal state . as a result , upon entering the superconducting phase , a substantial modification of the coherent phonon states of the individual units will occur .
we shall , for a moment , make the assumption that the deformation of the medium surrounding the bosonic charge carriers is not influenced by boson - fermion charge fluctuations . in that case
the coupling of the bound electron - pairs to the lattice deformations can be eliminated , along the lines usually adopted in studies of the polaron problem , leaving us with an effective hamiltonian given by the first four terms of the hamiltonian ( 1 ) and with a renormalized hybridization coupling @xmath104 . such a hamiltonian has been studied in detail within different schemes , such as self consistent diagramatic techniques@xcite , dynamical mean - field theory@xcite and renormalization group procedures@xcite , and gives rise to the following physics . as the temperature is lowered below a certain characteristic temperature @xmath63 , pairing of the itinerant electrons sets in , leading to an emptying out of the single electron spectrum near the chemical potential and thus to the opening of a pseudogap . upon further decreasing the temperature
these electron - pairs become itinerant@xcite , eventually condensing into a macroscopic superfluid ground state of the form @xmath105|0)_i\ ] ] ( @xmath106 being some model dependent parameters ) and exhibiting soundwave - like collective excitations as a consequence of the rigidity of the local phases @xmath107 .    in order to study how such a superfluid state affects the surrounding lattice deformations
, we now investigate this problem on the basis of a mean - field treatment of the full hamiltonian with the coupling of the phonons to the bound pairs explicitly taken into account . to this purpose
, we introduce the following order parameters for the charge carriers@xcite @xmath108 which we assume to be homogeneous in space . the corresponding mean field hamiltonian is given by :    @xmath109 \nonumber \\ h_b & = & ( \delta_b-2\mu ) \sum_i \left ( \rho_i^z + \frac{1}{2 }
\right ) + v x \sum_i [ \rho^+_i \ ; + \rho_i^- ] - \hbar \omega_0 \alpha \sum_i \left ( \rho_i^z + \frac{1}{2 } \right ) ( a_i+a_i^{+ } ) + \hbar \omega_0 \sum_i a^{+}_i a_i \ ; .\end{aligned}\ ] ]    on the level of this approximation , the effect of the charge fluctuations between bound pairs and pairs of itinerant electrons is diminished , being taken into account in a global , spatially independent , fashion .
the itinerancy of the intrinsically localized bound electron pairs , induced by the exchange coupling @xcite is not contained in the present approximation and therefore we can not account for collective soundwave - like excitations and the spatially correlated lattice deformations which result from them .
nevertheless , the present study does permit to investigate local excitations of the lattice , controlled by the local fluctuations of the order parameters to which they are coupled , and which can be related to a number of experiments .
considering itinerant bosons on a deformable lattice with weak boson - phonon coupling@xcite , leads to a hybridization of the bogoliubov modes ( coming from the repulsive interaction between the bosons ) and the einstein - like phonon modes . at long wavelengths this results in a modified sound velocity of the bogoliubov mode , with little spectral weight carried by the phonons , and a moderately renormalized einstein mode , principally of phonon character .
in such a context , what the present mean - field study permits us to do , is to investigate the effect of the superfluid state on such an einstein - like mode . due to the strong coupling between the bound electron pairs and the phonons ,
such an einstein mode will undergo a substantial renormalization , as it will be shown below .
the fermionic part @xmath110 of @xmath111 , being dynamically decoupled from the bosonic part @xmath112 and hence from the phonons , can be cast into the standard form @xmath113 .
the @xmath18 s and @xmath114 s denote the operators of bogoliubov quasi - particles having the dispersion @xmath115 which differs from the bare electron dispersion @xmath116 by showing a gap of size @xmath117 .
the corresponding eigenstates of @xmath110 are given by the standard bcs expression @xmath118| 0 \rangle \ ; .\ ] ] the bosonic part @xmath112 of @xmath111 , showing intrinsic coupling between the bosonic charge carriers and the phonons , is diagonalized in terms of a set of states of the form @xmath119 | 0 ) _
i|n>_i\ ] ] having eigenvalues @xmath120 .
the eigenstates of the mean - field hamiltonian , eq.(15 ) , are thus given by @xmath121    as in the atomic limit , the diagonalization of @xmath112 is carried out in a truncated hilbert space spanned by the eigenstates @xmath122 of the undeformed harmonic oscillator .
we ultimately arrive at the following set of self - consistent equations determining the temperature dependence of the order parameters @xmath123 \ ; , \end{aligned}\ ] ] @xmath124 $ ] denoting the partition function corresponding to @xmath112 .
this set of equations must be solved with the constraint that the total density of bosons and fermions remains fixed , i.e. @xmath125 , where @xmath126 \exp\,\left[-\beta e^{sc}_l(x)\right ] \ ; .\end{aligned}\ ] ]    let us now examine the modification of the composition of the local coherent phonon states on the individual lattice sites ( discussed above on the basis of the atomic limit for the normal state ) when a macroscopic phase - coherent state of the charge carriers emerges and the bound elelctron - pairs condense into a superfluid state . the probability distribution of the phonons in the superconducting phase is then described by @xmath127 \left [ ( u^{sc}_{ln})^2 + ( v^{sc } _ { ln})^2\right ] \ ; .\ ] ] the evaluation of @xmath128 in the superfluid phase follows closely the calculation for @xmath81 in the atomic limit .
in fact , the secular equation determining the coefficients @xmath129 and @xmath130 is formally the same as in the atomic limit , with the states @xmath131 and @xmath132 and @xmath133 , respectively .
the only modifications which then appear in this set of equations are ( @xmath4 ) a renormalized hybridization constant @xmath134 , instead of @xmath19 in the atomic limit , and ( @xmath135 ) an energy level for the diagonal elements involving the states @xmath136 being equal to @xmath137 , instead of @xmath138 for the diagonal elements involving the states @xmath139 .
this renormalization of @xmath19 is important since it reduces the bare @xmath19 by roughly an order of magnitude ( for the set of parameters chosen throughout this paper ) , expressing the fact that the superfluidity leads to a considerable reduction of the on - site correlation between the configurations with bound electron - pairs present and with such pairs being absent .
this has significant consequences on the composition of the local phonon clouds surrounding the individual sites and thus on their corresponding pdfs .    just below the superconducting @xmath140
( @xmath141 0.0067 in our present case ) , @xmath128 is determined by an incoherent superposition of essentially the two lowest eigenstates of the mean - field hamiltonian , eq.(15 ) , one giving rise to a broad peak at some finite @xmath82 and the other to a sharp maximum at @xmath85 . at low temperatures , however , the ground state alone is sufficient to describe @xmath128 .
the contribution to @xmath128 leading to a peak at some finite value of @xmath82 ( associated with the coefficients @xmath142 ) arises from the deformation induced by the presence of the bound electron pairs , while the contribution to @xmath128 leading to a maximum for @xmath85 ( associated with the coefficients @xmath143 ) arises from the absence of such an induced deformation for configurations where no bound pairs are present , with the oscillators being essentially undeformed .
the coherent superposition of deformed and undeformed oscillator states in the ground state is a direct consequence of the macroscopic phase locking of the composite bosonic states in the superconducting phase , resulting from the variational ground state of @xmath112 of the form @xmath144 the operators @xmath145 are defined as @xmath146 and differ significantly from the usual polaron shift operators @xmath147 , characterizing the deformations of a simply shifted oscillator state .
_ macroscopic _ phase locking and homogeneity imply that the coefficients @xmath148 and @xmath149 as well as phases @xmath107 are the same for all sites .
eq.(22 ) is a direct generalization of the superfluid ground state of a system of hard - core bosons on a lattice , which is of the form given in eq.(13 ) , and thus shows explicitly the locking together of the charge @xmath150 and lattice ( @xmath145 ) degrees of freedom on a local level .
( squares ) in the superconducting phase in terms of the contributions @xmath151 coming from the ground state ( triangles ) and @xmath152 coming from the first excited state ( circles ) of @xmath112.,width=264 ]    we shall next address ourselves to the question of how such a _ macroscopic _ coherent quantum state of the phonons can be tracked experimentally and , in particular , what difference there is to be expected with respect to the normal state features of the pdf , examined in the previous section .
the pdf in the superconducting phase is given by @xmath153     in the superconducting phase ( full lines ) into the contributions @xmath154 coming from the ground state ( dashed ) and @xmath155 coming from the first excited state ( pointed ) of @xmath112.,width=264 ]    in fig.[pdf0 ] we compare the pdf in the normal state and in the superconducting phase for different temperatures .
the features examined above on the basis of the composition of @xmath128 are perfectly reflected in the pdf . in
the normal state one observes a single peak at some finite @xmath94 , indicative of well defined shifted oscillator states . in the superconducting phase , on the contrary , one observes a weakly temperature dependent double - peak structure , which , as we can see from fig.[pdf2 ] , arises essentially from two contributions associated with the ground state and the first excited state of @xmath112 . at temperatures just below @xmath140 this double - peak structure results from an incoherent superposition of these two states , corresponding to a situation where localized bipolarons are randomly distributed over the lattice .
the statistical distribution of empty and doubly occupied sites will thus give rise to that double - peak structure .
in the limit of low temperature , on the contrary , the pdf arises @xmath156 from the ground state , with a peak centered near @xmath101 ( being determined by the coefficients @xmath143 ) and a peak centered around a finite value of @xmath94 ( being determined by the coefficients @xmath142 ) .
the appearence of a double - peak structure in the pdf has long since been recognized as a signature of delocalization of the charge carriers in polaronic systems@xcite .
the present study shows that , as the superconducting state sets in and is strengthened upon lowering the temperature , the _ macroscopic _ phase coherence of the system not only involves the charge carriers but also the deformations which surround them .
this is only possible if at @xmath84 these deformations are controlled by a single quantum state , i.e. , the ground state @xmath157 .
the superconductivity - induced delocalization of the self - trapped entities which characterize the normal state , is caused by the locking together of the local lattice deformations in a coherent superposition of deformed and essentially undeformed oscillator states , resulting in the developement of a double - peak structure in the pdf .
at least in the scenario presented here , this result does not bear out any undressings@xcite of the charge carriers , stripping off partially the phonon clouds , as one enters the superconducting phase .
such an effect would show up in a pdf which , in the superconducting phase , would be given by essentially a single - peak structure centered around @xmath101 .
the discussion in the previous section demonstrated the qualitative difference which exists between some basic lattice properties in the normal and the superconducting phase , and largely being controlled by the induced electron pairing in the subsystem of the intrinsically uncorrelated itinerant electrons .
what results is a competition between the local pairing , characterizing the normal - state phase , and the * k*-space pairing in the superconducting phase , expected to be particularly relevant around the superconducting transition where the spatial phase fluctuations play a major role and residual local electron pairing@xcite persists in the superconducting phase .
the qualitative difference in the structure of the pdf in the two phases leads us to expect squeezing effects of the local coherent phonon quantum states which are stronger in the pseudogap phase than in the superconducting phase . associated with it
, one should expect larger fluctuations of the local deformations of the environment and of the momenta of the ions representing it .
these fluctuations are given by @xmath158 with @xmath159    the expectation values in eq.([xp ] ) have to be calculated with respect to the eigenvectors of @xmath59 and @xmath160 for , respectively , for the normal and the superconducting phase .
the results for these fluctuations in the two phases are illustrated in fig.[dxdp ] .
in particular , the product @xmath161 clearly shows that the coherent phonon states in the normal state are considerably squeezed , i.e. , this product is for @xmath162 very close to the lower limit of the heisenberg uncertainty principle , equal to @xmath163 . in the superconducting phase , on the contrary , we observe practically no squeezing effect with decreasing temperature , the product @xmath164 being always bigger than what would be expected in the normal state in this temperature regime .
a possible experimental verification of these features comes from neutron absorption measurements@xcite which measure the kinetic energy @xmath165 of the local lattice vibrations .
indeed , the experimental results shown that , as the temperature is reduced , @xmath166 saturates in the superconducting phase and is distinctly above the values expected in the normal state at the same temperature .
this fact is born out in our study , as can be seen by inspection of fig.[dxdp](b ) .
( @xmath167 ) and with @xmath168 .
( c ) degree of the squeezing of the coherent phonon states as a function of @xmath169.,width=264 ]    finally , an issue which ought to be of crucial interest in studying the lattice properties connected with the onset of the superconducting order or of the electron pairing at @xmath63 in the normal state , is the question of the changes in the phonon spectral properties upon going from one phase to the other .
strong changes in the frequency of specific local modes , induced by strong electron - phonon coupling , have been well established by photoinduced infrared absorption measurements@xcite and been related to strong local displacements of the ions . in this context , the appearence of new additional modes occuring with changes from the normal into the superconducting phase have also been reported , although this is still a matter of debate@xcite .
more established are the strong intensity changes of certain raman active modes and the strong frequency renormalizations and lifetime broadening occuring at such transitions@xcite .
such effects can best be studied by energy loss spectra in inelastic neutron scattering or raman spectroscopy where phonons are emitted .
the cross - section for it is given by @xmath170 with @xmath171        here @xmath172 and @xmath173 are the eigenvalues and eigenvectors of the hamiltonians @xmath59 and @xmath112 , respectively .
we illustrate the results of such phonon emission spectra for the two phases in figs .
[ ph1 ] and [ ph2 ] .
quite generally , one observes a quasi - elastic contribution at frequency zero , which is the standard signature of coherent phonon states .
as the temperature is lowered below either @xmath63 or @xmath140 and the itinerant electrons get correlated by local pairing in the normal state or by cooper pairing in the superconducting phase , the local phonon mode first splits into several modes and finally evolves into a well defined mode which is softer in the pseudogap phase and harder in the superconducting phase .
these features of the phonon spectral properties are expected to hold qualitatively true .
a quantitative link - up between the normal phase and the superconducting phase as the temperature is decreased , would require a full description of the local boson - fermion exchange correlations , incorporating the electron itinerancy in the normal state and allowing for local electron pairing in the superconducting phase .
these are presently largely unresolved fundamental problems in such system where amplitude and phase order occur independently , as in the high - temperature superconductors with its pseudogap phase above @xmath140 where local electron pairing persists in an incoherent fashion .
the present study is an attempt to examine the condensed state of a system of composite particles such as charge carriers surrounded by clouds of bosonic excitations of the neighboring environment to which they are coupled and to see to what extent a macroscopic coherent quantum state can be induced in the system of these bosonic excitations . as a specific example we chose a system of electrons ,
moderately strongly coupled to local lattice deformations , which can be described in terms of resonant local - pair states ( localized bipolarons ) inside a fermi sea of essentially uncorrelated itinerant electrons with which they interact via a charge exchange term .
we find noticeable qualitative differences in the structure of the local coherent quantum states of the phonons in the superconducting phase and in the phase above it which corresponds to a normal state characterized by a pseudogap in the dos of the single - particle spectrum .
several measurable quantities related to static and dynamic lattice properties were examined , such as the pair distribution function , the fluctuations in the lattice positions and the corresponding ion kinetic energies of the atoms making up such a lattice and , finally , the phonon spectral properties given by energy - loss spectroscopy ( measurable by inelastic neutron scattering or raman spectroscopy ) .
the difference between local electron pairing in the normal state and * k*-space pairing in the superconducting phase has noticeable consequences on such measurable quantities , resulting in ( @xmath4 ) a single - peak structure of the pdf in the normal state and a double - peak structure in the superconducting phase , ( @xmath135 ) distinctively larger positional fluctuations and kinetic energy of the atomic clusters ( making up the lattice ) in the superconducting phase than in the normal state and , ( @xmath174 ) associated with it , a more effective squeezing of the phonon coherent quantum states in that latter .
for the superconducting phase , the present mean - field study is restricted to local phase fluctuations of the superconducting state given by a generalization of a phase - locked state of bosons on a lattice , eqs.(17,18,22 ) , in which these bosons correspond to composite particles in form of bipolarons [ see eq.(5 ) ] .
the examination of soundwave - like collective excitations of such a system is planned for some future work . there
, the question arises if the local dynamical deformations do get spatially correlated upon entering the superconducting state , as possibly suggested by ion channeling experiments@xcite .
a possible approach to solve this problem could be to integrate out the fermionic sector of the bfm hamiltonian eq.(1 ) , resulting in a gross - pitaevskii lagrangian@xcite for the effective itinerant bosonic charge carriers coupled to the phonons , followed by a treatment similar to that adopted for bosons on a defromable lattice@xcite .
alexandrov and j. ranninger , phys .
b , * 23 * , 1796 ( 1981 ) .
t. holstein , ann .
( ny ) * 8 * , 343 ( 1959 ) .
e. de mello and j. ranninger , phys .
b * 55 * , 14872 ( 1997 ) . c. zhang , e. jeckelmann and s. white , phys . rev .
b * 60 * , 14092 ( 1999 ) .
j. ranninger and s. robaszkiewicz , physica b * 135 * , 468 ( 1985 ) .
e. altman and a. auerbach , cond - mat/0108087 .
e. kochetov and m. mierzejewski , cond - mat/0204420 .
r. friedberg , t.d .
lee and h.k .
ren , phys .
b * 50 * , 10190 ( 1994 ) .
v.a . yurovski and a. ben - reuven , cond - mat/0205267 .
h. fehske , j. loos and g. wellein , z. phys .
b * 104 * , 619 ( 1997 ) .
w. stephan , phys .
b , * 54 * , 8981 ( 1997 ) .
romero , d.w .
brown and k. lindenberg , phys . rev .
b * 59 * , 13728 ( 1999 ) .
anderson , mater .
* 8 * , 153 ( 1976 ) .
s. lakkis , c. schlenker , b.k .
chakraverty , r. buder and m. marezio , phys . rev .
b * 14 * , 1429 ( 1976 ) .
chakraverty , j.m .
sienko and j. bonnerot , phys .
b * 17 * , 3781 ( 1978 ) .
o.f . schirmer and e. salje , j. phys .
c * 13 * , l1067 ( 1980 ) .
chernick , s.n .
lykov and n.i .
grechko , sov .
solid state * 24 * , 1661 ( 1982 ) .
j. ranninger and a. romano , phys .
lett . * 80 * , 5643 ( 1998 ) .
t. domanski , molecular physics reports * 34 * , 39 ( 2001 ) .
t. egami and s.j.l .
billinge , in _ physical properties of high temperature superconductors _ , edited by d. ginsberg ( world scientific , singapore , 1996 ) vol .
v , p. 265 ; j. rhler _ et al .
_ , in _ workshop in high-@xmath140 superconductivity 1996 : ten years after the discovery _ , edited by e. kaldis
_ et al . _
( kluwer , dordrecht , 1997 ) p. 469 .
j. ranninger , j.m .
robin and m. eschrig , phys .
lett . * 74 * , 4027 ( 1995 ) .
robin , a. romano and j. ranninger , phys .
lett . * 81 * , 2755 ( 1998 ) . t. domanski and j. ranninger , phys .
b * 63 * , 134505 ( 2001 ) . p. devillard and j. ranninger , phys .
lett . * 84 * , 5200 ( 2000 ) .
j. ranninger and l. tripodi , solid state commun . * 112 * , 349 ( 1999 ) .
g. jackeli and j. ranninger , phys .
b * 63 * , 184521 ( 2001 ) ; _ ibid . _ * 64 * , 104513 ( 2001 ) .
shore and l.m .
sander , phys .
b * 7 * , 4537 ( 1973 ) . j. hirsch , phys .
b * 62 * , 14487 ( 2000 ) .
i. kosztin , q. chen , b. janko and k. levin , phys .
b * 58 * , r5936 ( 1998 ) . c. taliani _ et al .
_ , in _ electronic properties of high-@xmath140 superconductors and related compounds _ , edited by h. kuzmany , m mehring and j. fink ( springer - verlag , berlin , series of solid state science , 1990 ) , p. 208
r. j. mcqueeney _
et al . _ ,
phys . rev .
* 82 * , 628 ( 1999 ) .
l. pintschovius and m. braden , phys .
b * 60 * , r15039 ( 1999 ) .
g. ruani and p. ricci , phys .
b * 55 * , 93 ( 1997 ) .
et al . _ ,
b * 59 * , 11495 ( 1999 ) .
et al . _ ,
phys . rev .
lett . * 65 * , 2712 ( 1990 ) .
sharma , t. venkatesan , z.h .
zhang , l.r .
liu , r. chu and w.k .
chu , phys .
* 77 * , 4624 ( 1996 ) .
geschkenbein , l.b .
ioffe and a.i .
larkin , phys .
b * 55 * , 3173 ( 1997 ) .