A new transfer-matrix algorithm for exact enumerations: Self-avoiding polygons on the square lattice

We present a new and more efficient implementation of transfer-matrix methods for exact enumerations of lattice objects. The new method is illustrated by an application to the enumeration of self-avoiding polygons on the square lattice. A detailed comparison with the previous best algorithm shows significant improvement in the running time of the algorithm. The new algorithm is used to extend the enumeration of polygons to length 130 from the previous record of 110.Comment: 17 pages, 8 figures, IoP style file

Series expansions of the percolation probability on the directed triangular lattice

We have derived long series expansions of the percolation probability for site, bond and site-bond percolation on the directed triangular lattice. For the bond problem we have extended the series from order 12 to 51 and for the site problem from order 12 to 35. For the site-bond problem, which has not been studied before, we have derived the series to order 32. Our estimates of the critical exponent $\beta$ are in full agreement with results for similar problems on the square lattice, confirming expectations of universality. For the critical probability and exponent we find in the site case: $q_c = 0.4043528 \pm 0.0000010$ and $\beta = 0.27645 \pm 0.00010$; in the bond case: $q_c = 0.52198\pm 0.00001$ and $\beta = 0.2769\pm 0.0010$; and in the site-bond case: $q_c = 0.264173 \pm 0.000003$ and $\beta = 0.2766 \pm 0.0003$. In addition we have obtained accurate estimates for the critical amplitudes. In all cases we find that the leading correction to scaling term is analytic, i.e., the confluent exponent $\Delta = 1$.Comment: 26 pages, LaTeX. To appear in J. Phys.

Directed percolation near a wall

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.

Complex-Temperature Singularities in the d=2d=2 Ising Model. III. Honeycomb Lattice

We study complex-temperature properties of the uniform and staggered susceptibilities $\chi$ and $\chi^{(a)}$ of the Ising model on the honeycomb lattice. From an analysis of low-temperature series expansions, we find evidence that $\chi$ and $\chi^{(a)}$ both have divergent singularities at the point $z=-1 \equiv z_{\ell}$ (where $z=e^{-2K}$), with exponents $\gamma_{\ell}'= \gamma_{\ell,a}'=5/2$. The critical amplitudes at this singularity are calculated. Using exact results, we extract the behaviour of the magnetisation $M$ and specific heat $C$ at complex-temperature singularities. We find that, in addition to its zero at the physical critical point, $M$ diverges at $z=-1$ with exponent $\beta_{\ell}=-1/4$, vanishes continuously at $z=\pm i$ with exponent $\beta_s=3/8$, and vanishes discontinuously elsewhere along the boundary of the complex-temperature ferromagnetic phase. $C$ diverges at $z=-1$ with exponent $\alpha_{\ell}'=2$ and at $v=\pm i/\sqrt{3}$ (where $v = \tanh K$) with exponent $\alpha_e=1$, and diverges logarithmically at $z=\pm i$. We find that the exponent relation $\alpha'+2\beta+\gamma'=2$ is violated at $z=-1$; the right-hand side is 4 rather than 2. The connections of these results with complex-temperature properties of the Ising model on the triangular lattice are discussed.Comment: 22 pages, latex, figures appended after the end of the text as a compressed, uuencoded postscript fil

Partially directed paths in a wedge

The enumeration of lattice paths in wedges poses unique mathematical challenges. These models are not translationally invariant, and the absence of this symmetry complicates both the derivation of a functional recurrence for the generating function, and solving for it. In this paper we consider a model of partially directed walks from the origin in the square lattice confined to both a symmetric wedge defined by $Y = \pm pX$, and an asymmetric wedge defined by the lines $Y= pX$ and Y=0, where $p > 0$ is an integer. We prove that the growth constant for all these models is equal to $1+\sqrt{2}$, independent of the angle of the wedge. We derive functional recursions for both models, and obtain explicit expressions for the generating functions when $p=1$. From these we find asymptotic formulas for the number of partially directed paths of length $n$ in a wedge when $p=1$. The functional recurrences are solved by a variation of the kernel method, which we call the ``iterated kernel method''. This method appears to be similar to the obstinate kernel method used by Bousquet-Melou. This method requires us to consider iterated compositions of the roots of the kernel. These compositions turn out to be surprisingly tractable, and we are able to find simple explicit expressions for them. However, in spite of this, the generating functions turn out to be similar in form to Jacobi $\theta$-functions, and have natural boundaries on the unit circle.Comment: 26 pages, 5 figures. Submitted to JCT

Force induced triple point for interacting polymers

We show the existence of a force induced triple point in an interacting polymer problem that allows two zero-force thermal phase transitions. The phase diagrams for three different models of mutually attracting but self avoiding polymers are presented. One of these models has an intermediate phase and it shows a triple point but not the others. A general phase diagram with multicritical points in an extended parameter space is also discussed.Comment: 4 pages, 8 figures, revtex

Unexpected Spin-Off from Quantum Gravity

We propose a novel way of investigating the universal properties of spin systems by coupling them to an ensemble of causal dynamically triangulated lattices, instead of studying them on a fixed regular or random lattice. Somewhat surprisingly, graph-counting methods to extract high- or low-temperature series expansions can be adapted to this case. For the two-dimensional Ising model, we present evidence that this ameliorates the singularity structure of thermodynamic functions in the complex plane, and improves the convergence of the power series.Comment: 10 pages, 4 figures; final, slightly amended version, to appear in Physica

Self-avoiding walks crossing a square

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 refereein

Test of Guttmann and Enting's conjecture in the eight-vertex model

We investigate the analyticity property of the partially resummed series expansion(PRSE) of the partition function for the eight-vertex model. Developing a graphical technique, we have obtained a first few terms of the PRSE and found that these terms have a pole only at one point in the complex plane of the coupling constant. This result supports the conjecture proposed by Guttmann and Enting concerning the ``solvability'' in statistical mechanical lattice models.Comment: 15 pages, 3 figures, RevTe

Effects of Eye-phase in DNA unzipping

The onset of an "eye-phase" and its role during the DNA unzipping is studied when a force is applied to the interior of the chain. The directionality of the hydrogen bond introduced here shows oscillations in force-extension curve similar to a "saw-tooth" kind of oscillations seen in the protein unfolding experiments. The effects of intermediates (hairpins) and stacking energies on the melting profile have also been discussed.Comment: RevTeX v4, 9 pages with 7 eps figure