 As one can check easily, one cannot obtain an exact closed set of dynamical equations, since the time derivatives of two-operator products involve three-operator products, whose time derivatives contain terms being the products of four operators, etc. Because of this hierarchy, we have to restrict ourselves to expectation values and – at some point – introduce a factorization that is based on physical considerations.

The simplest Heisenberg equations of motion read: $$\begin{aligned}
\frac{d}{dt}\sigma_x&=&\omega_0 \sigma_y, \label{der1}\\
\frac{d}{dt}\sigma_y&=&-\omega_0 \sigma_x+\sigma_z \left[\Omega(t)+\sum_n \Omega_n \left(a_n+a_n^{\dagger}\right)\right], \label{der2}\\
\frac{d}{dt}\sigma_z&=&-\sigma_y \left[ \Omega(t) + \sum_n \Omega_n (a_n+a_n^{\dagger})\right], \label{der3}\\
\frac{d}{dt}a_n&=&-i\left(\tilde{\omega}_n a_n+\frac{\Omega_n}{2}\sigma_x\right),\label{der4} \\
\frac{d}{dt}N_n&=&\frac{d}{dt} a_n^{\dagger}a_n= i\frac{\Omega_n}{2}\left(a_n-a_n^{\dagger}\right)\sigma_x \label{der5}.
\label{Heis}\end{aligned}$$ One can assume already at this point that the expectation values of the products of atomic an field operators factorize (e.g., $\langle\sigma_z a_n\rangle=\langle\sigma_z\rangle \langle a_n\rangle$), but in this case the influence of the photon number operators on the atomic system’s time evolution – which should be weak – is exactly zero. Additionally, since the photon annihilation operators’ dynamics are governed by a “classical” quantity (the expectation value of the dipole moment), the resulting photon statistics is necessarily Poissonian (i.e., the high harmonic modes are in coherent states).

![Central panel: the time evolution of the expectation values of the photon number operators $\langle N_n \rangle$ as a function of the frequency of the modes
As one can check easily, one cannot obtain an exact shut hardening of dynamical equations, since the meter derivative of two - operator products necessitate three - hustler products, whose time derivatives contain terms being the product of four operators, etc. Because of this hierarchy, we have to qualify ourselves to expectation values and – at some point – introduce a factorization that is based on physical considerations. 

 The simplest Heisenberg equation of motion read: $ $ \begin{aligned } 
 \frac{d}{dt}\sigma_x&=&\omega_0 \sigma_y, \label{der1}\\ 
 \frac{d}{dt}\sigma_y&=&-\omega_0 \sigma_x+\sigma_z \left[\Omega(t)+\sum_n \Omega_n \left(a_n+a_n^{\dagger}\right)\right ], \label{der2}\\ 
 \frac{d}{dt}\sigma_z&=&-\sigma_y \left [ \Omega(t) + \sum_n \Omega_n (a_n+a_n^{\dagger})\right ], \label{der3}\\ 
 \frac{d}{dt}a_n&=&-i\left(\tilde{\omega}_n a_n+\frac{\Omega_n}{2}\sigma_x\right),\label{der4 } \\ 
 \frac{d}{dt}N_n&=&\frac{d}{dt } a_n^{\dagger}a_n= i\frac{\Omega_n}{2}\left(a_n - a_n^{\dagger}\right)\sigma_x \label{der5 }. 
 \label{Heis}\end{aligned}$$ One can assume already at this point that the arithmetic mean values of the products of atomic an field operator factorize (e.g., $ \langle\sigma_z a_n\rangle=\langle\sigma_z\rangle \langle a_n\rangle$), but in this case the influence of the photon number operators on the atomic arrangement ’s time evolution – which should be weak – is exactly zero. Additionally, since the photon annihilation operators ’ dynamics are governed by a “ classical ” quantity (the expectation value of the dipole moment), the resulting photon statistics is necessarily Poissonian (i.e., the high harmonic modes are in coherent state). 

! [ Central panel: the prison term development of the arithmetic mean values of the photon number operator $ \langle N_n \rangle$ as a affair of the frequency of the modes
 As one can check easily, ont cannot obtain au exact closes set of dynamical equations, since tie tume dtgivatives of two-operagor produbts involce tiree-operator provhcts, whose tijc dernvetives contain jerms being dhe products ox wobr operators, etc. Because of this hiewarchy, ee have to restryct ptrsemnew to expectation values and – zt some point – introcuce a factorization that ls bwsed on physical clnsideratiobs.

Thq simplest Heksenberg equations of jotion read: $$\begin{aligned}
\frac{d}{dt}\rigma_r&=&\omega_0 \sigmq_y, \lahgl{der1}\\
\frac{d}{dt}\wigma_j&=&-\omega_0 \sigma_x+\sigma_z \laft[\Omegs(t)+\sum_n \Omega_n \kefv(a_n+a_b^{\dagger}\right)\right], \labxl{der2}\\
\frac{d}{dt}\sigma_z&=&-\sidma_y \left[ \Ojega(t) + \sum_n \Omega_b (q_n+a_n^{\ddgges})\rigfr], \lxbem{dxr3}\\
\fdac{d}{dt}w_n&=&-i\meft(\tilde{\ojega}_n a_n+\frax{\Omega_n}{2}\sigma_x\right),\lsbqo{der4} \\
\frac{d}{dt}N_h&=&\frac{d}{qt} a_n^{\dagger}a_n= i\frac{\Omega_n}{2}\left(a_n-a_n^{\dagger}\richt)\aigma_x \label{der5}.
\label{Heiw}\end{aligned}$$ One can adsume alrqady at this point that the expectation values of the 'ruduebf if atomic an field operators factorize (e.g., $\langlq\aibms_z a_n\rangle=\lannle\sigma_z\rangle \lamgpe w_n\rangle$), but kn this czse the influence lf the khoton number okeratprs on the atomic system’s tume evolutiou – qhich should be weck – is exactuy zgro. Adcitionally, since the phuton annihilatiln operatkfs’ dynamics are eovvrneg by a “classical” quantity (ehe expecvatiou value uf tne dipjle moment), the vasulting photon stwtistncs iv necessarlly Poissonian (i.e., the high harmonic modes are it cmherent ftatex).

![Central panej: the time evokution jf thd expectatpon valuev of the phjton number o[grators $\langlx N_n \rangje$ aw a dunctiov of the frequemcy of thv iodew
As one can check easily, one cannot exact set of equations, since the involve products, whose time contain terms being products of four operators, etc. Because this hierarchy, we have to restrict ourselves to expectation values and – at point – introduce a factorization that is based on physical considerations. The simplest equations motion $$\begin{aligned} \sigma_y, \label{der1}\\ \frac{d}{dt}\sigma_y&=&-\omega_0 \sigma_x+\sigma_z \left[\Omega(t)+\sum_n \Omega_n \left(a_n+a_n^{\dagger}\right)\right], \label{der2}\\ \frac{d}{dt}\sigma_z&=&-\sigma_y \left[ \Omega(t) + \sum_n \Omega_n (a_n+a_n^{\dagger})\right], \label{der3}\\ a_n+\frac{\Omega_n}{2}\sigma_x\right),\label{der4} \\ \frac{d}{dt}N_n&=&\frac{d}{dt} a_n^{\dagger}a_n= i\frac{\Omega_n}{2}\left(a_n-a_n^{\dagger}\right)\sigma_x \label{der5}. \label{Heis}\end{aligned}$$ One assume already at this that the expectation values of products atomic an operators (e.g., a_n\rangle=\langle\sigma_z\rangle \langle a_n\rangle$), in this case the influence of the photon number operators on the atomic system’s time evolution – should be is exactly Additionally, the annihilation operators’ dynamics by a “classical” quantity (the expectation dipole moment), the resulting photon statistics is necessarily (i.e., the harmonic modes are in coherent states). panel: the time evolution of the expectation values the photon number operators $\langle N_n \rangle$ as a function of the frequency of the
 As one can check easily, one canNot obtain aN exacT clOseD sEt of DynaMical equations, SInce The time derivatives of twO-operAtOR proDUcTs invOlve thrEE-oPERatOr PrOduCtS, WhOse tiMe dErivatiVes contain TerMs Being the prodUCtS of four opeRatOrs, etc. BecausE of This hiErArcHY, we haVe tO restRict ouRSelves To expectaTiON valueS And – at soME PoInt – iNtroduce a factorizATiON that is based on PhysicAl COnSIDerAtiOns.

The simpLeSt HeiSEnberg eQUaTIONs oF Motion read: $$\begIn{aligned}
\frAC{d}{dT}\sigma_X&=&\oMegA_0 \Sigma_y, \Label{DeR1}\\
\FraC{d}{dt}\sigma_y&=&-\oMega_0 \Sigma_x+\sigMa_z \lefT[\omega(t)+\sUM_n \Omega_N \left(a_N+a_n^{\DagGer}\rIGhT)\rIghT], \lABel{DEr2}\\
\FraC{D}{dt}\Sigma_z&=&-\siGmA_y \Left[ \OMega(T) + \SUM_N \OmeGa_n (A_n+a_n^{\DaggeR})\right], \label{deR3}\\
\frAc{d}{dT}A_n&=&-i\Left(\tIlde{\oMega}_N a_N+\frac{\omega_n}{2}\Sigma_X\rIght),\label{der4} \\
\fraC{d}{dt}n_n&=&\frac{d}{dt} A_n^{\dAgGer}A_n= I\frac{\oMega_n}{2}\lEft(A_n-a_N^{\dagger}\Right)\siGMa_x \LaBEL{DeR5}.
\label{Heis}\end{alignEd}$$ oNE cAn assume AlreadY At ThIS point thAt The ExpeCTAtion ValuES oF the prodUcts of AToMiC an fielD oPeratoRs FacTorIze (e.g., $\LAnglE\sigma_Z a_n\ranglE=\langLE\sigma_z\rangle \lANgle a_n\rangle$), bUT iN THiS Case The Influence of The pHOton NumbER oPerATors oN the aToMIc SYstem’s time evolution – WhIch shoUld be Weak – is exactly Zero. AdditiONALly, since The pHOtON annihilation oPeratOrs’ dynamicS Are goverNed by A “classicAl” quantitY (THe expectAtiOn vAluE of THE dIpole moment), thE REsulTiNg photoN stAtisticS is NecEssAriLy poissoniaN (i.e., the hiGh HaRmOnIc mOdes aRE in coherEnT stAtEs).

![CEntraL Panel: tHe timE evoLuTiON of The expeCTaTIOn vaLuEs Of thE phOtOn numBer oPEraTors $\lanGle N_n \rangLe$ aS A funCtIoN of the fRequency of the MoDes
 As one can check easily,one cannot  obta inanex actclos ed set of dyna m ical  equations, since thetimede r ivat i ve s oftwo-ope r at o r  pr od uc tsin v ol ve th ree -operat or product s,wh ose time der i va tives cont ain  terms being  th e prod uc tso f fou r o perat ors, e t c. Bec ause of t hi s  hiera r chy, we h av e to  restrict ourselv e st o expectationvalues  a n d– atsom e point –in trodu c e a fac t or i z a tio n  that is base d on physic a l c onside ra tio n s.

Th e sim pl e stHeisenbergequa tions ofmotion read: $ $ \begin{ aligne d}\fr ac{d } {d t} \si gm a _x& = &\ ome g a_0  \sigma_ y,  \ label {der 1 } \ \ 
\fr ac{ d}{d t}\si gma_y&=&-\ome ga_ 0 \s i gma _x+\s igma_ z \l ef t[\Om ega(t) +\sum _n  \Omega_n \left (a_n +a_n^{\da gge r} \ri gh t)\ri g ht], \ lab el{ der2}\\ 
\frac{ d }{d t} \ s i gm a_z&=&-\sigma_y \l ef t [  \ Omega(t)  + \su m _n  \ O mega_n ( a_ n+a _n^{ \ d agger })\r i gh t], \lab el{der 3 }\ \\frac{d }{ dt}a_n &= &-i \le ft(\t i lde{ \omega }_n a_n+ \frac { \Omega_n}{2}\s i gma_x\right), \ la b e l{ d er4}  \\ 
\frac{d}{d t}N_ n &=&\ frac { d} {dt }  a_n^ {\dag ge r }a _ n= i\frac{\Omega_n} {2 }\left (a_n- a_n^{\dagger} \right)\si g m a _x \labe l{de r 5} . 
\label{Heis}\ end{a ligned}$$O ne can a ssume  already  at thisp o int that  th e e xpe cta t i on  values of th e prod uc ts of a tom ic an f iel d o per ato rs  factoriz e (e.g.,  $ \l an gl e\s igma_ z  a_n\ran gl e=\ la ngl e\sig m a_z\ra ngle\lan gl ea _n\ rangle$ ) ,b u t in  t hi s ca seth e inf luen c e o f the p hoton num ber oper at or s on th e atomic syst em ’s time ev ol uti on – w h i ch shoul d be weak – is exactlyz ero. Ad dit ional ly,since the  ph oton a nni h ilatio n oper ators ’dyn a m ics a r e  g ove rn ed by a “c l a ssi cal”qu anti ty (the  expectation value ofthe dipole mo men t),t h eres u lt i ngph o ton s tatistics is ne cessarilyPo i ss onian (i.e . , t he  high h armonic  mode s  are in  coherent  states). 
![Ce n t ral  panel: th e time e volutiono f the ex pecta tio n valu es  of  thephoton num ber o perato rs  $\lan gle N _n  \rangle $ as a function of thefreque ncy o f t he modes
 As_one can_check easily, one cannot_obtain an_exact_closed set_of_dynamical equations, since_the time derivatives_of two-operator products involve_three-operator products, whose_time_derivatives contain terms being the products of four operators, etc. Because of this hierarchy,_we_have to_restrict_ourselves_to expectation values and –_at some point – introduce_a factorization_that is based on physical considerations.

The simplest Heisenberg_equations_of motion read:_$$\begin{aligned}
\frac{d}{dt}\sigma_x&=&\omega_0 \sigma_y, \label{der1}\\
\frac{d}{dt}\sigma_y&=&-\omega_0 \sigma_x+\sigma_z \left[\Omega(t)+\sum_n \Omega_n \left(a_n+a_n^{\dagger}\right)\right], \label{der2}\\
\frac{d}{dt}\sigma_z&=&-\sigma_y \left[_\Omega(t) + \sum_n \Omega_n (a_n+a_n^{\dagger})\right], \label{der3}\\
\frac{d}{dt}a_n&=&-i\left(\tilde{\omega}_n_a_n+\frac{\Omega_n}{2}\sigma_x\right),\label{der4} \\
\frac{d}{dt}N_n&=&\frac{d}{dt} a_n^{\dagger}a_n=_i\frac{\Omega_n}{2}\left(a_n-a_n^{\dagger}\right)\sigma_x_\label{der5}.
\label{Heis}\end{aligned}$$_One can assume already_at this point that the expectation_values of the products of atomic_an field operators factorize (e.g., $\langle\sigma_z a_n\rangle=\langle\sigma_z\rangle_\langle a_n\rangle$), but in this case_the influence of the photon_number operators_on the atomic system’s time_evolution – which_should be_weak – is_exactly zero. Additionally, since the photon_annihilation operators’ dynamics_are governed by a “classical” quantity_(the_expectation value of_the_dipole_moment), the_resulting photon statistics_is_necessarily Poissonian_(i.e.,_the high harmonic modes are in_coherent_states).

![Central panel: the time evolution of the_expectation values of the_photon_number operators $\langle N_n_\rangle$ as a function of_the frequency of the modes