many efficient sources of quasi - monochromatic hard radiation exploit transmission of ultrarelativistic electrons through periodic structures ( crystals or undulators ) @xcite . for such so - called coherent sources ,
high radiation brightness is relatively easy to achieve by increasing the periodic structure length .
the price to pay , however , is that as the photon emission probability reaches the order of unity , the description of the source performance must regard the possibility of creation of a few photons and electron - positron pairs per passing electron , i.e. , essentially electromagnetic cascading @xcite .
the proper procedure for calculation of electromagnetic multiple particle production at ultra - relativistic energies is via a system of kinetic equations for sequential photon and @xmath0 pair creation , allowing for energy redistribution at each branching . in a non - trivial field of the radiator , this complete procedure
is involved , and generally requires numerical simulation .
fortunately , in a quite typical case when typical photon energies are inferior to the incident electron energy , the calculations may appreciably simplify .
first of all , the probability of pair production is not enhanced so strongly as that of radiation @xcite , and thus can be neglected in the first approximation .
secondly , the negligibility of photon recoils allows treating the electron current as fully determined by the electron s initial conditions , entailing statistical independence of multiple photon emission acts @xcite .
altogether , that opens the possibility for semiclassical description of the cascading process .
however , an extra impediment is that when emitted photon energies belong to gamma - range and @xcite ) . ]
, the measured radiation spectra begin to depend on the photon detection method .
if it was feasible experimentally , in spite of high radiation intensity , to distinguish arrivals of individual photons , the detector signal under soft photon emission condition would be described merely by formula from classical electrodynamics for the radiation energy spectrum .
but narrow beaming of the radiation from ultra - relativistic electrons causes pileup of reactions from different photons in the downstream detector volume , thus affecting the detector signal .
spectrometers for gamma - quanta @xcite today are predominantly based on calorimetry ( measurement of the total energy deposition in a detector ) @xcite , but to capture most of the energy of electromagnetic shower created by the hard photon , the calorimeter transverse size must be in excess of the typical shower lateral spread ( molire radius ) , which amounts a few centimeters @xcite .
that in turn impedes angular resolution of the radiation within its lorentz - contracted emission cone , which might otherwise be employed to eliminate photon pileups .
scheme of gamma - radiation spectrum measurement yielding a multiphoton spectrum .
the spent electron is deflected by a magnet , so only the forward - flying photons hit the calorimeter . ]    in view of appreciable difficulties with eliminating photon pileup effects , a common attitude at laboratory measurements of gamma - radiation spectra is not to pursue the objective of photon counting at all .
then , a single electromagnetic calorimeter is used to measure the _ total _ energy deposited by _ all _ the @xmath1-quanta per event of electron passage through the radiator ( see fig .  [
fig : calorim - scheme ] ) . clearly , the latter method is equivalent to measuring the electron s radiative energy loss spectrum .
furthermore , sometimes a simplified procedure of radiation spectrum measurement is adopted , when spent electron energies are measured , which then truly yields the spectrum of full energy losses by electrons in the radiator .
even though the calorimetric method for spectrum measurement hinders applicability of formulas of classical electrodynamics , those can still be useful for the spectrum calculation .
indeed , the theory predicts that multiple soft photon emission probabilities factorize , expressing through the same classically calculated spectral intensity . to obtain the multiphoton spectrum ,
one thus needs to resum all @xmath2-photon probabilities while holding the total irradiated energy at a preselected value .
technically , the problem is similar to that of fast charged particle energy loss under statistically independent atom ionization acts in the medium @xcite , and some other problems @xcite , whose solutions were attained by integral transform methods . for generic energy loss spectrum by a fast particle in the medium , a laplace transform solution
was known since 1944 @xcite . in application to radiation , it was later utilized in @xcite , with the purpose of investigating soft ( photon - dominated ) electromagnetic showers in crystals and polycrystals , and checking the possibility of the shower length reduction due to coherent radiation effects .
relatively recently , in @xcite this method was used for studying the shape of the soft part of the radiation spectrum from electrons in an amorphous medium .
model solutions for shapes of synchrotron - like radiation spectra in crystals were discussed in @xcite .
there were also several investigations of shapes of radiation spectra in crystals , based on numerical simulation ( usually by monte - carlo ) , with @xcite or without @xcite neglect of @xmath0 pair production as a first approximation . notwithstanding the pretty long history of the problem , there still remain many open questions concerning multiphoton effects on coherent radiation spectrum shapes .
for example , since coherent radiation spectra are typically sharply peaked and discontinuous , the prime task is to investigate how do such features modify under multiphoton emission conditions . in so doing ,
it is also important to incorporate an incoherent bremsstrahlung component , which may require some special treatment , because of its divergent total probability .
next , it is important to give proper study for the behavior of the resummed coherent radiation spectrum in the high intensity ( or photon multiplicity ) limit .
there , one generally expects the spectrum to tend to a gaussian form , but corrections to this form may be significant at practice , and demand evaluation .
besides the abovementioned most urgent problems , there is a number of developments which are desirable to make in the theory .
first , it is expedient to supplement the notion of the resummed spectrum by that of photon multiplicity spectrum , which was proven measurable in experiments at cern @xcite .
secondly , it would be valuable for various purposes to find a possibility to reconstruct the single - photon spectrum from a calorimetrically measured one .
other issues include correspondence with the perturbation theory , separation of energetically ordered and anti - ordered contributions , etc .
the goal of the present paper is to provide a comprehensive formulation for the theory of resummed soft coherent radiation spectra , and with its aid , to tackle problems raised above .
although conditions of coherent radiation in different sources may vary , we concentrate on one of the most typical cases  dipole radiation ( for overview of its conditions see , e.g. , @xcite ) , under the simplifying assumption of approximately harmonic transverse oscillation of the radiating electron ( the neglect of higher harmonics ) .
we incorporate also an incoherent radiation component , which , fortunately , has a simple structure , free from model assumptions . the interference between coherent and incoherent spectral components may be neglected in the first approximation .
therewith , there emerges a fairly model - independent description of radiation , which can embrace diverse practical situations .
this paper is laid out as follows . in sec .
[ sec:2 ] we establish main relations of the theory of soft multiphoton spectra . along with the resummed spectrum itself
, the photon multiplicity spectrum is defined . in sec .
[ sec : coh+incoh ] , the developed generic techniques are applied to the study of coherent and incoherent contributions to the radiation .
modification of various spectral features due to multiphoton effects , as well as interplay of coherent and incoherent radiation components , is investigated .
[ sec : high - intensity ] explores the high - intensity limit , benefiting from the analogy between the multiphoton energy loss and a random walk .
having calculated the multiphoton spectrum in the central region , we also examine the spectrum behavior beyond it , and assess long - range effects due to the incoherent bremsstrahlung component .
[ sec : reconstr ] works out the principles of reconstruction of the single - photon radiation spectrum from the calorimetrically measured one .
the summary and some outlook is given in sec .
[ sec : summary ] .
a few technical issues are relegated to the appendices .
as was mentioned in the introduction , for radiation emitted by an ultra - relativistic electron ( or positron ) , the radiation angles are typically small ( inversely proportional to the electron s lorentz factor ) , wherefore they often remain experimentally unresolved . in the latter case , only the radiation spectrum is measured . for coherent radiation ,
the spectrum is concentrated on a characteristic energy scale @xmath3 , which is often much smaller than the initial electron energy @xmath4 : @xmath5 under those circumstances , the radiation recoils may be neglected , and the radiation process be treated as statistically independent emissions of photons by a classically moving charged particle .
the description of the multiphoton spectrum then proceeds as follows . ] .
let the classically evaluated radiation energy spectrum , integrated over relativistically small emission angles , equal @xmath6 .
the corresponding photon number density could be derived as @xmath7 that might as well be interpreted as the probability distribution for emission of a photon , but only provided its integral @xmath8 which is to correspond to the total emission probability , is @xmath9 .
sometimes the latter condition is fulfilled thanks to the smallness of electron coupling with the electromagnetic field , @xmath10 .
on the other hand , in a long radiator , @xmath11 can accumulate and get arbitrarily large .
as it becomes sizable , proper probabilistic treatment must take into account emission of an arbitrary number of photons .
assuming that the radiation process remains soft and intrinsically semiclassical , photon emission acts can be proven to be statistically independent @xcite , thereby obeying poisson statisticsa word of caution must be sound that poisson statistics holds only for a completely prescribed classical electromagnetic current . at practice , the radiation spectrum can yet depend on the electron impact parameters and momentum dispersion within the initial beam , fluctuations of the intra - crystal field due to thermal vibrations and lattice defects , etc .
therefore , the poisson distributions must be averaged over the beam and target ensembles , but only at the last step of the calculation ( cf .
, e.g. , @xcite ) . in the present paper
, we will refrain from implementation of the averaging procedures , just indicating which of them are most relevant for specific coherent radiation sources . fortunately , in a number of cases we consider , the ensemble averaging effects are negligible . ] : @xmath12 here factors @xmath13 and @xmath14 are consequences of probability conservation during photon generation is conventionally attributed to photon equivalence . ] . from the normalization of the total probability to unity , @xmath15 one readily derives a relation of @xmath13 with the single - photon spectrum : @xmath16 since @xmath13 may be regarded as the 0-th term of the series , it is to be interpreted as the probability of photon non - emission .
quantity @xmath11 also admits an independent physical meaning , being equal to the mean number of emitted photons ( the photon multiplicity ) : @xmath17    from the experimental point of view , completely exclusive spectrum @xmath18 is arguably too detailed . if we ( idealistically ) limit ourselves to a single - inclusive spectrum , where the energy of only one of the photons is measured , the corresponding observable would equal @xmath19 thus reproducing the photon number spectrum , connected with the classical energy spectrum via eq .
( [ primary ] ) .    for the case of calorimetric
measurement , however , the situation is different .
there , the probability of emission of an arbitrary number of photons with definite aggregate energy @xmath20 is measured , taking the form - functions in eqs .
( [ class ] ) or ( [ generic - sum - int ] ) must be replaced by a function of the detector response ( see , e.g. , @xcite ) .
that would be equivalent to ensemble averaging over the final state , which is left beyond the scope of this paper ( see footnote  [ exp - before - aver ] ) . at practice ,
though , the energy resolution of calorimeters is usually high enough for the @xmath21-function approximation in ( [ generic - sum - int ] ) to work satisfactorily . ]
@xmath22 the latter distribution ) are known as compound ( or generalized ) poisson distributions @xcite .
they are usually defined in terms of self - convolutions of the initial probability density @xmath23 it is easy to see that in eq .
( [ generic - sum - int ] ) , @xmath24 is an @xmath2-fold self - convolution of the single - photon spectrum . ]
is already properly normalized : @xmath25 as illustrated in fig .
[ fig : integr - domain ] , in the limit of low radiation intensity , distributions ( [ class ] ) and ( [ generic - sum - int ] ) match , since terms @xmath26 in their series are identical ( fig .
[ fig : integr - domain]a ) .
but the rest of the terms , with @xmath27 , do differ ( cf . fig .
[ fig : integr - domain]b ) .
distribution ( [ generic - sum - int ] ) is often called multiphoton spectrum , and in contrast to that , ( [ primary ] ) is called single - photon spectrum .     integration domains for partial probabilities in eqs .  ( [ class ] ) and ( [ generic - sum - int ] ) . ( a ) component @xmath26 . only the point
@xmath28 is selected by perfect detector of either kind .
( b ) component @xmath29 .
dashed lines mark the region selected by a perfect photon counter and spectrometer .
solid line , the region selected by a perfect calorimeter . ]
another way of representing the results is via a multiphoton energy spectrum @xmath30 in this paper , though , it will be advantageous for us to work with the probability spectral density .
so , in what follows , the term ` spectrum ' will mostly be used in relation to quantities such as @xmath31 and @xmath32 .    before proceeding with the resummation of series ( [ generic - sum - int ] ) , it is expedient to inspect its structure more closely .
first of all , in ( [ generic - sum - int ] ) , upper integration limits for individual photon energies may actually be lowered down to @xmath20 , granted that if energy of any of the emitted photons exceeds @xmath20 , the energy - conserving @xmath21-function will definitely yield zero . at @xmath26 , one
yet needs to ensure that the @xmath21-function falls into the integration interval completely , so the integration upper limit should be written more precisely as @xmath33 . with the forthcoming resummation procedure in mind , it is convenient to let the upper integration limit be the same everywhere , writing : @xmath34    next , even though in eq .
( [ generic - sum - int - upper - omega ] ) upper integration limits for partial probabilities are rendered finite , in integral ( [ w1-def ] ) entering @xmath13 through eq .
( [ w0-def ] ) , the integration extends over an infinite range of @xmath35 .
that fact needs care , since at @xmath36 , quantum corrections due to photon recoils invalidate the statistical independence of photon emission acts . for pure coherent radiation , whose spectrum is strongly suppressed beyond a sufficiently low energy @xmath3 , such a problem may be absent .
however , there often exists an incoherent bremsstrahlung contribution to @xmath31 , scaling as @xmath37 .
therewith , integral @xmath11 logarithmically diverges on the upper limit , whereby factor @xmath38 would nullify the multiphoton radiation spectrum for any finite @xmath20 , if there was no physical end of the bremsstrahlung spectrum , situated at @xmath39 . within the leading logarithmic ( ll ) accuracy ,
the upper integration limit could be merely replaced by @xmath4 : @xmath40 however , the usage of the leading - log accuracy in eq .
( [ poisson ] ) , where @xmath11 is exponentiated , may be unsatisfactory , as long as it introduces an indeterminate overall factor .
thence , it is worth promoting it to the next - to - leading logarithmic ( nll ) accuracy .
that is equivalent to effectively replacing eq .
( [ w1-regee ] ) by @xmath41 where parameter @xmath42 must be ab initio calculable in nll .
a study delivered in appendix  [ app : kappa ] ( though still under the neglect of pair production ) suggests that @xmath43 .
greater accuracy may be unnecessary , since @xmath44 only enters in terms a power @xmath45 with @xmath46 ( see sec .
[ subsec : incoh - bremsstr ] ) .
the incoherent bremsstrahlung contribution also dominates in the opposite extreme @xmath47 , where its behavior @xmath48 makes all the integrals present in eq .
( [ generic - sum - int - upper - omega ] ) diverge on the lower limit .
fortunately , though , the infrared ( ir ) divergence cancels for energy - resummed distributions by virtue of the bloch - nordsieck theorem @xcite , as will be demonstrated in the next subsection .
resummation of series of type ( [ generic - sum - int ] ) proceeds by applying a laplace transform , which reduces the multiple convolution to a simple product of equal single integrals , and the resulting power series resums to an exponential : @xmath49 [ prior to the resummation , we implemented a proper cutoff parameter @xmath50 , in the manner of eq .
( [ w1-nll ] ) ] . an immediate observation from eq .
( [ generic - exponentiation ] ) is that by virtue of factor @xmath51 tending to zero as @xmath47 , the infrared singularity brought by @xmath31 is suppressed , providing the integral convergence on the lower limit if @xmath52 . if @xmath53 is put to zero , eq .
( [ generic - exponentiation ] ) reproduces normalization condition ( [ total - probab-1 ] ) .
the laplace transform is further inverted by computing a mellin integral in the complex @xmath53-plane : @xmath54 here constant @xmath55 is required to be greater than real parts of all singularities of the integrand , otherwise being arbitrary . combining eqs .
( [ generic - mellin - int ] ) and ( [ generic - exponentiation ] ) yields the generic integral representation for the multiphoton radiation spectrum : @xmath56 the appearance in eq .
( [ generic - eq ] ) of a negative term proportional to dirac @xmath21-function does not imply negativity of the resummed spectrum at @xmath57 .
it just cancels the corresponding positive singularity in the contour integral , arising because its integrand tends to @xmath58 as @xmath59 . at @xmath60 ,
the @xmath21-term does not contribute , anyway , but it proves important at derivation of various integral relations involving @xmath32 [ in particular , maintaining the correct normalization ( [ total - probab-1 ] ) ] .
the integral on the r.h.s . of ( [ generic - eq ] )
, @xmath61 has an independent physical meaning , representing the distribution function for the radiating electrons , which is singular at @xmath57 , but is normalized to unity : ) , the integral over @xmath20 extends to infinity , because even if the single - photon spectrum terminates at @xmath62 , the spectrum resummed without account of electron energy degradation after photon emission will actually extend beyond the limit @xmath39 ( see discussion in sec .
[ subsec : incoh - bremsstr ] ) . ] @xmath63 representation ( [ generic - eq ] ) without the @xmath21-term , i.e. holding for function @xmath64 [ or pertaining to case @xmath65 , @xmath66 , was derived by landau @xcite on the basis of a kinetic equation .
next , let us note that contour integral representation ( [ generic - eq ] ) needs not be the only possible one .
the example of representation ( [ generic - sum - int - upper - omega ] ) shows that in the exponentiated laplace transform of @xmath31 , the upper integration limit may be arbitrary , provided it exceeds @xmath20 .
for instance , lowering it right down to @xmath20 leads to the integral representation @xmath67 indeed , expanding here @xmath68 to taylor series and integrating over @xmath53 termwise leads to decomposition ( [ generic - sum - int - upper - omega ] ) .
relation ( [ generic - eq - upperlimitomega ] ) was first noticed in @xcite , wherein it was regarded as only an _ approximate _ consequence of representation ( [ generic - eq ] ) ( in the sense that in the integral over @xmath53 , the contributing @xmath53 are @xmath69 , and hence in @xmath70 the contributing @xmath35 are @xmath71 ) . but here , by deriving both representations ( [ generic - eq ] ) and ( [ generic - eq - upperlimitomega ] ) from the same series , we reveal their exact equivalence .
although the derivation of representation ( [ generic - eq - upperlimitomega ] ) given above is illuminating , it may be desirable sometimes to avoid the use of power series , which may be beset by ir divergences . in that case , one may utilize the following generic statement :    [ lemma1 ] let @xmath72 be real parameters , and @xmath73 be a positive function of real variable , integrable on interval @xmath74 .
also let @xmath75 be an analytic function of complex variable @xmath76 in domain @xmath77 , and such that @xmath78 grows with @xmath79 slower than exponentially .
then , contour integral @xmath80 where @xmath81 is greater than real parts of all singularities of the integrand in the @xmath53-plane , is independent of @xmath82 on the semiaxis @xmath83 .
differentiation of eq .
( [ lemma ] ) by parameter @xmath82 gives : @xmath84 the condition that function @xmath85 grows with @xmath86 slower than exponentially guarantees the exponential decrease of the integrand of ( [ lemma - diff ] ) when @xmath83 .
thus , if the integration contour , defined to pass on the right of all singularities of the integrand , is withdrawn rightwards to infinity , the whole integral exponentially vanishes .
therewith , derivative ( [ lemma - diff ] ) is identically zero , which entails independence of integral ( [ lemma ] ) on @xmath82 under the conditions assumed .    in our case ( [ generic - eq ] ) , ( [ generic - eq - upperlimitomega ] ) , function @xmath75 has to be identified with an exponential @xmath87 ( where @xmath88 furnishes an ir cutoff , if necessary ) .
function ( [ fff ] ) satisfies the conditions of our lemma liberally : since @xmath89 increases at @xmath90 at most logarithmically , @xmath91|$ ] in this limit can only increase as some power of @xmath53 . as a whole , @xmath75 is actually ir - finite , so after unifying @xmath35-integrals in the exponent of @xmath92 , it is safe to put @xmath93 . by choosing in expression ( [ lemma ] ) @xmath82 arbitrarily small or large
, one can deduce eq .
( [ generic - eq - upperlimitomega ] ) from eq .
( [ generic - eq ] ) , or vice versa .
we will also employ lemma [ lemma1 ] for simplification of the exponentiated integrals with specific profiles of @xmath94 in the next section .
( [ generic - eq - upperlimitomega ] ) has an appealing property that apart from the factor @xmath95 the multiphoton spectrum involves only contributions from the single - photon spectrum with @xmath96 .
that looks natural in view of the positivity of contributing photon energies . in this sense
, it seems natural to term eq .
( [ generic - eq - upperlimitomega ] ) ` energetically ordered ' form [ to distinguish it from ` non - ordered ' form ( [ generic - eq ] ) ] .
still , one should mind the existence of factor ( [ w0omega - def ] ) , which depends on contributions from @xmath97 , and is nothing but the probability of non - emission of any photon with energy greater than @xmath20 @xcite ( check : @xmath98 with @xmath99 standing for the heaviside unit step function ) .
so , the effect of contributions from @xmath97 is ` non - dynamical ' , resulting in a rather trivial explicit factor , but nonetheless , due to this factor , the energetic ordering does not hold strictly for multiphoton probability spectrum .
remarkably , in sec .  [ sec : reconstr ] we will encounter representations which may be regarded even as energetically anti - ordered
.      the manifestation of multiphoton effects may be studied in terms of non - linearity of the radiation spectrum dependence on the radiator length ( or the crystalline target thickness ) , @xmath100 ( cf .
@xcite ) . in the simplest approximation ( reasonable for long radiators )
, one may assume that @xmath101 .
then , @xmath100-dependence of @xmath32 may be studied based on expansion ( [ generic - sum - int ] ) , where factor @xmath13 is yet a non - linear function of @xmath11 . expanding the latter factor to series in @xmath11 leads to a double series for @xmath32 .
in fact , a single power series in @xmath100 can be obtained using the following trick .
issuing from contour integral representation ( [ generic - eq ] ) , integrate in the exponent by parts : @xmath102 next , expand the @xmath31-dependent part of the exponential to power series : @xmath103 with @xmath104 now , it is straightforward to integrate over @xmath53 termwise with the use of the identity @xmath105 which yields @xmath106 the latter multiple integral converges even if the mean photon number @xmath11 diverges logarithmically ; its multiple derivatives always exist , as well .
in contrast to eq .
( [ generic - sum - int ] ) , however , terms of series ( [ dwn - pert - ser ] ) need not be everywhere positive .
that is chained to the fact that multiple photon emission leads to a redistribution of the spectrum , not only to pileup of events . from the standpoint of qed , negative terms in the perturbative expansion
typically arise from interference of loop ( photon reabsorption ) diagrams with lower - order ones @xcite .
( [ dwn - pert - ser ] ) indicates that for an observable such as the calorimetrically measured spectrum , the structure of those corrections assumes a simple generic form .
correspondence of ( [ dwn - pert - ser ] ) with expansion ( [ generic - sum - int ] ) may be established by differentiating in ( [ dwn - pert - ser ] ) the delta - function under the integral sign by formula @xmath107 and subsequently integrating by parts over all @xmath108 .
it should be minded thereat that the endpoint terms are non - vanishing , producing the powers of the total single - photon emission probability @xmath11 present in eq .
( [ generic - sum - int ] ) in the exponent .
however , if @xmath11 diverges , the integration by parts is no longer possible , wherewith series ( [ generic - sum - int ] ) does not exist , but eq .
( [ dwn - pert - ser ] ) remains valid , anyway .
characterization of compact probability distributions is often conducted in terms of their moments .
spectrum @xmath32 , however , is not normalized to unity , so it can not be used straightforwardly as a weighting distribution . in capacity of a normalized probability distribution
, one should either use electron distribution function ( [ = pi ] ) , or the rescaled radiation spectrum @xmath109 .
we will adopt the first option , which leads to simpler final results .
the mean emitted photon energy for our resummed spectrum is computed , e.g. , by termwise integration in eq .
( [ generic - sum - int ] ) : @xmath110\equiv\int_0^\infty d\omega \omega\frac{dw}{d\omega}\nonumber\\ & = & e^{-w_1}\sum_{n=1}^\infty\frac{n}{n!}w_1^{n-1}\overline{\omega_1\!}=\overline{\omega_1\!}.\end{aligned}\ ] ] thus , it coincides with that for the single - photon spectrum , @xmath111 that is quite natural , inasmuch as once photon energies in all the events are summed over , it no longer matters whether they were measured by photon counters or calorimeters .
[ in particular , that justifies the use of classical formulas for the rate of radiative energy losses ( see , e.g. , @xcite ) under conditions of multiple photon emission . ]    to further characterize the spectrum width , asymmetry , and so on , one can refer to higher moments about the mean value ( central moments ) , presuming their existence .
those are conveniently calculated based on representation ( [ generic - eq ] ) : @xmath112\qquad\qquad\nonumber\\ & = & \left(-\frac{\partial}{\partial s}-\overline{\omega_1\!}\right)^n e^{\int_0^{e } d\omega_1 \frac{dw_1}{d\omega_1}\left(e^{-s\omega_1}-1\right ) } \bigg|_{s=0}\nonumber\\ & \equiv&(-1)^n\frac{\partial^n}{\partial s^n}e^{\int_0^{e } d\omega_1 \frac{dw_1}{d\omega_1}\left(e^{-s\omega_1}-1+s\omega_1\right ) } \bigg|_{s=0}.\,\,\end{aligned}\ ] ] the above representation implies that exponential @xmath113 serves as a generating function @xcite for central moments , while @xmath114 offers a generating function for ordinary moments . applying eq .
( [ higher - moments ] ) for lowest @xmath2 , one finds @xmath115 @xmath116 @xmath117 etc . relations ( [ mean - omega ] ) , ( [ mean - omega2-alt ] ) ( historically known as campbell s theorem @xcite ) were discussed in @xcite with the objective of extracting information about the single - photon distribution from the resummed one .
( [ mean - omega3 ] ) , ( [ mean - omega4 ] ) may serve for the same purpose .
in addition , in sec .  [ sec : reconstr ] we will formulate a reconstruction procedure for the complete spectrum , not resorting to the notion of moments .
a predicament with the use of the above quoted moments is that with the account of an incoherent bremsstrahlung component ( sec .
[ sec : coh+incoh ] ) , they all diverge in the ultraviolet . therefore , they may seem to be only useful in the context of pure coherent radiation .
nonetheless , later on we will show that treatments of coherent and incoherent radiation components can always be separated , permitting the usage of spectral moments for pure coherent component regardless of the presence of an incoherent one .      as was mentioned in the introduction
, it is rather difficult in a multiphoton event to practically pin down individual energies of all the photons , i.e. , to measure joint probability distribution @xmath18 .
but there exist relatively compact experimental methods allowing to receive partial information about the photon number content .
for instance , placement of a thin converter and an ionizing particle counter upstream the calorimeter allows measuring the mean number of photons ( the photon multiplicity ) as a function of the total energy @xmath20 deposited in the calorimeter
see fig .
[ fig : multty - measur - scheme ] and , e.g. , @xcite . in this subsection
, we will extend our resummation procedure to description of the photon multiplicity spectrum .    in the spirit of eqs .
( [ mean - n - def ] ) and ( [ generic - sum - int ] ) , to construct the photon multiplicity at a given @xmath20 , one must incorporate in ( [ generic - sum - int ] ) a weighting factor @xmath2 equal to the number of emitted photons , whereupon divide the resulting sum by the unweighted expression : @xmath118 since all the weights @xmath2 in this sum are greater than unity , so must be their mean number : @xmath119 in particular , when @xmath120 , integration phase space volumes for @xmath27 terms shrink to zero , wherefore only term @xmath26 survives .
this term cancels with the correspondent term in the denominator , yielding @xmath121    the sum entering eq .
( [ mean - n - omega - def ] ) can again be evaluated by the laplace transform method , as follows :    [ mean - n - omega - contour ] @xmath122    [ the latter equality may be inferred from lemma  [ lemma1 ] with @xmath123 .
this may also be expressed as a convolution of the single - photon spectrum with the resummed one : @xmath124 here the last term stems from the @xmath21-function term of eq .
( [ generic - eq - upperlimitomega ] ) .    unlike @xmath32 , however , the multiplicity spectrum is not an ir - safe quantity .
but fortunately , the ir divergence of @xmath11 affects only the constant , @xmath20-independent part of @xmath125 .
indeed , the contribution from vicinity of the lower integration limit in the convolution term in ( [ mean - n - omega - convol ] ) is @xmath126 , and so is the l.h.s .
hence , at @xmath127 , @xmath20-dependent factors in both sides of eq .
( [ mean - n - omega - convol ] ) cancel .    to isolate the ir divergence of @xmath128 explicitly
, one can rewrite eq .
( [ mean - n - omega - convol ] ) , e.g. , as @xmath129 here the r.h.s .
is ir - safe ( thence we put there @xmath93 ) , while in the l.h.s .
, the dependence of @xmath128 on the ir cutoff @xmath88 must be additive , to cancel with the additive dependence on ir cutoff of the integral term .
on the other hand , the dependence of on @xmath50 for @xmath128 must disappear , as will be demonstrated in the next section .
generic resummation formulas established in the previous section are valid for radiation characterized by an arbitrary single - photon profile @xmath94 . in
what follows , we will mostly concentrate on a category of coherent radiation spectra involving an admixture of incoherent radiation . in that case ,
the salient features shared by the single - photon spectra are :    * spectral discontinuities at finite @xmath20 in the coherent radiation component ( coherent emission edges ) ; * a ` tail ' towards large @xmath20 brought by the incoherent bremsstrahlung component ; * an infrared divergence , originating from the incoherent bremsstrahlung , or from non - zero net deflection by coherent fields ( see , e.g. , @xcite ) . for simplicity ,
in this paper we will be only considering the first case .
our first task is to investigate how do the above mentioned features modify under multiphoton emission conditions .
since multiphoton effects interplay nonlinearly , it is expedient first to analyze them for pure coherent and incoherent components separately . combining them
, we can thereupon distinguish effects of their mutual influence .
let us first scrutinize the case of purely incoherent bremsstrahlung .
its single - photon spectrum must be proportional to bethe - heitler s spectrum of radiation at a single collisionthis assumption admittedly neglects dense - medium effects such as transition radiation due to the change of the dielectric density , and lpm radiation suppression by multiple scattering @xcite .
their neglect is justified under typical conditions @xmath130  mm , @xmath131  gev , @xmath132  mev , although some of those conditions can be relaxed at the expense of the others .
the influence of forward transition radiation ( noticeable when @xmath133 mm ) , and of lpm suppression on resummed radiation spectra was studied in @xcite . ] @xcite .
moreover , in the soft radiation limit we are restricted to , the single - photon spectrum reduces to the well - known semiclassical form @xcite @xmath134 where @xmath135 is proportional to the target thickness @xmath100 . in an amorphous matter , constant @xmath135 equals @xmath136 where the radiation length @xmath137 depends on the atomic density @xmath138 and ( net ) nucleus charge @xmath139 @xcite .
a similar estimate applies also for the incoherent component of coherent bremsstrahlung in a crystal . for the case of positron channeling , if their close collisions with atomic nuclei are sufficiently rare , factor @xmath140 should be replaced by the net interplanar electron density @xmath141 . finally , in an undulator , @xmath135 is proportional to the density of the residual air , though in high vacuum that can be negligible .    as was argued in the end of sec .
[ subsec : stat - indep - emission ] , the integral of ( [ dw1inc ] ) , which is logarithmically divergent in the uv , must be cut off at some photon energy commensurable with the electron energy : @xmath142 , with @xmath143 .
we implement this cutoff directly for approximation ( [ dw1inc ] ) , adopting a simplified model : @xmath144 where @xmath145 is the heaviside function .
integration contours for resummed spectra .
line @xmath146 , used in eq .  (
[ generic - eq ] ) , ( [ generic - eq - upperlimitomega ] ) .
deformed contour @xmath147 , used in eqs .
( [ eulergamma - def ] ) , ( [ dw - first - interv ] ) .
closed loop contour @xmath148 , applicable for pure coherent radiation ( see sec .  [ subsubsec : coh - in - fund - interv ] ) . ]    a convenient starting point for evaluation of the resummed spectrum corresponding to single - photon spectrum ( [ dw1incoh - cutoff ] ) is the energetically ordered contour integral representation , eq .
( [ generic - eq - upperlimitomega ] ) .
therein , the @xmath21-term is suppressed , since @xmath149 , and for the present case it will be neglected . once
( [ dw1incoh - cutoff ] ) is inserted to eq .
( [ generic - eq - upperlimitomega ] ) , ) . that depends on how large are relative fluctuations of parameter @xmath135 per electron passage .
if the radiation formation length @xmath150 is @xmath151 , which is plausible under conditions of footnote [ foot : no - medium - effects ] , many electron scatterings in the target contribute to the radiation independently , so fluctuations of @xmath135 are relatively small , justifying the use of resummation for the averaged spectrum .
averaging procedure for multiphoton radiation in a thin target will be considered elsewhere . ]
all @xmath20- and @xmath50-dependencies factor out after a simple rescaling of the integration variables , giving @xmath152 with @xmath153 to evaluate @xmath154 , one may regularize @xmath155 on the lower limit , then split it in two integrals , and by lemma  [ lemma1 ] , extend the upper limit of the integral involving @xmath156 to infinity .
thereafter , the ir regulator may again be sent to zero , and the sum of @xmath157-integrals evaluates : @xmath158 the integrand of the latter contour integral is vanishing as @xmath159 , so at complex infinity the contour can be turned to the left from the imaginary axis ( contour @xmath147 in fig .
[ fig : contours ] ) . using representation @xmath160 for the euler gamma function , and definition @xmath161 for euler s
constant , we get @xmath162 function ( [ fa ] ) is manifestly positive and monotonously decreasing ( see fig .
[ fig : phi , nu ] ) , its derivative in the origin being zero : @xmath163 other models for the single - photon spectrum suppression at large @xmath35 yield equivalent results ) one employs a slow exponential factor , @xmath164 with @xmath165 , the insertion of this form to representation ( [ generic - eq ] ) , with the upper limit of @xmath35-integral replaced by infinity , gives @xmath166 apart from factor @xmath167 ( which is close to unity when @xmath168 ) , this structure coincides with eq .
( [ 13 ] ) , once one identifies @xmath169 .
there , the single - photon spectrum is expressed by a gamma - distribution @xcite of index 0 , while the resummed spectrum is again a gamma - distribution , but of index @xmath135 .
that is a well - known example of functional stability more general than lvy stability ( see , e.g. , @xcite ) .
one can as well arrive at result ( [ 13 ] ) from the side of lvy - stable densities @xcite , for which the single - photon distribution is strictly scale - invariant : @xmath170 ; however , lvy distributions at finite @xmath21 , in contrast to ( [ 13 ] ) or ( [ gamma - distr - a ] ) , do not reduce to a product of functions of a single variable . ] .    regarding the shape of the resulting distribution , two remarks are in order .
first , spectrum ( [ 13 ] ) is suppressed compared to single - photon spectrum ( [ dw1inc ] ) everywhere in the region @xmath171 . at @xmath172
, there formally develops an enhancement of the spectrum , but it is not accurately described by eq .
( [ 13 ] ) , insofar as at @xmath173 quantum effects enhance , while at @xmath39 , the spectrum must strictly terminate at all . hence , there ansatz ( [ dw1incoh - cutoff ] ) for the single - photon spectrum breaks down , so our model description , though self - consistent formally , is inadequate in domain @xmath173 . with this reservation ,
the uniform suppression of the resummed spectrum at @xmath171 does not contradict , e.g. , the conservation of the mean photon energy at resummation [ eq .
( [ mean - omega ] ) ] .
graphs of functions @xmath154 [ eq .  ( [ fa ] ) ] and @xmath174 [ eq .  ( [ nu - through - digamma ] ) ] .
our physical assumptions require that @xmath175 ( see text ) . ]
secondly , the dependence of the resummed spectrum on @xmath44 enters eq .
( [ 13 ] ) through parameter @xmath142 , resulting in an overall factor @xmath176 .
thus , in order for result ( [ 13 ] ) not to be strongly sensitive to the actual value of @xmath44 , one needs fulfillment of condition @xmath177 ( otherwise , full kinetic equations for the electromagnetic cascade will need to be solved ) . since the accuracy of eq .
( [ 13 ] ) is @xmath178^a=1+\mathcal{o}(a)$ ] , function @xmath179 may be safely approximated by @xmath180 .
then , one readily checks that @xmath181 in accord with eq .
( [ total - probab-1 ] ) ( where , again , @xmath127 , owing to ir divergence of @xmath182 ) . property ( [ incoh - norm ] ) implies that despite the smallness of @xmath135 , function ( [ 13 ] ) is not uniformly small , but rather is similar to a @xmath21-function , peaking at @xmath120 where function ( [ 13 ] ) blows up : @xmath183 to prove eq .  ( [ delta - lim - ato0 ] ) , and estimate the narrowness of @xmath184 as a function of parameter @xmath135 , it may be expedient to evaluate a median energy @xmath185 at which there accumulates half probability of the resummed incoherent radiation , i.e. @xmath186 that yields @xmath187^{-1/a}e\lll e_e,\ ] ] which at small @xmath135 is exponentially small , thereby validating relation ( [ delta - lim - ato0 ] ) .
finally , the ( regularized ) photon multiplicity spectrum can be evaluated for the present case via eq .
( [ mean - n - omega - convol - regul ] ) . inserting ( [ dw1inc ] ) and
( [ 13 ] ) to ( [ mean - n - omega - convol - regul ] ) , and evaluating the integral , we get : @xmath188 with    @xmath189\\ & = & 1+a\left[\psi(1)-\psi(1+a)\right]\label{nu - through - digamma}\\ & = & 1-\frac{\pi^2}6 a^2+\mathcal{o}(a^3),\end{aligned}\ ] ]    involving digamma - function @xmath190 .
hence , photon multiplicity spectrum grows with @xmath20 strictly logarithmically , and proves to be independent of the ultraviolet ( uv ) cutoff @xmath50 .
the graph of function @xmath174 is shown in fig .
[ fig : phi , nu ] . at low @xmath135 , it is close to unity , which is natural , because no more than one photon is typically emitted per passing electron .
next we consider the case of pure coherent radiation .
it is principally different from the incoherent radiation case considered above , for as it involves an integrable ( ir- and uv - safe ) , although discontinuous single - photon spectrum .
physically , that may correspond to ( gamma - ray ) undulator radiation @xcite , or channeling radiation @xcite . as for coherent bremsstrahlung , usually containing an admixture of incoherent radiation component , too , it is more pertinent to the combined radiation case discussed in the next subsection .    in what concerns application to channeling radiation , though
, it should be stressed that the spectrum of the latter strongly depends on impact parameters of the charged particles , even when non - channeled passage events are rejected .
hence , at passing to multiphoton spectra , at first we have to perform the resummation for a definite impact parameter value , and only at the final step average the result over impact parameters .
we can not indulge into such specialized procedures here , leaving them for future studies .
so , in application to channeling radiation our results in this paper will only be preliminary .
what we wish to reflect in our present study is the transverse oscillatory motion of the radiating particle with respect to the direction of high longitudinal momentum , taking place in all the abovementioned cases .
thereat , the radiation spectrum shape depends on the particle oscillation harmonicity . in the simplest case of purely harmonic and small - amplitude oscillatory motion ,
the electromagnetic radiation spectrum has the structure ( see appendix  [ app : coh - rad - sources ] ) @xmath191 where coefficient @xmath192 in a straight crystal / undulator is proportional to the radiator length @xmath100 , and the profile function reads @xcite    [ f - def ] @xmath193    [ the latter form demonstrates the function positivity and symmetry with respect to midpoint @xmath194 . in case if aperture collimation is imposed on the photon beam , in view of the unambiguous relation between the photon energy @xmath35 and its emission angle [ see eq .
( [ omega - theta - corresp ] ) ] , function @xmath195 in eq .
( [ dw1coh - real ] ) will change , but still remain finite everywhere .    if the transverse oscillatory motion of the electron / positron happens to be highly anharmonic , or is already relativistic , higher harmonics in the single - photon spectrum develop , which can generally be incorporated as @xmath196 in particular , with the increase of the positron energy beyond @xmath197 gev , higher harmonics of channeling radiation in a crystal proliferate , and ultimately render the spectrum a synchrotron - like appearance .
this case may also be regarded as universal ; it was dealt with in papers @xcite .
henceforth , we will restrict ourselves to the simplest case of a single - harmonic dipole radiation described by eq .
( [ dw1coh - real ] ) . while in our equations
the shape of @xmath198 will be handled as generic ( to allow for the possibility of aperture collimation ) , in graphical illustrations specific form ( [ f - def ] ) will be used throughout .
( a ) multiphoton probability spectra computed for dipole harmonic coherent radiation defined through eqs .
( [ dw1coh - real][f - def ] ) ( solid curves ) .
red curve , corresponding to intensity parameter @xmath199 , @xmath200 [ small deviations from the single - photon spectrum ( [ f - def ] ) ] .
green curve , the same for @xmath201 ( the highest summit reached by the fundamental maximum ) .
blue curves , @xmath202 ( onset of high - intensity regime ) . dashed blue curve , gaussian approximation ( [ gauss ] ) .
dot - dashed blue curve , corrected gaussian approximation ( [ chebyshev ] ) .
( b ) multiphoton energy spectrum @xmath203 for the same single - photon spectrum as in ( a ) , and intensity parameter values @xmath204 , 2.5 , 6 .
the development of secondary maxima for the energy spectrum is more spectacular than for the probability spectrum @xmath205 .
( c ) photon multiplicity spectrum described by eqs .
( [ mean - n - omega - contour ] ) , ( [ dw1coh - real][f - def ] ) .
solid curves correspond to same parameters as in ( a ) .
dashed blue line is high - intensity approximation ( [ mean - n - at - gauss ] ) . ]    inserting entries ( [ dw1coh - real][f - def ] ) to mellin integrals ( [ generic - eq ] ) and ( [ mean - n - omega - contour - non - casual ] ) , and computing them numerically , one can trace the evolution of the resummed spectrum and the multiplicity spectrum with the increase of the intensity .
this is illustrated in fig .
[ fig : pure - coh - spectrum ] for a few values of @xmath206 .
the trends observed in those figures are : smoothening of the spectral shape with the increase of the intensity , the existence of an upper bound for the probability distribution ( fig .
[ fig : pure - coh - spectrum]a ) , the appearance of a second maximum in the energy spectrum at moderate intensity ( fig .
[ fig : pure - coh - spectrum]b , green curve ) , and monotonic increase of the photon multiplicity spectrum with @xmath20 ( fig .
[ fig : pure - coh - spectrum]c ) .
explaining the origin of those features will occupy us next .      when the single - photon spectrum is described by eq .
( [ dw1coh - real ] ) , at sufficiently low radiation intensity , the multiphoton spectrum will as well be concentrated within the interval @xmath207 ( termed fundamental )  see fig .
[ fig : pure - coh - spectrum]a , red and green curves .
so , it seems natural to begin with evaluation of the multiphoton spectrum in the fundamental interval .
the calculation is straightforward based on representation ( [ generic - eq ] ) , written for finite @xmath208 as @xmath209 here we dropped the @xmath21-term for @xmath60 , and by virtue of lemma  [ lemma1 ] , the upper limit in the integral over @xmath35 was replaced by infinity , without the necessity to employ the heaviside step function , which therefore was omitted . if @xmath195 is a polynomial , like ( [ f - def ] ) , the integral in the exponent of ( [ wcoh-1stinterv ] ) evaluates as a finite - order polynomial in @xmath210 : @xmath211 where @xmath212 stands for the @xmath213-th derivative of @xmath198 w.r.t . its argument . inserting ( [ polynom - s ] ) to ( [ wcoh-1stinterv ] ) , and rescaling the integration variable to @xmath214
, yields @xmath215 benefiting from the exponential decrease of the integrand at @xmath159 , along with the absence of cuts in the complex @xmath216-plane , we replaced the infinite integration contour by a closed loop encircling the origin , where the only singularity of the integrand resides ( contour @xmath148 in fig .
[ fig : contours ] ) .
[ [ infrared - limit . ] ] infrared limit .
+ + + + + + + + + + + + + + +    despite the closed integration contour , and representation of the exponent in the integrand by only a few power terms , integral  ( [ dwc - oint ] ) generally does not permit classification in terms of basic special functions .
, i.e. , function @xmath198 was constant
. such a model may be suitable for quick estimates , but generally is not acceptable numerically . ] but it is instructive at least to obtain its expansion about endpoint @xmath217 . there
, one of the exponentials tends to unity , @xmath218 , and if the second exponential is expanded into laurent series @xmath219 the only term surviving after the loop integration is @xmath220 , leaving the result @xmath221 this equation shows that in the limit @xmath120 , the single - photon probability dominates , albeit in conjunction with photon non - emission probability @xmath222 . for the photon multiplicity spectrum , eq .
( [ coh - omegato0 ] ) , when inserted to eq .
( [ mean - n - omega - convol ] ) leads to eq .
( [ mean - n - to1 ] ) .    to derive an @xmath223 correction to eq .
( [ coh - omegato0 ] ) , one needs to retain also a linear term in the expansion of @xmath224 : @xmath225\!\!+\!\mathcal{o}(\omega^2)\!\right\}\!.\ ] ] for our specific example ( [ f - def ] ) , @xmath226 .
but for a sufficiently high intensity , specifically @xmath227 the coefficient at the linear term in ( [ coh - omegato0-corr ] ) turns positive , wherewith the resummed spectrum in the whole fundamental interval becomes monotonously increasing
. that may be associated with the inception of high - intensity regime , which will be examined in the next section .    from eqs .
( [ coh - omegato0-corr ] ) and ( [ mean - n - omega - convol ] ) , one similarly infers an expression for the photon multiplicity spectrum @xmath228 this shows that the slope in the origin does not depend on @xmath229 , and is always positive , in accord with eq .
( [ mean - n - gtr1 ] ) .
the increase of the slope in the origin with the increase of parameter @xmath192 agrees with fig .  [
fig : pure - coh - spectrum]b .
[ [ par : fund - max ] ] fundamental maximum .
+ + + + + + + + + + + + + + + + + + + +    the opposite endpoint of the fundamental interval , @xmath230 , is of prime practical significance , as long as it represents the global spectral maximum at moderate radiation intensity . substituting @xmath231 in eq .
( [ dwc - oint ] ) leaves @xmath232 effectively , this is a function of a single dimensionless parameter @xmath233 , inasmuch as all @xmath234 . in fig .
[ fig : intensity - in - max ] , the dependence of ( [ max - int - vs - w1c ] ) on @xmath206 [ with @xmath198 defined by eq .
( [ f - def ] ) ] is displayed by the solid gray curve .
naturally , with the increase of @xmath192 , at first it rises proportionally , but eventually the rise halts and ends up with a decrease , formally owing to factor @xmath222 .
the ultimate suppression of the resumed spectrum at any fixed energy @xmath20 is not surprising , given the saturation of total probability ( [ total - probab-1 ] ) on one hand , and the spread of the spectrum to higher @xmath20 on the other hand .    to estimate the maximal summit of function ( [ max - int - vs - w1c ] )
, one can adopt the following approach . to capture the first bend of the nonlinear @xmath192-dependence ,
expand the contour integral in ( [ max - int - vs - w1c ] ) to second order in @xmath192 ( which is equivalent to keeping in eq .
( [ generic - sum - int ] ) the first two terms ) : @xmath235\nonumber\\ % \qquad\equiv e^{-\omega_0b\int_0 ^ 1dzp(z)}\left[bp(1)+\omega_0b^2\frac12\int_0 ^ 1dzp(z)p(1-z)\right]\\ & = & p(1)e^{-b\omega_0\xi}\left(b\omega_0+\frac{b^2\omega_0 ^ 2}2\eta\right),\,\,\label{40}\end{aligned}\ ] ] where @xmath236 @xmath237 the behavior of approximation ( [ 40 ] ) is illustrated in fig .
[ fig : intensity - in - max ] by the dashed gray curve .
the maximum of ( [ 40 ] ) can be explicitly found by differentiation by @xmath206 : it is situated at @xmath238 for @xmath198 described by eq .
( [ f - def ] ) ( wherewith @xmath239 , @xmath240 ) , eq .
( [ loc - max ] ) yields @xmath241 , at which @xmath242 . for comparison , the maximum of exact expression ( [ max - int - vs - w1c ] )
is achieved at @xmath243 , and amounts 0.84 ( see fig .
[ fig : intensity - in - max ] , solid gray curve ) .
hence , at any radiation intensity , the probability spectrum @xmath32 of multiphoton coherent radiation at practice never exceeds unity is such that @xmath157 is exceptionally small , while @xmath244 is sizable .
those are typical conditions for strong aperture collimation .
but that case is rather trivial , because multiphoton effects are then suppressed as a whole . ]  in marked contrast with the unlimited growth of the single - photon spectrum .
[ at that , the energy spectrum @xmath245 can indefinitely grow , @xmath246  see sec .
[ subsubsec : centr ] , eq .
( [ gauss ] ) ] .     the spectral intensity at the fundamental maximum ( @xmath230 ) vs. the radiator length ( in units of @xmath206 ) , for a single - photon spectrum described by eq .
( [ f - def ] ) .
solid gray curve , pure coherent radiation [ eq .
( [ max - int - vs - w1c ] ) ] . dashed gray curve , approximation  ( [ 40 ] ) .
solid black curve , the combination of coherent and incoherent radiation [ eq .
( [ max - dw ] ) ] in proportion @xmath247 , for @xmath248 ( @xmath249 ) .
dashed black curve , approximation  ( [ corr-1fund - max ] ) . ]      besides peaking at @xmath231 , the coherent radiation spectrum encounters at this point a discontinuity , which is actually easier to evaluate . for our simplified case ( [ f - def ] ) , the discontinuity of the single - photon spectrum ( defined to be a positive quantity ) just equals the height of its maximum : @xmath250 with the account of emission of several photons , the spectrum beyond @xmath3 deviates from zero , so higher terms of series ( [ generic - sum - int ] ) demand consideration . but
all terms with @xmath27 prove to be continuous - th term of ( [ generic - sum - int ] ) the integrand function @xmath251 is discontinuous on every face of an @xmath2-cubic domain , but the entering @xmath21-function restricts the integration domain to a slicing plane @xmath252 , which is oblique and nowhere parallel to any face of the cube for @xmath27 ( see fig .  [
fig : integr - domain]b ) .
thereby , for @xmath253 no discontinuities in the @xmath20-dependence can arise . ] .
hence , the discontinuity of the total spectrum just equals that of the first term .
this term yet involves the total photon non - emission probability factor @xmath254 with integer @xmath27 , the resummed spectrum remains continuous , but its @xmath255-th derivative encounters a discontinuity ( as can be noticed already from figs .
[ fig : pure - coh - spectrum ] ) .
that also follows from the analysis of the phase space available for integration for terms of series ( [ generic - sum - int]a , b ) .
furthermore , even in case if function @xmath198 is vanishing at @xmath256 , in representation ( [ generic - sum - int ] ) the function discontinuity at @xmath231 cancels , but discontinuities of its derivatives persist .
we leave the proof of those statements to the reader .    as regards the discontinuities of @xmath257 , fig .
[ fig : pure - coh - spectrum]c shows that they are similar to those of @xmath205 , but have opposite sign [ because @xmath257 contains @xmath258 in the denominator ] . at low radiation intensity
, the multiplicity spectrum exhibits a step - like behavior , approximately amounting to the smallest integer greater than @xmath259 .
that is traced to the fact that in an interval @xmath260 @xmath261 , the @xmath2-photon component dominates . as the intensity increases , function @xmath257 smoothens out , given that a @xmath2-photon component becomes competitive with fewer - photon components even in intervals where those components are not yet extinct .
having explored separately the shapes of resummed coherent and incoherent radiation spectra , we are now in a position to examine their nonlinear interplay .
the physical example when those components are both significant is coherent bremsstrahlung , occurring when an electron or positron crosses a family of crystalline planes at an above - critical ( although small ) angle .
thereat , the continuous interplanar potential acts on the electron periodically , evoking coherent radiation @xcite . at the same time
, the incoherent scattering on atomic nuclei in the planes causes incoherent bremsstrahlung , similar to that in amorphous matter . in a satisfactory approximation
, the single - photon radiation spectrum may be expressed as a sum of two non - interfering , because it is ir- and uv - safe , but that would complicate the structure of @xmath198 . on the other hand , if there is a coherent contribution to the electron r.m.s .
deflection angle , it ought to be included to @xmath262 . ] parts : @xmath263 where @xmath262 is given by eq .
( [ dw1inc ] ) . in graphical illustrations
, we will continue using eqs .
( [ dw1coh - real][f - def ] ) for @xmath264 , i.e. , neglect higher harmonics . in some cases
, those can be small indeed even for coherent bremsstrahlung ( ` one - point ' spectra , see @xcite ) . at the same time
, we saw in the previous section that the second maximum in the multiphoton spectrum can be generated even in the absence of secondary harmonics in the single - photon spectrum , so it will be expedient to check whether their effect is anyhow affected by incoherent radiation . for what concerns the ratio of magnitudes of coherent and incoherent radiation components , for coherent bremsstrahlung it obeys a relation @xmath265 with @xmath266 the misalignment angle between the electron momentum and the co - oriented family of atomic planes , and @xmath267 [ cf .
( [ a - amorph ] ) and ( [ dw - cb - fmax ] ) ] . at typical @xmath268 , ratio ( [ a / bomega0 ] ) is pretty small .
nevertheless , effects of incoherent radiation in the spectrum may be quite noticeable .
that is illustrated by fig .
[ fig : high - intensity ] , where in spite of rather large initial ratio @xmath269 , the incoherent bremsstrahlung plateau in the multiphoton energy spectrum @xmath245 is only a few times lower than the spectrum height in the maximum ( fig .
[ fig : high - intensity]b ) .
handling the effects of incoherent radiation proves to be most convenient with the aid of a convolution representation derived below .
once decomposition ( [ coh+inc ] ) is inserted to eq .
( [ generic - eq ] ) , the exponential in the integrand splits into a product of two factors , one depending on @xmath270 , and the other on @xmath271 . expressing then each of those factors via eq .
( [ generic - exponentiation ] ) through laplace transforms of the corresponding resummed radiation spectra , and doing the contour integral , we arrive at a convolution relation . ] : @xmath272 the term in the second line vanishes if @xmath182 divergesat practice , when the incoherent component receives a physical ir cutoff , factor @xmath273 may be not really small , so the last term in ( [ wcoh - times - winc ] ) may remain sizable . to be more quantitative , @xmath274 with @xmath88 the ir cutoff value .
but for cutoffs of different nature , and practical energies @xmath4 , ratio @xmath275 can hardly exceed @xmath276 , wherewith @xmath277 .
hence , @xmath273 can really be vanishing only provided @xmath278 .
we will restrict our consideration to the idealized case when this condition is met .
the extension to finite @xmath182 is straightforward , and basically resembles the situation described in the previous subsection . ] , whereupon expression ( [ wcoh - times - winc ] ) becomes linear in @xmath279 .
still , the resulting convolution is non - vanishing at @xmath280 , because , as we saw in sec .
[ subsec : incoh - bremsstr ] , the multiphoton incoherent radiation spectrum in this limit tends to a @xmath21-function , not to zero [ see eq .
( [ delta - lim - ato0 ] ) ] .
more precisely , if median @xmath185 defined by eq .
( [ omegaincoh12 ] ) proves to be much smaller than @xmath3 , convolution ( [ wcoh - times - winc ] ) will be close to @xmath205 .    observing further that in eq .
( [ wcoh - times - winc ] ) @xmath281 is a smooth function , whereas the coherent radiation spectrum @xmath282 is not , it may be beneficial to integrate in ( [ wcoh - times - winc ] ) by parts : @xmath283 invoking here eq .
( [ coh - omegato0 ] ) along with identity @xmath284 , the convolution relation recasts @xmath285\qquad\qquad\quad\nonumber\\ -\phi(a)\int_0^{\omega}d\omega ' \left(\frac{\omega'}{e}\right)^a\frac{d}{d\omega'}\left(\frac{dw_\text{c}}{d\omega''}\bigg|_{\omega''=\omega-\omega'}\right).\end{aligned}\ ] ] the merit of the latter representation is that it picks up discrete contributions from points at which @xmath282 encounters discontinuities , transforming to @xmath21-function - like terms for its derivative .
( a ) multiphoton radiation spectrum involving coherent component with magnitude @xmath201 , and incoherent component with @xmath286 .
the uv cutoff parameter is chosen at @xmath248 .
solid green curve , the exact result obtained by numerical evaluation of the contour integral .
long - dashed black curve , corrected low-@xmath20 approximation ( [ dw - small - omega - nlo ] ) . dotted black curve , approximation next to the coherent emission edge , eq .
( [ smearing - of - disc ] ) .
dot - dashed black curve , large-@xmath20 asymptotics ( [ simple - large - omega ] ) .
( b ) the energy spectrum for same conditions as in ( a ) .
short - dashed black curve , low-@xmath20 approximation ( [ dw - small - omega ] ) .
dot - dashed black curve , large-@xmath20 asymptotics ( [ simple - large - omega ] ) .
( c ) the photon multiplicity spectrum for same parameters as in ( a , b ) , and ir cutoff @xmath287 . solid green curve , exact result .
dot - dashed black curve , large-@xmath20 approximation ( [ n - omega - asympt - large ] ) . ]    in the limit @xmath120 , the integral term in eq .
( [ convol - byparts ] ) , vanishes due to the integration interval shrinkage , leaving @xmath288.\ ] ] clearly , here the leading term @xmath289 represents the contribution due to incoherent bremsstrahlung , but it receives an extra suppression by factor @xmath222 depending on the coherent radiation .
hence , the incoherent component can not be asserted to dominate in this limit alone , in contrast to the case of large @xmath20 ( see sec .
[ subsubsec : large - omega - asympt ] below ) .
the behavior of asymptotic approximation ( [ dw - small - omega ] ) is illustrated in fig .  [ fig : high - intensity]b by the short - dashed black curve .    in a similar fashion
, one can derive a convolution relation for the multiplicity spectrum , which reads @xmath290\\ + e^{-w_{1\text{c}}}\bar n_\text{i}(\omega)\frac{dw_\text{i}}{d\omega } + \underset{\hookrightarrow0}{e^{-w_{1\text{i}}}}\bar n_\text{c}(\omega)\frac{dw_\text{c}}{d\omega},\nonumber\end{aligned}\ ] ] with @xmath291 and @xmath292 the corresponding multiplicity spectra for pure coherent and incoherent components .
in particular , from eqs .
( [ convol - for - n - omega ] ) and ( [ wcoh - times - winc ] ) one observes that both @xmath32 , and @xmath293 involve @xmath50-dependent factors @xmath294 , which cancel between l.h.s . and
r.h.s . of eq .  ( [ convol - for - n - omega ] ) .
thus , @xmath128 is independent of @xmath50 , as anticipated .      for an efficient coherent radiator
, the coherent spectral component must prevail in the fundamental interval @xmath295 , anyway .
it is thus instructive to repeat the analysis of sec .
[ subsubsec : coh - in - fund - interv ] for a mixture of coherent and incoherent radiation components in the fundamental interval .
we will barely sketch it here , emphasizing the distinctions from the pure coherent radiation case .
once expressions ( [ coh+inc ] ) , ( [ dw1inc ] ) , ( [ dw1coh - real ] ) are inserted to eq .
( [ generic - eq - upperlimitomega ] ) , manipulations similar to those used at derivation of eq .
( [ dwc - oint ] ) yield @xmath296 here exponential decrease of the integrand at @xmath159 justifies deformation of the integration contour to shape @xmath147 ( see fig .
[ fig : contours ] ) , but its complete enclosure is impossible because of the presence of branching factor @xmath297 .
[ [ infrared - limit.-1 ] ] infrared limit . + + + + + + + + + + + + + + +    when @xmath298 , contour integral ( [ dw - first - interv ] ) is dominated by large @xmath216 . then , the leading - order term in asymptotic expansion ( [ exp - to - laurent ] ) reproduces result ( [ dw - small - omega ] )
. retaining next - to - leading order terms yields an @xmath223 correction : @xmath299+\mathcal{o}(\omega^2)\bigg\}.\end{aligned}\ ] ] compared to ( [ coh - omegato0-corr ] ) , apart from the overall factor @xmath300 , there emerges a term @xmath301 , whereas other terms remain essentially the same .
the term in the second line of eq .
( [ dw - small - omega - nlo ] ) changes its sign under the same condition ( [ monot - after4 ] ) , so inception of high - intensity regime is essentially independent of @xmath135 .
the behavior of asymptotic approximation ( [ dw - small - omega - nlo ] ) is illustrated in fig .
[ fig : high - intensity]a by the gray long - dashed curve .
[ [ reduction - of - the - fundamental - maximum - by - hard - incoherent - bremsstrahlung . ] ] reduction of the fundamental maximum by hard incoherent bremsstrahlung .
+ + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +    in the point of the spectral fundamental maximum , @xmath302 , integral ( [ dw - first - interv ] ) reduces to @xmath303 since @xmath135 and @xmath192 are both proportional to the target thickness , the increase of the latter at fixed @xmath266 [ cf .
( [ a / bomega0 ] ) ] corresponds to a simultaneous increase of @xmath192 and @xmath135 , at ratio @xmath304 held fixed . for an exemplary ratio @xmath305 , the spectral intensity in the fundamental maximum
is illustrated in fig .
[ fig : intensity - in - max ] by solid gray curve . compared to the pure coherent radiation spectrum
, its maximum is achieved at a slightly lower value of @xmath192 , and its height is somewhat lower , too .
this owes primarily to factor @xmath306 , where @xmath307 is a small number , but there can also be other appreciable effects of @xmath135 .    to assess changes in the height and location of the fundamental maximum due to incoherent radiation , one can , as in sec .
[ subsubsec : coh - in - fund - interv ] , expand the contour integral in ( [ max - dw ] ) to second order in @xmath192 and @xmath135 : @xmath308.\end{aligned}\ ] ] it may suffice to account for corrections in @xmath135 within the leading logarithmic accuracy , which gives
@xmath309\,\,\,\nonumber\\ + \frac{b^2\omega_0}2\int_0 ^ 1dzp(z)p(1-z ) \bigg\}\quad\end{aligned}\ ] ] ( see fig .
[ fig : intensity - in - max ] , dashed black curve ) . here
term @xmath310 can be sizable due to the large logarithm , despite the condition @xmath175 .
owing to its negative sign , that correction reduces the height of the fundamental coherent radiation maximum , in addition to factor @xmath222 .    in the spirit of sec .
[ subsubsec : coh - in - fund - interv ] , for approximation ( [ corr-1fund - max ] ) one can determine the target thickness at which the spectral radiation intensity in the fundamental maximum reaches its highest summit .
it is still given by formula ( [ loc - max ] ) , with parameter @xmath157 unaltered , but parameter @xmath311 modifies to @xmath312 corrections in the numerator and denominator of ( [ eta - coh+incoh ] ) both result in lowering of @xmath311 , and therethrough in an earlier turnover of the fundamental maximum , as confirmed by fig .
[ fig : intensity - in - max ] .
the largest effect , of course , stems from the logarithmic term in ( [ eta - coh+incoh ] ) .      a principally interesting effect of incoherent radiation concerns discontinuities of the coherent radiation spectrum .
as we saw in sec .
[ subsubsec : disc - dwc ] , the discontinuity at @xmath231 is damped by factor @xmath313 . with the addition of the incoherent bremsstrahlung component , this factor actually turns to zero , wherefore the discontinuity must be nullified . ] .
our analysis of the spectrum behavior then needs extension to a non - vanishing vicinity of point @xmath231 .    in vicinity of a singularity of the coherent part of the spectrum , it is convenient to use representation ( [ convol - byparts ] ) .
when @xmath20 slightly exceeds @xmath3 , the integrand of ( [ convol - byparts ] ) receives an extra contribution proportional to a @xmath21-function @xmath314 where @xmath315 is given by eq .
( [ disc ] ) .
integration of the @xmath21-singularity yields a sizable contribution even at an infinitesimal extension of @xmath20 beyond @xmath3 : @xmath316 eq .
( [ smearing - of - disc ] ) shows that for @xmath317 , factor @xmath318 vanishes as @xmath319 , indeed nullifying the discontinuity brought by factor @xmath320 , and making the whole spectrum everywhere continuous .
but the derivative of factor @xmath318 in point @xmath231 diverges , wherewith the resulting spectrum features a sharp spike . the behavior of approximation ( [ smearing - of - disc ] ) is shown in fig .  [ fig : high - intensity]a by the dotted black curve .
it should be realized , though , that the described subtle effect should be obscured in experiments with limited energy resolution and finite angular divergence of the initial beam , which smear the fundamental peak . besides that , the filling of the dip adjacent to the coherent radiation maximum is also partially provided by the pure coherent multiphoton spectrum , whose derivative beyond @xmath3 is negative ( see fig .
[ fig : pure - coh - spectrum]a ) .
thus , verification of threshold behavior ( [ smearing - of - disc ] ) may be feasible only with a rather perfect beam , and in a sufficiently thick target , when @xmath135 becomes sizable .      to accomplish our analysis in different regions , let us consider the large-@xmath20 asymptotics @xmath321 ( while still obeying condition @xmath322 ) .
in this limit , the incoherent bremsstrahlung spectrum decreases by a power law , whereas the multiphoton coherent radiation spectrum , in general , decreases exponentially ( see appendix  [ app : periph ] ) .
therefore , if eq .
( [ wcoh - times - winc ] ) is rewritten as @xmath323 the main contribution to the integral over @xmath324 comes from vicinity of the lower integration limit .
thus , we can expand @xmath325 to taylor series about point @xmath326 , and replace the upper integration limit by infinity : @xmath327 the integral in the second line vanishes due to eq .
( [ mean - omega ] ) , leaving @xmath328.\ ] ] furthermore , in the limit @xmath329 , eq .  ( [ large - omega - corr ] ) boils down to @xmath330,\ ] ] which is shown in figs .
[ fig : high - intensity]a , b by dot - dashed black curves .
( [ simple - large - omega ] ) implies that the combined spectrum ultimately decreases by the same power law as its ( resummed ) incoherent bremsstrahlung component . that property is in contrast with the infrared limit ( [ dw - small - omega ] ) , which yet involves factor @xmath222 of coherent photon non - emission probability [ note the difference between green and black long - dashed curves in fig .
[ fig : high - intensity]a at small @xmath20 ] . finally , it is straightforward to derive large-@xmath20 asymptotics for the photon multiplicity spectrum , with the aid of eq .
( [ convol - for - n - omega ] ) . due to factor @xmath332 in the integrand of ( [ convol - for - n - omega ] ) , again , the dominant contribution comes from the lower integration limit , where it is acceptable to replace @xmath333 $ ] , @xmath334 $ ] , and replace the upper integration limit by infinity .
therewith , eq .
( [ convol - for - n - omega ] ) yields @xmath335,\end{aligned}\ ] ] where @xmath336 is given by eq .
( [ mean - n - omega - incoh ] ) .
( terms involving @xmath222 canceled mutually . )
the integral entering ( [ 83 ] ) equals @xmath337 , as can be checked with the aid of eq .
( [ mean - n - omega - contour - non - casual ] ) . canceling in ( [ 83 ] )
the overall factors according to eq .
( [ simple - large - omega ] ) , we are left with @xmath338\left[1+\mathcal{o}\left(\frac{\overline{\omega_{1\!}}_{\text{c}}}{\omega}\right)\right].\ ] ]    the behavior of the photon multiplicity spectrum is shown in fig .
[ fig : high - intensity]c by the solid green curve . at large @xmath20 , it flattens out , and ultimately enters the regime of logarithmic growth as predicted by eq .
( [ n - omega - asympt - large ] ) ( dot - dashed black curve ) . yet before entering the asymptotic regime , as fig .
[ fig : high - intensity]b shows , the multiplicity spectrum first achieves a maximum , then passes through a shallow minimum .
the theory for such a behavior will be provided in sec .
[ subsec : conv - inc ] .
the preceding analysis revealed that multiphoton emission effects redistribute the coherent radiation spectrum , elevating it beyond the fundamental energy @xmath3 , at the expense of reduction below @xmath3 . at practice ,
though , one should yet be aware of existence of other physical processes which can produce superficially similar effects . in conclusion of this section
, we will briefly discuss such effects , too .
first of all , next to the coherent emission edge @xmath3 , the spectrum may receive an enhancement not only from multiphoton effects , but also from higher harmonics ( [ dw1coh - higher - harm ] ) . at that , higher harmonics produce also 3rd , 4th maxima , etc .
, whereas multiphoton effects can only give rise to a relatively narrow second maximum between @xmath3 and @xmath339 .
so , if no maxima are observed beyond @xmath339 , the simplest way to judge about the origin of the secondary peak would be to refer to the photon multiplicity spectrum : the greater its discontinuity at @xmath340 , the greater the significance of multiphoton effects .
if measurement of @xmath128 is unfeasible , an indirect criterion may be used : the relative contribution of higher harmonics does not vary with the extent of the radiator , whereas that of multiphoton effects does .
therefore , comparing the spectrum shapes at two different radiator lengths , and assessing the length - dependence nonlinearity ( cf .
[ subsec : ln ] ) , one can judge about the origin of the spectrum enhancement beyond @xmath3 .
another effect is the difference between asymptotics of the resummed energy spectrum @xmath245 at @xmath120 and at @xmath341 [ see fig .  [
fig : high - intensity]b , and eqs .
( [ simple - large - omega ] ) and ( [ dw - small - omega ] ) ] .
this may either owe to multiphoton effects , or to incomplete lpm - suppression of radiation in a finite - thickness target [ existing as well in an amorphous matter ( tsf - effect ) , see @xcite and refs . therein ] .
again , unambiguous discrimination between those effects may rely on the photon multiplicity spectrum behavior : if within the suppression region @xmath128 rises significantly above unity ( cf . , e.g. , fig .  [
fig : high - intensity]c ) , the suppression origin ought to be attributed to multiphoton effects , otherwise , it is more likely to be due to lpm - like effects .
discrimination criteria based on the target thickness dependence appear to be nonlinear both for multiphoton and incomplete lpm suppressions , and therefore are to be used cautiously .
the measure of significance of multiphoton effects in a resummed radiation spectrum is given by the mean photon multiplicity , defined by eqs .
( [ mean - n - def ] ) , ( [ w1-def ] ) . for the dominant coherent radiation component ( [ dw1coh - real][f - def ] ) , it estimates as @xmath342 from the practical viewpoint , it is important further to estimate how large this parameter can be for 3 basic coherent gamma - radiation source types : coherent bremsstrahlung , channeling radiation , and undulator radiation .
generic expressions for product @xmath206 for mentioned cases are quoted in appendix [ app : coh - rad - sources ] , so it is now left to assess the entering parameters for conditions of former and future experiments .    for coherent bremsstrahlung , parameter @xmath206
is described by eq .
( [ dw - cb - fmax ] ) .
early coherent bremsstrahlung spectrum measurements operated with relatively thin crystals ( @xmath343 mm ) and relatively large misalignment angles @xmath344 rad , at which @xmath208 was comfortably small . but with the advent of more practical radiation sources having @xmath345 cm and @xmath346 rad , formidable values @xmath347 were reached ( corresponding to an extensive electromagnetic shower ) .    for channeling radiation , the reference equation is ( [ dw - chann - fin ] ) .
it tells that at multi - gev positron energies , and @xmath348 , the photon multiplicity must achieve values @xmath349 .
multiphoton effects in channeling radiation experiments were found appreciable already when dealing with moderately high energies @xmath350 gev , and moderate - thickness targets , @xmath351 , corresponding to @xmath352 @xcite . in more recent cern channeling experiments @xcite , with @xmath353 gev and @xmath354 ,
photon multiplicities reached the order of 5 , and the measured spectra were apparently loosing the features characteristic of coherent radiation ( see also @xcite ) .
it should be realized , though , that at @xmath197 gev , channeling radiation becomes non - dipole , wherewith simple form ( [ dw1coh - real][f - def ] ) for the single - photon spectrum does not apply [ although contour integral representations ( [ generic - eq ] ) , ( [ generic - eq - upperlimitomega ] ) hold generally ] .    for undulator radiation
, we must refer to eq .
( [ dw - undul ] ) . for the newest and forthcoming undulators characterized by @xmath355 and @xmath356 @xcite
, parameter @xmath206 must be @xmath357 .
the only question is whether @xmath3 belongs to the gamma - range , for which the calorimetric method of spectrum measurement is pertinent . in the slac experiment e-166 @xcite
, @xmath3 reached @xmath358 mev , though only under moderate photon multiplicity @xmath359 . for the tesla design ( @xmath360 gev , @xmath361 m , @xmath362 , @xmath363 )
, @xmath3 must rise to 25 mev , and correspond to @xmath364 @xcite .
such photon energies are admittedly well into the gamma domain .
the present estimates show that virtually for all practical coherent radiation sources it appears both feasible and beneficial to enter the regime @xmath349 . in some cases
, it may be optimal to confine to values @xmath365 , when the spectrum is not strongly multiphoton yet ( see fig .
[ fig : intensity - in - max ] ) .
in other cases , such as positron sources @xcite , @xmath208 is demanded to be large indeed . in either case , exploring the asymptotic limit of high @xmath208 is of fundamental interest .
its study will occupy us in the remainder of this section .
to get feel of the trends for the resummed radiation spectrum behavior at high photon multiplicity , let us first take a look at spectral moments introduced in sec .
[ subsec : moments ] .
relations ( [ mean - omega ] ) , ( [ mean - omega2-alt ] ) suggest that with the increase of the radiator length , multiphoton spectrum moments grow proportionally : @xmath366 @xmath367 ( and normally @xmath368 ) .
therewith , the ratio of the width to mean value decreases : @xmath369 that decrease implies that the multiphoton spectrum becomes more sharply peaked .
furthermore , from eq .
( [ mean - omega4 ] ) one infers that the so - called skewness @xcite @xmath370 decreases as well , while the kurtosis @xmath371\ ] ] tends to a constant value 3 ( characteristic of a gaussian ) . ) , ratio @xmath372 appears to be always greater than 3 , and diverges at small @xmath192 , even though the underlying single photon spectrum may well be leptokurtic [ like that defined by eqs .
( [ dw1coh - real][f - def ] ) ] . that owes to our definition of the moments including term @xmath373 in the weighting distribution , which makes any distribution platikurtic . in the high - intensity limit , where @xmath149 exponentially , this definition suits us , anyway . ]
the above observations look natural from the point that the process of statistically independent photon emission is a kind of a random walk @xcite in the energy space . for our case ,
the walk is one - sided , with continuously distributed step size , yet the total probability of each step is less than unity . to substantiate the analogy further , @xmath94 in eq .
( [ generic - eq ] ) may be thought of as depending on the target thickness @xmath100 .
then , differentiating both sides of the equation by @xmath100 yields a kinetic equation    [ kinetic - eq ] @xmath374 where kernel @xmath375 stands for the radiative cross - section .
if we take into account that in the first term of ( [ kinetic - eq1 ] ) @xmath376 equals zero for @xmath97 , we may also rewrite the equation as @xmath377    the initial condition for eq .
( [ kinetic - eq1 ] ) or ( [ kinetic - eq2 ] ) is @xmath378    it is noteworthy that the last term in eqs .
( [ kinetic - eq ] ) makes them inhomogeneous with respect to @xmath32 .
that term is necessary to describe the change of the normalization with @xmath100 according to eq .
( [ total - probab-1 ] ) . at small @xmath100 , when @xmath379 , the inhomogeneous term dominates , but at large @xmath100
, it vanishes exponentially . for the electron distribution function @xmath380 ( see sec .  [ subsec : cont - int - repr ] ) , the kinetic equation will be strictly homogeneous :    [ kinetic - eq - for - pie ] @xmath381 \\
\equiv\int_0^\omega d\omega_1 \pi(e_e-\omega+\omega_1)\frac{\partial}{\partial l}\frac{dw_1}{d\omega_1}\quad \nonumber\\ -\pi(e_e-\omega)\int_0^e d\omega_1 \frac{\partial}{\partial l}\frac{dw_1}{d\omega_1}.\qquad\label{kinetic - eq - for - pie2}\end{aligned}\ ] ]    with initial condition @xmath382 in the literature , the kinetic equation is usually quoted in the form of electromagnetic cascade equation ( [ kinetic - eq - for - pie]b ) , but the function entering thereto is sometimes also called the multiphoton spectrum ( which may be misleading ) . in the presence of significant incoherent radiation , though , the inhomogeneous term should vanish , anyway , insofar as @xmath127 .    note , incidentally , that integro - differential equation ( [ kinetic - eq - for - pie ] ) is well - suited for numerical solution by monte - carlo method , allowing one to simulate a multiphoton radiation spectrum without actually computing separate @xmath2-photon components .
alternative forms ( [ kinetic - eq]a , b ) are suitable for that purpose , too , but granted that their initial condition ( [ kinetic - eq - init - cond-0 ] ) is non - singular , they must be solvable as well by non - random finite - difference methods .
in fact , knowledge of @xmath383 may be sufficient to extract also the photon multiplicity spectrum : assuming that the single - photon spectrum is proportional to @xmath100 , i.e. @xmath384 , eq .  ( [ kinetic - eq2 ] ) rewrites @xmath385 then , noticing that the first two terms on the r.h.s . of ( [ rewrit ] )
coincide with the r.h.s . of eq .
( [ mean - n - omega - convol ] ) , one arrives at a relation @xmath386 ( where for ir regularization , @xmath387 may need to be replaced by @xmath388 ) .
the linear homogeneous 1st - order integro - differential form of equation ( [ kinetic - eq - for - pie ] ) provides the correspondence between our multiphoton emission process and continuous random walks ( a rather general category of markov processes ) in principle may depend on @xmath100 , i.e. , on the electron history , in an arbitrary way , but its further promotion to multiphoton spectrum @xmath389 is sought as solution of a markovian differential equation with respect to increase of @xmath390 , or simply @xmath100 . ] .
( [ kinetic - eq ] ) actually describe a _ driven _ random walk , although it is equivalent to a free one , and belongs to the markovian process category , anyway .
now , let us turn to the case of large @xmath100 .
then , @xmath13 becomes small , and the radiation process definitely enters the free random walk regime .
intuitively , for long random walks , detail of the single - step distribution should fade away after a large number of steps .
the central limit theorem @xcite asserts that for a sufficiently long random walk , the particle probability distribution tends to a gaussian , involving only two parameters : the mean value and the variance of the single - step distribution ( single - photon spectrum ) .
it must be realized , however , that under the presence of incoherent bremsstrahlung component , the first moments defining the gaussian diverge .
that is a familiar situation when the gaussianity breaks down , and the diffusion becomes anomalous @xcite .
for such a case , there exists a generalized central limit theorem @xcite , but it is not really relevant for our conditions , as long as the incoherent component intensity , quantified by parameter @xmath135 , is not supposed to become large ( @xmath391 ) .    a reasonable way out of the present situation may be to treat the incoherent radiation component as a single - jump process , rather than a multistep random walk . to this end , one may first cope with the pure coherent radiation in the high multiplicity regime independently , where the problem of divergent moments is not encountered ( that may also be of intrinsic interest in application to undulator radiation , where @xmath135 is often negligible ) . at the final step of the calculation
, the contribution from the incoherent radiation component can be restored via convolution relation ( [ wcoh - times - winc ] ) .
following this approach , in the next subsection we deal with the pure coherent radiation case .      since typical @xmath20 in the coherent radiation spectrum at high intensity
are @xmath392 , energetically ordered representations offer in this case no advantage .
it is thus simpler to issue from contour integral representation ( [ generic - eq ] ) , @xmath393 where we explicitly introduced the coherent emission edge @xmath3 as the finite upper limit for @xmath35-integration .
obviously , as @xmath264 grows large , the exponential integrand of ( [ cont - int - for - wc ] ) steepens as a function of @xmath53 , peaking somewhere between the integration contour endpoints .
such a situation suggests application of the steepest descent method , which is a common tool in the random walk theory . for our case , though , the corresponding procedure yet involves certain subtleties which will be highlighted below .
the steepest descent method ( see , e.g. , @xcite ) determines the asymptotics of a contour integral by deforming its integration contour in the complex plane so that it passes through a point , to both sides from which the integrand decreases by its absolute value , without significant oscillations . for an analytic integrand
, this must be a regular extremum ( saddle point in the complex plane ) , where the derivative vanishes .
hence the equation for saddle point @xmath394 of the integrand of ( [ cont - int - for - wc ] ) emerges as ) involves no preexponential , the equation for the saddle point is basically unequivocal . ]
@xmath395\bigg|_{s = s_0}\nonumber\\ = \omega-\int_0^{\omega_0 } d\omega_1 \omega_1\frac{dw_{1\text{c}}}{d\omega_1}e^{-s_0\omega_1}=0.\end{aligned}\ ] ]    once solution @xmath396 to eq .
( [ saddle - point - eq ] ) is found , the exponent in eq .
( [ cont - int - for - wc ] ) is further expanded around this point to second order : @xmath397+\frac12 b[s_0(\omega)](s - s_0)^2%\nonumber\\ + \mathcal{o}\left[(s - s_0)^3\right]\!.\quad\end{aligned}\ ] ] to express coefficients @xmath398 and @xmath399 , we write in the exponent of ( [ cont - int - for - wc ] ) @xmath400,\end{aligned}\ ] ] and notice that term @xmath401 , when convolved with @xmath264 , cancels with term @xmath402 by virtue of eq .
( [ saddle - point - eq ] ) . that produces the expected structure ( [ hd ] ) , wherefrom we read off @xmath403=\int_0^{\omega_0 } d\omega_1 \frac{dw_{1\text{c}}}{d\omega_1}\left[e^{-s_0\omega_1}(1+s_0\omega_1)-1\right],\ ] ] and @xmath404=\int_0^{\omega_0 } d\omega_1 \omega_1 ^ 2\frac{dw_{1\text{c}}}{d\omega_1}e^{-s_0\omega_1}=-\frac1{s_0}\frac{\partial a}{\partial s_0}.\ ] ] the neglect in the exponent of ( [ cont - int - for - wc ] ) of taylor terms higher than quadratic reduces the integrand to a gaussian , and the contour integral evaluates    @xmath405
one should yet be aware that the saddle - point equation in the complex plane may have several solutions .
for our edgy single - photon spectrum ( [ dw1coh - real ] ) , the number of solutions is actually infinite . to demonstrate that
, the distribution of saddle points in the complex @xmath53-plane , and their motion with the variation of @xmath20 , is portrayed in fig .
[ fig : re - plot3d ] for specific single - photon spectrum profile ( [ dw1coh - real][f - def ] ) . there ,
black solid curves correspond to the @xmath20-independent imaginary part of eq .
( [ saddle - point - eq ] ) , giving saddle points trajectories .
dashed color curves correspond to the real part of eq .
( [ saddle - point - eq ] ) ; their intersections with the solid curves give saddle point locations , which depend on @xmath20 .
the figure suggests that for any @xmath20 , the number of saddle points must be infinite , so instead of eq .
( [ dwdomega - saddle ] ) one ougth to write @xmath406    where the summation runs over all the saddle points , and @xmath407 , @xmath408 .
( a ) the distribution of saddle points in the complex plane of @xmath53 .
solid black curves are solutions of equation @xmath409 ( saddle - point trajectories ) .
dashed curves are solutions of equation @xmath410 ( magenta , @xmath411 ; green , @xmath412 ; cyan , @xmath413 )
. intersections of solid curves with dashed ones give saddle points at a given @xmath414 . the density plot , function @xmath415 , for @xmath416 given by eq .
( [ h ] ) ( brighter regions correspond to greater values of @xmath417 ) .
( b ) values of function @xmath415 in the saddle points .
[ for @xmath413 ( cyan ) , the plotted function is multiplied by 0.02 . ] ]    nevertheless , in the limit @xmath418 , series ( [ dwdomega - saddle - sum ] ) must be dominated by a term with the greatest @xmath417 . the diagram in fig .
[ fig : re - plot3d]b shows the values of @xmath419 in saddle points . for any @xmath20 , the saddle point at @xmath420 appears to be the highest , and therefore asymptotically dominant .
thus , in the limit @xmath421 , eq .
( [ dwdomega - saddle ] ) is recovered , with @xmath422 being the unique _ real _ solution of eq .
( [ saddle - point - eq ] ) .
another issue is the robustness of the saddle point approximation itself .
that approximation breaks down in vicinity of the point where coefficient @xmath399 vanishes , i.e. at small @xmath20 corresponding to large @xmath394 . to determine the condition of applicability of the saddle - point approximation , one needs to assess higher terms in @xmath423 in eq .
( [ hd ] ) .
that leads to the inequality @xmath424 the l.h.s .
thereof admits interpretation as an order of magnitude of the number of emitted photons with energies below @xmath20 .
thus , ineq .
( [ saddle - point - cond ] ) implies that for the diffusion mode to develop , the same ` destination ' ( total energy ) @xmath20 must be reachable by many radiative scenarios with commensurable and sizable probability .
if all the premises listed above are fulfilled , eq .  ( [ dwdomega - saddle ] ) offers an explicit solution for the resummed spectrum , provided the real solution to equation ( [ saddle - point - eq ] ) is found and substituted to eqs .
( [ h][d ] ) .
an impediment is that even with simplest shapes of @xmath264 , like ( [ f - def ] ) or a uniform distribution , equation ( [ saddle - point - eq ] ) is transcendental . on the other hand , in the crudest approximation
, it may suffice to know its behavior around single point corresponding to the maximum of expression ( [ dwdomega - saddle ] ) .
the latter region , called central , will be scrutinized in the next paragraph .
the problem of finding the maximum of distribution ( [ dwdomega - saddle ] ) is essentially equivalent to tracing that for exponent @xmath398 , as long as the exponential should be an asymptotically steeper function than its preexponential factor .
inspection of integral ( [ h ] ) reveals that for any positive kernel @xmath264 , the integral achieves a maximum in point @xmath425 , since there the @xmath394-dependent factor @xmath426 $ ] turns to zero , while being negative elsewhere . from eq .
( [ saddle - point - eq ] ) , one finds that value @xmath425 corresponds to @xmath427 .
expansion the of r.h.s . of eq .
( [ h ] ) by deviations of @xmath394 from zero ( beginning from a quadratic term ) yields @xmath428= -\frac12 \overline{\omega^2_{1\!}}_{\text{c } } s_0 ^ 2+\mathcal{o}\left(\overline{\omega^3_{1\!}}_{\text{c}}s_0 ^ 3\right).\ ] ] if eq .
( [ saddle - point - eq ] ) in turn is linearized in @xmath394 , @xmath429 it solves as @xmath430 combining eqs .
( [ hcoh - quadr ] ) and ( [ saddle - point - eq - lin - soln ] ) , we obtain the explicit asymptotic expression for the probability spectrum in vicinity of the maximum ( central region ) : @xmath431 with @xmath432 as anticipated , this limiting distribution is merely a gaussian , with the mean value and variance complying with eqs .
( [ mean - omega][mean - omega2-alt ] ) , for an arbitrary shape of @xmath271 ( see , e.g. , @xcite ) . since dependencies @xmath433 for different shapes and magnitudes of @xmath264 and intensities are related by a linear change of the variables ( i.e.
, the multiphoton spectrum shape is invariant under changes of @xmath264 ) , the central region is also called scaling one @xcite . so far , the approach to gaussian asymptotics remains largely unexplored experimentally , although spectra resembling gaussians were obtained , e.g. , in .
approximation ( [ gauss ] ) holds as long as in eq .
( [ saddle - point - eq - lin ] ) nonlinear terms in @xmath394 keep relatively small , i.e. @xmath434 , the gaussian approximation must work in a sufficiently broad interval about its maximum .
the corresponding photon multiplicity spectrum @xmath435 is further evaluated by referring to representation ( [ mean - n - omega - convol ] ) .
its second term vanishes due to practically vanishing @xmath13 , or due to vanishing @xmath264 for @xmath436 , while in the first term , the integration is effectively limited to interval @xmath437 , where factor @xmath264 differs from zero . noting that within such an interval , the second factor @xmath438 varies slowly
, one can taylor - expand it to first order , say , about the origin : @xmath439 wherewith @xmath440 inserting here expression ( [ gauss ] ) , we obtain an utterly simple result : @xmath441 the linear function in the r.h.s
. of ( [ mean - n - at - gauss ] ) has positive intercept @xmath442 , and slope , but at small @xmath20 , its slope is @xmath443 [ see eq .
( [ n0=1 + ] ) ] . ]
@xmath444 ; the latter implies that @xmath257 increases by more than unity per each interval @xmath445 .
it is also evident that in vicinity of the gaussian maximum , @xmath446 which , according to ( [ mean - n - def ] ) , amounts the mean photon multiplicity for the whole coherent radiation process .
an approximately linear increase of @xmath257 in the range of typical @xmath20 was observed in experiment @xcite .
there was also observed a turnover of @xmath257 , one of the possible reasons for which will be explained later .    in fig .
[ fig : pure - coh - spectrum ] , the exact resummed spectrum shape for radiation intensity parameter @xmath202 ( blue solid curve ) is compared with the behavior of gaussian approximation ( [ gauss ] ) , ( [ mean - n - at - gauss ] ) ( blue dashed curve ) . at the chosen value of @xmath206 , a perfect convergence for the resummed spectrum ( fig .
[ fig : pure - coh - spectrum]a ) is not reached yet . at the same time , for @xmath257 ( fig .
[ fig : pure - coh - spectrum]c ) , the agreement is already excellent . remarkably , the linearity of the latter function holds well even far away from the central region , where the spectrum is no longer gaussian [ @xmath422 dependence gets highly nonlinear , and variation of the preexponent @xmath399 becomes significant , as well ] . to improve the accuracy of the diffusive approximation ,
one may be prompted to derive a correction to eq .
( [ gauss ] ) . to this end , note that eq .
( [ gauss ] ) corresponds to expanding the integral in the exponent in eq .
( [ cont - int - for - wc ] ) to second order in @xmath53 .
a natural improvement thus results if in the expansion of the exponentiated integral a cubic term in @xmath53 is retained : @xmath447 the latter contour integral is reducible to an airy function , but exponentiation of a small correction may be an excess of precision .
if we simply linearize the exponential , @xmath448 and evaluate the contour integral , it leads to an additive correction dating back to chebyshev ( see , e.g. , @xcite ) : @xmath449,\ ] ] with @xmath450 defined by eq .
( [ skewness ] ) , and @xmath451 the hermite polynomial of order 3 @xcite . the behavior of approximation ( [ chebyshev ] )
is illustrated in fig .
[ fig : pure - coh - spectrum ] and fig .
[ fig : logplot ] ( for a greater @xmath206 ) by dot - dashed curves , and is appreciably closer to the exact distribution than a simple gaussian ( dashed curves ) .
the evaluated skewness correction mildly breaks the scaling , and serves to describe the residual asymmetry of the spectrum .
one of its consequences is the redshift of the spectral maximum location to the value @xmath452 where the correction term is independent of the radiation intensity [ for the spectral shape ( [ f - def ] ) , @xmath453 .
however , the correction to the height of the maximum in this approximation is zero : @xmath454.\ ] ]    fig .
[ fig : logplot]a also indicates that even though the gaussian approximation , especially with skewness correction ( [ chebyshev ] ) , works fairly well , away from the central region the falloff law differs from the gaussian , and the scaling property gets strongly violated . in principle , this deviation can still be calculated based on the steepest descent method of sec .
[ subsubsec : steepest - descent ] , provided one copes with the saddle point equation at @xmath455 .
approximations derived along these lines for sub- and trans - central regions are described in appendix  [ app : periph ] .
let us finally inspect how does the multiphoton spectrum shape modify upon incorporation of an incoherent component . as in sec .
[ subsec : coh+incoh ] , this is accomplished by convolving the evaluated pure coherent radiation contribution with pure incoherent component ( [ 13 ] ) via eq .
( [ wcoh - times - winc ] ) . given the concentration of the coherent radiation probability in the central region
, it may be reasonable to employ in ( [ wcoh - times - winc ] ) the simplest gaussian form ( [ gauss ] ) . that leaves us with the expression @xmath456 ( term @xmath457 has been neglected as exponentially small at large @xmath208 . )     ( a ) log plot of the resummed pure coherent radiation spectrum , at intensity parameter @xmath458 .
solid purple curve , exact distribution evaluated by eqs .
( [ cont - int - for - wc ] ) , ( [ dw1coh - real][f - def ] ) .
dashed black curve , gaussian approximation ( [ gauss ] ) .
dot - dashed black curve , gaussian approximation with chebyshev correction ( [ chebyshev ] ) .
( b ) photon multiplicity spectrum under same conditions .
dashed black line , approximation ( [ mean - n - at - gauss ] ) . ]
the encountered integral still depends on rather many parameters , but under the high - intensity condition @xmath459 , their number can be effectively reduced .
the upper limit of the integral is always on the right of the maximum of the gaussian entering the integrand , hence at the upper integration limit , the gaussian is decreasing .
granted that this decrease is rapid , one can replace this integration limit by infinity , unless the integration interval is too short .
more precisely , the gaussian varies across the integration interval significantly if @xmath460 i.e. , @xmath461 condition ( [ small - omega - range ] ) is not very restrictive as long as its r.h.s . is independent of the radiation intensity , while its l.h.s .
is typically @xmath462 .
thus , in the high photon multiplicity limit , it appears legitimate virtually at all @xmath20 to replace the upper integration limit in eq .
( [ 1st - conv ] ) by infinity .
therewith , utilizing the value of the integral @xmath463 defining parabolic cylinder function @xmath464 @xcite , we arrive at the closed - form representation @xmath465 since function ( [ parab - cyl ] ) , similarly to ( [ gauss ] ) , depends on @xmath20 only through a linearly related variable @xmath466 , this may be regarded as extension of the scaling property beyond the vicinity of the maximum .
but it should be remembered that the spectrum shape also depends on parameter @xmath135 , so here we have a 2-parameter scaling .
the corresponding functions , though , differ from lvy distributions , rather being intermediate between lvy distributions and gaussian distributions .
the accuracy of approximation ( [ parab - cyl ] ) may be assessed from fig .
[ fig : loglogplot - a]a .
there one observes that the parabolic cylinder approximation ( dashed black curve ) holds fairly well almost everywhere , except the small-@xmath20 limit .
based on the obtained compact formula ( [ parab - cyl ] ) , we may now assess the deviation from the gaussian behavior in the central region . in our equations ,
parameter @xmath135 must always remain small , and it is straightforward to check from definition ( [ d - def ] ) that @xmath467 , wherefore , at @xmath280 distribution ( [ parab - cyl ] ) returns to gaussian form ( [ gauss ] ) . in the next approximation ,
the maximum of function ( [ d - def ] ) satisfies the equation @xmath468 at small @xmath135 , the solution of eq .
( [ eq - for - zmax ] ) is small , too , and in the linear approximation it equals @xmath469 in terms of @xmath20 , it casts as @xmath470 one may as well add here the skewness correction ( [ omega - max - skew ] ) ( provided both corrections are small ) , getting : @xmath471 note that those corrections have opposite signs , i.e. , the incoherent radiation component gives a blueshift , while skewness a redshift , and they partially compensate each other . evaluating the second derivative in the maximum , one can also show that the effect of @xmath135 is to broaden the spectrum .
( a ) resummed radiation spectrum for coherent component described by eqs .  ( [ dw1coh - real][f - def ] ) with @xmath458 , and incoherent radiation component with @xmath472 , @xmath473 . solid purple curve , exact distribution ; dashed black curve , parabolic cylinder approximation ( [ parab - cyl ] ) .
( b ) the photon multiplicity spectrum for same conditions as ( a ) , and @xmath474 . dashed black curve , parabolic cylinder approximation ( [ convol - for - n - omega2 ] ) .
dot - dashed black curve , high-@xmath20 asymptotics ( [ n - omega - asympt - large ] ) . ]    at large @xmath20 , the spectrum deviates from a gaussian strongly , developing a power - law ` tail ' .
the parabolic cylinder function at large values of its argument has asymptotics @xmath475+\mathcal{o}\!\left(\!e^{-\rho^2\!/2}\right)\!,\ ] ] entailing for the spectrum @xmath476 that agrees with eq .
( [ large - omega - corr ] ) modulo exponentially small terms .
the calculation of @xmath128 proceeds along the same lines when based on convolution relation ( [ convol - for - n - omega ] ) .
employing there gaussian approximations ( [ gauss ] ) , ( [ mean - n - at - gauss ] ) along with power - law and logarithmic expressions ( [ 13 ] ) , ( [ mean - n - omega - incoh ] ) for the pure incoherent component , we are led to the integral representation @xmath477.\quad\end{aligned}\ ] ] by the same reasoning as for eq .
( [ 1st - conv ] ) , the upper integration limit here may be extended to infinity .
evaluation of the integral gives the result @xmath478 where @xmath479 and @xmath480     ( a ) plots of function @xmath481 , for parameter @xmath482 , 0.03 , 0.1 , 0.3 ( top to bottom ) .
( b ) plots of @xmath483 for same values of @xmath135 . ]
shapes of functions ( [ r - def ] ) , ( [ s - def ] ) for several values of parameter @xmath135 are illustrated in figs .
[ fig : ra , sa]a , b .
function @xmath484 appears to be monotonous , though step - like , whereas @xmath485 features a hump .
the maximum of @xmath485 rises with the decrease of parameter @xmath135 , though only logarithmically .
note that function @xmath486 enters eq .
( [ convol - for - n - omega2 ] ) with a large coefficient , thence the shape of the maximum of @xmath128 is primarily described by @xmath485 .
an additive correction comes from @xmath484 , whose unit jump in the turnover region of @xmath485 partially compensates the drop of the photon multiplicity spectrum beyond its maximum .
the analysis of behavior of functions ( [ r - def ] ) , ( [ s - def ] ) is alleviated by the knowledge of their values in the origin , @xmath487 @xmath488,\ ] ] and asymptotics at large @xmath489 : @xmath490 @xmath491 @xmath492 @xmath493 with @xmath174 defined by eq .
( [ nu - def ] ) . at @xmath494 , i.e. asymptotic trans - central @xmath20 ,
function @xmath484 dominates , and along with terms @xmath495 reproduces result ( [ n - omega - asympt - large ] ) .    finally , we recall that the humpy structure of @xmath128 manifested itself already at moderate radiation intensity ( cf . fig .  [
fig : high - intensity]c ) . that may be non - accidental , since according to ineq .
( [ saddle - point - cond ] ) , the accuracy of the saddle - point approximation improves at large @xmath20 , where the maximum of @xmath128 is achieved .
although measurements of radiation spectra by the calorimetric method are ubiquitous , for further applications of the extracted photon beam it is often preferable to describe it in terms of photon number spectrum ( [ class ] ) , anyway .
needless to say , the rate of any reaction ( such as @xmath0 creation or nuclear reactions ) caused by a photon beam in a target is proportional to the photon number per energy interval , irrespective of temporal correlations between photons . besides that
, if the radiation spectrum is used for diagnostics of the radiator , the single - photon spectrum is more directly related to the classical radiation spectrum and the charged particle trajectory , and hence to the radiator parameters , than the multiphoton one .
thus , there arises a need for reconstructing single - photon spectra from calorimetrically measured spectra .
such a problem had been addressed some time ago in @xcite , where a procedure for reconstructing first moments ( the mean and the variance ) of the single - photon spectra from those of a multiphoton spectra was proposed . here
we will develop a reconstruction procedure for the complete spectrum .
the prerequisite , however , is fair independence of the spectrum on initial state parameters to be averaged over , wherewith the outcoming multiphoton spectrum is representable purely by contour integral ( [ generic - eq ] ) . as was argued in sec .
[ subsec : pure - coh ] , such conditions can be met for coherent bremsstrahlung or undulator radiation , but are unlikely for channeling radiation .
let us begin with the observation that integral relation ( [ generic - exponentiation ] ) , with the aid of eq .
( [ total - probab-1 ] ) can be cast in form where the l.h.s .
depends only on @xmath31 , while the r.h.s . only on @xmath32 : @xmath496 ( as long as @xmath32 is supposed to be inferred from experimental data , the upper limit of integration over @xmath20 is set to @xmath4 , which is the physical end of the spectrum . ) to further obtain a linear integral equation for the single - photon spectrum , it suffices to take the logarithm of both sides of eq .
( [ 90 ] ) :    @xmath497.\end{aligned}\ ] ]    after integration in the l.h.s . by parts
, this can also be written @xmath498.\,\,\end{aligned}\ ] ]    if @xmath11 is ir - finite , one can deal with eq .
( [ simpler ] ) arranged with the aid of relation ( [ total - probab-1 ] ) in a simple form @xmath499 since here the l.h.s .
is merely a laplace transform of @xmath31 , it inverts in the standard way :    [ 110 ] @xmath500 @xmath501 now , for @xmath502 ( when the @xmath21-term does not contribute ) , ir divergences do not appear in eq .
( [ reconstr ] ) , so the initial assumption of ir finiteness of @xmath11 may be relaxed .
( [ 110 ] ) might as well be derived from eq .
( [ finite ] ) , which is a bit more lengthy ( solving a 1st - order differential equation ) , but does not need presuming ir finiteness of @xmath11 . if @xmath13 is finite , integral ( [ log(1 + ... ) ] ) converges better , since the logarithm in the integrand vanishes at infinite integration limits . on the contrary , if @xmath503 , representation ( [ log(1 + ... ) ] ) is untenable , and ( [ reconstr ] ) ( without the @xmath21-term ) is the representation of choice .    from our experience with the resummation procedure in sec .
[ sec:2 ] , one may expect ambiguity in the contour integral representation for the reconstructed spectrum , as well .
indeed , applying lemma  [ lemma1 ] with @xmath504 , or @xmath505 , , the argument of the logarithm , @xmath506 , tends to zero as @xmath507 .
thereat , @xmath91|=\frac1{|z(s)|}$ ] increases , but only by a power law , @xmath508 , still satisfying the conditions of lemma  [ lemma1 ] . ]
one can see that the upper integration limit for the laplace transform of @xmath32 in ( [ 110 ] ) may be an arbitrary number greater than @xmath35 , e.g. equal @xmath509 : @xmath510    relation @xmath97 in integral ( [ reconstr - caus ] ) may seem to be physically counterintuitive , because it is energetically _ anti_-ordered .
but as we emphasized in sec .
[ subsec : cont - int - repr ] , the energetic ordering property does not hold exactly for our distributions , anyway .    another way to arrive at eq .
( [ reconstr - caus ] ) is to expand the logarithm in eq .
( [ log(1 + ... ) ] ) to power series , and evaluate the contour integral termwise :    [ reconstr - power - series - together ] @xmath511 here , as in passing from eq .
( [ generic - sum - int ] ) to ( [ generic - sum - int - upper - omega ] ) , one may replace upper limits @xmath4 by @xmath512 : @xmath513    in contrast to eq .
( [ generic - sum - int ] ) , however , series ( [ reconstr - power - series - together ] ) involves no factorial in the denominator , whence it must have a finite convergence radius , such that @xmath514 .
the latter property is natural , as long as resummed spectra are always bounded from above ( see sec .
[ par : fund - max ] , [ subsubsec : comb - in - fund - int ] ) .
series ( [ reconstr - power - series - together ] ) at sufficiently small intensities of @xmath11 may prove to be actually even more convenient than contour integral representations ( [ 110 ] ) .
in particular , in the limit @xmath47 , only the first term of ( [ reconstr - power - series - caus ] ) survives , giving @xmath515 which is equivalent to eq .
( [ coh - omegato0 ] ) . keeping @xmath516 terms
, one can derive an analog of eq .
( [ coh - omegato0-corr ] ) , etc .    for application of contour integral representations ( [ 110 ] ) , it is also highly desirable to map singularities . ] , as well as sectors of growth and decrease of the integrand , to which the integration contour can be adjusted .
as long as the radiation spectrum obeys eq .
( [ 108 ] ) , the integrand singularities in eqs .
( [ 110 ] ) must coincide with those of laplace transform of the single - photon spectrum .
but the latter transform may not have singularities in the right half - plane , since @xmath517 converges at any positive @xmath86 , so it is an analytic function of @xmath53 at @xmath518 . in case
if with the increase of @xmath35 , @xmath31 decreases faster than exponentially [ e.g. , strictly vanishes beyond a coherent emission edge @xmath3 , like in eq .
( [ dw1coh - real ] ) ] , the l.h.s . of ( [ 108 ] ) exists at any @xmath53 , so no singularities in the @xmath53 plane can emerge at all . in case
if @xmath31 decreases with @xmath35 exponentially , its laplace transform possesses singularities at finite negative @xmath86 . finally , if @xmath31 decreases by a power law , the l.h.s .
of ( [ 108 ] ) diverges at @xmath519 , so a singularity is situated at point @xmath520 .
as concerns the convergence sectors , the integrand of ( [ log(1 + ... ) ] ) obviously diverges at @xmath521 because of factor @xmath522 . at @xmath523 ,
the integrand may decrease only provided the argument of the logarithm , after analytic continuation to the negative @xmath86 half - plane , does not increase faster than exponentially . in the latter case ,
the integration contour may be deformed to shape @xmath147 ( see fig .
[ fig : contours ] ) .
but that appears to be impossible in the important case when @xmath31 is anticipated to decrease at @xmath524 faster than exponentially ( e.g. , pure single - harmonic coherent radiation ) .
then , as we noted above , the integrand of ( [ 110 ] ) may not have singularities in the whole complex plane , so if the integration contour might be deformed to @xmath147 , the result of the integration would identically equal zero , which is a clear contradiction . in this case , the only option for the integration contour in eqs .
( [ 110 ] ) is a line parallel to the imaginary axis ( @xmath146 ) .    at practice , when @xmath32 is measured as a number of events per bin of deposited energy , in eq .
( [ reconstr ] ) the integral involving @xmath32 must be replaced by the corresponding finite sum .
that is harmless for evaluation of the contour integrals , since the resulting function of @xmath53 remains analytic in the @xmath53 plane .
but it is mandatory that the spectrum measurements are absolute , providing correct normalization of the spectrum .    at high intensity , when any resummed coherent spectrum tends to gaussian form
, our reconstruction method will ultimately become unsustainable . to appreciate the encountered difficulties ,
consider a toy model for the single - photon spectrum : @xmath525 ( with @xmath526 ) .
the corresponding resummed spectrum evaluates in terms of a modified bessel function : @xmath527 laplace transform thereof equals @xmath528 now suppose that a calorimetric spectrum measurement gave instead of @xmath32 some slightly deviating distribution @xmath529 . inserting this to eq .
( [ log(1 + ... ) ] ) leads to the reconstructed underlying spectrum @xmath530 clearly , as @xmath192 increases , signal term @xmath531 becomes exponentially small , and ultimately inferior to the error term ; hence , for sufficiently large @xmath532 the method will loose efficiency . the best one can do in this situation , probably , is to extract first spectral moments as in ref .
@xcite , which may appear insufficient for complete characterization of a highly structured single - photon spectrum .    in the latter straitened circumstances
, it may also be worth engaging information contained in the photon multiplicity spectrum @xmath128 .
examination of structure of integral ( [ mean - n - omega - contour - non - casual ] ) shows that it is more complicated than ( [ generic - eq ] ) , and does not allow expressing @xmath94 explicitly .
but in the high - multiplicity limit , when the reconstruction method based on @xmath258 is obstructed , one can at least extract from @xmath128 an additional parameter @xmath208 , say , through eq .
( [ n - to - wi1c ] )
. its knowledge , along with the knowledge of moments @xmath533 and @xmath534 ( and possibly @xmath535 ) extracted from @xmath32 , may constrain the shape of @xmath94 more tightly .
the techniques developed in the present paper can be further employed to calculate spectra of specific coherent radiation sources , including channeling radiation . for undulator radiation , and partially for coherent bremsstrahlung , the results presented in this paper may already be applicable directly .
let us summarize the main physical lessons learned at this stage :    * the fundamental maximum of the multiphoton spectrum does not rise indefinitely with the increase of the radiator length @xmath100 , but eventually saturates , and subsequently decreases . instead , there is an elevation of the spectrum beyond the fundamental maximum , and its general spread ( sec .
[ sec : coh+incoh ] ) .
the high-@xmath20 spectral tail due to incoherent bremsstrahlung is the least affected by multiphoton effects , and its ratio w.r.t .
the fundamental maximum thereby increases .
* multiphoton effects alone can give rise to a second maximum in the coherent radiation spectrum , at @xmath536 , even when secondary harmonics in the single - photon spectrum are absent .
* the low-@xmath20 part of the multiphoton spectrum , as well as discontinuities of the spectrum for coherent radiation , are suppressed by the factor of photon non - emission probability ( sec .
[ sec : coh+incoh ] ) .
therein , manifestations of the incoherent radiation component in the spectrum are particularly pronounced .
* at high radiation intensity ( in a long radiator ) , the resummed spectrum of pure coheernt radiation tends to a gaussian form .
but under presence of an incoherent component , it features a power - law tail towards large @xmath20 ( sec .  [ sec : high - intensity ] ) .
such a diffusion regime may be regarded as _ weakly _ anomalous ( characterized by @xmath537 ) .
the corresponding limiting multiphoton spectrum is described by a parabolic cylinder function , being an intermediate asymptotic case between gaussian and lvy distributions . *
the photon multiplicity spectrum @xmath128 is an informative observable complementary to @xmath32 , and is measurable with the modern state of the technology .
the growth of @xmath128 with @xmath20 at low radiation intensity is stepwise , whereas at high intensity it smoothens out and tends to a linear law . in the region where the incoherent contribution overtakes the coherent one , @xmath128 saturates , and slightly subsides
, whereupon continues rising , but only logarithmically ( secs .
[ sec : coh+incoh ] , [ sec : high - intensity ] ) .
promising are also the formulas for reconstruction of the single - photon spectrum from the calorimetrically measured one , derived in sec .
[ sec : reconstr ] .
the applications thereof will be considered elsewhere .
it should also be remembered that although the present framework is well suited for fast evaluation of multiphoton emission effects , it is limited to conditions of soft radiation , which do not yet cover all the variety of interesting physical problems ( e.g. , hard coherent electromagnetic cascades , with deposited energies close to @xmath4 , investigated in experiments @xcite ) . nonetheless , for practical coherent radiation sources , as discussed in appendix  [ app : coh - rad - sources ] and the beginning of sec .
[ sec : high - intensity ] , conditions of soft radiation are quite typical . at the same time , albeit many experiments in the soft radiation domain found multiphoton effects to be appreciable , none of those experiments addressed multiphoton effects systematically .
the need for such dedicated experimental studies thus remains acute .
in sec .  [ subsec : incoh - bremsstr ] , we introduced parameter @xmath142 as an effective uv cutoff for the single - photon scale - invariant incoherent bremsstrahlung spectrum .
the value of this parameter can be inferred from a theory accurately describing acts of large radiative energy loss . for our purposes , we need only the nll soft asymptotics of this theory . as long as we deal solely with spectrum @xmath32 ( and are not concerned with the photon multiplicity spectrum ) , the coherent radiation component at large @xmath20 can be entirely neglected .
thereby , we return to the well - known electromagnetic cascade theory in an amorphous medium , but with the due account of the electron spin .    for the single - photon bremsstrahlung spectrum , instead of ( [ dw1incoh - cutoff ] )
we take the ultra - relativistic and completely screened approximation of the bethe - heitler formula : @xmath538 where @xcite @xmath539 for simplicity , the pair production process will be neglected , as before .
then , the kinetic equation for electron energy distribution @xmath540 remains one - component : @xmath541 where dimensionless parameter @xmath135 is proportional to the target thickness [ see eq .
( [ a - amorph ] ) ] . integro - differential equation ( [ shower - eq ] ) is endowed with the initial condition ( [ kinetic - eq - init - cond - for - pie ] ) : @xmath542    once the solution of eq .
( [ shower - eq][shower - init - cond ] ) is attained , the multiphoton radiation spectrum , equal to the radiative energy loss spectrum , is determined by the correspondence rule ( [ def - pie ] ) : @xmath543,\ ] ] and for this function we are interested in the limit @xmath544 .
note that by a change of integration variables , eq .
( [ shower - eq ] ) may can be recast @xmath545 which turns to recoilless kinetic equation ( [ kinetic - eq2 ] ) in the limit @xmath546 , when @xmath547
. however , form ( [ shower - eq ] ) is more convenient for solution .
solution of eq .
( [ shower - eq ] ) is obtained by applying mellin transform @xmath548 in terms of which the integro - differential equation becomes ordinary differential : @xmath549 with @xmath135-independent coefficient @xmath550 and initial condition @xmath551 integrating eq .
( [ shower - eq - mell ] ) , and applying the inverse mellin transform @xmath552 , leads to the solution @xmath553 by correspondence rule ( [ rel - with - multiphot - spectr ] ) , the multiphoton bremsstrahlung spectrum ensues @xmath554 in the last line , we linearized @xmath555 for small @xmath556 we are concerned with , and changed the integration variable to @xmath557 .
function @xmath558 at large arguments needed for us increases logarithmically : @xmath559\qquad\qquad\nonumber\\ = \ln s+\gamma_{\text{e}}+\!\int_0 ^ 1\!\frac{dz'}{z'}[w_{\text{bh}}(z')\!-\!w_{\text{bh}}(0)]+\mathcal{o}\!\left(s^{-1}\right)\!.\nonumber\\\end{aligned}\ ] ] inserting this approximation to eq .
( [ pi - solution ] ) , and doing the elementary contour integral , from comparison with eq .
( [ 13 ] , [ fa ] ) we ultimately infer the value for @xmath44 : @xmath560}=e^{-5/8}\approx0.5.\ ] ] according to eq .
( [ kappa58 ] ) , the effective uv cutoff must be set at about half the initial electron energy . with this value for @xmath44 , our eqs .
( [ 13 ] , [ fa ] ) agree with eqs .
( 2.11 ) , ( 2.19 ) of @xcite .
in this appendix , we will collect equations for spectral parameters of main coherent radiation sources , needed for numerical estimates of significance of multiphoton effects in sec .
[ sec : high - intensity ] .    any coherent radiation source we consider may be viewed as an ultrarelativistic electron ( or positron ) radiating under the influence of periodic force @xmath561 ( having period @xmath562 ) , directed predominantly transverse to the electron s large mean momentum ( which is taken parallel to @xmath563 ) . if the radiation is of dipole type , i.e. produced by a weakly accelerated electron remaining non - relativistic in its average rest frame ( e.r.f . ) , the emitted electromagnetic radiation in this frame will be monochromatic . in e.r.f . ,
period @xmath562 is lorentz - contracted to @xmath564 ( @xmath1 stands for the electron lorentz factor ) , corresponding to frequency , where frequency and mass have dimension of energy , while the electron charge is dimensionless . ]
the radiation in the e.r.f . is not strictly isotropic , as long as the driving force is known to be transverse w.r.t .
@xmath563 . for our purposes , we need only the angular distribution averaged over the azimuthal angles , which will thus obey a @xmath566 law , where @xmath567 is the angle between @xmath563 and the photon emission direction in e.r.f .
this information suffices to write the angle - differential radiation intensity : @xmath568 where @xmath569 and @xmath570 are the electron mass and charge .
[ to check the numerical coefficient , note that integration of ( [ i - dipole ] ) over photon emission angles with @xmath571 , and identification of @xmath572 with the electron acceleration , leads to the famous larmor formula . ]
( [ i - dipole ] ) gives an average ( over the oscillation period ) energy emitted by the electron in its rest frame per unit time .
dividing this by the photon energy @xmath573 gives the probability density , and the total probability is obtained by multiplying by the period value and the number of periods @xmath574 : @xmath575 the latter quantity must be lorentz invariant , which is opportune for expressing the desired radiation spectrum in the lab .    the advantage of referring to e.r.f .
variables is that @xmath576 is related to the photon momentum projection on the light front in @xmath76-direction , and therethrough to the photon energy in the lab frame : @xmath577 where @xmath3 is the maximal photon energy in the lab , corresponding to @xmath578 . expressing @xmath576 from ( [ omega - theta - corresp ] ) , inserting it to ( [ i - dipole ] ) , and differentiating , we convert the angular distribution in e.r.f . to the spectral distribution in the lab : @xmath579 with @xmath198 being function ( [ f - def ] ) , thus explaining its origin .
product @xmath580 entering the time integral is lorentz invariant , too , so it can be written as @xmath581 , where @xmath582 is the force , and @xmath583 the time in the lab . comparing ( [ 152 ] ) with eq .
( [ dw1coh - real ] ) , we infer @xmath584 where @xmath562 is the motion period in the lab , and obviously , @xmath585 .    although product @xmath206 is all we need for estimates of the photon multiplicity , it is also expedient to evaluate the photon energy scale @xmath3 independently .
the latter is obtained from e.r.f .
value @xmath586 by boosting along @xmath563 with the lorentz factor @xmath1 .
that results in multiplication by @xmath587 , so the maximal photon energy in the lab equals @xmath588 exceeding by factor of @xmath589 the driving force frequency in this frame .
the analysis of photon multiplicity for specific radiation sources further reduces to estimating in eq .
( [ dw - lab ] ) the driving force magnitude and period .
note that generally , the square of the integral in ( [ dw - lab ] ) is @xmath590 , while @xmath591 .
thus , in total , the radiation source brightness is proportional to its period . on the other hand , according to eq .
( [ 2gamma2 ] ) , the photon energy scales as @xmath592 , though that is generally compensable by an appropriate increase in @xmath1 .      for an electron crossing a family of atomic planes by a near - straight trajectory at a small misalignment angle @xmath266 , the time in eq .
( [ dw - lab ] ) expresses as @xmath593 , where @xmath594 is the coordinate transverse to the planes , and the oscillation period in the lab equals @xmath595 @xmath596 being the plane spacing .
therewith , eq .
( [ dw - lab ] ) becomes @xmath597 where @xmath582 as a function of time was replaced by the corresponding function of @xmath594 . if @xmath598 within a period describes roughly as @xmath599 evaluation of the integral in eq .
( [ dw - cb ] ) gives @xmath600}{\chi},\ ] ] where we used numerical estimates @xmath601 , and @xmath602 .
parameters @xmath100 and @xmath266 are not completely independent : they must obey inequality ( cf .
@xcite ) @xmath603 in order for the multiple scattering spread @xmath604 not to blur the coherent peak . at the same time , @xmath266 correlates with @xmath4 , since under assumption @xmath605 , eqs .
( [ 2gamma2 ] , [ t - cb ] ) yield @xmath606 for practical misalignment angles
@xmath607 , this constrains the electron energy to @xmath608 gev .
( [ cb - ineq ] ) and ( [ echi - ll20 ] ) are only compatible provided @xmath609 , which is essentially the condition ( [ all1 ] ) adopted in our paper .      for the case of channeling ,
the oscillation period depends on the particle energy : @xmath610 for a harmonic potential [ corresponding to force ( [ fx ] ) ] , the period will be constant , but with the account of anharmonicity , it may acquire also some dependence on the particle transverse energy in the channel .
for our estimates , it will suffice to adopt the harmonic approximation , letting @xmath611 where @xmath612 depends on the particle impact parameters , in the spirit of eq .
( [ fx ] ) . as an order of magnitude estimate
, we may use @xmath613 . substituting ( [ f - chann ] ) to eq .
( [ dw - lab ] ) , we get    @xmath614    at practice , the length of a crystal for channeling is limited by dechanneling effects . for positrons , the dechanneling length may be estimated crudely as @xcite @xmath615 in terms of which eq .
( [ dw - chann ] ) assumes the form @xmath616\right)^{3/2}\frac{l}{l_{\text{d}}}.\ ] ]      in an undulator , the particle is moving in a fixed channel , similarly to channeling , but the motion period , determined by the magnet spacing , is energy - independent , like in coherent bremsstrahlung . since @xmath562 is macroscopically large , according to eqs .
( [ dw - lab][2gamma2 ] ) , the source must be bright , but the photon energy be relatively low .
nonetheless , gamma - range undulators were created recently , with periods @xmath617 cm , and electron energies
@xmath618 gev @xcite .    for description of the undulator field strength ,
it is customary to introduce parameter @xmath619t[\text{cm}],\ ] ] where @xmath620 is the field amplitude , @xmath621 . in terms of @xmath622 , eq .
( [ dw - lab ] ) reads @xmath623 at practice , parameter @xmath622 can be @xmath624 ( somewhat violating the dipole radiation condition @xmath625 ) , whilst the number of periods can amount several hundreds or thousands .
the latter number is particularly large for undulator - based positron sources , where @xmath626 @xcite .
as had been pointed out in sec .
[ subsubsec : centr ] , beyond the central region specified by condition ( [ central - region ] ) , the gaussian approximation breaks down , and the multiphoton spectrum becomes dependent on all the detail of the single - photon spectrum .
even though in those regions the radiation spectrum is relatively faint ( and at high @xmath20 it may be overwhelmed by the incoherent bremsstrahlung ` tail ' ) , it may still prove measurable under high absolute radiation intensity .
therefore , it would be useful to find some approximations for peripheral regions , as well
.    suitable approximations for coherent radiation in peripheral regions can be built based on the same saddle - point approximation described in sec .
[ subsubsec : steepest - descent ] , provided transcendental saddle - point equation ( [ saddle - point - eq ] ) is solved at least approximately .
that requires development of a special approach on each side from @xmath533 .
let us consider the region @xmath627 first . from the viewpoint of eq .
( [ saddle - point - eq ] ) , this corresponds to large and negative @xmath394 ( see fig .  [
fig : s0c - vs - omega ] ) .
it is therefore beneficial to rewrite eq .
( [ saddle - point - eq ] ) as @xmath628,\ ] ] where the r.h.s . is smaller than any of the terms on the l.h.s .
( granted that the exponential within the integration domain is @xmath629 , achieving unity only at the endpoint @xmath630 ) .
thence , the r.h.s . can be treated as a perturbation .
graphical real solution of saddle - point equation ( [ saddle - point - eq ] ) , for @xmath264 given by eqs .
( [ dw1coh - real][f - def ] ) .
solid purple curve , exact solution .
dotted black line , central region approximation ( [ saddle - point - eq - lin - soln ] ) .
purple point , the central solution .
long - dashed black curve , large-@xmath20 approximation ( [ s0c - log-1 ] ) .
short - dashed black curve , improved large-@xmath20 approximation ( [ s0c - log-2 ] ) .
dot - dashed black curve , low-@xmath20 approximation ( [ quadr - eq ] ) . ]    to solve eq .
( [ saddle - point - eq - large - omega ] ) iteratively , first , we express @xmath631 from equating the l.h.s .
to zero ( leading log approximation ) : @xmath632 substituting this to the r.h.s . of ( [ saddle - point - eq - large - omega ] ) gives a nll approximation for @xmath631 : @xmath633+\mathcal{o}\left(\ln\ln\ln\frac{\omega}{b\omega_0 ^ 2}\right).\ ] ] its behavior is illustrated in fig .
[ fig : s0c - vs - omega ] by the long - dashed black curve . it does not seem to be accurate enough , which compels us to proceed to nnll : substitute @xmath631 in the r.h.s . of eq .
( [ saddle - point - eq - large - omega ] ) by that expressed from the l.h.s .
, whereupon replace the remaining @xmath631 in the r.h.s . by ll expression ( [ s0c - log ] ) : @xmath634^{1-z}\bigg\}.\end{aligned}\ ] ] the latter solution is to be inserted to eqs .
( [ h][d ] ) , which are then inserted to ( [ dwdomega - saddle ] ) .
the resulting expression is rather bulky , and will be omitted .
its accuracy can be visualized from fig .  [
fig : s0c - vs - omega ] , and may be regarded as satisfactory . the qualitative corollary from eq .
( [ s0c - log-2 ] ) is that the multiphoton pure coherent spectrum decreases by the law close to linear exponential [ since in ( [ dwdomega - saddle ] ) @xmath635 .
the opposite peripheral region @xmath636 corresponds to positive @xmath394 ( see fig .  [
fig : s0c - vs - omega ] ) .
there , the main contribution to the integral entering eq .
( [ saddle - point - eq ] ) stems from the lower integration limit .
asymptotically , at large @xmath394 , @xmath637 but the latter asymptotic solution becomes accurate at @xmath638 only , whereas we need to cover the whole region @xmath639 . therefore , eq .  ( [ saddle - point - eq ] ) needs to be solved accurately on the whole interval @xmath640 , which is challenging insofar as this equation is transcendental .    in the present situation , one of the simplest approaches may be to interpolate the solution between asymptotes ( [ 1/s2 ] ) and ( [ saddle - point - eq - lin ] ) , e.g. , writing it in form @xmath641 the behavior of the latter interpolation is illustrated in fig .
[ fig : s0c - vs - omega ] by the dot - dashed black curve , which is rather close to the exact solution on the interval of interest . expressing @xmath394 ( or @xmath642 ) from quadratic equation ( [ quadr - eq ] ) and inserting to eq .
( [ h][d ] ) , we get an approximation holding up to the central region ) . qualitatively , at large @xmath394 , @xmath643 , which may be interpreted in the sense that the function rises faster than a gaussian .
h. motz , p.  luchini .
_ undulators and free - electron lasers_. oxford : university press , 1990 ; p. rullhusen , x. artru , p.  dhez .
_ novel radiation sources using relativistic electrons_. world scientific , 1998 ;