the phenomenon of quantum interference can be used to control physical and chemical processes .
one approach , the ` @xmath3 ' scheme , uses a two - color light field to interfere @xmath4- and @xmath5-photon transitions @xcite .
interference between one- and two - photon transitions , for example , allows controllable polar asymmetry of photoelectrons in atomic ionization @xcite , controllable dissociation of hd@xmath6 @xcite , and controllable photocurrent injection in unbiased solids due to free carrier absorption @xcite , impurity - band transitions @xcite , and quantum well bound - continuum intersubband transitions @xcite . in biased asymmetric semiconductor double wells , ` @xmath7 ' interference allows control of carrier population and thz emission @xcite .
interband ` @xmath7 ' interference in unbiased semiconductors , which is our interest here , allows independent control of electrical current injection @xcite and spin current injection @xcite .
furthermore , in noncentrosymmetric semiconductors , it allows independent control of carrier populations ( i.e. , absorption ) @xcite and carrier spin polarization @xcite . in each scenario ,
the experimenter can control the interference by adjusting the phases of the two colors .
the controllable optical phases are not the only source of phase between the transition amplitudes . in general , there is also a material - dependent intrinsic phase @xcite .
phenomenologically , the intrinsic phase appears as a phase shift in the dependence of the process on a relative phase parameter of the optical fields .
additionally , selectivity between two processes is possible when their intrinsic phases differ @xcite ; for example , ` 1 + 3 ' experiments on diatomic molecules have controlled the branching ratio of ionization and dissociation channels @xcite .
the intrinsic phase can be strongly frequency dependent near resonances @xcite , and the hope that it might serve as a new spectroscopic observable @xcite has led to efforts to understand its physical origin .    whereas a resonance is necessary for a phase shift to a ` 1 + 3 ' process @xcite , it is not necessary for a phase shift to a ` 1 + 2 ' process .
for example , an intrinsic phase in the ` 1 + 2 ' photoionization of atoms is predicted from the simple model of a delta function potential @xcite .    nevertheless , microscopic theories for the interband ` @xmath7 ' processes in bulk semiconductors have until now predicted trivial intrinsic phases @xcite .
the photocurrent , for example , was predicted to be proportional to the sine of the relative phase parameter for all final energies @xcite .
each of these theories use the independent particle approximation , in which the coulomb attraction between the injected electron and hole is neglected .
that approximation is expected to be good for final energies well above the band gap , since in this case the electron and hole travel quickly away from each other .
however , close to the band gap , one generally expects to see deviations from the independent particle approximation . in the one - photon absorption spectrum , for instance
, it is well known that the electron - hole attraction is responsible for exciton lines below the band gap , and for an enhancement of the absorption above the gap known as _
sommerfeld _ or _
coulomb enhancement _
@xcite .
the effect of the electron - hole attraction on one - photon absorption has been studied with various degrees of sophistication . on the one hand , modern _ ab initio _ calculations that include coulomb effects
have recently given very good quantitative agreement with experimental spectra @xcite , although at the cost of significant computational overhead .
on the other hand , simple models of wannier excitons can describe coulomb effects near the band edge of many direct gap semiconductors .
these excitonic effects have long been understood qualitatively on the basis of a simple two - band model in the effective mass approximation @xcite , which is even quantitatively accurate for typical semiconductors @xcite .
excitonic effects on nonlinear optical properties of bulk semiconductors have also been studied in the effective mass , wannier exciton approximation @xcite and only recently with _ ab initio _
methods @xcite .
the two - photon absorption spectrum shows a different set of exciton lines and a coulomb enhancement that is weaker than its one - photon counterpart .
one- and two - photon absorption spectra have been measured sufficiently often that excitonic effects on them are well established .
in contrast , semiconductor ` 1 + 2 ' interference experiments have been done typically at only a single energy and typically many exciton binding energies above the band - gap . moreover
, these initial experiments lacked an absolute calibration of the relative optical phase , and thus were insensitive to the intrinsic phase .
such a calibration is possible @xcite , and could be used to verify the predictions we present here .
a nontrivial intrinsic phase would have implications for the use of ` @xmath7 ' current injection as an absolute measurement of the carrier - envelope phase of an ultrashort optical pulse @xcite .
in the present work , we extend the theory of ` 1 + 2 ' coherent control of bulk semiconductor interband transitions beyond the independent particle approximation , employing a set of approximations that are valid close to the @xmath8 point of a direct gap bulk semiconductor .
our investigation is limited to a perturbative treatment in the fields . in this limit of low photoinjected carrier density
, the only inter - particle interaction of importance is that between a single electron and hole . we show that , due to the electron - hole attraction , a nontrivial phase shift does in fact occur in the control of current and spin current , but not in the control of carrier population or spin polarization .
the intrinsic phase can be understood in terms of partial wave phase shifts due to the coulomb potential between electron and hole .
in addition , we find an enhancement of each process , and relate it to the coulomb enhancements of one- and two - photon absorption .    in the next section
, we establish notation necessary to describe the ` 1 + 2 ' processes in terms of one- and two - photon transition amplitudes . in section [ sec : model ] , we present the microscopic model .
the transition amplitudes are worked out in section [ sec : transition_amplitudes ] .
the final expressions for the ` 1 + 2 ' effects are given in section [ sec : results ] , and numerical results for gaas are presented . in section [ sec : discussion ] , further understanding of the enhancement and intrinsic phase is discussed , and we examine the ratios often used as figures of merit for ` 1 + 2 ' effects .
intermediate state coulomb enhancement is examined in appendix [ appx_intstateenhancement ] . for reference ,
the current injection tensor is worked out for parabolic bands in appendix [ appx_etafree ] .
the rate of photocurrent injection into an unbiased bulk semiconductor by a two color light field @xmath9 can be written @xmath10 where superscripts denote cartesian components and repeated indices are to be summed over @xcite .
the fourth rank tensor @xmath11 , called the current injection tensor , describes the material response .
it is purely imaginary in the independent particle approximation @xcite , but can be complex in general .
we define the intrinsic phase @xmath12 of the component @xmath13 as @xmath14 so that it is zero or @xmath15 in the independent particle approximation .
when the electron - hole interaction is included in the set of approximations used below , all the components of @xmath16 have the same phase .
that is , @xmath17 the intrinsic phase @xmath18 appears as a phase shift in the dependence of the current injection on the phase of the optical fields . for co -
linearly polarized fields ( @xmath19 and @xmath20 ) , for example , the current injection is @xmath21 we are ignoring scattering processes through which the current would relax to a steady state value under continuous illumination , or would decay to zero following pulsed excitation .
the current injection discussed here can be used as a source term in hydrodynamic equations that treat the scattering phenomenologically @xcite or in microscopic transport equations @xcite .
coulomb effects other than the excitonic effects we consider here play a role in scattering , especially at high densities of excited carriers .
such coulomb effects are outside the scope of this paper
.    the current injection tensor can be written in terms of one- and two - photon transition amplitudes .
we take as the initial state a clean , cold semiconductor . in a fermi s golden rule calculation for the ballistic current
, light produces transitions to final states @xmath22 with velocity @xmath23 , probability amplitude @xmath24 , and energy @xmath25 above the ground state .
the final states will be specified in detail in the next section ; here the label @xmath4 represents the set of quantum numbers for either an interacting or independent electron - hole pair .
thus , @xmath26 where @xmath27 is the electron charge ( negative ) , @xmath28 is a normalization volume , and @xmath29 is the amplitude for an _
i_-photon transition .
the @xmath30 take the form @xmath31 where the vector @xmath32 and second rank tensor @xmath33 depend only on properties of the material .
we consider them in detail in section [ sec : transition_amplitudes ] .    in the last equation of ( [ currentmicro ] ) ,
the first term is ( one - photon ) one - color current injection , the second term describes two - photon one - color current injection , while the interference term describes the ` 1 + 2 ' current . because of their different dependencies on the electric field amplitudes , these three terms can in principle be separated experimentally .
but the first two terms vanish for centrosymmetric crystals , and the first vanishes even for zincblende crystals ; we neglect them here .
excitonic effects on the first term were studied by shelest and ntin @xcite .
the third term survives in all materials . comparing its expression with the phenomenological form ( [ current phenomenological ] )
, we have @xmath34    even if scattering from impurities and phonons is neglected , the injection current described by ( [ currentmicro ] ) does not capture the full current density ; there are also optical rectification and ` shift ' contributions to the current @xcite .
the one - color varieties of these have been studied in some detail @xcite , but the ` 1 + 2 ' varieties have not . however
, the different time dependencies of the three current contributions allows for their separate examination experimentally , at least in principle , and rough order - of - magnitude estimates indicate that typically the injection current will be the largest ; it is the only contribution to the current we treat here .    as in ( [ currentmicro ] )
, the carrier injection rate can be written as @xmath35 where @xmath4 is the number density of injected electron - hole pairs .
the first two terms in ( [ carrier injection rate ] ) are the usual one- and two - photon absorption rates , which we denote by @xmath36 and @xmath37 respectively . in terms of one- and two - photon coefficients @xmath38 and @xmath39 that depend only on the properties of the material , they can be written as @xmath40 and @xmath41 @xcite . from ( [ carrier injection rate ] ) , @xmath42 the third term in ( [ carrier injection rate ] ) , denoted @xmath43 , allows population control as discussed and observed by fraser _
et al _ @xcite .
it can be written in terms of a third rank tensor @xmath44 as @xcite @xmath45 where @xmath46 it is purely real in the independent particle approximation @xcite .
expressions such as ( [ currentmicro ] ) and ( [ carrier injection rate ] ) can also be written for carrier spin polarization and spin current @xcite .
the interference terms of these describe ` 1 + 2 ' spin control @xcite and ` 1 + 2 ' spin current injection @xcite , which can be written in terms of material response pseudotensors @xmath47 @xcite and @xmath48 @xcite , respectively . in the independent particle approximation , @xmath47 is purely imaginary @xcite , while @xmath48 is purely real @xcite .
the phases of the material response tensors @xmath16 , @xmath49 , @xmath50 , and @xmath51 are related to the phases of the one- and two - photon matrix elements , @xmath52 and @xmath33 .
the one- and two - photon matrix elements also appear , respectively , in the one- and two - photon absorption coefficients , as can be seen from ( [ one_photon_schematic ] ) and ( [ two photon schematic ] ) .
there have been many theoretical investigations of one- and two - photon absorption near the direct gap of bulk semiconductors that include excitonic effects @xcite .
however , since one- and two - photon absorption are insensitive to the phases of the transition amplitudes , those calculations took no care to get the phases of the transition amplitudes correct . in the next two sections we find the transition amplitudes with the correct phases including excitonic effects
we first review the two - band effective mass model of wannier excitons ; the two bands are nondegenerate conduction and valence bands that are parabolic and isotropic with a direct gap @xmath53 at @xmath54 ( the @xmath8 point ) @xcite .
it has been used to study excitonic effects on one - photon absorption @xcite , two - photon absorption @xcite , and other nonlinear optical processes @xcite .
we subsequently describe a generalization that accounts for degeneracy and multiple bands .
it has been used for two - photon absorption @xcite , and has been implied whenever two - band results have been applied to actual semiconductors .    the total hamiltonian of the system can be written in the form @xmath55 . here
, @xmath56 is the field - free hamiltonian made up of the single - particle part @xmath57 and the part due to the coulomb interaction between carriers @xmath58 . using the minimal coupling hamiltonian
, the optical perturbation takes the form @xmath59 , where @xmath60 is the vector potential associated with the maxwell electric field and @xmath61 is the velocity operator associated with @xmath62 . in the long wavelength limit , the position dependence of @xmath60 is neglected , and thus the second term in @xmath63 may be neglected , since it can simply be absorbed in an overall time - dependent phase of the full system ket and hence can not cause any transitions between states of @xmath62 .
many approximate approaches to band structure calculation , including most pseudopotentials , and the truncation to a finite number of bands , implicitly assume an underlying field - free hamiltonian that is nonlocal ; there is then a correction to the interaction hamiltonian @xmath63 in the velocity gauge @xcite .
however , we neglect such nonlocal corrections , as has been the practice in previous calculations of coherent control effects @xcite .
the initial state is the `` vacuum '' @xmath64 ; it corresponds to completely filled valence bands and empty conduction bands .
if the coulomb interaction were neglected in a two - band model consisting of valence ( @xmath65 ) and conduction ( @xmath66 ) bands , the final states would be of the form @xmath67 , where the operator @xmath68 creates an electron in an eigenstate of @xmath57 , a bloch state @xmath69 with band index @xmath4 and wavevector @xmath0 .
the photon momentum has been neglected , consistent with the long wavelength approximation .
the coulomb interaction couples states at different @xmath0 ; thus including @xmath58 the final states are of the form @xmath70 where @xmath71 labels the state ; its physical meaning is given below .
in the effective mass wannier exciton approximation , the fourier transform @xmath72 which is the wavefunction for the relative coordinate between electron and hole , is a hydrogenic wavefunction satisfying @xmath73 where @xmath74 is the reduced mass in terms of the ( positive ) conduction and valence band effective masses , and @xmath75 is the coulomb potential , @xmath76 , screened by the static dielectric constant @xmath77 @xcite .
the state has energy @xmath78 we choose the states to be normalized over the volume @xmath28 by @xmath79 and @xmath80 ; as a result @xmath81 is unitless , having the normalization @xmath82 .
our focus in this paper is on the unbound solutions to ( [ exciton_schrodinger ] ) ; bound exciton states lack relative velocity between the electron and hole , and hence do not contribute to the ballistic current or spin current . for a fermi s golden rule calculation of the current or spin current , the unbound state must behave asymptotically like an outgoing plane wave in the relative coordinate between electron and hole ; @xmath83 is the wavevector of the outgoing plane wave . specifically , we must use `` ionization states '' rather than scattering states @xcite , as was done for atomic ` 1 + 2 ' ionization @xcite .
they are related by @xmath84 ^{*}$ ] @xcite .
calculations of one- or two - photon absorption are insensitive to an error in this choice of boundary condition , but the present calculation is not , since it is sensitive to the relative phase of the transition amplitudes .    the ionization state wavefunctions that solve ( [ exciton_schrodinger ] ) can be expressed as a sum over partial waves , @xmath85 where @xmath86 is the exciton bohr radius , and @xmath87 and @xmath88 mean @xmath89 and @xmath90@xcite .
the @xmath91 are legendre polynomials , @xmath92 is a confluent hypergeometric function , and @xmath93 is the gamma function .    such a two - band model of wannier excitons is useful for the description of many optical properties . however , near the band gap at the @xmath8 point of a typical zincblende semiconductor there are , counting spin , eight bands : two each of conduction ( @xmath66 ) , heavy hole ( @xmath94 ) , light hole ( @xmath95 ) and split - off hole ( @xmath96 ) .
other bands , especially the next higher conduction bands , can also be important for some processes , especially for population and spin control .
the existence of multiple bands and band degeneracy modifies the exciton hamiltonian , i.e. , the operator acting on @xmath81 in the left side of ( [ exciton_schrodinger ] ) . in the effective mass approximation , using a basis of @xmath8 point states , the kinetic part of the wannier exciton hamiltonian has a matrix structure @xcite .
even though this approximation neglects band warping , nonparabolicity , and inversion asymmetry , the hamiltonian lacks analytic eigenstates @xcite .
this is essentially due to the degeneracy of the @xmath94 and @xmath95 bands at the @xmath93 point . as a result of the difference between @xmath97 and @xmath98
there is ` envelope - hole coupling , ' @xcite which is a spin - orbit - like coupling between the orbital angular momentum of the exciton envelope function and the total angular momentum of the valence band @xmath8 point bloch functions @xcite .
baldereschi and lipari split the effective mass hamiltonian into a sum of terms based on symmetry , and showed that in a spherical approximation the envelope - hole coupling could be treated as a perturbation to the diagonal part , which has analytic , hydrogenic eigenstates @xcite . in order to extract the main physics , while preserving the simplicity of the two - band model , we neglect envelope - hole coupling entirely . in this approximation , ( [ exciton_schrodinger ] )
remains valid for each conduction - valence band pair , however one must use ` average ' effective masses for degenerate bands . specifically , the effective mass of the valence bands @xmath94 , @xmath95 , and @xmath96 is @xmath99 , where @xmath5 is the free electron mass , and @xmath100 is one of the luttinger parameters @xcite .
the upper conduction bands have a different average effective mass .
note that @xmath101 is independent of @xmath66 and @xmath65 within the set of bands @xmath102 .
the effect of envelope - hole coupling has been studied for exciton bound states @xcite , but not for optical processes involving unbound excitons in the continuum .    even within this model ,
the presence of multiple bands provides two types of terms in the sum over intermediate states in the two - photon amplitude : two - band terms , in which the intermediate and final states are in the same exciton series [ i.e. , two states of the form ( [ final state ] ) with the same @xmath66 and @xmath65 ] , and three - band terms , in which the intermediate and final states are in different series .
three - band terms are important for some processes but not for others . for current control
, three - band terms are important for cross - linearly polarized fields @xcite , and for spin current control , they are important for the spin current due to co - linearly polarized fields @xcite .
three - band terms are essential for population and spin control @xcite .
the velocity matrix elements involving the state @xmath104 are @xmath105 @xmath106 where @xmath107 is the velocity matrix element between bloch states .
in the independent particle approximation , the transition amplitudes are @xmath108 and @xmath109 , where @xmath110    with excitonic effects included , using the perturbation @xmath111 to second order gives the transition amplitudes @xmath112 and @xmath113 where the sum over intermediate states is over both bound and free excitons .
the two - photon transition amplitude is more difficult to deal with due to the sum over intermediate states ; however , in our set of approximations it has been done exactly @xcite .    in order to proceed analytically
, it is common to use ( [ vac_pairvelocity ] ) and ( [ twotwovelocity ] ) , and then make an expansion in @xmath0 of the velocity matrix elements @xmath114 about the @xmath8 point @xcite .
however , due to the degeneracy at the @xmath8 point , the coefficients of such an expansion can depend on the direction of @xmath0 @xcite . to proceed , we note that wannier excitons have large spatial extent and hence only a small region of wavevectors is important for them , i.e. , @xmath115 is peaked in the region of @xmath0 near @xmath83 @xcite .
this is especially true for final states with energies above the band gap .
thus , we expand @xmath114 about the @xmath93 point , approached in the direction @xmath116 , @xmath117 where @xmath118 and @xmath119 .
optical transitions due to the first term in ( [ bloch_v_expansion ] ) are called ` allowed ' , while those due to the second term are called ` forbidden ' .
we restrict ourselves to materials for which the ` allowed ' valence to conduction band transition does not vanish . keeping only the ` allowed ' term in ( [ vac_pairvelocity]),@xcite @xmath120 for the intravalence and intraconduction band transitions , the first two terms of ( [ bloch_v_expansion ] ) in ( [ twotwovelocity ] ) give @xcite @xmath121 \int \frac{d^{3}r}{l^{3}}\left ( \psi _ { cv}^{\bm{\kappa } } \left ( \mathbf{r}\right ) \right ) ^{*}\psi _ { c^{\prime } v^{\prime } } ^{\bm{\kappa } ^{\prime } } ( \mathbf{r } ) \\ & -\left [ \delta _ { v , v^{\prime } } \bm{\nabla}_{\mathbf{k}}v_{cc^{\prime } } ^{i}\left ( \mathbf{\hat{\bm{\kappa } } } \right ) -\delta _ { c , c^{\prime } } \bm{\nabla } _ { \mathbf{k}}v_{v^{\prime } v}^{i}\left ( \mathbf{\hat{\bm{\kappa } } } \right ) \right ] i\int \frac{d^{3}r}{l^{3}}\left ( \psi _ { cv}^{\bm{\kappa } } \left ( \mathbf{r}\right ) \right ) ^{*}\bm{\nabla } \psi _ { c^{\prime } v^{\prime } } ^{\bm{\kappa } ^{\prime } } ( \mathbf{r } ) . \end{split}\ ] ] in particular @xcite , @xmath122    for ge and for simple models of zincblende semiconductors that neglect lack of inversion symmetry , the first term in ( [ v_intra ] ) always vanishes .
this means that there are only ` allowed - forbidden ' two - photon transitions ( the interband transition is allowed , while the intraband transition is forbidden ) .
when the first term is nonvanishing , there are ` allowed - allowed ' two - photon transitions . in principle
, for materials that lack a center of inversion , there is also a small contribution to the ` allowed - forbidden ' two - photon transition from the first term in ( [ v_intra ] ) and the term in @xmath123 that comes from the second term in ( [ bloch_v_expansion ] ) ; we neglect it in what follows , but note that when compared to the dominant ` allowed - forbidden ' contribution that we consider here , it has a different coulomb enhancement but the same intrinsic phase ( see eq .  2.32 of rustagi _ et al_.@xcite ) .
we write @xmath124 , and discuss the ` allowed - forbidden ' and ` allowed - allowed ' transitions separately .    using ( [ ionization state ] ) , ( [ allowed_interband_v ] ) , and ( [ one photon rate ] ) , the one - photon transition amplitude is @xcite @xmath125 where only the ` allowed ' term is kept in @xmath126 .
the transition is to an @xmath127-wave .
the one - photon absorption coefficient is proportional to the norm of @xmath128 [ see ( [ carrier injection rate ] ) and ( [ one_photon_schematic])])@xcite .    for ` allowed - forbidden ' two - photon transitions , substituting ( [ allowed_interband_v ] ) and the second term of ( [ v_intra ] ) into ( [ two photon rate ] ) , @xmath129 where @xmath130 and we have used the coulomb green function , @xmath131 which is known analytically @xcite . in particular , @xmath132 where we define @xmath133 @xmath134 is the exciton binding energy , and
@xmath135 is a whittaker function with the integral representation @xmath136 since the green function depends only on the magnitude of @xmath137 , only the _ p_-wave of the final state survives the integral in eq .
( [ w_defined ] ) over the angles of @xmath137 , @xmath138 .
the integral over @xmath88 can be done using @xcite @xmath139 the final result is @xmath140 where only the ` allowed - forbidden ' term is kept in @xmath141 , and @xmath142    for ` allowed - allowed ' two - photon transitions , substituting ( [ allowed_interband_v ] ) and the first term of ( [ v_intra ] ) into ( [ two photon rate ] ) , @xmath143 \mathbf{e}_{\omega } \cdot \mathbf{v}_{c^{\prime } v^{\prime } } \left ( \mathbf{\hat{\bm{\kappa } } } \right ) \\ & \times \hbar \int d^{3}r\left ( \psi _ { cv}^{\bm{\kappa } } \left ( \mathbf{r}\right ) \right ) ^{*}g_{c^{\prime } v^{\prime } } \left ( \mathbf{r},\mathbf{0};\hbar \omega -e_{c^{\prime } v^{\prime } } ^{g}\right ) .
\end{split}\ ] ] since @xmath144 depends only on the magnitude of @xmath137 [ eq . ( [ coulombgreenfunction ] ) ] , only the @xmath127 part of the final state will survive the integration over angles of @xmath137
. again we use ( [ whittaker ] ) for the whittaker function . the integral over the magnitude of @xmath137
can be done using an identity obtained by taking a derivative with respect to @xmath145 of both sides of ( [ integral identity ] ) . finally , @xmath146 where only the ` allowed - allowed ' term is kept in @xmath141 , and @xmath147 this agrees with eq .  2.28 of rustagi @xcite , but note that we have defined @xmath148 , where @xmath149 is given , with a typographical error , in eq .  2.25 of that paper .    the factors @xmath150 and @xmath151 , which appear in ( [ n_af ] ) and ( [ n_aa ] ) are the enhancements due to the coulomb interaction in the intermediate states ; they are discussed further in appendix [ appx_intstateenhancement ] .
the one- and two - photon transition amplitudes were presented in the previous section on the basis of an expansion in @xmath0 of the bloch state velocity matrix elements .
the ` allowed ' one - photon transition amplitude @xmath152 is in ( [ one photon allowed ] ) , the ` allowed - forbidden ' two - photon transition amplitude is in ( [ two photon a - f ] ) , and the ` allowed - allowed ' two - photon transition amplitude is in ( [ two photon a - a final ] ) . from them , @xmath153 and @xmath154 may be extracted by comparison with the definitions in ( [ d1definition ] ) and ( [ d2definition ] ) .      the ` 1 + 2 ' current injection is dominated by interference of ` allowed ' one - photon transitions and ` allowed - forbidden ' two - photon transitions @xcite . substituting these into ( [ etaschematic ] ) , and using the gamma function identities @xmath155 and @xmath156 yields our final result for the current injection tensor @xmath157 where @xmath158 @xmath159 and @xmath160 with @xmath161 and @xmath162 where @xmath163 and @xmath164 note that only the ` allowed ' part of ( [ d1free ] ) and the ` allowed - forbidden ' part of ( [ d2free ] ) should be retained for a consistent solution .
we have written ( [ eta_result ] ) to separate the parts due to the electron - hole interaction . in the independent particle approximation ,
the current injection tensor @xmath165 is @xcite @xmath166 it is evaluated for parabolic bands in appendix [ appx_etafree ] .    ) ] and without [ eq.([eta_result_free ] ) ] excitonic effects .
the grey dotted and dash - dotted lines are based on a parabolic band calculation of @xmath167 that only includes two - band terms ; the dotted line includes excitonic effects , while the dash - dotted line does not .
the black solid and dashed lines are based on @xmath167 calculated with a nonperturbative solution of the @xmath168 @xmath169 hamiltonian ; the solid line includes excitonic effects , while the dashed line does not .
[ fig_eta],width=307 ]    for gaas , we present in fig .
[ fig_eta ] the magnitude of @xmath170 , based on @xmath167 calculated by two methods .
the first method , described in appendix [ appx_etafree ] , uses isotropic parabolic bands and includes only two - band terms ; it uses effective mass ratios for conduction , heavy hole , light hole , and split - off bands of @xmath171 , @xmath172 , @xmath173 , and @xmath174 respectively , @xmath175ev , the fundamental band gap @xmath176 is 1.519ev , and valence band spin - orbit splitting is 0.341ev @xcite .
the second method solves the @xmath177 @xmath169 hamiltonian including remote band effects , but in a spherical approximation with warping and spin - splitting neglected by replacing @xmath178 and @xmath179 with @xmath180 @xcite ; the calculation is nonperturbative in @xmath0 ( hence it includes band nonparabolicity ) and it includes both two- and three - band terms in the two - photon amplitude @xcite . the solid and dotted lines in fig.[fig_eta ] are calculated with ( [ eta_result ] ) , and hence include excitonic effects ; the coulomb enhancement part of the calculation uses @xmath181mev @xcite and the band parameters listed above . note that the solid black line in fig.[fig_eta ] is inconsistent in the sense that the coulomb enhancement is based on an expansion in @xmath0 , whereas the free - particle result that it enhances is nonperturbative in @xmath0 ; nevertheless , such an approach has given good agreement with experiments for one- and two - photon absorption @xcite .
the coulomb enhancement of @xmath16 can clearly be seen in fig.[fig_eta ] .
there is a kink in each curve at excess photon energy 341 mev corresponding to the onset of transitions from the @xmath96 band . at higher energies ,
the coulomb enhancement of @xmath96 transitions is larger than the coulomb enhancements of @xmath94 and @xmath95 transitions , since the former transitions are to conduction band states with lower energy .
hence , the kink in @xmath16 is enhanced by excitonic effects .    ) in gaas due to excitonic effects .
the solid line is calculated with eqs .
( [ delta explicitly defined ] ) and ( [ eta_result ] ) , and the dotted line is calculated with eq.([phase shift ] ) .
the inset is eq .
( [ phase shift ] ) plotted in scaled units.,width=307 ]    we extract the intrinsic phase of @xmath170 using ( [ delta explicitly defined ] ) .
the solid line in fig.[fig_delta ] is the intrinsic phase of @xmath182 calculated for gaas with the nonperturbative @xmath168 @xmath169 hamiltonian band model ; the result for the parabolic band model is nearly identical .
since we have used a spherical exciton model , the intrinsic phase is the same for all components of @xmath183 .
the intrinsic phase has its maximum value of @xmath1 at the band edge , and goes to zero as the light frequency increases .
the decrease is smooth except for a small kink at the onset of transitions from the @xmath96 band .
in fact , for excess photon energies less than the split - off energy , the intrinsic phase has the simple analytic form @xmath184 equation ( [ phase shift ] ) is plotted as the dotted line in fig.[fig_delta ] ; compared to the solid line , it is identical below the onset of @xmath96 transitions , and it makes a good approximation above the the onset of @xmath96 transitions . since ( [ phase shift ] ) only depends on the excess photon energy scaled by the exciton binding energy , we plot it as a function of this scaled energy in the inset of fig .
[ fig_delta ] ; it is useful for finding the intrinsic phase of materials other than gaas .    in @xmath11
( [ eta_result ] ) ] , the two- and three - band terms have different intermediate state coulomb enhancement @xmath185 . for many materials , however , @xmath185 is approximately equal for all the terms @xmath186 that contribute significantly to the total @xmath165 , as shown in appendix [ appx_intstateenhancement ] for gaas .
thus , at photon energies for which transitions from the heavy- and light - hole bands dominate @xmath187 , the coulomb enhancement becomes approximately independent of the sum over bands and we can make the simplification @xmath188 where the intrinsic phase is given by ( [ phase shift ] ) , and @xmath189     [ eq .  ( [ enhancement ] ) ] with @xmath190 , and the dotted line , applicable to carrier population and spin control is @xmath191 [ eq.([enhancement_aa ] ) ] with @xmath192,width=307 ]    the coulomb enhancement factor @xmath193 is plotted in fig .
[ fig_f ] with the approximation that @xmath190 ( see appendix [ appx_intstateenhancement ] ) .      the ` 1 + 2 ' carrier population control is dominated by interference of ` allowed ' one - photon transitions and ` allowed - allowed ' two - photon transitions @xcite .
substituting these into ( [ xischematic ] ) , and using the gamma function identity ( [ gamma identity ] ) , we find @xmath194 where @xmath195 and @xmath196 and @xmath197 are given by ( [ d2free ] ) and ( [ d1free ] ) .
note that only the ` allowed ' part of ( [ d1free ] ) and the ` allowed - allowed ' part of ( [ d2free ] ) should be retained for a consistent solution . in the independent particle approximation , @xmath198
thus , population control has a coulomb enhancement due to excitonic effects , but no phase shift .
note that ( [ population control result ] ) gives the population control tensor at final energies above the band edge .
there can also be population control of bound excitons when both one- and two - photon transitions are to the same excitonic state .
this can occur , for example , at @xmath127 excitons due to allowed - allowed two - photon transitions @xcite interfering with allowed one - photon transitions .    if @xmath199 is approximately the same for all the terms that significantly contribute to @xmath200 , then , at photon energies for which transitions from the heavy and light hole bands dominate @xmath49 , the coulomb enhancement becomes approximately independent of the sum over bands , and we can make the simplification , @xmath201 where @xmath202 the coulomb enhancement factor @xmath203 is plotted in fig .  [ fig_f ] with the approximation that @xmath192 ( see appendix [ appx_intstateenhancement ] ) .
the ` 1 + 2 ' spin current is dominated by interference of ` allowed ' one - photon transitions and ` allowed - forbidden ' two - photon transitions , whereas ` 1 + 2 ' spin control is dominated by interference of ` allowed ' one - photon transitions and ` allowed - allowed ' two - photon transitions @xcite . under the approximations that led to ( [ eta_c_simple ] ) and ( [ population control simple ] ) ,
the spin current injection pseudotensor is @xmath204 where @xmath205 is given by ( [ enhancement ] ) , @xmath206 is given by ( [ phase shift ] ) , and @xmath207 is the spin current injection pseudotensor in the independent particle approximation .
under similar approximations , the spin control pseudotensor is @xmath208 where @xmath191 is given by ( [ enhancement_aa ] ) , and @xmath209 is the spin control pseudotensor in the independent particle approximation .
spin control , like carrier population control , has a coulomb enhancement but no phase shift .
there can also be spin control of bound excitons , but it has not been included in ( [ zeta ] ) .
we now examine the relationship between the coulomb enhancements of the ` 1 + 2 ' processes and of one- and two - photon absorption ; the latter are denoted by @xmath210 and @xmath211 so that for @xmath212 , @xmath213 .
the relationship is particularly simple at photon energies for which transitions from the heavy- and light - hole bands are dominant and intermediate state coulomb enhancement is the same for each significant term in the sum over intermediate states .
the coulomb enhancements for the ` 1 + 2 ' processes are then given by ( [ enhancement ] ) and ( [ enhancement_aa ] ) .
for one - photon absorption , @xmath214 @xcite . in noncentrosymmetric semiconductors ,
two - photon absorption is dominated by allowed - allowed transitions just above the band gap , and by allowed - forbidden transitions at higher final energies ; the cross - over point in gaas is a few mev above the band gap @xcite . at photon energies for which allowed - allowed transitions dominate two - photon absorption , from ( [ two photon a - a final ] ) , @xmath215 and thus @xmath216 where @xmath217 , while at photon energies for which allowed - forbidden transitions dominate two - photon absorption , from ( [ two photon a - f ] ) , @xmath218 and thus @xmath219 note that , based on appendix [ appx_intstateenhancement ] , @xmath220 , which is the ratio of the two curves in fig .
[ fig_f ] . in centrosymmetric semiconductors ,
there are no allowed - allowed transitions , and only ( [ f_geometric_mean_af ] ) applies .
the ` 1 + 2 ' processes are often described by ratios .
for example , a useful quantity to describe the current is the swarm velocity @xcite , defined as the average velocity per injected electron - hole pair @xmath221 the swarm velocity is a maximum when the relative intensities of the two colors are chosen such that @xmath222 ; returning to ( [ currentmicro ] ) , if one associates the one- and two - photon amplitudes with the arms of an effective interferometer , this condition corresponds to balancing that interferometer . for fields co - linearly polarized along @xmath223 ,
the maximum swarm speed is @xmath224 a useful quantity to describe pure spin currents is the maximum spin separation distance @xcite ; it is proportional to @xmath225 . as a consequence of ( [ f_geometric_mean_af ] ) , the maximum swarm speed , and the maximum spin separation distance , are _ unaffected _ by excitonic effects when allowed - forbidden transitions dominate two - photon absorption @xcite .
however , close to the band edge , where allowed - allowed transitions dominate two - photon absorption , excitonic effects increase these ratios by a factor @xmath226 over their value in the independent particle approximation . in contrast , as a consequence of ( [ f_geometric_mean_af ] ) , excitonic effects do not affect the maximum control ratio for population and spin control ( @xmath227 and @xmath228 respectively @xcite ) close to the band edge and decrease them by a factor @xmath226 at higher photon energies for which allowed - forbidden transitions dominate two - photon absorption .    in the terminology of seideman @xcite , the excitonic phase shift of the ` 1 + 2 ' current and spin current is a direct phase shift .
it is due to the complex nature of the final state as it appears in the transition amplitudes .
thus it can be understood in terms of the partial wave phase shifts of the final state caused by the coulomb potential between electron and hole .
the coulomb interaction is rather unique due to its long range nature , so we first suppose the potential between the electron and hole falls off more rapidly than @xmath229 . in that simpler problem ,
the final state wave function is written as @xmath230 where the @xmath231 are real @xcite . if the potential between the particles is ignored , then the partial wave phase shifts , @xmath232 are zero .
the allowed one - photon pathway reaches an @xmath127 wave , while the allowed - forbidden two - photon pathway reaches a @xmath145 wave . substituting this form for the wave function into the one- and two - photon transition amplitudes , yields @xmath233 for the one - photon rate , where @xmath234 is real and depends on @xmath235 , and @xmath236 for the two - photon rate , where @xmath237 is real and depends on @xmath235 and @xmath238 . here
@xmath239 is the _
i_-photon transition amplitude when the potential between the particles is ignored .
it is then straightforward to see from ( [ etaschematic ] ) that the relative shift of the partial waves is responsible for the phase shift of the current and spin current .
that is , @xmath240 the use of ionization states as opposed to scattering states was important to get the correct sign of the intrinsic phase . with scattering states
, one would find @xmath241 .
due to the long - range nature of the coulomb potential , the partial wave phase shifts have a logarithmic @xmath88 dependent part , but it is the same for all partial waves and thus does not appear in the relative phase . the part of the coulomb partial wave phase shift @xmath242 that does not depend on @xmath88 is @xmath243 @xcite ; when inserted into ( [ related to partial waves ] ) , this reproduces ( [ phase shift ] ) .
in contrast , the allowed - allowed two - photon pathway reaches an @xmath127 wave and thus there is no phase shift for population control or spin control .    the expression ( [ related to partial waves ] ) for the intrinsic phase in terms of the scattering phases is particularly simple , since each pathway connects to only a single parity .
this contrasts with ` 1 + 2 ' ionization from an atomic @xmath127 state , for which the one - photon transition is to a @xmath145 wave and the two - photon transition is to both @xmath127 and @xmath244 waves ; the intrinsic phase is thus a weighting of the @xmath145-@xmath127 and @xmath145-@xmath244 partial wave shifts @xcite .
materials for which the first term in ( [ bloch_v_expansion ] ) is forbidden ( cu@xmath245o is an example ) have these same selection rules @xcite ; hence , they will have an intrinsic phase with a similar weighting .
the absence of a phase shift in population control can be connected to a symmetry of the second order nonlinear optical susceptibility . from considerations of energy transfer and macroscopic electrodynamics
, @xmath49 is related to the nonlinear susceptibility @xmath246 by @xmath247.\ ] ] in the independent particle approximation @xcite , @xmath248      ^{*},\ ] ] which is a generalization of overall permutation symmetry to resonant absorption . as a result of
( [ gopr ] ) , fraser _ et al_.  showed that @xmath249 is proportional to @xmath250 , and is thus purely real @xcite .
our result that @xmath49 remains real when excitonic effects are included suggests that ( [ gopr ] ) holds more generally .
in fact , it can be shown that ( [ gopr ] ) holds for any hamiltonian symmetric under time - reversal so long as @xmath251 is not resonant .
we have extended the theory of interband ` 1 + 2 ' processes in bulk semiconductors to include the electron - hole interaction . following previous theories
@xcite , we have used a framework based on ( i ) a separation of the initial carrier photoinjection and the subsequent carrier scattering , and ( ii ) a perturbative expansion in the optical field amplitudes , with injection rates obtained in a fermi s golden rule limit for the bichromatic field .
the injection rates for carrier population control , spin control , current injection , and spin current injection , have been described phenomenologically by tensors @xmath49 , @xmath51 , @xmath16 , and @xmath50 , respectively @xcite . like previous theories ,
we have used the long - wavelength limit , and neglect nonlocal corrections to the interaction hamiltonian .
but whereas previous theories of ` 1 + 2 ' photoinjection used the independent particle approximation , we have included excitonic effects .
we have shown that excitonic effects cause ( i ) an enhancement of each ` 1 + 2 ' process , and ( ii ) a phase shift for current injection and spin current injection .
our main results , the modifications of the aforementioned tensors relative to the independent particle approximation are given in ( [ eta_c_simple ] ) , ( [ population control simple ] ) , ( [ mu ] ) , and ( [ zeta ] ) .
these particularly simple results are valid at photon energies for which transitions from the heavy- and light - hole bands are dominant ; more general results are given for @xmath16 and @xmath49 in ( [ eta_result ] ) and ( [ population control result ] ) .
our results are based on the effective mass model of wannier excitons ; degenerate bands are included , but we use a spherical approximation to the exciton hamiltonian , and we neglect envelope - hole coupling .
this is a good approximation for many typical semiconductors , including gaas , since the electron - hole envelope function extends over many unit cells due to the screening of the coulomb interaction by the static dielectric constant @xcite . as a consequence of making the spherical approximation , the phase shifts and coulomb enhancements we find in this paper
are independent of crystal orientation .    also
, our results are limited to low excess photon energy since ( i ) the wannier exciton hamiltonian assumes parabolic bloch bands , and ( ii ) we have truncated the expansion in @xmath0 of the bloch state velocity matrix elements , which is the basis of the transition amplitude expansion . by comparing the black dashed line and grey dash - dotted line in fig .
[ fig_eta ] , one sees that higher order terms in @xmath0 ( for both bands and velocities ) are important in gaas for excess photon energies greater than about 200 mev .
this can then be considered the limit of validity of our calculation . however , combining the coulomb enhancement calculated assuming parabolic bands with the nonperturbative independent particle approximation result ( as was done for the solid black line in fig.[fig_eta ] ) likely gives a good approximation for a few hundred more mev ; this was the case for one- and two - photon absorption @xcite .
it is interesting to ask if there are other sources of intrinsic phases to the current ( or spin current ) besides the one that we have identified here , as these may produce spectral features in the intrinsic phase .
one possibility is the coupling between bound @xmath96-@xmath66 excitons and the unbound @xmath94-@xmath66 or @xmath95-@xmath66 excitons , since it is known that the intrinsic phase can show spectral features near a resonance @xcite .
another possibility is the envelope - hole coupling between the continua of unbound @xmath94-@xmath66 and @xmath95-@xmath66 excitons that was neglected in our treatment .
the effect of coupled continua in the general theory of ` @xmath3 ' phase shifts has not been considered to date .
finally , we note that the intrinsic phase and coulomb enhancement may be greater in reduced dimensional systems , which have greater exciton binding energies .
the carrier - carrier coulomb interaction was included in the theory for ` @xmath7 ' control of electrons in biased asymmetric quantum wells , although the intrinsic phase was not studied @xcite .
consider the functions @xmath252 and @xmath253 , which appear in ( [ n_af ] ) and ( [ n_aa ] ) ; we refer to them collectively as @xmath254 .
first , note that due to the energy conserving delta function in ( [ etaschematic ] ) , @xmath255 will be equal to @xmath256 [ see eq .
( [ energy conservation ] ) ] , and thus @xmath254 is a function only of @xmath257 , @xmath53 , @xmath258 , @xmath259 , and @xmath260 .
second , note that @xmath254 is defined so that if the electron - hole attraction is turned off , for example by letting @xmath261 , then @xmath262 @xcite .
this allows @xmath254 to be identified as part of the coulomb enhancement .
in particular , @xmath254 is the enhancement due to the coulomb interaction in the intermediate states ; if the coulomb interaction is neglected for the intermediate states , @xmath263 @xcite . [
note that lee and fan @xcite did not allow for @xmath264 in @xmath254 ( related to @xmath265 in their notation ) . ]    since the integrand is smooth for the parameter range of interest , numerical integration of @xmath254 is straightforward ; however , it need not be undertaken .
further simplification is possible since the parameter @xmath266 can be considered to be much less than one .
since most materials have an exciton binding energy that is much smaller than the band gap , @xmath251 is detuned from the band edge by many exciton binding energies at photon energies consistent with the approximations made here . in gaas , for example , when @xmath267 is within 500mev of the gap , @xmath266 is at most @xmath268 .
an expansion of @xmath252 for small @xmath266 , @xmath269 where @xmath270 , shows that @xmath252 is approximately @xmath271 and nearly constant as a function of @xmath257 .
the same is true of @xmath253 , which has the expansion @xmath272 where , with @xmath273 , @xmath274 in fact , when @xmath275 , @xmath276 even to fourth order in @xmath277 .
fig.[fig_n ] shows a numerical integration of @xmath252 using the parameters of gaas .
[ intermediate state coulomb enhancement , see ( [ two photon a - f ] ) ] for gaas .
the first ( second ) subscript is @xmath65 ( @xmath278 ) ; subscripts for @xmath279 are not shown ; @xmath280 denotes the heavy and light hole bands , and @xmath281 denotes the split off band.,width=307 ]
the tensor @xmath282 , defined in ( [ eta_ccvv ] ) , can be used to calculate the current injection tensor with or without excitonic effects using ( [ eta_result ] ) or ( [ eta_result_free ] ) .
it can be evaluated analytically in the approximation of parabolic bands .
part of the result for the 8 band kane model has been given before but without the split - off band as an initial or intermediate state @xcite .
we here give more detail , but only for the two - band terms .
we denote the bands by a double index @xmath4 and @xmath127 , where @xmath4 is one of @xmath283 , and @xmath127 runs over the two spin states for each band .
since the coulomb corrections to @xmath16 in ( [ eta_result ] ) do not depend on the spin index , we can include the sum over spin indices from ( [ eta_result ] ) in @xmath284 . and
, since we are only calculating two - band terms , we set @xmath285 and @xmath286 .
thus @xmath287 where @xmath288 for @xmath289 or @xmath94 , @xmath290 , while for @xmath291 , @xmath292 , where @xmath176 is the fundamental band gap , and @xmath293 is the spin - orbit splitting .
the interband velocity matrix elements are approximated by their value at @xmath294 , but still depend on the direction of @xmath0 . in terms of the orthogonal triple of unit vectors
@xmath295 , @xmath296 and @xmath297 , @xmath298 these matrix elements are : @xmath299 _ { s , s^{\prime } } \\
\mathbf{v}_{c , s;lh , s^{\prime } } & = & -\frac{1}{2}\sqrt{\frac{e_{p}}{3 m } } \left [ 2\mathbf{\hat{k}}\mathbb{i}+i\mathbf{\hat{l}}\sigma _ { y}-i\mathbf{\hat{m } } \sigma _ { x}\right ] _ { s , s^{\prime } }
\\ \mathbf{v}_{c , s;so , s^{\prime } } & = & \sqrt{\frac{e_{p}}{6m}}\left [ \mathbf{\hat{k}}\mathbb{i}-i\mathbf{\hat{l } } \sigma_{y}+i\mathbf{\hat{m}}\sigma_{x}\right]_{s , s^{\prime } } , \end{aligned}\ ] ] where @xmath300 is the kane energy @xcite . here
, @xmath301 is the @xmath302 identity matrix and @xmath303 are the pauli spin matrices .
of course , for parabolic bands , the intraband matrix elements are @xmath304 , where @xmath305 is the effective mass of band @xmath4 .
( in the proper kane model , the effective masses are given in terms of the parameters @xmath176 , @xmath293 , and @xmath300 , but we treat them as additional parameters , which is equivalent to including remote band effects on the effective masses . ) the sums over spin then yield @xmath306 the remaining angular integrals can be done using @xmath307    the result for @xmath308 is @xmath309 , \label{etafree parabolic}\ ] ] where the tensor properties are @xmath310 and the term involving @xmath311 should not be included if @xmath312 .
the free particle current injection tensor has also been investigated for parabolic bands , but with a simple three - band model @xcite .
that model does not have the matrix elements @xmath313 and @xmath314 , and thus differs from our result for @xmath315 .
this work was financially supported by the natural science and engineering research council , photonics research ontario , and the us defense advanced research projects agency .
we gratefully acknowledge many stimulating discussions with daniel ct , james fraser , ali najmaie , fred nastos , eugene sherman , art smirl , marty stevens , and henry van driel .


low energy nucleon - antinucleon annihilation is a fertile area for studying hadron dynamics and particularly pionization .
experiments give information on pion numbers , spectra , and correlations and seem to show little dependence on initial energy from annihilation at rest to kinetic energies of a few hundred mev .
results do depend on the iso - spin and spin quantum numbers of the annihilating pair .
understanding these results is a serious challenge to theory .    in this paper
we develop a description of nucleon - antinucleon annihilation into pions based on the creation of a coherent classical pion wave in the annihilation process and the quantization of this wave using the method of coherent states .
about 40% of annihilations go through meson resonances , mostly rho and omega mesons , and then to pions @xcite .
these channels do not lead to prompt pions and are not part of our treatment .
our first studies show our picture accounts for the major phenomenological features of the pions seen in low energy annihilation .
our picture divides the process into two steps .
the first is the dynamics of the nucleon and antinucleon .
we have little to say about this except to assume that once these two  touch " and annihilation begins it proceeds very rapidly .
this rapid annihilation leads to a pion pulse or pion wave that forms the basis of our coherent state .
it is the description and consequences of that state and the way in which its features reflect the spatial and temporal evolution of the annihilation region as well as its quantum numbers that is our principal concern @xcite .
our description of annihilation as occurring by the rapid radiation of pions , draws its inspiration from recent studies of annihilation in the skyrme model .
et al _ @xcite have studied numerically the collision of two skyrmions , one of baryon number one and the other of baryon number minus one .
they find that once there is significant overlap of the two , they disappear into pion radiation at nearly the causal limit .
et al _ @xcite studied this rapid decay in a schematic model .
they considered an initial well localized  blob " of skyrmionic matter with baryon number zero and total energy about twice the single skyrmion energy , and studied its time evolution .
they found that it decays into a pion wave nearly at the causal limit .
that pion wave is a propagating , but well localized pulse .
this picture of the annihilation process in terms of a rapid coherent radiation pulse is as far from the thermodynamic one as a broadcasting radio antenna is from a black body .
the skyrmionic description is classical , and is meant as purely suggestive .
does a classical starting point make sense for annihilation ?
let us assume that the rapid conversion to pion radiation seen in the skyrme treatments is a feature of the real annihilation .
then the pions are radiated from this  bright flash " as a coherent pulse , or coherent pion quantum wave .
since the mean number of pions radiated is not small , and the total energy to be radiated is large compared with the pion mass , @xmath0 , we can envisage this wave to be , in first approximation , a * classical coherent wave*. the remaining problem is then to quantize this wave .
this we do with the well known method of coherent states . our picture ,
then , has two central assumptions .
first we assume that rapid annihilation leads to a classical pion wave pulse .
this is a consequence of the rapid time scale of the annihilation and the large energy ( compared to the pion mass ) released .
second we assume that this classical wave can be quantized using coherent states .
this is basically the first term in a description with expansion parameter @xmath1 , where @xmath2 is the mean number of pions .
the error in our second assumption can therefore be estimated .
it should be noted that our starting point arises in classical qcd , which is the non - perturbative domain of the theory .
this seems to us the correct starting point for the study of annihilation .
it may be that other strongly interacting systems are also best approached from this direction .
it has recently been shown that that is true for the nuclear force @xcite . in such treatments
it is the reintroduction of quantum aspects that is the problem . for radiation fields ,
coherent states are constructed to solve that problem .    in the next section
, we develop the coherent state formalism for the pions .
we begin with a description of classical pion waves from a source , introduce coherent states , discuss the constraints of iso - spin and of energy and momentum conservation and give some statistical features of the coherent state . in section 3
we briefly investigate the consequences of our description for the annihilation data .
we only touch on what could be done , but we do see that our approach generally agrees with experiment .
finally section 4 discusses some directions for further work and presents a brief summary .
we imagine nucleon antinucleon annihilation as occurring very rapidly from a small region of high energy density . in first approximation
, we take this region to radiate a classical coherent pion wave .
if we further neglect @xmath3 interactions and the back reaction of the radiated pions on the source , we can write the classical wave equation for the pion wave @xmath4 , as @xmath5 here @xmath6 is the source of the pion field , @xmath7 .
this source will be assumed to be sharply localized in space and time . for the moment we neglect the isospin degree of freedom and will return to it later .
we also leave out the factor of @xmath8 in front of the source that is familiar from electrodynamics , since it adds nothing here .
eqn ( 1 ) gives the form of the pion wave as a classical wave radiated by the source @xmath9 .
the details of the annihilation mechanism are contained in this source function .    the equation for @xmath7 is easily solved in terms of @xmath9 by taking fourier transforms .
we define the fourier transform of @xmath7 by @xmath10 and a corresponding function , @xmath11 as the fourier transform of @xmath6 . the pion field
is then given by @xmath12 the pion radiation field comes entirely from the poles of the denominator in the integral .
we must specify how those poles are to be treated in order to get causal radiation . doing so and using the pole contributions only to evaluate the @xmath13 integral in ( 3 ) , we obtain the approximate expression for the pion field , for @xmath14 , @xmath15 where @xmath16 .
this form for the pion field is exact in the radiation region , but neglects non - radiating parts of the field coming from the @xmath13 dependence of @xmath17 . in the radiation zone
there is no source and @xmath7 given in ( 4 ) satisfies the source free wave equation .
it is this form for @xmath7 that we will take to represent the radiated pions of nucleon antinucleon annihilation .
it is not the classical pion wave that is detected in annihilation experiments but rather the quanta of that field , the pions .
we must therefore extract the quantum content of this classical wave .
this is exactly what the coherent state method achieves .
pioneered by glauber @xcite for optics in situations where the photon granularity is significant but not overwhelming , the coherent state formalism allows one to construct a quantum state that corresponds , in the large number of field quanta limit , to the classical wave field .
coherent states for pions were studied by horn and silver @xcite , and many of the arguments given here were considered by them .
they did not , however , study annihilation .
because the pion field satisfies the wave equation without sources , it can be decomposed into a linear superposition of free waves with wave number @xmath18 and corresponding energy @xmath19 .
that is the content of eqn ( 4 ) .
we can bring in quantum mechanics for this wave by introducing the standard creation and annihilation operators for the modes @xmath18 , @xmath20 and @xmath21 .
these satisfy the usual harmonic oscillator commutation relations , @xmath22 = \delta({\bf k}-{\bf k'}),\ ] ] with all other commutators zero .
the quantum field operator corresponding to the pion field is then @xmath23 the coherent state associated with a given classical state , @xmath24 , is the quantum state that is an eigenstate of the positive frequency part of @xmath25 , @xmath26 , with the corresponding positive energy part of the classical wave as its eigenvalue . recall that the positive frequency part of the state has the @xmath27 dependence . in the quantum operator for the field , this is the part that goes with the annihilation operator .
since the modes for each @xmath18 are independent , finding the eigenstate is equivalent to finding the eigenstate of one single annihilation operator .
the full state can then be built up by superposing the independent modes .
the normalized eigenstate , @xmath28 , with eigenvalue , @xmath29 , of a single mode annihilation operator , @xmath30 , ( @xmath31 ) , is given by @xmath32 where @xmath33 is the vacuum state .
this equation for the eigenstate of @xmath30 is the key to the entire coherent state formalism .
consider now the normalized quantum state @xmath34 defined by @xmath35 this is just a coherent state in which each mode @xmath18 carries weight @xmath36 .
it is clear , for the positive frequency part of the field , that @xmath37 with the eigenfunction @xmath38 given by @xmath39 this @xmath38 satisfies the source free wave equation for any function @xmath40 .
we then take @xmath34 to be the normalized quantum state corresponding to the classical wave @xmath38 .
note that @xmath34 does not contain a fixed number of quanta , but is rather a coherent superposition of states of different numbers of quanta .
the meaning of a coherent state arising from a classical source can be clarified by relating it to the standard s - matrix description of a field generated by a classical external source .
the classical source , @xmath9 , coupled to the field , @xmath7 , creates from the vacuum the quantum state @xcite @xmath41 upon four dimensional fourier transformation and integration over @xmath42 , this reduces , up to a normalization , to eqn ( 8) . the exponent in ( 11 ) is just the s - matrix in the interaction representation , and makes clear both the connection of the s - matrix and coherent state approaches and why the resulting state has an indefinite number of quanta .    in the electromagnetic case it seems natural to define states with indefinite numbers of photons , but even there
if the wave has a fixed total energy only a fixed number of photons above some frequency can be detected . with pions
the description in terms of a coherent mixture of states with different numbers of pions seems less comfortable , particularly in view of the finite pion mass .
no one would suggest representing a pion wave with total energy near the pion rest mass as a classical object . in annihilation , however , the total energy released is about 13 times the pion rest mass , and the approximation of representing it by a classical coherent wave seems at least plausible .
the purpose of this paper is to show that that description is not only plausible but useful and gives dynamical insight into the annihilation process and may have application in a number of other processes in which many pions are created .    to make the pion coherent state correspond to the radiating solution @xmath24 of eqn ( 4 ) , that is to make @xmath38 in ( 9 ) or ( 10 ) correspond to the positive frequency part of @xmath7
, we must take in ( 8) @xmath43 with this choice , our quantum coherent state represents the classical radiating pion wave generated by the source @xmath6 .
the pion is an iso - vector and comes in three charge states . in principle
each charge state can have its independent source , corresponding classical wave , and coherent state .
this is most easily formulated by making @xmath24 and @xmath6 of eqn ( 1 ) each iso - vectors .
then each charge state could have a source density that varies arbitrarily in space and time .
in fact nucleon antinucleon annihilation occurs in a state of fixed total iso - spin , 0 or 1 .
the full pion wave must have the same total iso - spin .
this condition of total iso - spin conservation imposes a * global * condition on the source and on the coherent state .
since the different @xmath18 modes of the wave are independent , they must each combine to the correct total iso - spin , which means that the iso - spin source density can not depend on @xmath44 or @xmath45 . in the coherent state @xmath34 of ( 8) we expect in general that @xmath46 is an iso - vector dotted into the now iso - vector creation operator .
but if we are to construct states of fixed iso - spin for any @xmath18 configuration , the iso - vector dependence of @xmath40 must be independent of @xmath18 . in that case @xmath40
is a constant vector in iso - spin space .
the direction of that iso - vector depends on the source .    for
the moment let us consider a coherent state specified by a fixed unit vector in iso - spin space @xmath47 and an iso - scalar function @xmath46 .
call that state @xmath48 . in place of (
8) we write @xmath49 that is a coherent state constructed from a particular superposition of pions , that superposition made up of pions all pointing in the iso - direction of @xmath47 .
this is not yet a state of fixed iso - spin .
in fact it is easy to see that it contains states of all iso - spin .
a state of fixed iso - spin must be projected from it .
call the projected state of total iso - spin @xmath50 with z - component @xmath51 , @xmath52 . the projection is given by @xcite @xmath53 where @xmath54 is the usual spherical harmonic , and @xmath55 is a normalization constant . the state defined in ( 14 ) is orthogonal to states of different total iso - spin or z - component projection . to see this and to calculate the normalization let us study the inner product of two projected states .
this is given by @xmath56 it is easy to show that @xmath57 for a strong field , large @xmath40 , this inner product is very sharply peaked around @xmath58 .
we will show in the next section that @xmath59 where @xmath2 is the mean number of pions emitted in the annihilation .
we will use this notation here both to save writing and to remind us that this integral , or @xmath2 , is large . to show the orthogonality in ( 16 ) write @xmath60 from this it is easy to show that @xmath61 from which the normalization can be read off . in particular for large @xmath2 , we can use the asymptotic form of the bessel function to yield for the normalization @xmath62    we should note that our iso - spin projection has the interesting property that states of even iso - spin contain only even numbers of pions and states of odd iso - spin only odd numbers . to escape this restriction , as the data does , we need to include the channels for annihilation into meson resonances .
we would now like to examine the properties of the coherent pion state generated by the assumed initial source @xmath6 .
we begin with a discussion of the case without iso - spin .
this coherent state is given by ( 8) with @xmath46 given in terms of the fourier transform of @xmath9 as in ( 12 ) .
although the state does not contain a fixed number of pions , we can ask what is the average number , @xmath2 . this is just the expectation value of the total pion number operator in the coherent state .
we need @xmath63 ( note that we are working in the heisenberg representation so that the states are time independent while the operators in general carry a time dependence .
however , the number operator is time independent . )
this expectation value is easily evaluated using the special properties of the coherent states to give ( 17 ) .
thus we see that the integral that enters to normalize the coherent state in ( 8) , is just @xmath2 .
it is large @xmath2 that corresponds to large field ( @xmath40 or @xmath9 large ) and thus to a good approximate correspondence between the classical and quantum descriptions .
it is also clear from this that the average single pion momentum distribution is given by @xmath64 in terms of the square of the fourier transform of the source density .    in the same way one obtains for the total energy released ,
the expectation of the hamiltonian , or @xmath65 in a quantum process , this energy is sharp and should be put equal to the total energy released in the annihilation .
this fact serves as a convenient way to normalize @xmath40 .
the sharpness of the state can be seen by calculating the dispersion in the number of pions emitted .
one finds @xmath66 which shows that the fractional dispersion in @xmath67 is @xmath68 .
this result also follows from the fact that the probability of finding a state of @xmath69 pions is given by a poisson distribution , as we shall show below .
one can use the coherent state to calculate the amplitude for finding some particular quantum configuration .
for example the amplitude for finding @xmath69 pions of momenta @xmath70 is given by @xmath71 the probability of finding @xmath69 pions of any momentum is then @xmath72 which is the poisson distribution .
it is clear that the coherent state does not
know " not to emit more pions than energy conservation permits .
we will turn to these constraints of energy and momentum conservation below .
our discussion here serves to emphasize that the coherent state approach is better at giving average information over the ensemble of annihilations than it is for particular quantum states .
if we include iso - spin projection as in section 2.3 , we can also probe questions about iso - spin populations . for example
let us study the average number of pions of a particular type in a projected coherent state with total iso - spin @xmath50 and z - component @xmath51 .
this is given by @xmath73 where @xmath0 is the pion type ( @xmath74 ) and there is no sum over @xmath0 implied . using the forms of the projected states ,
this becomes @xmath75 using ( 18 ) and expressing the @xmath47 in terms of spherical harmonics , the integrals can be calculated in terms of standard identities .
we continue to use ( 17 ) to express the integral over @xmath40 in terms of @xmath2 .
we now use ( 17 ) to define @xmath2 but as we shall see it is still the average number of pions summed over all pion types .
we then find @xmath76 in terms of a sum over clebsch - gordan coefficients . in obtaining ( 29 )
we have used the asymptotic forms of the bessel functions ( large @xmath2 ) so that there is no @xmath77 dependence in the normalization .    for nucleon - antinucleon annihilation
there are only two choices for @xmath50 , 0 or 1 . for the case of @xmath78 ( 29 )
reduces to @xmath79 for any @xmath0 . as we expect for a state of iso - spin zero the average number of pions is type independent , and each is one third of @xmath2 , making @xmath2 the true average number of pions summed over types , as promised . for @xmath80
there are three cases . for @xmath81
one finds @xmath82 we see that again the sum of averages over types is @xmath2 .
note the dominance of @xmath83 , an important feature of the data . for @xmath84
one finds @xmath85 with the same result for @xmath86 .
it is no surprise that for the @xmath78 case or for the @xmath80 , @xmath81 case the average number of @xmath87 and @xmath88 is the same .
it is a bit more puzzling for @xmath80 , @xmath89 . what happened to charge conservation ?
the initial @xmath89 state has charge @xmath90 .
however , @xmath91 is a number of order @xmath2 . to that order
all three @xmath51 values do correspond to states of average charge zero .
charge conservation is a @xmath1 effect , and we are neglecting terms of that order .
we will return to this approximation below .
the same charge ratios can be obtained in ( 28 ) by making the overlap of ( 16 ) proportional to @xmath92 .
we have already commented that this overlap is strongly peaked around @xmath58 for large @xmath2 .
this delta function approximation is correct in the large @xmath2 limit and it makes calculations of the average pion number by type and of higher moments of the distribution very simple .
we illustrate this delta function method in the calculation of isospin number correlations .
for example the joint average of @xmath0 type and @xmath55 type pions in a state of total iso - spin @xmath50 with z - component @xmath51 is @xmath93 in the large @xmath2 limit where the overlap may be replaced by a delta function , we obtain , @xmath94 in this expression the components of @xmath47 are , in the usual spherical coordinates , @xmath95 higher order forms are obtained from the generalization of ( 33 ) and ( 34 ) .
they come from higher powers of the number operator for the @xmath96th type in the expectation value with the coherent state replacement @xmath97 the integrals that appear in ( 34 ) are elementary , and hence these higher moments are easily calculated .
greiner , gong and muller @xcite have recently introduced polynomial coherent states in a discussion of pion condensates in heavy ion collisions .
they consider only @xmath78 states and calculate the iso - spin pair correlations defined as @xmath98 we obtain the same values as they do for @xmath99 for large @xmath2 , but our method is simpler .
we can also use the method to calculate the pair correlations in other states .
for example we find for @xmath80 , @xmath81 , @xmath100 and for @xmath80 , @xmath84 , @xmath101 .
other moments and other amplitudes are also easily calculated and agree with those of greiner , gong and muller where calculated by them .
thus we see that the coherent state approach makes definite predictions about average charge ratios in annihilations and can also be used to calculate higher moments of those distributions .
it also predicts that the single pion momentum distribution should be independent of pion type .
so far the coherent state we have defined does not have definite total energy or momentum .
for example it contains contributions from states with only one pion or with arbitrary numbers of pions , neither of which is allowed by conservation of total energy and momentum .
this fault is not serious for large @xmath2 , but is significant for the intermediate @xmath2 of annihilation .
energy and momentum conservation can be imposed on the coherent state , as was shown by horn and silver @xcite , at the expense of some further complication of the formalism .
let us begin by imposing energy and momentum conservation on the coherent state without iso - spin .
the method for imposing energy and momentum on the coherent state is formally the same as the projection method for imposing fixed iso - spin .
the projections commute , but taking them together complicates the physics . for iso - spin one first constructs a state with all the pions pointing in the same direction in iso - spin and then averages over that direction with the appropriate @xmath102 eigenfunction weight . to construct a state of fixed energy and momentum ( or four - momentum )
one first constructs a state of all the pions at a fixed place in space and time and averages over those places and times with the appropriate eigenfunctions of definite energy and momentum . to do this define the operator @xmath103 by @xmath104 then the ( un - normalized ) coherent state at @xmath44 and @xmath45 is @xmath105 it is more concise as well as more physical to use four dimensional notation .
we call @xmath106 the four - vector of position with space part @xmath44 and time part @xmath45 . thus @xmath107 becomes @xmath108 .
then the state of fixed total four - momentum , @xmath109 is given by @xmath110 to see how this works and to calculate the normalization , consider the overlap of two states of different four - momentum .
one easily finds @xmath111 where @xmath112 and where the fourth component of @xmath113 under the integral is @xmath114 .
a simple change of variable reduces the overlap to @xmath115 this clearly shows that states of different four - momentum are orthogonal , but the remaining integral is singular due to the undamped large @xmath106 behavior of the integrand . horn and silver
@xcite also pointed out this problem .
it arises from the first term in the expansion of the exponential .
this term , a  one " , corresponds to the no pion contribution to the annihilation process .
since there is in fact no such contribution that conserves energy and momentum , we can subtract it out .
similarly there is no one pion state that can contribute to annihilation and conserve four - momentum , thus we may also subtract out the one meson state and define a coherent state of total four - momentum @xmath109 that starts with two pions .
this is given in terms of @xmath116 by @xmath117 we then find @xmath118 the integral now is convergent and we have lost nothing by the subtractions since we have only removed states that can not contribute physically .
however the integral in ( 46 ) is still quite delicate and difficult to evaluate numerically because it contains delta functions coming from the @xmath119 integral .
we can get around this problem and shed more physical light on the expression by expanding the exponential and interchanging the order of integration .
let us call the integral in ( 46 ) @xmath120 .
we can write @xmath121 using the definition of @xmath96 and interchanging integrations one finds @xmath122 to be @xmath123
. integrals of this form are easily done by monte carlo methods following a program given by barger and phillips .
the individual terms in the @xmath69 sum in ( 47 ) represent the relative contributions to annihilation into @xmath69 pions . for fixed total energy
the sum must terminate .
for example , nucleon - antinucleon annihilation at rest can not go into more than 13 pions , and in that case @xmath124 must be zero for @xmath125 .
in fact @xmath126 is just the fraction of annihilations that go into @xmath69 pions .
thus the mean pion number is given by @xmath127 this is equal to the expression that one arrives at by taking the expectation value of the number operator in the state @xmath128 and normalizing .
this expression can also be written @xmath129 note that @xmath130 in terms of the old expression for @xmath2 , ( 17 ) .
for @xmath2 large , @xmath131 is large , and the integrands in ( 50 ) are dominated by the factors of @xmath132 and @xmath133 . the expression ( 50 )
can then be thought of in terms of a ratio of improper integrals @xmath134 these integrals can be approximated by the method of stationary phase . for annihilation at rest , the total energy released , @xmath135 , is large , and then the stationary phase point comes at @xmath136 .
thus the @xmath131 in the numerator of ( 51 ) , but not in the exponent , can be evaluated at @xmath137 to give that the new @xmath2 is approximately equal to the old one .
corrections are of order @xmath1 .
we have verified this in a number of numerical examples .
we see that four - momentum conservation can indeed be implemented for the coherent state , but for average quantities , and for @xmath2 large , that constraint is not important . as we will see below , it is important for the probability distribution and for higher moments of the distribution .
it is clear that one can combine the method outlined above for projecting onto states of good four - momentum with the method of the previous section for projecting onto states of good iso - spin .
we now turn to a discussion of how to do that .      combining the methods of previous sections
, we can construct a state of fixed total four - momentum @xmath109 and fixed iso - spin @xmath50 with z - component @xmath51 .
we first generalize @xmath138 of ( 38 ) to give @xmath139 then the appropriate generalization of ( 44 ) to include iso - spin is @xmath140 giving the state ( not normalized ) @xmath141 the inner product of two such states belonging to different four - momentum and iso - spin gives a four dimensional delta function on four - momentum and kronecker deltas on iso - spin and its z - component times the following normalization integral @xmath142 this integral can be calculated by the same expansion methods used in ( 46 ) and ( 47 ) . in terms of the @xmath143 defined in ( 47 )
the normalization integral of ( [ eq : normal ] ) is given by @xmath144 where @xmath145 this integral can be done using the identity in ( 18 ) by differentiating ( 18 ) with respect to @xmath2 @xmath69 times and then setting @xmath2 to zero .
one finds @xmath146 it should be recalled that sums over @xmath69 as in ( [ eq : fmsum ] ) have finite upper limits because of the constraint of energy conservation .
the expected number of pions of type @xmath0 corresponding to ( 27 ) is now @xmath147 the joint average of @xmath0 type and @xmath55 type pions as in ( 33 ) is given by @xmath148 in obtaining ( [ eq : corr ] ) , we have neglected a term from the commutator that is down by @xmath1 . in these expressions
, @xmath96 is still given as in ( 42 ) .
these expectation values can be evaluated , as before , by expanding .
once again they are expressed in terms of the @xmath124 .
we have @xmath149 and @xmath150 where @xmath151 and @xmath152    let us examine these forms in the cases of physical interest . for @xmath78 , it is easy to show that @xmath153 then from ( 56 ) and ( 61 ) we get @xmath154 as we expect . for @xmath80 , @xmath81 , we find @xmath155 these do not seem to give the same charge ratios we obtained before in this state , but those were calculated in the large @xmath2 limit . in that limit the dominant terms in the sums of ( 56 ) and ( 61 ) come from terms with @xmath156 .
thus the dominant terms in @xmath157 that matter have @xmath158 ( recall that @xmath77 is of order @xmath50 which is in turn of order @xmath159 ) . in this limit
it is easy to show that @xmath160 with corrections of order @xmath161 . in that limit ( 67 )
yields the results of ( 31 ) .
similarly for @xmath80 , @xmath84 , we find @xmath162 which yields ( 32 ) in the large @xmath2 limit , and shows how charge conservation enters for finite @xmath2 .
one can calculate the two pion averages in the same way .
the factors that enter become for @xmath80 , @xmath81 @xmath163 and for @xmath80 , @xmath84 @xmath164 which yield our previous results for the iso - spin pair correlations @xmath99 defined in sect .
2.5 in the large @xmath2 limit .
the evaluation of the averages now depend on the functional form of @xmath46 .
we turn to it in the next section .
as a first orientation into the phenomenological content of the coherent pion wave description of annihilation , let us study the relation between the average pion number @xmath2 of ( 17 ) and the energy released , ( 23 ) , for a very simple model of @xmath40 .
we do not take @xmath40 directly from the experimental pion spectrum since the pions from the meson resonance channels confuse the extraction .
rather suppose we assume that the pion field source turns on at @xmath165 and then decays exponentially in time , and that it has a spherically symmetric yukawa shape .
( spherical symmetry is appropriate for annihilation at rest . )
that is we take for @xmath6 of ( 1 ) , @xmath166 where @xmath167 is the source strength and where we have put in a factor of @xmath45 to make the time dependence continuous at @xmath165 . using ( 12 )
this leads to @xmath168 where @xmath169 is a strength and where we have multiplied @xmath40 by @xmath170 to model the p - wave nature of pion emission . with this form ,
the integrals for the mean number of pions and for the mean energy ( 17 ) and ( 23 ) can not both be done analytically , but they can be easily evaluated numerically . this requires a choice of the parameters @xmath171 , @xmath172 and @xmath169 .
the last can be fixed by requiring that the average energy be the energy released in annihilation . in units of the pion mass ( @xmath173 ) , and for annihilation at rest , this energy is 13.87 . for the range parameters , we take @xmath174 .
this corresponds to an annihilation region with a time and distance scale of half a pion compton wave length - a reasonable size .
it is not the purpose of this paper to make a careful study of the source function and to fit it to data .
rather we simply want to demonstrate that a quite reasonable choice for functional form and size yields a correspondingly reasonable account of the data .
thus we have not made an exhaustive parameter search . with our choice for parameters
we find that ( 17 ) gives for the average pion number @xmath175 . with the same choice of parameters , the average pion number calculated with four - momentum conservation imposed as in ( 50 ) gives @xmath176 .
this is certainly within @xmath1 of the unconstrained value , as promised .
also @xmath2 about 6 agrees with the data , for reviews see , @xcite , particularly when it is realized that we are only modeling that part of annihilation that goes into uncorrelated pions , and not the part that goes into other mesons that subsequently decay into pions . since these typically give few pions , the average number without them will be somewhat higher than the total average number .
although the average number from the poisson distribution and the distribution constrained by four - momentum conservation come out very close , the actual probability distributions in pion number are rather different . in fig .
1 we show the probability of finding @xmath69 pions as a function of @xmath69 for the poisson and for the  m - sum " cases ( as defined in ( 47 ) and ( 48 ) ) .
we see that the constrained case has a far narrower distribution than the poisson . for the poisson distribution
the @xmath177 ( @xmath178 is defined in ( 24 ) ) is @xmath179 . for the distribution constrained by four - momentum conservation
we find @xmath180 , which is quite close to the value deduced @xcite from experiment , @xmath181 .
this value has a statistical error of at least five percent and an additional unknown error from modeling the distribution of neutrals .
in fact the entire pion multiplicity distribution from the constrained calculation shown in fig . 1 is quite close to the experimental distribution @xcite , @xcite .
our constrained distribution is also indistinguishable ( on the scale of the figure ) from a gaussian distribution with the same average number and same variance .
that is what a statistical model would give .
one can also use the combination of iso - spin projection and four momentum projection of sect .
2.7 to calculate the pion number distributions .
they are shown in fig .
2a and 2b for @xmath78 and @xmath80 , @xmath81 respectively .
recall that for @xmath78 only an even number of pions is possible while for @xmath80 there is only an odd number .
the data is an appropriate average over these two , and again strongly resembles our figures @xcite , @xcite .
annihilation into vector mesons must be added to evade this odd - even effect .
let us now turn to a comparison of our picture with experiment for the average number of pions by iso - spin type .
we compare with the ratios reported in sedlak and simak .
to do this we need to estimate the relative population of @xmath78 and @xmath80 in proton antiproton annihilation at rest .
we do so following the work of locher and zou @xcite . then using the charge ratios in states of good iso - spin we found in sect .
2.5 , and for vanishing total third component of iso - spin , we find @xmath182 with an error of about ten percent to be compared with the experimental number reported as @xmath183 @xcite , .
this is remarkable agreement .
note that these ratios do not depend at all on the details of our parameters , but do depend on the structure of the coherent state picture .
in particular , it should be noted that the excess of @xmath83 s found experimentally comes out naturally in our treatment .
given @xmath184 , we can calculate charge averages including four - momentum conservation as in sect .
using the same parameters , we can find the average pion number by charge type and the variance in the various iso - spin states .
these are shown in table 1 .
the @xmath80 , @xmath84 case shows the effect of charge conservation , for @xmath84 , the average number of @xmath87 is one larger than the average number of @xmath88 .
we also see that with four - momentum conservation the ratio of @xmath83 to @xmath87 in the @xmath80 , @xmath81 case is somewhat larger than in the unrestricted case .    to make our formalism a real theory of the annihilation process we need to model annihilation in space and time as a classical source of pions .
we have presented a very simple example as an orientation .
more sophisticated models of the pion source could be made based on a dynamical theory of the annihilation process , for example the skyrme model calculation @xcite , or on nucleon - antinucleon potential model calculations .
these could include the spatial asymmetry seen in these annihilation calculations @xcite .
this asymmetry would lead to a different momentum distribution for the pions .
alternately one could use ( 22 ) to extract the form of @xmath46 from the single pion momentum distribution data .
we must also model the substantial part of annihilation that goes through meson resonances .
we plan to return to these ideas in later work .    given a form for @xmath46
, one can calculate not just single pion averages , but also correlations @xcite , @xcite and joint momentum distributions for multiple pions .
the great advantage of the pion coherent state approach is that it makes detailed , specific and rather restrictive predictions for these observables and thus is easily verifiable , or perhaps falsifiable .
clearly much work remains to be done to develop and test the coherent state approach developed here for pions radiated from nucleon antinucleon annihilation .
predictions for higher moments of pion distributions from nucleon antinucleon states of sharp iso - spin have to be confronted with data .
correlations have to be studied .
theories of nucleon antinucleon interaction and of the development of the annihilation region have to be studied and interpreted in terms of a source of the pion radiation . it may be possible , for example , to model annihilation classically in the skyrme model including the non - linear interactions of the pions , and the effects of the annihilation back on the source .
such a description would lead , far enough away from the annihilation region , to a free , non - interacting pion wave .
this wave could then be the source of the coherent state .
this approach is far simpler than a full quantum theory of the annihilation . beyond this study of the annihilation process , we must study the effect of other quantum numbers such as g - parity and angular momentum on annihilation in the coherent state description .
finally a better understanding of the relation between annihilation into pion radiation and into other mesons is needed .
we are exploring all these questions
. should these detailed descriptions fail , it would be interesting to go beyond the coherent state picture by including squeezed states .
one can also imagine other applications of the general methods advocated here .
the basic idea is to look for any fast coherent processes in which the total energy radiated is large compared with the mass of the radiated quanta and in which the interaction among those quanta or by the quanta back on the source can be neglected .
then the radiation process is described classically , and the subsequent wave quantized using coherent states . conserved quantum numbers have to be imposed on these states as we have imposed iso - spin .
one can then explore averages , moments , momentum distributions and their relation to the radiating source , all as we did for annihilation .
greiner , gong and muller @xcite have discussed one such possible application in heavy ion collisions , there are many more
. there may be similar opportunities in the description of hadronization in qcd jets @xcite .
other examples will occur to the reader .
it may even be possible to lift the restriction of no self interaction of the field or of the field back on the source so long as one can treat these classically and then use coherent states to quantize only in the radiation zone .
many qcd processes , like annihilation , may best be treated by starting in the classical non - perturbative domain of qcd .
it may also be possible to extend these ideas to heavy ion collisions in general @xcite , where one could take into account the statistical nature of the process by introducing a density matrix and averaging over an appropriate ensemble of coherent states .
such methods are well known in quantum optics @xcite , and may be an attractive alternative to the thermodynamic approach to heavy ion collisions .
in summary we have seen that the large energy released and the possibility of very short reaction time suggest that a classical , coherent pion wave pulse is radiated in nucleon antinucleon annihilation .
this wave can be quantized using the method of coherent states . from this
description averages and higher moments of pion distributions in annihilation can be calculated . in particular by imposing iso - spin and four - momentum conservation , interesting correlations among pion charge types
are easily calculated .
momentum distributions contain information on the spatial and temporal distribution of the pion wave source .
preliminary calculations show that all these features obtained from the coherent state picture agree with the data and therefore suggest that the coherent classical pion pulse picture of annihilation is valid for that part of annihilation that goes directly into pions .
we also speculate on other applications of classical meson waves and their subsequent quantization .
rda , fc , and jpd would like to thank the theory group of the division of nuclear and particle physics at the paul scherrer institute and for providing a very pleasant environment in which much of this work was done .
rda and bs are supported in part by the united states national science foundation .
the division de physique thorique is a research unit of the universities paris 11 and 6 associated to cnrs .
rda and bs would like to thank dr .
r. hollebeek for bringing reference to our attention .
hg c. amsler and f. myhrer , ann .
nucl . part .
sci . * 41 * , 219 ( 1991 ) .
+ a brief account of this work has been published in r.d .
amado , f. cannata , j - p .
dedonder , m.p . locher and b. shao ,
lett . * 72 * , 970 ( 1994 )
sommermann , r. seki , s. larson and s.e .
koonin , phys.rev . *
d45 * , 4303 ( 1992 ) .
+ b.shao , n.r.walet , and r.d.amado , phys.lett .
* b 303*,1 ( 1993 ) .
walet and r.d .
amado , phys.rev .
* c47 * , 498 ( 1993 ) . + r.j .
glauber , phys .
rev . * 131 * , 2767 ( 1963 ) + d. horn and r. silver , annals of physics * 66 * , 509 ( 1971 ) .
ginocchio and m.w .
kirson , phys .
letters * 44 * , 1744 ( 1980 ) + c. itzykson , j.b .
zuber , quantum field theory , ( mcgraw - hill , singapore ( 1985 ) ) , chapter 4 .
+ the iso - spin projection presented here was independently and long previously developed for coherent states by j.c.botke , d.j.scalapino , and r.l.sugar , phys.rev .
* d9 * , 813 ( 1974 ) .
we apologize to these authors for having missed this work in our earlier publication . v.
barger and r. phillips , collider physics , ( addison - wesley , redwood city , california ( 1987 ) ) , chapter 11 .
+ c. amsler and f. myhrer , ann .
nucl . part .
sci . * 41 * , 219 ( 1991 ) + j. sedlak and v. simak , sov .
* 19 * , 191 ( 1988 ) .
dover , t. gutsche , m. maruyama and a. faessler , prog .
* 29 * , 87 ( 1992 ) .
+ c. ghesquiere in : proc.int.symp .
on antinucleon - nucleon interaction , liblice - prague , 1974 , cern 74 - 18 , p. 436
boal , c.k .
gelbke and b.k .
jennings , rev.mod.phys . * 62 * , 553 ( 1990 ) + h.q .
song , b.s .
zou , m.p .
locher , j. riedlberger and p. troul , z. phys .
a - hadrons and nuclei * 342 * , 439 ( 1992 ) + m.p . locher and b.s .
zou , z.phys a * 341 * , 191 ( 1992 ) .
+ c. greiner , c. gong and b. muller , phys.lett .
b * 316 * , 226 ( 1993 ) .
dokshitzer , v.a .
khoze , s.i .
troyan and a.h .
mueller , rev .
* 60 * , 373 ( 1988 ) . + j - p .
blaizot and a. krzywicki , phys . rev .
* d46 * , 246 ( 1992 ) .
+    figure captions + figure 1 . the probability distribution , @xmath185 , for finding @xmath69 pions from nucleon antinucleon annihilation at rest .
the open circles are the poisson distribution and the solid squares are the distribution with the constraint of four - momentum conservation imposed . that distribution can not be distinguished ( on the scale of the figure ) from a gaussian of the same mean and variance .
+ figure 2 .
the probability distribution , @xmath185 , for finding @xmath69 pions from nucleon antinucleon annihilation at rest with iso - spin projection . in all cases
the constraint of four - momentum conservation has been imposed .
figure 2a is for @xmath78 and figure 2b for @xmath80 , @xmath81 .
the open circles are the probabilities with iso - spin projected and the solid squares without iso - spin projection .
the solid squares are the same as the solid squares in figure 1 .
note that in our model , for @xmath78 only an even number of pions is possible and for @xmath80 only an odd number .
this makes the iso - spin projected distributions appear peaked .
table 1 : the mean pion number and variance , and the mean pion number by charge type from nucleon antinucleon annihilation at rest in each of the three iso - spin states with the constraint of four - momentum conservation included .