Anisotropic Lipkin–Meshkov–Glick Model

• The anisotropic Lipkin–Meshkov–Glick model is a quantum many-body system featuring competing anisotropic spin interactions and distinctive phase transitions.
• Investigations reveal critical phenomena including both first-order and second-order quantum phase transitions driven by anisotropy and system size effects.
• Exact and mean-field methods provide actionable insights into entanglement dynamics, spectral structure, and applications in quantum technologies.

The anisotropic Lipkin–Meshkov–Glick (LMG) model is a quantum many-body system describing collective spins with competing, anisotropic infinite-range interactions in the $x$-$y$ plane, subject to a transverse field. Its anisotropic variant, realized by distinct couplings in $x$ and $y$ directions, exhibits a wealth of critical and dynamical phenomena, integrating perspectives from quantum phase transitions, integrability, magnetism, and non-equilibrium dynamics. The model is analytically tractable for small system sizes and solvable in the thermodynamic limit via mean-field, large-$N$, and exact (Richardson–Gaudin) techniques. It is pivotal both as a benchmark for critical phenomena in collective spin systems and as a theoretical framework for quantum technologies involving entanglement and collective dynamics.

1. Hamiltonian Forms and Anisotropy Structure

The general anisotropic LMG Hamiltonian for an ensemble of $N$ spin-$1/2$ particles (or, equivalently, a collective spin-$j=N/2$) is

$$H_{\mathrm{LMG}} = g\,J^z - \frac{\gamma_x}{2j} (J^x)^2 - \frac{\gamma_y}{2j}(J^y)^2$$

where $J^\alpha = \sum_{i=1}^N \sigma^\alpha_i/2$ are collective spin operators, $g$ is the transverse field, and $\gamma_x, \gamma_y$ control the interaction strengths in $x$ and $y$ directions respectively (Sorokin et al., 2014). Alternative forms involve $J_+^2 + J_-^2$ and $J_+J_- + J_-J_+$ or, in particle representation,

$$H_{\mathrm{LMG}} = E_0 J_z + \frac{V}{2}(J_+^2 + J_-^2) + \frac{W}{2}(J_+ J_- + J_- J_+)$$

with $V\neq W$ yielding $x$-$y$ anisotropy (Co' et al., 2018). The equivalence between matrix and collective spin representations allows direct computation of the full spectrum for modest $N$ and rigorous mean-field analysis for $N\to\infty$.

In the two-spin case ($N=2$), the Hamiltonian becomes

$$H = -\frac{J}{4}[\sigma^1_x \sigma^2_x + \gamma\,\sigma^1_y \sigma^2_y] - \frac{h}{2}(\sigma^1_z + \sigma^2_z) - \frac{J}{4}(1+\gamma)$$

where $\gamma$ encodes the $x/y$ anisotropy and $h$ is the transverse field (Çakmak et al., 2015).

2. Exact Solutions for Small-$N$ and Richardson–Gaudin Integrability

For $N=2$ or $N=3$, analytical diagonalization yields explicit eigenvalues and eigenvectors, showing how anisotropy drives spectral crossings. For $N=2$: $$E_\pm^{[2]}(\nu,\omega) = -\omega \pm \sqrt{1+\nu^2},\quad E_0^{[2]}(\nu,\omega) = -2\omega$$ where $\nu=|V|/E_0$, $\omega=|W|/E_0$ (Co' et al., 2018). Level crossings of the ground and first excited states occur at the critical anisotropy $\omega_c = \sqrt{1+\nu^2}$, a result that generalizes to arbitrary $N$ with $\omega_{2,1}^{[N]}(\nu) = \sqrt{\nu^2 + 1/(N-1)^2}$.

In large-$N$ or when formulating the model in terms of Schwinger bosons, the anisotropic LMG model becomes a specific case of the SU(1,1) Richardson–Gaudin integrable models, possessing $M$ "pairon" roots $\{e_\alpha\}$ determined by non-linear Bethe-equations. The exact eigenstates are constructed as

$$|\Psi\rangle = \prod_{\alpha=1}^M \left( \frac{a^\dagger a^\dagger}{e_\alpha + t } + \frac{b^\dagger b^\dagger}{e_\alpha - t} \right) |\nu_a, \nu_b\rangle$$

where $a, b$ are Schwinger bosons. The model displays first-, second-, and rare third-order quantum phase transitions mapped to the dynamics and topology changes of the pairons in the complex plane (Lerma et al., 2012).

3. Phase Diagram, Quantum Phase Transitions, and Criticality

The anisotropic LMG model features both continuous and discontinuous quantum phase transitions, depending on the symmetry and anisotropy. In the isotropic limit ($\gamma_x = \gamma_y$ or $V=W$), the system exhibits a second-order quantum phase transition at $g_c = \gamma_x$. For generic anisotropy and finite $N$, a sequence of ground-state level crossings occurs as $\omega$ increases—each corresponding to first-order transitions—with the critical anisotropy tending to zero as $N\to\infty$.

A mean-field analysis in the thermodynamic limit delineates three phases ("paramagnetic", "$x$-ferro", "$y$-ferro"), separated by lines $g_c^x = \gamma_x$, $g_c^y = \gamma_y + 2\kappa$ (for ring-coupled networks), featuring both exponential and algebraic ground-state degeneracies (Sorokin et al., 2014). At the critical points, the lowest excitation gap closes linearly or shows jumps, with the full excitation spectrum available from Bogoliubov analysis.

Classical (mean-field) dynamics additionally shows dynamical phase transitions: at separatrix crossings (saddles in the energy surface), the time-averaged magnetization becomes singular (logarithmic scaling), and in the anisotropic case, these bifurcations become richer, with distinct critical exponents for different observables (Yu, 2023).

4. Spectral Structure, Level Crossings, and Observables

Spectral features of the anisotropic LMG model are controlled by the interplay of anisotropy and system size:

• For $W=0$ (isotropic), the gap closes continuously at the critical point without degeneracies for finite $N$, and the transition is of second order.
• When $W\neq 0$, a cascade of exact level crossings occurs; these delimit genuine first-order transitions at finite $N$.
• In the thermodynamic limit, infinitesimal anisotropy ($\omega\to 0^+$) leads to a discontinuous change in the ground state.

For $N=2$, the full spectrum is available analytically, and the effect of anisotropy is explicit in the eigenvalues. For $N\gg 1$, spectral crossings, avoided crossings, and collapse of the gap can be tracked numerically or via Heine–Stieltjes ODE techniques (Co' et al., 2018, Lerma et al., 2012).

Observables such as ground-state magnetization and entanglement (concurrence) are computable in closed form, particularly for $N=2$, with entanglement peaking at the anisotropic crossing points and vanishing in the Ising or large field limits (Çakmak et al., 2015).

5. Classical Dynamics, Non-equilibrium Behavior, and Dynamical Phase Transitions

At the classical level, the system evolves according to coupled non-linear ordinary differential equations for $(S_x, S_y, S_z)$, exhibiting both librational (self-trapped) and rotational regimes—separated by dynamical phase transitions associated with saddle crossings on the Bloch sphere. The addition of two nonlinear interactions ($J_x \neq J_y$) causes the criticality in long-time averaged observables to display both logarithmic and non-logarithmic scaling, depending on the regime (Yu, 2023).

The exact solution of the classical dynamics in terms of Jacobi elliptic functions provides a benchmark for quantum dynamical phase transition (DPT) analysis and underpins semiclassical approximations (truncated Wigner, Holstein–Primakoff). DPT signatures, such as nonanalyticities in the Loschmidt echo and sudden jumps in infinite-time averages, are directly linked to classical separatrix crossings.

6. Applications: Single-Molecule Magnets, Quantum Engines, and Networks

The anisotropic LMG framework extends to the description of superparamagnetic single-molecule magnets (e.g., Fe$_8$), where macroscopic spins ($S\gg 1$) experience anisotropic energy barriers and tunneling (Campos et al., 2011). In this context, the Hamiltonian incorporates axial ($D S_z^2$) and rhombic ($E(S_x^2-S_y^2)$) anisotropies, leading to rich level-crossing and avoided-crossing phenomena, notably controlling quantum tunneling of magnetization rates.

As a working substance in quantum Otto cycles, the two-spin anisotropic LMG model reveals regimes of operation unattainable for noninteracting or even isotropic systems, allowing cooperative work enhancement, high-efficiency operation, and the approach to Carnot efficiency bounds by judicious tuning of anisotropy parameters (Çakmak et al., 2015).

In networked settings, anisotropic LMG models placed on a ring or more general graphs support collective quantum phases, exotic ground-state degeneracies, and modified critical fields dictated by network topology and interaction anisotropy (Sorokin et al., 2014).

7. Integrable Structure, Numerical Methods, and Quantum Information Perspectives

The integrability of the anisotropic LMG model via SU(1,1) Richardson–Gaudin mappings allows for exact characterization of the spectrum and phase diagram, especially in terms of the dynamics and complexification of the Bethe roots ("pairons"). This perspective highlights connections between avoided spectral crossings, true level crossings, higher-order QPTs, and topologically nontrivial many-body wavefunctions (Lerma et al., 2012).

The generalized Heine–Stieltjes ODE approach effectively computes the full set of Bethe roots for moderate system sizes, while large-scale diagonalization or mean-field/Bogoliubov techniques are used for larger $N$.

Current research also emphasizes the relevance of the anisotropic LMG model for engineered quantum platforms such as cold-atom ensembles, superconducting circuits, and quantum simulators, especially in realizing and probing dynamical criticality, collective entanglement, and fast quantum state switching mediated by tunable anisotropy (Yu, 2023).