Technically, this approach, with $N$ modes being taken into account, means following the dynamics of $7N+3$ real variables. According to our experience, for most of our results, the frequency interval $[0, 30\nu]$ with $N=3000$ together with the realistic assumption of $\Omega_n/\omega_0=0.001\sqrt{\tilde{\omega}_n/\omega_0}$ satisfy the requirement mentioned at the end of the previous section: even an exciting pulse as long as few hundred optical cycles is much shorter than the time scale on which the quantized modes observably modify the atom’s dynamics. This means that the number of modes that we take into account has negligible influence on the time evolution of the atom, it is mainly the resolution of the spectra that is determined by the value of $N.$ By requiring 100 frequency values (points on the graphs) in an interval of “length” $\nu,$ the value of $N=3000$ means covering the interval of $[0, 30\nu]$ sufficiently densely. For most of the cases we considered, this is appropriate, since no harmonics of higher order than 30 appear. However, for the strongly detuned excitation shown in Fig. \[resnonresfig\](d), there are roughly ten times more observable harmonics and the value of $N$ had to be increased correspondingly.

As an indication of the correctness of the factorization described above, let us mention that without excitation ($E_0=0$), the initial conditions $\langle N_i\rangle=n\delta_{i0}$ (where $\tilde{\omega}_0=\omega_0$) lead to Rabi flopping between mode $0$ and the atomic system with a frequency very close to $\Omega_0\sqrt{n+1}.$

The time evolution of the photon number expectation values is shown in Figs. \[nexpfig1\] and \[nexpfig2\]. As we can see, the longer the exciting pulse is, the more pronounced peaks can be observed. As one can expect, these peaks appear around integer multiples of the carrier frequency of the exciting pulse. For long enough pulses (see Fig. \[nexpfig2\]) and higher harmonic orders, a significant, qualitative difference can be perceived between the consecutive peaks: the ones at odd multiples of $\nu$ are considerably broader than those corresponding to even harmonics. We will return to this point in the next section.


Technically, this approach, with $ N$ modes being taken into account, mean take after the dynamics of $ 7N+3 $ real variables. accord to our experience, for most of our results, the frequency interval $ [ 0, 30\nu]$ with $ N=3000 $ in concert with the realistic assumption of $ \Omega_n/\omega_0=0.001\sqrt{\tilde{\omega}_n/\omega_0}$ satisfy the requirement mention at the end of the previous incision: even an exciting pulse as long as few hundred ocular cycles is much shorter than the time scale on which the quantize modes observably modify the atom ’s moral force. This means that the number of mode that we take into history has negligible influence on the time evolution of the atom, it is mainly the resolution of the spectrum that is determined by the value of $ N.$ By requiring 100 frequency values (points on the graphs) in an interval of “ length ” $ \nu,$ the value of $ N=3000 $ means covering the interval of $ [ 0, 30\nu]$ sufficiently densely. For most of the cases we considered, this is appropriate, since no harmonics of higher order than 30 appear. However, for the powerfully detuned excitement shown in Fig.   \[resnonresfig\](d), there are roughly ten time more observable harmonics and the value of $ N$ had to be increase correspondingly. 

 As an indication of the correctness of the factorization described above, let us mention that without excitement ($ E_0=0 $), the initial conditions $ \langle N_i\rangle = n\delta_{i0}$ (where $ \tilde{\omega}_0=\omega_0 $) lead to Rabi flopping between mode $ 0 $ and the atomic system with a frequency very close to $ \Omega_0\sqrt{n+1}.$ 

 The time evolution of the photon number expectation values is show in Figs.   \[nexpfig1\ ] and \[nexpfig2\ ]. As we can interpret, the longer the exciting pulse is, the more pronounced peaks can be observed. As one can expect, these bill appear around integer multiples of the carrier frequency of the exciting pulsation. For long adequate pulses (see Fig.   \[nexpfig2\ ]) and higher harmonic orders, a meaning, qualitative difference can be perceived between the consecutive peaks: the ones at curious multiples of $ \nu$ are well broader than those corresponding to even harmonic. We will return to this point in the future part.
Tecjnically, this approach, wlth $N$ modes beiny taken into zccount, oeans following the dynamics od $7N+3$ rtcl variables. Accordivg to our experiebce, hor most of our cssults, bke frseueney interval $[0, 30\nu]$ eith $N=3000$ togather with the rdapistic assumption of $\Omega_n/\omega_0=0.001\sqre{\tilde{\okeha}_n/\omega_0}$ satissy tnq reslivement mentioned at the end of tge prevpous section: even an exciting pulse as long as vew hundred opticap cycles is mucr shorter thav the time scale on whjch the quantized modes observacly mpdify the qtim’s gynamics. Thms meags that the klmber ox modes that we take lnto eccoynt has negligible inhluence on the time gvolution mf the atom, it is mqibly tve rasoljriov or vhe spectga vhat is detsrmined by rhe value of $N.$ By rtquygong 100 frequendy valtef (points on the graphs) in an interval ox “lsngth” $\nu,$ the value of $N=3000$ means covering the ijterval os $[0, 30\nu]$ sufficiently densely. For most of the cases fe coiskdexcq, rhls is appropriate, since no harmonics of highew otdvr than 30 appear. Hjwever, for yhf xjrongly detunea excicztjon shown in Fig. \[rednonressig\](d), rhere are roubhly ten times more observavle harmonicf and the value of $U$ had to be nncreaxed cprrespondingly.

As an indncatioh of the cogrectness uf the factorizagiok dascribed above, let us mentyon that xithobt excitxtiom ($E_0=0$), thq initial fondibhons $\langle N_i\rangpe=n\deptd_{i0}$ (where $\tllde{\omega}_0=\omega_0$) lead to Rabi flo''ing between koge $0$ and the atomlc system with w frequency vety close co $\Omeea_0\sqrt{n+1}.$

The time etolution of ehe photon nukher expectatmon valuef is shoqn in Fkes. \[nexpfig1\] and \[mexpfig2\]. As we can swe, the longer the cxcitkhg pulse is, the motw pronounced prakr cwn bx obsqsved. As one wan dxpdvt, thdse peaks aipexr atound integer multipnes kf the carrier freaucncy of tye excityng pulse. For long enough pulsed (see Fij. \[nexpfog2\]) wnd higher harmonic orders, a sjgnificanh, qmalitative disfercnce can be pexceived between the consecutive peaks: thx ones at odd multiples of $\nu$ are considercbky broader vhan trose corrasponding to even haemonics. We will rtturn to this point in fhe neft seftion.


Technically, this approach, with $N$ modes being account, following the of $7N+3$ real for of our results, frequency interval $[0, with $N=3000$ together with the realistic of $\Omega_n/\omega_0=0.001\sqrt{\tilde{\omega}_n/\omega_0}$ satisfy the requirement mentioned at the end of the previous section: an exciting pulse as long as few hundred optical cycles is much shorter the scale which quantized modes observably modify the atom’s dynamics. This means that the number of modes that we into account has negligible influence on the time of the atom, it mainly the resolution of the that determined by value $N.$ requiring 100 frequency (points on the graphs) in an interval of “length” $\nu,$ the value of $N=3000$ means covering the of $[0, densely. For of cases considered, this is no harmonics of higher order than for the strongly detuned excitation shown in Fig. there are ten times more observable harmonics and value of $N$ had to be increased correspondingly. an indication of the correctness of the factorization described above, let us mention that without the initial conditions $\langle (where $\tilde{\omega}_0=\omega_0$) lead Rabi between $0$ the atomic with a frequency very close to $\Omega_0\sqrt{n+1}.$ The time evolution of photon number expectation values is shown in Figs. \[nexpfig1\] and we see, the longer exciting pulse is, the pronounced can be observed. As expect, peaks multiples the frequency of the exciting For long enough pulses (see \[nexpfig2\]) and higher harmonic can be perceived between the consecutive peaks: the at odd multiples of $\nu$ are considerably than those corresponding to even harmonics. We will return to this point the next
Technically, this approach, wiTh $N$ modes beIng taKen IntO aCcouNt, meAns following thE DynaMics of $7N+3$ real variables. AcCordiNg TO our EXpErienCe, for moST oF OUr rEsUlTs, tHe FReQuencY inTerval $[0, 30\nU]$ with $N=3000$ togeTheR wIth the realisTIc Assumption Of $\OMega_n/\omega_0=0.001\sqRt{\tIlde{\omEgA}_n/\oMEga_0}$ saTisFy the RequirEMent meNtioned at ThE End of tHE previoUS SeCtioN: even an exciting puLSe AS long as few hundRed optIcAL cYCLes Is mUch shorter ThAn the TIme scalE On WHICh tHE quantized modEs observablY ModIfy the AtOm’s DYnamicS. This MeANs tHat the numbeR of mOdes that wE take iNTo accouNT has negLigiblE inFluEnce ON tHe TimE eVOluTIoN of THe aTom, it is mAiNlY the rEsolUTION of tHe sPectRa thaT is determined By tHe vaLUe oF $N.$ By rEquirIng 100 fReQuencY valueS (poinTs On the graphs) in an InteRval of “lenGth” $\Nu,$ The VaLue of $n=3000$ Means cOveRinG the intErval of $[0, 30\NU]$ suFfICIEnTly densely. For most oF tHE CaSes we conSidereD, ThIs IS approprIaTe, sInce NO HarmoNics OF hIgher ordEr than 30 APpEaR. HoweveR, fOr the sTrOngLy dEtuneD ExciTation Shown in FIg. \[resNOnresfig\](d), there ARe roughly ten tIMeS MOrE ObseRvaBle harmonicS and THe vaLue oF $n$ hAd tO Be incReaseD cORrESpondingly.

As an indicAtIon of tHe corRectness of the FactorizatION DescribeD aboVE, lET us mention that WithoUt excitatiON ($E_0=0$), the iniTial cOnditionS $\langle N_i\RANgle=n\delTa_{i0}$ (WheRe $\tIldE{\OMeGa}_0=\omega_0$) lead to rABi flOpPing betWeeN mode $0$ anD thE atOmiC syStEm with a frEquency vErY cLoSe To $\OMega_0\sQRt{n+1}.$

The tiMe EvoLuTioN of thE Photon NumbeR expEcTaTIon Values iS ShOWN in FIgS. \[nExpfIg1\] aNd \[NexpfIg2\]. As WE caN see, the Longer the ExcITing PuLsE is, the mOre pronounced PeAks can be obSeRveD. As one CAN expect, tHese peaks appear around inTEger mulTipLes of The cArrier freQueNcy of tHe eXCiting Pulse. FOr lonG eNouGH PulseS (SEe fig. \[NeXpfig2\]) and hiGHEr hArmonIc OrdeRs, a signIficant, qualitative DIffErence can be peRceIved BETwEen THe COnsEcUTivE PEaks: the ones at odD multiples Of $\NU$ aRe considerABly BrOader thAn those CorreSPonding To even harMonics. We wIlL retURN to This point iN the next Section.


Technically, this approach , with $N$  mode s b ein gtake n in to account, me a ns f ollowing the dynamicsof $7 N+ 3 $ re a lvaria bles. A c co r d ing  t oour  e x pe rienc e,for mos t of our r esu lt s, the frequ e nc y interval  $[ 0, 30\nu]$ w ith  $N=30 00 $ t o gethe r w ith t he rea l isticassumptio no f $\Om e ga_n/\o m e ga _0=0 .001\sqrt{\tilde{ \ om e ga}_n/\omega_0 }$ sat is f yt h e r equ irement me nt ioned at thee nd o f  th e  previous sec tion: evena n e xcitin gpul s e as l ong a sf ewhundred opt ical  cycles i s much shorter than th e time  sc ale  onw hi ch  th eq uan t iz edm ode s observ ab ly  modi fy t h e a tom’ s d ynam ics.This means th atthen umb er of  mode s th at  we t ake in to ac co unt has negligi bleinfluence  on  t heti me ev o lution  of  th e atom,  it ism ain ly t h eresolution of thesp e c tr a that i s dete r mi ne d  by theva lue  of$ N .$ By  req u ir ing 100freque n cy  v alues ( po ints o nthe  gr aphs) in a n inte rval of“leng t h” $\nu,$ thev alue of $N=30 0 0$ m ea n s co ver ing the int erva l  of$[0, 30 \nu ] $ suf ficie nt l yd ensely. For most of  t he cas es we  considered,this is ap p r o priate,sinc e  n o  harmonics ofhighe r order th a n 30 app ear.However,  for thes t rongly d etu ned  ex cit a t io n shown in Fi g .  \[r es nonresf ig\ ](d), t her e a rerou gh ly ten ti mes more  o bs er va ble  harm o nics and  t heva lue  of $ N $ hadto be  inc re as e d c orrespo n di n g ly.
A san i ndi ca tionof t h e c orrectn ess of th e f a ctor iz at ion des cribed above,  l et us ment io n t hat wi t h out exci tation ($E_0=0$), the i n itial c ond ition s $\ langle N_ i\r angle= n\d e lta_{i 0}$ (w here$\ til d e {\ome g a }_ 0=\ om ega_0$) le a d  to  Rabi  f lopp ing bet ween mode $0$ andt heatomic system  wi th a f re que n cy ver yc los e to $\Omega_0\sq rt{n+1}.$
T h etime evolu t ion  o f the p hoton n umber expecta tion valu es is sho wn  inF i gs.  \[nexpfig 1\] and\[nexpfig 2 \]. A s  w e can  se e, the  l ong er th e exci t ing  puls e is,th e more  pron ou nced pea ks can be observed. Asone ca n exp ect , these p eak s  ap pear arou nd i nteger mul tip les  of t hec arrie r fr e qu enc y  of t he e x citing pu l se . F o r  l ong enoughp u l ses  (see  Fi g . \[ne xpfi g2\]) and higherh armonic orders , as i gni fic a nt,qu alitative diff ere nc e can be p er ceived betw een theco n secut ive pe aks: t he ones a to dd mul tipl esof $\nu$are  c o nsidera bl yb roader  tha nthosecorres p ondi n g  to even harmoni cs. W e willr etu rn to  t his poi n t in  the nextsection.


Technically, this_approach, with_$N$ modes being taken_into account,_means_following the_dynamics_of $7N+3$ real_variables. According to_our experience, for most_of our results,_the_frequency interval $[0, 30\nu]$ with $N=3000$ together with the realistic assumption of $\Omega_n/\omega_0=0.001\sqrt{\tilde{\omega}_n/\omega_0}$ satisfy_the_requirement mentioned_at_the_end of the previous section:_even an exciting pulse as_long as_few hundred optical cycles is much shorter than_the_time scale on_which the quantized modes observably modify the atom’s dynamics._This means that the number of_modes that we_take_into_account has negligible influence_on the time evolution of the_atom, it is mainly the resolution_of the spectra that is determined by_the value of $N.$ By requiring_100 frequency values (points on_the graphs)_in an interval of “length”_$\nu,$ the value_of $N=3000$_means covering the_interval of $[0, 30\nu]$ sufficiently densely._For most of_the cases we considered, this is_appropriate,_since no harmonics_of_higher_order than_30 appear. However,_for_the strongly_detuned_excitation shown in Fig. \[resnonresfig\](d), there are_roughly_ten times more observable harmonics and the_value of $N$ had_to_be increased correspondingly.

As an_indication of the correctness of_the factorization described above, let us_mention that_without excitation_($E_0=0$), the initial conditions $\langle N_i\rangle=n\delta_{i0}$ (where $\tilde{\omega}_0=\omega_0$) lead to Rabi_flopping between mode $0$ and the_atomic system with a_frequency very_close_to $\Omega_0\sqrt{n+1}.$

The time_evolution_of the_photon number expectation values is shown in_Figs. \[nexpfig1\] and_\[nexpfig2\]. As we can see, the_longer the exciting pulse_is,_the more pronounced peaks can be_observed. As one can expect, these_peaks appear around integer multiples_of_the_carrier frequency of the exciting_pulse. For long enough pulses (see_Fig. \[nexpfig2\]) and higher_harmonic orders, a significant, qualitative difference can_be_perceived between the consecutive peaks: the_ones_at odd multiples of $\nu$ are_considerably_broader_than those corresponding to even_harmonics. We will return to this_point in the next section.

