 matrices $\sigma_i$ ($\tau_i$) operate in spin (particle-hole) space. Choosing the impurity spin ${\bf S}$ to point along the $z$ direction, the $4\times 4$ Hamiltonian separates into independent $2\times2$ Hamiltonians $${\cal H}_\pm = \xi_{\bf p} \tau_z (V\tau_z \pm J S) \delta({\bf r}) + \Delta \tau_x.
 \label{Hpm}$$ A standard calculation [@SPientka2013] shows that $H_+$ ($H_-$) has one subgap solution $[u_\epsilon({\bf r}),v_\epsilon({\bf r}]$ ($[u_{-\epsilon}({\bf r}),v_{-\epsilon}({\bf r})]$), whose energy we denote by $\epsilon$ ($-\epsilon$). Since the BdG formalism doubles the degrees of freedom, it is sufficient to consider only the solutions of one of these Hamiltonians, say $H_-$. Its subgap (YSR) state starts out at positive energies at small exchange couplings and crosses to negative energies for strong coupling. Keeping only the contribution of the subgap state with Bogoliubov operator $$\gamma_\epsilon = \int d{\bf r} [u^*_\epsilon({\bf r})  \psi_\downarrow({\bf r}) - v^*_\epsilon({\bf r})  \psi_\uparrow^\dagger({\bf r}) ]    ,$$ we can write the electron operators as $$\psi_\downarrow({\bf r}) = u_\epsilon({\bf r}) \gamma_\epsilon+\ldots$$ and $$\psi_\uparrow({\bf r}) = -v^*_\epsilon({\bf r})\gamma_\epsilon^\dagger + \ldots,$$ where the ellipses indicate the contributions of above-gap quasiparticles. The even parity state satisfies $\gamma_\epsilon|{\rm even}\rangle =0$ and the odd-parity state is $|{\rm odd}\rangle = \gamma_\epsilon^\dagger|{\rm even}\rangle$. Note that $\gamma_\epsilon^\dagger$ removes an electron spin of 1/2 from the electron system, so that the spin of the odd parity state is lower by 1/2 compared to the (spinless) even parity state.

In the regime of single-electron tunneling, the differential conductance at $eV=\pm (\Delta_t+\epsilon)$ (for a superconducting tip with gap $\Delta_t$) provides access to the Bogoliubov
matrices $ \sigma_i$ ($ \tau_i$) operate in spin (particle - hole) quad. choose the impurity spin $ { \bf S}$ to point along the $ z$ direction, the $ 4\times 4 $ Hamiltonian separates into autonomous $ 2\times2 $ Hamiltonians $ $ { \cal H}_\pm = \xi_{\bf p } \tau_z (V\tau_z \pm J S) \delta({\bf r }) + \Delta \tau_x. 
  \label{Hpm}$$ A standard calculation [ @SPientka2013 ] shows that $ H_+$ ($ H_-$) has one subgap solution $ [ u_\epsilon({\bf r}),v_\epsilon({\bf r}]$ ($ [ u_{-\epsilon}({\bf r}),v_{-\epsilon}({\bf r})]$), whose department of energy we announce by $ \epsilon$ ($ -\epsilon$). Since the BdG formalism doubles the degrees of freedom, it is sufficient to consider only the solution of one of these Hamiltonians, say $ H_-$. Its subgap (YSR) state start out at positivist energies at humble exchange couplings and crosse to negative energies for strong yoke. Keeping only the contribution of the subgap state of matter with Bogoliubov operator $ $ \gamma_\epsilon = \int d{\bf roentgen } [ u^*_\epsilon({\bf r })   \psi_\downarrow({\bf r }) - v^*_\epsilon({\bf r })   \psi_\uparrow^\dagger({\bf r }) ]    , $ $ we can write the electron operators as $ $ \psi_\downarrow({\bf radius }) = u_\epsilon({\bf r }) \gamma_\epsilon+\ldots$$ and $ $ \psi_\uparrow({\bf r }) = -v^*_\epsilon({\bf r})\gamma_\epsilon^\dagger + \ldots,$$ where the ellipses indicate the contributions of above - gap quasiparticles. The even parity state satisfies $ \gamma_\epsilon|{\rm even}\rangle = 0 $ and the odd - parity state is $ |{\rm odd}\rangle = \gamma_\epsilon^\dagger|{\rm even}\rangle$. Note that $ \gamma_\epsilon^\dagger$ removes an electron spin of 1/2 from the electron system, so that the tailspin of the curious parity department of state is lower by 1/2 compared to the (spinless) even parity state. 

 In the government of single - electron tunneling, the differential conductance at $ eV=\pm (\Delta_t+\epsilon)$ (for a superconducting tip with gap $ \Delta_t$) provides entree to the Bogoliubov
 mahrices $\sigma_i$ ($\tau_i$) operaue in spin (particle-hole) vpace. Dhoosing the impurity spin ${\bf S}$ to plibt aling the $z$ direction, thd $4\times 4$ Jamiltonuan wwparates iifo indeizndenf $2\timzs2$ Hamiltonians $${\gal H}_\pm = \xi_{\tf p} \tau_z (V\tau_s \om J S) \delta({\bf r}) + \Delta \tau_x.
 \label{Hpm}$$ A stancagd calculation [@SPitntha2013] sgows that $H_+$ ($H_-$) has one subgap solutjon $[u_\epvilon({\bf r}),v_\epsolon({\bf r}]$ ($[u_{-\epsilon}({\bf r}),v_{-\epsipon}({\bv r})]$), whose energy wf denote by $\epsyoon$ ($-\epsilon$). Rince the BdG formalisj doubles the degrees of freedoo, it ns sufficiebt to wonsider onoy thv solutions on one ox these Hamiltonians, xay $H_-$. Uts subgap (YSR) state vtarts out at posijive energhea at small exchante coupningv ana crussts vo hegatige xnergies fod strong coypling. Keeping only trv contributioh of tre subgap state with Bogoliubov operator $$\camja_\epsilon = \int d{\bf r} [u^*_\epsilon({\bf r})  \psi_\downartow({\bf r}) - v^*_\qpsilon({\bf r})  \psi_\uparrow^\dagger({\bf r}) ]    ,$$ we can wride thx dleebvon upegators as $$\psi_\downarrow({\bf r}) = u_\epsilon({\bf r}) \gammw_\spxikon+\ldots$$ and $$\pfi_\uparrow({\bf r}) = -f^*_\gpsilon({\bf r})\gamoa_\epsilon^\sagger + \ldots,$$ wherf the ejlipsws indicaue thr contributions of above-gap quasiparticjws. The even parity state satirfiex $\gamka_\epsilon|{\rm even}\rangle =0$ and the odd-parlty state ks $|{\rm odd}\rangle = \gakmd_\epsilon^\dagger|{\rm even}\ranglq$. Note thet $\gakma_\epsiuon^\dsgger$ wemoves an elecbson spin of 1/2 from hhe epewtron systfm, so that the spin of the odd 'erity state ix nowvr by 1/2 coiparec to the (spinjess) even parijy state.

Iu the fegime of aingle-enectron tungeling, the dixverential coiductance at $wV=\pm (\Delta_t+\dosilon)$ (for a sipercondubtnng tip wuth gap $\Delta_t$) proyides zccess to the Bigooiubov
matrices $\sigma_i$ ($\tau_i$) operate in spin (particle-hole) the spin ${\bf to point along 4$ separates into independent Hamiltonians $${\cal H}_\pm \xi_{\bf p} \tau_z (V\tau_z \pm J \delta({\bf r}) + \Delta \tau_x. \label{Hpm}$$ A standard calculation [@SPientka2013] shows that $H_+$ has one subgap solution $[u_\epsilon({\bf r}),v_\epsilon({\bf r}]$ ($[u_{-\epsilon}({\bf r}),v_{-\epsilon}({\bf r})]$), whose energy we by ($-\epsilon$). the formalism doubles the degrees of freedom, it is sufficient to consider only the solutions of one these Hamiltonians, say $H_-$. Its subgap (YSR) state out at positive energies small exchange couplings and crosses negative for strong Keeping the of the subgap with Bogoliubov operator $$\gamma_\epsilon = \int d{\bf r} [u^*_\epsilon({\bf r}) \psi_\downarrow({\bf r}) - v^*_\epsilon({\bf r}) \psi_\uparrow^\dagger({\bf r}) ,$$ we the electron as r}) u_\epsilon({\bf r}) \gamma_\epsilon+\ldots$$ r}) = -v^*_\epsilon({\bf r})\gamma_\epsilon^\dagger + \ldots,$$ indicate the contributions of above-gap quasiparticles. The even state satisfies even}\rangle =0$ and the odd-parity state $|{\rm odd}\rangle = \gamma_\epsilon^\dagger|{\rm even}\rangle$. Note that $\gamma_\epsilon^\dagger$ an electron spin of 1/2 from the electron system, so that the spin of the state is lower by compared to the even state. the of single-electron the differential conductance at $eV=\pm (\Delta_t+\epsilon)$ (for a superconducting tip with $\Delta_t$) provides access to the Bogoliubov
 matrices $\sigma_i$ ($\tau_i$) operate In spin (partIcle-hOle) SpaCe. chooSing The impurity spiN ${\Bf S}$ tO point along the $z$ directiOn, the $4\TiMEs 4$ HaMIlToniaN separaTEs INTo iNdEpEndEnT $2\TiMes2$ HaMilTonians $${\Cal H}_\pm = \xi_{\bf P} \taU_z (v\tau_z \pm J S) \delTA({\bF r}) + \Delta \tau_X.
 \laBel{Hpm}$$ A standArd CalculAtIon [@spientKa2013] sHows tHat $H_+$ ($H_-$) hAS one suBgap solutIoN $[U_\epsilON({\bf r}),v_\epSILoN({\bf r}]$ ($[U_{-\epsilon}({\bf r}),v_{-\epsilON}({\bF R})]$), whose energy we Denote By $\EPsILOn$ ($-\ePsiLon$). Since thE BDG forMAlism doUBlES THe dEGrees of freedoM, it is sufficIEnt To consIdEr oNLy the sOlutiOnS Of oNe of these HaMiltOnians, say $h_-$. Its suBGap (YSR) sTAte starTs out aT poSitIve eNErGiEs aT sMAll EXcHanGE coUplings aNd CrOsses To neGATIVe enErgIes fOr strOng coupling. KeEpiNg onLY thE contRibutIon oF tHe subGap staTe witH BOgoliubov operatOr $$\gaMma_\epsiloN = \inT d{\Bf r} [U^*_\ePsiloN({\Bf r})  \psi_\DowNarRow({\bf r}) - v^*_\Epsilon({\BF r})  \pSi_\UPARrOw^\dagger({\bf r}) ]    ,$$ we can wrItE THe Electron OperatORs As $$\PSi_\downarRoW({\bf R}) = u_\epSILon({\bf R}) \gamMA_\ePsilon+\ldOts$$ and $$\PSi_\UpArrow({\bf R}) = -v^*_\EpsiloN({\bF r})\gAmmA_\epsiLOn^\daGger + \ldOts,$$ where The elLIpses indicate tHE contributionS Of ABOvE-Gap qUasIparticles. THe evEN parIty sTAtE saTIsfieS $\gammA_\ePSiLOn|{\rm even}\rangle =0$ and thE oDd-pariTy staTe is $|{\rm odd}\rangLe = \gamma_\epsILON^\dagger|{\rM eveN}\RaNGle$. Note that $\gamMa_\epsIlon^\dagger$ REmoves an ElectRon spin oF 1/2 from the eLECtron sysTem, So tHat The SPIn Of the odd paritY STate Is Lower by 1/2 ComPared to The (SpiNleSs) eVeN parity stAte.

In the ReGiMe Of SinGle-elECtron tunNeLinG, tHe dIfferENtial cOnducTancE aT $ev=\Pm (\DElta_t+\epSIlON)$ (For a SuPeRconDucTiNg tip With GAp $\DElta_t$) prOvides accEss TO the boGoLiubov
 matrices $\sigma_i$ ($\ta u_i$) oper ate i n s pin  ( part icle -hole) space.C hoos ing the impurity spin${\bf  S } $ to po int a long th e  $ z $  di re ct ion ,t he  $4\t ime s 4$ Ha miltoniansep ar ates into in d ep endent $2\ tim es2$ Hamilto nia ns $${ \c alH }_\pm  =\xi_{ \bf p} \tau_z  (V\tau_z  \ p m J S) \delta( { \ bf  r})  + \Delta \tau_x.  \ l abel{Hpm}$$ Astanda rd ca l c ula tio n [@SPient ka 2013] shows t h at $ H _+$ ($H_-$) has o ne subgap s o lut ion $[ u_ \ep s ilon({ \bf r }) , v_\ epsilon({\b f r} ]$ ($[u_{ -\epsi l on}({\b f  r}),v_ {-\eps ilo n}( {\bf r} )] $),  w h ose en erg y  we  denoteby  $ \epsi lon$ ( $ - \eps ilo n$).  Sinc e the BdG for mal ismd oub les t he de gree sof fr eedom,  it i ssufficient to c onsi der onlythe  s olu ti ons o f  one o f t hes e Hamil tonians ,  sa y$ H _ -$ . Its subgap (YSR)  s t a te  startsout at po si t ive ener gi esat s m a ll ex chan g ecoupling s andc ro ss es to n eg ativeen erg ies  fors tron g coup ling. Ke eping only the contr i bution of the su b g ap stat e w ith Bogoliu bovo pera tor$ $\ gam m a_\ep silon  = \i n t d{\bf r} [u^*_\ep si lon({\ bf r} )  \psi_\down arrow({\bf r } ) - v^*_ \eps i lo n ({\bf r})  \ps i_\up arrow^\dag g er({\bfr}) ]     ,$$we can wr i t e the el ect ron  op era t o rs  as $$\psi_\d o w narr ow ({\bf r })= u_\ep sil on( {\b f r })  \gamma_\ epsilon+ \l do ts $$  an d $$\ p si_\upar ro w({ \b f r }) =- v^*_\e psilo n({\ bf  r } )\g amma_\e p si l o n^\d ag ge r +\ld ot s,$$wher e  th e ellip ses indic ate theco nt ributio ns of above-g ap  quasipart ic les . Thee v en parit y state satisfies $\gam m a_\epsi lon |{\rm  eve n}\rangle  =0 $ andthe odd-pa rity s tateis  $| { \ rm od d } \r ang le  = \gamma_ \ e psi lon^\ da gger |{\rm e ven}\rangle$. Note tha t $\gamma_\ep sil on^\ d a gg er$ re m ove sa n e l e ctron spin of 1 /2 from th ee le ctron syst e m,so  that t he spin  of t h e odd p arity sta te is low er  by1 / 2 c ompared to  the (sp inless) e v en pa r it y sta te. 

In t he  re gimeof sin g le- elect ron tu nn eling,  thedi fferenti al conductance at $eV=\ pm (\D elta_ t+\ epsilon)$  (f o r a  supercon duct ing tip wi thgap  $\De lta _ t$) p rovi d es  ac c ess t o th e  Bogoliub o v
 matrices_$\sigma_i$ ($\tau_i$)_operate in spin (particle-hole)_space. Choosing_the_impurity spin_${\bf_S}$ to point_along the $z$_direction, the $4\times 4$_Hamiltonian separates into_independent_$2\times2$ Hamiltonians $${\cal H}_\pm = \xi_{\bf p} \tau_z (V\tau_z \pm J S) \delta({\bf r})_+_\Delta \tau_x.
_\label{Hpm}$$_A_standard calculation [@SPientka2013] shows that_$H_+$ ($H_-$) has one subgap_solution $[u_\epsilon({\bf_r}),v_\epsilon({\bf r}]$ ($[u_{-\epsilon}({\bf r}),v_{-\epsilon}({\bf r})]$), whose energy we_denote_by $\epsilon$ ($-\epsilon$)._Since the BdG formalism doubles the degrees of freedom,_it is sufficient to consider only_the solutions of_one_of_these Hamiltonians, say $H_-$._Its subgap (YSR) state starts out_at positive energies at small exchange_couplings and crosses to negative energies for_strong coupling. Keeping only the contribution_of the subgap state with_Bogoliubov operator_$$\gamma_\epsilon = \int d{\bf r}_[u^*_\epsilon({\bf r}) _\psi_\downarrow({\bf r})_- v^*_\epsilon({\bf r})_ \psi_\uparrow^\dagger({\bf r}) ]  _ ,$$ we_can write the electron operators as_$$\psi_\downarrow({\bf_r}) = u_\epsilon({\bf_r})_\gamma_\epsilon+\ldots$$_and $$\psi_\uparrow({\bf_r}) = -v^*_\epsilon({\bf_r})\gamma_\epsilon^\dagger_+ \ldots,$$_where_the ellipses indicate the contributions of_above-gap_quasiparticles. The even parity state satisfies $\gamma_\epsilon|{\rm_even}\rangle =0$ and the_odd-parity_state is $|{\rm odd}\rangle_= \gamma_\epsilon^\dagger|{\rm even}\rangle$. Note that_$\gamma_\epsilon^\dagger$ removes an electron spin of_1/2 from_the electron_system, so that the spin of the odd parity state is_lower by 1/2 compared to the_(spinless) even parity state.

In_the regime_of_single-electron tunneling, the_differential_conductance at_$eV=\pm (\Delta_t+\epsilon)$ (for a superconducting tip with_gap $\Delta_t$)_provides access to the Bogoliubov