mathbf{C}}_{11}^{-1} = \left( \breve{\mathbf{C}}_{12} \breve{\mathbf{C}}_{22}^{-1} \right)^\mathsf{T},\end{aligned}$$ holds by considering the orthogonality of the eigenvector matrix $\breve{\mathbf{C}}$, $$\begin{aligned}
\mathbf{I} = \breve{\mathbf{C}}^\mathsf{T} \breve{\mathbf{C}}.\end{aligned}$$ This implies that subsystem and environment eigenvectors are orthogonal, leading to the expression, $$\begin{aligned}
\mathbf{0} 
&= \breve{\mathbf{C}}_\mathcal{E}^\mathsf{T} \breve{\mathbf{C}}_\mathcal{S}\\
&= \breve{\mathbf{C}}_{12}^\mathsf{T} \breve{\mathbf{C}}_{11} + \breve{\mathbf{C}}_{22}^\mathsf{T} \breve{\mathbf{C}}_{21} \\
&= \breve{\mathbf{C}}_{22}^{-\mathsf{T}} \breve{\mathbf{C}}_{12}^\mathsf{T} + \breve{\mathbf{C}}_{21} \breve{\mathbf{C}}_{11}^{-1} \\
&= \left( \breve{\mathbf{C}}_{12} \breve{\mathbf{C}}_{22}^{-1} \right)^\mathsf{T} + \breve{\mathbf{C}}_{21} \breve{\mathbf{C}}_{11}^{-1},\end{aligned}$$ from which Eq.  immediately follows.

Representation of the block-diagonalization matrix $\mathbf{Q}$ {#app:q}
===============================================================

In the literature,[@Cederbaum89; @Sikkema09; @Peng12; @Seino12] there exist multiple representations of the block-diagonalization matrix $\mathbf{Q}$. These representations are very similar to the expression in Eq.  and it is trivial to show that they achieve the exact same block-diagonalization. The different representations arise from the following different definitions:

1.  The definition of the block-diagonalization as either $$\begin{aligned}
    \tilde{\mathbf{F}} = \mathbf{Q} \mathbf{F} \mathbf{Q}^\mathsf{T},\end{aligned}$$ or $$\begin{aligned}
    \tilde{\mathbf{F}} = \
mathbf{C}}_{11}^{-1 } = \left (\breve{\mathbf{C}}_{12 } \breve{\mathbf{C}}_{22}^{-1 } \right)^\mathsf{T},\end{aligned}$$ holds by considering the orthogonality of the eigenvector matrix $ \breve{\mathbf{C}}$, $ $ \begin{aligned } 
 \mathbf{I } = \breve{\mathbf{C}}^\mathsf{T } \breve{\mathbf{C}}.\end{aligned}$$ This implies that subsystem and environment eigenvectors are orthogonal, lead to the formula, $ $ \begin{aligned } 
 \mathbf{0 } 
 & = \breve{\mathbf{C}}_\mathcal{E}^\mathsf{T } \breve{\mathbf{C}}_\mathcal{S}\\ 
 & = \breve{\mathbf{C}}_{12}^\mathsf{T } \breve{\mathbf{C}}_{11 } + \breve{\mathbf{C}}_{22}^\mathsf{T } \breve{\mathbf{C}}_{21 } \\ 
 & = \breve{\mathbf{C}}_{22}^{-\mathsf{T } } \breve{\mathbf{C}}_{12}^\mathsf{T } + \breve{\mathbf{C}}_{21 } \breve{\mathbf{C}}_{11}^{-1 } \\ 
 & = \left (\breve{\mathbf{C}}_{12 } \breve{\mathbf{C}}_{22}^{-1 } \right)^\mathsf{T } + \breve{\mathbf{C}}_{21 } \breve{\mathbf{C}}_{11}^{-1},\end{aligned}$$ from which Eq.    immediately follows. 

 Representation of the block - diagonalization matrix $ \mathbf{Q}$ { # app: q } 
 = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = 

 In the literature,[@Cederbaum89; @Sikkema09; @Peng12; @Seino12 ] there exist multiple representations of the engine block - diagonalization matrix $ \mathbf{Q}$. These representations are very similar to the construction in Eq.    and it is superficial to show that they achieve the accurate like block - diagonalization. The different representation arise from the following different definitions: 

 1.   The definition of the block - diagonalization as either $ $ \begin{aligned } 
     \tilde{\mathbf{F } } = \mathbf{Q } \mathbf{F } \mathbf{Q}^\mathsf{T},\end{aligned}$$ or $ $ \begin{aligned } 
     \tilde{\mathbf{F } } = \
matjbf{C}}_{11}^{-1} = \left( \breve{\mathbf{C}}_{12} \nreve{\mathbf{C}}_{22}^{-1} \rigkr)^\mathsh{T},\end{aljgned}$$ houds by considering the orthojonaoity if the eigenvector matfix $\breve{\lathbf{C}}$, $$\vegii{aligned}
\mathbf{I} = \breve{\mabkbf{C}}^\mzbhsf{T} \ureve{\mathbf{C}}.\end{sligned}$$ Thhs implies thad ruysystem and environment eigenvectors are oryhlgonal, leading to uhe expdvswion, $$\begin{aligned}
\mathbf{0} 
&= \brevs{\mathbf{B}}_\mathcal{E}^\mathsf{T} \nreve{\mathbf{C}}_\mathcal{S}\\
&= \breve{\lathhf{C}}_{12}^\mathsf{T} \breve{\mahhbf{C}}_{11} + \brevg{\jatrvf{C}}_{22}^\mathsf{T} \bfeve{\mathbf{B}}_{21} \\
&= \breve{\matgbf{C}}_{22}^{-\mathsf{T}} \breve{\mathbf{C}}_{12}^\mathsf{T} + \brere{\mathbf{C}}_{21} \btzce{\mwjhbf{C}}_{11}^{-1} \\
&= \left( \ureve{\mwthbf{C}}_{12} \breve{\mathbf{C}}_{22}^{-1} \sight)^\mayhsf{T} + \breve{\mabhbf{C}}_{21} \brwve{\mathbf{C}}_{11}^{-1},\end{aligned}$$ hrom which Eq.  immediwtely folnocs.

Representation of tye blocn-diaconauuzagioh kafrix $\mwthuf{Q}$ {#app:q}
===============================================================

In fhe literatyre,[@Cederbaum89; @Sikkems09; @[vmg12; @Seino12] thede exift multiple representations of the block-dpagohalization matrix $\mathbd{Q}$. These representatilns are vqry similar to the expression in Eq.  and it is trieial vo shiw gyah they achieve the exact same block-diagonalizwfipn. The different representayilnx arise from tfe followjng different defijitions:

1.  The definitijn og the block-diagonalization qs either $$\beyin{qligned}
    \tilde{\machbf{F}} = \mathby{Q} \matnbf{F} \kathbf{Q}^\mathsf{T},\end{aligneb}$$ or $$\bsgin{aligned}
    \tilde{\jxthbf{F}} = \
mathbf{C}}_{11}^{-1} = \left( \breve{\mathbf{C}}_{12} \breve{\mathbf{C}}_{22}^{-1} \right)^\mathsf{T},\end{aligned}$$ holds the of the matrix $\breve{\mathbf{C}}$, $$\begin{aligned} implies subsystem and environment are orthogonal, leading the expression, $$\begin{aligned} \mathbf{0} &= \breve{\mathbf{C}}_\mathcal{E}^\mathsf{T} &= \breve{\mathbf{C}}_{12}^\mathsf{T} \breve{\mathbf{C}}_{11} + \breve{\mathbf{C}}_{22}^\mathsf{T} \breve{\mathbf{C}}_{21} \\ &= \breve{\mathbf{C}}_{22}^{-\mathsf{T}} \breve{\mathbf{C}}_{12}^\mathsf{T} + \breve{\mathbf{C}}_{21} \breve{\mathbf{C}}_{11}^{-1} &= \left( \breve{\mathbf{C}}_{12} \breve{\mathbf{C}}_{22}^{-1} \right)^\mathsf{T} + \breve{\mathbf{C}}_{21} \breve{\mathbf{C}}_{11}^{-1},\end{aligned}$$ from which Eq. immediately follows. of block-diagonalization $\mathbf{Q}$ =============================================================== In the literature,[@Cederbaum89; @Sikkema09; @Peng12; @Seino12] there exist multiple representations of the block-diagonalization matrix $\mathbf{Q}$. representations are very similar to the expression in and it is trivial show that they achieve the same The different arise the different definitions: 1. definition of the block-diagonalization as either $$\begin{aligned} \tilde{\mathbf{F}} = \mathbf{Q} \mathbf{F} \mathbf{Q}^\mathsf{T},\end{aligned}$$ or $$\begin{aligned} \tilde{\mathbf{F}} = \
mathbf{C}}_{11}^{-1} = \left( \breve{\mathbf{C}}_{12} \brEve{\mathbf{C}}_{22}^{-1} \Right)^\MatHsf{t},\eNd{alIgneD}$$ holds by considERing The orthogonality of the eIgenvEcTOr maTRiX $\brevE{\mathbf{c}}$, $$\BeGIN{alIgNeD}
\maThBF{I} = \Breve{\MatHbf{C}}^\matHsf{T} \breve{\mAthBf{c}}.\end{aligned}$$ THIs Implies thaT suBsystem and enVirOnment EiGenVEctorS arE orthOgonal, LEading To the exprEsSIon, $$\begIN{aligneD}
\MAtHbf{0} 
&= \bReve{\mathbf{C}}_\mathcaL{e}^\mAThsf{T} \breve{\mathBf{C}}_\matHcAL{S}\\
&= \BREve{\MatHbf{C}}_{12}^\mathsf{t} \bReve{\mAThbf{C}}_{11} + \brEVe{\MATHbf{c}}_{22}^\Mathsf{T} \breve{\mAthbf{C}}_{21} \\
&= \breve{\MAthBf{C}}_{22}^{-\matHsF{T}} \bREve{\matHbf{C}}_{12}^\mAtHSf{T} + \Breve{\mathbf{c}}_{21} \breVe{\mathbf{C}}_{11}^{-1} \\
&= \Left( \brEVe{\mathbF{c}}_{12} \breve{\mAthbf{C}}_{22}^{-1} \RigHt)^\mAthsF{t} + \bReVe{\mAtHBf{C}}_{21} \BReVe{\mAThbF{C}}_{11}^{-1},\end{aliGnEd}$$ From wHich eQ.  IMMediAteLy foLlows.

representatioN of The bLOck-DiagoNalizAtioN mAtrix $\Mathbf{q}$ {#app:q}
===============================================================

in The literature,[@CeDerbAum89; @SikkemA09; @PeNg12; @seiNo12] There EXist muLtiPle RepreseNtationS Of tHe BLOCk-Diagonalization matRiX $\MAtHbf{Q}$. ThesE repreSEnTaTIons are vErY siMilaR TO the eXpreSSiOn in Eq.  anD it is tRIvIaL to show ThAt they AcHieVe tHe exaCT samE block-DiagonalIzatiON. The different rEPresentations ARiSE FrOM the FolLowing diffeRent DEfinItioNS:

1.  THe dEFinitIon of ThE BlOCk-diagonalization as EiTher $$\beGin{alIgned}
    \tilde{\matHbf{F}} = \mathbf{q} \MAThbf{F} \matHbf{Q}^\MAtHSf{T},\end{aligned}$$ oR $$\begiN{aligned}
    \tiLDe{\mathbf{f}} = \
mathbf{C}}_{11}^{-1} = \le ft( \breve {\mat hbf {C} }_ {12}  \br eve{\mathbf{C} } _{22 }^{-1} \right)^\mathsf {T},\ en d {ali g ne d}$$holds b y  c o n sid er in g t he or thogo nal ity ofthe eigenv ect or  matrix $\br e ve {\mathbf{C }}$ , $$\begin{a lig ned}
\ ma thb f {I} =  \b reve{ \mathb f {C}}^\ mathsf{T}  \ b reve{\ m athbf{C } } .\ end{ aligned}$$ This i m pl i es that subsys tem an de nv i r onm ent  eigenvect or s are orthogo n al , l ead i ng to the exp ression, $$ \ beg in{ali gn ed} \mathb f{0}
& =  \b reve{\mathb f{C} }_\mathca l{E}^\ m athsf{T }  \breve {\math bf{ C}} _\ma t hc al {S} \\ &=\ br eve { \ma thbf{C}} _{ 12 }^\ma thsf { T } \bre ve{ \mat hbf{C }}_{11} + \br eve {\ma t hbf {C}}_ {22}^ \mat hs f{T}\breve {\mat hb f{C}}_{21} \\
& = \b reve{\mat hbf {C }}_ {2 2}^{- \ mathsf {T} } \ breve{\ mathbf{ C }}_ {1 2 } ^ \m athsf{T} + \breve{ \m a t hb f{C}}_{2 1} \br e ve {\ m athbf{C} }_ {11 }^{- 1 }  \\
& = \l e ft ( \breve {\math b f{ C} }_{12}\b reve{\ ma thb f{C }}_{2 2 }^{- 1} \ri ght)^\ma thsf{ T } + \breve{\ma t hbf{C}}_{21}\ br e v e{ \ math bf{ C}}_{11}^{- 1},\ e nd{a lign e d} $$f rom w hichEq .  i mmediately follows. 
Repres entat ion of the bl ock-diagon a l i zation m atri x  $ \ mathbf{Q}$ {#a pp:q} 
========= = ======== ===== ======== ========= = = ======== === === === === = 
I n the literat u r e,[@ Ce derbaum 89;  @Sikke ma0 9;@Pe ng1 2;  @Seino12 ] thereex is tmu lti ple r e presenta ti ons  o f t he bl o ck-dia gonal izat io nm atr ix $\ma t hb f { Q}$.  T he se r epr es entat ions are  very s imilar to  th e  exp re ss ion inEq.  and it i strivial to  s how  thatt h ey achie ve the exact same block - diagona liz ation . Th e differe ntrepres ent a tionsarisefromth e f o l lowin g di ffe re nt definit i o ns: 

1. T he d efiniti on of the block-di a gon alization aseit her$ $ \b egi n {a l ign ed } 
 \tilde{\mathbf{ F}} = \mat hb f {Q } \mathbf{ F } \ ma thbf{Q} ^\maths f{T}, \ end{ali gned}$$ o r $$\begi n{ alig n e d}    \tilde {\mathbf {F}} = \
mathbf{C}}_{11}^{-1} =_\left( \breve{\mathbf{C}}_{12}_\breve{\mathbf{C}}_{22}^{-1} \right)^\mathsf{T},\end{aligned}$$ holds by_considering the_orthogonality_of the_eigenvector_matrix $\breve{\mathbf{C}}$, $$\begin{aligned}
\mathbf{I}_= \breve{\mathbf{C}}^\mathsf{T} \breve{\mathbf{C}}.\end{aligned}$$_This implies that subsystem_and environment eigenvectors_are_orthogonal, leading to the expression, $$\begin{aligned}
\mathbf{0} 
&= \breve{\mathbf{C}}_\mathcal{E}^\mathsf{T} \breve{\mathbf{C}}_\mathcal{S}\\
&= \breve{\mathbf{C}}_{12}^\mathsf{T} \breve{\mathbf{C}}_{11} + \breve{\mathbf{C}}_{22}^\mathsf{T} \breve{\mathbf{C}}_{21}_\\
&=_\breve{\mathbf{C}}_{22}^{-\mathsf{T}} \breve{\mathbf{C}}_{12}^\mathsf{T}_+_\breve{\mathbf{C}}_{21}_\breve{\mathbf{C}}_{11}^{-1} \\
&= \left( \breve{\mathbf{C}}_{12} \breve{\mathbf{C}}_{22}^{-1}_\right)^\mathsf{T} + \breve{\mathbf{C}}_{21} \breve{\mathbf{C}}_{11}^{-1},\end{aligned}$$ from_which Eq. _immediately follows.

Representation of the block-diagonalization matrix $\mathbf{Q}$ {#app:q}
===============================================================

In_the_literature,[@Cederbaum89; @Sikkema09; @Peng12;_@Seino12] there exist multiple representations of the block-diagonalization matrix_$\mathbf{Q}$. These representations are very similar_to the expression_in_Eq. _and it is trivial_to show that they achieve the_exact same block-diagonalization. The different representations_arise from the following different definitions:

1. _The definition of the block-diagonalization as_either $$\begin{aligned}
   _\tilde{\mathbf{F}} =_\mathbf{Q} \mathbf{F} \mathbf{Q}^\mathsf{T},\end{aligned}$$ or $$\begin{aligned}
_  _\tilde{\mathbf{F}} =_\