
As shown in the Fig. REF  and REF , the numerically obtained lattice relaxation is dominated by the solenoidal part, \(\nabla \times \left( \hat{z} \varepsilon ^U(x) \right)\) , for the lattice relaxation models of Ref. [1]} and [2]}, with their choice for parameters listed in Table REF . The resulting displacement field  [3]} is in qualitative agreement with the Bragg interferometry imaging of the strain fields in twisted bilayer graphene [4]}. In the following calculations, we will neglect the small irrotational part, and include only the solenoidal part.
<FIGURE>