}=[E{\bf S} - {\bf H} - {\bm \Sigma}_L - {\bm \Sigma}_R]^{-1}$, where in the local orbital basis $\{ \phi_i \}$ Hamiltonian matrix ${\bf H}$ is composed of elements $H_{ij} = \langle \phi_i |\hat{H}_{\rm KS}| \phi_{j} \rangle$ ($\hat{H}_{\rm KS}$ is the effective Kohn-Scham Hamiltonian obtained from the DFT self-consistent loop) and the overlap matrix ${\bf S}$ has elements $S_{ij} = \langle \phi_i | \phi_j \rangle$.

The transmission function Eq. (\[eq:telectron\]) obtained within the NEGF-DFT framework allows us to compute the following integrals [@Esfarjani2006] $$\label{eq:kintegral}
K_n(\mu) = \frac{2}{h} \int\limits_{-\infty}^{\infty} dE\, \mathcal{T}_{\rm el}(E)  (E - \mu)^n \left(-\frac{\partial f(E,\mu)}{\partial E} \right),$$ where is the Fermi-Dirac distribution function at the chemical potential $\mu$. The knowledge of $K_n(\mu)$ finally yields all electronic quantities in the expression for $ZT$: $G=e^2K_0(\mu)$; $S=K_1(\mu)/[eTK_0(\mu)]$; and $\kappa_{\rm el} = \{K_2(\mu) - [K_1(\mu)]^2/K_0(\mu)\}/T$.

The phonon thermal conductance is obtained from the phonon transmission function $\mathcal{T}_{\rm ph}(\omega)$ using the corresponding Landauer-type formula for the scattering region (molecule + portion of electrodes) attached to two semi-infinite electrodes: $$\label{eq:tphonon}
\kappa_{\rm ph} = \frac{\hbar^2}{2\pi k_B T^2} \int\limits_{0}^{\infty} d\omega\, \omega^2 \mathcal{T}_{\rm ph}(\omega) \frac{ e^{\hbar\omega/k_BT}}{(e^{\hbar\omega/k_BT}-1)^2}.$$ The phonon transmission function $\mathcal{T}_{\rm ph} (\omega)$ can be calculated
} =[ E{\bf S } - { \bf H } - { \bm \Sigma}_L - { \bm \Sigma}_R]^{-1}$, where in the local orbital basis $ \ { \phi_i \}$ Hamiltonian matrix $ { \bf H}$ is composed of elements $ H_{ij } = \langle \phi_i |\hat{H}_{\rm KS}| \phi_{j } \rangle$ ($ \hat{H}_{\rm KS}$ is the effective Kohn - Scham Hamiltonian obtained from the DFT self - reproducible loop topology) and the overlap matrix $ { \bf S}$ has elements $ S_{ij } = \langle \phi_i | \phi_j \rangle$. 

 The infection function Eq.   (\[eq: telectron\ ]) obtained within the NEGF - DFT model allows us to compute the following integral   [ @Esfarjani2006 ] $ $ \label{eq: kintegral } 
 K_n(\mu) = \frac{2}{h } \int\limits_{-\infty}^{\infty } dE\, \mathcal{T}_{\rm el}(E)   (E - \mu)^n \left(-\frac{\partial f(E,\mu)}{\partial E } \right),$$ where is the Fermi - Dirac distribution affair at the chemical potential $ \mu$. The cognition of $ K_n(\mu)$ finally yields all electronic quantities in the expression for $ ZT$: $ G = e^2K_0(\mu)$; $ mho = K_1(\mu)/[eTK_0(\mu)]$; and $ \kappa_{\rm el } = \{K_2(\mu) - [ K_1(\mu)]^2 / K_0(\mu)\}/T$. 

 The phonon thermal conductance is obtained from the phonon transmittance function $ \mathcal{T}_{\rm ph}(\omega)$ using the comparable Landauer - type formula for the scattering area (molecule + portion of electrodes) attached to two semi - infinite electrode: $ $ \label{eq: tphonon } 
 \kappa_{\rm ph } = \frac{\hbar^2}{2\pi k_B T^2 } \int\limits_{0}^{\infty } d\omega\, \omega^2 \mathcal{T}_{\rm ph}(\omega) \frac { e^{\hbar\omega / k_BT}}{(e^{\hbar\omega / k_BT}-1)^2}.$$ The phonon transmission routine $ \mathcal{T}_{\rm ph } (\omega)$ can be calculated
}=[E{\bf S} - {\bf H} - {\bm \Sigma}_L - {\bm \Rigma}_R]^{-1}$, where in the loral orbjtal basks $\{ \phi_i \}$ Hamiltonian matrix ${\bd H}$ iw composed of elements $H_{ij} = \lanhle \phi_i |\hat{I}_{\rm KS}| \phi_{j} \rangle$ ($\hat{H}_{\rm KS}$ ia the xffective Kohn-Sgham Hamiltmnian obtained ffol the DFT self-consistent loop) and tre overkao matrix ${\bf S}$ ras tleienta $S_{ij} = \langle \phi_i | \phi_j \rangle$.

The transmpssion function Ea. (\[eq:telectron\]) obtained withln tje NEGF-DFT framewogk allows uw to xompute the wollowing putegrals [@Esfzrjani2006] $$\label{eq:kintegral}
K_n(\mu) = \frxc{2}{h} \iut\limits_{-\infjv}^{\unfhi} dE\, \mathcal{V}_{\rm el}(V)  (E - \mu)^n \lefb(-\grac{\pastial f(R,\mu)}{\partial E} \rlght),$$ xherw is the Fermi-Dirac dmstribution function at the cvejical potential $\my$. Rhe ktowladge if $Y_n(\mh)$ hinzlly ylelvs all elecfronic quanrities in the exprexsyin for $ZT$: $G=e^2K_0(\ju)$; $S=K_1(\mt)/[eEK_0(\mu)]$; and $\kappa_{\rm el} = \{K_2(\mu) - [K_1(\mu)]^2/K_0(\mu)\}/T$.

The phmnoh thermal conductance iw obtained from the pjonon tragsmission function $\mathcal{T}_{\rm ph}(\omega)$ using the cmrres'ovdiun Uqnfauer-type formula for the scattering region (mjmevuke + portion of electrodes) ahtsshed to two sgmi-infiujts electrodes: $$\label{fq:tphonjn}
\kappa_{\rm ph} = \srac{\nbar^2}{2\pi k_B T^2} \int\limits_{0}^{\infty} e\omega\, \omega^2 \nathcal{T}_{\rm ph}(\omega) \frac{ e^{\hbar\umegs/k_BT}}{(e^{\nbar\omega/k_BT}-1)^2}.$$ The phonon trahsmission fknction $\mzghcal{T}_{\rm ph} (\omegx)$ csn be calculated
}=[E{\bf S} - {\bf H} - {\bm {\bm where in local orbital basis ${\bf is composed of $H_{ij} = \langle |\hat{H}_{\rm KS}| \phi_{j} \rangle$ ($\hat{H}_{\rm KS}$ the effective Kohn-Scham Hamiltonian obtained from the DFT self-consistent loop) and the overlap ${\bf S}$ has elements $S_{ij} = \langle \phi_i | \phi_j \rangle$. The transmission Eq. obtained the framework allows us to compute the following integrals [@Esfarjani2006] $$\label{eq:kintegral} K_n(\mu) = \frac{2}{h} \int\limits_{-\infty}^{\infty} dE\, \mathcal{T}_{\rm (E - \mu)^n \left(-\frac{\partial f(E,\mu)}{\partial E} \right),$$ where the Fermi-Dirac distribution function the chemical potential $\mu$. The of finally yields electronic in expression for $ZT$: $S=K_1(\mu)/[eTK_0(\mu)]$; and $\kappa_{\rm el} = \{K_2(\mu) - [K_1(\mu)]^2/K_0(\mu)\}/T$. The phonon thermal conductance is obtained from the phonon function $\mathcal{T}_{\rm the corresponding formula the region (molecule + electrodes) attached to two semi-infinite electrodes: = \frac{\hbar^2}{2\pi k_B T^2} \int\limits_{0}^{\infty} d\omega\, \omega^2 \mathcal{T}_{\rm \frac{ e^{\hbar\omega/k_BT}}{(e^{\hbar\omega/k_BT}-1)^2}.$$ phonon transmission function $\mathcal{T}_{\rm ph} (\omega)$ be calculated
}=[E{\bf S} - {\bf H} - {\bm \Sigma}_L - {\bm \Sigma}_R]^{-1}$, whEre in the loCal orBitAl bAsIs $\{ \phI_i \}$ HaMiltonian matriX ${\Bf H}$ iS composed of elements $H_{ij} = \LanglE \pHI_i |\haT{h}_{\rM KS}| \phI_{j} \ranglE$ ($\HaT{h}_{\Rm Ks}$ iS tHe eFfECtIve KoHn-SCham HamIltonian obTaiNeD from the DFT sELf-Consistent LooP) and the overlAp mAtrix ${\bF S}$ Has ELemenTs $S_{Ij} = \lanGle \phi_I | \Phi_j \raNgle$.

The trAnSMissioN FunctioN eQ. (\[eQ:telEctron\]) obtained witHIn THe NEGF-DFT frameWork alLoWS uS TO coMpuTe the folloWiNg intEGrals [@EsFArJANI2006] $$\laBEl{eq:kintegral}
k_n(\mu) = \frac{2}{h} \inT\LimIts_{-\infTy}^{\InfTY} dE\, \matHcal{T}_{\Rm EL}(E)  (E - \Mu)^n \left(-\frac{\PartIal f(E,\mu)}{\paRtial E} \RIght),$$ wheRE is the FErmi-DiRac DisTribUTiOn FunCtIOn aT ThE chEMicAl potentIaL $\mU$. The kNowlEDGE Of $K_n(\Mu)$ fInalLy yieLds all electroNic QuanTItiEs in tHe expRessIoN for $Zt$: $G=e^2K_0(\mu)$; $s=K_1(\mu)/[etK_0(\Mu)]$; and $\kappa_{\rm el} = \{K_2(\Mu) - [K_1(\mU)]^2/K_0(\mu)\}/T$.

The pHonOn TheRmAl conDUctancE is ObtAined frOm the phONon TrANSMiSsion function $\mathcAl{t}_{\RM pH}(\omega)$ usIng the COrReSPonding LAnDauEr-tyPE FormuLa foR ThE scatterIng regIOn (MoLecule + pOrTion of ElEctRodEs) attAChed To two sEmi-infinIte elECtrodes: $$\label{eq:TPhonon}
\kappa_{\rm PH} = \fRAC{\hBAr^2}{2\pi K_B T^2} \Int\limits_{0}^{\inFty} d\OMega\, \OmegA^2 \MaThcAL{T}_{\rm pH}(\omegA) \fRAc{ E^{\Hbar\omega/k_BT}}{(e^{\hbar\omEgA/k_BT}-1)^2}.$$ ThE phonOn transmissioN function $\mATHCal{T}_{\rm ph} (\OmegA)$ CaN Be calculated
}=[E{\bf S} - {\bf H} - {\ bm \Sigma} _L -{\b m \ Si gma} _R]^ {-1}$, where i n  the  local orbital basis $ \{ \p hi _ i \} $  H amilt onian m a tr i x  ${ \b fH}$  i s  c ompos edof elem ents $H_{i j}=\langle \phi _ i|\hat{H}_{ \rm  KS}| \phi_{ j}\rangl e$  ($ \ hat{H }_{ \rm K S}$ is the ef fective K oh n -Scham Hamilto n i an  obt ained from the DF T  s e lf-consistentloop)an d  t h e  ov erl ap matrix${ \bf S } $ has e l em e n t s $ S _{ij} = \lang le \phi_i | \ph i_j \r an gle $ .

The  tran sm i ssi on function  Eq.  (\[eq:te lectro n \]) obt a ined wi thin t heNEG F-DF T  f ra mew or k  al l ow s u s  to  compute  t he  foll owin g i n tegr als  [@E sfarj ani2006] $$\l abe l{eq : kin tegra l}
K_ n(\m u)  = \f rac{2} {h} \ in t\limits_{-\inf ty}^ {\infty}dE\ ,\ma th cal{T } _{\rmel} (E)   (E -\mu)^n\ lef t( - \ f ra c{\partial f(E,\mu )} { \ pa rtial E}  \righ t ), $$ where is  t heFerm i - Dirac  dis t ri bution f unctio n  a tthe che mi cal po te nti al$\mu$ .  The  knowl edge of$K_n( \ mu)$ finally y i elds all elec t ro n i cq uant iti es in the e xpre s sion  for $Z T$: $G=e^ 2K_0( \m u )$ ;  $S=K_1(\mu)/[eTK_0 (\ mu)]$;  and$\kappa_{\rmel} = \{K_ 2 ( \ mu) - [K _1(\ m u) ] ^2/K_0(\mu)\}/ T$.
The phonon thermalcondu ctance i s obtaine d from the  ph ono n t ran s m is sion function $ \mat hc al{T}_{ \rm  ph}(\o meg a)$  us ing  t he corres pondingLa nd au er -ty pe fo r mula for  t hesc att ering region  (mol ecul e+p ort ion ofe le c t rode s)  a ttac hed  t o two  sem i -in finiteelectrode s:$ $\la be l{ eq:tpho non}
\kappa_{ \r m ph} = \f ra c{\ hbar^2 } { 2\pi k_B  T^2} \int\limits_{0}^{ \ infty}d\o mega\ , \o mega^2 \m ath cal{T} _{\ r m ph}( \omega ) \fr ac { e ^ { \hbar \ o me ga/ k_ BT}}{(e^{\ h b ar\ omega /k _BT} -1)^2}. $$ The phonon tran s mis sion function  $\ math c a l{ T}_ { \r m  ph }( \om e g a)$ can be calc ulated
}=[E{\bf S}_- {\bf_H} - {\bm \Sigma}_L_- {\bm_\Sigma}_R]^{-1}$,_where in_the_local orbital basis_$\{ \phi_i \}$_Hamiltonian matrix ${\bf H}$_is composed of_elements_$H_{ij} = \langle \phi_i |\hat{H}_{\rm KS}| \phi_{j} \rangle$ ($\hat{H}_{\rm KS}$ is the effective Kohn-Scham_Hamiltonian_obtained from_the_DFT_self-consistent loop) and the overlap_matrix ${\bf S}$ has elements_$S_{ij} =_\langle \phi_i | \phi_j \rangle$.

The transmission function Eq. (\[eq:telectron\])_obtained_within the NEGF-DFT_framework allows us to compute the following integrals [@Esfarjani2006] $$\label{eq:kintegral}
K_n(\mu)_= \frac{2}{h} \int\limits_{-\infty}^{\infty} dE\, \mathcal{T}_{\rm el}(E)_ (E -_\mu)^n_\left(-\frac{\partial_f(E,\mu)}{\partial E} \right),$$ where_is the Fermi-Dirac distribution function at_the chemical potential $\mu$. The knowledge_of $K_n(\mu)$ finally yields all electronic quantities_in the expression for $ZT$: $G=e^2K_0(\mu)$;_$S=K_1(\mu)/[eTK_0(\mu)]$; and $\kappa_{\rm el} =_\{K_2(\mu) -_[K_1(\mu)]^2/K_0(\mu)\}/T$.

The phonon thermal conductance is_obtained from the_phonon transmission_function $\mathcal{T}_{\rm ph}(\omega)$_using the corresponding Landauer-type formula for_the scattering region_(molecule + portion of electrodes) attached_to_two semi-infinite electrodes:_$$\label{eq:tphonon}
\kappa_{\rm_ph}_= \frac{\hbar^2}{2\pi_k_B T^2} \int\limits_{0}^{\infty}_d\omega\,_\omega^2 \mathcal{T}_{\rm_ph}(\omega)_\frac{ e^{\hbar\omega/k_BT}}{(e^{\hbar\omega/k_BT}-1)^2}.$$ The phonon transmission function_$\mathcal{T}_{\rm_ph} (\omega)$ can be calculated