Nature of the collapse transition in interacting self-avoiding trails

We study the interacting self-avoiding trail (ISAT) model on a Bethe lattice of general coordination q and on a Husimi lattice built with squares and coordination q=4. The exact grand-canonical solutions of the model are obtained, considering that up to K monomers can be placed on a site and associating a weight ω_{i} with an i-fold visited site. Very rich phase diagrams are found with nonpolymerized, regular polymerized, and dense polymerized phases separated by lines (or surfaces) of continuous and discontinuous transitions. For a Bethe lattice with q=4 and K=2, the collapse transition is identified with a bicritical point and the collapsed phase is associated with the dense polymerized (solidlike) phase instead of the regular polymerized (liquidlike) phase. A similar result is found for the Husimi lattice, which may explain the difference between the collapse transition for ISATs and for interacting self-avoiding walks on the square lattice. For q=6 and K=3 (studied on the Bethe lattice only), a more complex phase diagram is found, with two critical planes and two coexistence surfaces, separated by two tricritical and two critical end-point lines meeting at a multicritical point. The mapping of the phase diagrams in the canonical ensemble is discussed and compared with simulational results for regular lattices.

Transfer-matrix study of a hard-square lattice gas with two kinds of particles and density anomaly

Using transfer matrix and finite-size scaling methods, we study the thermodynamic behavior of a lattice gas with two kinds of particles on the square lattice. Only excluded volume interactions are considered, so that the model is athermal. Large particles exclude the site they occupy and its four first neighbors, while small particles exclude only their site. Two thermodynamic phases are found: a disordered phase where large particles occupy both sublattices with the same probability and an ordered phase where one of the two sublattices is preferentially occupied by them. The transition between these phases is continuous at small concentrations of the small particles and discontinuous at larger concentrations, both transitions are separated by a tricritical point. Estimates of the central charge suggest that the critical line is in the Ising universality class, while the tricritical point has tricritical Ising (Blume-Emery-Griffiths) exponents. The isobaric curves of the total density as functions of the fugacity of small or large particles display a minimum in the disordered phase.

Large scale behavior of a two-dimensional model of anisotropic branched polymers

We study critical properties of anisotropic branched polymers modeled by semi-directed lattice animals on a triangular lattice. Using the exact transfer-matrix approach on strips of quite large widths and phenomenological renormalization group analysis, we obtained pretty good estimates of various critical exponents in the whole high-temperature region, including the point of collapse transition. Our numerical results suggest that this collapse transition belongs to the universality class of directed percolation.

Collapse transition of randomly branched polymers: renormalized field theory

We present a minimal dynamical model for randomly branched isotropic polymers, and we study this model in the framework of renormalized field theory. For the swollen phase, we show that our model provides a route to understand the well-established dimensional-reduction results from a different angle. For the collapse θ transition, we uncover a hidden Becchi-Rouet-Stora supersymmetry, signaling the sole relevance of tree configurations. We correct the long-standing one-loop results for the critical exponents, and we push these results on to two-loop order. For the collapse θ' transition, we find a runaway of the renormalization group flow, which lends credence to the possibility that this transition is a fluctuation-induced first-order transition. Our dynamical model allows us to calculate for the first time the fractal dimension of the shortest path on randomly branched polymers in the swollen phase as well as at the collapse transition and related fractal dimensions.

Renormalized field theory of collapsing directed randomly branched polymers

We present a dynamical field theory for directed randomly branched polymers and in particular their collapse transition. We develop a phenomenological model in the form of a stochastic response functional that allows us to address several interesting problems such as the scaling behavior of the swollen phase and the collapse transition. For the swollen phase, we find that by choosing model parameters appropriately, our stochastic functional reduces to the one describing the relaxation dynamics near the Yang-Lee singularity edge. This corroborates that the scaling behavior of swollen branched polymers is governed by the Yang-Lee universality class as has been known for a long time. The main focus of our paper lies on the collapse transition of directed branched polymers. We show to arbitrary order in renormalized perturbation theory with epsilon expansion that this transition belongs to the same universality class as directed percolation.

Thermodynamic behavior of a polymer with interacting bonds on a square lattice

Using the transfer matrix technique, finite-size scaling, phenomenological renormalization group, and conformal invariance ideas, the thermodynamic behavior of a polymer with interacting bonds on a square lattice has been studied. In this model, one monomer that belongs to the polymer has an activity x=e(beta(mu)), while the interactions between bonds of the polymer that are located on opposite edges of elementary squares of the lattice have a statistical weight y=e(-beta(epsilon)), where epsilon is the interaction energy. Next, the phase diagram of the model in the (x,y) plane was found, which shows three phases, two of them being polymerized. Furthermore, the densities of occupied sites and of bond interactions in each phase were calculated, in order to determine the nature of the transitions between the phases. The results obtained are consistent with a second-order transition line between the nonpolymerized and the regular polymerized phase and a first-order transition between the nonpolymerized and the dense polymerized phase. The boundary between both polymerized phases may be of first or second order, and thus evidence for a tricritical point is found.

Simple model for the DNA denaturation transition

We study pairs of interacting self-avoiding walks ¿omega(1), omega(2)¿ on the 3d simple cubic lattice. They have a common origin omega(1)(0)=omega 20, and are allowed to overlap only at the same monomer position along the chain: omega(1)(i) not equal omega(2)(j) for i not equal j, while omega(1)(i)=omega(2)(i) is allowed. The latter overlaps are indeed favored by an energetic gain epsilon. This is inspired by a model introduced long ago by Poland and Sheraga [J. Chem. Phys. 45, 1464 (1966)] for the denaturation transition in DNA where, however, self avoidance was not fully taken into account. For both models, there exists a temperature T(m) above which the entropic advantage to open up overcomes the energy gained by forming tightly bound two-stranded structures. Numerical simulations of our model indicate that the transition is of first order (the energy density is discontinuous), but the analog of the surface tension vanishes and the scaling laws near the transition point are exactly those of a second-order transition with crossover exponent straight phi=1. Numerical and exact analytic results show that the transition is second order in modified models where the self-avoidance is partially or completely neglected.

Side-chain entropy and packing in proteins

What role does side-chain packing play in protein stability and structure? To address this question, we compare a lattice model with side chains (SCM) to a linear lattice model without side chains (LCM). Self-avoiding configurations are enumerated in 2 and 3 dimensions exhaustively for short chains and by Monte Carlo sampling for chains up to 50 main-chain monomers long. This comparison shows that (1) side-chain degrees of freedom increase the entropy of open conformations, but side-chain steric exclusion decreases the entropy of compact conformations, thus producing a substantial entropy that opposes folding; (2) there is side-chain "freezing" or ordering, i.e., a sharp decrease in entropy, near maximum compactness; and (3) the different types of contacts among side chains (s) and main-chain elements (m) have different frequencies, and the frequencies have different dependencies on compactness. mm contacts contribute significantly only at high densities, suggesting that main-chain hydrogen bonding in proteins may be promoted by compactness. The distributions of mm, ms, and ss contacts in compact SCM configurations are similar to the distributions in protein structures in the Brookhaven Protein Data Bank. We propose that packing in proteins is more like the packing of nuts and bolts in a jar than like the pairwise matching of jigsaw puzzle pieces.