Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2610.05057

Onset of a quantum Fisher-Selke sequence in the transverse-field ANNNI model and its bilayer

Tuan-Vu Truong, Vu-Quoc-Minh Nguyen, Hoang-Long Nguyen, Thanh-Tien Nguyen

quant-phcond-mat.str-elSubmitted Oct 4, 2026• Announced Oct 6, 2026• 0 citations
arXivHTMLPDF

Abstract

We study the transverse-field axial next-nearest-neighbor Ising (ANNNI) model on the square lattice and on a bilayer. At the classical multiphase point $κ=J_2/J_0=1/2$, every configuration of straight domain walls that separates domains of at least two sites is a ground state. The walls are rigid, so perturbation theory in the transverse field $g=Γ/J_0$ can be carried out for each configuration separately, and we take it to order $g^8$. Quantum fluctuations select $\langle 3\rangle$ ($\uparrow\uparrow\uparrow\downarrow\downarrow\downarrow$) at order $g^2$, because spins next to a wall gain the most energy from virtual flips. The structures $\langle 23\rangle$, $\langle 223\rangle$ and $\langle 2223\rangle$ then open between $\langle 3\rangle$ and the antiphase $\langle 2\rangle$ at orders $g^4$, $g^6$ and $g^8$, in the order of the Fisher-Selke sequence of the classical three-dimensional model. A second layer suppresses each successive phase more strongly, beyond the stiffening of the local fields. We test the two widest phases with stochastic series expansion quantum Monte Carlo, comparing the free energies of the competing ordered states by thermodynamic integration, since annealing freezes in the modulation that forms at the ordering transition. The simulations reproduce the $\langle 3\rangle$ wedge to about 1%, and in a pre-registered test $\langle 23\rangle$ is stable at six of seven points, with the seventh resolved by post-hoc runs.

Related papers

Ground state properties of the one-dimensional axial next-nearest-neighbor Ising model in a transverse field

Yun-Tong Yang, Chen, Fu-Zhou, Luo, Hong-Gang

Describing and understanding the consequences of competing interactions remains profoundly challenging in both classical and quantum systems, as it is difficult to identify suitable order parameters, thereby hindering the characterization of certain phases such as the floating phase found in the one-dimensional axial next-nearest-neighbor Ising (ANNNI) model in the presence of the frustration interactions between the nearest and next-nearest neighbor sites. In this work, we employ the pattern picture to explore the frustration physics in such a model. This picture has been comprehensively detailed in our previous article [Yang and Luo, Phys. Rev. E \textbf{112}, 044102 (2025)]. Here, we apply it to the ANNNI model with periodic boundary conditions, considering system sizes ranging from $L=16$ to $L=128$. Our results demonstrate that the ground state of the system comprises four phases: ferromagnetic, paramagnetic, floating, and $\langle 2,2 \rangle$ antiphase. The transition from the ferromagnetic to the paramagnetic phase is continuous, analogous to that in the transverse-field Ising model, while the transition from the floating phase to the antiphase is first order. Furthermore, the floating phase exhibits particularly intriguing characteristics: states with distinct domain structures emerge successively as the frustration parameter increases. As the system size grows, this succession becomes progressively denser, leading to the reasonable inference that it eventually approaches a continuous variation in the thermodynamic limit. To validate the effectiveness of our picture, we computed the second derivative of the ground-state energy, which exhibits multiple dips within the floating phase$-$consistent with the pattern language.

arXiv:2301.08891•Jan 21, 2023•2 citations

Floating phase in the one-dimensional transverse axial next-nearest-neighbor Ising model

Anjan Kumar Chandra, Subinay Dasgupta

To study the ground state of an axial next-nearest-neighbor Ising chain under transverse field as a function of frustration parameter kappa and field strength Gamma, we present here two different perturbative analyses. In one, we consider the (known) ground state at kappa=0.5 and Gamma=0 as the unperturbed state and treat an increase of the field from 0 to Gamma coupled with an increase of kappa from 0.5 to 0.5+rGamma/J as perturbation. The first-order perturbation correction to eigenvalue can be calculated exactly and we could conclude that there are only two phase-transition lines emanating from the point kappa=0.5, Gamma=0. In the second perturbation scheme, we consider the number of domains of length 1 as the perturbation and obtain the zeroth-order eigenfunction for the perturbed ground state. From the longitudinal spin-spin correlation, we conclude that floating phase exists for small values of transverse field over the entire region intermediate between the ferromagnetic phase and antiphase.

arXiv:cond-mat/0612144•Physical Review E•Feb 7, 2007•28 citations

Kosterlitz-Thouless phase transition of the axial next-nearest-neighbor Ising model in two dimensions

Takayuki Shirakura, Fumitaka Matsubara, Naoki Suzuki

The spin structure of an axial next-nearest-neighbor Ising (ANNNI) model in two dimensions (2D) is a renewed problem because different Monte Carlo (MC) simulation methods predicted different spin orderings. The usual equilibrium simulation predicts the occurrence of a floating incommensurate (IC) Kosterlitz-Thouless (KT) type phase, which never emerges in non-equilibrium relaxation (NER) simulations. In this paper, we first examine previously published results of both methods, and then investigate a higher transition temperature ${T}_{c1}$ between the IC and paramagnetic phases. In the usual equilibrium simulation, we calculate the chain magnetization on larger lattices (up to $512\ifmmode\times\else\texttimes\fi{}512$ sites) and estimate ${T}_{c1}\ensuremath{\approx}1.16J$ with frustration ratio $\ensuremath{\kappa}(\ensuremath{\equiv}\ensuremath{-}{J}_{2}/{J}_{1})=0.6$. We examine the nature of the phase transition in terms of the Binder ratio ${g}_{L}$ of spin overlap functions and the correlation-length ratio $\ensuremath{\xi}/L$. In the NER simulation, we observe the spin dynamics in equilibrium states by means of an autocorrelation function and also observe the chain magnetization relaxations from the ground and disordered states. These quantities exhibit an algebraic decay at $T\ensuremath{\lesssim}1.17J$. We conclude that the two-dimensional ANNNI model actually admits an IC phase transition of the KT type.

arXiv:1406.6464•Physical Review B•Oct 8, 2014•8 citations

Ising models with competing axial interactions in transverse fields

Parongama Sen, Bikas K. Chakrabarti

The mean-field phase diagram of a three-dimensional axial-next-nearest-neighbor Ising (ANNNI) model in a transverse field has been studied here. The dynamical modes are found by using the random-phase approximation and the phase diagram is also obtained from the softening condition of these modes. For a one-dimensional ANNNI model in a transverse field, these dynamical modes are obtained here more accurately using a Jordan-Wigner-type transformation (with some approximations). Softening of these modes gives the zero-temperature phase diagram in one dimension.

Physical review. B, Condensed matter•Jul 1, 1989•18 citations

Quantum fidelity approach to the ground-state properties of the one-dimensional axial next-nearest-neighbor Ising model in a transverse field

O. F. de Alcantara Bonfim, B. Boechat, J. Florencio

In this work we analyze the ground-state properties of the s=1/2 one-dimensional axial next-nearest-neighbor Ising model in a transverse field using the quantum fidelity approach. We numerically determined the fidelity susceptibility as a function of the transverse field B_{x} and the strength of the next-nearest-neighbor interaction J_{2}, for systems of up to 24 spins. We also examine the ground-state vector with respect to the spatial ordering of the spins. The ground-state phase diagram shows ferromagnetic, floating, and 〈2,2〉 phases, and we predict an infinite number of modulated phases in the thermodynamic limit (L→∞). Paramagnetism only occurs for larger magnetic fields. The transition lines separating the modulated phases seem to be of second order, whereas the line between the floating and the 〈2,2〉 phases is possibly of first order.

arXiv:1705.09679•Physical review. E•Oct 17, 2017•10 citations

Quantum Monte Carlo study of the bond- and site-diluted transverse-field Ising model

Calvin Krämer, Max Hörmann, Kai Phillip Schmidt

We study the transverse-field Ising model on a square lattice with bond- and site-dilution at zero temperature by stochastic series expansion quantum Monte Carlo simulations. Tuning the transverse field $h$ and the dilution $p$, the quantum phase diagram of both models is explored. Both quantum phase diagrams show long-range order for small $h$ and small $p$. The ordered phase of each is separated from the disordered (quantum) Griffiths phase by second-order phase transitions on two critical lines touching at a multi-critical point. Using Binder ratios we locate quantum critical points with high accuracy. The order-parameter critical exponent $β$ and the average correlation-length exponent $ν_{\mathrm{av}}$ are determined along the critical lines and at the multi-critical points via finite-size scaling. We find three internally consistent sets of critical exponents and compare them with potentially connected universality classes. The quantum Griffiths phase in the vicinity of the phase transition lines is analyzed through the local susceptibility. Our results indicate that activated scaling occurs not only at the percolation transition, but also at the phase transition line for $p$ smaller than the percolation threshold $p_c$.

arXiv:2505.07627•Physical review. B./Physical review. B•Jul 8, 2025•1 citation
We study the transverse-field axial next-nearest-neighbor Ising (ANNNI) model on the square lattice and on a bilayer. At the classical multiphase point κ=J2/J0=1/2κ=J_2/J_0=1/2κ=J2​/J0​=1/2, every configuration of straight domain walls that separates domains of at least two sites is a ground state. The walls are rigid, so perturbation theory in the transverse field g=Γ/J0g=Γ/J_0g=Γ/J0​ can be carried out for each configuration separately, and we take it to order g8g^8g8. Quantum fluctuations select ⟨3⟩\langle 3\rangle⟨3⟩ (↑↑↑↓↓↓\uparrow\uparrow\uparrow\downarrow\downarrow\downarrow↑↑↑↓↓↓) at order g2g^2g2, because spins next to a wall gain the most energy from virtual flips. The structures ⟨23⟩\langle 23\rangle⟨23⟩, ⟨223⟩\langle 223\rangle⟨223⟩ and ⟨2223⟩\langle 2223\rangle⟨2223⟩ then open between ⟨3⟩\langle 3\rangle⟨3⟩ and the antiphase ⟨2⟩\langle 2\rangle⟨2⟩ at orders g4g^4g4, g6g^6g6 and g8g^8g8, in the order of the Fisher-Selke sequence of the classical three-dimensional model. A second layer suppresses each successive phase more strongly, beyond the stiffening of the local fields. We test the two widest phases with stochastic series expansion quantum Monte Carlo, comparing the free energies of the competing ordered states by thermodynamic integration, since annealing freezes in the modulation that forms at the ordering transition. The simulations reproduce the ⟨3⟩\langle 3\rangle⟨3⟩ wedge to about 1%, and in a pre-registered test ⟨23⟩\langle 23\rangle⟨23⟩ is stable at six of seven points, with the seventh resolved by post-hoc runs.
Describing and understanding the consequences of competing interactions remains profoundly challenging in both classical and quantum systems, as it is difficult to identify suitable order parameters, thereby hindering the characterization of certain phases such as the floating phase found in the one-dimensional axial next-nearest-neighbor Ising (ANNNI) model in the presence of the frustration interactions between the nearest and next-nearest neighbor sites. In this work, we employ the pattern picture to explore the frustration physics in such a model. This picture has been comprehensively detailed in our previous article [Yang and Luo, Phys. Rev. E \textbf{112}, 044102 (2025)]. Here, we apply it to the ANNNI model with periodic boundary conditions, considering system sizes ranging from L=16L=16L=16 to L=128L=128L=128. Our results demonstrate that the ground state of the system comprises four phases: ferromagnetic, paramagnetic, floating, and ⟨2,2⟩\langle 2,2 \rangle⟨2,2⟩ antiphase. The transition from the ferromagnetic to the paramagnetic phase is continuous, analogous to that in the transverse-field Ising model, while the transition from the floating phase to the antiphase is first order. Furthermore, the floating phase exhibits particularly intriguing characteristics: states with distinct domain structures emerge successively as the frustration parameter increases. As the system size grows, this succession becomes progressively denser, leading to the reasonable inference that it eventually approaches a continuous variation in the thermodynamic limit. To validate the effectiveness of our picture, we computed the second derivative of the ground-state energy, which exhibits multiple dips within the floating phase−-−consistent with the pattern language.
To study the ground state of an axial next-nearest-neighbor Ising chain under transverse field as a function of frustration parameter kappa and field strength Gamma, we present here two different perturbative analyses. In one, we consider the (known) ground state at kappa=0.5 and Gamma=0 as the unperturbed state and treat an increase of the field from 0 to Gamma coupled with an increase of kappa from 0.5 to 0.5+rGamma/J as perturbation. The first-order perturbation correction to eigenvalue can be calculated exactly and we could conclude that there are only two phase-transition lines emanating from the point kappa=0.5, Gamma=0. In the second perturbation scheme, we consider the number of domains of length 1 as the perturbation and obtain the zeroth-order eigenfunction for the perturbed ground state. From the longitudinal spin-spin correlation, we conclude that floating phase exists for small values of transverse field over the entire region intermediate between the ferromagnetic phase and antiphase.
The spin structure of an axial next-nearest-neighbor Ising (ANNNI) model in two dimensions (2D) is a renewed problem because different Monte Carlo (MC) simulation methods predicted different spin orderings. The usual equilibrium simulation predicts the occurrence of a floating incommensurate (IC) Kosterlitz-Thouless (KT) type phase, which never emerges in non-equilibrium relaxation (NER) simulations. In this paper, we first examine previously published results of both methods, and then investigate a higher transition temperature Tc1{T}_{c1}Tc1​ between the IC and paramagnetic phases. In the usual equilibrium simulation, we calculate the chain magnetization on larger lattices (up to $512\ifmmode\times\else\texttimes\fi{}512$ sites) and estimate ${T}_{c1}\ensuremath{\approx}1.16J$ with frustration ratio $\ensuremath{\kappa}(\ensuremath{\equiv}\ensuremath{-}{J}_{2}/{J}_{1})=0.6$. We examine the nature of the phase transition in terms of the Binder ratio gL{g}_{L}gL​ of spin overlap functions and the correlation-length ratio $\ensuremath{\xi}/L$. In the NER simulation, we observe the spin dynamics in equilibrium states by means of an autocorrelation function and also observe the chain magnetization relaxations from the ground and disordered states. These quantities exhibit an algebraic decay at $T\ensuremath{\lesssim}1.17J$. We conclude that the two-dimensional ANNNI model actually admits an IC phase transition of the KT type.
The mean-field phase diagram of a three-dimensional axial-next-nearest-neighbor Ising (ANNNI) model in a transverse field has been studied here. The dynamical modes are found by using the random-phase approximation and the phase diagram is also obtained from the softening condition of these modes. For a one-dimensional ANNNI model in a transverse field, these dynamical modes are obtained here more accurately using a Jordan-Wigner-type transformation (with some approximations). Softening of these modes gives the zero-temperature phase diagram in one dimension.
In this work we analyze the ground-state properties of the s=1/2 one-dimensional axial next-nearest-neighbor Ising model in a transverse field using the quantum fidelity approach. We numerically determined the fidelity susceptibility as a function of the transverse field B_{x} and the strength of the next-nearest-neighbor interaction J_{2}, for systems of up to 24 spins. We also examine the ground-state vector with respect to the spatial ordering of the spins. The ground-state phase diagram shows ferromagnetic, floating, and 〈2,2〉 phases, and we predict an infinite number of modulated phases in the thermodynamic limit (L→∞). Paramagnetism only occurs for larger magnetic fields. The transition lines separating the modulated phases seem to be of second order, whereas the line between the floating and the 〈2,2〉 phases is possibly of first order.
We study the transverse-field Ising model on a square lattice with bond- and site-dilution at zero temperature by stochastic series expansion quantum Monte Carlo simulations. Tuning the transverse field hhh and the dilution ppp, the quantum phase diagram of both models is explored. Both quantum phase diagrams show long-range order for small hhh and small ppp. The ordered phase of each is separated from the disordered (quantum) Griffiths phase by second-order phase transitions on two critical lines touching at a multi-critical point. Using Binder ratios we locate quantum critical points with high accuracy. The order-parameter critical exponent βββ and the average correlation-length exponent νavν_{\mathrm{av}}νav​ are determined along the critical lines and at the multi-critical points via finite-size scaling. We find three internally consistent sets of critical exponents and compare them with potentially connected universality classes. The quantum Griffiths phase in the vicinity of the phase transition lines is analyzed through the local susceptibility. Our results indicate that activated scaling occurs not only at the percolation transition, but also at the phase transition line for ppp smaller than the percolation threshold pcp_cpc​.