Source author record

I. Jensen

I. Jensen appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math-ph math.MP cond-mat.stat-mech math.CO cond-mat hep-lat hep-th nlin.SI physics.comp-ph

Catalog footprint

What is connected

9works

9topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2016arXiv

Is the full susceptibility of the square-lattice Ising model a differentially algebraic function?

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$.

preprint2014arXiv

Integrability vs non-integrability: Hard hexagons and hard squares compared

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.

preprint2013arXiv

Hard hexagon partition function for complex fugacity

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.

preprint2010arXiv

Square lattice Ising model $\tildeχ^{(5)}$ ODE in exact arithmetic

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}(Φ^{(5)}) = E^{(5)}$, where $Φ^{(5)} =\tildeχ^{(5)}-\tildeχ^{(3)}/2+\tildeχ^{(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χ^{(5)}$ at the ferromagnetic critical point and to obtain a limited number of analytic continuations of $\tildeχ^{(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χ^{(5)}$ is singular at $w=1/2$ on an infinite number of branches.

preprint2009arXiv

High order Fuchsian equations for the square lattice Ising model: $\tildeχ^{(5)}$

We consider the Fuchsian linear differential equation obtained (modulo a prime) for $\tildeχ^{(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χ^{(1)}$ and $\tildeχ^{(3)}$ can be removed from $\tildeχ^{(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χ^{(3)}$ and $\tildeχ^{(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χ^{(n)}$'s.

preprint2009arXiv

High order Fuchsian equations for the square lattice Ising model: $χ^{(6)}$

This paper deals with $\tildeχ^{(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χ^{(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 $Φ^{(6)}=\tildeχ^{(6)} - {2 \over 3} \tildeχ^{(4)} + {2 \over 45} \tildeχ^{(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χ^{(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.

preprint2008arXiv

Exact generating function for 2-convex polygons

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.

preprint2005arXiv

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 $λ^{L^2+o(L^2)}$. We estimate $λ= 1.744550 \pm 0.000005$ as well as obtaining strict upper and lower bounds, $1.628 < λ< 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/μ$ the average length of a SAW grows as $L$, while for $x > 1/μ$ it grows as $L^2$. Here $μ$ is the growth constant of unconstrained SAW in ${\mathbb Z}^2$. For $x = 1/μ$ 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 $τ^{L^2+o(L^2)}$ on the same $L \times L$ lattice. We give precise estimates for $τ$ as well as upper and lower bounds, and prove that $τ< λ.$

preprint1994arXiv

Low-Temperature Series Expansions for the Spin-1 Ising Model

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 $Δ_{1} = 1.1 \pm 0.1$ is found.