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

Physical review. B, Condensed matter

Anisotropic three-dimensionalXYmodel and vortex-loop scaling

Subodh R. Shenoy, Biplab Chattopadhyay

Published Apr 1, 1995• 57 citations
Publisher

Abstract

Scaling equations are obtained for thermally excited vortex loops in the anisotropic three-dimensional (3D) XY model with interplane/intraplane coupling ratio ${\mathit{K}}_{\mathrm{\ensuremath{\perp}}}$/${\mathit{K}}_{\mathrm{\ensuremath{\parallel}}}$=${\ensuremath{\gamma}}_{0}^{\mathrm{\ensuremath{-}}2}$. For high-${\mathit{T}}_{\mathit{c}}$ superconductors, ${\ensuremath{\gamma}}_{0}^{\mathrm{\ensuremath{-}}2}$, related to the Ginzburg-Landau masses, is in a strong anisotropy regime, ${\ensuremath{\gamma}}_{0}^{\mathrm{\ensuremath{-}}1}$0.5, and a length scale ${\mathit{r}}_{0}$=${\ensuremath{\gamma}}_{0}$${\mathit{a}}_{0}$ naturally arises (${\mathit{a}}_{0}$= lattice constant). (i) For loop major axes a>${\mathit{r}}_{0}$, the dominant excitations are 3D elliptical vortex loops cutting multiple planes, with a\ensuremath{\sim}${\ensuremath{\xi}}_{\mathrm{\ensuremath{-}}}$(T)=${\mathit{a}}_{0}$\ensuremath{\Vert}\ensuremath{\varepsilon}${\mathrm{\ensuremath{\Vert}}}^{\mathrm{\ensuremath{-}}\ensuremath{\nu}}$ where \ensuremath{\Vert}\ensuremath{\varepsilon}\ensuremath{\Vert}\ensuremath{\equiv}\ensuremath{\Vert}(T-${\mathit{T}}_{\mathit{c}}$)/${\mathit{T}}_{\mathit{c}}$\ensuremath{\Vert}, and \ensuremath{\nu}=0.67. The vorticity segment components (\ensuremath{\mu}=x^,y^,z^) are ${\mathit{J}}_{\mathrm{\ensuremath{\mu}}}$(r)=0,\ifmmode\pm\else\textpm\fi{}1, and interact via a Biot-Savart-like law. The renormalized anisotropy ${\ensuremath{\gamma}}_{\mathit{l}}^{\mathrm{\ensuremath{-}}1}$\ensuremath{\rightarrow}1, asymptotically isotropic, as l=ln(a/${\mathit{a}}_{0}$)\ensuremath{\rightarrow}\ensuremath{\infty}: anisotropy is irrelevant. (ii) For loop scales ${\mathit{r}}_{0}$>a>2${\mathit{a}}_{0}$ the dominant excitations are quasi-2D rectangular loops with short sides ${\mathit{J}}_{\mathit{z}}$(r)=\ifmmode\pm\else\textpm\fi{}1 cutting single planes, that are thus effectively decoupled at finite scales ${\mathit{r}}_{0}$. The ${\mathit{J}}_{\mathit{z}}$=\ifmmode\pm\else\textpm\fi{}1 vortex components, of in-plane separation \ensuremath{\Vert}R${\ensuremath{\rightarrow}}_{\mathrm{W}}$mN\ensuremath{\Vert}${\mathit{r}}_{0}$, interact via a logarithmic [ln(\ensuremath{\Vert}R${\ensuremath{\rightarrow}}_{\mathrm{\ensuremath{\parallel}}}$\ensuremath{\Vert}/${\mathit{r}}_{0}$)] plus linear [${\mathit{Q}}_{0}$(\ensuremath{\Vert}R${\ensuremath{\rightarrow}}_{\mathrm{\ensuremath{\parallel}}}$\ensuremath{\Vert}/${\mathit{a}}_{0}$-1)] potential. As ${\ensuremath{\gamma}}_{0}^{\mathrm{\ensuremath{-}}1}$\ensuremath{\rightarrow}0, the coefficient ${\mathit{Q}}_{0}$\ensuremath{\rightarrow}0, and the Kosterlitz-Thouless limit is recovered. (iii)The 3D transition temperature ${\mathit{T}}_{\mathit{c}}$ versus ${\mathit{K}}_{\mathrm{\ensuremath{\perp}}}$/${\mathit{K}}_{\mathrm{\ensuremath{\parallel}}}$, calculated for strong anisotropies from ``2D'' (3D) scaling for scales ${\mathit{r}}_{0}$ (>${\mathit{r}}_{0}$), matches existing Monte Carlo data well. Critical regions \ensuremath{\Vert}${\mathrm{\ensuremath{\varepsilon}}}_{\mathit{c}}$\ensuremath{\Vert}\ensuremath{\le}${\ensuremath{\gamma}}_{0}^{\mathrm{\ensuremath{-}}1/\ensuremath{\nu}}$ are estimated. Contact is made with ideas of an intrinsic critical current arising from linear \ensuremath{\sim}${\ensuremath{\gamma}}_{0}^{\mathrm{\ensuremath{-}}1}$${\mathit{R}}_{\mathrm{\ensuremath{\parallel}}}$ vortex segment effective potentials.

DOI
10.1103/physrevb.51.9129

Related papers

Monte Carlo analysis of the two-dimensionalXYmodel. II. Comparison with the Kosterlitz renormalization-group equations

Peter Olsson

The temperature scale of relevance for the vortex excitations in the two-dimensional XY model is identified. The size dependence of the helicity modulus at constant vortex temperature is then obtained by Monte Carlo simulations. The Monte Carlo data are then compared with the results for the length-dependent dielectric function from the Kosterlitz renormalization-group equations. The agreement is excellent for temperatures up to and slightly above ${\mathit{T}}_{\mathit{c}}$. The critical temperature is determined to be ${\mathit{T}}_{\mathit{c}}$=0.892 13(10). The temperature dependence of the characteristic length right above ${\mathit{T}}_{\mathit{c}}$ is also obtained and found to be in agreement with recent results from high-temperature series expansions.

Physical review. B, Condensed matter•Aug 1, 1995•110 citations

Critical behavior of the two-dimensionalXYmodel

Rajan Gupta, CLIVE F. BAILLIE

We present detailed Monte Carlo results for the susceptibility \ensuremath{\chi}, correlation length \ensuremath{\xi}, and specific heat ${\mathit{C}}_{\mathit{v}}$ for the XY model. The simulations are done on ${64}^{2}$, ${128}^{2}$, ${256}^{2}$, and ${512}^{2}$ lattices over the temperature range 0.98\ensuremath{\le}T\ensuremath{\le}1.43 corresponding to 270. Fits to \ensuremath{\chi} and \ensuremath{\xi} data favor a Kosterlitz-Thouless (KT) singularity over a second-order transition; however, unconstrained four-parameter KT fits do not confirm the predicted values \ensuremath{\nu}=0.5 and \ensuremath{\eta}=0.25. Our best estimate, ${\mathit{T}}_{\mathit{c}}$=0.894(5), is obtained using KT fits with \ensuremath{\nu} fixed at 0.5. The exponent \ensuremath{\eta} is calculated as a function of temperature in the spin-wave phase using Monte Carlo renormalization-group and finite-size-scaling methods. Both methods give consistent results and we find \ensuremath{\eta}\ensuremath{\approxeq}0.235 at T=0.894. We also present results for the behavior of vortex density across the transition and exhibit how the dilute-gas approximation breaks down.

Physical review. B, Condensed matter•Feb 1, 1992•196 citations

A Monte Carlo Study of the Criticality of the 3D XY Model

Sindre Vaskinn

In this master thesis the criticality of the 3D XY model with cubic anisotropy has been examined by the means of Monte Carlo simulation. The 3D XY model is a special case of the $O(n)$ spin-vector model which describe numerous phenomena's. Here, $n$ is the number of vector dimensions and for the XY model $n=2$. The $O(n)$-model is invariant for global rotations. This means that it is arbitrary which direction in the symmetry breaking occurs. In real system one should always expect perturbations from a perfect system, here it is done by adding a cubic anisotropy (order $p=4$) then the symmetry of the Hamiltonian is lowered and it will be advantageous for the spins to be ordered in directions $90^o$ apart. This symmetry is called $C_{4}$, which is a part of the cyclic group, and it restricts the symmetry breaking to one of these four directions. \n\nThe criticality of this model is when the spins suddenly align as the temperature is lowered to a critical temperature. Sufficiently close to the critical temperature quantities like the magnetization, which is the sum of spins, will show a power law growth characterised by critical exponents specified by the class of universality of the transition. The critical exponents can be obtained by finite-size scaling. By applying renormalization theory on the critical region one can show that some of parameters in the Hamiltonian are: relevant and effect the scaling, irrelevant and does not effect the scaling and dangerously irrelevant meaning that they only affect some of the critical exponents. Previous Monte Carlo and renormalization calculus show that anisotropy which is relevant for $p<4$ becomes dangerously irrelevant for $p\\geq4$. The question of interest is what kind of effect the dangerously irrelevant parameter in the Hamiltonian has on the scaling. This is interesting because the effect of dangerously irrelevant parameters is hardly studied and anisotropy arise in many crystal systems because of the lattice geometry. \n\nThe results are obtained by Monte Carlo simulations and the critical slowing down close to the phase transition is solved by the method of parallel tempering. This method was chosen because it is simple to implement in a parallel code.\n\nThe Binder cumulant for different quantities and the peak of the magnetic susceptibility was used support a critical temperature $T_c = 2.2019 \\ \\textnormal k$ for no field applied and no anisotropy. The known critical exponents for the no field and no anisotropy was used. The scaling made a good collapse, meaning that the right $T_c$ and the correct critical exponents where used. The scaling did not take into account a magnetic field and therefore a good collapse of curves where not obtained in this case. When the cubic anisotropy was added the critical exponent for the magnetic susceptibility and the cubic order parameter changed, while the other exponents where unchanged. This supports the view that cubic anisotropy is dangerously irrelevant for the 3D XY model. The scaling of the magnetic susceptibility showed a possible new kind of scaling, a scaling which is different on the two sides of the transition. Although a more accurate determination of the critical exponent in the ordered phase require more extensive simulations. It appears that this exponent is larger than the exponent in the disordered phase, where one recovered the isotropic exponent. The difference in the critical exponents has been predicted to increase with $p$ and an investigation of the hexagonal anisotropy $p=6$ has been started.

52•Jan 1, 2015•0 citations

Finite-size scaling study of the vortex-free three-dimensionalXYmodel

Ronald Fisch

Monte Carlo simulations have been used to study a version of the XY model without vortex loops for J/T=0 on simple cubic and face-centered-cubic lattices. Finite-size scaling of data from L\ifmmode\times\else\texttimes\fi{}L\ifmmode\times\else\texttimes\fi{}L lattices with L up to 32 is used to obtain values for the order parameter \ensuremath{\Vert}M\ensuremath{\Vert} and a 1/L finite-size correction on each lattice. The order-parameter susceptibility is found to be diverging as ${\mathit{L}}^{0.8}$. We then study the crossover from three-dimensional to one-dimensional behavior on L\ifmmode\times\else\texttimes\fi{}W\ifmmode\times\else\texttimes\fi{}W simple cubic lattices, with W=4 and 8, and L up to 1024.

Physical review. B, Condensed matter•Nov 1, 1995•9 citations

Crossover and Critical Behavior in the Layered XY Model

Roman Kracht, Andrea Trombettoni, Ilaria Maccari +1 more

Motivated by the interplay between 2D and 3D scaling signatures observed in unconventional layered superconductors, we present a systematic Monte Carlo study of the three-dimensional classical XY model with anisotropic in-plane $J_\parallel$ and inter-plane $J_\perp$ couplings. Our study includes very small values of the system anisotropy $Δ=J_\perp /J_\parallel$ not studied before, and focuses on characterizing the crossover from quasi-2D topological scaling to genuine 3D critical behavior. The numerical results for the critical temperature unambiguously reveal a logarithmic scaling with $Δ$, directly related to the topological scaling in the 2D limit. Despite the 3D nature of the layered XY criticality, topological scaling signatures survive up to system sizes comparable to the crossover length $\ell_J$, which diverges at small $Δ$ with a scaling behavior reminiscent of the Berezinskii-Kosterlitz-Thouless (BKT) transition. This shows that genuine 3D symmetry-breaking behavior emerges only at exceedingly large system sizes when the anisotropy is very strong. Our results indicate that new experimental evidence is required to clarify the extent to which the critical signatures observed in layered strongly correlated materials are shaped by their pronounced anisotropy.

arXiv:2603.19351•Mar 19, 2026•0 citations

Criticality and novel quantum liquid phases in Ginzburg--Landau theories with compact and non-compact gauge fields

J. Smiseth

We have studied the critical properties of three-dimensional U(1)-symmetric lattice gauge theories. The models apply to various physical systems such as insulating phases of strongly correlated electron systems as well as superconducting and superfluid states of liquid metallic hydrogen under extreme pressures. This thesis contains an introductory part and a collection of research papers of which seven are published works and one is submitted for publication.\n\nPaper I: Critical properties of the 2+1-dimensional compact abelian Higgs model with gauge charge q=2 are studied. We introduce a novel method of computing the third moment M3 of the action which allows us to extract correlation length and specific heat critical exponents ν and α without invoking hyperscaling. Finite-size scaling analysis of M3 yields the ratio (1+α)/ν and 1/ν separately. We find that α and ν vary along the critical line of the theory, which however exhibits a remarkable resilience of Z2 criticality. We conclude that the model is a fixed-line theory, which we propose to characterize the zero temperature quantum phase transition from a Mott-Hubbard insulator to a charge fractionalized insulator in two spatial dimensions.\n\nPaper II: Large scale Monte Carlo simulations are employed to study phase transitions in the three-dimensional compact abelian Higgs model in adjoint representations of the matter field, labeled by an integer q, for q=2,3,4,5. We also study various limiting cases of the model, such as the Zq lattice gauge theory, dual to the 3DZq spin model, and the 3D xy spin model which is dual to the Zq lattice gauge theory in the limit q → ∞. In addition, for benchmark purposes, we study the 2D square lattice 8-vertex model, which is exactly solvable and features non-universal critical exponents. The critical exponents α and ν are calculated from finite size scaling of the third moment of the action, and the method is tested thoroughly on models with known values for these exponents. We have found that for q=3, the three-dimensional compact abelian Higgs model exhibits a second order phase transition line which joins a first order phase transition line at a tricritical point. The results for q=2 in Paper I are reported with a higher lever of detail.\n\nPaper III: This paper is based on a talk by F. S. Nogueira in the Aachen HEP 2003 conference where a review of the results for the compact abelian Higgs model from Paper I and Paper II was presented, as well as the results for the q=1 case studied by F. S. Nogueira, H. Kleinert and A. Sudbø.\n\nPaper IV: We study the effects of a Chern-Simons (CS) term in the phase structure of two different abelian gauge theories in three dimensions. By duality transformations we show how the compact U(1) gauge theory with a CS term for certain values of the CS coupling can be written as a gas of vortex loops interacting through steric repulsion. This theory is known to exhibit a phase transition governed by proliferation of vortex loops. We also employ Monte Carlo simulations to study the non-compact U(1) abelian Higgs model with a CS term. Finite size scaling of the third moment of the action yields critical exponents α and ν that vary continuously with the strength of the CS term, and a comparison with available analytical results is made.\n\nPaper V: The critical properties of N-component Ginzburg-Landau theory are studied in d=2+1 dimensions. The model is dualized to a theory of N vortex fields interacting through a Coulomb and a screened potential. The model with N=2 shows two anomalies in the specific heat. From Monte Carlo simulations we calculate the critical exponents α and ν and the mass of the gauge field. We conclude that one anomaly corresponds to an inverted 3D xy fixed point, while the other corresponds to a 3D xy fixed point. There are N fixed points, namely one corresponding to an inverted 3D xy fixed point, and N-1corresponding to neutral 3D xy fixed points. Applications are briefly discussed.\n\nPaper VI: The phase diagram and critical properties of the N-component London superconductor are studied both analytically and through large-scale Monte-Carlo simulations in d=2+1 dimensions. The model with different bare phase stiffnesses for each flavor is a model of superconductivity which should arise out of metallic phases of light atoms under extreme pressure. A projected mixture of electronic and protonic condensates in liquid metallic hydrogen under extreme pressure is the simplest example, corresponding to N=2 with individually conserved matter fields. We compute critical exponents α and ν for N=2 and N=3. The results from Paper V are presented at a higher level of detail. For the arbitrary N case, there are N fixed points,namely one charged inverted 3D xy fixed point, and N-1 neutral 3D xy fixed points. We explicitly identify one charged vortex mode and N-1 neutral vortex modes. The model for N=2 and equal bare phase stiffnesses corresponds to a field theoretical description of an easy-plane quantum antiferromagnet. In this case, the critical exponents are computed and found to be non 3D xy values. Furthermore, we study the model in an external magnetic field, and find a novel feature, namely N-1 superfluid phases arising out of N charged condensates. In particular, for N=2 we point out the possibility of two novel types of field-induced phase transitions in ordered quantum fluids: i) A phase transition from a superconductor to a superfluid or vice versa, driven by tuning an external magnetic field. This identifies the superconducting phase of liquid metallic hydrogen as a novel quantum fluid. ii) A phase transition corresponding to a quantum fluid analogue of sublattice melting, where a composite field-induced Abrikosov vortex lattice is decomposed and disorders the phases of the constituent condensate with lowest bare phase stiffness. Both transitions belong to the 3D xy universality class.\n\nPaper VII: We consider the vortex superconductor with two individually conserved condensates in a finite magnetic field. The ground state is a lattice of cocentered vortices in both order parameters. We find two novel phase transitions when temperature is increased at fixed magnetic field. i) A "vortex sublattice melting" transition where vortices in the field with lowest phase stiffness ("light vortices") loose cocentricity with the vortices with large phase stiffness ("heavy vortices"), entering a liquid state (the structure factor of the light vortex sublattice vanishes continuously.) This transition is in the 3D xy universality class. ii) A first order melting transition of the lattice of heavy vortices in a liquid of light vortices.\n\nPaper VIII: We report on large-scale Monte Carlo simulations of a novel type of a vortex matter phase transition which should take place in a three dimensional two-component superconductor. We identify the regime where first, at a certain temperature a field-induced lattice of co-centered vortices of both order parameters melts, causing the system to loose superconductivity. In this state the two-gap system retains a broken composite symmetry and we observe that at a higher temperature it undergoes an extra phase transition where the disordered composite one-flux-quantum vortex lines are "ionized" into a "plasma" of constituent fractional flux vortex lines in individual order parameters. This is the hallmark of the superconductor-to-superfluid-to-normal fluid phase transitions projected to occur in e.g. liquid metallic hydrogen.

•Jan 1, 2005•0 citations
Scaling equations are obtained for thermally excited vortex loops in the anisotropic three-dimensional (3D) XY model with interplane/intraplane coupling ratio K\ensuremath⊥{\mathit{K}}_{\mathrm{\ensuremath{\perp}}}K\ensuremath⊥​/K\ensuremath∥{\mathit{K}}_{\mathrm{\ensuremath{\parallel}}}K\ensuremath∥​=\ensuremathγ0\ensuremath−2{\ensuremath{\gamma}}_{0}^{\mathrm{\ensuremath{-}}2}\ensuremathγ0\ensuremath−2​. For high-Tc{\mathit{T}}_{\mathit{c}}Tc​ superconductors, \ensuremathγ0\ensuremath−2{\ensuremath{\gamma}}_{0}^{\mathrm{\ensuremath{-}}2}\ensuremathγ0\ensuremath−2​, related to the Ginzburg-Landau masses, is in a strong anisotropy regime, \ensuremathγ0\ensuremath−1{\ensuremath{\gamma}}_{0}^{\mathrm{\ensuremath{-}}1}\ensuremathγ0\ensuremath−1​0.5, and a length scale r0{\mathit{r}}_{0}r0​={\ensuremath{\gamma}}_{0}$${\mathit{a}}_{0} naturally arises (a0{\mathit{a}}_{0}a0​= lattice constant). (i) For loop major axes a>r0{\mathit{r}}_{0}r0​, the dominant excitations are 3D elliptical vortex loops cutting multiple planes, with a\ensuremath{\sim}\ensuremathξ\ensuremath−{\ensuremath{\xi}}_{\mathrm{\ensuremath{-}}}\ensuremathξ\ensuremath−​(T)=a0{\mathit{a}}_{0}a0​\ensuremath{\Vert}\ensuremath{\varepsilon}\ensuremath∥\ensuremath−\ensuremathν{\mathrm{\ensuremath{\Vert}}}^{\mathrm{\ensuremath{-}}\ensuremath{\nu}}\ensuremath∥\ensuremath−\ensuremathν where \ensuremath{\Vert}\ensuremath{\varepsilon}\ensuremath{\Vert}\ensuremath{\equiv}\ensuremath{\Vert}(T-Tc{\mathit{T}}_{\mathit{c}}Tc​)/Tc{\mathit{T}}_{\mathit{c}}Tc​\ensuremath{\Vert}, and \ensuremath{\nu}=0.67. The vorticity segment components (\ensuremath{\mu}=x^,y^,z^) are J\ensuremathμ{\mathit{J}}_{\mathrm{\ensuremath{\mu}}}J\ensuremathμ​(r)=0,\ifmmode\pm\else\textpm\fi{}1, and interact via a Biot-Savart-like law. The renormalized anisotropy \ensuremathγl\ensuremath−1{\ensuremath{\gamma}}_{\mathit{l}}^{\mathrm{\ensuremath{-}}1}\ensuremathγl\ensuremath−1​\ensuremath{\rightarrow}1, asymptotically isotropic, as l=ln(a/a0{\mathit{a}}_{0}a0​)\ensuremath{\rightarrow}\ensuremath{\infty}: anisotropy is irrelevant. (ii) For loop scales r0{\mathit{r}}_{0}r0​>a>2a0{\mathit{a}}_{0}a0​ the dominant excitations are quasi-2D rectangular loops with short sides Jz{\mathit{J}}_{\mathit{z}}Jz​(r)=\ifmmode\pm\else\textpm\fi{}1 cutting single planes, that are thus effectively decoupled at finite scales r0{\mathit{r}}_{0}r0​. The Jz{\mathit{J}}_{\mathit{z}}Jz​=\ifmmode\pm\else\textpm\fi{}1 vortex components, of in-plane separation \ensuremath{\Vert}R\ensuremath→W{\ensuremath{\rightarrow}}_{\mathrm{W}}\ensuremath→W​mN\ensuremath{\Vert}r0{\mathit{r}}_{0}r0​, interact via a logarithmic [ln(\ensuremath{\Vert}R\ensuremath→\ensuremath∥{\ensuremath{\rightarrow}}_{\mathrm{\ensuremath{\parallel}}}\ensuremath→\ensuremath∥​\ensuremath{\Vert}/r0{\mathit{r}}_{0}r0​)] plus linear [Q0{\mathit{Q}}_{0}Q0​(\ensuremath{\Vert}R\ensuremath→\ensuremath∥{\ensuremath{\rightarrow}}_{\mathrm{\ensuremath{\parallel}}}\ensuremath→\ensuremath∥​\ensuremath{\Vert}/a0{\mathit{a}}_{0}a0​-1)] potential. As \ensuremathγ0\ensuremath−1{\ensuremath{\gamma}}_{0}^{\mathrm{\ensuremath{-}}1}\ensuremathγ0\ensuremath−1​\ensuremath{\rightarrow}0, the coefficient Q0{\mathit{Q}}_{0}Q0​\ensuremath{\rightarrow}0, and the Kosterlitz-Thouless limit is recovered. (iii)The 3D transition temperature Tc{\mathit{T}}_{\mathit{c}}Tc​ versus K\ensuremath⊥{\mathit{K}}_{\mathrm{\ensuremath{\perp}}}K\ensuremath⊥​/K\ensuremath∥{\mathit{K}}_{\mathrm{\ensuremath{\parallel}}}K\ensuremath∥​, calculated for strong anisotropies from ``2D'' (3D) scaling for scales r0{\mathit{r}}_{0}r0​ (>r0{\mathit{r}}_{0}r0​), matches existing Monte Carlo data well. Critical regions \ensuremath{\Vert}\ensuremathεc{\mathrm{\ensuremath{\varepsilon}}}_{\mathit{c}}\ensuremathεc​\ensuremath{\Vert}\ensuremath{\le}\ensuremathγ0\ensuremath−1/\ensuremathν{\ensuremath{\gamma}}_{0}^{\mathrm{\ensuremath{-}}1/\ensuremath{\nu}}\ensuremathγ0\ensuremath−1/\ensuremathν​ are estimated. Contact is made with ideas of an intrinsic critical current arising from linear \ensuremath{\sim}{\ensuremath{\gamma}}_{0}^{\mathrm{\ensuremath{-}}1}$${\mathit{R}}_{\mathrm{\ensuremath{\parallel}}} vortex segment effective potentials.
The temperature scale of relevance for the vortex excitations in the two-dimensional XY model is identified. The size dependence of the helicity modulus at constant vortex temperature is then obtained by Monte Carlo simulations. The Monte Carlo data are then compared with the results for the length-dependent dielectric function from the Kosterlitz renormalization-group equations. The agreement is excellent for temperatures up to and slightly above Tc{\mathit{T}}_{\mathit{c}}Tc​. The critical temperature is determined to be Tc{\mathit{T}}_{\mathit{c}}Tc​=0.892 13(10). The temperature dependence of the characteristic length right above Tc{\mathit{T}}_{\mathit{c}}Tc​ is also obtained and found to be in agreement with recent results from high-temperature series expansions.
We present detailed Monte Carlo results for the susceptibility \ensuremath{\chi}, correlation length \ensuremath{\xi}, and specific heat Cv{\mathit{C}}_{\mathit{v}}Cv​ for the XY model. The simulations are done on 642{64}^{2}642, 1282{128}^{2}1282, 2562{256}^{2}2562, and 5122{512}^{2}5122 lattices over the temperature range 0.98\ensuremath{\le}T\ensuremath{\le}1.43 corresponding to 270. Fits to \ensuremath{\chi} and \ensuremath{\xi} data favor a Kosterlitz-Thouless (KT) singularity over a second-order transition; however, unconstrained four-parameter KT fits do not confirm the predicted values \ensuremath{\nu}=0.5 and \ensuremath{\eta}=0.25. Our best estimate, Tc{\mathit{T}}_{\mathit{c}}Tc​=0.894(5), is obtained using KT fits with \ensuremath{\nu} fixed at 0.5. The exponent \ensuremath{\eta} is calculated as a function of temperature in the spin-wave phase using Monte Carlo renormalization-group and finite-size-scaling methods. Both methods give consistent results and we find \ensuremath{\eta}\ensuremath{\approxeq}0.235 at T=0.894. We also present results for the behavior of vortex density across the transition and exhibit how the dilute-gas approximation breaks down.
In this master thesis the criticality of the 3D XY model with cubic anisotropy has been examined by the means of Monte Carlo simulation. The 3D XY model is a special case of the O(n)O(n)O(n) spin-vector model which describe numerous phenomena's. Here, nnn is the number of vector dimensions and for the XY model n=2n=2n=2. The O(n)O(n)O(n)-model is invariant for global rotations. This means that it is arbitrary which direction in the symmetry breaking occurs. In real system one should always expect perturbations from a perfect system, here it is done by adding a cubic anisotropy (order p=4p=4p=4) then the symmetry of the Hamiltonian is lowered and it will be advantageous for the spins to be ordered in directions 90o90^o90o apart. This symmetry is called C4C_{4}C4​, which is a part of the cyclic group, and it restricts the symmetry breaking to one of these four directions. \n\nThe criticality of this model is when the spins suddenly align as the temperature is lowered to a critical temperature. Sufficiently close to the critical temperature quantities like the magnetization, which is the sum of spins, will show a power law growth characterised by critical exponents specified by the class of universality of the transition. The critical exponents can be obtained by finite-size scaling. By applying renormalization theory on the critical region one can show that some of parameters in the Hamiltonian are: relevant and effect the scaling, irrelevant and does not effect the scaling and dangerously irrelevant meaning that they only affect some of the critical exponents. Previous Monte Carlo and renormalization calculus show that anisotropy which is relevant for p<4p<4p<4 becomes dangerously irrelevant for pgeq4p\\geq4pgeq4. The question of interest is what kind of effect the dangerously irrelevant parameter in the Hamiltonian has on the scaling. This is interesting because the effect of dangerously irrelevant parameters is hardly studied and anisotropy arise in many crystal systems because of the lattice geometry. \n\nThe results are obtained by Monte Carlo simulations and the critical slowing down close to the phase transition is solved by the method of parallel tempering. This method was chosen because it is simple to implement in a parallel code.\n\nThe Binder cumulant for different quantities and the peak of the magnetic susceptibility was used support a critical temperature Tc=2.2019textnormalkT_c = 2.2019 \\ \\textnormal kTc​=2.2019textnormalk for no field applied and no anisotropy. The known critical exponents for the no field and no anisotropy was used. The scaling made a good collapse, meaning that the right TcT_cTc​ and the correct critical exponents where used. The scaling did not take into account a magnetic field and therefore a good collapse of curves where not obtained in this case. When the cubic anisotropy was added the critical exponent for the magnetic susceptibility and the cubic order parameter changed, while the other exponents where unchanged. This supports the view that cubic anisotropy is dangerously irrelevant for the 3D XY model. The scaling of the magnetic susceptibility showed a possible new kind of scaling, a scaling which is different on the two sides of the transition. Although a more accurate determination of the critical exponent in the ordered phase require more extensive simulations. It appears that this exponent is larger than the exponent in the disordered phase, where one recovered the isotropic exponent. The difference in the critical exponents has been predicted to increase with ppp and an investigation of the hexagonal anisotropy p=6p=6p=6 has been started.
Monte Carlo simulations have been used to study a version of the XY model without vortex loops for J/T=0 on simple cubic and face-centered-cubic lattices. Finite-size scaling of data from L\ifmmode\times\else\texttimes\fi{}L\ifmmode\times\else\texttimes\fi{}L lattices with L up to 32 is used to obtain values for the order parameter \ensuremath{\Vert}M\ensuremath{\Vert} and a 1/L finite-size correction on each lattice. The order-parameter susceptibility is found to be diverging as L0.8{\mathit{L}}^{0.8}L0.8. We then study the crossover from three-dimensional to one-dimensional behavior on L\ifmmode\times\else\texttimes\fi{}W\ifmmode\times\else\texttimes\fi{}W simple cubic lattices, with W=4 and 8, and L up to 1024.
Motivated by the interplay between 2D and 3D scaling signatures observed in unconventional layered superconductors, we present a systematic Monte Carlo study of the three-dimensional classical XY model with anisotropic in-plane J∥J_\parallelJ∥​ and inter-plane J⊥J_\perpJ⊥​ couplings. Our study includes very small values of the system anisotropy Δ=J⊥/J∥Δ=J_\perp /J_\parallelΔ=J⊥​/J∥​ not studied before, and focuses on characterizing the crossover from quasi-2D topological scaling to genuine 3D critical behavior. The numerical results for the critical temperature unambiguously reveal a logarithmic scaling with ΔΔΔ, directly related to the topological scaling in the 2D limit. Despite the 3D nature of the layered XY criticality, topological scaling signatures survive up to system sizes comparable to the crossover length ℓJ\ell_JℓJ​, which diverges at small ΔΔΔ with a scaling behavior reminiscent of the Berezinskii-Kosterlitz-Thouless (BKT) transition. This shows that genuine 3D symmetry-breaking behavior emerges only at exceedingly large system sizes when the anisotropy is very strong. Our results indicate that new experimental evidence is required to clarify the extent to which the critical signatures observed in layered strongly correlated materials are shaped by their pronounced anisotropy.
We have studied the critical properties of three-dimensional U(1)-symmetric lattice gauge theories. The models apply to various physical systems such as insulating phases of strongly correlated electron systems as well as superconducting and superfluid states of liquid metallic hydrogen under extreme pressures. This thesis contains an introductory part and a collection of research papers of which seven are published works and one is submitted for publication.\n\nPaper I: Critical properties of the 2+1-dimensional compact abelian Higgs model with gauge charge q=2 are studied. We introduce a novel method of computing the third moment M3 of the action which allows us to extract correlation length and specific heat critical exponents ν and α without invoking hyperscaling. Finite-size scaling analysis of M3 yields the ratio (1+α)/ν and 1/ν separately. We find that α and ν vary along the critical line of the theory, which however exhibits a remarkable resilience of Z2 criticality. We conclude that the model is a fixed-line theory, which we propose to characterize the zero temperature quantum phase transition from a Mott-Hubbard insulator to a charge fractionalized insulator in two spatial dimensions.\n\nPaper II: Large scale Monte Carlo simulations are employed to study phase transitions in the three-dimensional compact abelian Higgs model in adjoint representations of the matter field, labeled by an integer q, for q=2,3,4,5. We also study various limiting cases of the model, such as the Zq lattice gauge theory, dual to the 3DZq spin model, and the 3D xy spin model which is dual to the Zq lattice gauge theory in the limit q → ∞. In addition, for benchmark purposes, we study the 2D square lattice 8-vertex model, which is exactly solvable and features non-universal critical exponents. The critical exponents α and ν are calculated from finite size scaling of the third moment of the action, and the method is tested thoroughly on models with known values for these exponents. We have found that for q=3, the three-dimensional compact abelian Higgs model exhibits a second order phase transition line which joins a first order phase transition line at a tricritical point. The results for q=2 in Paper I are reported with a higher lever of detail.\n\nPaper III: This paper is based on a talk by F. S. Nogueira in the Aachen HEP 2003 conference where a review of the results for the compact abelian Higgs model from Paper I and Paper II was presented, as well as the results for the q=1 case studied by F. S. Nogueira, H. Kleinert and A. Sudbø.\n\nPaper IV: We study the effects of a Chern-Simons (CS) term in the phase structure of two different abelian gauge theories in three dimensions. By duality transformations we show how the compact U(1) gauge theory with a CS term for certain values of the CS coupling can be written as a gas of vortex loops interacting through steric repulsion. This theory is known to exhibit a phase transition governed by proliferation of vortex loops. We also employ Monte Carlo simulations to study the non-compact U(1) abelian Higgs model with a CS term. Finite size scaling of the third moment of the action yields critical exponents α and ν that vary continuously with the strength of the CS term, and a comparison with available analytical results is made.\n\nPaper V: The critical properties of N-component Ginzburg-Landau theory are studied in d=2+1 dimensions. The model is dualized to a theory of N vortex fields interacting through a Coulomb and a screened potential. The model with N=2 shows two anomalies in the specific heat. From Monte Carlo simulations we calculate the critical exponents α and ν and the mass of the gauge field. We conclude that one anomaly corresponds to an inverted 3D xy fixed point, while the other corresponds to a 3D xy fixed point. There are N fixed points, namely one corresponding to an inverted 3D xy fixed point, and N-1corresponding to neutral 3D xy fixed points. Applications are briefly discussed.\n\nPaper VI: The phase diagram and critical properties of the N-component London superconductor are studied both analytically and through large-scale Monte-Carlo simulations in d=2+1 dimensions. The model with different bare phase stiffnesses for each flavor is a model of superconductivity which should arise out of metallic phases of light atoms under extreme pressure. A projected mixture of electronic and protonic condensates in liquid metallic hydrogen under extreme pressure is the simplest example, corresponding to N=2 with individually conserved matter fields. We compute critical exponents α and ν for N=2 and N=3. The results from Paper V are presented at a higher level of detail. For the arbitrary N case, there are N fixed points,namely one charged inverted 3D xy fixed point, and N-1 neutral 3D xy fixed points. We explicitly identify one charged vortex mode and N-1 neutral vortex modes. The model for N=2 and equal bare phase stiffnesses corresponds to a field theoretical description of an easy-plane quantum antiferromagnet. In this case, the critical exponents are computed and found to be non 3D xy values. Furthermore, we study the model in an external magnetic field, and find a novel feature, namely N-1 superfluid phases arising out of N charged condensates. In particular, for N=2 we point out the possibility of two novel types of field-induced phase transitions in ordered quantum fluids: i) A phase transition from a superconductor to a superfluid or vice versa, driven by tuning an external magnetic field. This identifies the superconducting phase of liquid metallic hydrogen as a novel quantum fluid. ii) A phase transition corresponding to a quantum fluid analogue of sublattice melting, where a composite field-induced Abrikosov vortex lattice is decomposed and disorders the phases of the constituent condensate with lowest bare phase stiffness. Both transitions belong to the 3D xy universality class.\n\nPaper VII: We consider the vortex superconductor with two individually conserved condensates in a finite magnetic field. The ground state is a lattice of cocentered vortices in both order parameters. We find two novel phase transitions when temperature is increased at fixed magnetic field. i) A "vortex sublattice melting" transition where vortices in the field with lowest phase stiffness ("light vortices") loose cocentricity with the vortices with large phase stiffness ("heavy vortices"), entering a liquid state (the structure factor of the light vortex sublattice vanishes continuously.) This transition is in the 3D xy universality class. ii) A first order melting transition of the lattice of heavy vortices in a liquid of light vortices.\n\nPaper VIII: We report on large-scale Monte Carlo simulations of a novel type of a vortex matter phase transition which should take place in a three dimensional two-component superconductor. We identify the regime where first, at a certain temperature a field-induced lattice of co-centered vortices of both order parameters melts, causing the system to loose superconductivity. In this state the two-gap system retains a broken composite symmetry and we observe that at a higher temperature it undergoes an extra phase transition where the disordered composite one-flux-quantum vortex lines are "ionized" into a "plasma" of constituent fractional flux vortex lines in individual order parameters. This is the hallmark of the superconductor-to-superfluid-to-normal fluid phase transitions projected to occur in e.g. liquid metallic hydrogen.