ansatz2}
\rho_{\rm d}=  {\rm diag}\left(\frac{1-p}{2}, \frac{1-p}{2},p\right),$$ expressed in the basis $\{ |100\rangle,|010\rangle,|001\rangle \}$ with QD occupation states $|n_1,n_2,n_3\rangle$ for $N_{\rm d}=1$. Note that since we assumed $\epsilon_1=\epsilon_2$, the occupation probabilities of QDs 1 and 2 are equal. The occupation probability $0<p\ll 1$ refers to the energetically highest QD 3. Equation is consistent with our numerical analysis of the Lindblad equation for $\rho(t)$, where we again find a factorized density matrix at long times, $\rho(t)\simeq \rho_{\rm M}(t)\otimes \rho_d$. We note that the dark space turns out to be independent of the concrete value of $p$.

Going through the corresponding steps in Sec. \[sec2c\], we arrive at a Lindblad equation for $\rho_{\rm M}(t)$, $$\label{LindbladTwoMBQ}
    \partial_t\rho_{\rm M}(t)=-i[\tilde H_{\rm L},\rho_{\rm M}(t)]+
    \sum_{a=1}^6 \tilde \Gamma_a \mathcal{L}[K_a] \rho_{\rm M}(t).$$ The six jump operators are denoted by $K_a$, with the respective dissipative transition rates $\tilde\Gamma_a$. With $\lambda_{LR}\equiv t_{LR}/E_C \ll 1$, we obtain $$\begin{aligned}
\nonumber
  K_1^{} &=& K_4^\dagger = ie^{i(\beta_3-\beta_1)}\frac{|\lambda_{1,1R}\lambda_{3,3R}|}{\lambda_{LR}} X_R \\
  &&\qquad - \,e^{i\beta_3} |\lambda_{1,3L}\lambda_{3,3R}|  Z_L Y_R,\nonumber\\ \label{jumpK}
  K_2^{} &=& K_5^\dagger =-ie^{i(\beta_3-\beta_2)}\frac{|\lambda_{2,4R}\lambda_{3,3R}|}{\lambda_{LR}}Z_
ansatz2 } 
 \rho_{\rm d}=   { \rm diag}\left(\frac{1 - p}{2 }, \frac{1 - p}{2},p\right),$$ expressed in the basis $ \ { |100\rangle,|010\rangle,|001\rangle \}$ with QD occupation states $ |n_1,n_2,n_3\rangle$ for $ N_{\rm d}=1$. Note that since we assumed $ \epsilon_1=\epsilon_2 $, the occupation probability of QDs 1 and 2 are adequate. The occupation probability $ 0 < p\ll 1 $ refers to the energetically highest QD 3. equality is consistent with our numeral analysis of the Lindblad equation for $ \rho(t)$, where we again discover a factorized density matrix at retentive times, $ \rho(t)\simeq \rho_{\rm M}(t)\otimes \rho_d$. We note that the colored space turns out to be autonomous of the concrete value of $ p$. 

 Going through the corresponding steps in Sec.   \[sec2c\ ], we arrive at a Lindblad equation for $ \rho_{\rm M}(t)$, $ $ \label{LindbladTwoMBQ } 
     \partial_t\rho_{\rm M}(t)=-i[\tilde H_{\rm L},\rho_{\rm M}(t)]+ 
     \sum_{a=1}^6 \tilde \Gamma_a \mathcal{L}[K_a ] \rho_{\rm M}(t).$$ The six jump operators are denoted by $ K_a$, with the respective dissipative passage rates $ \tilde\Gamma_a$. With $ \lambda_{LR}\equiv t_{LR}/E_C \ll 1 $, we obtain $ $ \begin{aligned } 
 \nonumber 
   K_1^ { } & = & K_4^\dagger = ie^{i(\beta_3-\beta_1)}\frac{|\lambda_{1,1R}\lambda_{3,3R}|}{\lambda_{LR } } X_R \\ 
   & & \qquad - \,e^{i\beta_3 } |\lambda_{1,3L}\lambda_{3,3R}|   Z_L Y_R,\nonumber\\ \label{jumpK } 
   K_2^ { } & = & K_5^\dagger = -ie^{i(\beta_3-\beta_2)}\frac{|\lambda_{2,4R}\lambda_{3,3R}|}{\lambda_{LR}}Z _
answtz2}
\rho_{\rm d}=  {\rm diag}\left(\fvac{1-p}{2}, \frac{1-p}{2},p\right),$$ expresved in the basks $\{ |100\rangle,|010\rangle,|001\rangle \}$ with QE occypation states $|n_1,n_2,n_3\rangue$ for $N_{\rl d}=1$. Note thau since we assumev $\epsilok_1=\zpsilkk_2$, the iccupation pronabilities mf QDs 1 and 2 ase eeual. The occupation probability $0<p\ll 1$ referx ho the energetycalkr hifhest QD 3. Equation is consistent wjth our numerical anslysis of the Lindblad equwtioj for $\rho(t)$, where wf again fine a sqctorized devsity matrpr at long tjmes, $\rho(t)\simeq \rho_{\rm M}(t)\otimes \rfo_d$. Wz note that tye fdrk space tnrns olt to be indeivndent mf the voncrete value of $p$.

Giing through the corrxsponding steps in Sgc. \[sec2c\], we drxive at a Lindblad eqyarion xor $\sho_{\ro M}(t)$, $$\lagek{LjndblafTwkMBQ}
    \padtial_t\rho_{\rm M}(t)=-i[\tilde H_{\rm L},\rho_{\rk I}(n)]+
    \sum_{a=1}^6 \tilse \Gamia_w \mathcal{L}[K_a] \rho_{\rm M}(t).$$ The six jump operdtods are denoted by $K_a$, wirh the respective disdipative eransition rates $\tilde\Gamma_a$. With $\lambda_{LR}\equiv t_{NR}/E_C \mu 1$, cc uvtwin $$\begin{aligned}
\nonumber
  K_1^{} &=& K_4^\dagger = ie^{i(\beta_3-\fsts_1)}\fgac{|\lambda_{1,1R}\lambda_{3,3R}|}{\jambda_{LR}} X_R \\
  &&\qaoad - \,e^{i\beta_3} |\laobda_{1,3L}\lcjbsa_{3,3R}|  Z_L Y_R,\nonumber\\ \label{jompK}
  J_2^{} &=& K_5^\daggew =-ie^{o(\beta_3-\beta_2)}\frac{|\lambda_{2,4R}\lambda_{3,3R}|}{\oambda_{LR}}Z_
ansatz2} \rho_{\rm d}= {\rm diag}\left(\frac{1-p}{2}, \frac{1-p}{2},p\right),$$ expressed basis |100\rangle,|010\rangle,|001\rangle \}$ QD occupation states that we assumed $\epsilon_1=\epsilon_2$, occupation probabilities of 1 and 2 are equal. The probability $0<p\ll 1$ refers to the energetically highest QD 3. Equation is consistent our numerical analysis of the Lindblad equation for $\rho(t)$, where we again find factorized matrix long $\rho(t)\simeq \rho_{\rm M}(t)\otimes \rho_d$. We note that the dark space turns out to be independent of concrete value of $p$. Going through the corresponding in Sec. \[sec2c\], we at a Lindblad equation for M}(t)$, \partial_t\rho_{\rm M}(t)=-i[\tilde L},\rho_{\rm \sum_{a=1}^6 \Gamma_a \mathcal{L}[K_a] \rho_{\rm The six jump operators are denoted by $K_a$, with the respective dissipative transition rates $\tilde\Gamma_a$. With $\lambda_{LR}\equiv \ll 1$, $$\begin{aligned} \nonumber &=& = X_R \\ &&\qquad |\lambda_{1,3L}\lambda_{3,3R}| Z_L Y_R,\nonumber\\ \label{jumpK} K_2^{} &=&
ansatz2}
\rho_{\rm d}=  {\rm diag}\left(\fraC{1-p}{2}, \frac{1-p}{2},p\riGht),$$ exPreSseD iN the BasiS $\{ |100\rangle,|010\rangle,|001\rANgle \}$ With QD occupation states $|N_1,n_2,n_3\raNgLE$ for $n_{\Rm D}=1$. Note That sinCE wE ASsuMeD $\ePsiLoN_1=\EpSilon_2$, The OccupatIon probabiLitIeS of QDs 1 and 2 are EQuAl. The occupAtiOn probabilitY $0<p\lL 1$ referS tO thE EnergEtiCally HighesT qD 3. EquaTion is conSiSTent wiTH our numERIcAl anAlysis of the LindblAD eQUation for $\rho(t)$, wHere we AgAIn FINd a FacTorized denSiTy matRIx at lonG TiMES, $\Rho(T)\Simeq \rho_{\rm M}(t)\oTimes \rho_d$. We NOte That thE dArk SPace tuRns ouT tO Be iNdependent oF the Concrete vAlue of $P$.

going thROugh the CorresPonDinG stePS iN SEc. \[sEc2C\], We aRRiVe aT A LiNdblad eqUaTiOn for $\Rho_{\rM m}(T)$, $$\LAbel{linDblaDTwoMbQ}
    \partial_t\rho_{\Rm M}(T)=-i[\tiLDe H_{\Rm L},\rhO_{\rm M}(t)]+
    \Sum_{a=1}^6 \TiLde \GaMma_a \maThcal{l}[K_A] \rho_{\rm M}(t).$$ The six jUmp oPerators aRe dEnOteD bY $K_a$, wiTH the reSpeCtiVe dissiPative tRAnsItION RaTes $\tilde\Gamma_a$. With $\LaMBDa_{lR}\equiv t_{lR}/E_C \ll 1$, WE oBtAIn $$\begin{aLiGneD}
\nonUMBer
  K_1^{} &=& K_4^\DaggER = iE^{i(\beta_3-\beTa_1)}\frac{|\LAmBdA_{1,1R}\lambdA_{3,3R}|}{\Lambda_{lR}} x_R \\
  &&\qQuaD - \,e^{i\beTA_3} |\lamBda_{1,3L}\laMbda_{3,3R}|  Z_L Y_r,\nonuMBer\\ \label{jumpK}
  K_2^{} &=& k_5^\Dagger =-ie^{i(\beta_3-\BEtA_2)}\FRaC{|\LambDa_{2,4R}\Lambda_{3,3R}|}{\lambDa_{LR}}z_
ansatz2}
\rho_{\rm d}=  {\ rm diag}\l eft(\ fra c{1 -p }{2} , \f rac{1-p}{2},p\ r ight ),$$ expressed in thebasis  $ \ { |1 0 0\ rangl e,|010\ r an g l e,| 00 1\ ran gl e  \ }$ wi thQD occu pation sta tes  $ |n_1,n_2,n_3 \ ra ngle$ for$N_ {\rm d}=1$.Not e that  s inc e  we a ssu med $ \epsil o n_1=\e psilon_2$ ,t he occ u pationp r ob abil ities of QDs 1 an d  2 are equal. The  occup at i on p rob abi lity $0<p\ ll  1$ r e fers to th e e ner g etically high est QD 3. E q uat ion is  c ons i stentwithou r  nu merical ana lysi s of theLindbl a d equat i on for$\rho( t)$ , w here we  a gai nf ind afac t ori zed dens it ymatri x at l o n g ti mes , $\ rho(t )\simeq \rho_ {\r m M} ( t)\ otime s \rh o_d$ .We no te tha t the  d ark space turns  out  to be in dep en den tof th e  concr ete  va lue of$p$.

G o ing  t h r o ug h the correspondin gs t ep s in Sec . \[se c 2c \] ,  we arri ve  at  a L i n dblad  equ a ti on for $ \rho_{ \ rm  M }(t)$,$$ \label {L ind bla dTwoM B Q}
   \pa rtial_t\ rho_{ \ rm M}(t)=-i[\t i lde H_{\rm L} , \r h o _{ \ rm M }(t )]+
    \su m_{a = 1}^6  \ti l de  \G a mma_a  \mat hc a l{ L }[K_a] \rho_{\rm M} (t ).$$ T he si x jump operat ors are de n o t ed by $K _a$, wi t h the respecti ve di ssipativet ransitio n rat es $\til de\Gamma_ a $ . With $ \la mbd a_{ LR} \ e qu iv t_{LR}/E_C \ ll 1 $,  we obt ain  $$\beg in{ ali gne d}\n onumber
 K_1^{}&= &K_ 4^ \da gger=  ie^{i(\ be ta_ 3- \be ta_1) } \frac{ |\lam bda_ {1 ,1 R }\l ambda_{ 3 ,3 R } |}{\ la mb da_{ LR} }X_R \ \
 & &\q quad -\,e^{i\be ta_ 3 } |\ la mb da_{1,3 L}\lambda_{3, 3R }|  Z_L Y_ R, \no number \ \  \label{ jumpK}
  K_2^{} &=& K_5 ^ \dagger  =- ie^{i (\be ta_3-\bet a_2 )}\fra c{| \ lambda _{2,4R }\lam bd a_{ 3 , 3R}|} { \ la mbd a_ {LR}}Z_
ansatz2}
\rho_{\rm d}=_ {\rm_diag}\left(\frac{1-p}{2}, \frac{1-p}{2},p\right),$$ expressed in_the basis_$\{_|100\rangle,|010\rangle,|001\rangle \}$_with_QD occupation states_$|n_1,n_2,n_3\rangle$ for $N_{\rm_d}=1$. Note that since_we assumed $\epsilon_1=\epsilon_2$,_the_occupation probabilities of QDs 1 and 2 are equal. The occupation probability $0<p\ll 1$_refers_to the_energetically_highest_QD 3. Equation is consistent_with our numerical analysis of_the Lindblad_equation for $\rho(t)$, where we again find a_factorized_density matrix at_long times, $\rho(t)\simeq \rho_{\rm M}(t)\otimes \rho_d$. We note that_the dark space turns out to_be independent of_the_concrete_value of $p$.

Going through_the corresponding steps in Sec. \[sec2c\], we_arrive at a Lindblad equation for_$\rho_{\rm M}(t)$, $$\label{LindbladTwoMBQ}
    \partial_t\rho_{\rm_M}(t)=-i[\tilde H_{\rm L},\rho_{\rm M}(t)]+
  _ \sum_{a=1}^6 \tilde \Gamma_a \mathcal{L}[K_a]_\rho_{\rm M}(t).$$_The six jump operators are_denoted by $K_a$,_with the_respective dissipative transition_rates $\tilde\Gamma_a$. With $\lambda_{LR}\equiv t_{LR}/E_C \ll_1$, we obtain_$$\begin{aligned}
\nonumber
  K_1^{} &=& K_4^\dagger =_ie^{i(\beta_3-\beta_1)}\frac{|\lambda_{1,1R}\lambda_{3,3R}|}{\lambda_{LR}}_X_R \\
 _&&\qquad_-_\,e^{i\beta_3} |\lambda_{1,3L}\lambda_{3,3R}|_ Z_L Y_R,\nonumber\\_\label{jumpK}
_ K_2^{}_&=&_K_5^\dagger =-ie^{i(\beta_3-\beta_2)}\frac{|\lambda_{2,4R}\lambda_{3,3R}|}{\lambda_{LR}}Z_