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:2312.04468 · Physical review. E

Helicity modulus and chiral symmetry breaking for boundary conditions with finite twist

Gaurav Khairnar, Thomas Vojta

Published Feb 12, 2025• 5 citations
arXivHTMLPDF

Abstract

We study the response of a two-dimensional classical XY model to a finite (noninfinitesimal) twist of the boundary conditions. We use Monte Carlo simulations to evaluate the free energy difference between periodic and twisted-periodic boundary conditions and find deviations from the expected quadratic dependence on the twist angle. Consequently, the helicity modulus (spin stiffness) shows a nontrivial dependence on the twist angle. We show that the deviation from the expected behavior arises because of the mixing of states with opposite chirality which leads to an additional entropy contribution in the quasi-long-range ordered phase. We give an improved prescription for the numerical evaluation of the helicity modulus for a finite twist, and we discuss the spontaneous breaking of the chiral symmetry for the antiperiodic boundary conditions. We also discuss applications to discrete spin systems and some experimental scenarios where boundary conditions with finite twist are necessary.

DOI
10.1103/physreve.111.024114

Related papers

Monte Carlo studies of the finite three-dimensionalXY-model helicity modulus

K. K. Mon

A new logarithmic term in the finite-size helicity modulus has been predicted recently for the three-dimensional classical XY-model universality class under nonperiodic boundary conditions. These theoretical ideas have important implications for the analysis of superfluid-density measurements in confined geometries. We present numerical results from Monte Carlo simulations of the helicity modulus for the three-dimensional classical XY model on a finite simple-cubic lattice, with pinned-surface-spin boundary conditions. Our data are consistent with the presence of an additional logarithmic term and an estimate for the amplitude is given.

Physical review. B, Condensed matter•Aug 1, 1991•2 citations

Critical behavior of the helicity modulus for the classical Heisenberg model

Robert Goodell Brown, Mikael Ciftan

The critical scaling of the helicity modulus of the classical O(3) $3d$ Heisenberg ferromagnet is studied directly. Monte Carlo methods that impose either an antiperiodic boundary condition or a finite twist of definite handedness across otherwise periodic boundaries in one lattice direction are used to measure scale-dependent enthalpy variations in a simple cubic lattice at the ferromagnetic critical temperature. Finite-size scaling is then used to determine the critical exponents ${v}_{E}$ and ${v}_{F}$ for helicity and, by evaluating three independent hyperscaling-linked pairs of $\ensuremath{\nu}$ and $\ensuremath{\alpha}$, to test hyperscaling for this model. It is observed that antiperiodic boundary conditions in particular constrain the lattice to have a nonzero topological charge, establishing a connection between topological charge and helicity in the model.

Physical Review B•Dec 14, 2006•12 citations

Response to a twist in systems with Z p symmetry: The two-dimensional p -state clock model

Yuta Kumano, Koji Hukushima, Yusuke Tomita +1 more

We study response to a twist in the two-dimensional $p$-state clock model, which has the discrete ${Z}_{p}$ symmetry. The response is measured in terms of helicity modulus, which is usually defined with respect to an infinitesimal twist. However, we demonstrate that such a definition is inappropriate for the clock model. The helicity modulus must be defined with respect to a finite, quantized twist which matches the discrete ${Z}_{p}$ symmetry of the model. Numerical study of the appropriately defined helicity modulus resolves controversy over the clock model, showing the existence of two Berezinskii-Kosterlitz-Thouless transitions for $p>4$.

arXiv:1301.6166•Physical Review B•Sep 30, 2013•52 citations

Chiral and spin order in the two-dimensional ±JXYspin glass: Domain-wall scaling analysis

Muriel Ney‐Nifle, H. J. Hilhorst

This is an analytic study of the two-dimensional XY spin glass with \ifmmode\pm\else\textpm\fi{}J disorder. The Hamiltonian has a continuous spin symmetry and a discrete chiral symmetry, and therefore possesses, potentially, two different order parameters and correlation lengths. The cost of breaking the symmetries is probed by comparing the ground-state energy under periodic (P) boundary conditions with the one under antiperiodic (AP) and under reflecting (R) boundary conditions. Two energy differences (domain-wall energies) appear, \ensuremath{\Delta}${\mathit{E}}^{\mathrm{AP}}$ and \ensuremath{\Delta}${\mathit{E}}^{\mathit{R}}$, whose scaling behavior with system size is nontrivially related to the correlation length exponents. For a specific distribution of the \ifmmode\pm\else\textpm\fi{}J disorder we show that the chiral, and spin correlation lengths diverge with the same exponent as T goes down to 0. The common exponent has a common cause, viz. the low reversal energy of domains of chiral variables. For general disorder we give a heuristic argumentation in terms of droplet excitations that leads to spin ordering on a longer, or equal, scale than chiral ordering. These results are in contrast with intepretations of Monte Carlo simulations.

arXiv:cond-mat/9409049•Physical review. B, Condensed matter•Apr 1, 1995•14 citations

Symmetries and boundary conditions with a twist

D'Amico, Irene, de Oliveira, L. N., Zawadzki, K.

Interest in finite-size systems has risen in the last decades, due to the focus on nanotechnological applications and because they are convenient for numerical treatment that can subsequently be extrapolated to infinite lattices.Independently of the envisioned application, special attention must be given to boundary condition, which may or may not preserve the symmetry of the infinite lattice.Here we present a detailed study of the compatibility between boundary conditions and conservation laws.The conflict between open boundary conditions and momentum conservation is well understood, but we examine other symmetries, as well: we discuss gauge invariance, inversion, spin, and particle-hole symmetry and their compatibility with open, periodic, and twisted boundary conditions.In the interest of clarity, we develop the reasoning in the framework of the one-dimensional half-filled Hubbard model, whose Hamiltonian displays a variety of symmetries.Our discussion includes analytical and numerical results.Our analytical survey shows that, as a rule, boundary conditions break one or more symmetries of the infinite-lattice Hamiltonian.The exception is twisted boundary condition with the special torsion Θ = πL/2, where L is the lattice size.Our numerical results for the ground-state energy at half-filling and the energy gap for L = 2-7 show how the breaking of symmetry affects the convergence to the L → ∞ limit.We compare the computed energies and gaps with the exact results for the infinite lattice drawn from the Bethe-Ansatz solution.The deviations are boundary-condition dependent.The special torsion yields more rapid convergence than open or periodic boundary conditions.For sizes as small as L = 7, the numerical results for twisted condition are very close to the L → ∞ limit.We also discuss the ground-state electronic density and magnetization at half filling under the three boundary conditions.

•Jul 17, 2017•12 citations

Helicity modulus in the bilayer XY model by the Monte Carlo worm algorithm

A. Masini, Alessandro Cuccoli, A. Rettori +2 more

The behavior of the helicity modulus has been frequently employed to investigate the onset of the topological order characterizing the low-temperature phase of the two-dimensional XY model. We here present how the analysis based on the use of this key quantity can be applied to the study of the properties of coupled layers. To this aim, we first discuss how to extend the popular worm algorithm to a layered sample, and in particular to the evaluation of the longitudinal helicity, that we introduce taking care of the fact that the virtual twist representing the elastic deformation one applies to properly define the helicity modulus can act on a single layer or on all of them. We then apply the method to investigate the bilayer XY model, showing how the helicity modulus can be used to determine the phase diagram of the model as a function of temperature and interlayer coupling strength.

arXiv:2407.11507•Physical review. B./Physical review. B•Mar 11, 2025•2 citations
We study the response of a two-dimensional classical XY model to a finite (noninfinitesimal) twist of the boundary conditions. We use Monte Carlo simulations to evaluate the free energy difference between periodic and twisted-periodic boundary conditions and find deviations from the expected quadratic dependence on the twist angle. Consequently, the helicity modulus (spin stiffness) shows a nontrivial dependence on the twist angle. We show that the deviation from the expected behavior arises because of the mixing of states with opposite chirality which leads to an additional entropy contribution in the quasi-long-range ordered phase. We give an improved prescription for the numerical evaluation of the helicity modulus for a finite twist, and we discuss the spontaneous breaking of the chiral symmetry for the antiperiodic boundary conditions. We also discuss applications to discrete spin systems and some experimental scenarios where boundary conditions with finite twist are necessary.
A new logarithmic term in the finite-size helicity modulus has been predicted recently for the three-dimensional classical XY-model universality class under nonperiodic boundary conditions. These theoretical ideas have important implications for the analysis of superfluid-density measurements in confined geometries. We present numerical results from Monte Carlo simulations of the helicity modulus for the three-dimensional classical XY model on a finite simple-cubic lattice, with pinned-surface-spin boundary conditions. Our data are consistent with the presence of an additional logarithmic term and an estimate for the amplitude is given.
The critical scaling of the helicity modulus of the classical O(3) 3d3d3d Heisenberg ferromagnet is studied directly. Monte Carlo methods that impose either an antiperiodic boundary condition or a finite twist of definite handedness across otherwise periodic boundaries in one lattice direction are used to measure scale-dependent enthalpy variations in a simple cubic lattice at the ferromagnetic critical temperature. Finite-size scaling is then used to determine the critical exponents vE{v}_{E}vE​ and vF{v}_{F}vF​ for helicity and, by evaluating three independent hyperscaling-linked pairs of $\ensuremath{\nu}$ and $\ensuremath{\alpha}$, to test hyperscaling for this model. It is observed that antiperiodic boundary conditions in particular constrain the lattice to have a nonzero topological charge, establishing a connection between topological charge and helicity in the model.
We study response to a twist in the two-dimensional ppp-state clock model, which has the discrete Zp{Z}_{p}Zp​ symmetry. The response is measured in terms of helicity modulus, which is usually defined with respect to an infinitesimal twist. However, we demonstrate that such a definition is inappropriate for the clock model. The helicity modulus must be defined with respect to a finite, quantized twist which matches the discrete Zp{Z}_{p}Zp​ symmetry of the model. Numerical study of the appropriately defined helicity modulus resolves controversy over the clock model, showing the existence of two Berezinskii-Kosterlitz-Thouless transitions for p>4p>4p>4.
This is an analytic study of the two-dimensional XY spin glass with \ifmmode\pm\else\textpm\fi{}J disorder. The Hamiltonian has a continuous spin symmetry and a discrete chiral symmetry, and therefore possesses, potentially, two different order parameters and correlation lengths. The cost of breaking the symmetries is probed by comparing the ground-state energy under periodic (P) boundary conditions with the one under antiperiodic (AP) and under reflecting (R) boundary conditions. Two energy differences (domain-wall energies) appear, \ensuremath{\Delta}EAP{\mathit{E}}^{\mathrm{AP}}EAP and \ensuremath{\Delta}ER{\mathit{E}}^{\mathit{R}}ER, whose scaling behavior with system size is nontrivially related to the correlation length exponents. For a specific distribution of the \ifmmode\pm\else\textpm\fi{}J disorder we show that the chiral, and spin correlation lengths diverge with the same exponent as T goes down to 0. The common exponent has a common cause, viz. the low reversal energy of domains of chiral variables. For general disorder we give a heuristic argumentation in terms of droplet excitations that leads to spin ordering on a longer, or equal, scale than chiral ordering. These results are in contrast with intepretations of Monte Carlo simulations.
Interest in finite-size systems has risen in the last decades, due to the focus on nanotechnological applications and because they are convenient for numerical treatment that can subsequently be extrapolated to infinite lattices.Independently of the envisioned application, special attention must be given to boundary condition, which may or may not preserve the symmetry of the infinite lattice.Here we present a detailed study of the compatibility between boundary conditions and conservation laws.The conflict between open boundary conditions and momentum conservation is well understood, but we examine other symmetries, as well: we discuss gauge invariance, inversion, spin, and particle-hole symmetry and their compatibility with open, periodic, and twisted boundary conditions.In the interest of clarity, we develop the reasoning in the framework of the one-dimensional half-filled Hubbard model, whose Hamiltonian displays a variety of symmetries.Our discussion includes analytical and numerical results.Our analytical survey shows that, as a rule, boundary conditions break one or more symmetries of the infinite-lattice Hamiltonian.The exception is twisted boundary condition with the special torsion Θ = πL/2, where L is the lattice size.Our numerical results for the ground-state energy at half-filling and the energy gap for L = 2-7 show how the breaking of symmetry affects the convergence to the L → ∞ limit.We compare the computed energies and gaps with the exact results for the infinite lattice drawn from the Bethe-Ansatz solution.The deviations are boundary-condition dependent.The special torsion yields more rapid convergence than open or periodic boundary conditions.For sizes as small as L = 7, the numerical results for twisted condition are very close to the L → ∞ limit.We also discuss the ground-state electronic density and magnetization at half filling under the three boundary conditions.
The behavior of the helicity modulus has been frequently employed to investigate the onset of the topological order characterizing the low-temperature phase of the two-dimensional XY model. We here present how the analysis based on the use of this key quantity can be applied to the study of the properties of coupled layers. To this aim, we first discuss how to extend the popular worm algorithm to a layered sample, and in particular to the evaluation of the longitudinal helicity, that we introduce taking care of the fact that the virtual twist representing the elastic deformation one applies to properly define the helicity modulus can act on a single layer or on all of them. We then apply the method to investigate the bilayer XY model, showing how the helicity modulus can be used to determine the phase diagram of the model as a function of temperature and interlayer coupling strength.