 $\epsilon_j$ for QD $j=1,2$, one arrives at $$\label{HDots}
 H_{\rm d}=\sum_{j=1,2} \epsilon_j d_j^{\dagger} d^{}_j,$$ see Ref. [@Karzig2017] for details. The energies $\epsilon_j$ can be controlled by variation of the gate voltage parameter $n_g$. Without loss of generality, we take $\epsilon_2>\epsilon_1$ throughout, where both energies should satisfy $|\epsilon_j|\ll {\rm min}\{E_C,\Delta\}$. In addition, we employ a time-dependent electromagnetic driving field which can pump single electrons between the two QDs via a tunnel link. To that end, a suitable AC voltage can be applied to a gate electrode located near this link. The respective Hamiltonian contribution is given by [@Platero2004] $$\label{Hdriv}
H_{{\rm drive}}(t) = w(t) d_1^{\dagger} d_2^{}+{\rm h.c.},\quad
w(t) = t_{12}+2A \cos\left(\omega_0 t\right),$$ where $\omega_0$ denotes the drive frequency and $A$ the drive amplitude. In what follows, we assume that the static contribution vanishes, $t_{12}=0$, because a small coupling $t_{12}\ne 0$ will not affect the dissipator in the Lindblad equation, see Eq.  below, and thus does not change the physics in a qualitative manner.

In this work, we consider the Coulomb valley regime where the total charge on the box is fixed by the charging term in Eq.  on time scales $\delta t>1/E_C$ [@Romito2014]. The total particle number on the QDs, $N_{\rm d}=\sum_j d_j^\dagger d_j^{}$, is therefore also conserved on such time scales. For even $N_{\rm d}\in\{0,2\}$, the inter-QD dynamics is effectively frozen out. We here mainly focus on the case $N_{\rm d}=1$, where the pair of QDs forms a spin-1/2 degree of freedom corresponding to Pauli operators $\tau_{x,y,z}$ with $\tau_\pm=(\tau_x\pm i \tau_y)/2$, $$\tau_+ = \tau_-^\dagger = d_1
$ \epsilon_j$ for QD $ j=1,2 $, one arrives at $ $ \label{HDots } 
  H_{\rm d}=\sum_{j=1,2 } \epsilon_j d_j^{\dagger } d^{}_j,$$ see Ref.   [ @Karzig2017 ] for details. The energies $ \epsilon_j$ can be operate by mutant of the gate voltage parameter $ n_g$. Without loss of generalization, we take $ \epsilon_2>\epsilon_1 $ throughout, where both energies should satisfy $ |\epsilon_j|\ll { \rm min}\{E_C,\Delta\}$. In accession, we employ a time - pendent electromagnetic driving field which can pump individual electrons between the two QDs via a tunnel liaison. To that end, a suitable AC voltage can be applied to a gate electrode located near this connection. The respective Hamiltonian contribution is given by [ @Platero2004 ] $ $ \label{Hdriv } 
 H_{{\rm drive}}(t) = w(t) d_1^{\dagger } d_2^{}+{\rm h.c.},\quad 
 w(t) = t_{12}+2A \cos\left(\omega_0 t\right),$$ where $ \omega_0 $ denotes the drive frequency and $ A$ the drive amplitude. In what be, we assume that the static contribution vanishes, $ t_{12}=0 $, because a small coupling $ t_{12}\ne 0 $ will not involve the dissipator in the Lindblad equation, see Eq.    below, and thus does not change the purgative in a qualitative manner. 

 In this work, we consider the Coulomb valley regime where the total charge on the box is fix by the charging term in Eq.    on time scales $ \delta t>1 / E_C$ [ @Romito2014 ]. The total particle number on the QDs, $ N_{\rm d}=\sum_j d_j^\dagger d_j^{}$, is therefore also conserved on such clock time scales. For even $ N_{\rm d}\in\{0,2\}$, the inter - QD dynamics is efficaciously freeze out. We here mainly focus on the case $ N_{\rm d}=1 $, where the couple of QDs forms a spin-1/2 academic degree of freedom corresponding to Pauli operators $ \tau_{x, y, z}$ with $ \tau_\pm=(\tau_x\pm i \tau_y)/2 $, $ $ \tau_+ = \tau_-^\dagger = d_1
 $\epdilon_j$ for QD $j=1,2$, one arriyes at $$\label{HDots}
 H_{\rm d}=\vum_{j=1,2} \elsilon_j a_j^{\dagger} d^{}_j,$$ see Ref. [@Karzig2017] foc derails. The energies $\epsilon_j$ can be clntrollee by cariation of the gate vombage 'acameter $n_g$. Withput loss ox generality, wa gane $\epsilon_2>\epsilon_1$ throughout, where foth enrrhies should sajisfy $|\qpsimon_j|\ll {\rm min}\{E_C,\Delta\}$. In addition, ws emploj a time-dependent electromagnetic driving fleld which can pump sijgle electrins fwtween the tdo QDs via a tunnel lihk. To that end, a suitable AC voutage can be apkljef to a gate xlectrjde located kvar thiv link. Yhe respective Hakilronian contribution iv given by [@Platero2004] $$\label{Hdrhv}
K_{{\rm drive}}(t) = w(t) d_1^{\daggee} e_2^{}+{\rm h.w.},\quag
w(t) = r_{12}+2A \zos\meht(\ojega_0 t\gigit),$$ where $\omsga_0$ denotes the drive frequencu wbd $A$ the drivs amplyttde. In what follows, we assume that the vtafic contribution vanishws, $t_{12}=0$, because a small foupling $e_{12}\ne 0$ will not affect the dissipator in the Lindbldd eqnagiou, see Dw.  helow, and thus does not change the physics in z auslitative manncr.

In this work, we vojsoqer the Couloob valley regime where the hotal crarge on the bjx ix fixed by the charging tern in Eq.  on tpme wcales $\delta t>1/E_C$ [@Rlmito2014]. The tutal partocle number on the QDs, $U_{\rm d}=\shm_j d_j^\daggeg d_j^{}$, is tgdrefore also conrerned mn such ukme scales. For evqn $N_{\rm d}\ii\{0,2\}$, the inter-QA dymamics is effectlvely frozen out. We here mainpy focus on hhe case $N_{\rm d}=1$, where the pair oh QDs forms a s[in-1/2 degree jf frcedom corresponqing to Pauli pperatoxs $\tau_{b,y,z}$ with $\tzu_\pm=(\tau_e\pm i \tau_y)/2$, $$\twu_+ = \tau_-^\dagger = d_1
$\epsilon_j$ for QD $j=1,2$, one arrives at d}=\sum_{j=1,2} d_j^{\dagger} d^{}_j,$$ Ref. [@Karzig2017] for be by variation of gate voltage parameter Without loss of generality, we take throughout, where both energies should satisfy $|\epsilon_j|\ll {\rm min}\{E_C,\Delta\}$. In addition, we employ time-dependent electromagnetic driving field which can pump single electrons between the two QDs a link. that a suitable AC voltage can be applied to a gate electrode located near this link. The Hamiltonian contribution is given by [@Platero2004] $$\label{Hdriv} H_{{\rm = w(t) d_1^{\dagger} d_2^{}+{\rm w(t) = t_{12}+2A \cos\left(\omega_0 t\right),$$ $\omega_0$ the drive and the amplitude. In what we assume that the static contribution vanishes, $t_{12}=0$, because a small coupling $t_{12}\ne 0$ will not affect dissipator in equation, see below, thus not change the a qualitative manner. In this work, Coulomb valley regime where the total charge on box is by the charging term in Eq. time scales $\delta t>1/E_C$ [@Romito2014]. The total particle on the QDs, $N_{\rm d}=\sum_j d_j^\dagger d_j^{}$, is therefore also conserved on such time scales. $N_{\rm d}\in\{0,2\}$, the inter-QD is effectively frozen We mainly on case $N_{\rm where the pair of QDs forms a spin-1/2 degree of freedom to Pauli operators $\tau_{x,y,z}$ with $\tau_\pm=(\tau_x\pm i \tau_y)/2$, $$\tau_+ = d_1
 $\epsilon_j$ for QD $j=1,2$, one arrives aT $$\label{HDotS}
 H_{\rm d}=\Sum_{J=1,2} \epSiLon_j D_j^{\daGger} d^{}_j,$$ see Ref. [@KaRZig2017] fOr details. The energies $\epSilon_J$ cAN be cONtRolleD by variATiON Of tHe GaTe vOlTAgE paraMetEr $n_g$. WitHout loss of GenErAlity, we take $\ePSiLon_2>\epsilon_1$ ThrOughout, where BotH energIeS shOUld saTisFy $|\epsIlon_j|\lL {\Rm min}\{E_c,\Delta\}$. In aDdITion, we EMploy a tIME-dEpenDent electromagnetIC dRIving field whicH can puMp SInGLE elEctRons betweeN tHe two qds via a tUNnEL LInk. tO that end, a suitAble AC voltaGE caN be appLiEd tO A gate eLectrOdE LocAted near thiS linK. The respeCtive HAMiltoniAN contriBution Is gIveN by [@PLAtErO2004] $$\laBeL{hdrIV}
H_{{\Rm dRIve}}(T) = w(t) d_1^{\daggEr} D_2^{}+{\rM h.c.},\quAd
w(t) = T_{12}+2a \COS\lefT(\omEga_0 t\Right),$$ Where $\omega_0$ denOteS the DRivE freqUency And $A$ ThE drivE ampliTude. IN wHat follows, we assUme tHat the staTic CoNtrIbUtion VAnisheS, $t_{12}=0$, bEcaUse a smaLl couplINg $t_{12}\Ne 0$ WILL nOt affect the dissipaToR IN tHe LindblAd equaTIoN, sEE Eq.  below, AnD thUs doES Not chAnge THe Physics iN a qualITaTiVe manneR.

IN this wOrK, we ConSider THe CoUlomb vAlley regIme whERe the total charGE on the box is fiXEd BY ThE CharGinG term in Eq.  on Time SCaleS $\delTA t>1/e_C$ [@ROMito2014]. THe totAl PArTIcle number on the QDs, $N_{\Rm D}=\sum_j d_J^\daggEr d_j^{}$, is therefoRe also consERVEd on such Time SCaLEs. For even $N_{\rm d}\iN\{0,2\}$, the iNter-QD dynaMIcs is effEctivEly frozeN out. We herE MAinly focUs oN thE caSe $N_{\RM D}=1$, wHere the pair of qdS forMs A spin-1/2 deGreE of freeDom CorResPonDiNg to Pauli OperatorS $\tAu_{X,y,Z}$ wIth $\Tau_\pm=(\TAu_x\pm i \taU_y)/2$, $$\Tau_+ = \TaU_-^\daGger = d_1
 $\epsilon_j$ for QD $j=1, 2$, one ar rives  at  $$ \l abel {HDo ts}
 H_{\rm d} = \sum _{j=1,2} \epsilon_j d_ j^{\d ag g er}d ^{ }_j,$ $ see R e f. [ @Ka rz ig 201 7] fo r det ail s. Theenergies $ \ep si lon_j$ can b e  c ontrolledbyvariation of  th e gate  v olt a ge pa ram eter$n_g$. Withou t loss of  g e nerali t y, we t a k e$\ep silon_2>\epsilon_ 1 $t hroughout, whe re bot he ne r g ies  sh ould satis fy  $|\e p silon_j | \l l { \rm min}\{E_C,\De lta\}$. Ina ddi tion,we  em p loy atime- de p end ent electro magn etic driv ing fi e ld whic h  can pu mp sin gle  el ectr o ns  b etw ee n  th e  t woQ Dsvia a tu nn el  link . To t h a t en d,a su itabl e AC voltagecan  bea ppl ied t o a g ateel ectro de loc atedne ar this link. T he r espective  Ha mi lto ni an co n tribut ion  is  givenby [@Pl a ter o2 0 0 4 ]$$\label{Hdriv}
H_ {{ \ r mdrive}}( t) = w ( t)  d _ 1^{\dagg er } d _2^{ } + {\rmh.c. } ,\ quad
w(t ) = t_ { 12 }+ 2A \cos \l eft(\o me ga_ 0 t \righ t ),$$  where  $\omega _0$ d e notes the driv e  frequency an d  $ A $  t h e dr ive  amplitude.  Inw hatfoll o ws , w e  assu me th at th e  static contributio nvanish es, $ t_{12}=0$, be cause a sm a l l  couplin g $t _ {1 2 }\ne 0$ will n ot af fect the d i ssipator  in t he Lindb lad equat i o n, see E q. be low , a n d  t hus does notc h ange  t he phys ics  in a q ual ita tiv e m an ner.

Inthis wor k,  w eco nsi der t h e Coulom bval le y r egime wherethe t otal  c ha r geon theb ox i s fi xe dby t hech argin g te r m i n Eq. on time s cal e s $\ de lt a t>1/E _C$ [@Romito2 01 4]. The to ta l p articl e number o n the QDs, $N_{\rm d}=\ s um_j d_ j^\ dagge r d_ j^{}$, is  th erefor e a l so con served  on s uc h t i m e sca l e s.  Fo reven $N_{\ r m  d} \in\{ 0, 2\}$ , the i nter-QD dynamics i s  ef fectively fro zen  out . We  he r em ain ly foc u s  on the case $N _{\rm d}=1 $, wh ere the pa i r o fQDs for ms a sp in-1/ 2  degree  of freed om corres po ndin g toPauli oper ators $\ tau_{x,y, z }$ wi t h$\tau _\p m=(\ta u_ x\p m i \ tau_y) / 2$,  $$\t au_+ =  \ tau_-^ \dagg er  = d_1
 $\epsilon_j$_for QD_$j=1,2$, one arrives at_$$\label{HDots}
 H_{\rm_d}=\sum_{j=1,2}_\epsilon_j d_j^{\dagger}_d^{}_j,$$_see Ref. [@Karzig2017] for_details. The energies_$\epsilon_j$ can be controlled_by variation of_the_gate voltage parameter $n_g$. Without loss of generality, we take $\epsilon_2>\epsilon_1$ throughout, where both_energies_should satisfy_$|\epsilon_j|\ll_{\rm_min}\{E_C,\Delta\}$. In addition, we employ_a time-dependent electromagnetic driving field_which can_pump single electrons between the two QDs via_a_tunnel link. To_that end, a suitable AC voltage can be applied_to a gate electrode located near_this link. The_respective_Hamiltonian_contribution is given by_[@Platero2004] $$\label{Hdriv}
H_{{\rm drive}}(t) = w(t) d_1^{\dagger}_d_2^{}+{\rm h.c.},\quad
w(t) = t_{12}+2A \cos\left(\omega_0 t\right),$$_where $\omega_0$ denotes the drive frequency and_$A$ the drive amplitude. In what_follows, we assume that the_static contribution_vanishes, $t_{12}=0$, because a small_coupling $t_{12}\ne 0$_will not_affect the dissipator_in the Lindblad equation, see Eq. _below, and thus_does not change the physics in_a_qualitative manner.

In this_work,_we_consider the_Coulomb valley regime_where_the total_charge_on the box is fixed by_the_charging term in Eq.  on time scales_$\delta t>1/E_C$ [@Romito2014]. The_total_particle number on the_QDs, $N_{\rm d}=\sum_j d_j^\dagger d_j^{}$,_is therefore also conserved on such_time scales._For even_$N_{\rm d}\in\{0,2\}$, the inter-QD dynamics is effectively frozen out. We here_mainly focus on the case $N_{\rm_d}=1$, where the pair_of QDs_forms_a spin-1/2 degree_of_freedom corresponding_to Pauli operators $\tau_{x,y,z}$ with $\tau_\pm=(\tau_x\pm i_\tau_y)/2$, $$\tau_+_= \tau_-^\dagger = d_1