BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Alan D. Sokal appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
I give a combinatorial interpretation of the multiple Laguerre polynomials of the first kind of type II, generalizing the digraph model found by Foata and Strehl for the ordinary Laguerre polynomials. I also give an explicit integral representation for these polynomials, which shows that they form a multidimensional Stieltjes moment sequence whenever $x \le 0$.
We solve analytically the differential equations for a skier on a circular hill and for a particle on a loop-the-loop track when the hill or track is endowed with a coefficient of kinetic friction $μ$. For each problem, we determine the exact "phase diagram" in the two-dimensional parameter plane.
We find Stieltjes-type and Jacobi-type continued fractions for some "master polynomials" that enumerate permutations, set partitions or perfect matchings with a large (sometimes infinite) number of simultaneous statistics. Our results contain many previously obtained identities as special cases, providing a common refinement of all of them.
I show that a hypergeometric function ${}_{p}F_q(a_1,\ldots,a_p;b_1,\ldots,b_q;\,\cdot\,)$ with $p \le q$ belongs to the Laguerre--Pólya class $LP^+$ for arbitrarily large $b_{p+1},\ldots,b_q > 0$ if and only if, after a possible reordering, the differences $a_i - b_i$ are nonnegative integers. This result arises as an easy corollary of the case $p=q$ proven two decades ago by Ki and Kim. I also give explicit examples for the case ${}_{1}F_2$.
I give a short and completely elementary proof of Takagi's 1921 theorem on the zeros of a composite polynomial $f(d/dz) \, g(z)$.
We show that, for Hankel matrices, total nonnegativity (resp. total positivity) of order r is preserved by sum, Hadamard product, and Hadamard power with real exponent t \ge r-2. We give examples to show that our results are sharp relative to matrix size and structure (general, symmetric or Hankel). Some of these examples also resolve the Hadamard critical-exponent problem for totally positive and totally nonnegative matrices.
Two decades ago, Chauve, Dulucq and Guibert showed that the number of rooted trees on the vertex set $[n+1]$ in which exactly $k$ children of the root are lower-numbered than the root is $\binom{n}{k} \, n^{n-k}$. Here I give a simpler proof of this result.
We define an infinite sequence of generalizations, parametrized by an integer $m \ge 1$, of the Stieltjes--Rogers and Thron--Rogers polynomials; they arise as the power-series expansions of some branched continued fractions, and as the generating polynomials for $m$-Dyck and $m$-Schröder paths with height-dependent weights. We prove that all of these sequences of polynomials are coefficientwise Hankel-totally positive, jointly in all the (infinitely many) indeterminates. We then apply this theory to prove the coefficientwise Hankel-total positivity for combinatorially interesting sequences of polynomials. Enumeration of unlabeled ordered trees and forests gives rise to multivariate Fuss--Narayana polynomials and Fuss--Narayana symmetric functions. Enumeration of increasing (labeled) ordered trees and forests gives rise to multivariate Eulerian polynomials and Eulerian symmetric functions, which include the univariate $m$th-order Eulerian polynomials as specializations. We also find branched continued fractions for ratios of contiguous hypergeometric series ${}_r \! F_s$ for arbitrary $r$ and $s$, which generalize Gauss' continued fraction for ratios of contiguous ${}_2 \! F_1$; and for $s=0$ we prove the coefficientwise Hankel-total positivity. Finally, we extend the branched continued fractions to ratios of contiguous basic hypergeometric series ${}_r \! ϕ_s$.
We introduce the generic Lah polynomials $L_{n,k}(ϕ)$, which enumerate unordered forests of increasing ordered trees with a weight $ϕ_i$ for each vertex with $i$ children. We show that, if the weight sequence $ϕ$ is Toeplitz-totally positive, then the triangular array of generic Lah polynomials is totally positive and the sequence of row-generating polynomials $L_n(ϕ,y)$ is coefficientwise Hankel-totally positive. Upon specialization we obtain results for the Lah symmetric functions and multivariate Lah polynomials of positive and negative type. The multivariate Lah polynomials of positive type are also given by a branched continued fraction. Our proofs use mainly the method of production matrices; the production matrix is obtained by a bijection from ordered forests of increasing ordered trees to labeled partial Lukasiewicz paths. We also give a second proof of the continued fraction using the Euler--Gauss recurrence method.
We find a Thron-type continued fraction (T-fraction) for the ordinary generating function of the Ward polynomials, as well as for some generalizations employing a large (indeed infinite) family of independent indeterminates. Our proof is based on a bijection between super-augmented perfect matchings and labeled Schröder paths, which generalizes Flajolet's bijection between perfect matchings and labeled Dyck paths.
I give an elementary proof of Wall's continued-fraction characterization of Hausdorff moment sequences.
We consider two different interpretations of the Chu--Vandermonde identity: as an identity for polynomials, and as an identity for infinite matrices. Each interpretation leads to a class of possible generalizations, and in both cases we obtain a complete characterization of the solutions.
I study the sequences of Euler and Springer numbers from the point of view of the classical moment problem.
We prove the complete monotonicity on $(0,\infty)^n$ for suitable inverse powers of the spanning-tree polynomials of graphs and, more generally, of the basis generating polynomials of certain classes of matroids. This generalizes a result of Szego and answers, among other things, a long-standing question of Lewy and Askey concerning the positivity of Taylor coefficients for certain rational functions. Our proofs are based on two_ab initio_ methods for proving that $P^{-β}$ is completely monotone on a convex cone $C$: the determinantal method and the quadratic-form method. These methods are closely connected with harmonic analysis on Euclidean Jordan algebras (or equivalently on symmetric cones). We furthermore have a variety of constructions that, given such polynomials, can create other ones with the same property: among these are algebraic analogues of the matroid operations of deletion, contraction, direct sum, parallel connection, series connection and 2-sum. The complete monotonicity of $P^{-β}$ for some $β> 0$ can be viewed as a strong quantitative version of the half-plane property (Hurwitz stability) for $P$, and is also related to the Rayleigh property for matroids.
We find zero-free regions in the complex plane at large |q| for the multivariate Tutte polynomial (also known in statistical mechanics as the Potts-model partition function) Z_G(q,w) of a graph G with general complex edge weights w = {w_e}. This generalizes a result of Sokal (cond-mat/9904146) that applies only within the complex antiferromagnetic regime |1+w_e| \le 1. Our proof uses the polymer-gas representation of the multivariate Tutte polynomial together with the Penrose identity.
We introduce a new graph invariant Λ(G) that we call maxmaxflow, and put it in the context of some other well-known graph invariants, notably maximum degree and its relatives. We prove the equivalence of two "dual" definitions of maxmaxflow: one in terms of flows, the other in terms of cocycle bases. We then show how to bound the total number (or more generally, total weight) of various classes of subgraphs of G in terms of either maximum degree or maxmaxflow. Our results are motivated by a conjecture that the modulus of the roots of the chromatic polynomial of G can be bounded above by a function of Λ(G).
This is a response to five comments [American Psychologist 69, 626-629 and 632-635 (2014)] on our article arXiv:1307.7006.
This is a response to Barbara Fredrickson's comment [American Psychologist 68, 814-822 (2013)] on our article arXiv:1307.7006. We analyze critically the renewed claims made by Fredrickson (2013) concerning positivity ratios and "flourishing", and attempt to disentangle some conceptual confusions; we also address the alleged empirical evidence for nonlinear effects. We conclude that there is no evidence whatsoever for the existence of any "tipping points", and only weak evidence for the existence of any nonlinearity of any kind. Our original concern, that the application of advanced mathematical techniques in psychology and related disciplines may not always be appropriate, remains undiminished.
We prove the existence of long-range order at sufficiently low temperatures, including zero temperature, for the three-state Potts antiferromagnet on a class of quasi-transitive plane quadrangulations, including the diced lattice. More precisely, we show the existence of (at least) three infinite-volume Gibbs measures, which exhibit spontaneous magnetization in the sense that vertices in one sublattice have a higher probability to be in one state than in either of the other two states. For the special case of the diced lattice, we give a good rigorous lower bound on this probability, based on computer-assisted calculations that are not available for the other lattices.
We prove that the (real or complex) chromatic roots of a series-parallel graph with maxmaxflow Lambda lie in the disc |q-1| < (Lambda-1)/log 2. More generally, the same bound holds for the (real or complex) roots of the multivariate Tutte polynomial when the edge weights lie in the "real antiferromagnetic regime" -1 \le v_e \le 0. This result is within a factor 1/log 2 \approx 1.442695 of being sharp
We examine critically the claims made by Fredrickson and Losada (2005) concerning the construct known as the "positivity ratio". We find no theoretical or empirical justification for the use of differential equations drawn from fluid dynamics, a subfield of physics, to describe changes in human emotions over time; furthermore, we demonstrate that the purported application of these equations contains numerous fundamental conceptual and mathematical errors. The lack of relevance of these equations and their incorrect application lead us to conclude that Fredrickson and Losada's claim to have demonstrated the existence of a critical minimum positivity ratio of 2.9013 is entirely unfounded. More generally, we urge future researchers to exercise caution in the use of advanced mathematical tools such as nonlinear dynamics and in particular to verify that the elementary conditions for their valid application have been met.
We exhibit infinite families of two-dimensional lattices (some of which are triangulations or quadrangulations of the plane) on which the q-state Potts antiferromagnet has a finite-temperature phase transition at arbitrarily large values of q. This unexpected result is proven rigorously by using a Peierls argument to measure the entropic advantage of sublattice long-range order. Additional numerical data are obtained using transfer matrices, Monte Carlo simulation, and a high-precision graph-theoretic method.
The classic Cayley identity states that \det(\partial) (\det X)^s = s(s+1)...(s+n-1) (\det X)^{s-1} where X=(x_{ij}) is an n-by-n matrix of indeterminates and \partial=(\partial/\partial x_{ij}) is the corresponding matrix of partial derivatives. In this paper we present straightforward combinatorial proofs of a variety of Cayley-type identities, both old and new. The most powerful of these proofs employ Grassmann algebra (= exterior algebra) and Grassmann-Berezin integration. Among the new identities proven here are a pair of "diagonal-parametrized" Cayley identities, a pair of "Laplacian-parametrized" Cayley identities, and the "product-parametrized" and "border-parametrized" rectangular Cayley identities.
I study the leading root x_0(y) of the partial theta function Θ_0(x,y) = \sum_{n=0}^\infty x^n y^{n(n-1)/2}, considered as a formal power series. I prove that all the coefficients of -x_0(y) are strictly positive. Indeed, I prove the stronger results that all the coefficients of -1/x_0(y) after the constant term 1 are strictly negative, and all the coefficients of 1/x_0(y)^2 after the constant term 1 are strictly negative except for the vanishing coefficient of y^3.
Let R_αbe the Riesz distribution on a simple Euclidean Jordan algebra, parametrized by the complex number α. I give an elementary proof of the necessary and sufficient condition for R_αto be a locally finite complex measure (= complex Radon measure).
We study, via Monte Carlo simulation, the dynamic critical behavior of the Chayes-Machta dynamics for the Fortuin-Kasteleyn random-cluster model, which generalizes the Swendsen-Wang dynamics for the q-state Potts ferromagnet to non-integer q \ge 1. We consider spatial dimension d=2 and 1.25 \le q \le 4 in steps of 0.25, on lattices up to 1024^2, and obtain estimates for the dynamic critical exponent z_{CM}. We present evidence that when 1 \le q \lesssim 1.95 the Ossola-Sokal conjecture z_{CM} \ge β/νis violated, though we also present plausible fits compatible with this conjecture. We show that the Li-Sokal bound z_{CM} \ge α/νis close to being sharp over the entire range 1 \le q \le 4, but is probably non-sharp by a power. As a byproduct of our work, we also obtain evidence concerning the corrections to scaling in static observables.
We argue that the 4-state Potts antiferromagnet has a finite-temperature phase transition on any Eulerian plane triangulation in which one sublattice consists of vertices of degree 4. We furthermore predict the universality class of this transition. We then present transfer-matrix and Monte Carlo data confirming these predictions for the cases of the union-jack and bisected hexagonal lattices.
We study, using transfer-matrix methods, the partition-function zeros of the square-lattice q-state Potts antiferromagnet at zero temperature (= square-lattice chromatic polynomial) for the special boundary conditions that are obtained from an m x n grid with free boundary conditions by adjoining one new vertex adjacent to all the sites in the leftmost column and a second new vertex adjacent to all the sites in the rightmost column. We provide numerical evidence that the partition-function zeros are becoming dense everywhere in the complex q-plane outside the limiting curve B_\infty(sq) for this model with ordinary (e.g. free or cylindrical) boundary conditions. Despite this, the infinite-volume free energy is perfectly analytic in this region.
I give a proof of the uniform boundedness theorem that is elementary (i.e. does not use any version of the Baire category theorem) and also extremely simple.
We introduce several infinite families of new critical exponents for the random-cluster model and present scaling arguments relating them to the k-arm exponents. We then present Monte Carlo simulations confirming these predictions. These new exponents provide a convenient way to determine k-arm exponents from Monte Carlo simulations. An understanding of these exponents also leads to a radically improved implementation of the Sweeny Monte Carlo algorithm. In addition, our Monte Carlo data allow us to conjecture an exact expression for the shortest-path fractal dimension d_min in two dimensions: d_min = (g+2)(g+18)/(32g) where g is the Coulomb-gas coupling, related to the cluster fugacity q via q = 2 + 2 cos(gπ/2) with 2 \le g \le 4.
We prove, by simple manipulation of commutators, two noncommutative generalizations of the Cauchy-Binet formula for the determinant of a product. As special cases we obtain elementary proofs of the Capelli identity from classical invariant theory and of Turnbull's Capelli-type identities for symmetric and antisymmetric matrices.
We obtain a characterization of generalized Stieltjes functions of any order λ> 0 in terms of inequalities for their derivatives on (0,\infty). When λ=1, this provides a new and simple proof of a characterization of Stieltjes functions first obtained by Widder in 1938.
I reconsider the problem of a raindrop falling through mist, collecting mass, and generalize it to allow an arbitrary power-law form for the accretion rate. I show that the coupled differential equations can be solved by the simple trick of temporarily eliminating time (t) in favor of the raindrop's mass (m) as the independent variable
We derive some new structural results for the transfer matrix of square-lattice Potts models with free and cylindrical boundary conditions. In particular, we obtain explicit closed-form expressions for the dominant (at large |q|) diagonal entry in the transfer matrix, for arbitrary widths m, as the solution of a special one-dimensional polymer model. We also obtain the large-q expansion of the bulk and surface (resp. corner) free energies for the zero-temperature antiferromagnet (= chromatic polynomial) through order q^{-47} (resp. q^{-46}). Finally, we compute chromatic roots for strips of widths 9 <= m <= 12 with free boundary conditions and locate roughly the limiting curves.
We study the correction-to-scaling exponents for the two-dimensional self-avoiding walk, using a combination of series-extrapolation and Monte Carlo methods. We enumerate all self-avoiding walks up to 59 steps on the square lattice, and up to 40 steps on the triangular lattice, measuring the mean-square end-to-end distance, the mean-square radius of gyration and the mean-square distance of a monomer from the endpoints. The complete endpoint distribution is also calculated for self-avoiding walks up to 32 steps (square) and up to 22 steps (triangular). We also generate self-avoiding walks on the square lattice by Monte Carlo, using the pivot algorithm, obtaining the mean-square radii to ~0.01% accuracy up to N = 4000. We give compelling evidence that the first non-analytic correction term for two-dimensional self-avoiding walks is Delta_1 = 3/2. We compute several moments of the endpoint distribution function, finding good agreement with the field-theoretic predictions. Finally, we study a particular invariant ratio that can be shown, by conformal-field-theory arguments, to vanish asymptotically, and we find the cancellation of the leading analytic correction.
We study the q-state Potts model with nearest-neighbor coupling v=e^{βJ}-1 in the limit q,v \to 0 with the ratio w = v/q held fixed. Combinatorially, this limit gives rise to the generating polynomial of spanning forests; physically, it provides information about the Potts-model phase diagram in the neighborhood of (q,v) = (0,0). We have studied this model on the square and triangular lattices, using a transfer-matrix approach at both real and complex values of w. For both lattices, we have computed the symbolic transfer matrices for cylindrical strips of widths 2 \le L \le 10, as well as the limiting curves of partition-function zeros in the complex w-plane. For real w, we find two distinct phases separated by a transition point w=w_0, where w_0 = -1/4 (resp. w_0 = -0.1753 \pm 0.0002) for the square (resp. triangular) lattice. For w > w_0 we find a non-critical disordered phase, while for w < w_0 our results are compatible with a massless Berker-Kadanoff phase with conformal charge c = -2 and leading thermal scaling dimension x_{T,1} = 2 (marginal operator). At w = w_0 we find a "first-order critical point": the first derivative of the free energy is discontinuous at w_0, while the correlation length diverges as w \downarrow w_0 (and is infinite at w = w_0). The critical behavior at w = w_0 seems to be the same for both lattices and it differs from that of the Berker-Kadanoff phase: our results suggest that the conformal charge is c = -1, the leading thermal scaling dimension is x_{T,1} = 0, and the critical exponents are ν= 1/d = 1/2 and α= 1.
The multivariate Tutte polynomial (known to physicists as the Potts-model partition function) can be defined on an arbitrary finite graph G, or more generally on an arbitrary matroid M, and encodes much important combinatorial information about the graph (indeed, in the matroid case it encodes the full structure of the matroid). It contains as a special case the familiar two-variable Tutte polynomial -- and therefore also its one-variable specializations such as the chromatic polynomial, the flow polynomial and the reliability polynomial -- but is considerably more flexible. I begin by giving an introduction to all these problems, stressing the advantages of working with the multivariate version. I then discuss some questions concerning the complex zeros of the multivariate Tutte polynomial, along with their physical interpretations in statistical mechanics (in connection with the Yang--Lee approach to phase transitions) and electrical circuit theory. Along the way I mention numerous open problems. This survey is intended to be understandable to mathematicians with no prior knowledge of physics.
We elucidate the close connection between the repulsive lattice gas in equilibrium statistical mechanics and the Lovasz local lemma in probabilistic combinatorics. We show that the conclusion of the Lovasz local lemma holds for dependency graph G and probabilities {p_x} if and only if the independent-set polynomial for G is nonvanishing in the polydisc of radii {p_x}. Furthermore, we show that the usual proof of the Lovasz local lemma -- which provides a sufficient condition for this to occur -- corresponds to a simple inductive argument for the nonvanishing of the independent-set polynomial in a polydisc, which was discovered implicitly by Shearer and explicitly by Dobrushin. We also present some refinements and extensions of both arguments, including a generalization of the Lovasz local lemma that allows for "soft" dependencies. In addition, we prove some general properties of the partition function of a repulsive lattice gas, most of which are consequences of the alternating-sign property for the Mayer coefficients. We conclude with a brief discussion of the repulsive lattice gas on countably infinite graphs.
I show that the zeros of the chromatic polynomials P_G(q) for the generalized theta graphs Θ^{(s,p)} are, taken together, dense in the whole complex plane with the possible exception of the disc |q-1| < 1. The same holds for their dichromatic polynomials (alias Tutte polynomials, alias Potts-model partition functions) Z_G(q,v) outside the disc |q+v| < |v|. An immediate corollary is that the chromatic zeros of not-necessarily-planar graphs are dense in the whole complex plane. The main technical tool in the proof of these results is the Beraha-Kahane-Weiss theorem on the limit sets of zeros for certain sequences of analytic functions, for which I give a new and simpler proof.
We study the chromatic polynomial P_G(q) for m \times n triangular-lattice strips of widths m <= 12_P, 9_F (with periodic or free transverse boundary conditions, respectively) and arbitrary lengths n (with free longitudinal boundary conditions). The chromatic polynomial gives the zero-temperature limit of the partition function for the q-state Potts antiferromagnet. We compute the transfer matrix for such strips in the Fortuin--Kasteleyn representation and obtain the corresponding accumulation sets of chromatic zeros in the complex q-plane in the limit n\to\infty. We recompute the limiting curve obtained by Baxter in the thermodynamic limit m,n\to\infty and find new interesting features with possible physical consequences. Finally, we analyze the isolated limiting points and their relation with the Beraha numbers.
We study the chromatic polynomials (= zero-temperature antiferromagnetic Potts-model partition functions) P_G(q) for m \times n rectangular subsets of the square lattice, with m \le 8 (free or periodic transverse boundary conditions) and n arbitrary (free longitudinal boundary conditions), using a transfer matrix in the Fortuin-Kasteleyn representation. In particular, we extract the limiting curves of partition-function zeros when n \to\infty, which arise from the crossing in modulus of dominant eigenvalues (Beraha-Kahane-Weiss theorem). We also provide evidence that the Beraha numbers B_2,B_3,B_4,B_5 are limiting points of partition-function zeros as n \to\infty whenever the strip width m is \ge 7 (periodic transverse b.c.) or \ge 8 (free transverse b.c.). Along the way, we prove that a noninteger Beraha number (except perhaps B_{10}) cannot be a chromatic root of any graph.
At zero temperature, the 3-state antiferromagnetic Potts model on a square lattice maps exactly onto a point of the 6-vertex model whose long-distance behavior is equivalent to that of a free scalar boson. We point out that at nonzero temperature there are two distinct types of excitation: vortices, which are relevant with renormalization-group eigenvalue 1/2; and non-vortex unsatisfied bonds, which are strictly marginal and serve only to renormalize the stiffness coefficient of the underlying free boson. Together these excitations lead to an unusual form for the corrections to scaling: for example, the correlation length diverges as β\equiv J/kT \to \infty according to ξ\sim A e^{2β} (1 + bβe^{-β} + ...), where b is a nonuniversal constant that may nevertheless be determined independently. A similar result holds for the staggered susceptibility. These results are shown to be consistent with the anomalous behavior found in the Monte Carlo simulations of Ferreira and Sokal.
I show that there exist universal constants $C(r) < \infty$ such that, for all loopless graphs $G$ of maximum degree $\le r$, the zeros (real or complex) of the chromatic polynomial $P_G(q)$ lie in the disc $|q| < C(r)$. Furthermore, $C(r) \le 7.963906... r$. This result is a corollary of a more general result on the zeros of the Potts-model partition function $Z_G(q, {v_e})$ in the complex antiferromagnetic regime $|1 + v_e| \le 1$. The proof is based on a transformation of the Whitney-Tutte-Fortuin-Kasteleyn representation of $Z_G(q, {v_e})$ to a polymer gas, followed by verification of the Dobrushin-Kotecký-Preiss condition for nonvanishing of a polymer-model partition function. I also show that, for all loopless graphs $G$ of second-largest degree $\le r$, the zeros of $P_G(q)$ lie in the disc $|q| < C(r) + 1$. Along the way, I give a simple proof of a generalized (multivariate) Brown-Colbourn conjecture on the zeros of the reliability polynomial for the special case of series-parallel graphs.
The generalized theta graph Θ_{s_1,...,s_k} consists of a pair of endvertices joined by k internally disjoint paths of lengths s_1,...,s_k \ge 1. We prove that the roots of the chromatic polynomial $pi(Θ_{s_1,...,s_k},z) of a k-ary generalized theta graph all lie in the disc |z-1| \le [1 + o(1)] k/\log k, uniformly in the path lengths s_i. Moreover, we prove that Θ_{2,...,2} \simeq K_{2,k} indeed has a chromatic root of modulus [1 + o(1)] k/\log k. Finally, for k \le 8 we prove that the generalized theta graph with a chromatic root that maximizes |z-1| is the one with all path lengths equal to 2; we conjecture that this holds for all k.
We study the antiferromagnetic q-state Potts model on the square lattice for q=3 and q=4, using the Wang-Swendsen-Kotecky (WSK) Monte Carlo algorithm and a powerful finite-size-scaling extrapolation method. For q=3 we obtain good control up to correlation length $ξ\sim 5000$; the data are consistent with $ξ(β) = A e^{2β} β^p (1 + a_1 e^{-β} + ...)$ as $β\to\infty$, with $p \approx 1$. The staggered susceptibility behaves as $χ_{stagg} \sim ξ^{5/3}$. For q=4 the model is disordered ($ξ\ltapprox 2$) even at zero temperature. In appendices we prove a correlation inequality for Potts antiferromagnets on a bipartite lattice, and we prove ergodicity of the WSK algorithm at zero temperature for Potts antiferromagnets on a bipartite lattice.
Using results from conformal field theory, we compute several universal amplitude ratios for the two-dimensional Ising model at criticality on a symmetric torus. These include the correlation-length ratio x^\star = \lim_{L\to\infty} ξ(L)/L and the first four magnetization moment ratios V_{2n} = <{\cal M}^{2n}>/<{\cal M}^2>^n. As a corollary we get the first four renormalized 2n-point coupling constants for the massless theory on a symmetric torus, G_{2n}^*. We confirm these predictions by a high-precision Monte Carlo simulation.
We have performed a high-precision Monte Carlo study of the dynamic critical behavior of the Swendsen-Wang algorithm for the two-dimensional 3-state Potts model. We find that the Li-Sokal bound ($τ_{int,E} \geq const \times C_H$) is almost but not quite sharp. The ratio $τ_{int,E} / C_H$ seems to diverge either as a small power ($\approx 0.08$) or as a logarithm.
We introduce a variant of the multi-grid Monte Carlo (MGMC) method, based on the embedding of an $XY$ model into the target model, and we study its mathematical properties for a variety of nonlinear $σ$-models. We then apply the method to the two-dimensional $O(N)$-symmetric nonlinear $σ$-models (also called $N$-vector models) with $N=3,4,8$ and study its dynamic critical behavior. Using lattices up to $256 \times 256$, we find dynamic critical exponents $z_{int,{\cal M}^2} \approx 0.70 \pm 0.08$, $0.60 \pm 0.07$, $0.52 \pm 0.10$ for $N=3,4,8$, respectively (subjective 68\% confidence intervals). Thus, for these asymptotically free models, critical slowing-down is greatly reduced compared to local algorithms, but not completely eliminated; and the dynamic critical exponent does apparently vary with $N$. We also analyze the static data for $N=8$ using a finite-size-scaling extrapolation method. The correlation length $ξ$ agrees with the four-loop asymptotic-freedom prediction to within $\approx 1\%$ over the interval $12 < ξ< 650$.
We carry out a high-precision simulation of the two-dimensional $SU(3)$ principal chiral model at correlation lengths $ξ$ up to $\sim 4 \times 10^5$, using a multi-grid Monte Carlo (MGMC) algorithm and approximately one year of Cray C-90 CPU time. We extrapolate the finite-volume Monte Carlo data to infinite volume using finite-size-scaling theory, and we discuss carefully the systematic and statistical errors in this extrapolation. We then compare the extrapolated data to the renormalization-group predictions. The deviation from asymptotic scaling, which is $\approx 12%$ at $ξ\sim 25$, decreases to $\approx 2%$ at $ξ\sim 4 \times 10^5$. We also analyze the dynamic critical behavior of the MGMC algorithm using lattices up to $256 \times 256$, finding the dynamic critical exponent $z_{int,{\cal M}^2} \approx 0.45 \pm 0.02$ (subjective 68% confidence interval). Thus, for this asymptotically free model, critical slowing-down is greatly reduced compared to local algorithms, but not completely eliminated.
We solve exactly the general one-dimensional $O(N)$-invariant spin model taking values in the sphere $S^{N-1}$, with nearest-neighbor interactions, in finite volume with periodic boundary conditions, by an expansion in hyperspherical harmonics. The possible continuum limits are discussed for a general one-parameter family of interactions, and an infinite number of universality classes is found. For these classes we compute the finite-size-scaling functions and the leading corrections to finite-size scaling. A special two-parameter family of interactions (which includes the mixed isovector/isotensor model) is also treated, and no additional universality classes appear. In the appendices we give new formulae for the Clebsch-Gordan coefficients and 6--$j$ symbols of the $O(N)$ group, and some new generalizations of the Poisson summation formula; these may be of independent interest.
We study the antiferromagnetic $q$-state Potts model on the square lattice for $q=3$ and $q=4$, using the Wang-Swendsen-Kotecký Monte Carlo algorithm and a new finite-size-scaling extrapolation method. For $q=3$ we obtain good control up to correlation length $ξ\sim 5000$; the data are consistent with $ξ(β) = A e^{2β} β(1 + a_1 e^{-β} + \ldots)$ as $β\to\infty$. For $q=4$ the model is disordered even at zero temperature.
We prove that for any Monte Carlo algorithm of Metropolis type, the autocorrelation time of a suitable ``energy''-like observable is bounded below by a multiple of the corresponding ``specific heat''. This bound does not depend on whether the proposed moves are local or non-local; it depends only on the distance between the desired probability distribution $π$ and the probability distribution $π^{(0)}$ for which the proposal matrix satisfies detailed balance. We show, with several examples, that this result is particularly powerful when applied to non-local algorithms.
I make three conceptual points regarding multigrid methods for computing propagators in lattice gauge theory: 1) The class of operators handled by the algorithm must be stable under coarsening. 2) Problems related by symmetry should have solution methods related by symmetry. 3) It is crucial to distinguish the vector space $V$ from its dual space $V^*$. All the existing algorithms violate one or more of these principles.
Previous theories of dilute polymer solutions have failed to distinguish clearly between two very different ways of taking the long-chain limit: (I) $N \to\infty$ at fixed temperature $T$, and (II) $N \to\infty$, $T \to T_θ$ with $x \equiv N^ϕ(T-T_θ)$ fixed. I argue that the modern two-parameter theory (continuum Edwards model) applies to case II --- not case I --- and in fact gives exactly the crossover scaling functions for $x \ge 0$ modulo two nonuniversal scale factors. A Wilson-type renormalization group clarifies the connection between crossover scaling functions and continuum field theories. [Also contains a general discussion of the connection between the Wilson and field-theoretic renormalization groups. Comments solicited.]