The ordering of charges on half-filled hypercubic lattices is investigated numerically, where electroneutrality is ensured by background charges. This system is equivalent to the
s=1/2 Ising lattice model with antiferromagnetic
1/r interaction. The temperature dependences of specific heat, mean staggered occupation, and of a generalized susceptibility indicate continuous order-disorder phase transitions at finite temperatures in two- and three-dimensional systems. In contrast, the susceptibility of the one-dimensional system exhibits singular behavior at vanishing temperature. For the two- and three-dimensional cases, the critical exponents are obtained by means of a finite-size scaling analysis. Their values are consistent with those of the Ising model with short-range interaction, and they imply that the studied model cannot belong to any other known universality class. Samples of up to 1400,
1122, and
223 sites are considered for dimensions 1 to 3, respectively.
The percolation, Ising, and O(
n) models constitute fundamental systems in statistical and condensed matter physics. For short-range-interacting cases, the nature of their phase transitions is well established by renormalization-group theory. However, the universality of the transitions in these models remains elusive when algebraically decaying long-range interactions
∼1/rd+σ are introduced, where
d is the dimensionality and
σ is the decay exponent. Building upon insights from Lévy flight, i.e., long-range simple random walk, we propose three universality diagrams in the
(d,σ) plane for the percolation model, the O(
n) model, and the Fortuin-Kasteleyn Ising model, respectively. The conjectured universality diagrams are consistent with recent high-precision numerical studies and rigorous mathematical results, offering a unified perspective on critical phenomena in systems with long-range interactions.
Determining the threshold value
σ∗ that separates the short-range (SR) and long-range (LR) universality classes in phase transitions remains a controversial issue. While Sak's criterion,
σ∗=2−ηSR, has been widely accepted, recent studies of two-dimensional (2D) models with long-range interactions have challenged it. In this work, we focus on the crossover between LR and SR criticality in several classical 2D statistical models, including the XY, Heisenberg, percolation, and Ising models, whose interactions decay as
1/r2+σ. Our previous simulations for the XY, Heisenberg, and percolation models consistently indicate a universal boundary at
σ∗=2. Here, we complete the picture by performing large-scale Monte Carlo simulations of the 2D LR-Ising model, reaching lattice sizes up to
L=8192. By analyzing the Fortuin-Kasteleyn critical polynomial
Rp, the Binder ratio
Qm, and the anomalous dimension
η, we obtain convergent and self-consistent evidence that the universality class already changes sharply at
σ=2. Taken together, these results establish a unified scenario for LR interacting systems: across all studied models, the crossover from LR to SR universality occurs at
σ∗=2.