🔬 condensed matter

Thermal phase transitions in a mixed-spin Ising model on the Lieb lattice: Exact results beyond zero magnetic field

This paper presents an exact analytical solution and Monte Carlo verification for a mixed-spin Ising model on a Lieb lattice, revealing a rich phase diagram with ferrimagnetic, disordered, and ferromagnetic phases, including dome-shaped discontinuous transitions, Ising-type critical lines, and reentrant behavior, all exactly solvable at finite magnetic fields when the effective field vanishes.

Original authors: Jozef Strecka, Katarina Karlova

Published 2026-07-14

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Jozef Strecka, Katarina Karlova

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a microscopic city built on a special grid called the Lieb lattice. In this city, there are two types of residents living in a very specific arrangement. The "nodal" residents are simple spinners that can only face Up or Down (like a coin on its edge). The "decorating" residents are a bit more complex; they can face Up, face Down, or simply sit still (zero).

Scientists Jozef Strečka and Katarína Karl'ová decided to throw a party for these residents by turning up the heat (temperature) and blowing a strong wind (magnetic field). Their goal? To see how the residents rearrange themselves and whether they suddenly snap into a new order.

The Magic Map

Usually, figuring out how these mixed-up residents behave is like trying to solve a giant, tangled knot of string. But the authors found a magic map (a mathematical trick called the "generalized decoration-iteration transformation"). This map translates the complicated city of mixed residents into a simpler, well-known city of just the simple spinners on a square grid.

Here is the catch: This map only works perfectly if the "effective wind" inside the simplified city is exactly zero. It sounds like a strict rule, but the authors discovered something wild: even if you blow a real, strong wind on the original city, there are specific combinations of wind strength and a special "personal space" force (called uniaxial single-ion anisotropy) where the effective wind inside the map cancels out to zero. When this happens, the complex problem becomes exactly solvable, like turning a locked door into an open hallway.

The Three Neighborhoods

By using this map, the authors mapped out the "ground state"—what the city looks like when it's freezing cold. They found three distinct neighborhoods:

The Ferrimagnetic (FRI) Neighborhood: Here, the simple spinners face one way, and the complex residents face the opposite way. It's a bit of a tug-of-war, but they are locked in a stable, ordered dance.
The Disordered (DP) Neighborhood: The complex residents have decided to "sit still" (zero spin) because the "personal space" force is too strong. The simple spinners are all facing the wind, but the complex ones are just chilling.
The Ferromagnetic (FM) Neighborhood: The wind is so strong that everyone—simple and complex residents alike—faces the same direction. Total alignment!

The paper shows that the boundaries between these neighborhoods are sharp lines. If you cross them, the residents snap instantly from one arrangement to another.

The Dome of Change

When the authors turned up the heat, they didn't just see a smooth slide from one neighborhood to another. They found a dome-shaped wall of sudden changes.

Imagine a hill made of sudden jumps. If you walk up this hill (increasing temperature), you might suddenly jump from the "Disordered" neighborhood to the "Ferrimagnetic" one, and then jump right back out again. This is called a discontinuous phase transition. It's like stepping off a curb and landing on a different floor instantly, rather than walking up a ramp.

This dome is capped by a critical line. If you walk along the very top edge of this dome, the sudden jumps disappear and turn into a smooth, continuous flow. The authors proved that these smooth transitions belong to a famous family of behaviors known as the Ising universality class, the same family as the classic models Onsager solved decades ago.

The "Re-Entrance" Surprise

The most playful part of the discovery happens in a very narrow, tricky zone where the "personal space" force is just slightly less than 2.0 times the interaction strength ( $D/J \lesssim 2$ ).

In this specific spot, the authors found a re-entrant phase transition. This is a mouthful, but the behavior is bizarre:

Start at a low temperature: The city is in the Disordered (DP) state.
Warm it up a little: Suddenly, it snaps into the Ferrimagnetic (FRI) state.
Warm it up more: It snaps back into the Disordered (DP) state.

It's like putting on a jacket because it's cold, taking it off because it's getting warmer, and then putting it back on because it's too warm? No, that doesn't make sense in real life, but in this magnetic city, the residents literally leave their ordered dance, join the party, and then leave the party again just because the temperature changed. The authors identified two consecutive jumps: DP $\rightarrow$ FRI $\rightarrow$ DP.

How Do We Know?

The authors didn't just guess this.

The Math: They used the magic map to derive exact analytical predictions. This means they solved the equations perfectly, not approximately. They proved that the re-entrance and the critical points exist mathematically.
The Simulation: To make sure their math wasn't playing tricks, they ran classical Monte Carlo simulations. Think of this as a supercomputer playing the game millions of times to see what actually happens. They simulated a city with 14,400 simple spins (which corresponds to 43,200 total spins in the mixed model).
The Match: The simulation results matched the exact math perfectly. The "re-entrant" jumps appeared in the simulation exactly where the math said they would.

What They Didn't Find

It's important to note what the paper doesn't say. They didn't find a smooth, gradual change everywhere. They explicitly showed that the transition between the Ferrimagnetic and Ferromagnetic phases is not a sudden jump; it's a smooth crossover. The "dome" of sudden jumps only exists between the Ferrimagnetic and Disordered phases. Also, they didn't find these re-entrant jumps everywhere; they only happen in that tiny, specific slice of the parameter space where the anisotropy is just under 2.0.

The Big Picture

The authors also noticed a funny coincidence. A completely different model, involving quantum spins on a diamond-shaped grid, behaves almost exactly like this one. Even though the microscopic residents are different (one uses quantum mechanics, the other uses classical spins), their "magic maps" lead to the exact same simplified city. This suggests that the weird re-entrant behavior isn't just a fluke of one specific model, but a fundamental feature of how these types of magnetic cities behave.

In short, by using a clever mathematical trick, the authors proved that a mixed-spin magnetic system can undergo a "re-entrant" dance, jumping into order and then back out of it, all while being exactly solvable even when a magnetic field is present. And thanks to the simulations, we know this isn't just a theory—it's a real, predicted phenomenon in this mathematical world.

Technical Summary: Thermal Phase Transitions in a Mixed-Spin Ising Model on the Lieb Lattice

Problem Statement
The paper investigates the ground-state and finite-temperature properties of a mixed spin-1/2 and spin-1 Ising model situated on a decorated square (Lieb) lattice. The system incorporates uniaxial single-ion anisotropy ( $D$ ) and an external magnetic field ( $h$ ). While the zero-field limit of this model was previously solved, the behavior under finite magnetic fields remains a challenge for exact analytical treatment. The study aims to determine the exact phase diagram, specifically focusing on the interplay between the magnetic field and anisotropy, and to identify the nature of thermal phase transitions (continuous vs. discontinuous) and potential reentrant phenomena.

Methodology
The authors employ the generalized decoration-iteration transformation to map the mixed-spin model on the Lieb lattice exactly onto an effective spin-1/2 Ising model on a square lattice.

Mapping: The partition function of the original mixed-spin system is factorized, allowing the summation over the spin-1 decorating sites to be performed exactly. This results in an effective Hamiltonian for the nodal spin-1/2 sites characterized by an effective nearest-neighbor interaction ( $J_{eff}$ ) and an effective magnetic field ( $h_{eff}$ ).
Exact Solvability Condition: The model becomes exactly solvable for finite applied magnetic fields only within a specific parameter subspace where the effective field vanishes ( $h_{eff} = 0$ ). In this regime, the system maps to the exactly solved zero-field Ising model on a square lattice.
Validation: The exact analytical predictions are independently verified using classical Monte Carlo simulations of the effective spin-1/2 model. The simulations utilize the Metropolis algorithm on a square lattice of 14,400 spins, with results mapped back to the original mixed-spin variables using exact statistical relations derived from the transformation.

Key Contributions and Results

Ground-State Phase Diagram:
The zero-temperature analysis identifies three distinct phases:
- Ferrimagnetic (FRI): Characterized by antiparallel alignment of decorating spins ( $S_i=1$ ) and nodal spins ( $\sigma_i=-1/2$ ).
- Disordered (DP): Characterized by nonmagnetic decorating spins ( $S_i=0$ ) and fully polarized nodal spins ( $\sigma_i=1/2$ ).
- Ferromagnetic (FM): Characterized by full alignment of all spins with the external field ( $S_i=1, \sigma_i=1/2$ ).
  The boundaries between these phases are determined by comparing ground-state energies, meeting at a triple point at $[D/J, h/J] = [3, 4]$ .
Finite-Temperature Phase Transitions:
The study reveals that finite-temperature phase transitions occur exclusively along the ground-state coexistence line between the FRI and DP phases, where the condition $h_{eff} = 0$ is satisfied.
- Discontinuous Transitions: A dome-shaped surface of discontinuous thermal phase transitions separates the FRI and DP phases.
- Continuous Transitions: This dome is terminated by a line of Ising-type critical points where the transitions become continuous. These critical points belong to the universality class of the two-dimensional Ising model.
- Critical Condition: The exact critical temperature is derived analytically for the subspace where $h_{eff}=0$ .
Reentrant Phase Transitions:
A significant finding is the identification of a narrow parameter region ( $D/J \lesssim 2$ ) where thermally induced reentrant phase transitions occur. In this regime, as temperature increases, the system undergoes a sequence of transitions: Disordered Phase (DP) $\to$ Ferrimagnetic Phase (FRI) $\to$ Disordered Phase (DP). This implies the FRI phase is stabilized only within an intermediate temperature range, enclosed by the DP phase at both lower and higher temperatures.
Monte Carlo Verification:
The Monte Carlo simulations confirm the analytical predictions, reproducing the magnetization jumps associated with discontinuous transitions, the critical behavior near the terminal points, and the specific temperature dependencies of the reentrant transitions.

Significance and Claims
The paper claims to provide exact analytical results for a mixed-spin Ising model on a Lieb lattice in the presence of a finite magnetic field, a regime where general exact solutions are typically unavailable. The authors emphasize that the exact solvability is achieved by identifying the specific parameter manifold where the effective field vanishes, despite the presence of a non-zero external field.

The work highlights the utility of the decoration-iteration transformation in solving complex mixed-spin systems. Furthermore, the authors draw a quantitative similarity between their results and those of the spin-1/2 Ising-Heisenberg model on a diamond-decorated square lattice. They argue that despite microscopic differences (classical spin-1 vs. quantum Heisenberg dimers), both models map to nearly identical effective spin-1/2 Ising models, explaining the shared occurrence of reentrant transitions and similar phase diagrams. The study suggests that the mixed-spin Ising model on the Lieb lattice serves as a valuable, exactly solvable benchmark for understanding the complex critical behavior found in more intricate hybrid classical-quantum spin systems.