} -1}) G_0^{\rm R} \partial_{X^0} G_0^{{\rm R} -1}]
  + {{\rm c.c.}}\label{eq:mp3b} \\
  & \mp \frac{i \hbar^2}{4} \delta_{\hat k}^{\phantom{\hat k} {\hat \imath}} \int \frac{d^D \pi}{(2 \pi \hbar)^D} \int d X^0
  \operatorname{tr}[{\hat \Sigma}_1 {\hat G}_0]^< + {{\rm c.c.}}\label{eq:mp3c} \\
  & \pm \frac{i \hbar^2}{4} \int \frac{d^D \pi}{(2 \pi \hbar)^D} \int d X^0 f^{\prime}(-\pi_{\hat 0}) \pi_{\hat k}
  \operatorname{tr}[\Sigma_1^{< (1)} (G_0^{\rm R} - G_0^{\rm A}) \partial_{\pi_{\hat \imath}} G_0^{{\rm R} -1} G_0^{\rm R}] + {{\rm c.c.}}, \label{eq:mp3d}\end{aligned}$$

\[eq:mp3\]

where $\epsilon_{ABC}$ in Eq.  is the antisymmetric tensor with $\epsilon_{\pi_{\hat 0} \pi_{\hat \imath} X^0} = 1$. $f(\omega) = (e^{\beta \omega} \mp 1)^{-1}$ is the distribution function for boson or fermion.

Below we focus on the clean and noninteracting limit ${{\hat {\bm \Sigma}}} = 0$ to derive the Berry-phase formula in the Bloch basis. The retarded Green’s function is given by $G_0^{\rm R} = [(-\pi_{\hat 0}) - {\cal H}(X^0) + \mu + i \eta]^{-1}$ with $\eta \to +0$, leading to $\partial_{\pi_{\hat 0}} G_0^{{\rm R} -1} = -1$, $\partial_{\pi_{\hat \imath}} G_0^{{\rm R} -1} = -v^{\hat \imath}$, and $\partial_{X^0} G_0^{{\rm R} -
} -1 }) G_0^{\rm R } \partial_{X^0 } G_0^{{\rm R } -1 } ] 
   + { { \rm c.c.}}\label{eq: mp3b } \\ 
   & \mp \frac{i \hbar^2}{4 } \delta_{\hat k}^{\phantom{\hat k } { \hat \imath } } \int \frac{d^D \pi}{(2 \pi \hbar)^D } \int d X^0 
   \operatorname{tr}[{\hat \Sigma}_1 { \hat G}_0]^ < + { { \rm c.c.}}\label{eq: mp3c } \\ 
   & \pm \frac{i \hbar^2}{4 } \int \frac{d^D \pi}{(2 \pi \hbar)^D } \int d X^0 f^{\prime}(-\pi_{\hat 0 }) \pi_{\hat k } 
   \operatorname{tr}[\Sigma_1^ { < (1) } (G_0^{\rm radius } - G_0^{\rm A }) \partial_{\pi_{\hat \imath } } G_0^{{\rm R } -1 } G_0^{\rm radius } ] + { { \rm c.c. } }, \label{eq: mp3d}\end{aligned}$$ 

 \[eq: mp3\ ] 

 where $ \epsilon_{ABC}$ in Eq.    is the antisymmetric tensor with $ \epsilon_{\pi_{\hat 0 } \pi_{\hat \imath } X^0 } = 1$. $ f(\omega) = (e^{\beta \omega } \mp 1)^{-1}$ is the distribution function for boson or fermion. 

 Below we focus on the clean and noninteracting terminus ad quem $ { { \hat { \bm \Sigma } } } = 0 $ to derive the Berry - phase formula in the Bloch basis. The retarded Green ’s affair is given by $ G_0^{\rm R } = [ (-\pi_{\hat 0 }) - { \cal H}(X^0) + \mu + i \eta]^{-1}$ with $ \eta \to +0 $, contribute to $ \partial_{\pi_{\hat 0 } } G_0^{{\rm R } -1 } = -1 $, $ \partial_{\pi_{\hat \imath } } G_0^{{\rm R } -1 } = -v^{\hat \imath}$, and $ \partial_{X^0 } G_0^{{\rm R } -
} -1}) G_0^{\gm R} \partial_{X^0} G_0^{{\rm R} -1}]
  + {{\ro c.c.}}\label{eq:mp3b} \\
  & \mp \fcac{i \hbzr^2}{4} \delta_{\fat k}^{\phantom{\hat k} {\hat \imath}} \mnt \drac{d^E \pi}{(2 \pi \hbar)^D} \int d X^0
  \operatorjame{tr}[{\har \Sijma}_1 {\hat G}_0]^< + {{\rm c.c.}}\label{eq:mi3e} \\
  & \ll \frcc{m \hbar^2}{4} \int \frac{c^D \pi}{(2 \pi \htar)^D} \int d X^0 f^{\[rkmz}(-\pi_{\hat 0}) \pi_{\hat k}
  \operatorname{tr}[\Sigma_1^{< (1)} (G_0^{\rm R} - H_0^{\rm A}) \partial_{\py_{\hat \ymatg}} G_0^{{\rm R} -1} G_0^{\rm R}] + {{\rm c.c.}}, \label{eq:mp3d}\ens{alignev}$$

\[eq:mp3\]

where $\epsikon_{ABC}$ in Eq.  is the antisylmetgic tensor with $\epdilon_{\pi_{\hat 0} \pi_{\hwr \imath} X^0} = 1$. $w(\omega) = (e^{\bttc \omega} \mp 1)^{-1}$ is the distribution function fur boxon or ferniin.

Bfnow we focuw on nhe clean and nonintesacting limit ${{\hat {\bm \Xigka}}} = 0$ to derive the Berry-'hase formula in the Bloch bavia. The retarded Grweb’s futctimn ir gixen bb $G_0^{\dm R} = [(-\oi_{\het 0}) - {\cal H}(X^0) + \mu + i \eta]^{-1}$ with $\eta \to +0$, leadimg np $\partial_{\pi_{\hzt 0}} G_0^{{\ri W} -1} = -1$, $\partial_{\pi_{\hat \imath}} G_0^{{\rm R} -1} = -v^{\hat \imdth}$, and $\partial_{X^0} G_0^{{\rm R} -
} -1}) G_0^{\rm R} \partial_{X^0} G_0^{{\rm R} {{\rm \\ & \frac{i \hbar^2}{4} \delta_{\hat \frac{d^D \pi \hbar)^D} \int X^0 \operatorname{tr}[{\hat \Sigma}_1 G}_0]^< + {{\rm c.c.}}\label{eq:mp3c} \\ & \frac{i \hbar^2}{4} \int \frac{d^D \pi}{(2 \pi \hbar)^D} \int d X^0 f^{\prime}(-\pi_{\hat 0}) \pi_{\hat \operatorname{tr}[\Sigma_1^{< (1)} (G_0^{\rm R} - G_0^{\rm A}) \partial_{\pi_{\hat \imath}} G_0^{{\rm R} -1} G_0^{\rm + c.c.}}, \[eq:mp3\] $\epsilon_{ABC}$ in Eq. is the antisymmetric tensor with $\epsilon_{\pi_{\hat 0} \pi_{\hat \imath} X^0} = 1$. $f(\omega) (e^{\beta \omega} \mp 1)^{-1}$ is the distribution function boson or fermion. Below focus on the clean and limit {\bm \Sigma}}} 0$ derive Berry-phase formula in Bloch basis. The retarded Green’s function is given by $G_0^{\rm R} = [(-\pi_{\hat 0}) - {\cal H}(X^0) \mu + with $\eta +0$, to 0}} G_0^{{\rm R} -1$, $\partial_{\pi_{\hat \imath}} G_0^{{\rm R} -1} and $\partial_{X^0} G_0^{{\rm R} -
} -1}) G_0^{\rm R} \partial_{X^0} G_0^{{\rm R} -1}]
  + {{\rm c.c.}}\label{Eq:mp3b} \\
  & \mp \fraC{i \hbaR^2}{4} \deLta_{\HaT k}^{\phAntoM{\hat k} {\hat \imath}} \iNT \fraC{d^D \pi}{(2 \pi \hbar)^D} \int d X^0
  \operaTornaMe{TR}[{\hat \sIgMa}_1 {\hat g}_0]^< + {{\rm c.c.}}\laBEl{EQ:Mp3c} \\
  & \Pm \FrAc{i \HbAR^2}{4} \iNt \fraC{d^D \Pi}{(2 \pi \hbaR)^D} \int d X^0 f^{\prIme}(-\Pi_{\Hat 0}) \pi_{\hat k}
  \opeRAtOrname{tr}[\SiGma_1^{< (1)} (g_0^{\rm R} - G_0^{\rm A}) \partIal_{\Pi_{\hat \iMaTh}} G_0^{{\RM R} -1} G_0^{\rm r}] + {{\rm C.c.}}, \labEl{eq:mp3D}\End{aliGned}$$

\[eq:mp3\]

wHeRE $\epsilON_{ABC}$ in EQ.  IS tHe anTisymmetric tensor WItH $\Epsilon_{\pi_{\hat 0} \pi_{\Hat \imaTh} x^0} = 1$. $F(\oMEGa) = (e^{\BetA \omega} \mp 1)^{-1}$ is ThE distRIbution FUnCTIOn fOR boson or fermiOn.

Below we foCUs oN the clEaN anD NonintEractInG LimIt ${{\hat {\bm \SigmA}}} = 0$ to dErive the BErry-phASe formuLA in the BLoch baSis. the RetaRDeD GReeN’s FUncTIoN is GIveN by $G_0^{\rm R} = [(-\pI_{\hAt 0}) - {\Cal H}(X^0) + \Mu + i \eTA]^{-1}$ WITh $\etA \to +0$, LeadIng to $\Partial_{\pi_{\hat 0}} G_0^{{\Rm R} -1} = -1$, $\PartIAl_{\pI_{\hat \iMath}} G_0^{{\Rm R} -1} = -v^{\HaT \imatH}$, and $\paRtial_{x^0} G_0^{{\Rm R} -
} -1}) G_0^{\rm R} \partia l_{X^0} G_ 0^{{\ rmR}-1 }]
 + { {\rm c.c.}}\la b el{e q:mp3b} \\
  & \mp \fr ac{i\h b ar^2 } {4 } \de lta_{\h a tk } ^{\ ph an tom {\ h at  k} { \ha t \imat h}} \int \ fra c{ d^D \pi}{(2\ pi  \hbar)^D}  \i nt d X^0
  \ ope ratorn am e{t r }[{\h at\Sigm a}_1 { \ hat G} _0]^< + { {\ r m c.c. } }\label { e q: mp3c } \\
  & \pm \fra c {i \hbar^2}{4} \i nt \fr ac { d^ D \pi }{( 2 \pi \hba r) ^D} \ i nt d X^ 0  f ^ { \ pri m e}(-\pi_{\hat  0}) \pi_{\ h atk}
  \ op era t orname {tr}[ \S i gma _1^{< (1)}(G_0 ^{\rm R}- G_0^ { \rm A}) \partia l_{\pi _{\ hat  \im a th }}  G_ 0^ { {\r m  R } - 1 } G _0^{\rmR} ]+ {{\ rm c . c . } }, \ lab el{e q:mp3 d}\end{aligne d}$ $

\ [ eq: mp3\] 

whe re $ \e psilo n_{ABC }$ in  E q.  is the anti symm etric ten sor  w ith  $ \epsi l on_{\p i_{ \ha t 0} \p i_{\hat \im at h } X^ 0} = 1$. $f(\omega )= (e ^{\beta\omega }  \ mp 1)^{-1}$  i s t he d i s tribu tion fu nction f or bos o nor  fermio n. 

Belo wwefoc us on thecleanand noni ntera c ting limit ${{ \ hat {\bm \Sig m a} } }  = 0$ t o d erive the B erry - phas e fo r mu lai n the  Bloc hb as i s. The retarded Gre en ’s fun ction  is given by$G_0^{\rmR } = [(-\pi _{\h a t0 }) - {\cal H}( X^0)+ \mu + i\ eta]^{-1 }$ wi th $\eta  \to +0$, l eading t o $ \pa rti al_ { \ pi _{\hat 0}} G_ 0 ^ {{\r mR} -1}= - 1$, $\p art ial _{\ pi_ {\ hat \imat h}} G_0^ {{ \r mR}  -1 } = - v ^{\hat \ im ath }$ , a nd $\ p artial _{X^0 } G_ 0^ {{ \ rmR} -
} -1})_G_0^{\rm R}_\partial_{X^0} G_0^{{\rm R} -1}]
_ +_{{\rm_c.c.}}\label{eq:mp3b} \\
__& \mp \frac{i_\hbar^2}{4} \delta_{\hat k}^{\phantom{\hat_k} {\hat \imath}} \int_\frac{d^D \pi}{(2 \pi_\hbar)^D}_\int d X^0
  \operatorname{tr}[{\hat \Sigma}_1 {\hat G}_0]^< + {{\rm c.c.}}\label{eq:mp3c} \\
  &_\pm_\frac{i \hbar^2}{4}_\int_\frac{d^D_\pi}{(2 \pi \hbar)^D} \int d_X^0 f^{\prime}(-\pi_{\hat 0}) \pi_{\hat k}
_ \operatorname{tr}[\Sigma_1^{<_(1)} (G_0^{\rm R} - G_0^{\rm A}) \partial_{\pi_{\hat \imath}}_G_0^{{\rm_R} -1} G_0^{\rm_R}] + {{\rm c.c.}}, \label{eq:mp3d}\end{aligned}$$

\[eq:mp3\]

where $\epsilon_{ABC}$ in Eq.  is_the antisymmetric tensor with $\epsilon_{\pi_{\hat 0}_\pi_{\hat \imath} X^0}_=_1$._$f(\omega) = (e^{\beta \omega}_\mp 1)^{-1}$ is the distribution function_for boson or fermion.

Below we focus_on the clean and noninteracting limit ${{\hat_{\bm \Sigma}}} = 0$ to derive_the Berry-phase formula in the_Bloch basis._The retarded Green’s function is_given by $G_0^{\rm_R} =_[(-\pi_{\hat 0}) -_{\cal H}(X^0) + \mu + i_\eta]^{-1}$ with $\eta_\to +0$, leading to $\partial_{\pi_{\hat 0}}_G_0^{{\rm_R} -1} =_-1$,_$\partial_{\pi_{\hat_\imath}} G_0^{{\rm_R} -1} =_-v^{\hat_\imath}$, and_$\partial_{X^0}_G_0^{{\rm R} -