We revisit an integrable lattice model of polymer collapse using numerical simulations. This model was first studied by Blöte and Nienhuis (1989 J. Phys. A: Math. Gen. 22 1415) and it describes polymers with some attraction, providing thus a model for the polymer collapse transition. At a particular set of Boltzmann weights the model is integrable and the exponents ν = 12/23 ≈ 0.522 and γ = 53/46 ≈ 1.152 have been computed via identification of the scaling dimensions x t = 1/12 and x h = −5/48. We directly investigate the polymer scaling exponents via Monte Carlo simulations using the pruned-enriched Rosenbluth method algorithm. By simulating this polymer model for walks up to length 4096 we find ν = 0.576(6) and γ = 1.045(5), which are clearly different from the predicted values. Our estimate for the exponent ν is compatible with the known θ-point value of 4/7 and in agreement with very recent numerical evaluation by Foster and Pinettes (2012 J. Phys. A: Math. Theor. 45 505003).
We revisit an integrable lattice model of polymer collapse using numerical simulations. This model was first studied by Blöte and Nienhuis (1989 J. Phys. A: Math. Gen. 22 1415) and it describes polymers with some attraction, providing thus a model for the polymer collapse transition. At a particular set of Boltzmann weights the model is integrable and the exponents ν = 12/23 ≈ 0.522 and γ = 53/46 ≈ 1.152 have been computed via identification of the scaling dimensions x t = 1/12 and x h = −5/48. We directly investigate the polymer scaling exponents via Monte Carlo simulations using the pruned-enriched Rosenbluth method algorithm. By simulating this polymer model for walks up to length 4096 we find ν = 0.576(6) and γ = 1.045(5), which are clearly different from the predicted values. Our estimate for the exponent ν is compatible with the known θ-point value of 4/7 and in agreement with very recent numerical evaluation by Foster and Pinettes (2012 J. Phys. A: Math. Theor. 45 505003).
We report a numerical investigation of the tricritical behaviour of the dilute equilibrium polymerization model with annealed dilution on the square lattice using transfer matrix techniques. In general, the tricritical exponent ν t is found to be compatible with the value for the tricritical n = 1 model. In the limit where the dilute model is equivalent to the interacting polymer model (whose tricritical point is known as the Θ-point), we observe a crossover of ν t to the tricritical n → 0 model value.
We study the collapse of two-dimensional polymers, via an O(n) model on the square lattice that allows for dilution, bending rigidity and short-range monomer attractions. This model contains two candidates for the theta point, ΘBN and ΘDS, both exactly solvable. The relative stability of these points, and the question of which one describes the `generic' theta point, have been the source of a long-standing debate. Moreover, the analytically predicted exponents of ΘBN have never been convincingly observed in numerical simulations. In the present paper, we shed a new light on this confusing situation. We show in particular that the continuum limit of ΘBN is an unusual conformal field theory, made in fact of a simple dense polymer decorated with {\sl non-compact degrees of freedom}. This implies in particular that the critical exponents take continuous rather than discrete values, and that corrections to scaling lead to an unusual integral form. Furthermore, discrete states may emerge from the continuum, but the latter are only normalizable---and hence observable---for appropriate values of the model's parameters. We check these findings numerically. We also probe the non-compact degrees of freedom in various ways, and establish that they are related to fluctuations of the density of monomers. Finally, we construct a field theoretic model of the vicinity of ΘBN and examine the flow along the multicritical line between ΘBN and ΘDS.
The nature of the θ point for a polymer in two dimensions has long been debated, with a variety of candidates put forward for the critical exponents. This includes those derived by Duplantier and Saleur for an exactly solvable model. We use a representation of the problem via the CP^{N-1}σ model in the limit N→1 to determine the stability of this critical point. First we prove that the Duplantier-Saleur (DS) critical exponents are robust, so long as the polymer does not cross itself: They can arise in a generic lattice model and do not require fine-tuning. This resolves a longstanding theoretical question. We also address an apparent paradox: Two different lattice models, apparently both in the DS universality class, show different numbers of relevant perturbations, apparently leading to contradictory conclusions about the stability of the DS exponents. We explain this in terms of subtle differences between the two models, one of which is fine-tuned (and not strictly in the DS universality class). Next we allow the polymer to cross itself, as appropriate, e.g., to the quasi-two-dimensional case. This introduces an additional independent relevant perturbation, so we do not expect the DS exponents to apply. The exponents in the case with crossings will be those of the generic tricritical O(n) model at n=0 and different from the case without crossings. We also discuss interesting features of the operator content of the CP^{N-1} model. Simple geometrical arguments show that two operators in this field theory, with very different symmetry properties, have the same scaling dimension for any value of N (or, equivalently, any value of the loop fugacity). Also we argue that for any value of N the CP^{N-1} model has a marginal odd-parity operator that is related to the winding angle.
By extending a previous analysis of bulk properties, accurate Monte Carlo enumerations of SAW with attractive n.n. interactions on the square lattice are used to obtain γ 1 t = 0.57 ± 0.09 and γ 11 t = -0.53 ± 0.10 for the entropic boundary exponents at the Θ-point. This gives new indications concerning the possible validity of recent conjectures about the Θ-point universality class. In particular, under plausible assumptions, these results appear compatible with a field-theoretic description in terms of a nonunitary superconformal field theory with central charge c = 1/2, for which γ 1 t = 15/28 and γ 11 t = -4/7 should be expected.
The nature of the theta point for a polymer in two dimensions has long been debated, with a variety of candidates put forward for the critical exponents. This includes those derived by Duplantier and Saleur (DS) for an exactly solvable model. We use a representation of the problem in terms of the CPN−1 sigma model in the limit N→1 to determine the stability of this critical point. First we prove that the DS critical exponents are robust so long as the polymer does not cross itself: they can arise in a generic lattice model, and do not require fine tuning. However there is an apparent paradox: two different lattice models, apparently both in the DS universality class, show different numbers of relevant perturbations, apparently leading to contradictory conclusions about the stability of the DS exponents. We explain this in terms of subtle differences between the two models, one of which is fine-tuned. Next, we allow polymer crossings, as appropriate e.g. to the quasi-2D case. This introduces an additional independent relevant perturbation, so we do not expect the DS exponents to apply: the exponents in the case with crossings will be those of the generic tricritical O(n) model at n=0, and different to the case without crossings. We also discuss interesting features of the operator content of the CPN−1 model. Simple geometrical arguments show that two operators in this field theory, with very different symmetry properties, have the same scaling dimension for any value of N. Also we argue that for any N there is a marginal parity-odd operator which is related to the winding angle.
Using a cluster-flipping Monte Carlo algorithm combined with a generalization of the histogram reweighting scheme of Ferrenberg and Swendsen [Phys. Rev. Lett. 61, 2635 (1988); 63, 1195 (1989)], we have studied the equilibrium properties of the thermal random-field Ising model on a cubic lattice in three dimensions. We have equilibrated systems of L\ifmmode\times\else\texttimes\fi{}L\ifmmode\times\else\texttimes\fi{}L spins, with values of L up to 32, and for these systems the cluster-flipping method appears to a large extent to overcome the slow equilibration seen in single-spin-flip methods. From the results of our simulations we have extracted values for the critical exponents and the critical temperature and randomness of the model by finite size scaling. For the exponents we find \ensuremath{\nu}=1.02\ifmmode\pm\else\textpm\fi{}0.06, \ensuremath{\beta}=0.06\ifmmode\pm\else\textpm\fi{}0.07, \ensuremath{\gamma}=1.9\ifmmode\pm\else\textpm\fi{}0.2, and \ensuremath{\gamma}\ifmmode\bar\else\textasciimacron\fi{}=2.9\ifmmode\pm\else\textpm\fi{}0.2, where \ensuremath{\nu}, \ensuremath{\beta}, \ensuremath{\gamma}, and \ensuremath{\gamma}\ifmmode\bar\else\textasciimacron\fi{} govern the critical singularities in the correlation length, magnetization, and connected and disconnected susceptibilities, respectively. \textcopyright{} 1996 The American Physical Society.