# Scaling behaviour of lattice animals at the upper critical dimension

###### Abstract

We perform numerical simulations of the lattice-animal problem at the upper critical dimension on hypercubic lattices in order to investigate logarithmic corrections to scaling there. Our stochastic sampling method is based on the pruned-enriched Rosenbluth method (PERM), appropriate to linear polymers, and yields high statistics with animals comprised of up to 8000 sites. We estimate both the partition sums (number of different animals) and the radii of gyration. We re-verify the Parisi-Sourlas prediction for the leading exponents and compare the logarithmic-correction exponents to two partially differing sets of predictions from the literature. Finally, we propose, and test, a new Parisi-Sourlas-type scaling relation appropriate for the logarithmic-correction exponents.

## 1 Introduction

A lattice animal is a cluster of connected sites on a regular lattice. The enumeration of such objects – also called polyominoes – is a combinatorial problem of interest to mathematicians [1], while in physics, they are closely linked to the problems of percolation [2] and clustering in spin models [3]. In chemistry they form a basis for models of randomly branched polymers in good solvents [4]. Lattice animals linked by translations are considered as belonging to the same equivalence class, and as such are considered to be essentially the same. Of interest is , the number of distinct animals containing sites. A related objective is the calculation of the radius of gyration , related to the average distance of occupied sites from the centre of mass of the lattice animal. A number of variants of the lattice animal are studied: bond lattice animals, which are clusters of connected bonds; weakly embedded and strongly embedded trees. It is believed that all these different models belong to the same universality class [5].

It is now established that the number of lattice animals and the radius of gyration behave, to leading order in , as [4]

(1) | |||||

(2) |

Here is related to the rate of growth of and where is the Hausdorff (fractal) dimension of the lattice animals.

It is useful to define the generating function (or grand-canonical partition function) for the model as

(3) |

The growth constant is related to the critical fugacity defining the radius of convergence of the generating function (3). To leading order,

(4) |

giving the identification . The growth constant is related to the critical fugacity of the corresponding field theory and depends on the lattice coordination number, while the entropic exponent and the correlation-length exponent are universal [6].

Lattice animals may be viewed as the graphs arising from high-temperature expansions of related magnetic models, in particular the high temperature expansion of (the derivative of) the free energy of the Potts models, which in turn is related to the percolation problem. Parisi and Sourlas related the problem of branched polymers, and hence lattice animals, in dimension with the Yang-Lee edge problem in dimension, and predicted that and are related by [7]

(5) |

This relation was re-derived in an interesting fashion [8]: identifying as the high-temperature expansion of a magnetic model leads to the identification of with the magnetic and thermal scaling fields of the related magnetic model, indicating that the lattice animal is controlled by a single scaling field. This leads to the relations . The exponents and are the usual critical indices related to the divergence of the specific heat and susceptibility of the related Potts model. Substituting the mean-field exponents into the usual hyperscaling relation would lead to an upper critical dimension of , whereas the correct upper critical dimension is , indicating that there is an anomalous scaling and that hyperscaling is modified, with replaced by :

(6) |

and Eq.(5) is recovered.

In Ref. [5], the Parisi-Sourlas predictions for the leading behaviour for both Eqs.(1) and (2) were verified in dimensions to using a high-statistics numerical study with lattice animals with up to several thousand sites in each case. The measured values of and were compatible with the Parisi-Sourlas scaling relation.

Although experimentally inaccessible, a complete understanding of the lattice-animal problem includes the upper critical dimension . At and above this dimension, the critical exponents take on the mean-field values and [9]. In eight dimensions, the scaling forms (1) and (2) are modified by multiplicative logarithmic corrections. Indeed, in the high-precision study of Ref. [5], very large corrections to Eqs.(1) and (2) were reported in eight dimensions. While it was presumed that these corrections are logarithmic in nature, no attempt at a detailed fit to them was made because the authors were unaware of theoretical predictions beyond the leading order, and because of the notorious difficulty in fitting to such logarithms. It is expected that at the upper critical dimension , and scale as

(7) | |||||

(8) |

with and . The values of the logarithmic correction exponents and are the subject of the present article.

The mean-field exponents for the lattice animal model correspond to the exponents calculated from a theory with reduced temperature and where the reduced magnetic field is used as a temperature-like variable. This is consistent with the realisation, stated above, that in this model there is only one scaling field, linked to the magnetic field of the underlying magnetic model. The full set of mean-field exponents are

(9) |

These exponents are related by Fisher renormalisation to the standard mean-field exponents obtained setting and varying . These Fisher-renormalised exponents are

(10) |

Whilst in the lattice animal model there is only a single scaling field (), in the equivalent Yang-Lee model it is possible to vary both and independently.

This pairing via Fisher renormalisation and scaling relations have permitted new analytic predictions for the logarithmic corrections [10]: , consistent with the prediction of Ref. [11] and . The latter differs from a previous renormalisation-group based prediction: [11]. We therefore considered it worthwhile to revisit the problem of lattice animals in dimensions in an attempt to discern whether the numerics support either of these analytic predictions for logarithmic corrections.

## 2 Scaling at the upper critical dimension

While the leading exponents and in Eqs.(7) and (8) are not in doubt, there are two sets of predictions in the literature for their logarithmic-correction counterparts and . In Ref. [11], Ruiz-Lorenzo analytically studied these and other logarithmic corrections for a generic scalar field theory at its upper critical dimension . This theory, with imaginary coupling, is known to describe the Yang-Lee problem [12]. The latter originates from the study of the Yang-Lee edge singularity, which may be regarded as a critical or pseudo-critical point.

Parisi and Sourlas advanced a relationship between the Yang-Lee singularity in dimensions and the lattice-animal problem in dimensions [7]. Recently, an exact mapping between the two problems established this relationship on a rigorous footing [13]. The renormalisation-group calculation of Ref. [11] for the Yang-Lee problem ( theory with imaginary coupling) in dimensions yields a free energy as a function of the magnetic field, the singular part of which is of the form

(11) |

The grand canonical partition function for the lattice animals is

(12) |

where, with given by Eq.(8), scales as

(13) |

The Parisi-Sourlas mapping, then, identifies Eq.(11) for the Yang-Lee problem with Eq.(13) for the lattice animals, with the magnetic field in the former case being replaced by the fugacity in the latter. This leads to the predictions and for the lattice-animal problem in dimensions [11], which are supported by previous direct calculations [4]. The -approach also leads to analytic predictions for the correlation length,

(14) |

which translates to and for lattice animals.

The lattice animal and Yang-Lee problems can be considered either with the field or the order parameter held constant, with constant field being the more natural in field theory [14]. The corresponding two sets of critical exponents are linked via Fisher renormalization [15]. The Fisher renormalization scheme for logarithmic-corrections was recently established in Ref. [10]. In Ref. [11], Ruiz-Lorenzo has also determined the constant-order-parameter critical exponents and, in particular, the constant-order-parameter logarithmic-correction exponent for the correlation-length is given as [11]. Applying the Fisher-renormalization relations from Ref. [10] with this value yields the prediction , which is different to the constant-field estimate quoted above.

To summarize, the renormalization-group approach of Ref. [11] yields constant-field estimates for the leading entopic and Flory critical exponents which agree with mean-field, have been checked numerically in Ref. [5] and which are not in doubt. The corresponding logarithmic exponents are and . In Ref. [11], constant-order-parameter estimates are also given, which, when Fisher renormalized also yield but . while all leading exponents and Ruiz-Lorenzo’s RG calculations for the theory agree with mean-field predictions, and while calculations for the logarithmic-correction exponents agree with all previous estimates where they exist (besides the Yang-Lee problem, these include for spin glasses and for percolation in six dimensions), the disparity between the estimates for requires further investigation and we chose a non-perturbative, numerical approach. It is also necessary to check if is supported numerically, as this has not been tested non-perturbatively before.

## 3 Numerical approach

The numerical data was obtained using the prune-enriched Rosenbluth method (PERM) which is a variant of the Rosenbluth-Rosenbluth Monte Carlo method for self-avoiding walks designed to avoid the ensembles being dominated by a few high weight clusters and to avoid undue time being used calculating clusters with small weights [16]. This is avoided by reducing the width of the weights distribution by pruning low-weight configurations while cloning high-weight ones. To apply this approach to lattice animals, we have to estimate the cluster weight while it is still growing. The approach, which is discussed in detail in Ref. [5], generates independent clusters from different Monte Carlo tours and therefore leads to straightforward estimates for the errors in the raw data for the partition sums and gyration radii.

In Fig. 1 is plotted against . The constant is set to the best fit value . The asymptotic value of the slope is consistent with the expected value of . To investigate the logarithmic corrections, the leading scaling behaviour is subtracted out, and is plotted against in Fig. 2. The middle line shows the plot corresponding to the best fit value of , whilst the other two lines correspond to the upper and lower error bounds (upper line) and (lower line). The insert shows the last section of the curve, which can be seen to have a slope consistent with . (Note that although the horizontal axis has relatively short range, since it is on a log-log scale it corresponds to a wide range of animals sizes, from to .)

The slope in the sizes of animals calculated is sensitive to the precise value of . One may attempt to eliminate by using two values of , as

(15) |

In Fig. 3 we plot against , and the larger- value of the slope appears closer to than to the expected . However, this method relies heavily on the first half of the data, which corresponds to relatively small values and is far from asymptotic. Indeed, the pre-asymptotic portion of the best-fit curve in Fig. 2 is also well fitted by a straight line of slope , and we therefore consider the asymptotic regime not to have been reached in Fig. 3.

The situation for is less clear. Although the leading behaviour is again well verified, as shown by Fig. 4, the curve for the logarithmic correction is far from having reached its asymptote. The insert in Fig. 4 compares the last portion of the graph with the prediction . Whilst tending towards the correct value in the insert, the true asymptotic value remains to be determined. However, the graph does indicate that the logarithmic correction exponent is likely positive, and so supports more than .

## 4 Discussion

The scaling relation , introduced in Ref. [7], is essentially hyperscaling with the dimension replaced by for lattice animals corresponding to dimensional reduction arising in the mapping from lattice animals to the Yang-Lee model. In Ref. [17], a set of scaling relations for logarithmic corrections were developed, which included the corresponding hyperscaling relation in which is the correction exponent for the specific heat or free energy and is a logarithmic-correction exponent for the finite-size scaling of the correlation length.

In the case of lattice animals, may be identified with , from Eq.(13). Then, reducing the dimensionality appropriately, we find

(16) |

This is the logarithmic counterpart to the Parisi-Sourlas equation.

The value was proposed in Ref. [17] for lattice animals, also on the basis of scaling relations. Together with the estimates and , Eq.(16) holds in the present case. Indeed, the value fits the full set of scaling relations for logarithmic corrections proposed in Ref. [17]. On the other hand, does not satisfy the scaling relations when used in conjunction with the other exponent values known and reported in Ref. [11].

We have revisited the problem of lattice animals at the upper critical dimension and re-verified that the universal exponents and take their mean-field values there. We also provide numerical evidence in support of the (uncontested) logarithmic counterpart to the entopic index, and give an estimate for the growth non-universal constant . Regarding the logarithmic counterpart of the exponent, there are two candidate values in the literature. One of these is a direct constant-field calculation in Ref. [11], and the other is a Fisher renormalized version of the constant-order-parameter value also determined in Ref. [11]. Only the latter is consistent with the scaling relations, including a logarithmic version of the Parisi-Sourlas relation proposed herein. Our numerical evidence also indicates that the latter value is more likely to be the correct one, but, because of difficulties in achieving the asymptotic scaling regime, does not absolutely confirm the particular value.

## References

- [1] S. Golomb, Polyominoes: Puzzles, Patterns, Problems and Packings (Princeton Univ. Press, Princeton, N.J. 1994).
- [2] D. Stauffer and A. Aharony, Introduction to Percolation Theory, 2nd ed. (Taylor and Francis, London, 1992).
- [3] C.M. Fortuin and P.W. Kasteleyn, Physica 57, 536 (1972).
- [4] T. C. Lubensky and J. Isaacson, Phys. Rev. Lett. 41, 829 (1978); 42, 410(E) (1979); Phys. Rev. A 20, 2130 (1979).
- [5] H.-P. Hsu, W. Nadler and P. Grassberger, J. Phys. A: Math. Gen. 38, 775 (2005); Computer Physics Communications 169, 114 (2005).
- [6] T.C. Lubensky and J. Isaacson, Phys. Rev. Lett 41, 829 (1978); ibid 42, 410 (1979) (erratum).
- [7] G. Parisi and N. Sourlas, Phys. Rev. Lett. 46, 871 (1981).
- [8] F. Family, J. Phys. A 15, L583 (1982)
- [9] J. Adler, Y. Meir, A.B. Harris, A. Aharony and J.A.M.S. Duarté, Phys. Rev. B 38, 4941 (1988).
- [10] R. Kenna, H.-P Hsu and C. von Ferber, J. Stat. Mech. (JSTAT): Theor. Exp. (2008) L10002.
- [11] J.J. Ruiz-Lorenzo, J. Phys. A 31, 8773 (1998).
- [12] M.E. Fisher, Phys. Rev. Lett. 40, 1610 (1978).
- [13] D. Brydges and J.Z. Imbrie, Annals of Mathematics 158, 1019 (2003); J. Stat. Phys. 110, 503 (2003); J.Z. Imbrie, Annales Henri Poincaré 4, S445 (2003); J. Phys. A 37, L137 (2004).
- [14] T.C. Lubensky and A.J. McKane, J. Physique Lettres 42, L331 (1981).
- [15] M.E. Fisher, Phys. Rev. 176, 257 (1968).
- [16] C.M. Care, Phys. Rev. E 56, 1181 (1997); C.M. Care and R. Ettelaie, Phys. Rev. E 62, 1397 (2000).
- [17] R. Kenna, D.A. Johnston, and W. Janke, Phys. Rev. Lett. 96, 115701 (2006); ibid. 97, 155702 (2006).