On the finite-size behavior of systems with asymptotically large critical shift

Exact results of the finite-size behavior of the susceptibility in three-dimensional mean spherical model films under Dirichlet-Dirichlet, Dirichlet-Neumann and Neumann-Neumann boundary conditions are presented. The corresponding scaling functions are explicitly derived and their asymptotics close to, above and below the bulk critical temperature $T_c$ are obtained. The results can be incorporated in the framework of the finite-size scaling theory where the exponent $\lambda$ characterizing the shift of the finite-size critical temperature with respect to $T_c$ is smaller than $1/\nu$, with $\nu$ being the critical exponent of the bulk correlation length.Comment: 24 pages, late

Monte Carlo simulation of ice models

We propose a number of Monte Carlo algorithms for the simulation of ice models and compare their efficiency. One of them, a cluster algorithm for the equivalent three colour model, appears to have a dynamic exponent close to zero, making it particularly useful for simulations of critical ice models. We have performed extensive simulations using our algorithms to determine a number of critical exponents for the square ice and F models.Comment: 32 pages including 15 postscript figures, typeset in LaTeX2e using the Elsevier macro package elsart.cl

Nature of Phase Transitions of Superconducting Wire Networks in a Magnetic Field

We study $I$-$V$ characteristics of periodic square Nb wire networks as a function of temperature in a transverse magnetic field, with a focus on three fillings 2/5, 1/2, and 0.618 that represent very different levels of incommensurability. For all three fillings, a scaling behavior of $I$-$V$ characteristics is found, suggesting a finite temperature continuous superconducting phase transition. The low-temperature $I$-$V$ characteristics are found to have an exponential form, indicative of the domain-wall excitations.Comment: 5 pages, also available at http://www.neci.nj.nec.com/homepages/tang.htm

Correlated percolation and the correlated resistor network

We present some exact results on percolation properties of the Ising model, when the range of the percolating bonds is larger than nearest-neighbors. We show that for a percolation range to next-nearest neighbors the percolation threshold Tp is still equal to the Ising critical temperature Tc, and present the phase diagram for this type of percolation. In addition, we present Monte Carlo calculations of the finite size behavior of the correlated resistor network defined on the Ising model. The thermal exponent t of the conductivity that follows from it is found to be t = 0.2000 +- 0.0007. We observe no corrections to scaling in its finite size behavior.Comment: 16 pages, REVTeX, 6 figures include

Finite-size scaling and conformal anomaly of the Ising model in curved space

We study the finite-size scaling of the free energy of the Ising model on lattices with the topology of the tetrahedron and the octahedron. Our construction allows to perform changes in the length scale of the model without altering the distribution of the curvature in the space. We show that the subleading contribution to the free energy follows a logarithmic dependence, in agreement with the conformal field theory prediction. The conformal anomaly is given by the sum of the contributions computed at each of the conical singularities of the space, except when perfect order of the spins is precluded by frustration in the model.Comment: 4 pages, 4 Postscript figure

Renormalization Group Study of the Intrinsic Finite Size Effect in 2D Superconductors

Vortices in a thin-film superconductor interact logarithmically out to a distance on the order of the two-dimensional (2D) magnetic penetration depth $\lambda_\perp$, at which point the interaction approaches a constant. Thus, because of the finite $\lambda_\perp$, the system exhibits what amounts to an {\it intrinsic} finite size effect. It is not described by the 2D Coulomb gas but rather by the 2D Yukawa gas (2DYG). To study the critical behavior of the 2DYG, we map the 2DYG to the massive sine-Gordon model and then perform a renormalization group study to derive the recursion relations and to verify that $\lambda_\perp$ is a relevant parameter. We solve the recursion relations to study important physical quantities for this system including the renormalized stiffness constant and the correlation length. We also address the effect of current on this system to explain why finite size effects are not more prevalent in experiments given that the 2D magnetic penetration depth is a relevant parameter.Comment: 8 pages inRevTex, 5 embedded EPS figure

Bloch-Wall Phase Transition in the Spherical Model

The temperature-induced second-order phase transition from Bloch to linear (Ising-like) domain walls in uniaxial ferromagnets is investigated for the model of D-component classical spin vectors in the limit D \to \infty. This exactly soluble model is equivalent to the standard spherical model in the homogeneous case, but deviates from it and is free from unphysical behavior in a general inhomogeneous situation. It is shown that the thermal fluctuations of the transverse magnetization in the wall (the Bloch-wall order parameter) result in the diminishing of the wall transition temperature T_B in comparison to its mean-field value, thus favouring the existence of linear walls. For finite values of T_B an additional anisotropy in the basis plane x,y is required; in purely uniaxial ferromagnets a domain wall behaves like a 2-dimensional system with a continuous spin symmetry and does not order into the Bloch one.Comment: 16 pages, 2 figure

Phase transitions in the antiferromagnetic XY model with a kagome lattice

The ground state of the antiferromagnetic XY model with a kagome lattice is characterized by a well developed accidental degeneracy. As a consequence the phase transition in this system consists in unbinding of pairs of fractional vortices. Addition of the next-to-nearest neighbors (NNN) interaction leads to stabilization of the long-range order in chirality (staggered chirality). We show that the phase transition, related with destruction of this long-range order, can happen as a separate phase transition below the temperature of the fractional vortex pairs unbinding only if the NNN coupling is extremely weak, and find how the temperature of this transition depends on coupling constants. We also demonstarte that the antiferromagnetic ordering of chiralities and, accordingly, the presence of the second phase transition are induced by the free energy of spin wave fluctuations even in absence of the NNN coupling.Comment: 10 pages (Revtex) + 8 figures (in 2 postscript files

Dynamic Critical Behavior of a Swendsen-Wang-Type Algorithm for the Ashkin-Teller Model

We study the dynamic critical behavior of a Swendsen-Wang-type algorithm for the Ashkin--Teller model. We find that the Li--Sokal bound on the autocorrelation time ($\tau_{{\rm int},{\cal E}} \ge {\rm const} \times C_H$) holds along the self-dual curve of the symmetric Ashkin--Teller model, and is almost but not quite sharp. The ratio $\tau_{{\rm int},{\cal E}} / C_H$ appears to tend to infinity either as a logarithm or as a small power ($0.05 \leq p \leq 0.12$). In an appendix we discuss the problem of extracting estimates of the exponential autocorrelation time.Comment: 59 pages including 3 figures, uuencoded g-compressed ps file. Postscript size = 799740 byte

Sine-Gordon mean field theory of a Coulomb Gas

Sine-Gordon field theory is used to investigate the phase diagram of a neutral Coulomb gas. A variational mean field free energy is constructed and the corresponding phase diagrams in two (2d) and three dimensions (3d) are obtained. When analyzed in terms of chemical potential, the Sine-Gordon theory predicts the phase diagram topologically identical with the Monte Carlo simulations and a recently developed Debye-H\"uckel-Bjerrum (DHBj) theory. In 2d we find that the infinite order Kosterlitz-Thouless line terminates in a tricritical point, after which the metal-insulator transition becomes first order. However, when the transformation from chemical potential to the density is made the whole of the insulating phase is mapped onto zero density.Comment: 5 pages, Revtex with twocolumn style, 2 Postscript figures. Submitted to PR