 \hat N_{\rm tip})}.\end{aligned}$$ Here, we have introduced different chemical potentials for each subsystem: $\mu_{S}$ for the sample system S and $\mu_{\rm tip}$ for the STM tip. The bias voltage $V$ enters through the difference $\mu_{\rm tip}-\mu_S=eV$. For convenience, we define the chemical potential $\mu_{\rm S}$ of the system S as a reference energy and absorb it into the definition of the single particle energy.

Evaluating the current up quadratic order in the tunneling matrix elements yields $$\begin{aligned}
\label{eq:13}
I(t)  &=&  i \int_{t_0}^t d t'  \left[
\langle \hat V(t')  \hat j_{I}(t)   \rangle_0
-
\langle \hat j_{I}(t) \hat V(t')   \rangle_0
\right]
\non
&&
+ O(t^3_s),
\label{eqn:8}\end{aligned}$$ where $t_s$ is a measure of the order of magnitude of the largest tunneling matrix element $t_{\mu\sigma}$. Note that $\langle \hat V(t') \hat j_{I}(t) \rangle_0=\langle  \hat j_{I}(t) \hat V(t') \rangle_0^*$, ensuring that the current is real.

Substituting the linear expansion of tunneling matrix elements in the displacements of the vibrational modes $\nu$, Eq. , into $\langle \hat V(t')  \hat j_{I}(t)   \rangle_0$ yields

$$\begin{aligned}
\label{eq:Vj}
\langle \hat V(t')  \hat j_{I}(t)    \rangle_0&=&
\sum_{\mu\sigma\mu'\sigma'}  
t_{\mu \sigma}t_{\mu' \sigma'}
\Big \langle (1 +\sum_\nu  \lambda^{\rm tip}_{\mu\nu}\hat X_\nu (t'))(1 +
\sum_{\nu'}  \lambda^{\rm tip}_{\mu'\nu'}\hat X_{\nu'} (t)  )
\left(
d^\dagger_{\mu\sigma}(t')  c_{0\sigma,\rm tip}(
\hat N_{\rm tip})}.\end{aligned}$$ Here, we have introduced different chemical potentials for each subsystem: $ \mu_{S}$ for the sample arrangement S and $ \mu_{\rm tip}$ for the STM peak. The bias voltage $ V$ enters through the difference $ \mu_{\rm tip}-\mu_S = eV$. For appliance, we define the chemical potential $ \mu_{\rm S}$ of the organization S as a reference book energy and steep it into the definition of the individual particle department of energy. 

 Evaluating the current up quadratic order in the tunneling matrix elements yields $ $ \begin{aligned } 
 \label{eq:13 } 
 I(t)   & = &   i \int_{t_0}^t d triiodothyronine'   \left [ 
 \langle \hat V(t')   \hat j_{I}(t)    \rangle_0 
 - 
 \langle \hat j_{I}(t) \hat V(t')    \rangle_0 
 \right ] 
 \non 
 & & 
 + O(t^3_s), 
 \label{eqn:8}\end{aligned}$$ where $ t_s$ is a measure of the decree of magnitude of the largest tunneling matrix component $ t_{\mu\sigma}$. Note that $ \langle \hat V(t') \hat j_{I}(t) \rangle_0=\langle   \hat j_{I}(t) \hat V(t') \rangle_0^*$, ensuring that the current is veridical. 

 Substituting the linear expansion of tunneling matrix elements in the supplanting of the vibrational modes $ \nu$, Eq.  , into $ \langle \hat V(t')   \hat j_{I}(t)    \rangle_0 $ yields 

 $ $ \begin{aligned } 
 \label{eq: Vj } 
 \langle \hat V(t')   \hat j_{I}(t)     \rangle_0&= & 
 \sum_{\mu\sigma\mu'\sigma' }  
 t_{\mu \sigma}t_{\mu' \sigma' } 
 \Big \langle (1 + \sum_\nu   \lambda^{\rm tip}_{\mu\nu}\hat X_\nu (t'))(1 + 
 \sum_{\nu' }   \lambda^{\rm tip}_{\mu'\nu'}\hat X_{\nu' } (t)  ) 
 \left (
 d^\dagger_{\mu\sigma}(t')   c_{0\sigma,\rm tip } (
 \hah N_{\rm tip})}.\end{aligned}$$ Here, we have introdocwd difherent dhemical potentials for each subsystxm: $\my_{S}$ foe the sample system S xnd $\mu_{\rm nip}$ for tye SUM tip. The bias voltage $V$ entera thrmngh the differekce $\mu_{\rm ti[}-\mu_S=eV$. For coneeviznce, we define the chemical potentiaj $\mu_{\rm X}$ lf the system F as w rervrtnce energy and absorb it into ths definption of the singke particle energy.

Evaluatijg tje current up quadgatic order in eye tunneling matrix eltmznts yields $$\begin{aligned}
\label{eq:13}
I(t)  &=&  i \int_{g_0}^t d c'  \left[
\langoe \hah V(t')  \hat j_{I}(v)   \raggle_0
-
\langle \hat j_{I}(t) \hdt V(t')   \rangle_0
\right]
\nok
&&
+ O(t^3_s),
\nabwl{eqn:8}\end{aligned}$$ where $t_s$ is a measure os the ordar of magnitude of rhw larcest tunvwlivg javris elemfnt $t_{\mu\sigma}$. Hote that $\lqngle \hat V(t') \hat j_{I}(u) \rwbgle_0=\langle  \hzt j_{I}(t) \hwt V(t') \rangle_0^*$, ensuring that the current ps rsal.

Substituting the linwar expansion of tunngling matryx elements in the displacements of the vibrationdl mover $\nb$, Eq. , ivro $\langle \hat V(t')  \hat j_{I}(t)   \rangle_0$ yields

$$\begin{ajjgmec}
\label{eq:Vj}
\langje \hat V(t')  \nah k_{Y}(t)    \rangle_0&=&
\som_{\mu\sigma\mh'\sigma'}  
t_{\mu \sigma}t_{\lu' \sigmw'}
\Big \oangle (1 +\stm_\nu  \lambda^{\rm tip}_{\mu\nu}\hat X_\nu (t'))(1 +
\sum_{\nu'}  \lambba^{\rn tip}_{\mu'\nu'}\hat X_{\nu'} (t)  )
\left(
d^\daggzr_{\mu\sibma}(t')  c_{0\sigma,\rm tip}(
\hat N_{\rm tip})}.\end{aligned}$$ Here, we have introduced potentials each subsystem: for the sample for STM tip. The voltage $V$ enters the difference $\mu_{\rm tip}-\mu_S=eV$. For convenience, define the chemical potential $\mu_{\rm S}$ of the system S as a reference and absorb it into the definition of the single particle energy. Evaluating the up order the matrix elements yields $$\begin{aligned} \label{eq:13} I(t) &=& i \int_{t_0}^t d t' \left[ \langle \hat V(t') \hat \rangle_0 - \langle \hat j_{I}(t) \hat V(t') \rangle_0 \non && + O(t^3_s), where $t_s$ is a measure the of magnitude the tunneling element $t_{\mu\sigma}$. Note $\langle \hat V(t') \hat j_{I}(t) \rangle_0=\langle \hat j_{I}(t) \hat V(t') \rangle_0^*$, ensuring that the current is real. the linear tunneling matrix in displacements the vibrational modes , into $\langle \hat V(t') \hat $$\begin{aligned} \label{eq:Vj} \langle \hat V(t') \hat j_{I}(t) \rangle_0&=& t_{\mu \sigma}t_{\mu' \Big \langle (1 +\sum_\nu \lambda^{\rm tip}_{\mu\nu}\hat (t'))(1 + \sum_{\nu'} \lambda^{\rm tip}_{\mu'\nu'}\hat X_{\nu'} (t) ) d^\dagger_{\mu\sigma}(t') c_{0\sigma,\rm tip}(
 \hat N_{\rm tip})}.\end{aligned}$$ Here, we Have introdUced dIffEreNt ChemIcal Potentials for eACh suBsystem: $\mu_{S}$ for the sample SysteM S ANd $\mu_{\RM tIp}$ for The STM tIP. THE BiaS vOlTagE $V$ ENtErs thRouGh the diFference $\mu_{\Rm tIp}-\Mu_S=eV$. For convENiEnce, we defiNe tHe chemical poTenTial $\mu_{\Rm s}$ of THe sysTem s as a rEferenCE energY and absorB iT Into thE DefinitION oF the Single particle eneRGy.

eValuating the cuRrent uP qUAdRATic OrdEr in the tunNeLing mATrix eleMEnTS YIelDS $$\begin{aligned}
\Label{eq:13}
I(t)  &=&  i \iNT_{t_0}^t D t'  \left[
\LaNglE \Hat V(t')  \hAt j_{I}(t)   \RaNGle_0
-
\Langle \hat j_{I}(T) \hat v(t')   \rangle_0
\rIght]
\noN
&&
+ o(t^3_s),
\labeL{Eqn:8}\end{aLigned}$$ WheRe $t_S$ is a MEaSuRe oF tHE orDEr Of mAGniTude of thE lArGest tUnneLING MatrIx eLemeNt $t_{\mu\Sigma}$. Note that $\LanGle \hAT V(t') \Hat j_{I}(T) \rangLe_0=\laNgLe  \hat J_{I}(t) \hat v(t') \ranGlE_0^*$, ensuring that thE curRent is reaL.

SuBsTitUtIng thE Linear ExpAnsIon of tuNneling MAtrIx ELEMeNts in the displacemeNtS OF tHe vibratIonal mODeS $\nU$, eq. , into $\laNgLe \hAt V(t')  \HAT j_{I}(t)   \rAnglE_0$ YiElds

$$\begiN{alignED}
\lAbEl{eq:Vj}
\lAnGle \hat v(t')  \Hat J_{I}(t)    \RanglE_0&=&
\Sum_{\mU\sigma\Mu'\sigma'}  
t_{\Mu \sigMA}t_{\mu' \sigma'}
\Big \laNGle (1 +\sum_\nu  \lambdA^{\Rm TIP}_{\mU\Nu}\haT X_\nU (t'))(1 +
\sum_{\nu'}  \lambDa^{\rm TIp}_{\mu'\Nu'}\haT x_{\nU'} (t)  )
\lEFt(
d^\daGger_{\mU\sIGmA}(T')  c_{0\sigma,\rm tip}(
 \hat N_{\rm tip})}.\end{a ligned}$$Here,  we  ha ve  int rodu ced differentc hemi cal potentials for eac h sub sy s tem: $\ mu_{S }$ fort he s amp le  s yst em Sand $ \mu _{\rm t ip}$ for t heST M tip. The b i as  voltage $ V$enters throu ghthe di ff ere n ce $\ mu_ {\rmtip}-\ m u_S=eV $. For co nv e nience ,  we def i n ethechemical potentia l  $ \ mu_{\rm S}$ of  the s ys t em S  as  areferenceen ergya nd abso r bi t int o  the definiti on of the s i ngl e part ic lee nergy. 

Eva lu a tin g the curre nt u p quadrat ic ord e r in th e  tunnel ing ma tri x e leme n ts  y iel ds $$\ b eg in{ a lig ned}
\la be l{ eq:13 }
I( t ) &=& i\int _{t_0 }^t d t'  \le ft[ 
\la n gle  \hat  V(t' )  \ ha t j_{ I}(t)  \ra ng le_0
-
\langle\hat  j_{I}(t)  \h at  V( t' )   \ r angle_ 0
\ rig ht]
\no n
&&
+O (t^ 3_ s ) , 
\ label{eqn:8}\end{a li g n ed }$$ wher e $t_s $  i sa  measure  o f t he o r d er of  mag n it ude of t he lar g es ttunneli ng  matri xele men t $t_ { \mu\ sigma} $. Notethat$ \langle \hat V ( t') \hat j_{I } (t ) \r a ngle _0= \langle  \h at j _ {I}( t) \ h at  V( t ') \r angle _0 ^ *$ ,  ensuring that thecu rrentis re al.

Substitu ting the l i n e ar expan sion of tunneling matr ix el ements int he displ aceme nts of t he vibrat i o nal mode s $ \nu $,Eq. ,  i nto $\langle\ h at V (t ')  \ha t j _{I}(t)   \ra ngl e_0 $yields

$ $\begin{ al ig ne d} 
\l abel{ e q:Vj}
\l an gle  \ hat  V(t' )   \hat  j_{I }(t)   \ ran gle_0&= & 
\ s u m_{\ mu \s igma \mu '\ sigma '} t_{ \mu \si gma}t_{\m u'\ sigm a' }\Big \l angle (1 +\su m_ \nu  \lamb da ^{\ rm tip } _ {\mu\nu} \hat X_\nu (t'))(1 +
\s u m_{\nu' } \lamb da^{ \rm tip}_ {\m u'\nu' }\h a t X_{\ nu'} ( t)  ) 
\ lef t ( 
d^\d a g ge r_{ \m u\sigma}(t ' )   c _{0\s ig ma,\ rm tip} (
 \hat_N_{\rm tip})}.\end{aligned}$$_Here, we have introduced_different chemical_potentials_for each_subsystem:_$\mu_{S}$ for the_sample system S_and $\mu_{\rm tip}$ for_the STM tip._The_bias voltage $V$ enters through the difference $\mu_{\rm tip}-\mu_S=eV$. For convenience, we define the_chemical_potential $\mu_{\rm_S}$_of_the system S as a_reference energy and absorb it_into the_definition of the single particle energy.

Evaluating the current_up_quadratic order in_the tunneling matrix elements yields $$\begin{aligned}
\label{eq:13}
I(t)  &=& _i \int_{t_0}^t d t'  \left[
\langle_\hat V(t') _\hat_j_{I}(t)_  \rangle_0
-
\langle \hat_j_{I}(t) \hat V(t')   \rangle_0
\right]
\non
&&
+_O(t^3_s),
\label{eqn:8}\end{aligned}$$ where $t_s$ is a measure_of the order of magnitude of the_largest tunneling matrix element $t_{\mu\sigma}$. Note_that $\langle \hat V(t') \hat_j_{I}(t) \rangle_0=\langle_ \hat j_{I}(t) \hat V(t')_\rangle_0^*$, ensuring that_the current_is real.

Substituting the_linear expansion of tunneling matrix elements_in the displacements_of the vibrational modes $\nu$, Eq. ,_into_$\langle \hat V(t')__\hat_j_{I}(t) _ \rangle_0$ yields

$$\begin{aligned}
\label{eq:Vj}
\langle_\hat_V(t') _\hat_j_{I}(t)    \rangle_0&=&
\sum_{\mu\sigma\mu'\sigma'} _
t_{\mu_\sigma}t_{\mu' \sigma'}
\Big \langle (1 +\sum_\nu  \lambda^{\rm_tip}_{\mu\nu}\hat X_\nu (t'))(1 +
\sum_{\nu'}__\lambda^{\rm tip}_{\mu'\nu'}\hat X_{\nu'} (t)_ )
\left(
d^\dagger_{\mu\sigma}(t')  c_{0\sigma,\rm tip}(