}}\omega_3 - {\mathrm{i}}\omega_4)
 \times \nonumber \\ & & \hspace{1em}
   \mathcal{G}_{\gamma\gamma'}({\mathrm{i}}\omega_3) \mathcal{G}_{\delta\delta'}({\mathrm{i}}\omega_4)
   V^{\mathrm{rn},-1}_{\gamma' i} V^{\mathrm{rn},-1}_{\delta' j}.
   \label{eq:TwoParticleGF:connected:frequency:initial}\end{aligned}$$ In order to solve the frequency integral analytically, we again transform into the eigenbasis of $H_0 + \Sigma$. Eq. (\[eq:TwoParticleGF:connected:frequency:initial\]) now reads $$\begin{aligned}
 & & \mathcal{G}^{(2,c)}_{ijij}({\mathrm{i}}\omega_1, {\mathrm{i}}\omega_2, {\mathrm{i}}\omega_3, {\mathrm{i}}\omega_4) =
   - 2\pi \sum_{\alpha\beta\gamma\delta} \sum_{\tilde\mu\tilde\nu\tilde\rho\tilde\sigma}
     V^{\mathrm{ri}}_{i\tilde\mu} V^{\mathrm{ri}}_{j\tilde\nu}
 \times \nonumber \\ & & \hspace{1em}
   \mathcal{G}_{\tilde\mu\tilde\mu}({\mathrm{i}}\omega_1) \mathcal{G}_{\tilde\nu\tilde\nu}({\mathrm{i}}\omega_2)
    V^{\mathrm{in}}_{\tilde\mu\alpha} V^{\mathrm{in}}_{\tilde\nu\beta}
 \times \nonumber \\ & & \hspace{1em}
   \Gamma_{\alpha\beta\gamma\delta} \delta({\mathrm{i}}\omega_1 + {\mathrm{i}}\omega_2 - {\mathrm{i}}\omega_3 - {\mathrm{i}}\omega_4)
 \times \nonumber \\ & & \hspace{1em}
   V^{\mathrm{in},-1}_{\gamma\tilde\rho} V^{\mathrm{in},-1}_{\delta\tilde\sigma}
   \mathcal{G}_{\tilde\rho\tilde\rho}({\mathrm{i}}\omega_3) \mathcal{G}_{\tilde\sigma\tilde\sigma}({\mathrm{i
} } \omega_3 - { \mathrm{i}}\omega_4) 
  \times \nonumber \\ & & \hspace{1em } 
    \mathcal{G}_{\gamma\gamma'}({\mathrm{i}}\omega_3) \mathcal{G}_{\delta\delta'}({\mathrm{i}}\omega_4) 
    V^{\mathrm{rn},-1}_{\gamma' i } V^{\mathrm{rn},-1}_{\delta' j }. 
    \label{eq: TwoParticleGF: connected: frequency: initial}\end{aligned}$$ In order to resolve the frequency built-in analytically, we again transform into the eigenbasis of $ H_0 + \Sigma$. Eq.   (\[eq: TwoParticleGF: connected: frequency: initial\ ]) nowadays read $ $ \begin{aligned } 
  & & \mathcal{G}^{(2,c)}_{ijij}({\mathrm{i}}\omega_1, { \mathrm{i}}\omega_2, { \mathrm{i}}\omega_3, { \mathrm{i}}\omega_4) = 
    - 2\pi \sum_{\alpha\beta\gamma\delta } \sum_{\tilde\mu\tilde\nu\tilde\rho\tilde\sigma } 
      V^{\mathrm{ri}}_{i\tilde\mu } V^{\mathrm{ri}}_{j\tilde\nu } 
  \times \nonumber \\ & & \hspace{1em } 
    \mathcal{G}_{\tilde\mu\tilde\mu}({\mathrm{i}}\omega_1) \mathcal{G}_{\tilde\nu\tilde\nu}({\mathrm{i}}\omega_2) 
     V^{\mathrm{in}}_{\tilde\mu\alpha } V^{\mathrm{in}}_{\tilde\nu\beta } 
  \times \nonumber \\ & & \hspace{1em } 
    \Gamma_{\alpha\beta\gamma\delta } \delta({\mathrm{i}}\omega_1 + { \mathrm{i}}\omega_2 - { \mathrm{i}}\omega_3 - { \mathrm{i}}\omega_4) 
  \times \nonumber \\ & & \hspace{1em } 
    V^{\mathrm{in},-1}_{\gamma\tilde\rho } V^{\mathrm{in},-1}_{\delta\tilde\sigma } 
    \mathcal{G}_{\tilde\rho\tilde\rho}({\mathrm{i}}\omega_3) \mathcal{G}_{\tilde\sigma\tilde\sigma}({\mathrm{i
}}\omeha_3 - {\mathrm{i}}\omega_4)
 \times \nunumber \\ & & \hspaew{1em}
   \kathcam{G}_{\gamma\gxmma'}({\mathrm{i}}\omega_3) \mathcal{G}_{\delva\deota'}({\maukrm{i}}\omega_4)
   V^{\mathrm{rv},-1}_{\gamma' i} N^{\mathrm{rn},-1}_{\eelte' j}.
   \label{eq:Two'zrticleNY:connsgted:fxewuency:initial}\ekd{aligned}$$ It order to solee tke frequency integral analytically, wq again tgansform into jhe epgqnbaaps of $H_0 + \Sigma$. Eq. (\[eq:TwoParticleGR:connecued:frequency:initiak\]) now reads $$\begin{aligned}
 & & \matjcal{G}^{(2,c)}_{ijij}({\mathrm{i}}\olega_1, {\mathrm{u}}\omedq_2, {\mathrm{i}}\omeea_3, {\mathrm{i}}\omega_4) =
   - 2\pi \sum_{\alpha\beta\gamma\delta} \sum_{\tildd\mu\tikde\nu\tilde\tki\tipge\sigma}
     V^{\matrrm{ri}}_{i\tilde\mm} V^{\mathsm{ri}}_{j\tikde\nu}
 \times \nokumbec \\ & & \hspace{1em}
   \mathcal{G}_{\tmlde\mu\tilde\mu}({\mathrm{i}}\jmega_1) \matvccl{G}_{\tilde\nu\tilde\nu}({\mathem{u}}\omegd_2)
    V^{\magyrm{kn}}_{\tjlve\mh\alpha} V^{\methrm{in}}_{\tilds\nu\beta}
 \timws \nonumber \\ & & \hspave{1qn}
   \Gamma_{\alphz\beta\gwmia\delta} \delta({\mathrm{i}}\omega_1 + {\mathrm{i}}\omega_2 - {\mzthrm{i}}\omega_3 - {\mathrm{i}}\ometa_4)
 \times \nonumber \\ & & \jspace{1em}
   V^{\mathrm{in},-1}_{\gamma\tilde\rho} V^{\mathrm{in},-1}_{\delta\tilde\sigma}
   \mavhzal{Y}_{\bllde\fyo\hilde\rho}({\mathrm{i}}\omega_3) \mathcal{G}_{\tilde\sigma\tilde\syfms}({\msthrm{i
}}\omega_3 - {\mathrm{i}}\omega_4) \times \nonumber \\ & \mathcal{G}_{\gamma\gamma'}({\mathrm{i}}\omega_3) V^{\mathrm{rn},-1}_{\gamma' i} j}. \label{eq:TwoParticleGF:connected:frequency:initial}\end{aligned}$$ In integral we again transform the eigenbasis of + \Sigma$. Eq. (\[eq:TwoParticleGF:connected:frequency:initial\]) now reads & & \mathcal{G}^{(2,c)}_{ijij}({\mathrm{i}}\omega_1, {\mathrm{i}}\omega_2, {\mathrm{i}}\omega_3, {\mathrm{i}}\omega_4) = - 2\pi \sum_{\alpha\beta\gamma\delta} \sum_{\tilde\mu\tilde\nu\tilde\rho\tilde\sigma} V^{\mathrm{ri}}_{i\tilde\mu} V^{\mathrm{ri}}_{j\tilde\nu} \nonumber \\ & & \hspace{1em} \mathcal{G}_{\tilde\mu\tilde\mu}({\mathrm{i}}\omega_1) \mathcal{G}_{\tilde\nu\tilde\nu}({\mathrm{i}}\omega_2) V^{\mathrm{in}}_{\tilde\mu\alpha} V^{\mathrm{in}}_{\tilde\nu\beta} \times \nonumber \\ & \hspace{1em} \delta({\mathrm{i}}\omega_1 {\mathrm{i}}\omega_2 {\mathrm{i}}\omega_3 - {\mathrm{i}}\omega_4) \times \nonumber \\ & & \hspace{1em} V^{\mathrm{in},-1}_{\gamma\tilde\rho} V^{\mathrm{in},-1}_{\delta\tilde\sigma} \mathcal{G}_{\tilde\rho\tilde\rho}({\mathrm{i}}\omega_3) \mathcal{G}_{\tilde\sigma\tilde\sigma}({\mathrm{i
}}\omega_3 - {\mathrm{i}}\omega_4)
 \times \nonUmber \\ & & \hspacE{1em}
   \maThcAl{G}_{\GaMma\gAmma'}({\Mathrm{i}}\omega_3) \maTHcal{g}_{\delta\delta'}({\mathrm{i}}\omegA_4)
   V^{\matHrM{Rn},-1}_{\gaMMa' I} V^{\matHrm{rn},-1}_{\deLTa' J}.
   \LAbeL{eQ:TWoPArTIcLeGF:cOnnEcted:frEquency:iniTiaL}\eNd{aligned}$$ In oRDeR to solve thE frEquency integRal AnalytIcAllY, We agaIn tRansfOrm intO The eigEnbasis of $h_0 + \SIGma$. Eq. (\[eQ:twoPartICLegF:coNnected:frequency:iNItIAl\]) now reads $$\begiN{alignEd}
 & & \MAtHCAl{G}^{(2,C)}_{ijIj}({\mathrm{i}}\oMeGa_1, {\matHRm{i}}\omegA_2, {\MaTHRM{i}}\oMEga_3, {\mathrm{i}}\omeGa_4) =
   - 2\pi \sum_{\alphA\BetA\gamma\DeLta} \SUm_{\tildE\mu\tiLdE\Nu\tIlde\rho\tildE\sigMa}
     V^{\mathrm{Ri}}_{i\tilDE\mu} V^{\matHRm{ri}}_{j\tiLde\nu}
 \tImeS \noNumbER \\ & & \hSpAce{1Em}
   \MAthCAl{g}_{\tiLDe\mU\tilde\mu}({\MaThRm{i}}\omEga_1) \mATHCAl{G}_{\tIldE\nu\tIlde\nU}({\mathrm{i}}\omega_2)
    v^{\maThrm{IN}}_{\tiLde\mu\Alpha} v^{\matHrM{in}}_{\tiLde\nu\bEta}
 \tiMeS \nonumber \\ & & \hspace{1Em}
   \GaMma_{\alpha\bEta\GaMma\DeLta} \deLTa({\mathRm{i}}\OmeGa_1 + {\mathrM{i}}\omega_2 - {\MAthRm{I}}\OMEgA_3 - {\mathrm{i}}\omega_4)
 \times \NoNUMbEr \\ & & \hspace{1Em}
   V^{\matHRm{In},-1}_{\GAmma\tildE\rHo} V^{\MathRM{In},-1}_{\delTa\tiLDe\Sigma}
   \matHcal{G}_{\tILdE\rHo\tilde\RhO}({\mathrM{i}}\OmeGa_3) \mAthcaL{g}_{\tilDe\sigmA\tilde\siGma}({\maTHrm{i
}}\omega_3 - {\mathrm{i}}\ omega_4)
\time s \ non um ber\\ &  & \hspace{1em } 
  \mathcal{G}_{\gamma\ga mma'} ({ \ math r m{ i}}\o mega_3) \m a t hca l{ G} _{\ de l ta \delt a'} ({\math rm{i}}\ome ga_ 4) 
   V^{\math r m{ rn},-1}_{\ gam ma' i} V^{\m ath rm{rn} ,- 1}_ { \delt a'j}.
  \lab e l{eq:T woParticl eG F :conne c ted:fre q u en cy:i nitial}\end{align e d} $ $ In order tosolveth e  f r e que ncy  integralan alyti c ally, w e  a g a i n t r ansform intothe eigenba s isof $H_ 0+ \ S igma$.  Eq.(\ [ eq: TwoParticle GF:c onnected: freque n cy:init i al\]) n ow rea ds$$\ begi n {a li gne d}  &&  \ mat h cal {G}^{(2, c) }_ {ijij }({\ m a t h rm{i }}\ omeg a_1,{\mathrm{i}}\ ome ga_2 ,  {\ mathr m{i}} \ome ga _3, { \mathr m{i}} \o mega_4) =
   -2\pi  \sum_{\a lph a\ bet a\ gamma \ delta}  \s um_ {\tilde \mu\til d e\n u\ t i l de \rho\tilde\sigma}  V ^{\mathr m{ri}} _ {i \t i lde\mu}V^ {\m athr m { ri}}_ {j\t i ld e\nu}
 \ times\ no nu mber \\  &  & \hs pa ce{ 1em }
  \ math cal{G} _{\tilde \mu\t i lde\mu}({\math r m{i}}\omega_1 )  \ m a th c al{G }_{ \tilde\nu\t ilde \ nu}( {\ma t hr m{i } }\ome ga_2) 
 V ^{\mathrm{in}}_{\ti ld e\mu\a lpha}  V^{\mathrm{i n}}_{\tild e \ n u\beta} \ti m es \nonumber \\ &  & \h space{1em}    \Gamm a_{\a lpha\bet a\gamma\d e l ta} \del ta( {\m ath rm{ i } }\ omega_1 + {\m a t hrm{ i} }\omega _2- {\mat hrm {i} }\o meg a_ 3 - {\mat hrm{i}}\ om eg a_ 4) 
 \ times \nonumbe r\\&& \ hspac e {1em}   V^ {\ma th rm { in} ,-1}_{\ g am m a \til de \r ho}V^{ \m athrm {in} , -1} _{\delt a\tilde\s igm a }
  \ ma thcal{G }_{\tilde\rho \t ilde\rho}( {\ mat hrm{i} } \ omega_3)  \mathcal{G}_{\tilde\si g ma\tild e\s igma} ({\m athrm{i
}}\omega_3 -_{\mathrm{i}}\omega_4)
 \times_\nonumber \\ & &_\hspace{1em}
 __\mathcal{G}_{\gamma\gamma'}({\mathrm{i}}\omega_3) \mathcal{G}_{\delta\delta'}({\mathrm{i}}\omega_4)
__ V^{\mathrm{rn},-1}_{\gamma' i}_V^{\mathrm{rn},-1}_{\delta' j}.
 _ \label{eq:TwoParticleGF:connected:frequency:initial}\end{aligned}$$ In order_to solve the_frequency_integral analytically, we again transform into the eigenbasis of $H_0 + \Sigma$. Eq. (\[eq:TwoParticleGF:connected:frequency:initial\]) now_reads_$$\begin{aligned}
 &_&_\mathcal{G}^{(2,c)}_{ijij}({\mathrm{i}}\omega_1,_{\mathrm{i}}\omega_2, {\mathrm{i}}\omega_3, {\mathrm{i}}\omega_4) =
 _ - 2\pi \sum_{\alpha\beta\gamma\delta} \sum_{\tilde\mu\tilde\nu\tilde\rho\tilde\sigma}
_ _  V^{\mathrm{ri}}_{i\tilde\mu} V^{\mathrm{ri}}_{j\tilde\nu}
 \times \nonumber \\ &_&_\hspace{1em}
  _\mathcal{G}_{\tilde\mu\tilde\mu}({\mathrm{i}}\omega_1) \mathcal{G}_{\tilde\nu\tilde\nu}({\mathrm{i}}\omega_2)
    V^{\mathrm{in}}_{\tilde\mu\alpha} V^{\mathrm{in}}_{\tilde\nu\beta}
 \times \nonumber_\\ & & \hspace{1em}
  _\Gamma_{\alpha\beta\gamma\delta} \delta({\mathrm{i}}\omega_1 +_{\mathrm{i}}\omega_2_-_{\mathrm{i}}\omega_3 - {\mathrm{i}}\omega_4)
 \times_\nonumber \\ & & \hspace{1em}
 _ V^{\mathrm{in},-1}_{\gamma\tilde\rho} V^{\mathrm{in},-1}_{\delta\tilde\sigma}
   \mathcal{G}_{\tilde\rho\tilde\rho}({\mathrm{i}}\omega_3)_\mathcal{G}_{\tilde\sigma\tilde\sigma}({\mathrm{i