Moiré magnetism in a bilayer Ising model

Moiré patterns in magnetic bilayers generate spatially modulated interlayer exchange interactions that can give rise to nonuniform magnetic textures. We study a minimal classical bilayer Ising model with a moiré-modulated interlayer coupling, generated either by relative twist or differential strain between the layers. Using large-scale classical Monte Carlo simulations, we show that the ordering transition remains in the conventional two-dimensional Ising universality class, even when the low-temperature state is domain-textured. At low temperatures, we find a smooth crossover between a uniform ferromagnet and domain-textured state, in which the spins locally follow the sign of the interlayer exchange. We demonstrate that there is no breaking of layer symmetry for twisted bilayers. The location of the crossover is determined by a simple geometric energy balance between bulk interlayer exchange and intralayer domain-wall costs. Our results provide a minimal framework for understanding how moiré-modulated magnetic textures can emerge from geometric energetics without requiring a thermodynamic phase transition.

💡 Research Summary

The authors investigate a minimal classical Ising model for a magnetic bilayer in which a moiré pattern—generated either by a relative twist between the layers or by differential strain—produces a spatially periodic modulation of the inter‑layer exchange coupling. The Hamiltonian consists of a ferromagnetic nearest‑neighbour interaction within each layer (J = 1) and a modulated inter‑layer term J′ Φ(u)σ₁σ₂, where Φ(u)=Φ₀+∑ₐcos(bₐ·u) captures the moiré‑induced alternation between ferromagnetic and antiferromagnetic exchange. The control parameters are the twist angle ϕ (or equivalently the strain ratio a) that sets the moiré unit‑cell size L_M, and the overall inter‑layer coupling strength J′. By varying Φ₀ the authors can bias the system toward overall ferromagnetism or keep it unbiased (Φ₀≈0).

Large‑scale Monte Carlo simulations are performed on square lattices with periodic boundary conditions. Both Metropolis single‑spin updates and Swendsen‑Wang/Wolff cluster moves are employed: cluster updates efficiently equilibrate the high‑temperature paramagnetic–ordered transition, while single‑spin flips are essential for low‑temperature relaxation of domain walls. The Binder cumulant U₂=3/2(1−⟨m⁴⟩/3⟨m²⟩²) is used to locate the critical temperature Tc for several system sizes expressed in terms of the number of moiré cells N_M = L/L_M. All data collapse onto a single crossing point, indicating a continuous transition from the paramagnet directly into an ordered phase. Finite‑size scaling of the slope dU₂/dT yields an effective exponent 1/ν≈0.97±0.02, in excellent agreement with the exact 2D Ising value ν=1. Thus, despite the presence of a long‑wavelength modulation and competing ferro‑ and antiferromagnetic inter‑layer bonds, the thermal transition remains in the conventional 2D Ising universality class.

Below Tc the ordered phase can assume two distinct textures. When Φ₀ is sufficiently positive or J′ is small, the bulk ferromagnetic exchange dominates and the system adopts a uniform ferromagnetic state with all spins aligned. When Φ₀≈0 and J′ is comparable to the intra‑layer coupling, the spatially varying Φ(u) forces spins in each moiré cell to follow the local sign of the inter‑layer exchange, producing a domain‑textured state. In this state, domain walls appear at the boundaries of the moiré cells, typically in only one of the two layers (especially for the strained case where the layers have different spin densities).

A key question is whether the transition between the uniform ferromagnet and the domain‑textured state constitutes a genuine thermodynamic phase transition. The authors examine the layer‑exchange symmetry σ_{i,1}↔σ_{i,2}. In the twisted bilayer this symmetry is exact, while in the strained bilayer it is explicitly broken by the unequal lattice constants. They define a layer‑polarization order parameter ⟨P²⟩ based on the static structure factor at the moiré reciprocal vectors. For twisted bilayers, ⟨P²⟩ scales as 1/N_M and vanishes in the thermodynamic limit, showing that no spontaneous layer symmetry breaking occurs. Consequently, the crossover between the two low‑temperature textures is not a phase transition but a smooth crossover governed by energetic considerations.

The authors construct a simple geometric energy balance: the cost of creating a domain wall per moiré cell is roughly E_wall≈2J(1−Φ₀), while the gain from aligning with the modulated inter‑layer exchange is E_inter≈J′Φ₀L_M². Equating these yields a linear relation J′∝(a−1)/a (or equivalently J′∝tan ϕ) that delineates the crossover line in the (J′,a) or (J′,ϕ) plane. This prediction matches the numerical crossover boundary shown in Fig. 4, confirming that the location of the crossover is set by a straightforward competition between bulk inter‑layer exchange and intralayer domain‑wall energy.

In summary, the study demonstrates that moiré‑induced spatial modulation of inter‑layer exchange does not alter the universality class of the magnetic ordering transition in a bilayer Ising system. The emergence of domain textures at low temperature is a consequence of geometric energy balance rather than a separate symmetry‑breaking transition. These findings provide a minimal, analytically tractable framework for understanding moiré magnetism in real two‑dimensional van‑der‑Waals magnets such as CrI₃ or Cr₂Ge₂Te₆, where twist or strain can be used to engineer magnetic textures without affecting the critical temperature.