In this subsection, we consider the model proposed in Ref. [1]}. Note that Ref. [1]} does not provide a general formula for the intra-layer hoppings, but lists its magnitude for a set of discrete hopping distance. Here, we fit the intra-layer hopping with the following formula:
\(t_{intra}(r) = t_0 e^{- \alpha _0 \bar{r}^2} \cos (\beta _0 \bar{r}) + t_1 \bar{r}^2 e^{- \alpha _1 (\bar{r} - r_1)^2} \) 
