
(\vec T' + \vec \tau_\beta)} \langle \phi_\beta (\vec T') | {\cal H} | \phi_\alpha (\vec T) \rangle e^{i \vec k \cdot (\vec T +
\vec \tau_\alpha)}\end{aligned}$$ Assuming that the inter-site tunneling amplitude depends on the layer-projected difference coordinate, the matrix element can be expressed $$\begin{aligned}
\langle \phi_\beta (\vec T') | {\cal H} | \phi_\alpha (\vec T) \rangle = \frac{1}{(2 \pi)^2} \int \, d^2 q \, f(\vec q) \, e^{i
\vec q \cdot (\vec T' + \vec \tau'_\beta - \vec T - \vec \tau_\alpha)}\end{aligned}$$ Carrying out the lattice sums in Eqn. 21 and expressing the momenta in terms of their differences from the respective zone corners, $\vec k = \vec K + \vec q$, one obtains an expression for the interlayer tunneling amplitude in terms of sums over the reciprocal lattices of the reference ($\vec G$) and rotated ($\vec G'$) layers $$\begin{aligned}
{\cal T}_{\beta \alpha} (\vec q',\vec q) = \frac{1}{\cal A} \sum_{\vec G, \vec G'} \, f(\vec q + \vec K + \vec G) \, e^{i \vec
G' \cdot
\vec \tau_\beta} e^{-i \vec  G \cdot \tau_\alpha} \, \nonumber\\
 \,\,\,\,\,\,\,\,\,\,\,\,\,\, \delta(\vec q' - \vec q +  \Delta K  + \vec G' -
\vec G)\end{aligned}$$ where ${\cal A}$ is the area of the unit cell.

When $q \ll G$ Eqn. 23 describes two distinct types of interlayer tunneling processes: (1) Direct interlayer terms conserve the [*crystal momentum*]{} $\vec k$ and occur when $\Delta K = |\vec K(\theta) - \vec K + \vec G' - \vec G|$. (Note that this occurs for $\vec G = \vec G' = 0$ and for all boosts by the reciprocal lattice vectors $\vec {\cal G} = \vec G' - \vec G$ that symmetr
(\vec T' + \vec \tau_\beta) } \langle \phi_\beta (\vec T') | { \cal H } | \phi_\alpha (\vec T) \rangle e^{i \vec k \cdot (\vec T + 
 \vec \tau_\alpha)}\end{aligned}$$ Assuming that the inter - site tunneling amplitude depends on the layer - projected remainder coordinate, the matrix chemical element can be expressed $ $ \begin{aligned } 
 \langle \phi_\beta (\vec T') | { \cal H } | \phi_\alpha (\vec T) \rangle = \frac{1}{(2 \pi)^2 } \int \, d^2 q \, f(\vec q) \, e^{i 
 \vec q \cdot (\vec T' + \vec \tau'_\beta - \vec T - \vec \tau_\alpha)}\end{aligned}$$ Carrying out the wicket sums in Eqn. 21 and expressing the momenta in condition of their differences from the respective zone corners, $ \vec k = \vec K + \vec q$, one obtain an expression for the interlayer tunneling amplitude in term of sums over the reciprocal lattices of the reference ($ \vec G$) and rotated ($ \vec G'$) layers $ $ \begin{aligned } 
 { \cal T}_{\beta \alpha } (\vec q',\vec q) = \frac{1}{\cal A } \sum_{\vec G, \vec G' } \, f(\vec q + \vec K + \vec G) \, e^{i \vec 
 G' \cdot 
 \vec \tau_\beta } e^{-i \vec   G \cdot \tau_\alpha } \, \nonumber\\ 
  \,\,\,\,\,\,\,\,\,\,\,\,\,\, \delta(\vec q' - \vec q +   \Delta K   + \vec gram' - 
 \vec G)\end{aligned}$$ where $ { \cal A}$ is the area of the unit cell. 

 When $ q \ll G$ Eqn. 23 identify two distinct types of interlayer tunneling processes: (1) Direct interlayer footing conserve the [ * crystal momentum * ] { } $ \vec k$ and occur when $ \Delta K = |\vec K(\theta) - \vec K + \vec G' - \vec G|$. (Note that this occurs for $ \vec G = \vec G' = 0 $ and for all boost by the reciprocal lattice vectors $ \vec { \cal G } = \vec G' - \vec G$ that symmetr

(\vec T' + \vec \tau_\beta)} \langle \pmi_\beta (\vec T') | {\cal H} | \phm_\alpha (\bec T) \ravgle e^{i \vec k \cdot (\vec T +
\vec \tqu_\alpya)}\end{aligned}$$ Assuming ghat the pnter-site tunieling amplitude depends on ths layzr-'rojected diffetence coorditate, the matrif dlzment can be expressed $$\begin{aligned}
\lwngle \pni_\heta (\vec T') | {\caj H} | \[hi_\amihc (\vec T) \rangle = \frac{1}{(2 \pi)^2} \int \, d^2 s \, f(\vec q) \, e^{i
\vec q \ccot (\vec T' + \vec \tau'_\beta - \vef T - \vec \tau_\alpha)}\end{allgned}$$ Carryung jyt the lattize sums in Eqn. 21 and exkressing the momenta in terms of theix differencgs fgmm the respxctive zone corners, $\vec k = \vec K + \vec q$, one obtsinv ab expression for the mnterlayer tunneling amplituda nn terms of sums over tye rewiprmcal oatgicts pf the rffecence ($\vec G$) and rotatee ($\vec G'$) layers $$\begim{ajpbned}
{\cal T}_{\betz \alphw} (\dec q',\vec q) = \frac{1}{\cal A} \sum_{\vec G, \vec G'} \, f(\nec s + \vec K + \vec G) \, e^{i \vec
T' \cdot
\vec \tau_\beta} e^{-i \gec  G \cdjt \tau_\alpha} \, \nonumber\\
 \,\,\,\,\,\,\,\,\,\,\,\,\,\, \delta(\vec q' - \vec q +  \Delta N  + \vxc G' -
\rcg G)\eve{apigned}$$ where ${\cal A}$ is the area of the unit cejm.

Wnek $q \ll G$ Eqn. 23 defcribes two dlsyynct types of interlaysr tunneling procedses: (1) Dyrect interlaytr tetms conserve the [*crystal monentum*]{} $\vec k$ qnd occur when $\Delca K = |\vec K(\tketa) - \fec K + \vec G' - \vec G|$. (Note thac this occurs for $\vec G = \vsz G' = 0$ and for alu bpovts by the reciprocal lattyce vectocs $\vee {\cal G} = \vec G' - \ves G$ that sjmmetv
(\vec T' + \vec \tau_\beta)} \langle \phi_\beta | H} | (\vec T) \rangle T \vec \tau_\alpha)}\end{aligned}$$ Assuming the inter-site tunneling depends on the layer-projected difference coordinate, matrix element can be expressed $$\begin{aligned} \langle \phi_\beta (\vec T') | {\cal H} \phi_\alpha (\vec T) \rangle = \frac{1}{(2 \pi)^2} \int \, d^2 q \, f(\vec \, \vec \cdot T' + \vec \tau'_\beta - \vec T - \vec \tau_\alpha)}\end{aligned}$$ Carrying out the lattice sums in 21 and expressing the momenta in terms of differences from the respective corners, $\vec k = \vec + q$, one an for interlayer tunneling amplitude terms of sums over the reciprocal lattices of the reference ($\vec G$) and rotated ($\vec G'$) layers {\cal T}_{\beta q',\vec q) \frac{1}{\cal \sum_{\vec \vec G'} \, + \vec K + \vec G) G' \cdot \vec \tau_\beta} e^{-i \vec G \cdot \, \nonumber\\ \delta(\vec q' - \vec q + K + \vec G' - \vec G)\end{aligned}$$ where A}$ is the area of the unit cell. When $q \ll G$ Eqn. 23 describes types of interlayer tunneling (1) Direct interlayer conserve [*crystal $\vec and occur $\Delta K = |\vec K(\theta) - \vec K + \vec G' \vec G|$. (Note that this occurs for $\vec G = = and for all by the reciprocal lattice $\vec G} = \vec G' G$ symmetr

(\vec T' + \vec \tau_\beta)} \langle \phi_\beTa (\vec T') | {\cal H} | \Phi_\alPha (\Vec t) \rAnglE e^{i \vEc k \cdot (\vec T +
\vec \TAu_\alPha)}\end{aligned}$$ Assuming tHat thE iNTer-sITe TunneLing ampLItUDE dePeNdS on ThE LaYer-prOjeCted difFerence cooRdiNaTe, the matrix eLEmEnt can be exPreSsed $$\begin{aliGneD}
\langlE \pHi_\bETa (\vec t') | {\caL H} | \phi_\Alpha (\vEC T) \rangLe = \frac{1}{(2 \pi)^2} \iNt \, D^2 Q \, f(\vec q) \, E^{I
\vec q \cdOT (\VeC T' + \veC \tau'_\beta - \vec T - \vec \taU_\AlPHa)}\end{aligned}$$ CaRrying OuT ThE LAttIce Sums in Eqn. 21 aNd ExpreSSing the MOmENTA in TErms of their diFferences frOM thE respeCtIve ZOne corNers, $\vEc K = \Vec k + \vec q$, one obtAins An expressIon for THe interLAyer tunNeling AmpLitUde iN TeRmS of SuMS ovER tHe rECipRocal latTiCeS of thE refERENCe ($\veC G$) aNd roTated ($\Vec G'$) layers $$\begIn{aLignED}
{\caL T}_{\betA \alphA} (\vec Q',\vEc q) = \frAc{1}{\cal A} \Sum_{\veC G, \Vec G'} \, f(\vec q + \vec K + \veC G) \, e^{i \Vec
G' \cdot
\vEc \tAu_\BetA} e^{-I \vec  G \CDot \tau_\AlpHa} \, \nOnumber\\
 \,\,\,\,\,\,\,\,\,\,\,\,\,\, \Delta(\veC Q' - \veC q +  \dELTa k  + \vec G' -
\vec G)\end{aligneD}$$ wHERe ${\Cal A}$ is thE area oF ThE uNIt cell.

WhEn $Q \ll g$ Eqn. 23 DEScribEs twO DiStinct tyPes of iNTeRlAyer tunNeLing prOcEssEs: (1) DIrect INterLayer tErms consErve tHE [*crystal momentUM*]{} $\vec k$ and occur WHeN $\dElTA K = |\veC K(\tHeta) - \vec K + \vec g' - \vec g|$. (note That THiS ocCUrs foR $\vec G = \VeC g' = 0$ aND for all boosts by the rEcIprocaL lattIce vectors $\vec {\Cal G} = \vec G' - \veC g$ THat symmeTr

(\vec T' + \vec \tau_\bet a)} \langl e \ph i_\ bet a(\ve c T' ) | {\cal H} | \phi _\alpha (\vec T) \rang le e^ {i \vec k\cdot  (\vecT  + \ vec  \ ta u_\ al p ha )}\en d{a ligned} $$ Assumin g t ha t the inter- s it e tunnelin g a mplitude dep end s on t he  la y er-pr oje cteddiffer e nce co ordinate,  t h e matr i x eleme n t  c an b e expressed $$\be g in { aligned}
\lang le \ph i_ \ be t a  (\ vec  T') | {\c al  H} | \phi_\a l ph a ( \ve c  T) \rangle =  \frac{1}{( 2  \p i)^2}\i nt\ , d^2q \,f( \ vec  q) \, e^{i 
\ve c q \cdot  (\vec T' + \v e c \tau' _\beta  -\ve c T-  \ ve c \ ta u _\a l ph a)} \ end {aligned }$ $Carry ingo u t thelat tice  sums  in Eqn. 21 a ndexpr e ssi ng th e mom enta  i n ter ms oftheir  d ifferences from  the  respecti vezo neco rners ,  $\vec  k= \ vec K +  \vec q $ , o ne o b ta ins an expressionfo r th e interl ayer t u nn el i ng ampli tu dein t e r ms of  sum s  o ver therecipr o ca llattice sof the  r efe ren ce ($ \ vecG$) an d rotate d ($\ v ec G'$) layers $$\begin{alig n ed } {\ c al T }_{ \beta \alph a} ( \ vecq',\ v ec  q) = \fr ac{1} {\ c al A} \sum_{\vec G, \v ec  G'} \ , f(\ vec q + \vecK + \vec G ) \ , e^{i \ vecG '\ cdot
\vec \tau _\bet a} e^{-i \ v ec  G \c dot \ tau_\alp ha} \, \n o n umber\\ \, \,\ ,\, \,\ , \ ,\ ,\,\,\,\,\,\, \ delt a( \vec q'  -\vec q+ \De lta  K +  \vec G'-
\vec G )\ en d{ al ign ed}$$ where ${ \c alA} $ i s the area o f the  uni tce l l.
When $ q  \ l l  G$Eq n.  23des cr ibestwod ist inct ty pes of in ter l ayer  t un nelingprocesses: (1 )Direct int er lay er ter m s  conserv e the [*crystal momentu m *]{} $\ vec  k$ a nd o ccur when  $\ DeltaK = |\vecK(\the ta) -  \ vec K  + \v e c  G ' -  \ vec G|$. ( N o tethatth is o ccurs f or $\vec G = \vecG ' =  0$ and for a llboos t s  b y t h er eci pr o cal l attice vectors$\vec {\ca lG }= \vec G'-  \v ec  G$ tha t symme tr

(\vec T'_+ \vec_\tau_\beta)} \langle \phi_\beta (\vec_T') |_{\cal_H} |_\phi_\alpha_(\vec T) \rangle_e^{i \vec k_\cdot (\vec T +
\vec_\tau_\alpha)}\end{aligned}$$ Assuming that_the_inter-site tunneling amplitude depends on the layer-projected difference coordinate, the matrix element can be_expressed_$$\begin{aligned}
\langle \phi_\beta_(\vec_T')_| {\cal H} | \phi_\alpha_(\vec T) \rangle = \frac{1}{(2_\pi)^2} \int_\, d^2 q \, f(\vec q) \, e^{i
\vec_q_\cdot (\vec T'_+ \vec \tau'_\beta - \vec T - \vec \tau_\alpha)}\end{aligned}$$_Carrying out the lattice sums in_Eqn. 21 and_expressing_the_momenta in terms of_their differences from the respective zone_corners, $\vec k = \vec K_+ \vec q$, one obtains an expression_for the interlayer tunneling amplitude in_terms of sums over the_reciprocal lattices_of the reference ($\vec G$)_and rotated ($\vec_G'$) layers_$$\begin{aligned}
{\cal T}_{\beta \alpha}_(\vec q',\vec q) = \frac{1}{\cal A}_\sum_{\vec G, \vec_G'} \, f(\vec q + \vec_K_+ \vec G)_\,_e^{i_\vec
G' \cdot
\vec_\tau_\beta} e^{-i \vec__G \cdot_\tau_\alpha}_\, \nonumber\\
 \,\,\,\,\,\,\,\,\,\,\,\,\,\, \delta(\vec q' -_\vec_q +  \Delta K  +_\vec G' -
\vec G)\end{aligned}$$_where_${\cal A}$ is the_area of the unit cell.

When_$q \ll G$ Eqn. 23 describes_two distinct_types of_interlayer tunneling processes: (1) Direct interlayer terms conserve the [*crystal momentum*]{}_$\vec k$ and occur when $\Delta_K = |\vec K(\theta)_- \vec_K_+ \vec G'_-_\vec G|$._(Note that this occurs for $\vec G_= \vec_G' = 0$ and for all_boosts by the reciprocal_lattice_vectors $\vec {\cal G} = \vec_G' - \vec G$ that symmetr