“…, where d is the inter-layer distance, r ij is the distance between sites i and j, λ is the parameter that controls the decay of the interlayer coupling, and t ⊥,0 is the largest possible inter-layer coupling realized at r ij = d. In the following we take t ⊥,0 = 0.36t, λ = 10/a, and d = a, where a is the lattice constant of the triangular lattice. From the computational point of view, we will use the twist scaling relation for computational convenience [65,66], and we compute the valley expectation V z = ±1 by means of the valley operator [49,[67][68][69][70].…”