 over the QD subsystem, we arrive at a Lindblad equation for the $2\times 2$ density matrix $\rho_{\rm M}(t)$ only, $$\label{LindbladMBQ}
\partial_t\rho_{\rm M}(t) = -i[\tilde H_{\rm L},\rho_{\rm M} ]+  \sum_{s=\pm}  \tilde\Gamma_s \mathcal{L}[\tilde J_s]\rho_{\rm M}(t),$$ where the jump operators $\tilde J_\pm$ have been defined in Eq. . The dissipative transition rates $\tilde \Gamma_\pm$ in Eq.  are given by $$\label{mod1}
  \tilde \Gamma_+ = p \Gamma_+,\quad \tilde \Gamma_- = (1-p)\Gamma_-,$$ cf. Eqs.  and. The coherent time evolution in Eq.  is governed by the Hamiltonian $$\label{mod22}
   \tilde H_{\rm L} =  2(1-2p)g_0 \tilde J_z + \sum_\pm \tilde h_\pm \tilde J_\pm^\dagger \tilde J_\pm^{},$$ where $\tilde J_z$ has been specified in Eq.  and the Lamb shifts $\tilde h_\pm$ are given by $$\label{mod2}
    \tilde h_+ = p h_+,\quad \tilde h_- = (1-p) h_-.$$ The drive amplitude $A$ then appears only implicitly through the dependence $p=p(A)$. We note that within the RWA, no contributions $\propto p_\perp$ appear in Eq. . Indeed, the RWA allows one to neglect terms $\propto \tau_+\rho\tau_+$ which stem from $p_\perp\ne 0$.

Importantly, apart from the initial transient behavior, all of our numerical results for the Majorana dynamics obtained from the full Lindblad equation for the combined QD-MBS system, Eq. , are quantitatively reproduced by using the simpler Lindblad equation. This statement is valid for arbitrary model parameters subject to Eqs.  and. We emphasize that the integration over the QD degrees of freedom as carried out above relies on the facts that (i) the convergence towards the target state is dictated by the Majorana sector, and that (ii) the QD and MBS degrees of freedom always decouple
over the QD subsystem, we arrive at a Lindblad equation for the $ 2\times 2 $ density matrix $ \rho_{\rm M}(t)$ only, $ $ \label{LindbladMBQ } 
 \partial_t\rho_{\rm M}(t) = -i[\tilde H_{\rm L},\rho_{\rm M } ] +   \sum_{s=\pm }   \tilde\Gamma_s \mathcal{L}[\tilde J_s]\rho_{\rm M}(t),$$ where the jump operators $ \tilde J_\pm$ have been define in Eq.  . The dissipative passage rates $ \tilde \Gamma_\pm$ in Eq.    are given by $ $ \label{mod1 } 
   \tilde \Gamma_+ = p \Gamma_+,\quad \tilde \Gamma_- = (1 - p)\Gamma_-,$$ cystic fibrosis.   Eqs.    and. The coherent time evolution in Eq.    is govern by the Hamiltonian $ $ \label{mod22 } 
    \tilde H_{\rm fifty } =   2(1 - 2p)g_0 \tilde J_z + \sum_\pm \tilde h_\pm \tilde J_\pm^\dagger \tilde J_\pm^{},$$ where $ \tilde J_z$ has been specified in Eq.    and the Lamb shifts $ \tilde h_\pm$ are given by $ $ \label{mod2 } 
     \tilde h_+ = p h_+,\quad \tilde h_- = (1 - phosphorus) h_-.$$ The drive amplitude $ A$ then appears only implicitly through the dependence $ phosphorus = p(A)$. We note that within the RWA, no contributions $ \propto p_\perp$ look in Eq.  . Indeed, the RWA allows one to fail term $ \propto \tau_+\rho\tau_+$ which stem from $ p_\perp\ne 0$. 

 Importantly, apart from the initial transient behavior, all of our numeral results for the Majorana dynamics obtain from the full Lindblad equality for the combined QD - bachelor of medicine system, Eq.  , are quantitatively reproduced by using the simple Lindblad equation. This statement is valid for arbitrary model parameters subject to Eqs.    and. We emphasize that the integration over the QD degrees of freedom as carried out above relies on the facts that (i) the convergence towards the target state is dictated by the Majorana sector, and that (ii) the QD and MBS degrees of exemption constantly decouple
 ovfr the QD subsystem, we avrive at a Lindblad equetion fkr the $2\tkmes 2$ density matrix $\rho_{\rm M}(v)$ onoy, $$\lavel{LindbladMBQ}
\partial_t\fho_{\rm M}(t) = -i[\tilde Y_{\rm O},\eho_{\rm M} ]+  \sum_{s=\pm}  \bnlde\Gzlma_s \nathcal{L}[\tilde M_s]\rho_{\rm M}(t),$$ fhere the jump ooexators $\tilde J_\pm$ have been defined ig Eq. . Thr fissipative trwnsiuiog rafvs $\tilde \Gamma_\pm$ in Eq.  are giveh by $$\lauel{mod1}
  \tilde \Gsmma_+ = p \Gamma_+,\quad \tilde \Galma_- = (1-p)\Gamma_-,$$ cf. Eqs.  and. Hhe coherenj tiiw evolution kn Eq.  is governed by tge Hamiltonian $$\label{mod22}
   \tilde H_{\rm K} =  2(1-2p)g_0 \tilee J_z + \sum_\pm \tildx h_\pm \nilde J_\pm^\daggcg \tilde J_\pm^{},$$ whrre $\tilde J_z$ hss ueen specified in Eq.  and vhe Lamb shifts $\tildg h_\pm$ are ciren by $$\label{mod2}
    \tiodw h_+ = k h_+,\qudd \tkode h_- = (1-p) h_-.$$ The dgivx amplitude $A$ then appwars only implicitlu eyrough the delendense $p=p(A)$. We note that within the RWA, no conuribufions $\propto p_\perp$ appeqr in Eq. . Indeed, the RAA allows one to neglect terms $\propto \tau_+\rho\tau_+$ which stem from $o_\pexi\ke 0$.

Iopogtantly, apart from the initial transient behadjot, sll of our numcrical results for tje Iajorana dynaoics oyfajned from the full Lindblwd eqyation fow thr combined QD-MBS system, Eq. , qre quantitaniveoy reproduced by uding the sioplet Lindnlad equation. This statzment js valid fog arbitradh model parametefs xutject to Eqs.  and. We emphasyze that vhe iutegratiun ofer thq QD degrefs of freedom as carried out cbove relies on the facts that (i) the convergenrx towards the tdrgvt state ns dicbated by the Matorana sector, cnd that (ii) tfe QD and JBS degcees of freeqom always dewluple
over the QD subsystem, we arrive at equation the $2\times density matrix $\rho_{\rm = H_{\rm L},\rho_{\rm M} \sum_{s=\pm} \tilde\Gamma_s \mathcal{L}[\tilde M}(t),$$ where the jump operators $\tilde have been defined in Eq. . The dissipative transition rates $\tilde \Gamma_\pm$ in are given by $$\label{mod1} \tilde \Gamma_+ = p \Gamma_+,\quad \tilde \Gamma_- = (1-p)\Gamma_-,$$ Eqs. The time in Eq. is governed by the Hamiltonian $$\label{mod22} \tilde H_{\rm L} = 2(1-2p)g_0 \tilde J_z + \tilde h_\pm \tilde J_\pm^\dagger \tilde J_\pm^{},$$ where $\tilde has been specified in and the Lamb shifts $\tilde are by $$\label{mod2} h_+ p \tilde h_- = h_-.$$ The drive amplitude $A$ then appears only implicitly through the dependence $p=p(A)$. We note that within RWA, no p_\perp$ appear Eq. Indeed, RWA allows one terms $\propto \tau_+\rho\tau_+$ which stem from apart from the initial transient behavior, all of numerical results the Majorana dynamics obtained from the Lindblad equation for the combined QD-MBS system, Eq. are quantitatively reproduced by using the simpler Lindblad equation. This statement is valid for arbitrary subject to Eqs. and. emphasize that the over QD of as carried above relies on the facts that (i) the convergence towards the state is dictated by the Majorana sector, and that (ii) and degrees of freedom decouple
 over the QD subsystem, we arrivE at a LindblAd equAtiOn fOr The $2\tImes 2$ Density matrix $\rHO_{\rm M}(T)$ only, $$\label{LindbladMBQ}
\pArtiaL_t\RHo_{\rm m}(T) = -i[\Tilde h_{\rm L},\rho_{\RM M} ]+  \SUM_{s=\pM}  \tIlDe\GAmMA_s \MathcAl{L}[\Tilde J_s]\Rho_{\rm M}(t),$$ wheRe tHe Jump operatorS $\TiLde J_\pm$ have BeeN defined in Eq. . the DissipAtIve TRansiTioN rateS $\tilde \gAmma_\pm$ In Eq.  are giVeN By $$\labeL{Mod1}
  \tildE \gAmMa_+ = p \GAmma_+,\quad \tilde \GammA_- = (1-P)\GAMma_-,$$ cf. Eqs.  and. The CohereNt TImE EVolUtiOn in Eq.  is goVeRned bY The HamiLToNIAN $$\laBEl{mod22}
   \tilde H_{\rm l} =  2(1-2p)g_0 \tilde J_z + \sUM_\pm \Tilde h_\Pm \TilDE J_\pm^\daGger \tIlDE J_\pM^{},$$ where $\tilde j_z$ haS been specIfied iN eq.  and thE lamb shiFts $\tilDe h_\Pm$ aRe giVEn By $$\LabEl{MOd2}
    \tILdE h_+ = p H_+,\QuaD \tilde h_- = (1-p) H_-.$$ THe Drive AmplITUDE $A$ thEn aPpeaRs onlY implicitly thRouGh thE DepEndenCe $p=p(A)$. we noTe That wIthin tHe RWA, No Contributions $\prOpto P_\perp$ appeAr iN EQ. . InDeEd, the rwA alloWs oNe tO neglecT terms $\pROptO \tAU_+\RHo\Tau_+$ which stem from $p_\pErP\NE 0$.

IMportantLy, aparT FrOm THe initiaL tRanSienT BEhaviOr, alL Of Our numerIcal reSUlTs For the MAjOrana dYnAmiCs oBtainED froM the fuLl LindblAd equATion for the combINed QD-MBS systeM, eq. , ARE qUAntiTatIvely reprodUced BY usiNg thE SiMplER LindBlad eQuATiON. This statement is valId For arbItrarY model parametErs subject TO eQs.  and. We eMphaSIzE That the integraTion oVer the QD deGRees of frEedom As carrieD out above RELies on thE faCts ThaT (i) tHE CoNvergence towaRDS the TaRget staTe iS dictatEd bY thE MaJorAnA sector, anD that (ii) tHe qD AnD MbS dEgreeS Of freedoM aLwaYs DecOuple
 over the QD subsystem, we  arrive at  a Li ndb lad  e quat ionfor the $2\tim e s 2$  density matrix $\rho_ {\rmM} ( t)$o nl y, $$ \label{ L in d b lad MB Q} 
\p ar t ia l_t\r ho_ {\rm M} (t) = -i[\ til de  H_{\rm L},\ r ho _{\rm M} ] + \sum_{s=\pm}   \ tilde\ Ga mma _ s \ma thc al{L} [\tild e  J_s]\ rho_{\rmM} ( t),$$w here th e ju mp o perators $\tildeJ _\ p m$ have been d efined  i n  E q .  .The  dissipati ve  tran s ition r a te s $ \ti l de \Gamma_\pm $ in Eq.  a r e g iven b y$$\ l abel{m od1} \ til de \Gamma_+  = p  \Gamma_+ ,\quad \tilde\ Gamma_-  = (1- p)\ Gam ma_- , $$  c f.Eq s . a nd . T h e c oherentti me  evol utio n i n  Eq.   i s go verne d by the Hami lto nian $$\ label {mod2 2}
 \tild e H_{\ rm L}  =   2(1-2p)g_0 \t ilde  J_z + \s um_ \p m \ ti lde h _ \pm \t ild e J _\pm^\d agger \ t ild eJ _ \ pm ^{},$$ where $\til de J _z $ has be en spe c if ie d  in Eq. a ndtheL a mb sh ifts $\ tilde h_ \pm$ a r egi ven by$$ \label {m od2 }
   \t i ldeh_+ =p h_+,\q uad \ t ilde h_- = (1- p ) h_-.$$ Thed ri v e  a m plit ude  $A$ then a ppea r s on ly i m pl ici t ly th rough  t h ed ependence $p=p(A)$.  W e note  that  within the R WA, no con t r i butions$\pr o pt o  p_\perp$ appe ar in  Eq. . Ind e ed, theRWA a llows on e to negl e c t terms$\p rop to\ta u _ +\ rho\tau_+$ wh i c h st em  from $ p_\ perp\ne  0$ .
Imp ort an tly, apar t from t he  i ni ti altrans i ent beha vi or,  a llof ou r  numer icalresu lt sf orthe Maj o ra n a  dyn am ic s ob tai ne d fro m th e  fu ll Lind blad equa tio n  for  t he  combin ed QD-MBS sys te m, Eq. , a re  qu antita t i vely rep roduced by using the si m pler Li ndb lad e quat ion. This  st atemen t i s  valid  for a rbitr ar y m o d el pa r a me ter ssubject to E qs.   and .We e mphasiz e that the integra t ion  over the QDdeg rees o ffre e do m  as  c a rri e d  out above reli es on thefa c ts  that (i)t heco nvergen ce towa rds t h e targe t state i s dictate dby t h e  Ma jorana sec tor, and  that (ii )  theQ Dand M BSdegree soffreed om alw a ysdecou ple
 over_the QD_subsystem, we arrive at_a Lindblad_equation_for the_$2\times_2$ density matrix_$\rho_{\rm M}(t)$ only,_$$\label{LindbladMBQ}
\partial_t\rho_{\rm M}(t) = -i[\tilde_H_{\rm L},\rho_{\rm M}_]+_ \sum_{s=\pm}  \tilde\Gamma_s \mathcal{L}[\tilde J_s]\rho_{\rm M}(t),$$ where the jump operators $\tilde J_\pm$ have_been_defined in_Eq. ._The_dissipative transition rates $\tilde \Gamma_\pm$_in Eq.  are given by_$$\label{mod1}
 _\tilde \Gamma_+ = p \Gamma_+,\quad \tilde \Gamma_- =_(1-p)\Gamma_-,$$_cf. Eqs.  and. The_coherent time evolution in Eq.  is governed by the_Hamiltonian $$\label{mod22}
   \tilde H_{\rm_L} = _2(1-2p)g_0_\tilde_J_z + \sum_\pm \tilde_h_\pm \tilde J_\pm^\dagger \tilde J_\pm^{},$$ where_$\tilde J_z$ has been specified in_Eq.  and the Lamb shifts $\tilde h_\pm$_are given by $$\label{mod2}
  _ \tilde h_+ = p_h_+,\quad \tilde_h_- = (1-p) h_-.$$ The_drive amplitude $A$_then appears_only implicitly through_the dependence $p=p(A)$. We note that_within the RWA,_no contributions $\propto p_\perp$ appear in_Eq. ._Indeed, the RWA_allows_one_to neglect_terms $\propto \tau_+\rho\tau_+$_which_stem from_$p_\perp\ne_0$.

Importantly, apart from the initial transient_behavior,_all of our numerical results for the_Majorana dynamics obtained from_the_full Lindblad equation for_the combined QD-MBS system, Eq. ,_are quantitatively reproduced by using the_simpler Lindblad_equation. This_statement is valid for arbitrary model parameters subject to Eqs.  and._We emphasize that the integration over_the QD degrees of_freedom as_carried_out above relies_on_the facts_that (i) the convergence towards the target_state is_dictated by the Majorana sector, and_that (ii) the QD_and_MBS degrees of freedom always decouple