 the helicity modulus $\Gamma$ (related to the superfluid density) has an increasingly abrupt change with $1/J_\tau$ for increasing $L$, indicative of a jump as at a BKT transition. If the transition out of the CSF is indeed a BKT transition, $\Gamma$ can be fit to the finite size scaling form $\Gamma (L) \! =\!  A \left( 1\! + \! \frac{1}{2 (\log L + C)} \right)$ (with fit parameters $A$,$C$) right at the transition point, with $A$ taking on the universal value of $2/\pi$, while $C$ is a non-universal constant [@weber_minnhagen; @olsson.prl1995]. Fitting $\Gamma(L)$ to this form, we find that the error to this fit shows a sharp minimum [@weber_minnhagen; @olsson.prl1995] at a certain $1/J_\tau$ (Fig.\[Fig:classical\_trans\] inset), with $A\! \approx\! 2/\pi$ at this dip. This not only allows us to precisely locate the transition out of the CSF state, but also [*confirms*]{} its BKT nature.

![(Color online) (A) Helicity modulus $\Gamma$ versus $1/J_\tau$ for different system sizes for $J_\perp\!=\!1$. (A-Inset) RMS error of fit to the BKT finite size scaling form of $\Gamma$ shows a deep minimum [@supp] at the transition, at $1/J_\tau \!=\! 0.887(1)$, and yields a jump $\Delta \Gamma\! \approx\! 0.637$, close to the BKT value $2/\pi$. (B) Binder cumulants for the staggered current versus $1/J_\tau$ (for different $L$ for $J_\perp\!=\!1$) intersecting at a continuous transition at $1/J_\tau\! =\! 0.981(4)$. (B-inset) Critical susceptibility versus $L$ gives the ratio of critical exponents $\gamma/\nu \! \approx\! 1.72$, very close to 2D Ising value $\gamma/\nu\!=\!7/4$. Error bars are smaller than the symbol sizes. []{data-label="Fig:classical_trans"}](xy.eps){width="2.8in"
the helicity modulus $ \Gamma$ (related to the superfluid density) has an increasingly abrupt change with $ 1 / J_\tau$ for increase $ L$, indicative of a leap as at a BKT transition. If the transition out of the CSF is indeed a BKT transition, $ \Gamma$ can be fit to the finite size scaling shape $ \Gamma (L) \! = \!   A \left (1\! + \! \frac{1}{2 (\log L + C) } \right)$ (with fit parameters $ A$, $ C$) right at the conversion compass point, with $ A$ taking on the cosmopolitan value of $ 2/\pi$, while $ C$ is a non - cosmopolitan constant [ @weber_minnhagen; @olsson.prl1995 ]. Fitting $ \Gamma(L)$ to this shape, we find that the error to this fit indicate a sharp minimum [ @weber_minnhagen; @olsson.prl1995 ] at a certain $ 1 / J_\tau$ (Fig.\[Fig: classical\_trans\ ] inset), with $ A\! \approx\! 2/\pi$ at this dip. This not only allow us to precisely locate the transition out of the CSF state, but besides [ * confirms * ] { } its BKT nature. 

! [ (Color online) (A) Helicity modulus $ \Gamma$ versus $ 1 / J_\tau$ for different arrangement sizes for $ J_\perp\!=\!1$. (A - Inset) RMS mistake of fit to the BKT finite size scaling form of $ \Gamma$ shows a deep minimum [ @supp ] at the transition, at $ 1 / J_\tau \!=\! 0.887(1)$, and yields a leap $ \Delta \Gamma\! \approx\! 0.637 $, close to the BKT value $ 2/\pi$. (B) Binder cumulants for the staggered current versus $ 1 / J_\tau$ (for different $ L$ for $ J_\perp\!=\!1 $) intersecting at a continuous transition at $ 1 / J_\tau\! = \! 0.981(4)$. (B - inset) Critical susceptibility versus $ L$ gives the ratio of critical exponents $ \gamma/\nu \! \approx\! 1.72 $, very close to 2D Ising value $ \gamma/\nu\!=\!7/4$. Error bars are smaller than the symbol sizes. [ ] { data - label="Fig: classical_trans"}](xy.eps){width="2.8 in "
 thf helicity modulus $\Gamma$ (related to the superfnuid dsnsity) hxs an increasingly abrupt chenge with $1/J_\tau$ for increasing $L$, indicatine of a jymp es at a BKT transition. In the bransntmon out of the GSF is indead a BKT transhtkou, $\Gamma$ can be fit to the finite sizq scalimg form $\Gamma (L) \! =\!  A \jeft( 1\! + \! \frac{1}{2 (\log L + C)} \right)$ (with fit pzrametegs $A$,$C$) right at thr transition point, with $A$ hakijg on the universap value of $2/\kj$, wrule $C$ is a nun-universal constant [@wgber_minnhagen; @olsson.prl1995]. Fitting $\Eamma(K)$ to this doem, ag find that vhe ergor to this fln shows a sharl minimum [@webev_minniageb; @olsson.prl1995] at a certein $1/J_\tau$ (Fig.\[Fig:classycal\_trans\] iuset), with $A\! \approx\! 2/\pi$ ar thiv di[. Thkw nut knky allowd ua to preciaely locate the transition out os the CSF stats, but wlfo [*confirms*]{} its BKT nature.

![(Color online) (D) Hslicity modulus $\Gamma$ vwrsus $1/J_\tau$ for differgnt system sizes for $J_\perp\!=\!1$. (A-Inset) RMS error of fit to the BNT fiiige wide szqllng form of $\Gamma$ shows a deep minimum [@supp] ae tne transition, at $1/J_\tau \!=\! 0.887(1)$, and ylekqs a jump $\Delja \Gammc\! \alprox\! 0.637$, close to thf BKT vwlue $2/\pi$. (B) Bindtr cukulants for the staggered cyrrent versuf $1/J_\tau$ (for differenc $L$ for $J_\per'\!=\!1$) intetsectimg at a continuous tranritikn at $1/J_\tau\! =\! 0.981(4)$. (B-inset) Dfitical suscepticilpty eersus $L$ gives the ratio os criticao ex'onents $\eamms/\nu \! \a[prox\! 1.72$, verj close to 2D Ising value $\gammc/\nu\!=\!7/4$. Esror bars wre smaller than the symbol sizxx. []{data-label="Fog:wlavsical_trcns"}](xy.eis){width="2.8in"
the helicity modulus $\Gamma$ (related to the has increasingly abrupt with $1/J_\tau$ for jump at a BKT If the transition of the CSF is indeed a transition, $\Gamma$ can be fit to the finite size scaling form $\Gamma (L) =\! A \left( 1\! + \! \frac{1}{2 (\log L + C)} \right)$ (with parameters right the point, with $A$ taking on the universal value of $2/\pi$, while $C$ is a non-universal constant @olsson.prl1995]. Fitting $\Gamma(L)$ to this form, we find the error to this shows a sharp minimum [@weber_minnhagen; at certain $1/J_\tau$ inset), $A\! 2/\pi$ at this This not only allows us to precisely locate the transition out of the CSF state, but also its BKT online) (A) modulus versus for different system $J_\perp\!=\!1$. (A-Inset) RMS error of fit finite size scaling form of $\Gamma$ shows a minimum [@supp] the transition, at $1/J_\tau \!=\! 0.887(1)$, yields a jump $\Delta \Gamma\! \approx\! 0.637$, close the BKT value $2/\pi$. (B) Binder cumulants for the staggered current versus $1/J_\tau$ (for different $J_\perp\!=\!1$) intersecting at a transition at $1/J_\tau\! 0.981(4)$. Critical versus gives the of critical exponents $\gamma/\nu \! \approx\! 1.72$, very close to 2D value $\gamma/\nu\!=\!7/4$. Error bars are smaller than the symbol sizes.
 the helicity modulus $\Gamma$ (reLated to the SuperFluId dEnSity) Has aN increasingly aBRupt Change with $1/J_\tau$ for increAsing $l$, iNDicaTIvE of a jUmp as at A bKt TRanSiTiOn. IF tHE tRansiTioN out of tHe CSF is indEed A BkT transition, $\gAmMa$ can be fit To tHe finite size ScaLing foRm $\gamMA (L) \! =\!  A \leFt( 1\! + \! \fRac{1}{2 (\loG L + C)} \rigHT)$ (with fIt parametErS $a$,$C$) righT At the trANSiTion Point, with $A$ taking oN ThE Universal value Of $2/\pi$, whIlE $c$ iS A Non-UniVersal consTaNt [@webER_minnhaGEn; @OLSSon.PRl1995]. Fitting $\GammA(L)$ to this forM, We fInd thaT tHe eRRor to tHis fiT sHOws A sharp minimUm [@weBer_minnhaGen; @olsSOn.prl1995] at A Certain $1/j_\tau$ (FiG.\[FiG:clAssiCAl\_TrAns\] InSEt), wITh $a\! \apPRox\! 2/\Pi$ at this DiP. THis noT onlY ALLOws uS to PrecIsely Locate the tranSitIon oUT of The CSf statE, but AlSo [*conFirms*]{} iTs BKT NaTure.

![(Color online) (a) HelIcity moduLus $\gaMma$ VeRsus $1/J_\TAu$ for dIffEreNt systeM sizes fOR $J_\pErP\!=\!1$. (a-iNsEt) RMS error of fit to tHe bkt fInite sizE scaliNG fOrM Of $\Gamma$ sHoWs a Deep MINimum [@Supp] AT tHe transiTion, at $1/j_\TaU \!=\! 0.887(1)$, aNd yieldS a Jump $\DeLtA \GaMma\! \ApproX\! 0.637$, ClosE to the bKT value $2/\Pi$. (B) BiNDer cumulants foR The staggered cURrENT vERsus $1/j_\taU$ (for differeNt $L$ fOR $J_\peRp\!=\!1$) inTErSecTIng at A contInUOuS Transition at $1/J_\tau\! =\! 0.981(4)$. (B-inSeT) CritiCal suSceptibility vErsus $L$ giveS THE ratio of CritICaL Exponents $\gamma/\Nu \! \appRox\! 1.72$, very cloSE to 2D IsinG valuE $\gamma/\nu\!=\!7/4$. error bars ARE smaller ThaN thE syMboL SIzEs. []{data-label="FiG:CLassIcAl_trans"}](Xy.ePs){width="2.8In"
 the helicity modulus $\Ga mma$ (rela ted t o t hesu perf luid  density) hasa n in creasingly abrupt chan ge wi th $1/J _ \t au$ f or incr e as i n g $ L$ ,ind ic a ti ve of  ajump as  at a BKTtra ns ition. If th e  t ransitionout  of the CSFisindeed  a  BK T  tran sit ion,$\Gamm a $ canbe fit to  t h e fini t e sizes c al ingform $\Gamma (L)\ != \!  A \left( 1 \! + \ !\ fr a c {1} {2(\log L +C) } \ri g ht)$ (w i th f i t p a rameters $A$, $C$) righta t t he tra ns iti o n poin t, wi th $A$  taking ontheuniversal  value of $2/\ p i$, whi le $C$  is  anon- u ni ve rsa lc ons t an t [ @ web er_minnh ag en ; @ol sson . p r l 1995 ].Fitt ing $ \Gamma(L)$ to  th is f o rm,  we f ind t hatth e err or tothisfi t shows a sharp  min imum [@we ber _m inn ha gen;@ olsson .pr l19 95] ata certa i n $ 1/ J _ \ ta u$ (Fig.\[Fig:clas si c a l\ _trans\]  inset ) ,wi t h $A\! \ ap pro x\!2 / \pi$at t h is  dip. Th is not on ly  allows  u s to p re cis ely  loca t e th e tran sition o ut of the CSF state, but also [*co n fi r m s* ] {} i tsBKT nature. 

![ ( Colo r on l in e)( A) He licit ym od u lus $\Gamma$ versus  $ 1/J_\t au$ f or differentsystem siz e s for $J_\ perp \ != \ !1$. (A-Inset)  RMSerror of f i t to the  BKTfinite s ize scali n g  form of  $\ Gam ma$  sh o w sa deep minimu m [@su pp ] at th e t ransiti on,  at  $1 /J_ \t au \!=\!0.887(1) $,  a nd  y iel ds aj ump $\De lt a \ Ga mma \! \a p prox\!  0.63 7$,cl os e  to  the BK T  v a l ue $ 2/ \p i$.(B)  B inder  cum u lan ts forthe stagg ere d  cur re nt  versus  $1/J_\tau$ ( fo r differen t$L$  for $ J _ \perp\!= \!1$) intersecting at a continu ous  tran siti on at $1/ J_\ tau\!=\! 0.981( 4)$. ( B-ins et ) C r i tical s us cep ti bility ver s u s $ L$ gi ve s th e ratio  of critical expon e nts  $\gamma/\nu\!\app r o x\ ! 1 . 72 $ , v er y  cl o s e to 2D Ising v alue $\gam ma / \n u\!=\!7/4$ .  Er ro r barsare sma llert han the  symbol s izes. []{ da ta-l a b el= "Fig:class ical_tra ns"}](xy. e ps){w i dt h="2. 8in "
 the_helicity modulus_$\Gamma$ (related to the_superfluid density)_has_an increasingly_abrupt_change with $1/J_\tau$_for increasing $L$,_indicative of a jump_as at a_BKT_transition. If the transition out of the CSF is indeed a BKT transition, $\Gamma$_can_be fit_to_the_finite size scaling form $\Gamma_(L) \! =\!  A_\left( 1\!_+ \! \frac{1}{2 (\log L + C)} \right)$_(with_fit parameters $A$,$C$)_right at the transition point, with $A$ taking on_the universal value of $2/\pi$, while_$C$ is a_non-universal_constant_[@weber_minnhagen; @olsson.prl1995]. Fitting $\Gamma(L)$_to this form, we find that_the error to this fit shows_a sharp minimum [@weber_minnhagen; @olsson.prl1995] at a_certain $1/J_\tau$ (Fig.\[Fig:classical\_trans\] inset), with $A\!_\approx\! 2/\pi$ at this dip._This not_only allows us to precisely_locate the transition_out of_the CSF state,_but also [*confirms*]{} its BKT nature.

![(Color_online) (A) Helicity_modulus $\Gamma$ versus $1/J_\tau$ for different_system_sizes for $J_\perp\!=\!1$._(A-Inset)_RMS_error of_fit to the_BKT_finite size_scaling_form of $\Gamma$ shows a deep_minimum_[@supp] at the transition, at $1/J_\tau \!=\!_0.887(1)$, and yields a_jump_$\Delta \Gamma\! \approx\! 0.637$,_close to the BKT value_$2/\pi$. (B) Binder cumulants for the_staggered current_versus $1/J_\tau$_(for different $L$ for $J_\perp\!=\!1$) intersecting at a continuous transition at_$1/J_\tau\! =\! 0.981(4)$. (B-inset) Critical susceptibility_versus $L$ gives the_ratio of_critical_exponents $\gamma/\nu \!_\approx\!_1.72$, very_close to 2D Ising value $\gamma/\nu\!=\!7/4$. Error_bars are_smaller than the symbol sizes. []{data-label="Fig:classical_trans"}](xy.eps){width="2.8in"