 \epsilon_0 \tilde{\chi}[n,N;\omega] \tilde{\chi}[N,m;\omega] }{1+i \epsilon_0 \tilde{\chi}[N,N;\omega]}.\end{aligned}$$ Since there is a driving force only on the first site, we just need to find the susceptibilities to a force on the first site: $$\begin{aligned}
\label{eq:chi_full}
    \tilde{\chi}_{\epsilon_0}[n,1;\omega]
    =
    i^{n}
    \frac{U_{N-n}(\frac{\omega}{2J})-\frac{\epsilon_0}{J} U_{N-1-n}(\frac{\omega}{2J})}
    {J U_N(\frac{\omega}{2J})+(i\frac{\kappa}{2}-\epsilon_0)U_{N-1}(\frac{\omega}{2J})-i\frac{\epsilon_0}{J} \frac{\kappa}{2}U_{N-2}(\frac{\omega}{2J})}\end{aligned}$$

We now compute the steady state total photon number $\bar{n}_{\rm tot}(\epsilon_0)$ when $\epsilon_0 \neq 0$. Recall we are interested in the regime where $\epsilon_0/\kappa \ll 1$ but $e^{A(N-1)} \epsilon_0/\kappa$ is not a priori small. The form of our susceptibilities Eqs.(\[eq:chi\_xx\_full\])-(\[eq:chi\_xp\_full\]) implies that $A$ doesn’t effect the spectrum, but just the residue of the poles as expected from our previous discussion. A non-zero value of $\epsilon_0$ changes both the coherent drive-induced photon number, in addition to drive-independent photons generated from input vacuum fluctuations. The leading order correction to the total photon number when $\epsilon_0 \neq 0$ is therefore $$\begin{aligned}
    \bar{n}_{\rm tot}(\epsilon_0)
    =
    \bar{n}_{\rm tot}(0)
    +
    (
    c
    \frac{\beta^2}{\kappa}
    +
    d
    )
    e^{4A(N-1)}
    (\frac{\epsilon_0}{\kappa})^2 
    +
    \mathcal{O}(

\epsilon_0 \tilde{\chi}[n, N;\omega ] \tilde{\chi}[N, m;\omega ] } { 1+i \epsilon_0 \tilde{\chi}[N, N;\omega]}.\end{aligned}$$ Since there is a driving force only on the inaugural web site, we just need to witness the susceptibility to a force on the first site: $ $ \begin{aligned } 
 \label{eq: chi_full } 
     \tilde{\chi}_{\epsilon_0}[n,1;\omega ] 
     = 
     i^{n } 
     \frac{U_{N - n}(\frac{\omega}{2J})-\frac{\epsilon_0}{J } U_{N-1 - n}(\frac{\omega}{2J }) } 
     { J U_N(\frac{\omega}{2J})+(i\frac{\kappa}{2}-\epsilon_0)U_{N-1}(\frac{\omega}{2J})-i\frac{\epsilon_0}{J } \frac{\kappa}{2}U_{N-2}(\frac{\omega}{2J})}\end{aligned}$$ 

 We nowadays compute the steady state of matter total photon number $ \bar{n}_{\rm tot}(\epsilon_0)$ when $ \epsilon_0 \neq 0$. Recall we are concerned in the regime where $ \epsilon_0/\kappa \ll 1 $ but $ e^{A(N-1) } \epsilon_0/\kappa$ is not a priori small. The form of our susceptibilities Eqs.(\[eq: chi\_xx\_full\])-(\[eq: chi\_xp\_full\ ]) implies that $ A$ does n’t effect the spectrum, but barely the residue of the poles as expected from our previous discussion. A non - zero value of $ \epsilon_0 $ changes both the coherent campaign - induced photon number, in addition to tug - independent photons generated from input signal vacuum fluctuations. The leading order correction to the entire photon number when $ \epsilon_0 \neq 0 $ is therefore $ $ \begin{aligned } 
     \bar{n}_{\rm tot}(\epsilon_0) 
     = 
     \bar{n}_{\rm tot}(0) 
     + 
     (
     c 
     \frac{\beta^2}{\kappa } 
     + 
     d 
    ) 
     e^{4A(N-1) } 
     (\frac{\epsilon_0}{\kappa})^2 
     + 
     \mathcal{O } (
 \epdilon_0 \tilde{\chi}[n,N;\omega] \tiude{\chi}[N,m;\omega] }{1+i \epsiloi_0 \tilde{\dhi}[N,N;\omeea]}.\end{aligned}$$ Since there is e druving force only on the firrt site, wv just newd ti find the susceptinnlitisd to e force on the nirst site: $$\tegin{aligned}
\lateu{ee:chi_full}
    \tilde{\chi}_{\epsilon_0}[n,1;\omega]
    =
    i^{m}
    \frac{U_{N-n}(\frac{\jmegs}{2T})-\frad{\vpwilon_0}{J} U_{N-1-n}(\frac{\omega}{2J})}
    {J U_N(\fdac{\omege}{2J})+(i\frac{\kappa}{2}-\epsolon_0)U_{N-1}(\frac{\omega}{2J})-i\frac{\epsilln_0}{J} \vrac{\kappa}{2}U_{N-2}(\frac{\omeha}{2J})}\end{aligngs}$$

We bow compute ghe steady state total photon number $\bar{n}_{\rm tot}(\epsilov_0)$ wheu $\epsilon_0 \ngq 0$. Ggcall we are integested in the regime fhere $\elsilon_0/\kappa \ll 1$ bnt $e^{Q(N-1)} \epsilon_0/\kappa$ is nov a priori small. The form of mux susceptibilities Eqw.(\[ew:chi\_xf\_fuln\])-(\[eq:cfu\_xp\_wulm\]) mmpmies tjat $A$ doesn’t sffect the wpectrum, but just tne grsidue of ths polef ws expected from our previous discussiot. A non-zero value of $\epsilin_0$ changes both the clherent dwive-induced photon number, in addition to drive-indapendxng pkitons tejerated from input vacuum fluctuations. The lewsimg order correctlon to the total pnohom number when $\gpsilon_0 \nes 0$ is therefore $$\behin{aligged}
    \bar{n}_{\rm uot}(\epxilon_0)
    =
    \bar{n}_{\rm tot}(0)
    +
    (
    c
    \drac{\beta^2}{\kappa}
    +
    d
    )
    e^{4A(N-1)}
    (\ftac{\epsilon_0}{\kappa})^2 
    +
    \matgcal{O}(

\epsilon_0 \tilde{\chi}[n,N;\omega] \tilde{\chi}[N,m;\omega] }{1+i \epsilon_0 \tilde{\chi}[N,N;\omega]}.\end{aligned}$$ Since a force only the first site, the to a force the first site: \label{eq:chi_full} \tilde{\chi}_{\epsilon_0}[n,1;\omega] = i^{n} \frac{U_{N-n}(\frac{\omega}{2J})-\frac{\epsilon_0}{J} U_{N-1-n}(\frac{\omega}{2J})} U_N(\frac{\omega}{2J})+(i\frac{\kappa}{2}-\epsilon_0)U_{N-1}(\frac{\omega}{2J})-i\frac{\epsilon_0}{J} \frac{\kappa}{2}U_{N-2}(\frac{\omega}{2J})}\end{aligned}$$ We now compute the steady state total photon number $\bar{n}_{\rm tot}(\epsilon_0)$ $\epsilon_0 \neq 0$. Recall we are interested in the regime where $\epsilon_0/\kappa \ll but \epsilon_0/\kappa$ not priori small. The form of our susceptibilities Eqs.(\[eq:chi\_xx\_full\])-(\[eq:chi\_xp\_full\]) implies that $A$ doesn’t effect the spectrum, but the residue of the poles as expected from previous discussion. A non-zero of $\epsilon_0$ changes both the drive-induced number, in to photons from input vacuum The leading order correction to the total photon number when $\epsilon_0 \neq 0$ is therefore $$\begin{aligned} \bar{n}_{\rm = \bar{n}_{\rm ( c + ) (\frac{\epsilon_0}{\kappa})^2 + \mathcal{O}(
 \epsilon_0 \tilde{\chi}[n,N;\omega] \tilDe{\chi}[N,m;\omeGa] }{1+i \epSilOn_0 \tIlDe{\chI}[N,N;\oMega]}.\end{aligned}$$ sInce There is a driving force onLy on tHe FIrst SItE, we juSt need tO FiND The SuScEptIbILiTies tO a fOrce on tHe first sitE: $$\beGiN{aligned}
\labeL{Eq:Chi_full}
    \tilDe{\cHi}_{\epsilon_0}[n,1;\omEga]
    =
    I^{n}
    \frac{u_{N-N}(\frAC{\omegA}{2J})-\fRac{\epSilon_0}{J} u_{n-1-n}(\frac{\Omega}{2J})}
    {J U_N(\FrAC{\omega}{2j})+(I\frac{\kaPPA}{2}-\ePsilOn_0)U_{N-1}(\frac{\omega}{2J})-i\frAC{\ePSilon_0}{J} \frac{\kappA}{2}U_{N-2}(\fraC{\oMEgA}{2j})}\End{AliGned}$$

We now cOmPute tHE steady STaTE TOtaL Photon number $\bAr{n}_{\rm tot}(\epsILon_0)$ When $\epSiLon_0 \NEq 0$. RecaLl we aRe INteRested in the RegiMe where $\epSilon_0/\kAPpa \ll 1$ buT $E^{A(N-1)} \epsiLon_0/\kapPa$ iS noT a prIOrI sMalL. THE foRM oF ouR SusCeptibilItIeS Eqs.(\[eQ:chi\_XX\_FULl\])-(\[eq:Chi\_Xp\_fuLl\]) impLies that $A$ doesN’t eFfecT The SpectRum, buT jusT tHe resIdue of The poLeS as expected from Our pRevious diScuSsIon. a nOn-zerO Value oF $\epSilOn_0$ changEs both tHE coHeRENT dRive-induced photon nUmBER, iN additioN to driVE-iNdEPendent pHoTonS genERAted fRom iNPuT vacuum fLuctuaTIoNs. the leadInG order CoRreCtiOn to tHE totAl photOn number When $\ePSilon_0 \neq 0$ is therEFore $$\begin{aligNEd}
    \BAR{n}_{\RM tot}(\EpsIlon_0)
    =
    \bar{n}_{\rm tOt}(0)
    +
    (
    c
    \fRAc{\beTa^2}{\kaPPa}
    +
    D
    )
    e^{4A(n-1)}
    (\Frac{\ePsiloN_0}{\kAPpA})^2 
    +
    \Mathcal{O}(

 \epsilon_0 \tilde{\chi}[n ,N;\omega]  \til de{ \ch i} [N,m ;\om ega] }{1+i \ep s ilon _0 \tilde{\chi}[N,N;\o mega] }. \ end{ a li gned} $$ Sinc e  t h e reis  a  dr iv i ng  forc e o nly onthe firstsit e,  we just nee d  t o find the  su sceptibiliti esto a f or ceo n the  fi rst s ite: $ $ \begin {aligned} 
\ l abel{e q :chi_fu l l }   \tilde{\chi}_{\ep s il o n_0}[n,1;\omeg a]
   =  i^{ n}    \frac{ U_ {N-n} ( \frac{\ o me g a } {2J } )-\frac{\epsi lon_0}{J} U _ {N- 1-n}(\ fr ac{ \ omega} {2J}) }  {J U_N(\fra c{\o mega}{2J} )+(i\f r ac{\kap p a}{2}-\ epsilo n_0 )U_ {N-1 } (\ fr ac{ \o m ega } {2 J}) - i\f rac{\eps il on _0}{J } \f r a c { \kap pa} {2}U _{N-2 }(\frac{\omeg a}{ 2J}) } \en d{ali gned} $$
We  nowcomput e the  s teady state tot al p hoton num ber  $ \ba r{ n}_{\ r m tot} (\e psi lon_0)$  when $ \ eps il o n _ 0\neq 0$. Recall we  a r e  i ntereste d in t h ere g ime wher e$\e psil o n _0/\k appa \l l 1$ but  $e^{A ( N- 1) } \epsi lo n_0/\k ap pa$  is  nota  pri ori sm all. The  form of our suscept i bilities Eqs. ( \[ e q :c h i\_x x\_ full\])-(\[ eq:c h i\_x p\_f u ll \]) impli es th at $A $  doesn’t effect the  s pectru m, bu t just the re sidue of t h e poles as  exp e ct e d from our pre vious  discussio n . A non- zerovalue of  $\epsilo n _ 0$ chang esbot h t hec o he rent drive-in d u cedph oton nu mbe r, in a ddi tio n t o d ri ve-indepe ndent ph ot on sge ner atedf rom inpu tvac uu m f luctu a tions.  Thelead in go rde r corre c ti o n  toth etota l p ho ton n umbe r  wh en $\ep silon_0 \ neq 0$ i sth erefore  $$\begin{ali gn ed}
    \b ar {n} _{\rmt o t}(\epsi lon_0)
    =
    \bar{n } _{\rm t ot} (0)
   + 
    (
   c 
    \ fra c {\beta ^2}{\k appa} 
  +     d   )   e^{4A(N- 1 ) }
   (\ fr ac{\ epsilon _0}{\kappa})^2 
  +    \mathcal{ O}( 

 \epsilon_0_\tilde{\chi}[n,N;\omega] \tilde{\chi}[N,m;\omega]_}{1+i \epsilon_0 \tilde{\chi}[N,N;\omega]}.\end{aligned}$$ Since_there is_a_driving force_only_on the first_site, we just_need to find the_susceptibilities to a_force_on the first site: $$\begin{aligned}
\label{eq:chi_full}
    \tilde{\chi}_{\epsilon_0}[n,1;\omega]
    =
 __ i^{n}
___ \frac{U_{N-n}(\frac{\omega}{2J})-\frac{\epsilon_0}{J} U_{N-1-n}(\frac{\omega}{2J})}
  _ {J U_N(\frac{\omega}{2J})+(i\frac{\kappa}{2}-\epsilon_0)U_{N-1}(\frac{\omega}{2J})-i\frac{\epsilon_0}{J} \frac{\kappa}{2}U_{N-2}(\frac{\omega}{2J})}\end{aligned}$$

We now_compute the_steady state total photon number $\bar{n}_{\rm tot}(\epsilon_0)$ when_$\epsilon_0_\neq 0$. Recall_we are interested in the regime where $\epsilon_0/\kappa \ll_1$ but $e^{A(N-1)} \epsilon_0/\kappa$ is not_a priori small._The_form_of our susceptibilities Eqs.(\[eq:chi\_xx\_full\])-(\[eq:chi\_xp\_full\])_implies that $A$ doesn’t effect the_spectrum, but just the residue of_the poles as expected from our previous_discussion. A non-zero value of $\epsilon_0$_changes both the coherent drive-induced_photon number,_in addition to drive-independent photons_generated from input_vacuum fluctuations._The leading order_correction to the total photon number_when $\epsilon_0 \neq_0$ is therefore $$\begin{aligned}
  __\bar{n}_{\rm tot}(\epsilon_0)
 ___=
 _  \bar{n}_{\rm_tot}(0)
_ __+
    (
 __ c
    \frac{\beta^2}{\kappa}
 _  +
 __ d
  _ )
   _e^{4A(N-1)}
    (\frac{\epsilon_0}{\kappa})^2 
_ _ +
_   \mathcal{O}(
