
\omega\right\vert \ll\left\vert \delta_{i}\right\vert $ such that $\delta
_{1}\approx\delta_{2}$. In the tight-binding limit, the external potential corresponds to the Hamiltonian $$H_{\text{e}}\left(  t\right)  =\hbar\Omega\sum_{l=0,1\text{; }q=\text{L,R}}\cos\left(  \varphi_{l,q}\mathbf{+}\omega t\right)  b_{l,q}^{\dag}b_{l,q},
\label{eq:He}$$ where $\varphi_{l,\text{L/R}}=\mp\frac{k_{x}\lambda_{\text{s}}}{4}+\frac{lk_{y}\lambda_{\text{s}}}{2}$. The wavevectors $k_{x}$ and $k_{y}$ can be represented as $k_{x}=k_{1}\cos\phi_{1}-k_{2}\cos\phi_{2}$ and $k_{y}=k_{1}\sin\phi_{1}-k_{2}\sin\phi_{2}$. To eliminate the time-dependent Hamiltonian $H_{\text{e}}\left(  t\right)  $, we perform a unitary transformation $U=\!\exp[-i\frac{\Omega}{\omega}\sum_{l=0,1\text{; }q=\text{L,R}}\sin\left(
\varphi_{l,q}\mathbf{+}\omega t\right)  b_{l,q}^{\dag}b_{l,q}]$, which yields the following effective Hamiltonian$$\begin{aligned}
H_{\text{eff}}=  &  \sum_{l}\frac{\Delta}{2}\left(  b_{l,\text{R}}^{\dag
}b_{l,\text{R}}-b_{l,\text{L}}^{\dag}b_{l,\text{L}}\right) \nonumber\\
&  -\sum_{\substack{l\neq-1,0,1\\q=\text{L,R}}}J_{y}b_{l+1,q}^{\dag}b_{l,q}-\sum_{l\neq0,1}J_{x}b_{l,\text{L}}^{\dag}b_{l,\text{R}}+\text{H.c.}\nonumber\\
&  -\sum
\omega\right\vert \ll\left\vert \delta_{i}\right\vert $ such that $ \delta 
 _ { 1}\approx\delta_{2}$. In the tight - binding limit, the external potential corresponds to the Hamiltonian $ $ H_{\text{e}}\left (  t\right)   = \hbar\Omega\sum_{l=0,1\text {; } q=\text{L, R}}\cos\left (  \varphi_{l, q}\mathbf{+}\omega t\right)   b_{l, q}^{\dag}b_{l, q }, 
 \label{eq: He}$$ where $ \varphi_{l,\text{L / R}}=\mp\frac{k_{x}\lambda_{\text{s}}}{4}+\frac{lk_{y}\lambda_{\text{s}}}{2}$. The wavevectors $ k_{x}$ and $ k_{y}$ can be represent as $ k_{x}=k_{1}\cos\phi_{1}-k_{2}\cos\phi_{2}$ and $ k_{y}=k_{1}\sin\phi_{1}-k_{2}\sin\phi_{2}$. To rule out the time - dependent Hamiltonian $ H_{\text{e}}\left (  t\right)   $, we perform a one transformation $ U=\!\exp[-i\frac{\Omega}{\omega}\sum_{l=0,1\text {; } q=\text{L, R}}\sin\left (
 \varphi_{l, q}\mathbf{+}\omega t\right)   b_{l, q}^{\dag}b_{l, q}]$, which yields the following effective Hamiltonian$$\begin{aligned } 
 H_{\text{eff}}=   &   \sum_{l}\frac{\Delta}{2}\left (  b_{l,\text{R}}^{\dag 
 } b_{l,\text{R}}-b_{l,\text{L}}^{\dag}b_{l,\text{L}}\right) \nonumber\\ 
 &   -\sum_{\substack{l\neq-1,0,1\\q=\text{L, R}}}J_{y}b_{l+1,q}^{\dag}b_{l, q}-\sum_{l\neq0,1}J_{x}b_{l,\text{L}}^{\dag}b_{l,\text{R}}+\text{H.c.}\nonumber\\ 
 &   -\sum

\omeha\right\vert \ll\left\vert \dtlta_{i}\right\vert $ socy that $\delta
_{1}\zpprox\deuta_{2}$. In the tight-binding limiv, thw exttgnal potential corresoonds to nhe Hamilronien $$H_{\text{e}}\left(  t\cjght)  =\hncr\Omefw\sum_{n=0,1\vext{; }q=\text{L,R}}\cos\keft(  \varpvi_{l,q}\mathbf{+}\omegd g\rnght)  b_{l,q}^{\dag}b_{l,q},
\label{eq:He}$$ where $\varphy_{l,\text{L/T}}=\mo\frac{k_{x}\lambda_{\tgxt{s}}}{4}+\fgas{lk_{y}\mambda_{\text{s}}}{2}$. The wavevectors $k_{x}$ and $k_{y}$ can be representrd as $k_{x}=k_{1}\cos\phi_{1}-k_{2}\cos\phi_{2}$ and $k_{y}=k_{1}\din\phi_{1}-k_{2}\sin\phi_{2}$. To epiminate thg tiiw-dependent Hxmiltonian $H_{\text{e}}\left(  t\right)  $, we perform a unitary gransyormation $U=\!\grp[-i\fgdc{\Omega}{\omege}\sum_{l=0,1\tvxt{; }q=\text{L,R}}\sik\keft(
\vasphi_{l,q}\msthbf{+}\omega t\rinht)  u_{l,q}^{\dqg}b_{l,q}]$, which yields thx following effectivg Hamiltonhau$$\begin{aligned}
H_{\text{eff}}=  &  \sum_{n}\fraw{\Delgq}{2}\lewt(  b_{k,\tsxt{R}}^{\dah
}b_{l,\vext{R}}-b_{l,\text{M}}^{\dag}b_{l,\text{L}}\eight) \nonumber\\
&  -\sum_{\xufwtack{l\neq-1,0,1\\q=\texf{L,R}}}J_{y}b_{j+1,q}^{\qag}b_{l,q}-\sum_{l\neq0,1}J_{x}b_{l,\text{L}}^{\dag}b_{l,\text{R}}+\text{H.c.}\nmnujber\\
&  -\sum
\omega\right\vert \ll\left\vert \delta_{i}\right\vert $ such that $\delta the limit, the potential corresponds to }q=\text{L,R}}\cos\left( t\right) b_{l,q}^{\dag}b_{l,q}, \label{eq:He}$$ $\varphi_{l,\text{L/R}}=\mp\frac{k_{x}\lambda_{\text{s}}}{4}+\frac{lk_{y}\lambda_{\text{s}}}{2}$. The wavevectors and $k_{y}$ can be represented as and $k_{y}=k_{1}\sin\phi_{1}-k_{2}\sin\phi_{2}$. To eliminate the time-dependent Hamiltonian $H_{\text{e}}\left( t\right) $, we perform a transformation $U=\!\exp[-i\frac{\Omega}{\omega}\sum_{l=0,1\text{; }q=\text{L,R}}\sin\left( \varphi_{l,q}\mathbf{+}\omega t\right) b_{l,q}^{\dag}b_{l,q}]$, which yields the following effective Hamiltonian$$\begin{aligned} H_{\text{eff}}= \sum_{l}\frac{\Delta}{2}\left( }b_{l,\text{R}}-b_{l,\text{L}}^{\dag}b_{l,\text{L}}\right) & & -\sum

\omega\right\vert \ll\left\vert \dElta_{i}\right\Vert $ sUch ThaT $\dElta
_{1}\ApprOx\delta_{2}$. In the tiGHt-biNding limit, the external pOtentIaL CorrESpOnds tO the HamILtONIan $$h_{\tExT{e}}\lEfT(  T\rIght)  =\hBar\omega\suM_{l=0,1\text{; }q=\texT{L,R}}\CoS\left(  \varphi_{l,Q}\MaThbf{+}\omega t\RigHt)  b_{l,q}^{\dag}b_{l,q},
\lAbeL{eq:He}$$ wHeRe $\vARphi_{l,\TexT{L/R}}=\mp\Frac{k_{x}\LAmbda_{\tExt{s}}}{4}+\frac{lK_{y}\LAmbda_{\tEXt{s}}}{2}$. The wAVEvEctoRs $k_{x}$ and $k_{y}$ can be repREsENted as $k_{x}=k_{1}\cos\phI_{1}-k_{2}\cos\pHi_{2}$ ANd $K_{Y}=K_{1}\siN\phI_{1}-k_{2}\sin\phi_{2}$. To ElIminaTE the timE-DePENDenT hamiltonian $H_{\tExt{e}}\left(  t\riGHt)  $, wE perfoRm A unITary trAnsfoRmATioN $U=\!\exp[-i\frac{\OMega}{\Omega}\sum_{l=0,1\Text{; }q=\tEXt{L,R}}\sin\LEft(
\varpHi_{l,q}\maThbF{+}\omEga t\RIgHt)  B_{l,q}^{\DaG}B_{l,q}]$, WHiCh yIEldS the follOwInG effeCtivE hAMIltoNiaN$$\begIn{aliGned}
H_{\text{eff}}=  &  \sUm_{l}\Frac{\dEltA}{2}\left(  B_{l,\texT{R}}^{\daG
}b_{L,\text{r}}-b_{l,\texT{L}}^{\dag}B_{l,\Text{L}}\right) \nonumBer\\
&  -\sUm_{\substacK{l\nEq-1,0,1\\Q=\teXt{l,R}}}J_{y}b_{L+1,Q}^{\dag}b_{l,Q}-\suM_{l\nEq0,1}J_{x}b_{l,\tExt{L}}^{\dag}B_{L,\teXt{r}}+\TEXt{h.c.}\nonumber\\
&  -\sum

\omega\right\vert \ll\lef t\vert \de lta_{ i}\ rig ht \ver t $such that $\de l ta
_ {1}\approx\delta_{2}$.  In t he tigh t -b indin g limit ,  t h e  ex te rn alpo t en tialcor respond s to the H ami lt onian $$H_{\ t ex t{e}}\left ( t\right)  =\ hba r\Omeg a\ sum _ {l=0, 1\t ext{;  }q=\t e xt{L,R }}\cos\le ft (   \var p hi_{l,q } \ ma thbf {+}\omega t\right )  b _{l,q}^{\dag}b _{l,q} ,\ la b e l{e q:H e}$$ where  $ \varp h i_{l,\t e xt { L / R}} = \mp\frac{k_{x }\lambda_{\ t ext {s}}}{ 4} +\f r ac{lk_ {y}\l am b da_ {\text{s}}} {2}$ . The wav evecto r s $k_{x } $ and $ k_{y}$  ca n b e re p re se nte da s $ k _{ x}= k _{1 }\cos\ph i_ {1 }-k_{ 2}\c o s \ p hi_{ 2}$  and  $k_{ y}=k_{1}\sin\ phi _{1} - k_{ 2}\si n\phi _{2} $.  To e limina te th etime-dependentHami ltonian $ H_{ \t ext {e }}\le f t(  t\ rig ht)   $, we  perfor m  aun i t a ry  transformation $U =\ ! \ ex p[-i\fra c{\Ome g a} {\ o mega}\su m_ {l= 0,1\ t e xt{;}q=\ t ex t{L,R}}\ sin\le f t( 
\ varphi_ {l ,q}\ma th bf{ +}\ omega t\ri ght) b_{l,q}^ {\dag } b_{l,q}]$, whi c h yields thef ol l o wi n g ef fec tive Hamilt onia n $$\b egin { al ign e d}
H_ {\tex t{ e ff } }=  &  \sum_{l}\fra c{ \Delta }{2}\ left(  b_{l,\ text{R}}^{ \ d a g
}b_{l, \tex t {R } }-b_{l,\text{L }}^{\ dag}b_{l,\ t ext{L}}\ right ) \nonum ber\\
& - \ sum_{\su bst ack {l\ neq - 1 ,0 ,1\\q=\text{L , R }}}J _{ y}b_{l+ 1,q }^{\dag }b_ {l, q}- \su m_ {l\neq0,1 }J_{x}b_ {l ,\ te xt {L} }^{\d a g}b_{l,\ te xt{ R} }+\ text{ H .c.}\n onumb er\\ 
&  - \su m

\omega\right\vert \ll\left\vert_\delta_{i}\right\vert $_such that $\delta
_{1}\approx\delta_{2}$. In_the tight-binding_limit,_the external_potential_corresponds to the_Hamiltonian $$H_{\text{e}}\left( _t\right)  =\hbar\Omega\sum_{l=0,1\text{; }q=\text{L,R}}\cos\left(_ \varphi_{l,q}\mathbf{+}\omega t\right)__b_{l,q}^{\dag}b_{l,q},
\label{eq:He}$$ where $\varphi_{l,\text{L/R}}=\mp\frac{k_{x}\lambda_{\text{s}}}{4}+\frac{lk_{y}\lambda_{\text{s}}}{2}$. The wavevectors $k_{x}$ and $k_{y}$ can be represented as $k_{x}=k_{1}\cos\phi_{1}-k_{2}\cos\phi_{2}$ and_$k_{y}=k_{1}\sin\phi_{1}-k_{2}\sin\phi_{2}$._To eliminate_the_time-dependent_Hamiltonian $H_{\text{e}}\left(  t\right) _$, we perform a unitary_transformation $U=\!\exp[-i\frac{\Omega}{\omega}\sum_{l=0,1\text{;_}q=\text{L,R}}\sin\left(
\varphi_{l,q}\mathbf{+}\omega t\right)  b_{l,q}^{\dag}b_{l,q}]$, which yields the following_effective_Hamiltonian$$\begin{aligned}
H_{\text{eff}}=  &_ \sum_{l}\frac{\Delta}{2}\left(  b_{l,\text{R}}^{\dag
}b_{l,\text{R}}-b_{l,\text{L}}^{\dag}b_{l,\text{L}}\right) \nonumber\\
&  -\sum_{\substack{l\neq-1,0,1\\q=\text{L,R}}}J_{y}b_{l+1,q}^{\dag}b_{l,q}-\sum_{l\neq0,1}J_{x}b_{l,\text{L}}^{\dag}b_{l,\text{R}}+\text{H.c.}\nonumber\\
&  -\sum