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

Approaching the Kosterlitz-Thouless transition for the classical XY model with tensor networks

Laurens Vanderstraeten, Bram Vanhecke, Andreas Martin Läuchli, Frank Verstraete

Published Dec 26, 2019• 39 citations
arXivHTMLPDF

Abstract

We apply variational tensor-network methods for simulating the Kosterlitz-Thouless phase transition in the classical two-dimensional XY model. In particular, using uniform matrix product states (MPS) with non-Abelian O(2) symmetry, we compute the universal drop in the spin stiffness at the critical point. In the critical low-temperature regime, we focus on the MPS entanglement spectrum to characterize the Luttinger-liquid phase. In the high-temperature phase, we confirm the exponential divergence of the correlation length and estimate the critical temperature with high precision. Our MPS approach can be used to study generic two-dimensional phase transitions with continuous symmetries.

DOI
10.1103/physreve.100.062136

Related papers

Hybrid Berezinskii-Kosterlitz-Thouless and Ising topological phase transition in the generalized two-dimensional XY model using tensor networks

Feng-Feng Song, Guang-Ming Zhang

In tensor network representation, the partition function of a generalized two-dimensional XY spin model with topological integer and half-integer vortex excitations is mapped to a tensor product of one-dimensional quantum transfer operator, whose eigenequation can be solved by an algorithm of variational uniform matrix product states. Using the singularities of the entanglement entropy, we accurately determine the complete phase diagram of this model. Both the integer vortex-antivortex binding and half-integer vortex-antivortex binding phases are separated from the disordered phase by the usual Berezinskii-Kosterlitz-Thouless (BKT) transitions, while a continuous topological phase transition exists between two different vortex binding phases, exhibiting a logarithmic divergence of the specific heat and exponential divergence of the spin correlation length. A hybrid BKT and Ising universality class of topological phase transition is thus established. We further prove that three phase transition lines meet at a multicritical point from which a deconfinement crossover line extends into the disordered phase.

arXiv:2011.05682•Physical review. B./Physical review. B•Jan 19, 2021•13 citations

Critical analysis of two-dimensional classical XY model

Raghav G Jha

Abstract We consider the two-dimensional classical XY model on a square lattice in the thermodynamic limit using tensor renormalization group and precisely determine the critical temperature corresponding to the Berezinskii–Kosterlitz–Thouless (BKT) phase transition to be 0.89290(5) which is an improvement compared to earlier studies using tensor network methods.

arXiv:2004.06314•Journal of Statistical Mechanics Theory and Experiment•Aug 1, 2020•22 citations

Contrasting pseudocriticality in the classical two-dimensional Heisenberg and RP2 models: Zero-temperature phase transition versus finite-temperature crossover

Lander Burgelman, Lukas Devos, Bram Vanhecke +2 more

Tensor-network methods are used to perform a comparative study of the two-dimensional classical Heisenberg and RP^{2} models. We demonstrate that uniform matrix product states (MPSs) with explicit SO(3) symmetry can probe correlation lengths up to O(10^{3}) sites accurately, and we study the scaling of entanglement entropy and universal features of MPS entanglement spectra. For the Heisenberg model, we find no signs of a finite-temperature phase transition, supporting the scenario of asymptotic freedom. For the RP^{2} model we observe an abrupt onset of scaling behavior, consistent with hints of a finite-temperature phase transition reported in previous studies. A careful analysis of the softening of the correlation length divergence, the scaling of the entanglement entropy, and the MPS entanglement spectra shows that our results are inconsistent with true criticality, but are rather in agreement with the scenario of a crossover to a pseudocritical region which exhibits strong signatures of nematic quasi-long-range order at length scales below the true correlation length. Our results reveal a fundamental difference in scaling behavior between the Heisenberg and RP^{2} models: Whereas the emergence of scaling in the former shifts to zero temperature if the bond dimension is increased, it occurs at a finite bond-dimension independent crossover temperature in the latter.

arXiv:2202.07597•Physical review. E•Jan 13, 2023•7 citations

Tensor network approach to the two-dimensional fully frustrated XY model and a chiral ordered phase

Feng-Feng Song, Guang-Ming Zhang

A general framework is proposed to solve the two-dimensional (2D) fully frustrated XY model for the Josephson junction arrays in a perpendicular magnetic field. The essential idea is to encode the ground-state local rules induced by frustrations in the local tensors of the partition function. The partition function is then expressed in terms of a product of a 1D transfer matrix operator, whose eigenequation can be solved by an algorithm of matrix product states rigorously. The singularity of the entanglement entropy for the 1D quantum analog provides a stringent criterion to distinguish various phase transitions without identifying a priori order parameter. Two very close phase transitions are determined at ${T}_{c1}\ensuremath{\approx}0.4459$ and ${T}_{c2}\ensuremath{\approx}0.4532$, respectively. The former corresponds to a Berezinskii-Kosterlitz-Thouless phase transition describing the phase coherence of XY spins, and the latter is an Ising-like continuous phase transition below which a chirality order with spontaneously broken ${Z}_{2}$ symmetry is established.

arXiv:2112.03550•Physical review. B./Physical review. B•Apr 20, 2022•22 citations

Tensor network approach to the fully frustrated XY model on a kagome lattice with a fractional vortex-antivortex pairing transition

Feng-Feng Song, Guang-Ming Zhang

We have developed a tensor network approach to the two-dimensional fully frustrated classical XY spin model on the kagome lattice, and clarified the nature of the possible phase transitions of various topological excitations. We find that the standard tensor network representation for the partition function does not work due to the strong frustrations in the low-temperature limit. To avoid the direct truncation of the Boltzmann weight, based on the duality transformation, we introduce a different representation to build the tensor network with local tensors lying on the centers of the elementary triangles of the kagome lattice. Then the partition function is expressed as a product of one-dimensional transfer matrix operators, whose eigenequation can be solved by the variational uniform matrix product state algorithm accurately. The singularity of the entanglement entropy for the one-dimensional quantum operator provides a stringent criterion for the possible phase transitions. Through a systematic numerical analysis of thermodynamic properties and correlation functions in the thermodynamic limit, we prove that the model exhibits a single Berezinskii-Kosterlitz-Thouless phase transition only, which is driven by the unbinding of $\frac{1}{3}$ fractional vortex-antivortex pairs determined at ${T}_{c}\ensuremath{\simeq}0.075{J}_{1}$ accurately. The absence of long-range order of chirality or quasi-long-range order of integer vortices has been verified in the whole finite-temperature range. Thus, the long-standing controversy about the phase transitions in this fully frustrated XY model on the kagome lattice is solved rigorously, which provides a plausible way to understand the charge-$6e$ superconducting phase observed experimentally in the two-dimensional kagome superconductors.

arXiv:2304.11944•Physical review. B./Physical review. B•Jul 20, 2023•12 citations

Berezinskii–Kosterlitz–Thouless Transition – A Universal Neural Network Study with Benchmarks

Y.-H. Tseng, Fu-Jiun Jiang

Using a supervised neural network (NN) trained once on a one-dimensional lattice of 200 sites, we calculate the Berezinskii–Kosterlitz–Thouless phase transitions of the two-dimensional (2D) classical XY and the 2D generalized classical XY models. In particular, both the bulk quantities Binder ratios and the spin states of the studied systems are employed to construct the needed configurations for the NN prediction. By applying semiempirical finite-size scaling to the relevant data, the critical points obtained by the NN approach agree well with the known results established in the literature. This implies that for each of the considered models, the determination of its various phases requires only a little information. The outcomes presented here demonstrate convincingly that the employed universal NN is not only valid for the symmetry breaking related phase transitions, but also works for calculating the critical points of the phase transitions associated with topology. The efficiency of the used NN in the computation is examined by carrying out several detailed benchmark calculations.

Results in Physics•Jan 1, 2021•7 citations
We apply variational tensor-network methods for simulating the Kosterlitz-Thouless phase transition in the classical two-dimensional XY model. In particular, using uniform matrix product states (MPS) with non-Abelian O(2) symmetry, we compute the universal drop in the spin stiffness at the critical point. In the critical low-temperature regime, we focus on the MPS entanglement spectrum to characterize the Luttinger-liquid phase. In the high-temperature phase, we confirm the exponential divergence of the correlation length and estimate the critical temperature with high precision. Our MPS approach can be used to study generic two-dimensional phase transitions with continuous symmetries.
In tensor network representation, the partition function of a generalized two-dimensional XY spin model with topological integer and half-integer vortex excitations is mapped to a tensor product of one-dimensional quantum transfer operator, whose eigenequation can be solved by an algorithm of variational uniform matrix product states. Using the singularities of the entanglement entropy, we accurately determine the complete phase diagram of this model. Both the integer vortex-antivortex binding and half-integer vortex-antivortex binding phases are separated from the disordered phase by the usual Berezinskii-Kosterlitz-Thouless (BKT) transitions, while a continuous topological phase transition exists between two different vortex binding phases, exhibiting a logarithmic divergence of the specific heat and exponential divergence of the spin correlation length. A hybrid BKT and Ising universality class of topological phase transition is thus established. We further prove that three phase transition lines meet at a multicritical point from which a deconfinement crossover line extends into the disordered phase.
Tensor-network methods are used to perform a comparative study of the two-dimensional classical Heisenberg and RP^{2} models. We demonstrate that uniform matrix product states (MPSs) with explicit SO(3) symmetry can probe correlation lengths up to O(10^{3}) sites accurately, and we study the scaling of entanglement entropy and universal features of MPS entanglement spectra. For the Heisenberg model, we find no signs of a finite-temperature phase transition, supporting the scenario of asymptotic freedom. For the RP^{2} model we observe an abrupt onset of scaling behavior, consistent with hints of a finite-temperature phase transition reported in previous studies. A careful analysis of the softening of the correlation length divergence, the scaling of the entanglement entropy, and the MPS entanglement spectra shows that our results are inconsistent with true criticality, but are rather in agreement with the scenario of a crossover to a pseudocritical region which exhibits strong signatures of nematic quasi-long-range order at length scales below the true correlation length. Our results reveal a fundamental difference in scaling behavior between the Heisenberg and RP^{2} models: Whereas the emergence of scaling in the former shifts to zero temperature if the bond dimension is increased, it occurs at a finite bond-dimension independent crossover temperature in the latter.
A general framework is proposed to solve the two-dimensional (2D) fully frustrated XY model for the Josephson junction arrays in a perpendicular magnetic field. The essential idea is to encode the ground-state local rules induced by frustrations in the local tensors of the partition function. The partition function is then expressed in terms of a product of a 1D transfer matrix operator, whose eigenequation can be solved by an algorithm of matrix product states rigorously. The singularity of the entanglement entropy for the 1D quantum analog provides a stringent criterion to distinguish various phase transitions without identifying a priori order parameter. Two very close phase transitions are determined at ${T}_{c1}\ensuremath{\approx}0.4459$ and ${T}_{c2}\ensuremath{\approx}0.4532$, respectively. The former corresponds to a Berezinskii-Kosterlitz-Thouless phase transition describing the phase coherence of XY spins, and the latter is an Ising-like continuous phase transition below which a chirality order with spontaneously broken Z2{Z}_{2}Z2​ symmetry is established.
We have developed a tensor network approach to the two-dimensional fully frustrated classical XY spin model on the kagome lattice, and clarified the nature of the possible phase transitions of various topological excitations. We find that the standard tensor network representation for the partition function does not work due to the strong frustrations in the low-temperature limit. To avoid the direct truncation of the Boltzmann weight, based on the duality transformation, we introduce a different representation to build the tensor network with local tensors lying on the centers of the elementary triangles of the kagome lattice. Then the partition function is expressed as a product of one-dimensional transfer matrix operators, whose eigenequation can be solved by the variational uniform matrix product state algorithm accurately. The singularity of the entanglement entropy for the one-dimensional quantum operator provides a stringent criterion for the possible phase transitions. Through a systematic numerical analysis of thermodynamic properties and correlation functions in the thermodynamic limit, we prove that the model exhibits a single Berezinskii-Kosterlitz-Thouless phase transition only, which is driven by the unbinding of 13\frac{1}{3}31​ fractional vortex-antivortex pairs determined at ${T}_{c}\ensuremath{\simeq}0.075{J}_{1}$ accurately. The absence of long-range order of chirality or quasi-long-range order of integer vortices has been verified in the whole finite-temperature range. Thus, the long-standing controversy about the phase transitions in this fully frustrated XY model on the kagome lattice is solved rigorously, which provides a plausible way to understand the charge-6e6e6e superconducting phase observed experimentally in the two-dimensional kagome superconductors.
Using a supervised neural network (NN) trained once on a one-dimensional lattice of 200 sites, we calculate the Berezinskii–Kosterlitz–Thouless phase transitions of the two-dimensional (2D) classical XY and the 2D generalized classical XY models. In particular, both the bulk quantities Binder ratios and the spin states of the studied systems are employed to construct the needed configurations for the NN prediction. By applying semiempirical finite-size scaling to the relevant data, the critical points obtained by the NN approach agree well with the known results established in the literature. This implies that for each of the considered models, the determination of its various phases requires only a little information. The outcomes presented here demonstrate convincingly that the employed universal NN is not only valid for the symmetry breaking related phase transitions, but also works for calculating the critical points of the phase transitions associated with topology. The efficiency of the used NN in the computation is examined by carrying out several detailed benchmark calculations.