the existence and nature of melting transitions in two - dimensional ( 2d ) systems remains a topic of interest and uncertainty in both theory and experiments .
in particular , the melting of the vortex lattice in a thin film superconductor has received much attention . a continuous melting transition from an abrikosov lattice to a vortex liquid via a phase of hexatic order , as in the kosterlitz - thouless - halperin - nelson - young ( kthny ) theory , was first proposed over fifteen years ago.@xcite a `` transition '' in the current - voltage response of thin film superconductors has been seen in several experiments.@xcite the vortex system becomes depinned above some transition temperature ( the depinned vortices will dissipate energy as they move , leading to an ohmic response ) and the rise of the resistivity of the film is presumed to be caused by the melting of the vortex lattice
. however , it seems that this transition temperature varies with the strength of the pinning centers in the thin films . for films of amorphous nbo
the pinning is weak and no features can be detected in the current - voltage characteristics which may be attributed to a `` transition '' .
this has been interpreted as being due to the absence of the abrikosov vortex lattice.@xcite    the superconducting state arises from macroscopic correlations in the phase of the cooper pair wave function , or equivalently the ginzburg - landau order parameter .
this phase coherence is usually a consequence of off - diagonal long range order ( odlro ) , i.e.  non - vanishing @xmath0 , where @xmath1 is the superconducting order parameter .
it has been predicted that thermal excitations of the shear modes of the vortex lattice should destroy odlro in the crystalline state and the length scale over which phase coherence persists has been estimated.@xcite    oneill and moore have performed monte carlo ( mc ) simulations of thin film superconductors within the lowest landau level ( lll ) approximation using a spherical geometry to minimize finite size boundary effects.@xcite no freezing transition out of the vortex liquid phase was found in the temperature region investigated . the correlation function associated with the positions of the vortices
was calculated and the associated length scale of these lattice - like correlations was found to diverge only in the zero temperature limit . it was noted that this happened on a similar length scale over which there was expected to be phase coherence from the arguments of refs .  and  .
in the work of oneill and moore , the superconducting film existed on the surface of a sphere , and this geometry is the one employed here .
the magnetic field perpendicular to the surface is generated by a magnetic monopole placed at the center of the sphere.@xcite the strength of the monopole must be quantized to have a well defined superconducting order parameter , but the strength of the field at the surface may be continuously varied by changing the radius of the sphere . the lll approximation , which assumes an infinite magnetic screening length , works well in the 2d case where the effective screening length is @xmath2 because of the magnetic field outside the superconductor.@xcite @xmath3 becomes large as the film thickness @xmath4 is reduced ; @xmath5 is the bulk magnetic penetration length .
it is often the case that @xmath3 is greater than the film dimensions for small @xmath4 .
more recently , strong evidence for a first order transition at a finite temperature in an alternative finite size approximation with quasi - periodic boundary conditions ( qp ) on an infinite plane has come from mc simulations.@xcite this model has a periodic amplitude of the order parameter and a fixed change in the phase as a unit cell containing @xmath6 vortices is crossed ( the change in phase depends on the choice of gauge ) .
the quasi - periodicity in the phase places some constraints on the motion of the vortices , but it is not obvious that this will change the thermodynamic properties of the system .
the observed transition temperature is within the range investigated on the sphere .
this raises the question of whether or not there really is no transition on the sphere .
if so , then what are the reasons that different boundary conditions give different results ?
if there is a freezing transition , at least on the plane , does this contradict the theoretical prediction of no odlro at finite temperatures ?
apparently not , as the work of sasik , stroud and tesanovic@xcite ( sst ) has concluded that in the thermodynamic limit even the low temperature state below the first order transition has no odlro , despite the presence of `` bragg peaks '' in the correlations of the order parameter density ( the peaks will not be delta - functions but the algebraic singularities expected for a 2d crystalline system@xcite ) .
these conclusions of sst are taken as support for the `` charge - density wave '' phase ( with a modulation in the cooper pair density , but no odlro ) proposed by tesanovic.@xcite    in this work we again fail to find a phase transition in the spherical geometry .
we have investigated the correlations associated with both the _ phase _ and the _ amplitude _ of the superconducting order parameter .
one advantage of the spherical geometry is that the liquid state may be investigated over a larger range of temperatures .
this has allowed us to compare the growing length scales of the phase and amplitude correlations over a large variation in temperature .
section  [ sec : method ] contains the formulation of our model . in section
[ sec : pert ] it is shown how the high temperature perturbation theory on the spherical geometry gives the same free energy per vortex as the conventional theory in the thermodynamic limit . in section
[ sec : ground ] we describe how the ground state of the vortex system on the sphere is affected by the necessity of containing twelve disclination defects within the triangular lattice .
we have previously concluded that these defects with long range strains will give a finite energy cost per vortex relative to the flat plane ground state even in the thermodynamic limit.@xcite however , the presence of dislocation defects were not considered , and these can screen the strains of the disclinations to recover the thermodynamic limit of the triangular lattice on a flat plane .
this screening is described in the framework of elasticity theory , and the numerical evidence for dislocations in the ground state is reexamined . in section  [ sec :
density ] our calculations on the four - point density correlation function of the vortex liquid are shown .
this is followed by our investigations on the phase coherence in section  [ sec : phase ] .
the spherical geometry is in principle realizable using the end of a long thin solenoid to approximate the field from a monopole and placing it inside a spherical superconducting film .
however there would always be a problem in the placement of the solenoid which would lead to a large scale inhomogeneity in the field strength at the spherical surface .
additional simulations are described in section  [ sec : inhom ] , where we add just such an inhomogeneous magnetic field .
finally , in the light of different conclusions from simulations of the lll system in different geometries , we discuss some of the problems in these geometries in section  [ sec : discussion ] and comment on the experimental situation .
some conclusions are drawn in section  [ sec : conclusions ] .
we model a thin film superconductor by a spherical shell of radius @xmath7 and thickness @xmath8 .
an external magnetic field @xmath9 , always directed perpendicular to the shell s surface is created by a magnetic monopole at the center of the sphere .
we neglect all fluctuations in the magnetic induction so that it always has the mean value @xmath10 on the superconducting surface .
dirac s quantization condition@xcite ( assuming that the basic unit of charge is that of a cooper pair ) requires that the total flux through the surface is an integer multiple of flux quanta @xmath11 .
our choice of gauge that satisfies @xmath12 is @xmath13 .
we also assume that the superconductor is described by a ginzburg - landau ( gl ) complex order parameter @xmath1 with a phenomenological free energy hamiltonian given by @xmath14= \int d^3r \left [ \alpha(t ) { |\psi|}^2                   + \frac{\beta}{2}{|\psi|}^4   + \frac{1}{2 m } \psi^ * { d}^2 \psi \right],\label{eq : glenergy}\ ] ] where @xmath15 and @xmath16 is the gauge invariant derivative operator .
the lll approximation consists in expanding this order parameter in terms of the eigenstates of @xmath17 , and then only keeping the degenerate level of states with the lowest eigenvalue .
if we only consider the two dimensions on the surface of the sphere then the subspace of this lll is spanned by @xmath18 orthonormal functions:@xcite @xmath19 with @xmath20 and @xmath21}^{1/2}$ ] .
we write the order parameter as @xmath22 with @xmath23 , and measure lengths in units of the magnetic length @xmath24 which fixes @xmath25 .
[ note that the lattice spacing @xmath26 in a triangular vortex lattice will be given by @xmath27 . ] the lll hamiltonian can then be written as @xmath28 , \nonumber\end{aligned}\ ] ] with @xmath29 .
the quartic coupling term is given by@xcite @xmath30 where @xmath31 and @xmath32 $ ] is the binomial coefficient .
our effective reduced temperature variable is @xmath33 = \frac{dq^2}{k_bt}\alpha ' \left[-1+\frac{t}{t_{c0}}+\frac{h}{h_{c2}(0)}\right].\ ] ] where @xmath34 and @xmath35 is the mean field transition temperature .
@xmath36 is the straight line extrapolation to zero temperature of the upper critical field , within mean field theory . for better computational efficiency , eq .
( [ eq : oneill ] ) maybe rewritten as@xcite @xmath37   , \ ] ] where @xmath38 and @xmath39 is the beta function .
@xmath40 is the heaviside step function .
an alternative description of the order parameter on the sphere in the lll approximation is possible if we write ( following ref .  )
@xmath41 , where @xmath42 .
the @xmath6th order polynomial may be written as a product of @xmath18 terms to give : @xmath43 therefore the order parameter is determined by the positions of the @xmath6 zeros ( vortices ) @xmath44 along with an overall complex amplitude @xmath45 .
it is shown in appendix  [ ap : inv ] that the gl free energy hamiltonian depends only on the _ relative _ positions of these vortices plus the overall complex amplitude , and not on the choice of coordinate system @xmath46 that arose from our choice of gauge .
this translational invariance is assumed throughout this paper .
note that there are no constraints on the allowed positions @xmath44 which is not the case with qp boundary conditions .
the properties of these systems are determined by the partition function @xmath47}$ ] .
we investigate these properties , in particular correlations in the `` density '' @xmath48 and the phase @xmath49 , with the metropolis algorithm which uses an mc method to sample the phase space.@xcite we have made measurements for runs up to @xmath50 mc steps for system sizes as large as @xmath51 .
we have mainly concentrated on the temperature range between @xmath52 and @xmath53 .
we have also used standard minimization routines to investigate the ground states of this model .
a perturbation expansion has been developed about the gaussian limit of ( [ eq : glenergy ] ) at @xmath54 to approximate various properties of the lll system.@xcite in two dimensions the geometry used has always been a flat plane .
the series may be represented by feynman diagrams and many quantities in two and three dimensions have been evaluated to high order.@xcite the series can extend to negative @xmath55 as long as the propagator is renormalized at least to the level of the hartree - fock approximation ( hf ) .
we consider the perturbation series for the free energy using the spherical geometry that leads to the hamiltonian ( [ eq : ham ] ) .
the usual treatment follows for the feynman diagrams in a @xmath56 theory , with directed lines representing the gaussian propagator @xmath57 connecting vertices representing the interaction @xmath58 , and internal lines are summed from @xmath59 to @xmath6 .
all diagrams with loops are removed by renormalizing the propagator in a self - consistent fashion .
this is the hf approximation , and the renormalized propagator , @xmath60 , must satisfy @xmath61 where @xmath62 . on the sphere we find : @xmath63 so that in the large @xmath6 limit on the sphere , the hf propagator approaches that on the flat plane .
exactly the same diagrams appear in the free energy expansion for the spherical geometry as those for the plane .
the contribution of each diagram with @xmath64 vertices to the coefficient of the @xmath64th order term of the expansion is a product of a combinatorial factor and the integral / sum over the `` momenta '' along each line of the interactions at each vertex .
the only differences on the sphere are these sums .
we now look at the lowest order feynman graph without loops , shown in fig .
[ fig : feyn]a , which has two vertices connected by four lines .
the sum that appears in this diagram is @xmath65 . on the plane
this is equal to @xmath66 , and we find on the sphere : @xmath67,\ ] ] which also gives the result on the plane with a @xmath68 correction .
generalizing our method to consider higher order terms , we calculate a particularly simple type of diagram to arbitrary order : the ring diagrams shown in fig .  [ fig : feyn]b .
for such a ring diagram of order @xmath64 , the relevant sum is @xmath69 .
we find on the sphere : @xmath70,\end{aligned}\ ] ] where the result on a plane from ref .
is @xmath71 . so the contribution to the free energy to @xmath64th order approaches that on the plane as long as @xmath64 is much less than the system size @xmath6 .
although we have not verified this for more complex diagrams , we expect a similar result .
this is because the interactions at each vertex have only @xmath68 ( and higher power ) corrections on the sphere from the interactions on the plane , so with @xmath64 vertices in a diagram , there should be a total correction of order @xmath72 .
therefore as long as the perturbation theory is useful ( ie . on the high temperature side of any phase transition ) the thermodynamic limit on the sphere should be the same as that for the flat plane . intuitively , it is also reasonable that the free energy per vortex in the liquid phase should be the same on the flat plane and on the sphere . in the liquid phase ,
correlation lengths are finite ( see section  [ sec : sim ] ) , and so long as they remain small compared to the sphere radius , then the liquid on the sphere is `` unaware '' of its curvature .
= 8 cm     +
the ground state of the lll vortex system on an infinite plane was long ago discovered to be the triangular lattice.@xcite this is also the ground state on the finite systems used in the qp simulations when commensurate boundary conditions are chosen . however , on a sphere the closest configuration the vortices can make to an ideal triangular lattice must contain twelve `` disclinations '' , i.e.  twelve vortices that only have five nearest neighbors .
we have previously studied the ground states by numerical minimization of the hamiltonian in eq .
( [ eq : ham ] ) for different system sizes.@xcite the following conclusions were made : there exist some values of @xmath6 ( magic numbers ) for which the ground state has a lower energy than nearby values , and the energies as @xmath73 do not seem to converge to the infinite plane value .
extrapolation of our numerical results gives an extra energy on the sphere in this large @xmath6 limit of @xmath74% of the ground state energy per vortex on the infinite plane .
the vortex configurations of these ground states make up a triangular network , but with twelve five - coordinated centers .
the result that as @xmath73 the ground states do not look the same as for an infinite plane may be explained by the presence of the disclinations ( five - coordinated centers ) , which are defects with long range strains.@xcite therefore there will always be a finite fraction of the vortex system that is distorted from an ideal triangular lattice .
this is supported by calculations within elasticity theory ( which assumes only small distortions ) that give a similar finite energy cost of the distortions on a sphere in the large @xmath6 limit when an icosahedral arrangement of the disclinations is assumed.@xcite    this result of an extra energy per vortex in the thermodynamic limit would be equally valid at finite temperatures as long as there is a non - zero shear modulus ( ie . in a solid phase ) .
if it were true there could be a paradox : imagine a 2d system with a continuous phase transition from a liquid to a crystalline solid phase . above the transition
the free energies per vortex have been shown in section  [ sec : pert ] to be the same on the sphere and the plane and so this transition must remain at the same temperature in each case .
as at the transition the free energy is continuous , then the free energies in the solid state must also be equal .    in fact , one can do a better job on the ground states on the sphere by having a number of dislocations ( which are bound pairs of five - fold and seven - fold disclinations ) in addition to the twelve free disclinations .
this means we have to take away the restriction that every vortex is six - fold coordinated other than the twelve five - fold centers used in the elasticity calculation of ref .  .
there is no topological problem with this as one can pull an opposite pair of dislocations out of a perfect lattice .
that extra dislocation defects may lower the energy may be seen in two different ways .
first , the interaction energy between a dislocation and a disclination may be negative , depending on the direction of the dislocation and this interaction energy grows faster with system size than the self energy of the dislocation .
alternatively , consider a region on a sphere bounded by three great circles that intersect at neighboring disclinations . even in the large @xmath6 limit this region will be distorted from a triangle ( all three of its angles will be @xmath75 ) and so a perfect triangular lattice in this region will be incommensurate with the boundaries , where a series of dislocations will occur .
so we have a picture of a ground state in the large @xmath6 limit made up of twenty regions of perfect lattice , but with sides containing dislocations , and corners containing disclinations .    there are still competing effects which discourage dislocations , and which may explain why they were not seen in ref .  .
first there is the curvature of the sphere which favors the distortions of a disclination , but which decreases as @xmath76 .
second , there is the positive self energy of a dislocation , which is generally expected to increase as the logarithm of the system size.@xcite third there is the interaction between two dislocations which is directional , and approximately proportional to the logarithm of the distance between them . to estimate the density of dislocations needed to reduce the strain energy of the disclinations we use the framework of elasticity theory .
the elasticity theory of twelve disclinations on a sphere is described in ref .  .
the method is briefly explained here and then generalized to allow the inclusion of dislocations .
we consider small strains from an ideal lattice , @xmath77 .
these strains are related to a stress field by hooke s law , @xmath78 . for the lll system
the elastic constant @xmath5 diverges ( the system is incompressible ) and the shear modulus @xmath79 is given by @xmath80 .
the stress tensor has zero divergence which allows the problem to be expressed in terms of the airy stress function , @xmath81 , defined by @xmath82 .
five - fold disclinations are defined by the change in bond angle by @xmath83 when a path encircles one such defect .
dislocations , on the other hand , are defined by their burgers vector , but they are also equivalent to a disclination dipole pair .
the stress field corresponding to a particular configuration of these topological defects is given by the solutions to the biharmonic equation on a sphere : @xmath84 where @xmath85 is the 2d young s modulus which in the lll system is given by @xmath86 .
the free disclinations are labeled by @xmath87 , and @xmath4 labels the dislocations with burgers vectors * b*@xmath88 . in the large @xmath7 limit where we neglect the energy of bending the lattice over the spherical curvature ,
the energy of a given configuration is @xmath89 where the integrals are over the surface of the sphere .
we will need to use coordinate frames rotated so that the reference axis passes through the defect .
for instance , if a disclination is at the position @xmath90 in our original coordinate frame , a general point @xmath46 is described in the rotated frame by the coordinates @xmath91 ( the polar angle from the disclination ) and @xmath92 ( the azimuthal angle in the rotated frame relative to the @xmath93 direction ) , where @xmath94 for the reference direction of the azimuthal angle @xmath92 in the new frame also to be rotated , a further transformation @xmath95 is required ( as in the standard treatment of euler angles@xcite ) .
we use this extra rotation for the dislocation defects , but for simplicity we do not define the transformation here .    the solution of eq .
( [ eq : biharm ] ) in the absence of dislocations was found to be @xmath96 where the @xmath87th disclination is at @xmath90 and @xmath97 + 1 \}.\ ] ] the elastic energy cost of the twelve disclinations is @xmath98 ^ 2 $ ] , which is proportional to the area of the system , and gives a finite energy cost per vortex .
if we now consider additional dislocation defects , we may write @xmath99 where @xmath100 by differentiating the solution for a disclination , we find : @xmath101 where we take the size of the burgers vector to be the lattice spacing @xmath26 , and @xmath102 is the azimuthal angle about the axis through the dislocation measured from the direction of the burgers vector .    within elasticity theory
, we can consider separately the contributions of the self energy of each defect , and the pairwise interaction energy between the defects .
the self energy of a dislocation on a sphere is @xmath103\right\ } + e_{core},\ ] ] where @xmath104 is a cutoff of the order of the lattice spacing ( within @xmath104 of the dislocation , the linear elasticity theory must break down completely ) and @xmath105 is the energy of the distortion inside this cutoff . as expected , this becomes proportional to @xmath106 in the limit of large @xmath7 , where : @xmath107 + e_{core}.\ ] ]    in the presence of more than one dislocation the total elastic energy will be @xmath108 which is a sum of the self energies plus pairwise interactions between every defect ( the interactions between disclinations are included in @xmath109 ) .    if we consider a disclination at @xmath110 and a dislocation at @xmath111 with a burgers vector pointing perpendicular to the @xmath112 line ( such that the five - seven pair in the dislocation point towards the disclination ) , we find an interaction energy : @xmath113\right].\ ] ] allowing the dislocation to rotate by an angle @xmath114 will give the above energy multiplied by a factor @xmath115 .
= 8 cm     +    the interaction energy between two dislocations , one at @xmath116 , the other at @xmath111 , with opposite burgers vectors perpendicular to the line joining the dislocations , is found to be @xmath117 } { ( 1+\cos{\theta_d } ) } \right].\ ] ] again as expected , the total energy of a dislocation pair , @xmath118 , in the large @xmath7 limit becomes proportional to the logarithm of the distance between the two opposite dislocations .
notice that the interaction energy between a dislocation and disclination increases more rapidly with @xmath7 than the self energy and the dislocation - dislocation interaction energy .
as @xmath119 may be negative ( depending on the orientation of the burgers vector ) , there is an instability to the formation of dislocations , for large enough system size .
if we consider the twelve disclinations to be fixed at the corners of an icosahedron , we find the energy landscape of a single dislocation to have minima along the lines that join the disclinations . the energy @xmath120 ( which is a sum of the interactions between the dislocation and the twelve disclinations )
is shown in fig .
[ fig : disloc ] as a function of position on one of these lines , with a burgers vector that is perpendicular to the line .
the position is denoted by the polar angle @xmath121 from one disclination , with a neighboring disclination located at @xmath122 $ ] . by symmetry
, the energy is zero exactly on a disclination and at the midpoint between two disclinations , and there is one negative stationary point .    for an estimate of
when dislocations should first start to appear in the ground states as @xmath7 increases , we consider the energy to put two dislocations of opposite burgers vectors on each line joining disclinations , at the positions that minimize @xmath119 .
such a configuration has an elastic energy cost of @xmath123 \nonumber\\ & & + 60e_{core},\end{aligned}\ ] ] where we assume @xmath124 and have only taken into account the dislocation - dislocation interactions between pairs on the same line . if we neglect @xmath105 we find that @xmath125 when @xmath126 which corresponds to @xmath127
however , with an estimate of @xmath128 , this critical value increases to @xmath129 .
the critical value of @xmath6 is sensitive to the estimated values of the core energy , and the neglect of extra dislocation pair interactions .
also , we have neglected the local curvature , which for small systems will favor the presence of the disclination strains , and therefore disfavor dislocations that screen these strains .
our estimates are compared with numerical results later in this section .
having established that dislocations may lower the total elastic energy cost on the sphere for large enough systems , it remains to show that large numbers of dislocations may screen out the strain field of the disclinations so as to cancel the @xmath130 dependence of the elastic energy .
= 8 cm     +    we consider series of dislocations along the lines between the disclinations .
the sum of the strains from each defect may be approximated at large distances by an integral .
for such a series with equally spaced dislocations a distance @xmath131 apart , the resulting strain fields are those from a positive and a negative disclination at each end of the line , but with a charge of size @xmath132 ( where @xmath133 for a five - fold disclination ) .
this charge arises from the ratio of the distance between the five - seven pair in a dislocation to the distance between dislocations in the series .
this is seen by writing the sum of the dislocation strain fields as an integral .
a series of dislocations running from @xmath134 to some point @xmath135 will have a total strain field of @xmath136 because the strain field @xmath137 of a dislocation at @xmath138 with burgers vector perpendicular to the @xmath139 direction is just the differential of the strain field @xmath140 of a disclination : @xmath141 $ ] .
if we imagine such series with @xmath142 between each disclination and the midpoints to neighboring disclinations , as is shown near one disclination in fig .
[ fig : screen ] , then the strain field of each disclination will be exactly screened away .
however , there will be an effective disclination of charge @xmath143 at each midpoint [ the midpoint makes a polar angle @xmath144 with the disclination , so that @xmath145 . what is needed is a variable spacing that starts from @xmath146 at the disclination , but diverges at the midpoint .
we may require that the spacing has the scaling form @xmath147 with the condition that @xmath148 ( this form is chosen as any higher power of @xmath7 would lead to a spacing that increases faster than the length of the series ) .    in the large @xmath6 limit an exact cancellation of the @xmath130 terms in the energy
is obtained .
the strain field from one series starting at @xmath149 and lying along @xmath139 is given by @xmath150^{\theta_m}_0 - \frac{l_0}{r}\int_0^{\theta_m }    \sigma^s_{\theta ' } ( \theta,\phi ) \frac{dg^{-1}}{d\theta'}d\theta'\nonumber\\ & = & -\frac{1}{5}\sigma^s_0 + \frac{1}{[5+(r / l_0)g(\theta_m ) ] } \sigma^s_{\theta_m } -\frac{l_0}{r}i(\theta,\phi),\end{aligned}\ ] ] where @xmath151 , which remains of order one as the system size increases .
in the large @xmath7 limit we write the strain field of the @xmath152th series that is connected to the disclination @xmath153 as @xmath154 where @xmath152 runs from one to five . the elastic energy cost due to all the defects , @xmath155 is @xmath156 which shows that all of the terms of order @xmath130 cancel , leaving terms of order one which depend on the gradients of @xmath157 and the value of @xmath158 .
our approximation in simplifying the effective charge at the midpoint to be proportional to @xmath159 ignores a correction of order @xmath76 to this charge . as the energy of a disclination is proportional to the area of the system ,
the energy cost of this ignored extra charge is also of order one .
therefore , the production of these series of dislocations of variable spacing , along the lines between disclinations , is an effective way of screening the diverging strains of the disclinations .
the ground state energy per vortex in the large @xmath6 limit is the same as that for an infinite flat plane apart from a correction proportional to @xmath160 due to the core energies of the dislocations .
we now consider the numerical evidence for the above conclusions .
it is sufficient to minimize the hamiltonian of eq .
( [ eq : ham ] ) by minimizing the abrikosov ratio @xmath161 the triangular lattice on the flat plane has the value@xcite @xmath162 . in ref .
, the ground states of system sizes up to @xmath163 were found numerically , and it was suggested that there was a finite energy cost per vortex in the large @xmath6 limit , consistent with the presence of twelve disclinations . the existence of magic numbers was observed and explained , although for large systems this only has a small effect on the total energy .
we have extended this investigation to larger systems ; this is motivated by the above elasticity results which suggest that above a critical size , dislocations will appear in the ground states and these will screen the strains of the disclinations to recover the large @xmath6 limit ground state energy of the flat plane . in fig .
[ fig : screen ] the sort of configuration is shown that we expect to find near each disclination .
this figure demonstrates just how large the spherical system must be to accomodate the dislocations and screen out the strains over a large region .    in fig .
[ fig : beta ] our numerical results are shown for ground states of systems larger than @xmath163 up to @xmath164 .
we also include results from ref .   for smaller system sizes .
it is important to note that , as the systems get large , there is an increasing number of metastable states that the minimization routines may fall into . in each case
we have taken the lowest energy after several hundred runs , but this is always just an upper bound on the true ground state energy .
= 8 cm     +    plotting the abrikosov ratio against @xmath160 the results for @xmath165 fall on a straight line that extrapolates to give a finite difference for the flat plane as @xmath73 .
however , this trend changes above @xmath166 , with energies falling below the straight line .
the results are explained by the reasoning that once systems are large enough for dislocations to overcome their self energy and mutual repulsion , they may screen the strains of the disclinations and reduce the total ground state energy . in the numerical ground states with @xmath167 we observe dislocations near the disclinations with the correct orientation of burgers vector (
the distribution of the dislocations may be irregular due to the difficulties in locating the true global ground state ) .
the elasticity calculation combined with this numerical work is strong evidence that the ground state energy per vortex in the thermodynamic limit is the same on the sphere as on the flat plane .
it should be mentioned however that for the system sizes we look at in mc simulations , the ground states have a higher energy per vortex by @xmath168-@xmath169% of the ground state energy on an infinite plane .
the jump in the energy at the transition observed in the qp systems is of a similar order of magnitude ( @xmath170 ) as the energy difference of the spherical ground states from the qp ground states .
it must be a possibility that the raising of the ground state energy inhibits a freezing transition to the low temperature state , which could explain the apparent absence of a phase transition in our simulations on the sphere .
we now go on to measure the correlations in the 2d vortex system using mc simulations . in the following , the system sizes we often use correspond to the magic numbers with particularly low energy ground states .
this is because a phase transition would be most likely where the ground state is least frustrated .
also , where we start the simulations in the ground state , we have most confidence that our numerical ground states are the correct ones in the magic number cases .
however , we fail to observe commensurability effects in the mc simulations .
the magic numbers are a property of the ground state , and do not affect the vortex liquid over the temperature range investigated .
oneill and moore have investigated the correlation function associated with vortex positions.@xcite this method is unwieldy as it requires finding the roots of an @xmath6th order polynomial at every measurement .
it is more convenient to look at the four point density correlation function ( briefly investigated by lee and moore@xcite and related to the correlator studied on the qp systems@xcite ) @xmath171 this is most revealing in reciprocal space . on a spherical surface
this means that we need to decompose the function into spherical harmonics : @xmath172 within the lll subspace , this correlator can be written as thermal averages over the coefficients @xmath173 : @xmath174,\end{aligned}\ ] ] with @xmath175 given in appendix  [ ap : struc ] and @xmath176 equal to the mean density @xmath177 . to make comparison with different temperatures we look at the function @xmath178 and
then normalize this function by its infinite temperature limit @xmath179 ( also given in appendix  [ ap : struc ] ) .
we average our measurements over @xmath180 , so our investigations center on the function @xmath181 $ ] with @xmath182 .
@xmath183 is the same correlator that has been looked at in a plane geometry by hu and macdonald@xcite as well as by yeo and moore,@xcite but modified for a sphere .
examples of this function for different @xmath55 and @xmath6 are shown in fig .
[ fig:1 ] .
we have concentrated our measurements on the form of the peak in @xmath184 at the first reciprocal lattice vector . in our units for a triangular lattice with a unit cell of area @xmath185 this is given by @xmath186 .
in the entire temperature range we look at ( @xmath187 ) we find that a lorentzian always gives an excellent fit to this peak ( see fig .
[ fig:1 ] ) .
this is consistent with an exponential decrease of the density correlations in real space .
the associated correlation length @xmath188 is easily extracted from the width @xmath189 of the peak at half its maximum , as @xmath190 . in fig .
[ fig:2 ] we plot @xmath191 against @xmath55 for different system sizes , which shows how the correlation length grows as the temperature is lowered , and reveals when finite - size effects become important . down to @xmath192 the results for the larger systems appear to be independent of system size . for @xmath193 , finite size effects for @xmath194 make its values of @xmath195 fall beneath the trend of the larger systems .
a linear fit at low temperatures is at least consistent with the results from vortex - vortex correlations@xcite that @xmath196 , as @xmath197 .    for small @xmath198 corrections to this scaling form are apparent , so that the range of this fit is small .
nevertheless , it is clear that there is no jump in the value of @xmath188 in the vicinity of the observed transition of qp simulations , or the predicted kthny transition . for the integral of @xmath184 under a given peak to remain constant
we would expect the general form @xmath199 in the vicinity of the peak.@xcite this implies that the peak height is proportional to the inverse width . using the parameters from our lorentzian fits ,
this relation is demonstrated in the inset of fig .
[ fig:1 ] .
= 8 cm     +    = 8 cm     +    these results should be compared to those of yeo and moore where @xmath184 was calculated within a high temperature approximation that uses a non - perturbative method to go to lower temperatures.@xcite in this work similar zero temperature scaling is seen down to @xmath192 , although the peak heights do not grow as fast .
in contrast , hu and macdonald,@xcite who calculate @xmath184 with an mc simulation on a plane with periodic boundary conditions , see an order of magnitude jump in the peaks at a transition temperature of @xmath200 . below this first order transition
, @xmath184 shows `` bragg peaks '' at the reciprocal lattice vectors .
however , in the temperature range above this freezing , the correlation function appears to be quantitatively the same as on the sphere at the same temperature .
this is as it should be as the liquid phase on the sphere has the same properties as the liquid on the flat plane as long as the system size is sufficiently large .
it is long range coherence in the phase of the order parameter ( odlro ) that gives rise to superconductivity . within the lll approximation
the phase can only be changed by the movement of the vortices , so one might expect a simple relation between vortex correlations and phase coherence .
sasik , stroud and tesanovic@xcite have defined a phase correlation function as @xmath201 . from measurements of this in their mc simulation with qp boundary conditions they found that in the thermodynamic limit @xmath202 fell to zero over the order of a lattice constant even when in the solid phase characterized by long range density correlations .    as our model is always in the liquid state
it is clear that @xmath202 will always be trivially short ranged simply due to the translations of the vortex system .
however we are interested in the phase correlations between points moving in the same patch of correlated vortices .
it is the length scale associated with this type of correlation that is physically relevant .
that is , when this length scale approaches the limiting scale of a given sample ( eg . the distance between pinning centers ) it will become superconducting . to investigate this length scale we need to measure the correlations in the phase from a fixed vortex .
we also need to fix the rotations about this vortex ( these would also destroy phase coherence ) . in practise
we can do this by ( a ) fixing one vortex at a pole ( @xmath149 ) by setting @xmath203 , and ( b ) fixing the nearest six vortices by keeping @xmath204 for @xmath205 .
@xmath206 are the coefficients for the ground state with a fixed orientation .
these conditions effectively fix the order parameter in a local region around the central fixed vortex .
we can then measure the correlation between the phase at a point @xmath207 and at @xmath208 : latexmath:[\[\sigma_{fix}(\hbox{\bf r}')= \langle \psi^*(\hbox{\bf 0 } ) \psi(\hbox{\bf r}')/    the difference between this and the correlation function of sst . as well as the fixing of the vortices
, our definition contains no part dependent on the amplitude of the order parameter . in terms of the mc variables ,
the phase at @xmath149 is given by @xmath210 where @xmath211 is the azimuthal angle of approach as @xmath212 .
fixing this angle @xmath211 to be zero , the phase correlation function may be written as @xmath213    the reader may notice a problem in measuring the phase of a complex function on the surface of a sphere .
if a path is followed that encircles one vortex , the phase of @xmath1 will change by @xmath214 .
however , an ambiguity exists in the definition of the inside or the outside of a closed loop in a closed surface and this could lead to an ill defined phase . in fact , it is our choice of gauge that places on one axis a string ( or infinitely long and thin solenoid ) terminating at the monopole that removes this ambiguity .
this choice effectively makes the south pole ( @xmath215 ) the `` outside '' of the system , and the north pole ( @xmath149 ) the center .
this results in a singularity in the phase of the order parameter at @xmath215 .
encircling the south pole    = 17.5 cm    2 without apparently encircling any vortex results in a phase change of @xmath216 ( this is the origin of the dirac quantization of the charge on a monopole ) .
this does not imply a singularity in the supercurrent , as the gauge part of the current will exactly cancel with this wave function contribution .
our measurements of @xmath217 are shown in fig .
[ fig:4 ] using the projection @xmath218 . for comparison
, we also show the mean density @xmath219 with the same coefficients fixed
. a growing range of phase coherence as the temperature is lowered is demonstrated , and coincides with the growth in density correlations .
the length scale @xmath220 is found by plotting @xmath221 against @xmath222 and fitting a gaussian over the region outside the fixed vortices , i.e.  for @xmath223 with @xmath224 .
we find that @xmath220 increases proportionally to @xmath198 in the temperature range @xmath225 ( see fig .  [ fig:6 ] ) , with again corrections to scaling visible for @xmath226 .
the similarity between figures  [ fig:2 ]  and  [ fig:6 ] indicates that the length scales for phase and density coherence grow proportionally to each other as the temperature is reduced .
the length scales are unchanged within the measured errors for the different system sizes studied . a relation between the growth of correlations in the density @xmath227 and in the phase of @xmath1
is clearly observed .
the scaling of phase coherence @xmath228 was predicted by oneill and moore from two distinct arguments .
one of these involved considering the thermal excitation of small amplitude shear modes from the ground state lattice to define a phase coherence length within the crystalline phase .
the second argument considered the free energy cost of phase fluctuations over a finite region of the vortex system .
it is this second argument that is directly relevant to the liquid phase with a finite length scale of crystalline order , and we repeat and amplify it here .
the free energy cost of phase fluctuations over a region with crystalline correlations was shown to be @xmath229 , where @xmath230 is the phase of the order parameter . from this result ,
the free energy cost of a phase change of @xmath231 over a region of linear extent @xmath232 will be @xmath233 .
when such excitations proliferate , phase coherence will be destroyed .
hence , there will be phase coherence only over the length scale @xmath232 when the free energy @xmath234 is of order @xmath235 .
the @xmath55 dependence of the mean - field elastic shear modulus @xmath236 is known in the large field limit , @xmath237 , with the simple result that @xmath238 .
= 8 cm    we believe that there is another kind of phase fluctuation within the liquid that can be identified with a topological defect a `` braid '' .
consider a ring of vortices of a given radius that move around the circumference of the ring .
( in three dimensions this vortex movement may take place along the field direction to form a topologically constrained defect hence the name , braid .
we have previously calculated the energies of such defects.@xcite ) such a movement may be created by a simple combination of changes to the basis coefficients . if we make the transformation @xmath239 for @xmath240 and make no change to the remaining coefficients @xmath241 then such a braid movement will occur in a ring of vortices that enclose approximately @xmath242 vortices ( see fig .
[ fig : braid ] ) .
= 6 cm    if @xmath243 the vortices in the ring will each hop to a neighbors position , with a phase change in the order parameter of @xmath214 inside the ring ( i.e.  no essential change ) .
excitations will involve values of @xmath114 less than @xmath214 which increase the phase at the center of the ring by @xmath114 while leaving the phase unchanged outside of the ring .
the energies of these braid excitations will define a length scale over which they proliferate , and thus destroy phase coherence .
the energy costs of braids in a crystal phase would increase with the size of the ring , and so phase coherence would not be destroyed by this mechanism in the crystalline state .
nevertheless , we believe that the excitations of such defects are present on the length scale @xmath232 in the liquid phase , and are a plausible mechanism for the loss of phase coherence seen in fig .
[ fig:4 ] .
notice that in fig .
[ fig:4 ] it is when the density correlations become ring - like ( as would happen due to braid excitations ) that the phase coherence is lost .
if an experiment using the spherical geometry was ever to be realized , there would always be a problem in placing the end of a solenoid ( the `` monopole '' ) exactly at the center of the sphere .
this will cause the external field to be inhomogeneous over the surface of the sphere .
such inhomogeneities are generic in experiments with any geometry , so it is worthwhile investigating these effects .
we consider displacing the monopole from the center of the sphere by a distance @xmath4 along the @xmath244-axis .
it is straightforward to show that for small @xmath245 , this induces a change in the vector potential of @xmath246 .
an inhomogeneity is introduced into the magnetic field at the surface which has a radial component of @xmath247 .
substituting this new field into the free energy gives an extra term to first order in @xmath248 of @xmath249 where @xmath250 .
before we look at our mc results with this extra term , we consider its effect within the high - temperature perturbation theory . with the addition of ( [ eq : extra ] ) to the hamiltonian ,
the hf equation becomes @xmath251 with @xmath252 .
the sum in this equation is no longer independent of @xmath180 .
however , in the large @xmath6 limit , @xmath253 has the form , @xmath254,\ ] ] where @xmath255^{1/2}$ ] .
the width of the interaction is of order @xmath256 over most of the sphere .
if the variation in the renormalized propagator @xmath257 is smooth over the range of @xmath256 ( which corresponds to changes in real space being smooth over the distance @xmath258 ) then the sum over @xmath259 is dominated by contributions when @xmath260 to give @xmath261 .
therefore we have @xmath18 independent hf equations @xmath262 .
each has the same solution , but with an @xmath180-dependent effective temperature : @xmath263 $ ] with @xmath264 . this will be true for any change in the hamiltonian of the form @xmath265 . where the effective temperature becomes @xmath266 .
such a change corresponds to changes in real space in one direction only over length scales longer than @xmath258 .
this procedure can be generalized to an arbitrary diagram in the perturbative scheme .
any diagram is dominated by contributions where the indices in and out of a vertex are within @xmath267 of each other .
therefore a given diagram is just the sum of the respective diagrams without the inhomogeneity but at the effective temperature @xmath268 .
for example , the free energy may be expanded in a power series of the propagator , @xmath269 . with the inhomogeneity this generalizes to @xmath270
this result may break down when variation in @xmath271 is on a scale less than @xmath256 . in the thermodynamic
limit our conclusion will hold as long as the perturbation expansion remains useful .
= 8 cm     +    interesting effects may occur when there is a phase transition in the pure system .
however , as our simulations fail to see such a transition , the above analysis is sufficient to explain our numerical results . in fig .
[ fig : homvs ] we plot the @xmath180-dependent propagator @xmath272 from mc simulations with @xmath273 . and compare to the zero inhomogeneity propagator with the @xmath180-dependent effective temperature .
the results confirm that the system behaves with the local properties of the effective temperature .
in this section we discuss some of the problems of different boundary conditions when interpreting the results of numerical simulations of superconducting thin films within the lll approximation .
we first mention some of the shortcomings of the spherical geometry used in our simulations .
we then go on to discuss the other main geometry used qp boundary conditions . in this case there exist problems not present in the spherical geometry due to the constraints on the phase of the order parameter .
we also consider a third geometry with a circular geometry , where initial work has not found a phase transition . finally , we look at the implications of recent experiments .
an important difficulty on the sphere is that the crystalline state free energy is raised due to the curvature and the necessary presence of twelve disclinations .
a liquid with correlation length less than the system size will be unaffected by the curvature .
the energy difference should disappear in the thermodynamic limit ( see section  [ sec : ground ] ) .
however , for the system sizes studied , this energy remains on the same order of magnitude of the jump in energy seen at the transition observed with qp boundary conditions .
the presence of the topologically required disclinations on the sphere is clearly an important feature not present in real thin films .
however , their effects will only become important as the correlation length of the system approaches the distance between disclinations .
this is a finite size effect as experienced by any geometry , although the effective system size will be smaller for a given @xmath6 than on a plane . for the largest system sizes studied ,
finite - size effects probably become important in the region of interest at @xmath274 .    the mobility of these disclinations , if large , would be a mechanism for destroying crystalline order in a finite - sized system . at low temperatures the disclinations should not be mobile , both because of their mutual repulsion ( as can be seen from the elasticity calculation ) and the topological barrier to their motion through the lattice . at higher temperatures where there is a finite density of dislocations ( eg . in a hexatic phase ) the dislocations provide a topological mechanism for the motion of the disclinations . in finite size systems , the strain field associated with the disclination may even create dislocations by pulling apart the dislocation pairs that are always present in a 2d crystal . from inspection of snapshots of the systems at finite temperatures we have observed that dislocations are attracted to the disclinations , and that the positions of the disclinations are not perfectly icosahedral as they are in the ground states .
our measurements so far have been insensitive to hexatic bond - angle order .
this order could never be globally defined on the sphere , and will only be a meaningful concept within the regions free of any disclinations , which are small in the system sizes we study .
it must be conceded , therefore , that our numerical evidence from the sphere for the non - existence of a phase - transition is not yet conclusive .
larger system size studies are clearly needed , but even increasing sizes by an order of magnitude may not resolve the problem .
our work on the ground states suggests that very large systems with @xmath275 may be needed to reach certain properties of the thermodynamic limit such as the screening of the disclinations by dislocations .      in the system investigated by refs .   a constant magnetic field * h * perpendicular to a thin - film superconductor
is represented by the landau gauge , @xmath276 . in this gauge
the lll states are @xmath277}$ ] , where @xmath278 is a `` momentum '' in the @xmath279 direction that must be fixed by the boundary conditions . for a general state in the lll subspace ,
@xmath280 , the system is periodic in the x direction with period @xmath281 .
quasi - periodicity in the y - direction may also be imposed by setting @xmath282 .
this condition leads to @xmath283 with @xmath284 .
so each principle region has an area @xmath285 and contains @xmath6 zeros of @xmath1 that correspond to vortices . to have a geometry commensurate with a triangular vortex lattice , @xmath286 and @xmath287
are fixed as integers .
( note the difference between the magnetic length @xmath258 and the lattice spacing @xmath26 both of which are defined in section  [ sec : method ] . )
the difference between this system and more conventional periodic boundary conditions comes from the constraints on the phase of the order parameter as the principal region is crossed . that such a constraint exists is apparent from the fact that there are only @xmath288 degrees of freedom available to describe the @xmath288 coordinates of the vortex positions as well as an overall phase and amplitude .
it has been shown that the phase constraints lead to the center of mass of the vortex positions to be fixed,@xcite which accounts for the missing two parameters .
a consequence of the fixed center of mass is that there are only @xmath6 distinct ground states , each a triangular lattice separated by a minimum discrete distance .    in mc simulations with qp boundary conditions ,
evidence is seen for a first order transition .
the observed melting temperature drops as system size is increased but extrapolation gives an estimate for the thermodynamic transition@xcite at @xmath289 . in the low temperature state dominated by the ground state there is a finite energy barrier between different ground states .
we now estimate what that energy barrier is .
= 8 cm     +    figure  [ fig : qp9 ] illustrates the sort of movement of vortices between ground states that obeys the center of mass constraint .
the example is for @xmath290 , and a short inspection will reveal that there are only 9 distinct ground states that may be related by such moves . in general the different ground states have the form @xmath291}$ ] for @xmath259 equal to an integer multiple of @xmath292 , and zero otherwise .
each ground state is labeled by @xmath293 , where @xmath294 must be an integer from @xmath169 to @xmath292 , and the condition @xmath295 quantizes the allowed values of @xmath296 to integers from @xmath169 to @xmath297 .
these ground states are exactly equivalent to the eilenberger@xcite states @xmath298 , with @xmath299 where @xmath300 and @xmath301 .
these are related by @xmath302 } \phi(\hbox{\bf r}+\hbox{\bf r}_0| \hbox{\bf 0})$ ] .
consider the change @xmath303 for any configuration of @xmath304 .
for an infinite system this just corresponds to a uniform translation of the vortices in the x direction by @xmath305 .
however , with qp boundary conditions this change will involve movement of vortices at the edges of the principal region in the opposite direction , thus conserving the center of mass .
( an analogous change in the spherical geometry gives the braid in section  [ sec : phase ] . )
it is exactly this sort of movement that is a path between different ground states .
for example , if we start from @xmath306 and apply this transformation with @xmath114 changing from @xmath59 to @xmath307 we arrive at the state @xmath308 .    what is the energy barrier for such a path ?
the free energy of the qp system is given by @xmath309   , \end{aligned}\ ] ] with @xmath310\sqrt{3}\pi/4n_x^2\}}$ ] .
changing the phase of the coefficients does not affect the quadratic energy .
also , the contributions to the quartic term @xmath311 will be zero when @xmath259 , @xmath312 , @xmath313 , and @xmath87 correspond to the same region ( ie .
they are all between @xmath59 and @xmath314 ) . as the width of the gaussian in the quartic interaction is of order @xmath292 , contributions will be negligible unless these variables are within @xmath315 of each other , and the edge of the principal region . therefore the energy barrier of this path
will be of order @xmath256
. this @xmath6 dependence has been observed numerically.@xcite    = 8 cm     +    a more quantitative estimate may be obtained by considering the movement of the vortices at the edge of the principal region relative to the vortices in the bulk of this region ( see fig  [ fig : inf ] ) . at any point on the path this represents a defect from the perfect ground state , the energy cost of which may be calculated using the methods of ref .  .
we find an energy barrier of @xmath316 which is proportional to @xmath292 as predicted .
this is an upper bound , as only the movement of a line of vortices has been allowed , without relaxation of the surrounding lattice .
however , as the interaction energy is short ranged,@xcite we do not expect the @xmath6 dependence to change .
it is also a zero temperature result , but it should still hold for a low temperature state with long correlations . we emphasize that this energy barrier grows with system size , despite the fact that the positions of the discrete ground states get closer to each other .
additional evidence for the size of these barriers comes from the simulations themselves . in ref .
a plot is shown of @xmath317 at a temperature corresponding to the extrapolated thermodynamic transition , where the modulations in this quantity show that the system is dominated by just one of the @xmath6 ground states during the length of the measurement , @xmath318 mc steps .
the measurements by sst of this low temperature state are clearly not equilibrated as not all of the ground states are being sampled ; there should always be a translationally invariant density for any _ finite _ system .
it is also a possibility that energy barriers between ground states affect the evidence for a first order phase transition .
for instance , the free energy of the liquid will be overestimated if the simulation is not run for long enough ( the standard time for these simulations has been @xmath319 mc steps ) which may alter the melting temperature .
this casts doubt on the extrapolation of the melting temperature to the thermodynamic limit .
even if a first order ( or continuous ) transition remains after full equilibration , the discrete ground state problem makes the qp system unsuitable for transport measurements .
an alternative set of boundary conditions may be constructed in the symmetric gauge,@xcite @xmath320 .
the lll states may be formed by the complete basis with circular symmetry @xmath321 $ ] , where @xmath322 . to perform numerical simulations
, one has to truncate the series of these states to a finite number of terms , say , @xmath323 to @xmath6 .
this truncated system has @xmath6 vortices , as @xmath324}\prod_{i=1}^n   ( \zeta-\zeta_i)$ ] , and the order parameter drops rapidly to zero outside of a radius @xmath325 from the origin .
little work has been done on this geometry as finite size effects are strong .
numerical minimizations on small systems have shown that the ground state energy per vortex is _ lower _ than for the infinite plane,@xcite and vortices near the edge of the system are attracted to the `` boundary''@xcite ( near where the order parameter is zero ) .
the finite size effects are greater than with spherical or qp boundary conditions where there are no edges of the finite systems and each vortex resides in an equivalent local environment .
the effect of the boundary is mainly in a fraction of @xmath326 of the bulk of the system , and should be unimportant in the thermodynamic limit .
a similar geometry has been used in simulations of the 2d one component plasma ( ocp ) , where the particles are confined to a circular disk.@xcite a first order transition was observed and this contrasts with the absence of a transition in simulations of the ocp on a spherical geometry.@xcite this has motivated us to perform simulations on the circular geometry in search of a phase transition . we have performed the real space phase and density correlation measurements as in section  [ sec : phase ] for systems up to @xmath51 and find results consistent with those on the sphere .
we have also measured the internal energy while heating from the ground state and cooling from high temperatures , for system sizes up to @xmath327 .
these energy measurements ( see fig  [ fig : circ ] ) have revealed no hysteresis on the scale observed in the qp system for temperatures greater than @xmath328 , and hence no evidence for a phase transition .
= 8 cm     +      it seems that from a theoretical point of view there are three possible scenarios for the behavior of thin film superconductors in the presence of a magnetic field : a ) the standard picture of kthny melting of a 2d crystal where a dislocation mediated continuous transition occurs between a low temperature state with quasi - crystalline order and a high temperature liquid or hexatic phase ; b ) a first order transition where the liquid free energy becomes lower than the solid phase before kthny melting takes place ( as is seen in qp simulations ) ; c ) the absence of a distinct low - temperature state , but where the correlation lengths of the liquid diverge as @xmath329 ( this is the prediction of high order perturbation theory from the high temperature side , and is consistent with the numerical simulations on the sphere ) .
we now discuss some experimental results and their implications for these possibilities .
a melting transition of the flux lattice in in / in0@xmath330 thin films was observed by gammel _
et al._@xcite the zero field limit of the transition was consistent with the kosterlitz - thouless transition in the absence of external fields , though the magnetic field dependence of the melting line did not satisfy the predictions of the kthny theory .
the measured transitions were mainly at low fields , so the relevance of our simulations here are doubtful .    in the experiments of berghuis _
et al._@xcite the current voltage response was measured of amorphous nb@xmath331ge@xmath330 films at different fields and temperatures .
a drop to zero in the a.c .
ohmic resistivity is seen , for example at fields @xmath332 at @xmath333 , though the d.c .
response remains finite ( ie .
the critical current is zero . ) .
this drop is attributed to a kthny transition consistent with the temperature and film thickness dependence of the melting field . a similar conclusion is made by yazdani _
et al._@xcite from a.c .
measurements on films of amorphous mo@xmath334ge@xmath335 .
it should be emphasized that these a.c .
measurements at a given frequency @xmath336 probe structural changes of the vortex system over a restricted length scale.@xcite while this is useful to measure properties independent of the presence of a large pinning length , one would expect even in the absence of a phase transition to see a large drop in the resistivity at fixed frequency when the relaxation times grow larger than @xmath337 . in numerical simulations on the sphere,@xcite the plastic time scale
was observed to increase as @xmath338 and we would expect a sharp drop when @xmath339 .
in contrast to the above experiments , measurements by nikulov _
et al._@xcite on films of amorphous nbo@xmath330 with weak pinning have found no signatures of a resistive transition down to fields @xmath340 .
this implies that the transitions observed previously are not universal properties of superconducting thin films .
while the lack of a resistive drop may be attributed to the total absence of pinning , one would expect in the presence of a first order transition , where there is a discontinuous change in @xmath341 , that this would be observed as a step in the flux flow resistivity .
thus the interpretation of these experiments remains controversial . while the kthny theory provides a useful framework to discuss flux lattice melting in two dimensions , its effects
have only been observed over short length scales.@xcite as far as we know there has been no experimental indication of a first order transition in superconducting thin films . in principle , the results of these experiments can still be explained simply by the growing length scales and relaxation times of a vortex liquid in the absence of a phase transition .
we have shown that the model with a spherical geometry remains in the vortex liquid phase , at least down to temperatures @xmath342 for system sizes up to @xmath343 .
this could be due to the raised energy of the ground states because of the necessary presence of twelve disclinations in a triangular network over a sphere .
however , we have made similar investigations on an alternative finite size model : considering the lll basis states on an infinite plane with a circular symmetric gauge and truncating the series to include only @xmath6 vortices.@xcite this system has a lower ground state energy per vortex than the infinite system and no disclinations , but we see the same length scales as in section  [ sec : phase ] for temperatures down to @xmath328 . the main evidence for a first order transition in the qp model is the double peak structure observed in the probability distribution of the internal energy at the transition .
our investigations on the internal energy for both the sphere and the plane with circular coordinates have revealed no such signatures .
the fact that the spherical system remains in the liquid phase down to low temperatures allows us to measure the variation of the correlations in the liquid state with temperature over a large regime .
measurements in the liquid phase have the advantages that the system reaches equilibrium over a reasonable time scale at low temperatures . our measurements of the density correlations are consistent with those in the liquid phase using qp boundary conditions.@xcite the extracted length scales demonstrate an asymptotic linear dependence on @xmath198 .
the scaling analysis of oneill and moore predicted that the _ phase coherence _
length should grow as @xmath344 .
our results for the phase coherence are consistent with this scaling form , and also demonstrate that in the liquid the phase and density correlations grow with identical temperature dependence .
our results for the scaling of the phase correlations contradict the conclusions of sst , who claim that in the thermodynamic limit @xmath73 the phase coherence has a range only of the order of the lattice spacing in their `` solid '' phase .
however , this discrepancy is really due to different definitions of what is meant by phase coherence .
sst do not predict the growth of the length scale of phase correlations as the temperature goes to zero .
although their results are claimed as evidence for the charge - density wave of tesanovic,@xcite a finite length scale that grows with decreasing temperature or even a power law form of odlro would still be consistent with this exotic phase .
it seems unlikely to us that the low - temperature phase below a first order transition would have a shorter range of phase coherence than the liquid phase .
if we measure the phase correlations without fixing any vortices , we also find no phase coherence beyond a length of one lattice spacing , due to the translational freedom of the liquid phase . in the solid phase of the qp model the vortex system appears to be fixed near one ground state due to the finite energy barrier between the @xmath6 ground states , and this allows phase coherence to be observed over a long range . as @xmath6 increases
, there is a decrease in the range of phase coherence .
these measurements by sst demonstrate just how strong the finite - size effects are in the qp simulations below the observed transition . without fixing the motion of any vortices ,
the phase coherence should always have a trivial length scale of the inter - vortex spacing . by restricting the motion of a small number of vortices
, we have been able to measure a more physically relevant phase correlation length that is consistent with the expected growth of order as the temperature is reduced .
an interesting question remains as to whether a similarly defined phase coherence in the solid phase with qp boundary conditions would result in true long range order , or another length scale that also exhibits zero temperature scaling .    despite the differences in results from spherical and qp boundary conditions
, our results from high temperature perturbation theory and on the ground states suggest that the thermodynamic limit will be the same on the sphere as the infinite flat plane .
this contradicts the earlier statement of ref .   where the effect of allowing additional dislocation defects into the ground state was not considered .
even so , the results on the ground states show how large system sizes must be for certain properties to approach the thermodynamic limit .
a definite drawback on the sphere is the small effective system size compared to the number of degrees of freedom .
the spherical geometry is perhaps not only relevant for simulational purposes .
it is possible that an experiment could be constructed with an amorphous type - ii superconductor sputtered over the surface of a sphere , and one end of a long and thin solenoid inserted into the center to approximate the field from a monopole .
in contrast an experimental realization of the qp boundary conditions is not possible ( these boundary conditions are primarily of interest in extrapolating to the thermodynamic limit of an infinite number of free vortices ) .
we feel that the spherical geometry , with its property of translational invariance , remains a useful system for numerical simulations , especially to measure properties of the unpinned vortex liquid .
an investigation into dynamical properties with a view to examining the length scales associated with non - local resistivity@xcite is planned , which would not be possible using qp boundary conditions due to the absence of full translational invariance of the vortex positions .
we would like to acknowledge useful interactions with gianni blatter , tom blum , vadim geshkenbein , roman sasik and joonhyun yeo .
mjwd would also like to thank han lee for the use of his source code and the epsrc for financial support .
in this appendix we will show that the gl free energy hamiltonian in eq .
( [ eq : ham ] ) depends only on the positions of the zeros in the order parameter on the sphere .
first , from eq .
( [ eq : prod ] ) we have @xmath345\nonumber\\ & = & d\prod_{i=1}^n[a_i\sin(\theta/2)e^{i\phi}-b_i\cos(\theta/2)]\nonumber\\ & = & d\prod_{i=1}^nf_i(\theta,\phi),\end{aligned}\ ] ] with @xmath346 and @xmath347 .
if @xmath348 when @xmath349 and @xmath350 then we may choose : @xmath351    now , the gl - lll free energy is a functional that depends only on the value of @xmath352 at each point on the sphere .
if we can show that @xmath353 depends only on the relative position of the @xmath152th zero then the hamiltonian will be invariant with respect to coordinate rotations of all the zeros : @xmath354\nonumber\\ & = & \frac{1}{2}[1-\cos{\theta'(i)}]\nonumber\\ & = & \sin^2[\theta'(i)/2]\end{aligned}\ ] ] where @xmath355 is the polar angle measured from the axis through the @xmath152th zero to the point @xmath46 , defined by an analogous expression to eq .
[ eq : theta ] .
this is independent of the coordinate frame and so we have shown that @xmath356 depends only on the relative positions of the vortices plus the overall amplitude .
to find the reciprocal space density correlator from thermal averages of the lll basis state coefficients ( see eq .
( [ eq : recip ] ) ) we need the integrals @xmath357 which may be written : @xmath358}^{1/2}\times\nonumber\\ & & \hspace{-1 cm } \times\frac{4\pi p!(n - p+m)!}{(n+m+1)!m ! } \,_3f_2\left [ \begin{array}{l }   m - l , m+l+1,p+1 \\ m+1,n+m+2   \end{array}\right ] , \end{aligned}\ ] ] where @xmath359 $ ] is a generalized hypergeometric series . to find the infinite temperature limit of @xmath360 we consider the case @xmath323 and calculate the averages @xmath361 and @xmath362 exactly in the limit @xmath54 ( this involves only gaussian integrals ) .
our result is @xmath363}^2 \nonumber\\ & = & \frac{(n!)^2}{(n - l)!(n+l+1)!}.\label{eq : infinite}\end{aligned}\ ] ] ( we are indebted to t. blum for this final result . ) this is consistent with the fact that in the limit @xmath364 , the infinite temperature correlation function should be identical with that for an infinite plane,@xcite and the expression in eq .
( [ eq : infinite ] ) should tend to @xmath365 .
999 d.  s.  fisher , phys .
b * 22 * , 1190 ( 1980 ) ; b.  a.  huberman and s.  doniach , phys .  rev
.  lett . * 43 * , 950 ( 1979 ) .
p. l. gammel , a. f. hebard and d. j. bishop , phys.rev .
* 60 * , 144 ( 1988 ) .
p.  berghuis and p.  h.  kes , phys
.  rev .
b * 47 * , 262 ( 1993 ) ; p.  berghuis , a.  l.  f.  van der slot and p.  h.  kes , phys .
65 * , 2583 ( 1990 ) .
a.  yazdani , w.  r.  white , m.  r.  hahn , m.  gabay , m.  r.  beasley and a.  kapitulnik , phys .
.  lett . * 70 * , 505 ( 1993 ) .
a.  v.  nikulov , d.  y.  remisov and v.  a.  oboznov , phys .
lett . * 75 * , 2586 ( 1995 ) .
m.  a.  moore , phys .
b * 39 * , 136 ( 1989 ) ; * 45 * , 7336 ( 1992 ) . j.  a.  oneill and m.  a.  moore , phys .  rev
.  lett . * 69 * , 2582 ( 1992 ) ; phys .
b * 48 * , 374 ( 1993 ) .
p.  a.  m.  dirac , proc .
r.  soc .
london , ser .
a * 133 * , 60 ( 1931 ) .
j.  pearl , appl .  phys .
lett . * 5 * , 65 ( 1964 ) .
z.  tesanovic and l.  xing , phys .
67 * , 22729 ( 1991 ) ; l.  xing and z.  tesanovic , physica c * 196 * , 241 ( 1992 ) . y.  kato and n.  nagaosa , phys .  rev .
b * 47 * , 2932 ( 1993 ) ; * 48 * , 7383 ( 1993 ) .
j.  hu and a.  h.  macdonald , phys .
rev .  lett . * 71 * , 432 ( 1993 ) ; phys .
b * 49 * , 15263 ( 1994 ) .
r.  sasik and d.  stroud , phys .
b * 49 * , 16074 ( 1993 ) .
r.  sasik , d.  stroud and z.  tesanovic , phys .
b * 51 * , 3042 ( 1995 ) .
d. r. nelson , in _ phase transitions and critical phenomena _ , edited by c. domb and j. l. lebowitz ( academic press , new york , 1983 ) , vol .  7 , p.  5 .
z.  tesanovic , physica c * 220 * , 303 ( 1994 ) .
m.  j.  w.  dodgson , j. phys .
a : math . gen . * 29 * , 2499 ( 1996 ) .
s.  m.  roy and v.  singh , phys .
lett . * 51 * , 2069 ( 1983 ) .
h.  h.  lee and m.  a.  moore , phys .
b * 49 * , 9240 ( 1994 ) .
n.  metropolis , a.  w.  rosenbluth , m.  n.  rosenbluth , a.  h.  teller and e.  teller , j.   chem .  phys . * 21 * , 1087 ( 1953 ) .
g. j. ruggeri and d. j. thouless , j. phys .
f : metal phys .
, * 6 * , 2063 ( 1976 ) . j. hu , a. h. macdonald and b. d. mckay , phys .
b * 49 * , 15263 ( 1994 ) ; e. brezin , a. fujita and a. hikami , phys . rev . lett . *
65 * , 1949 ( 1990 ) .
g. p. mccauley and d. j. thouless , j. phys .
f : metal phys . , * 6 * , 109 ( 1976 ) . w.  h.  kleiner , l.  m.  roth and s.  h.  autler , phys
.  rev . *
133 * , 1226 ( 1964 ) .
h.  s.  seung and d.  r.  nelson , phys .
a * 38 * , 1005 ( 1988 ) .
f. r. n. nabarro , _ theory of crystal dislocations _
( clarendon press , oxford , 1967 ) , pp .  53 - 7
. h. goldstein , _ classical mechanics _ ( addison - wesley , 1980 ) , pp .
j.  yeo and m.  a.  moore , phys .
lett . * 76 * , 1142 ( 1996 ) .
m.  j.  w.  dodgson and m.  a.  moore , phys .
b * 51 * , 11887 ( 1995 ) .
f. d. m. haldane and e. h. rezayi , phys .
b * 31 * , 2529 ( 1985 ) .
g. eilenberger , phys .
* 164 * , 628 ( 1967 ) .
r. sasik ( private communications ) .
j. a. oneill , phd thesis ( university of manchester ) .
choquard and j. clerouin , phys .
lett . * 50 * , 2086 ( 1983 ) .
j. m. caillol , d. levesque , j. j. weis and j. p. hansen , jour .
* 28 * , 325 ( 1982 ) .
a. k. kienappel and m. a. moore ( private communication ) .