the dynamics of relaxing quantum systems in sufficiently strong laser fields , which are neither stationary nor monochromatic , is an important aspect of the theory of interaction between atoms and photon fields .
this problem is of particular practical significance because of dramatic advances in the precision of measurements presently attainable in spectroscopic experimental studies of hydrogenic and few  particle atoms . among the simple atoms that attracted attention in the last few years one should mention : _ hydrogen _ and its isotopes @xcite , _ positronium _
@xcite , denoted @xmath2 , _ muonium _ ( see @xcite and references therein ) , denoted @xmath3 , and _ helium _ atom @xcite .
the measurements of the ground state hfs and the lamb shift are of particular interest for these systems .
high precision spectroscopy of low  lying transitions in hydrogenic atoms offers a unique opportunity to test qed calculations and to refine both the fundamental constants and various properties of respective nuclei .
this urges further theoretical developments intended to describe the well established photon - induced processes in simple atoms with much higher precision than techniques used so far are able to provide , in order to allow a proper interpretation of newly available experimental data . in the present paper the stepwise 3-photon ionization of muonium in a _
non  monochromatic _ laser field with a _ time
dependent _ amplitude is studied .
it has been recognized already a long time ago ( see @xcite for a more detailed discussion ) that the @xmath1 transition offers unique opportunities for high precision spectroscopy due to the narrow natural line width , @xmath4 , of the @xmath5-state . in the hydrogen atom , for example , @xmath6  hz which enables the quality factor of @xmath7 to be achieved already today , as well as suggests future resolutions of order of 1 part in @xmath8 .
experimentally , the @xmath1 transition can be induced doppler  free by absorbing two photons from two counter  propagating laser beams . in particular , this scheme is presently used at the rutherford appleton laboratory  @xcite for determination of the @xmath1 energy separation in muonium to the @xmath9  mhz accuracy .
this study offers an opportunity to improve the present knowledge of the muon / electron mass ratio and , thereby , of the muon mass itself . in the course of this investigation ,
two counter  propagating pulsed laser beams with almost identical intensities , @xmath10 , and the same photon energies , @xmath11 ( @xmath12 nm ) are used .
thus , the frequency @xmath13 is such that it is tuned into resonance with the 2photon doppler  free @xmath14 transition . unlike the hydrogen case however , the @xmath1 transition in muonium is hard to be observed directly , i.e. by detecting radiation emitted as a result of the @xmath5state deexcitation .
this is due to the fact that the appropriate line intensities in muonium happen to be weaker by several orders of magnitude , owing to much lower densities at which muonium atoms can be produced .
furthermore , the latter reason necessitates the use of higher laser intensities ( by a factor of @xmath15 , at least ) that are required for a reasonable signal strength .
this demands a _ pulsed _ rather than a cw laser source ( as in the case of hydrogen @xcite ) to be employed .
the @xmath5state of muonium is detected therefore via its photoionization by the third photon absorbed from one of the laser beams @xcite .
schematically , the experiment in muonium is presented in figure  [ fig1](a ) where @xmath16 stands for the resonant energy ( frequency ) which drives the 2photon @xmath1 transition .
despite the fact that the laser frequency can be calibrated , in principle , to a rather high accuracy , there are a number of systematic error sources each of which proves to be essential for precise determination of the @xmath1 energy separation in muonium . in this paper
we report the results of the theoretical study intended to allow for one of the most important among these effects , that is , the time - dependent frequency variation ( chirp ) of the laser field .
this phenomenon arises due to the rapid refraction index variation of the _ laser media _ and is hardly avoidable , especially with powerful pulsed lasers , unless some technical developments are made in order to compensate the frequency alteration  @xcite .
as will be demonstrated below , the chirped laser signal immediately makes its appearance in appreciable spurious shifting and broadening of spectral / photoionization lines , as well as leads to a noticeable distortion of the line shapes . to estimate the scale of the effect considered , one should note , for example , that the chirp  induced shift of the center of the ionization line is roughly equal to a characteristic magnitude of the chirp itself . for the @xmath1 experiment in muonium
, the latter can typically amount to @xmath17 , thus leading to the relative shift of the line s center by @xmath18 .
the effects of such an order are usually completely neglected in atomic physics .
nonetheless , these happen to be crucial for studies where an absolute accuracy of @xmath19 is anticipated .
it is therefore of primary practical importance to work out a relatively simple theoretical scheme which could allow one to describe certain types of multi - photon ionization phenomena in hydrogenic atoms , where an arbitrary time variation of the frequency of the laser field is accompanied by an arbitrary time modulation of its amplitude .
although the problem of the hydrogenic systems interaction with monochromatic resonant laser field with / without the amplitude modulation is investigated quite fully by now ( see @xcite and references therein ) , there have been much fewer studies where both the amplitude and the frequency of the laser signal vary with time .
even though , the authors confine themselves in most cases to either 1photon ( rather than 2photon ) resonant transition @xcite , or treat only somewhat special forms ( `` abrupt step '' , `` exponential field pulses '' etc . ) of the field modulation @xcite . under these circumstances , it seems very desirable to reconsider the problem as a whole , while adapting it specifically to the conditions of the @xmath1 experiment in muonium .
we consider , without a loss of generality , a set of three levels : @xmath20 and @xmath21 , of a hydrogenic atom with the charge @xmath22 of its nucleus , such that their one  particle energies satisfy the inequality : @xmath23 ; this rather general setup is shown in figure [ fig1](b ) .
it is supposed that the given 3level system is exposed to two counter  propagating ( along the @xmath24axis ) laser waves with equal time
dependent ( circular ) frequencies , @xmath25 , polarizations vectors , @xmath26 , and wave vectors , @xmath27 , where @xmath28 is the fine structure constant and @xmath29 denotes the speed of light .
the stationary part , @xmath30 , of the field frequency @xmath31 is assumed to be such that it is in the 2 and 1photon resonances with the pairs of states , @xmath32 and @xmath33 , respectively .
this implies that @xmath34 with @xmath35 being the differences of the one  particle energies .
in addition , it is supposed that the @xmath36transition is forbidden between @xmath37 and @xmath38 and is allowed between @xmath38 and @xmath21 .    under the resonance conditions assumed above , 3photon ionization of the system ,
that is , @xmath39 transition , occurs predominantly as a 2step process : ( i ) a 2photon resonant absorption from the state @xmath37 into the state @xmath38 , followed by ( ii ) a 1photon resonant transition between the levels @xmath38 and @xmath21 .
this simple physical picture is a consequence of the evident estimate of the probability for an atom to absorb 3 photons of equal energies ( the atomic units , @xmath40 , are used throughout the paper ; in this system , the speed of light is equal to @xmath41 ) : @xmath42 here , @xmath43 are the 2 and 1photon resonant transition probabilities , @xmath44 is the total width of the state @xmath38 , and @xmath45 denotes some characteristic mean energy which a two  fold summation over intermediate states in the exact expression to define @xmath46 ( see @xcite for a more detailed discussion ) can be reduced to . equation ( [ 3phot ] ) readily follows then from the estimate , @xmath47 , which is fulfilled for most states of atoms in moderately strong laser fields , except extraordinarily short  living ones .
a situation when @xmath37 stands for the ground , @xmath38 for the @xmath5 , and @xmath21 for the continuum @xmath48states , respectively , is of our primary concern here . in this case @xmath49 and @xmath50 since @xmath51 is dominated by the photoionization rate of the @xmath5level ( @xmath52 ) .
accordingly , we will assume that the natural line widths of the @xmath37- and @xmath38levels are both equal to zero : @xmath53 .
a straightforward generalization of the latter condition , that might be required for the treatment of _ excited _ states , can be achieved by adding appropriate imaginary parts to the energies of the levels  @xcite : @xmath54 .    in the semiclassical approximation , which happens to be accurate enough for our purposes ,
the electric fields inducing the 2 and 1photon transitions can be taken in the form @xmath55 here , the real functions @xmath56 and @xmath57 describe the ( macroscopic ) spatial inhomogeneity of the laser field and its time - dependent amplitude .
these are related direct to the laser intensity @xmath58 which , along with the chirped ( circular ) laser frequency @xmath31 , is obtainable direct from measurements . for the @xmath1 experiment with muonium , @xmath58 and @xmath31
can be well approximated as @xcite : @xmath59 where @xmath60  mj is the energy within one laser pulse ; @xmath61  ns , @xmath62  mm denote the temporal and spatial dispersions of the external laser field ( in the @xmath63 plane perpendicular to the direction of the beams propagation ) .
the full temporal width of the laser pulse is expressed then as @xmath64  ns .
parameter @xmath65 mhz defines the magnitude of the chirp @xmath66 , being actually detected in the experiment , and @xmath67 stands for the error function @xcite .    as
usually happens to be true in most practical situations , both the amplitude @xmath68 and the additonal phase @xmath69 vary much slower with time than @xmath70 and @xmath31(see equations ( [ intens])-([omega ] ) ) .
in addition , the magnitude of the chirp is normally much smaller than the `` base '' frequency @xmath30
. it is therefore permissible to assume that @xmath31 and @xmath68 are subject to the following general conditions @xcite : @xmath71 within the framework of the conventional time
dependent perturbation theory @xcite , the exact time
dependent wave function of the quasi
stationary state @xmath72 in the field is sought in the form @xmath73 where @xmath74 , being the amplitudes to find an atom in one of the unperturbed states @xmath75 , satisfy the following well known system of coupled differential equations : @xmath76 here , @xmath77 denotes the matrix element of @xmath78 being the operator of the `` particle  laser field '' interaction . due to
the given below reasons , @xmath78 can be taken in the form @xmath79 where @xmath80 is the dipole operator of an electron ( or a muon ) and @xmath81 denotes the electric field  ( [ fields ] ) seen by the atom in its rest frame .
for particular laser parameters of interest , the accuracy of the above approximation follows from the evident estimate : @xmath82 .
this implies that the contribution of the quadrupole terms @xmath83 , which originate from the spatial _ macroscopic _ inhomogeneity of @xmath81 , is expected to be a factor @xmath84 smaller than the contribution of the ( dipole ) terms in ( [ dip_appr ] ) .
almost the same estimate holds also for corrections arising due to the quadrupole component @xmath85 of the photon field itself .
this is despite the fact that the states @xmath37 and @xmath38 are coupled , in accord with the initial assumptions of our model and equation ( [ dip_appr ] ) , in only the second order of the perturbation theory in @xmath78 , i.e. by taking into account at least the 2photon absorption / emission . indeed , unlike the second order dipole coupling , the 1photon quadrupole @xmath86 in the amplitude .
if @xmath87 , as it is assumed throughout the paper , then the terms of such a type are known to be almost negligible @xcite . in our case , in particular , their relative contribution is of the order of @xmath88 , as can be shown by averaging the transition amplitude over some time interval : @xmath89 .
the same `` off - resonant '' argument is fully applicable to magnetic transitions as well , thus elucidating why the @xmath90 and @xmath91contributions are suppressed by a factor @xmath92 as compared to the dipole one .    by virtue of ( [ a])-([dip_appr ] ) , one can readily retrieve , after some conventional algebra , the following form for effective quasiclassical one  particle operators which couple two given states of reference @xcite : @xmath93 all summations are performed here over complete set of one - particle atomic orbitals
. equations ( [ vba])-([vab ] ) describe the real transitions , whereas ( [ vbb ] ) defines the _ ac stark _ shift of the level @xmath38 ( or @xmath37 , if @xmath94 ) which arises due to interaction of the atom with the incident radiation .
note that the sum in ( [ vbb ] ) coincides , by definition @xcite , with the dipole dynamic polarizability of the state @xmath38 at the frequency @xmath95 .
also , it should be pointed out that only the leading terms , as given by equations ( [ vba])-([vbb ] ) , were retained in the course of derivation .
the most important among neglected terms are of the form : @xmath96 } + c.c .
\right ) \right|^2 \simeq \frac{1}{\omega_{\rm l}^4 } \left|\dot{e}(t)\right|^2u^2(\bi{r}).\ ] ] apparently , this contribution is smaller than @xmath97 by a factor of @xmath98 , as follows from ( [ intens])-([omega ] ) , in accord with ( [ cond_omega])-([cond_e ] ) .
as was assumed above , there are no other discrete resonant atomic levels ( either 1 or 2photon ) except @xmath20 .
this enables the system ( [ a ] ) to be reduced to these two states only , and @xmath99 , is ignored , since we are mainly concerned here with the lowest @xmath100- and @xmath5states . ] while describing the 1st stage of the two
step 3photon ionization process under consideration .
indeed , an inclusion of the off - resonant states would give rise in the right - hand side of equations ( [ ala1])-([cc1 ] ) to the rapidly @xmath101varying terms @xmath102 .
after temporal averaging ( see above ) , these happen to be of order of @xmath103 , and can therefore be omitted to the accuracy adopted in this study .
on introducing the new functions , @xmath104 \right\ } , \label{alpha_c}\ ] ] and replacing the matrix elements in ( [ a ] ) by equations ( [ vba])-([vbb ] ) , one arrives at the following system of equations : @xmath105 , \\
\fl \label{cc1 } \dot{c}(t ) & = & - \rmi \:d_{b , a }   \frac{4\pi}{c }   i(\bi{r},t ) - \rmi \ : \left [ d_{b , b } \frac{4\pi}{c } i(\bi{r},t ) - \left ( 4\pi \delta \nu + \dot{\alpha}(t)+2 \dot{\phi}(t)\right ) \right ] c(t),\end{aligned}\ ] ] where @xmath106 is the time
independent frequency detuning off the resonance , @xmath107 is the intensity of the single laser beam and @xmath108 denote the effective matrix elements given by the sums in ( [ vba])-([vbb ] ) .
it has been taken into account here that the 3level system is actually probed by a superposition of two laser fields ( [ fields ] ) , rather than by a single one ; this results eventually in doubling all matrix elements .
in addition , we have neglected a contribution due to the 2photon absorption from _ each _ of two laser beams .
this contribution comes from the second  order interactions of an atom with only one of two light waves and makes its appearance finally in two doppler  broadened terms , i.e. depending on a @xmath109component of atomic velocity in the laboratory frame , in the expression to define ionization probabilities .
however , after averaging over @xmath109 , as would be required within the framework of a systematic approach , these terms turn out to be small compared , at the center of the photoionization line @xcite , with the doppler  free ones which originate from ( [ ala1])-([cc1 ] ) .
these arguments provide the ground for the approximation made , although the effects caused by atomic movement in the media are beyond the scope of our present consideration .
it is relevant to point out here that the matrix elements , @xmath108 , in ( [ ala1])-([cc1 ] ) are @xmath110dependent . in particular
, @xmath111 has a non  zero imaginary part .
the latter allows for a non
 zero photoionization rate of the @xmath38state , due to the action of the photon field , since @xmath30 is supposed to exceed the photoionization threshold @xmath112 .
indeed , as was noted above , @xmath111 is proportional to the dynamic ( tensor ) polarizability @xmath113 of the level @xmath38 , which is complex - valued if @xmath114 ( see our recent paper @xcite for a more detailed discussion ) .
@xmath115 and @xmath116 define then the _ ac stark _ shift of the level ( see equation ( [ vbb ] ) ) and its decay probability ( @xmath117 ) via the single photoionization .
this result is in agreement with the optical theorem  @xcite : @xmath118 , where @xmath119 is the fine structure constant and @xmath120 denotes the total photoionization cross section of the state @xmath38 .
numerical calculations show that the single photoionization happens to be the main channel of the @xmath38level depopulation .
it is therefore crucial for the correct description of the time
evolution of the 3level system that this depopulation mechanism is taken in fact into account in ( [ ala1])-([cc1 ] ) through the matrix elements @xmath108 .
an allowance for the 2nd step of the process considered can be made by augmenting the system ( [ ala1])-([cc1 ] ) with an additional equation that explicitly describes the @xmath121 coupling through the single photoionization of the state @xmath38 .
given by ( [ cont ] ) , this extra equation is nothing else but a differential form of the conservation law which determines the probability balance between the states of reference .
moreover , equations ( [ ala1])-([cc1 ] ) can be decoupled by substituting @xmath122 of ( [ ala1 ] ) into ( [ cc1 ] ) to finally get @xmath123 c(t )   \\
\lo+   \rmi d_{b , a } \frac{4\pi}{c }   i(\bi{r},t ) c^2(t ) , \nonumber   \\ \label{ala2 } \fl \dot{\alpha}(t ) =    \frac{4\pi}{c } i(\bi{r},t ) \left [ d_{a , a } + d_{b , a } c(t ) \right],\\ \label{cont }   \fl \dot{w}_c(t ) =   \frac{2}{\omega_{\rm l } }   \sigma^{(\gamma)}_b ( \omega_{\rm l } ) i(\bi{r},t ) \left| c(t)\right|^2   \exp \left\{2\im   \alpha(t ) \right\},\end{aligned}\ ] ] where @xmath124 is the probability for an atom to be ionized by the time @xmath101 , via absorption of 3 photons of equal energy . under the actual experimental conditions @xcite , initial population of the level @xmath38 , i.e. before the laser pulse is shot , is usually rather low .
this is because the major fraction of the muonium atoms is produced in the _ ground _ state @xmath37 where the population is proportional to @xmath125 . in view of this physical picture and by virtue of ( [ alpha_c ] )
, it is natural to impose the following initial conditions on the functions @xmath126 : @xmath127 these enable @xmath128)$ ] to be determined uniquely .
being of our primary concern , this quantity constitutes the required 3photon resonant ionization probability of an atom at infinitely large positive times , i.e. when the interaction between the single laser pulse and the system has already ceased . as a function of the laser frequency detuning @xmath129 , @xmath128)$ ] describes , and will be therefore used as a synonym of , the 3photon ionization line profile / shape . for brevity ,
the latter will be denoted as @xmath130)$ ] where the argument @xmath131 $ ] is introduced in order to indicate explicitly that the line profile depends on the form of the chirp . in particular ,
if @xmath132 then @xmath133)$ ] will be denoted as simply @xmath134 .
equations ( [ cc2])-([cont ] ) are well suited for numerical calculations with arbitrary functions defining the spatial and temporal distributions of the laser pulse . in addition
, the given system is rather convenient for analytic treatment , as will be demonstrated in subsection [ subsec : anal_cons ] . in particular
, one can develop further perturbative expansion of equations ( [ cc2])-([ala2 ] ) in terms of the intensity @xmath58 .
this yields analytic formulae which , after setting @xmath135 , agree with those formerly derived in @xcite within the framework of a 2level model , for _ unchirped _ laser signals . despite a relatively simple form of the appropriate equations ,
however , a not very straightforward numerical integration is still required in order to calculate , for example , photoionization line profiles in most practical situations , let alone the fact that the result of the work @xcite is valid for rather weak laser intensities only . note that equations ( [ cc2])-([cont ] ) are free from the latter restriction .
furthermore , unlike our approach employing the probability amplitudes , the problems similar to those considered in this paper are usually treated ( e.g. see @xcite ) in the formalism of the density matrix @xcite . under the physical conditions adopted here , these two approaches are equivalent . to facilitate the adequate comparison however , it is relevant to present here an alternative form of ( [ cc2])-([cont ] ) , by rewriting the system in terms of the density matrix elements : @xmath136 and @xmath137 , denoted as @xmath138 .
this yields : @xmath139 here , 5 unknown functions are subject to the following initial conditions : @xmath140 it should be pointed out that only real quantities enter equations ( [ r2s])-([r_cont ] ) .
this circumstance may be advantageous , especially for numerical calculations , since most standard numerical packages are not applicable direct to systems of ode involving complex  valued coefficients .
apart from that , ( [ cc2])-([cont ] ) and ( [ r2s])-([r_cont ] ) generate identical results , as has been proved by extensive numerical tests carried out by us for laser signals with various forms of the chirp , pulse envelope , and spatial distribution .
however , before we demonstrate appropriate results ( see section [ sec : appl ] ) , it is instructive to present here an approximate but reasonably accurate analytic solution for @xmath141 .
this consideration seems to be rather useful for a qualitative description of the process as a whole , since it reveals some key features which make their appearance due to the non - stationarity of both the frequency and the amplitude of the laser signal .
let us assume , in agreement with ( [ intens ] ) , that the intensity of the incident radiation is gaussian both in time and space and use the most general form of the chirp , @xmath142 ( cf . equation ( [ omega ] ) ) . to start with , consider riccati  type equation ( [ cc2 ] ) to determine the complex  valued function @xmath143 , equation ( [ cc2 ] ) takes the form @xmath144 c(x )   \nonumber \\ \label{eq : full } \lo+   \rmi d_{b , a } { \cal i } \e^{-\frac{\tau^2 x^2}{2 \sigma^2_{\rm t}}}=0,\end{aligned}\ ] ] where the following notations were introduced for brevity : @xmath145 note that @xmath146 depends on the time
independent parameters of the pulse only , but not on @xmath147 .    due to its non - linearity , equation ( [ eq : full ] )
admits only numerically solving , except for the trivial case , @xmath148 , where @xmath141 is of the form @xmath149 - \delta d \biggr\}. \label{c0}\end{aligned}\ ] ] however , the problem can be simplified considerably by linearizing equation ( [ eq : full ] ) .
this approximation , which the consideration developed below is stemmed on , is justifiable if @xmath150 holds uniformly for @xmath151 and arbitrary @xmath152 and @xmath153 . to satisfy the latter condition ,
it is sufficient ( but not necessary ) to assume that @xmath154 as readily follows from equation ( [ c : simp ] ) . if @xmath155 and @xmath156 , for example , then @xmath157
this imposes the following restriction on the pulse energy : @xmath158 , where particular values of the atomic matrix element , @xmath159 ( see table [ tabd ] ) , and the typical laser parameters have been used .
the above estimates are illustrated in figure [ fig1a ] where @xmath160 obtained by direct numerically solving equation ( [ eq : full ] ) is plotted versus the dimensionless variable @xmath147 , at several @xmath161values .
in addition , for each pulse energy , the appropriate @xmath162value defined by equation ( [ cond_c ] ) is indicated for reference as well .     as a function of @xmath147 , for various pulse energies , @xmath161 , and fixed values of the chirp , @xmath163 , and the frequency detuning , @xmath164 .
the laser field frequency is assumed to be of the form : @xmath165 ; the rest laser parameters used are : @xmath166.,width=302,height=302 ]    the linearized form of equation ( [ eq : full ] ) reads : @xmath167   \widetilde{c}(x ) + \rmi d_{b , a}\ : { \cal i } \e^{-\frac{\tau^2 x^2}{2 \sigma^2_{\rm t}}}=0.\end{aligned}\ ] ] its solution obeying the zero initial condition ( [ cond_inf ] ) at @xmath168 has the form @xmath169 to demonstrate the fair accuracy of this relation for @xmath170 and @xmath171 , for example , the graph of @xmath172 as a function of @xmath147 is presented , for several chirp values , in figure [ fig1b ] , along with appropriate graphs for numerical solution of the non - linear equation ( [ eq : full ] ) . together with plots in figure [ fig1a ] , these graphs provide instructive information about the magnitude of @xmath173 at various time moments .
in addition , both figures clearly demonstrate a rather peculiar way in which the values for the above ratio are ( strongly ) distorted due the presence of the chirp .     and @xmath174 as functions of @xmath147 , for fixed values of the pulse energy , @xmath175 , frequency detuning , @xmath176 , and for several chirp values : @xmath177 . the frequency of the non - monochromatic laser field
is assumed to be of the form : @xmath178 , and the rest laser parameters used are : @xmath179 .
+ @xmath180  numerical solution of equation ( [ eq : full ] ) ; @xmath181  analytic solution , @xmath182 , of equation ( [ eq : simp ] ) , as given by ( [ c : simp]).,width=302,height=302 ]    it is relevant to note that @xmath183 approximates the function @xmath141 qualitatively correctly for even those @xmath161 which violate condition ( [ cond_c ] ) .
an elaborate comparison with exact numerical solutions of ( [ eq : full ] ) shows that , for most applications , the region of the validity of equation ( [ c : simp ] ) can be safely extended up to @xmath184 , depending on the rest laser parameters . at higher pulse energies ,
the contribution of non  linear effects becomes so essential that the neglect of the @xmath185term is no longer permissible .
this can be seen , for example , in figure [ fig1a ] where , for @xmath186 and @xmath187 , appropriate non  linear contribution amounts to 20% .    as follows from the method of our derivation , equation ( [ c : simp ] ) can be readily generalized on a wide class of similar multi - photon ionization processes that are induced by not very intense laser fields and which involve the chirps and laser amplitudes , such that @xmath188 and @xmath189 sufficiently fast as @xmath190 .
note that the latter condition ensures the rapid convergence of the integral in ( [ c : simp ] ) .
in addition , it makes certain that the values of all physical quantities of interest are determined actually by only a somewhat narrow domain of @xmath147 within which the laser intensity is peaked .
equation ( [ c : simp ] ) is suitable to be used further for the analytical study of the 3photon ionization with muonium , where @xmath191 and @xmath192 . in particular , @xmath183 readily enables one to obtain analytic formulae for the probabilities : @xmath193 , @xmath194 , and @xmath195)$ ] . however , these expressions are not given explicitly here as they turn out to be too bulky . instead
, the formula ( [ c : simp ] ) will be used in the following subsection [ subsec : wood ] , while addressing an interesting complementary problem : a calculation of the integral of the 3photon ionization line profile over the entire range of the frequency detunings , @xmath196 .      in what follows ,
it is shown that , for given @xmath161 and @xmath153 , the integral of the 3photon ionization line profile , @xmath133)$ ] , over the entire range of the frequency detunings @xmath129 ,
@xmath197)\equiv \int_{-\infty}^{+\infty } w_c^{\infty}(\delta \nu , a_\omega , r;[\phi ] ) \rmd \delta \nu,\ ] ] is almost independent of particular form of the chirp , at least for low pulse energies . to be more precise , equation ( [ eq : rho_gen ] ) defines in fact a non - linear functional of @xmath69 whose `` strength '' depends , though , on the rest laser parameters . within certain range of these parameters
( @xmath161 , in the first instance ) , the functional relation ( [ eq : rho_gen ] ) turns out to be weak , so that @xmath198)$ ] plays in this case the role of a sort of adiabatic invariant for 3photon resonant ionization processes occurring under the action of the fields with slow  varying frequencies .
hence , an `` almost  conservation '' of @xmath198)$ ] provides a simple and easily verifiable approximate criterion which enables the photoionization line profiles corresponding to different values / forms of the chirp to be compared and normalized .
the problem discussed here was first addressed , to our knowledge , in @xcite , within the framework of the similar 3level model and for low pulse energies only .
in contrast with our present approach , however , the genuine photoionization rate of the level @xmath38 ( that is , @xmath5 ) was set in @xcite equal , respectively , to zero and to some time and energy
independent constant , while describing the dynamics of the 1st and the 2nd stages of the resonant 3photon ionization .
although this approximation provides acceptable results for @xmath199 , say , it is of interest , both experimental and theoretical , to reconsider the entire problem by releasing the above simplifying assumptions and getting , thereby , a deeper insight into the physics of the process .    [ [ section ] ]    let us consider first the simplest case of low @xmath161 , such that condition ( [ cond_c ] ) is satisfied . in addition , we will still assume , without a loss of generality , that the intensity of the incident radiation is gaussian both in time and space , without specifying , though , particular form of the chirp for a moment . by making use of equation ( [ cont ] ) ,
this allows @xmath198)$ ] to be expressed in terms of @xmath183 as @xmath200 ) = \frac{c}{4\pi }          \frac{2}{\omega_{\rm l } } \sigma^{(\gamma)}_b { \cal i } \int_{- \infty}^{+ \infty } \rmd \delta \nu   \int_{- \infty}^{+ \infty } \rmd x   \exp \left\{- \frac { \tau^2 x^2}{2 \sigma^2_{\rm t } } \right\ } \left|\widetilde{c}(x,\delta \nu , a_\omega , r)\right|^2 \nonumber \\ \label{eq : rho - cappr } \lo \times \exp \left\ { 2 d_{b , a } { \cal i } \im           \left [ \int_{-\infty}^{x }     \exp \left\{- \frac { \tau^2 u^2}{2 \sigma^2_{\rm t } } \right\ } \widetilde{c}(u,\delta \nu , a_\omega , r ) \rmd u          \right ] \right\}.\end{aligned}\ ] ] here , our former notations ( [ defin : iab ] ) have been employed . to this end , it must be noted that @xmath201 is satisfied for all @xmath129 and @xmath202 , except those @xmath203 which are in the vicinity of the lineshape s maximum , @xmath204 , and for @xmath161 such that @xmath205 . the latter condition happens to be more restrictive as compared to the one assumed here .
this is due the fact that the diagonal matrix elements are usually much bigger than the non - diagonal ones ( see table [ tabd ] for comparison ) . by making use of ( [ c : simp ] ) , ( [ cond_omega ] ) and presuming that the above condition is fulfilled
, one can show that @xmath206 \approx \nonumber \\ \label{eq : exp } \fl - \sqrt{2 \pi } d^2_{b , a } { \cal i}^3 \frac{\sigma_{\rm t}}{\tau }   \frac{\im(\delta d ) } { \left[4 \pi \tau \delta \nu + 2 \dot{\phi}(0 ) + { \cal i } \re ( \delta d)\right]^2 + { \cal i}^2 [ \im ( \delta d)]^2 } \left [ 1+{\rm erf } \left ( \frac{\tau x } { \sqrt{2}\sigma_{\rm t}}\right ) \right].\end{aligned}\ ] ] note that , for any @xmath147 , arbitrary laser parameters and @xmath207 , this expression is negative , as it should be , since @xmath208 . for small pulse energies considered here , the right  hand side of ( [ eq : exp ] ) is proportional to @xmath209 , which permits the appropriate exponent in ( [ eq : rho - cappr ] ) to be expanded in terms of @xmath161 . on retaining two leading terms in this expansion , this yields @xmath200 ) = \frac{c}{4\pi }          \frac{2}{\omega_{\rm l } } \sigma^{(\gamma)}_b ( \omega_{\rm l } ) { \cal i } \int_{- \infty}^{+ \infty } \rmd \delta \nu   \int_{- \infty}^{+ \infty } \rmd x   \exp \left\{- \frac { \tau^2 x^2}{2 \sigma^2_{\rm t } } \right\ } \left|\widetilde{c}(x,\delta \nu , a_\omega , r)\right|^2 \nonumber \\
\label{exp_ser } \fl \times \left\ { 1 - \frac{\sigma_{\rm t}}{\tau } \frac{\sqrt{2 \pi } d^2_{b , a } { \cal i}^3   \im(\delta d ) } { \left[4 \pi \tau \delta \nu + 2 \dot{\phi}(0 ) + { \cal i } \re ( \delta d)\right]^2 + { \cal i}^2 [ \im ( \delta d)]^2 } \left [ 1+{\rm erf } \left ( \frac{\tau x } { \sqrt{2}\sigma_{\rm t}}\right ) \right ] \right\}.\end{aligned}\ ] ] further calculations are straightforward but cumbersome .
hence , some technical details will be given below for the first term in the curly brackets in ( [ exp_ser ] ) only ; the corresponding contribution is denoted as @xmath210 ) $ ] . in this simplest case , by substituting ( [ c : simp ] ) for @xmath183 and interchanging the order of integration , one arrives at @xmath211 ) = \nonumber \\ \fl   \frac{c}{4\pi }   \frac{\sigma^{(\gamma)}_b ( \omega_{\rm l}){\cal i}^3}{\tau \omega_{\rm l } } d^2_{b , a }   \int_{-\infty}^{+\infty }   \exp\left\ { - \frac { \tau^2 x^2}{2 \sigma^2_{\rm t } } -   \sqrt{2 \pi } \frac{\sigma_{\rm t}}{\tau } { \cal i } \im(\delta d )   { \rm erf}\left ( \frac{\tau x } { \sqrt{2}\sigma_{\rm t}}\right )   \right\ } \rmd x \nonumber \\ \lo{\times } \int_{-\infty}^{x }   \exp\left\ { - \frac { \tau^2 u^2}{\sigma^2_{\rm t } } +   \sqrt{2 \pi } \frac{\sigma_{\rm t}}{\tau } { \cal i } \im(\delta d )
{ \rm erf}\left ( \frac{\tau u } { \sqrt{2}\sigma_{\rm t}}\right )   \right\ } \rmd u \nonumber \\
\label{eq : wood1 } \fl = \frac{c}{4\pi }   \frac{\sigma_{\rm t } \sigma^{(\gamma)}_b { \cal i}^2}{\sqrt{2 } \tau^2 \omega_{\rm l } }   \frac{d^2_{b , a}}{\im(\delta d ) }   \left\{\sqrt{\frac{\pi}{2 } } - \int_{-\infty}^{+\infty } \exp\left [ - 2 u^2 - \sqrt{2 \pi } \frac{\sigma_{\rm t}}{\tau } { \cal i } \im(\delta d )   { \rm erfc}(u )   \right ]   \rmd u \right\ } .\end{aligned}\ ] ] for two most important cases of interest , @xmath212 and @xmath213 , the integral in ( [ eq : wood1 ] ) can be evaluated , respectively , by expanding the integrand into the power series in @xmath214 and by using the method of the steepest descent .
this yields @xmath215   \rmd u   \nonumber \\ \label{step : desc } \lo{\asymp }
\left\ { \begin{array}[c]{ll }   \sqrt{\pi/2 } - \pi ( \sigma_{\rm t } / \tau ) \im(\delta d ) { \cal i }          \qquad & { \rm for } \ ;
|{\cal i } \im(\delta d ) | \ll 1   \\ \ms \left ( \sqrt{\pi }   \tau^2 / 2^{3/2 } \sigma^2_{\rm t } \right )                                   \left [ \im(\delta d )   { \cal i } \right]^{-2 }           \qquad & { \rm for }   \ ; |{\cal i } \im(\delta d ) | \gtrsim 1       \end{array }     \right . .\end{aligned}\ ] ] accordingly , we are left finally with @xmath211 )    \nonumber \\
\label{rho_0 } \fl = \frac{c}{4\pi } \sqrt{\frac{\pi}{2 } } \frac{\sigma_{\rm t } } { \tau^2 \omega_{\rm l } }   \sigma^{(\gamma)}_b ( \omega_{\rm l } ) d^2_{b , a }   \times \left\ {          \begin{array}[c]{ll }           \sqrt{\pi } ( \sigma_{\rm t}/\tau ) { \cal i}^3           \qquad & { \rm for } \ ;
|{\cal i } \im(\delta d ) | \ll 1 \\ \ms          [ \sqrt{2 } \im(\delta d ) ] ^{-1 } { \cal i}^2          \qquad & { \rm for }   \ ; |{\cal i } \im(\delta d ) | \gg 1          \end{array } \right .
.\end{aligned}\ ] ] note that in either limiting case , the result is independent , in accord with our initial statement , of particular form of the chirp @xmath216 . to clarify this point
, one should emphasize that the only two simplifications made in the course of derivation were : ( i ) the use of @xmath183 given by ( [ c : simp ] ) and ( ii ) the replacement of the appropriate exponent by a unity .
both approximations are fully justifiable for @xmath217 assumed here .
obviously , the latter condition is consistent with the _ first _ relation in ( [ rho_0 ] ) , which has been formerly obtained by means of a different technique ( and cast in a slightly different analytic form ) in @xcite .
apart from being chirp  independent , this result demonstrates also that @xmath210 ) \propto a^3_{\omega}$ ] , as opposed to the next limiting case , that is , @xmath218 , where @xmath210 ) \propto a^2_{\omega}$ ] . a gradual decrease of the exponent is clearly evidenced by numerical simulations carried out by us for muonium , while scanning sufficiently broad @xmath161domain .
some results of this numerical study are discussed in more detail in section  [ sec : appl ] ( see figure [ fig : power ] , for example ) .
these show , in particular , that equations ( [ rho_0 ] ) describe @xmath210)$ ] qualitatively correctly for both low and moderately high pulse energies .    in a similar manner
one can also obtain appropriate contribution , denoted @xmath219)$ ] , coming from the second term in the curly brackets in ( [ exp_ser ] ) .
after tedious calculations whose details will be given elsewhere , one ends up with the following expression for a sum of two contributions : @xmath211 )   + \rho^{(1)}(a_\omega , r;[\phi ] )   \nonumber \\
\label{rho_01 } \fl = \left\ {           \begin{array}[c]{ll }   \frac{c}{4\pi } \frac{\pi}{\sqrt{2 } } \frac{\sigma^2_{\rm t } } { \tau^3 \omega_{\rm l } }   \sigma^{(\gamma)}_b ( \omega_{\rm l } ) d^2_{b , a } { \cal i}^3   \left (    1 - \sqrt{2 } \pi \frac{\sigma^2_{\rm t } } { \tau^2 } d^2_{b , a } { \cal i}^2     \right )           \qquad & { \rm for } \ ;
|{\cal i } \im(\delta d ) | \ll 1 \\
\ms \frac{c}{4\pi } \frac{\sqrt{\pi}}{2 } \frac{\sigma_{\rm t } } { \tau^2 \omega_{\rm l } }   \sigma^{(\gamma)}_b ( \omega_{\rm l } )   \frac {   d^2_{b , a } } { \im(\delta d ) } { \cal i}^2 \left (    1 - \frac{4 \sqrt{2 } \pi } { 3 }   \frac{\sigma^3_{\rm t } } { \tau^3 }   d^2_{b , a } { \cal i}^2 \right )          \qquad & { \rm for }   \ ; |{\cal i } \im(\delta d ) | \gtrsim 1          \end{array } \right .
.\end{aligned}\ ] ] here , according to ( [ cond_omega ] ) , we have discarded the terms in @xmath219 ) $ ] containing the second and higher order derivatives of @xmath220 .
this simplification results in the fact that entire equation ( [ rho_01 ] ) turns out to be completely independent of the chirp .    the right  hand side of equation ( [ cond_omega ] ) depends ( through @xmath202 ) on the radius @xmath153 . to eliminate this dependence and to facilitate , thereby , an adequate comparison with the results of numerical calculations , the left and right  hand sides of ( [ cond_omega ] ) can be integrated over the entire @xmath221plane perpendicular to the direction of the laser beam s propagation ( @xmath24axis ) .
a clear physical meaning of such an averaging procedure is discussed in the comments following equation ( [ w_avar ] ) .    by virtue of the identity , @xmath222 which defines the spatially  independent part of @xmath202 , denoted as @xmath223 , this yields eventually the following @xmath153independent result : @xmath224 ) \approx   \int_{-\infty}^{+\infty}\!\!\int_{-\infty}^{+\infty } \left(\rho^{(0)}(a_\omega , r;[\phi ] )   + \rho^{(1)}(a_\omega , r;[\phi ] ) \right )
\rmd x \rmd
y \nonumber \\ \label{rho_fin } \fl = \left\ {           \begin{array}[c]{ll }   \frac{c}{4\pi } \frac{\sqrt{2 } \pi^2}{3 } \frac{\sigma^2_{\rm t } \sigma^2_{\rm r } } { \tau^3 \omega_{\rm l } }   \sigma^{(\gamma)}_b ( \omega_{\rm l } ) d^2_{b , a } \widetilde{{\cal i}}^3   \left (    1 - \frac { 3 \sqrt{2}\pi } { 5 } \frac{\sigma^2_{\rm t } } { \tau^2 } d^2_{b , a}\widetilde{{\cal i}}^2     \right )           \qquad & { \rm for } \ ; |\widetilde{{\cal i } } \im(\delta d ) | \ll 1 \\ \ms \frac{c}{4\pi } \frac{\pi^{3/2}}{2 } \frac{\sigma_{\rm t } \sigma^2_{\rm r } } { \tau^2 \omega_{\rm l } }   \sigma^{(\gamma)}_b ( \omega_{\rm l } )   \frac {   d^2_{b , a } } { \im(\delta d ) } \widetilde{{\cal i}}^2 \left (    1 - \frac{2 \sqrt{2 } \pi } { 3 }   \frac{\sigma^3_{\rm t } } { \tau^3 }      d^2_{b , a } \widetilde{{\cal i}}^2 \right )          \qquad & { \rm for }   \ ; |\widetilde{{\cal i } } \im(\delta d ) | \gtrsim 1          \end{array } \right . .\end{aligned}\ ] ]    , over the entire range of the frequency detunings @xmath225 , as a function of the chirp s magnitude , @xmath226 ; three particular pulse energies are considered : @xmath227 .
the chirped frequency of the laser field employed is of the form : @xmath228 .
the rest laser parameters used are : @xmath229.,width=283,height=283 ]    in order to check this equation , the numerical calculation of @xmath230)$ ] , in a wide range of @xmath161 and for typical values of the rest physical parameters used in the 1s2s experiment in muonium , have been carried out .
this was done by an extra numerical integration of the 3photon ionization line profiles @xmath231 ( see equation ( [ w_avar ] ) for definition ) over the entire @xmath129domain .
these results are shown in figure [ fig1c ] , for three particular pulse energies . at @xmath175 ,
being the only value displayed which satisfies condition ( [ cond_c ] ) , the plot demonstrates that @xmath230)$ ] is conserved to the 1% accuracy , within a wide range of the chirp s magnitude . the mean numerical value , @xmath232 ) \approx
0.080\:{\rm mhz}$ ] , is to be compared with that obtained by means of the _ first _ equation in ( [ rho_fin ] ) : @xmath233 ) \approx 0.1242 - 0.0292=0.095\:{\rm mhz}$ ] . to make these results comparable ,
the latter value includes also an additional weighting factor , @xmath234 . even though parameter @xmath162 is very close to unity : @xmath235 , for @xmath175 and @xmath236 used here , numerical and analytical data are consistent to 20% of relative accuracy .
interesting to note that the first / second relation in ( [ rho_fin ] ) overestimates / underestimate the true value of @xmath230 ) $ ] .
much better agreement can be anticipated ( and actually takes place ) for lower pulse energies and/or bigger spatial dispersions @xmath237 . for @xmath238 and @xmath239 (
apppropriate line profiles are displayed in figure [ fig2](a ) ) , for example , numerical and analytical results read : @xmath240 ) \approx 0.90\cdot 10^{-4}\:{\rm mhz}$ ] and @xmath241 , respectively .
[ [ section-1 ] ]    for higher pulse energies and/or smaller spatial dispersions @xmath237 , such that @xmath242 , equation ( [ rho_fin ] ) fails to describe adequately the true function @xmath230)$ ] .
this refers both to the absolute values of @xmath230)$ ] and , particularly , to its dependence on the magnitude of the chirp , which makes its appearance at high pulse energies .
it has been elucidated above that the main underlying reason for this failure originates from the fact that `` the linear approximation '' ceases to be valid any longer , due to corrections caused by the @xmath243term in equation ( [ eq : full ] ) .
this non  linear term does not allow any reasonably accurate analytic treatment to be developed , so that the numerical analysis of the problem must be used instead .
to simulate the high  energy regime , the results of our calculation are presented in figure [ fig1c ] , for two pulse energies , @xmath244 , and for typical values of the rest physical parameters used in the 1s2s experiment in muonium . as in the case of low @xmath161 ,
this was done by the numerical integration of the 3photon ionization line profiles @xmath245 over the entire @xmath129domain .
the plots demonstrate in fact a noticeable deviation of @xmath230)$ ] from a constant value , with an increase of the pulse energy .
for example , at @xmath187 and for @xmath226 in the range @xmath246 , this deviation ( that is , the dependence on the chirp s magnitude ) amounts already to 30% ; in absolute units , the data read : @xmath247)=2.265\:{\rm mhz}$ ] for @xmath248 and @xmath247)=2.90\:{\rm mhz}$ ] for @xmath249 .
furthermore , the gauge of the given deviation may be expected to be even more pronounced at higher pulse energies , thus preventing us actually from viewing @xmath230)$ ] as a conserving quantity in this case .
although such @xmath161 are of no practical interest at the moment for the @xmath100@xmath5 experimental study with muonium , the entire effect might be of relevance for future experimental studies where much higher laser intensities are involved .
let us apply the results obtained in preceding sections to particular laser parameters adopted in the 1s-2s experiment in muonium @xcite , by assuming that @xmath250 and setting @xmath251 .
it is convenient in this section to measure the energies in the units : @xmath252 , where @xmath253 stands for the reduced mass of the electron and muon . here
, the value for the electron to @xmath254 mass ratio , @xmath255 , was used .
the quantity @xmath256 stands for the unity of distance , with @xmath257 being the bohr radius .
accordingly , one  particle energies and dynamic polarizabilities , @xmath258 , are measured in the units of @xmath259 and @xmath260 .    as was mentioned above ,
3photon ionization probability , @xmath261 , can be obtained by accurate numerically solving either the system ( [ ala2])-([cont ] ) or ( [ r2s])-([r_cont ] ) . in both cases
this has been done by the stepwise time integration , starting from @xmath262 , while scanning sufficiently broad @xmath129domain centered at @xmath263 ( i.e. @xmath264 ) , for all @xmath265 on the @xmath153mesh : @xmath266 .
numerical values of the matrix elements , @xmath267 , which enter ( [ ala2])-([cont ] ) were calculated in @xcite ; these are compiled in table [ tabd ] , together with the value of the single photoionization cross section of the @xmath5level at @xmath268 .
the latter is given explicitly by ( see @xcite and references therein ) @xmath269 where @xmath119 is the fine structure constant and @xmath270 denotes the ionization potential of the @xmath5state .
@xmath271 [ tabd ]    at the final stage of our numerical procedure , a spatial averaging has been performed in order to obtain @xmath153-independent ionization profiles , @xmath272 .
this has been achieved by an extra integration of @xmath273 over entire @xmath221plane perpendicular to the direction of the laser beam s propagation : @xmath274 basically , this newly introduced quantity can be interpreted as an averaged ionization probability related to the entire laser beam spot , since @xmath245 is independent of the distance from the beam s axis .
alternatively , in view of the low values of the muon density in the media @xcite ( see below ) , @xmath275 can be called as `` the ionization probability per one ( muon ) atom in the beam '' .
an auxiliary factor , @xmath276 , has been introduced in ( [ w_avar ] ) in order to make @xmath277 dimensionless , as the original probability @xmath261 is .
this particular choice was adopted according to the so  called rule being inherent to various problems involving gaussian law . in the case
considered , the `` @xmath278''rule makes its appearance through the fact that the values of @xmath261 happen to be almost negligible for @xmath279 , unless the pulse energies higher than @xmath280 , say , are considered .
for such a high energy regime , the probability @xmath273 has a long  distance `` tail '' spreading out beyond the _ effective _ beam radius , @xmath281 , within which the values of @xmath261 may be very close to , yet less than , @xmath9 .
in this situation the proper normalization factor different from that chosen above must be used so as to prevent @xmath282 from being bigger than @xmath9 . as is clearly evidenced by our calculations , the values of @xmath277 defined by equation ( [ w_avar ] )
exceed 1 starting already from @xmath283 , depending on the rest laser parameters .
one should point out , however , that such @xmath161 are hardly attainable currently with proper pulsed lasers sources in the required frequency range .
this makes the high energy regime in the resonant 3photon ionization with muonium to be mostly of academic rather than practical interest at the moment , thus justifying the definition of equation ( [ w_avar ] ) . to avoid any misunderstanding ,
it must be clearly stated once again that the spurious effect mentioned here originates solely from the particular choice of the normalization factor in ( [ w_avar ] ) , and has , therefore , nothing to do with either the model developed in the work or a lack of the numerical accuracy employed by us .
the latter was chosen to be equal to @xmath84 which proves to be sufficient to ensure that @xmath284 holds ( as it should ! ) for any particular individual set of all parameters involved in the problem , including @xmath285 , and @xmath237 .
the averaged probability @xmath277 readily enables one to estimate the expected number of experimentally observed ionization events , @xmath286 , defined as the number of the muonium atoms which are ionized by a sequence of @xmath287 identical gaussian laser pulses ( [ intens ] ) within the time @xmath288 and detected eventually in apparatus .
such a formulation of the problem corresponds direct to the actual experimental situation where the gas media containing muonium atoms interacts with @xmath289 laser pulses per second . under these conditions ,
the required number of the ionization events to be detected during the observation time @xmath290 can be estimated by means of the following simple relation : @xmath291 here , @xmath292 is the spatial density of the muonium atoms in the chamber , @xmath293 and @xmath294 stand , respectively , for the _ effective _ cross section of the laser beam ( see above ) and the total path which this beam travels in the media , @xmath295 denotes an efficiency of registration of the muonium atoms in apparatus ; finally , @xmath296 is the value of the averaged 3photon ionization probability @xmath297 at its maximum , for particular values of the chirp and the energy of the pulse .
provided @xmath226 and @xmath161 are fixed , this maximum is achieved approximately at @xmath298 , almost independently of @xmath161 , as can be seen in figure [ fig2](b ) .
in particular , @xmath299 thus leading , for example , to the following expected number of events detectable during the observation time @xmath300 : @xmath301 . it is rather instructive that this value happens to be quite close to that found preliminary in the @xmath1 experiment in muonium being currently underway .
, of the ground state of muonium by the laser field with the chirped frequency @xmath302 , versus laser frequency detuning , @xmath303 .
two pairs of different values for the pulse energy , @xmath161 , and the spatial dispersion of the laser signal , @xmath237 , are used : ( a ) @xmath304 , and ( b ) @xmath305.,width=529,height=283 ]    , of the ground state of muonium by the laser field with the chirped frequency @xmath306 , versus laser frequency detuning , @xmath307 , for two pulse energies : ( a ) @xmath308 and ( b ) @xmath309 .
, width=529,height=283 ]    .
the chirped frequency of the laser field is assumed to be of the form : @xmath310 ; ( b ) - the power  dependence of the line profile s maximum value . in both cases
the rest laser parameters are the same as in figures [ fig3](a , b ) .
, width=529,height=264 ]    some results of our simulations are presented in figures [ fig2]-[fig5 ] , for various magnitudes of the chirp @xmath226 and various parameters of the laser pulse which are typically used in the measurements of the @xmath1 energy separation in muonium ; particular values employed are indicated in figures [ fig2],[fig3 ] .
both these graphs demonstrate appreciable red  shift of the maxima ( i.e. towards lower frequencies , relative to @xmath311 ) of the photoionization profiles with an increase of the chirp . it must be noted , however , that this shift is caused by a combination of two competing mechanisms : ( i ) the shift being due to the combined _ ac stark _
effect for the @xmath100- and @xmath5levels and ( ii ) that arising because of the chirp itself .
the curves show that the former mechanism turns out to be small compared with the latter at relatively small pulse energies , @xmath312 . at @xmath313 and @xmath314 , for instance , i.e. in the case of completely unchirped laser signal , the _ ac stark _ shift amounts approximately to @xmath19 .
one can anticipate , however , that the relative contribution of the _ ac stark _ shift will increase as the laser intensity increases , since its magnitude is linearly proportional to the laser intensity , whereas the chirp  induced shift is almost intensity  independent .
this behaviour can be seen in figures [ fig3](a , b ) and is summarized in figure [ fig4](a ) showing , in particular , that the _ ac stark _
shift becomes equals to @xmath315 at @xmath316 .
in addition to the shift of the maxima , the curves in figures [ fig2],[fig3 ] exhibit a strong dependence on both the spatial dispersion @xmath237 and the power @xmath161 of the laser signal . it should be noted that an appreciable distortion of the photoionization line shapes happens to be particularly enhanced at relatively large chirp values , where the line profiles become asymmetric .
generally speaking , this asymmetry is an intrinsic feature of the 3photon resonant photoionization occurring under the action of the _ chirped _ pulsed laser signal .
moreover , an asymmetry of the ionization line shapes happens to be one of the most pronounced manifestations of the non - monochromaticity of the laser field .
this phenomenon is basically due to the fact that the chirp , even though its relative magnitude amounts to only @xmath317 in our calculations , violates an equivalence of initial ( @xmath318 ) and final ( @xmath319 ) time moments .
indeed , the photoionization probability turns out to be either more efficient or suppressed , depending on whether the time  dependent chirped laser frequency is in the resonance with the 2photon @xmath320 transition or slightly off it .
the process as a whole becomes then somewhat `` time  sensitive '' , although the form of the laser pulse envelope used is time  invariant .
it should be noted , however , that the gauge of the above asymmetry depends strongly on the laser power , as well as particular forms of the laser amplitude and the chirp .
this dependence can be seen in figures [ fig3](a , b ) where higher @xmath161values as compared with figures [ fig2](a , b ) are used .
the former show that the asymmetry discussed is less pronounced as long as higher laser intensities are involved .
two essential characteristics of the line profiles maximum , that are , the shift and the magnitude , are plotted in figure [ fig4](a , b ) as functions of the laser power for various values of the chirp .
in particular , the graph [ fig4](b ) provides information about the dependence of the line shapes maximum value , @xmath321 , on the pulse energy @xmath161 .
it is quite natural to parameterize this relation by the following simple power law : @xmath322 .
the exponent here proves to be a slow  varying function of the pulse energy , such that @xmath323 . according to this condition
, @xmath324 can be estimated as @xmath325.\ ] ]    dependence of the exponent , @xmath324 , defined by the relation : + @xmath326.,width=264,height=264 ]    one should note , for example , an appreciable deviation of the @xmath324values from @xmath327 which must be expected for a 2step 3photon ionization of an atom in the state @xmath37 by a _ weak , monochromatic and spatially homogeneous laser field _ , without the relaxation of the intermediate levels . for @xmath314 , the @xmath161dependence of @xmath324 is displayed in figure  [ fig : power ] .
this graph demonstrates that @xmath328 gradually decreases as the pulse energy increases and tends , for @xmath329 , to an almost constant value , @xmath330 , thus indicating on the presence of saturation in the 3photon transition considered .
separation in muonium ; ( c ) - photoionization line profile measured in , versus laser frequency detuning , @xmath307 .
, width=566,height=264 ]    finally , it should be mentioned that the above consideration has been developed mainly in the attempt to simulate those line profiles of the two
step 3photon ionization probabilities that are supposed to be measured soon in the new @xmath1 experimental investigation of muonium .
essential physical parameters employed in this study are still not fixed completely and are subject to further changes .
it is therefore rather tempting to apply our present technique to appropriate data formerly collected within the framework of the previously employed experimental setup  @xcite .
these are shown in figures  [ fig5](a , b ) , along with the photoionization line profile obtained within the framework of the current theoretical approach , figure  [ fig5](c ) .
the fact that the corresponding experimental line shape ( not shown ) is virtually indistinguishable from the theoretical one provides an additional and independent check of the overall validity of the model developed in this work .
we have developed a simple theoretical scheme intended to describe , to the @xmath19 accuracy , the stepwise 3photon resonant photoionization in hydrogenic systems , induced by the chirped laser field with time
dependent amplitude .
it has been shown that such an accuracy can be achieved within the framework of a simple 3-level model , by taking into account ( i ) the _ ac stark _ shifts and ( ii ) non - zero ionization rates of the levels involved , together with ( iii ) a spatial inhomogeneity of the laser signal and ( iv ) arbitrary @xmath101-dependencies of its intensity ( pulse envelope ) , @xmath331 , and frequency , @xmath31 . the system of equations ( [ ala2])-([cont ] ) or , equivalently , equations ( [ r2s])-([r_cont ] ) is a key point of the method employed .
these are of independent significance as being not specifically restricted to particular states of reference , laser intensities , and chirps , so that the results obtained for @xmath332 and @xmath333states can be generalized on arbitrary @xmath334 levels of hydrogen  like atoms .
excited @xmath335levels , such that @xmath336 , are of particular current interest for ultra - high precision laser spectroscopy , and it should be expected that an adequate interpretation of experimental data should require an accurate theoretical account of a wide spectrum of light  induced effects .
some of these effects are beyond the scope of our present consideration and , in the first instance , comprise appropriate corrections due to inevitable motion of atoms in the media and associated _ second _ order doppler shifts , as being of major importance .
indeed , even for cw lasers where the signals are usually almost unchirped in the laboratory frame , the amplitude and the frequency of the pulsed laser signals become essentially time  dependent in the atomic rest frame , thus leading , as was demonstrated above , to the shifts and distortions of ionization / excitation lines .
the relative contribution of the effects discussed is estimated as @xmath337 with @xmath109 being the atomic velocity in the laboratory frame , and these should be incorporated in the theoretical scheme developed here , along with relativistic / radiative corrections @xmath338 to the operator of the particle  laser field interaction .
this work is currently in progress , and we consider it as a subject of forthcoming publications .
the authors are indebted to g. zu putlitz for his constant support and encouragement .
this work has been strongly motivated and influenced by discussions within the muonium 1s-2s collaboration , particularly with p.e.g .
baird , m.g .
boshier , p.g.h .
sandars , and w.t . toner .
one of us ( v.y . ) , wishes to acknowledge his gratitude to prof .
i.p . grant ( frs ) for his kind advice and to the volkswagen  stiftung and the royal society for financial support .
also , this work has been funded in part by the grants from nato , the royal society and the bundesminister fr bildung und forschung of germany .
99 nez f , plimmer m d and bourzeix s 1993 _ europhys .
lett . _ * 24 * 635 bourzeix s 1996 384 berkeland d , boshier m d and hinds e 1995 2470 schmidt - kaler f , leibfried d , seel s 1995 a * 51 * 2789 udem th , huber a 1997 2646 fee m s , chu s , mills a p 1993 a * 48 * 192 mills jr .
a p 1993 _ hyperfine interactions _ * 76 * 233 maas f e 1994 247 jungmann k 1995 _ physikalisches bltter _ * 51 * 1167 g. zu putlitz 1996 _ hyperfine interactions _ * 103 * 103 boshier m g , hughes v w , jungmann k and zu putlitz g 1996 _ comments atom . mol .
phys . _ * 33 * 17 eikema k s , ubachs u , vassen w and hogervorst w 1997 a * 55 * 1866 ; 1997 1690 giacobino e and cagnac b 1980 _ in progress in optics xvii _ ed e wolf ( amsterdam : north - holland ) pp 86 - 161 zon b a 1974 _ opt .
spectrosc . _ * 36 * 838 zon b a and katsnelson b g 1973 _ sov.phys.jetp _ * 65 * 947 zon b a , manakov n l and rapoport l p 1975 _ opt .
spectrosc . _ * 38 * 6 beausoleil g r and hnsch t w 1986 a * 33 * 1661 kaplan a e 1973 _ sov.phys.jetp _ * 65 * 1416 kaplan a e 1975 _ sov.phys.jetp _ * 68 * 823 lifshitz e m and pitaevskii l p 1974 _ relativistic quantum theory_. course of theoretical physics , part i , vol 4 ( oxford : pergamon ) landau l d and lifshitz e m 1974 _ nonrelativistic quantum theory_. course of theoretical physics , vol 3 ( oxford : pergamon ) bakule p , cornish p s , lane i , baird p e g , toner w t , sandars p g h , liu w , towrie m , _ private communication _ 1970 _ handbook of mathematical functions _ ed m abramowitz and i stegun ( dover , new york ) .
cohen - tannoudji c 1968 _ cargese lectures on physics _
vol 2 ( gordon and breach : n.y . )
levinson m d 1982 _ introduction to nonlinear spectroscopy _ ( academic : new york ) salzman w r 1971 220 vasilenko l s , chebotaev v p and shishaev a v 1970 _ pisma zh .
* 12 * 161 bloembergen n and levensen m d 1976 _ high resolution laser spectroscopy _ ed k shimoda ( springer , berlin ) p 315 sargent iii m , scully m o and lamb jr
. w e 1974 _ laser physics _ ( reading , mass . : addison - wesley ) bassani f , forney j j and quattropani a 1977 1070 yakhontov v and jungmann k 1996 _ _ d * 38 * 141 woodman j 1992 _ ph.d . thesis _ , university of oxford , clarendon laboratory ( unpublished )