Search

We discuss the implications of studies of partition function zeros and equimodular curves for the analytic properties of the Ising model on a square lattice in a magnetic field. In particular we consider the dense set of singularities in the susceptibility of the Ising model at $H=0$ found by Nickel and its relation to the analyticity of the field theory computations of Fonseca and Zamolodchikov., Comment: 21 pages, 13 figures

We consider self-avoiding lattice polygons, in the hypercubic lattice, as a model of a ring polymer adsorbed at a surface and either being desorbed by the action of a force, or pushed towards the surface. We show that, when there is no interaction with the surface, then the response of the polygon to the applied force is identical (in the thermodynamic limit) for two ways in which we apply the force. When the polygon is attracted to the surface then, when the dimension is at least 3, we have a complete characterization of the critical force--temperature curve in terms of the behaviour, (a) when there is no force, and, (b) when there is no surface interaction. For the 2-dimensional case we have upper and lower bounds on the free energy. We use both Monte Carlo and exact enumeration and series analysis methods to investigate the form of the phase diagram in two dimensions. We find evidence for the existence of a \emph{mixed phase} where the free energy depends on the strength of the interaction with the adsorbing line and on the applied force.

We study the class of non-holonomic power series with integer coefficients that reduce, modulo primes, or powers of primes, to algebraic functions. In particular we try to determine whether the susceptibility of the square-lattice Ising model belongs to this class, and more broadly whether the susceptibility is a solution of a differentially algebraic equation. Initial results on Tutte's non-linear ordinary differential equation (ODE) and other simple quadratic non-linear ODEs suggest that a large set of differentially algebraic power series solutions with integer coefficients might reduce to algebraic functions modulo primes, or powers of primes. Here we give several examples of series with integer coefficients and non-zero radius of convergence that reduce to algebraic functions modulo (almost) every prime (or power of a prime). These examples satisfy differentially algebraic equations with the encoding polynomial occasionally possessing quite high degree (and thus difficult to identify even with long series). Additionally, we have extended both the high- and low-temperature Ising square-lattice susceptibility series to 5043 coefficients. We find that even this long series is insufficient to determine whether it reduces to algebraic functions modulo $3$, $5$, etc. This negative result is in contrast to the comparatively easy confirmation that the corresponding series reduce to algebraic functions modulo powers of $2$., Comment: 36 pages

In this paper we compare the integrable hard hexagon model with the non-integrable hard squares model by means of partition function roots and transfer matrix eigenvalues. We consider partition functions for toroidal, cylindrical, and free-free boundary conditions up to sizes $40\times40$ and transfer matrices up to 30 sites. For all boundary conditions the hard squares roots are seen to lie in a bounded area of the complex fugacity plane along with the universal hard core line segment on the negative real fugacity axis. The density of roots on this line segment matches the derivative of the phase difference between the eigenvalues of largest (and equal) moduli and exhibits much greater structure than the corresponding density of hard hexagons. We also study the special point $z=-1$ of hard squares where all eigenvalues have unit modulus, and we give several conjectures for the value at $z=-1$ of the partition functions., Comment: 46 pages

We investigate the behaviour of the mean size of directed compact percolation clusters near a damp wall in the low-density region, where sites in the bulk are wet (occupied) with probability $p$ while sites on the wall are wet with probability $p_w$. Methods used to find the exact solution for the dry case ($p_w=0$) and the wet case ($p_w=1$) turn out to be inadequate for the damp case. Instead we use a series expansion for the $p_w=2p$ case to obtain a second order inhomogeneous differential equation satisfied by the mean size, which exhibits a critical exponent $\gamma=2$, in common with the wet wall result. For the more general case of $p_w=rp$, with $r$ rational, we use a modular arithmetic method of finding ODEs and obtain a fourth order homogeneous ODE satisfied by the series. The ODE is expressed exactly in terms of $r$. We find that in the damp region $0

We study the analyticity of the partition function of the hard hexagon model in the complex fugacity plane by computing zeros and transfer matrix eigenvalues for large finite size systems. We find that the partition function per site computed by Baxter in the thermodynamic limit for positive real values of the fugacity is not sufficient to describe the analyticity in the full complex fugacity plane. We also obtain a new algebraic equation for the low density partition function per site., Comment: 49 pages, IoP styles files, lots of figures (png mostly) so using PDFLaTeX. Some minor changes added to version 2 in response to referee reports

We obtain in exact arithmetic the order 24 linear differential operator $L_{24}$ and right hand side $E^{(5)}$ of the inhomogeneous equation$L_{24}(\Phi^{(5)}) = E^{(5)}$, where $\Phi^{(5)} =\tilde{\chi}^{(5)}-\tilde{\chi}^{(3)}/2+\tilde{\chi}^{(1)}/120$ is a linear combination of $n$-particle contributions to the susceptibility of the square lattice Ising model. In Bostan, et al. (J. Phys. A: Math. Theor. {\bf 42}, 275209 (2009)) the operator $L_{24}$ (modulo a prime) was shown to factorize into $L_{12}^{(\rm left)} \cdot L_{12}^{(\rm right)}$; here we prove that no further factorization of the order 12 operator $L_{12}^{(\rm left)}$ is possible. We use the exact ODE to obtain the behaviour of $\tilde{\chi}^{(5)}$ at the ferromagnetic critical point and to obtain a limited number of analytic continuations of $\tilde{\chi}^{(5)}$ beyond the principal disk defined by its high temperature series. Contrary to a speculation in Boukraa, et al (J. Phys. A: Math. Theor. {\bf 41} 455202 (2008)), we find that $\tilde{\chi}^{(5)}$ is singular at $w=1/2$ on an infinite number of branches., Comment: 25 pages, 2 figures, IoP style files

This paper deals with $\tilde{\chi}^{(6)}$, the six-particle contribution to the magnetic susceptibility of the square lattice Ising model. We have generated, modulo a prime, series coefficients for $\tilde{\chi}^{(6)}$. The length of the series is sufficient to produce the corresponding Fuchsian linear differential equation (modulo a prime). We obtain the Fuchsian linear differential equation that annihilates the "depleted" series $\Phi^{(6)}=\tilde{\chi}^{(6)} - {2 \over 3} \tilde{\chi}^{(4)} + {2 \over 45} \tilde{\chi}^{(2)}$. The factorization of the corresponding differential operator is performed using a method of factorization modulo a prime introduced in a previous paper. The "depleted" differential operator is shown to have a structure similar to the corresponding operator for $\tilde{\chi}^{(5)}$. It splits into factors of smaller orders, with the left-most factor of order six being equivalent to the symmetric fifth power of the linear differential operator corresponding to the elliptic integral $E$. The right-most factor has a direct sum structure, and using series calculated modulo several primes, all the factors in the direct sum have been reconstructed in exact arithmetics., Comment: 23 pages

We consider the Fuchsian linear differential equation obtained (modulo a prime) for $\tilde{\chi}^{(5)}$, the five-particle contribution to the susceptibility of the square lattice Ising model. We show that one can understand the factorization of the corresponding linear differential operator from calculations using just a single prime. A particular linear combination of $\tilde{\chi}^{(1)}$ and $\tilde{\chi}^{(3)}$ can be removed from $\tilde{\chi}^{(5)}$ and the resulting series is annihilated by a high order globally nilpotent linear ODE. The corresponding (minimal order) linear differential operator, of order 29, splits into factors of small orders. A fifth order linear differential operator occurs as the left-most factor of the "depleted" differential operator and it is shown to be equivalent to the symmetric fourth power of $L_E$, the linear differential operator corresponding to the elliptic integral $E$. This result generalizes what we have found for the lower order terms $\tilde{\chi}^{(3)}$ and $\tilde{\chi}^{(4)}$. We conjecture that a linear differential operator equivalent to a symmetric $(n-1)$-th power of $L_E$ occurs as a left-most factor in the minimal order linear differential operators for all $\tilde{\chi}^{(n)}$'s., Comment: 33 pages

We investigate a two-dimensional problem of an isolated self-interacting end-grafted polymer, pulled by one end. In the thermodynamic limit, we find that the model has only two different phases, namely a collapsed phase and a stretched phase. We show that the phase diagram obtained by Kumar {\it at al.\} [Phys. Rev. Lett. {\bf 98}, 128101 (2007)] for small systems, where differences between various statistical ensembles play an important role, differ from the phase diagram obtained here in the thermodynamic limit., Comment: 20 pages, 22 figures

We calculate very long low- and high-temperature series for the susceptibility $\chi$ of the square lattice Ising model as well as very long series for the five-particle contribution $\chi^{(5)}$ and six-particle contribution $\chi^{(6)}$. These calculations have been made possible by the use of highly optimized polynomial time modular algorithms and a total of more than 150000 CPU hours on computer clusters. For $\chi^{(5)}$ 10000 terms of the series are calculated {\it modulo} a single prime, and have been used to find the linear ODE satisfied by $\chi^{(5)}$ {\it modulo} a prime. A diff-Pad\'e analysis of 2000 terms series for $\chi^{(5)}$ and $\chi^{(6)}$ confirms to a very high degree of confidence previous conjectures about the location and strength of the singularities of the $n$-particle components of the susceptibility, up to a small set of ``additional'' singularities. We find the presence of singularities at $w=1/2$ for the linear ODE of $\chi^{(5)}$, and $w^2= 1/8$ for the ODE of $\chi^{(6)}$, which are {\it not} singularities of the ``physical'' $\chi^{(5)}$ and $\chi^{(6)},$ that is to say the series-solutions of the ODE's which are analytic at $w =0$. Furthermore, analysis of the long series for $\chi^{(5)}$ (and $\chi^{(6)}$) combined with the corresponding long series for the full susceptibility $\chi$ yields previously conjectured singularities in some $\chi^{(n)}$, $n \ge 7$. We also present a mechanism of resummation of the logarithmic singularities of the $\chi^{(n)}$ leading to the known power-law critical behaviour occurring in the full $\chi$, and perform a power spectrum analysis giving strong arguments in favor of the existence of a natural boundary for the full susceptibility $\chi$., Comment: 54 pages, 2 figures

Single molecule force spectroscopy provide details of the underlying energy surfaces of proteins which are essential to the understanding of their unfolding process. Recently, it has been observed experimentally that by pulling proteins in different directions relative to their secondary structure, one can gain a better understanding of the shape of the energy landscape. We consider simple lattice models which are anisotropic in nature to study the response of a force in unfolding of a polymer. Our analytical solution of the model, supported by extensive numerical calculations, reveal that the force temperature diagrams are very different depending on the direction of the applied force. We find that either unzipping or shearing kind transitions dominate the dynamics of the unfolding process depending solely on the direction of the applied force., Comment: 8 pages, 11 figures

Polygons are described as almost-convex if their perimeter differs from the perimeter of their minimum bounding rectangle by twice their `concavity index', $m$. Such polygons are called \emph{$m$-convex} polygons and are characterised by having up to $m$ indentations in their perimeter. We first describe how we conjectured the (isotropic) generating function for the case $m=2$ using a numerical procedure based on series expansions. We then proceed to prove this result for the more general case of the full anisotropic generating function, in which steps in the $x$ and $y$ direction are distinguished. In so doing, we develop tools that would allow for the case $m > 2$ to be studied. %In our proof we use a `divide and conquer' approach, factorising 2-convex %polygons by extending a line along the base of its indents. We then use %the inclusion-exclusion principle, the Hadamard product and extensions to %known methods to derive the generating functions for each case., Comment: 28 pages, 15 figures, IoP style

Polygons are described as almost-convex if their perimeter differs from the perimeter of their minimum bounding rectangle by twice their `concavity index', $m$. Such polygons are called \emph{$m$-convex} polygons and are characterised by having up to $m$ indentations in the side. We use a `divide and conquer' approach, factorising 2-convex polygons by extending a line along the base of its indents. We then use the inclusion-exclusion principle, the Hadamard product and extensions to known methods to derive the generating functions for each case., Comment: 53 pages

We study a restricted class of self-avoiding walks (SAW) which start at the origin (0, 0), end at $(L, L)$, and are entirely contained in the square $[0, L] \times [0, L]$ on the square lattice ${\mathbb Z}^2$. The number of distinct walks is known to grow as $\lambda^{L^2+o(L^2)}$. We estimate $\lambda = 1.744550 \pm 0.000005$ as well as obtaining strict upper and lower bounds, $1.628 < \lambda < 1.782.$ We give exact results for the number of SAW of length $2L + 2K$ for $K = 0, 1, 2$ and asymptotic results for $K = o(L^{1/3})$. We also consider the model in which a weight or {\em fugacity} $x$ is associated with each step of the walk. This gives rise to a canonical model of a phase transition. For $x < 1/\mu$ the average length of a SAW grows as $L$, while for $x > 1/\mu$ it grows as $L^2$. Here $\mu$ is the growth constant of unconstrained SAW in ${\mathbb Z}^2$. For $x = 1/\mu$ we provide numerical evidence, but no proof, that the average walk length grows as $L^{4/3}$. We also consider Hamiltonian walks under the same restriction. They are known to grow as $\tau^{L^2+o(L^2)}$ on the same $L \times L$ lattice. We give precise estimates for $\tau$ as well as upper and lower bounds, and prove that $\tau < \lambda.$, Comment: 27 pages, 9 figures. Paper updated and reorganised following refereeing

We analyse new exact enumeration data for self-avoiding polygons, counted by perimeter and area on the square, triangular and hexagonal lattices. In extending earlier analyses, we focus on the perimeter moments in the vicinity of the bicritical point. We also consider the shape of the critical curve near the bicritical point, which describes the crossover to the branched polymer phase. Our recently conjectured expression for the scaling function of rooted self-avoiding polygons is further supported. For (unrooted) self-avoiding polygons, the analysis reveals the presence of an additional additive term with a new universal amplitude. We conjecture the exact value of this amplitude., Comment: 17 pages, 3 figures

Using a technique based on production of ion Coulomb crystals, the isotope shifts of the 4s2 1S0 -> 4s5p 1P1 transition for all naturally occurring isotopes of calcium as well as the hyperfine splitting of the 4s5p 1P1 state in Ca-43 have been measured. The field shift and specific mass shift coefficients as well as the hyperfine structure constants for Ca-43 have been derived from the data., Comment: 6 pages,6 figures, to be published in Phys. Rev. A

Exactly solvable models of planar polygons, weighted by perimeter and area, have deepened our understanding of the critical behaviour of polygon models in recent years. Based on these results, we derive a conjecture for the exact form of the critical scaling function for planar self-avoiding polygons. The validity of this conjecture was recently tested numerically using exact enumeration data for small values of the perimeter on the square and triangular lattices. We have substantially extended these enumerations and also enumerated polygons on the hexagonal lattice. We also performed Monte-Carlo simulations of the model on the square lattice. Our analysis supports the conjecture that the scaling function is given by the logarithm of an Airy function., Comment: 11 pages, 5 figures

Percolation on two-dimensional small-world networks has been proposed as a model for the spread of plant diseases. In this paper we give an analytic solution of this model using a combination of generating function methods and high-order series expansion. Our solution gives accurate predictions for quantities such as the position of the percolation threshold and the typical size of disease outbreaks as a function of the density of "shortcuts" in the small-world network. Our results agree with scaling hypotheses and numerical simulations for the same model., Comment: 7 pages, 3 figures, 2 tables

We analyze new data for self-avoiding polygons, on the square and triangular lattices, enumerated by both perimeter and area, providing evidence that the scaling function is the logarithm of an Airy function. The results imply universal amplitude combinations for all area moments and suggest that rooted self-avoiding polygons may satisfy a $q$-algebraic functional equation., Comment: 9 pages

We derive low-temperature series (in the variable $u = \exp[-\beta J/S^2]$) for the spontaneous magnetisation, susceptibility and specific heat of the spin-$S$ Ising model on the square lattice for $S=\frac32$, 2, $\frac52$, and 3. We determine the location of the physical critical point and non-physical singularities. The number of non-physical singularities closer to the origin than the physical critical point grows quite rapidly with $S$. The critical exponents at the singularities which are closest to the origin and for which we have reasonably accurate estimates are independent of $S$. Due to the many non-physical singularities, the estimates for the physical critical point and exponents are poor for higher values of $S$, though consistent with universality., Comment: 14 pages, LaTeX with IOP style files (ioplppt.sty), epic.sty and eepic.sty. To appear in J. Phys. A

The finite-lattice method of series expansions has been combined with a new extrapolation procedure to extend the low-temperature series for the specific heat, spontaneous magnetisation and susceptibility of the spin-1 Ising model on the square lattice. The extended series were derived by directly calculating the series to order 99 (in the variable $u=\exp[-J/k_B T])$ and using the new extrapolation procedure to calculate an additional 13--14 terms., Comment: 12 pages, LaTeX with IOP style files (ioplppt.sty) and epsf.sty. 3 ps files included. To appear in J. Phys. A

Series expansion methods are used to study directed bond percolation clusters on the square lattice whose lateral growth is restricted by a wall parallel to the growth direction. The percolation threshold $p_c$ is found to be the same as that for the bulk. However the values of the critical exponents for the percolation probability and mean cluster size are quite different from those for the bulk and are estimated by $\beta_1 = 0.7338 \pm 0.0001$ and $\gamma_1 = 1.8207 \pm 0.0004$ respectively. On the other hand the exponent $\Delta_1=\beta_1 +\gamma_1$ characterising the scale of the cluster size distribution is found to be unchanged by the presence of the wall. The parallel connectedness length, which is the scale for the cluster length distribution, has an exponent which we estimate to be $\nu_{1\parallel} = 1.7337 \pm 0.0004$ and is also unchanged. The exponent $\tau_1$ of the mean cluster length is related to $\beta_1$ and $\nu_{1\parallel}$ by the scaling relation $\nu_{1\parallel} = \beta_1 + \tau_1$ and using the above estimates yields $\tau_1 = 1$ to within the accuracy of our results. We conjecture that this value of $\tau_1$ is exact and further support for the conjecture is provided by the direct series expansion estimate $\tau_1= 1.0002 \pm 0.0003$., Comment: 12pages LaTeX, ioplppt.sty, to appear in J. Phys. A

The finite lattice method of series expansion has been used to extend low-temperature series for the partition function, order parameter and susceptibility of the spin-1 Ising model on the square lattice. A new formalism is described that uses two distinct transfer matrix approaches in order to significantly reduce computer memory requirements and which permits the derivation of the series to 79th order. Subsequent analysis of the series clearly confirms that the spin-1 model has the same dominant critical exponents as the spin-$\frac{1}{2}$ Ising model. Accurate estimates for both the critical temperature and non-physical singularities are obtained. In addition, evidence for a non-analytic confluent correction with exponent $\Delta_{1} = 1.1 \pm 0.1$ is found., Comment: 24 pages (+7 uuencoded postscript figures), LaTeX