Resolution of the Two-Dimensional Ferromagnetic Spin-3/2 Ising Model via Cluster Growth

We propose a computational methodology based on a hierarchical cluster growth process to solve spin-3/2 Ising models efficiently. The method circumvents the exponential complexity ((4^{N})) of the canonical ensemble partition function by iteratively constructing finite magnetic clusters of size (N_g), where the effective spin state of a site in generation (g+1) is determined by the local magnetization of a cluster from generation (g). This approach, which shares conceptual ground with effective field theories, allows the study of systems of effectively very large size (N = N_0 (N_g)^{g}). We apply the formalism to the ferromagnetic spin-3/2 Ising model on a honeycomb lattice, modeling the monolayer CrI$3$, a prototypical two-dimensional Ising magnet. The model, calibrated using the experimental transition temperature ((T{c} \simeq 45) K), successfully reproduces key experimental features: the temperature dependence of the magnetization (m(T)), including its inflection point, and the broadened peak in the specific heat (c_v(T)). We also compute the entropy (s(T)), finding a finite residual value at low temperatures consistent with the system’s double degeneracy. Our results demonstrate that this hierarchical cluster method provides a quantitatively accurate and computationally efficient framework for studying complex magnetic systems.

💡 Research Summary

The authors introduce a hierarchical cluster‑growth scheme to tackle the notoriously difficult spin‑3/2 Ising model on a two‑dimensional lattice. Traditional canonical‑ensemble approaches require handling 4^N microstates, which quickly becomes intractable for realistic system sizes. By constructing a series of “generations” (g = 0, 1, 2, …) of clusters, the method replaces each spin in generation g + 1 with the average magnetization of a finite cluster from generation g. Consequently, only 4^{N_g} configurations need to be evaluated at any stage, while the effective total number of spins scales as N = N₀(N_g)^g, allowing the study of macroscopic domains with modest computational effort.

The zeroth generation (g = 0) is a microscopic cluster of N₀ spins (chosen as 4 sites). Its Hamiltonian is the usual nearest‑neighbour Ising form H(0)=‑J(0)∑⟨i,j⟩S_i^z S_j^z with S_i^z∈{±3/2, ±1/2}. For higher generations the spin variables are replaced by local magnetizations m_i(g) and the effective Hamiltonian becomes H(g)=‑J(g)∑⟨i,j⟩m_i(g)m_j(g). The dimensionless coupling K(g)=J(g)/k_BT is renormalized recursively via J(g)=a^g J(0), where the scale factor a is a tunable parameter that controls how the effective interaction evolves with the generation. The magnetization update rule is taken as F(m)=m/m_sat with m_sat=1.5 (the saturation magnetization of a spin‑3/2 ion), ensuring that each generation preserves the physical bound 0 ≤ |m| ≤ m_sat.

To connect the abstract cluster size to a real material, the authors equate the effective particle number N = N₀(N_g)^g_max with the number of magnetic ions in a given volume (n mol × N_A). For a monolayer of CrI₃, the experimentally measured nearest‑neighbour exchange is J(0)≈6.91 meV, which the authors translate into an energy scale γ(0)=J(0)/k_B≈10² K. Using N₀=N_g=4 and a maximum generation g_max=10, the effective system contains roughly 10⁶ spins, comparable to a macroscopic magnetic domain.

Numerical results show that the zeroth‑generation cluster displays only smooth finite‑size behavior, whereas generations g ≥ 1 exhibit a sharp drop in magnetization m(T) with a clear inflection point, reproducing the experimentally observed Curie transition near 45 K. The specific heat c_v(T) develops a broadened peak whose height and width depend on both a and g_max, reflecting finite‑size rounding of the critical singularity. The entropy s(T) approaches a residual value of k_B ln 2 as T→0, consistent with the double degeneracy of the spin‑3/2 ground state (±3/2 versus ±1/2). By locating the value a_v at which the internal energy u(a) shows an inflection, the authors identify a proxy for the critical point, analogous to a fixed point in renormalization‑group theory.

The paper also explores a generalized magnetization recurrence F_α(m)=(m/m_sat)^α. Physical constraints (F_α(0)=0, F_α(m_sat)=1, monotonicity) restrict α to the interval (0, 2); the linear case α = 1 is found to be the simplest and most consistent with experimental data.

Compared with Monte‑Carlo or conventional real‑space RG, the hierarchical cluster method does not integrate out degrees of freedom nor perform explicit length rescaling; it works directly with thermodynamic observables at each scale. This makes it computationally cheap while still capturing essential critical behavior, and it can be readily extended to include longer‑range interactions, anisotropies, or higher‑dimensional lattices.

In summary, the hierarchical cluster‑growth framework provides a novel, efficient, and quantitatively accurate tool for studying spin‑3/2 Ising systems and, by extension, other complex magnetic materials. Its successful application to monolayer CrI₃ demonstrates its ability to reproduce key experimental signatures—magnetization curve, specific‑heat peak, and low‑temperature entropy—while dramatically reducing the computational burden associated with the exponential state space of high‑spin Ising models.