Exact partition function zeros of a polymer on a simple cubic lattice

Polymers in Physics, Chemistry and Biology: Behavior of Linear Polymers in Fractal Structures

We start presenting an overview on recent applications of linear polymers and networks in condensed matter physics, chemistry and biology by briefly discussing selected papers (published within 2022-2024) in some detail. They are organized into three main subsections: polymers in physics (further subdivided into simulations of coarse-grained models and structural properties of materials), chemistry (quantum mechanical calculations, environmental issues and rheological properties of viscoelastic composites) and biology (macromolecules, proteins and biomedical applications). The core of the work is devoted to a review of theoretical aspects of linear polymers, with emphasis on self-avoiding walk (SAW) chains, in regular lattices and in both deterministic and random fractal structures. Values of critical exponents describing the structure of SAWs in different environments are updated whenever available. The case of random fractal structures is modeled by percolation clusters at criticality, and the issue of multifractality, which is typical of these complex systems, is illustrated. Applications of these models are suggested, and references to known results in the literature are provided. A detailed discussion of the reptation method and its many interesting applications are provided. The problem of protein folding and protein evolution are also considered, and the key issues and open questions are highlighted. We include an experimental section on polymers which introduces the most relevant aspects of linear polymers relevant to this work. The last two sections are dedicated to applications, one in materials science, such as fractal features of plasma-treated polymeric materials surfaces and the growth of polymer thin films, and a second one in biology, by considering among others long linear polymers, such as DNA, confined within a finite domain.

Universality class of the special adsorption point of two-dimensional lattice polymers

In recent work [Rodrigues et al., Phys. Rev. E 100, 022121 (2019)10.1103/PhysRevE.100.022121], evidence was found that the surface adsorption transition of interacting self-avoiding trails (ISATs) placed on the square lattice displays a nonuniversal behavior at the special adsorption point (SAP) where the collapsing polymers adsorb. In fact, different surface exponents ϕ^{(s)} and 1/δ^{(s)} were found at the SAP depending on whether the surface orientation is horizontal (HS) or diagonal (DS). Here, we revisit these systems and study other ones, through extensive Monte Carlo simulations, considering much longer trails than previous works. Importantly, we demonstrate that the different exponents observed in the reference above are due to the presence of a surface-attached-globule (SAG) phase in the DS system, which changes the multicritical nature of the SAP and is absent in the HS case. By considering a modified horizontal surface (mHS), on which the trails are forbidden from having two consecutive steps, resembling the DS situation, a stable SAG phase is found in the phase diagram, and both DS and mHS systems present similar 1/δ^{(s)} exponents at the SAP, namely, 1/δ^{(s)}≈0.44, whereas 1/δ^{(s)}≈0.34 in the HS case. Intriguingly, while ϕ^{(s)}≈1/δ^{(s)} is found for the DS and HS scenarios, as expected, in the mHS case ϕ^{(s)} is about 10% smaller than 1/δ^{(s)}. These results strongly indicate that at least two universality classes exist for the SAPs of adsorbing ISATs on the square lattice.

Solvent-Induced Negative Energetic Elasticity in a Lattice Polymer Chain

The negative internal energetic contribution to the elastic modulus (negative energetic elasticity) has been recently observed in polymer gels. This finding challenges the conventional notion that the elastic moduli of rubberlike materials are determined mainly by entropic elasticity. However, the microscopic origin of negative energetic elasticity has not yet been clarified. Here, we consider the n-step interacting self-avoiding walk on a cubic lattice as a model of a single polymer chain (a subchain of a network in a polymer gel) in a solvent. We theoretically demonstrate the emergence of negative energetic elasticity based on an exact enumeration up to n=20 and analytic expressions for arbitrary n in special cases. Furthermore, we demonstrate that the negative energetic elasticity of this model originates from the attractive polymer-solvent interaction, which locally stiffens the chain and conversely softens the stiffness of the entire chain. This model qualitatively reproduces the temperature dependence of negative energetic elasticity observed in the polymer-gel experiments, indicating that the analysis of a single chain can explain the properties of negative energetic elasticity in polymer gels.

Exact Enumeration Approach to Estimate the Theta Temperature of Interacting Self-Avoiding Walks on the Simple Cubic Lattice

We compute the exact root-mean-square end-to-end distance of the interacting self-avoiding walk (ISAW) up to 27 steps on the simple cubic lattice. These data are used to construct a fixed point equation to estimate the theta temperature of the collapse transition of the ISAW. With the Bulirsch-Stoer extrapolation method, we obtain accurate results that can be compared with large-scale long-chain simulations. The free parameter ω in extrapolation is precisely determined using a parity property of the ISAW. The systematic improvement of this approach is feasible by adopting the combination of exact enumeration and multicanonical simulations.

Partition-function-zero analysis of polymer adsorption for a continuum chain model

Polymer chains undergoing adsorption are expected to show universal critical behavior which may be investigated using partition function zeros. The focus of this work is the adsorption transition for a continuum chain, allowing for investigation of a continuous range of the attractive interaction and comparison with recent high-precision lattice model studies. The partition function (Fisher) zeros for a tangent-hard-sphere N-mer chain (monomer diameter σ) tethered to a flat wall with an attractive square-well potential (range λσ, depth ε) have been computed for chains up to N=1280 with 0.01≤λ≤2.0. In the complex-Boltzmann-factor plane these zeros are concentrated in an annular region, centered on the origin and open about the real axis. With increasing N, the leading zeros, w_{1}(N), approach the positive real axis as described by the asymptotic scaling law w_{1}(N)-y_{c}∼N^{-ϕ}, where y_{c}=e^{ε/k_{B}T_{c}} is the critical point and T_{c} is the critical temperature. In this work, we study the polymer adsorption transition by analyzing the trajectory of the leading zeros as they approach y_{c} in the complex plane. We use finite-size scaling (including corrections to scaling) to determine the critical point and the scaling exponent ϕ as well as the approach angle θ_{c}, between the real axis and the leading-zero trajectory. Variation of the interaction range λ moves the critical point, such that T_{c} decreases with λ, while the results for ϕ and θ_{c} are approximately independent of λ. Our values of ϕ=0.479(9) and θ_{c}=56.8(1.4)^{∘} are in agreement with the best lattice model results for polymer adsorption, further demonstrating the universality of these constants across both lattice and continuum models.

Transverse-size critical exponent of directed percolation from Yang-Lee zeros of survival probability

By using transfer-matrix method we compute survival probabilities for the directed percolation problem on strips of a square lattice, and get very precise estimates of their Yang-Lee zeros lying closest to the real axis in the complex plane of occupation probability. This allows us to get accurate values for transverse-size critical exponent and percolation threshold.

Physical mechanism for biopolymers to aggregate and maintain in non-equilibrium states

Many human or animal diseases are related to aggregation of proteins. A viable biological organism should maintain in non-equilibrium states. How protein aggregate and why biological organisms can maintain in non-equilibrium states are not well understood. As a first step to understand such complex systems problems, we consider simple model systems containing polymer chains and solvent particles. The strength of the spring to connect two neighboring monomers in a polymer chain is controlled by a parameter s with s → ∞ for rigid-bond. The strengths of bending and torsion angle dependent interactions are controlled by a parameter s A with s A  → -∞ corresponding to no bending and torsion angle dependent interactions. We find that for very small s A , polymer chains tend to aggregate spontaneously and the trend is independent of the strength of spring. For strong springs, the speed distribution of monomers in the parallel (along the direction of the spring to connect two neighboring monomers) and perpendicular directions have different effective temperatures and such systems are in non-equilibrium states.

Finite-size corrections and scaling for the dimer model on the checkerboard lattice

Lattice models are useful for understanding behaviors of interacting complex many-body systems. The lattice dimer model has been proposed to study the adsorption of diatomic molecules on a substrate. Here we analyze the partition function of the dimer model on a 2M×2N checkerboard lattice wrapped on a torus and derive the exact asymptotic expansion of the logarithm of the partition function. We find that the internal energy at the critical point is equal to zero. We also derive the exact finite-size corrections for the free energy, the internal energy, and the specific heat. Using the exact partition function and finite-size corrections for the dimer model on a finite checkerboard lattice, we obtain finite-size scaling functions for the free energy, the internal energy, and the specific heat of the dimer model. We investigate the properties of the specific heat near the critical point and find that the specific-heat pseudocritical point coincides with the critical point of the thermodynamic limit, which means that the specific-heat shift exponent λ is equal to ∞. We have also considered the limit N→∞ for which we obtain the expansion of the free energy for the dimer model on the infinitely long cylinder. From a finite-size analysis we have found that two conformal field theories with the central charges c=1 for the height function description and c=-2 for the construction using a mapping of spanning trees can be used to describe the dimer model on the checkerboard lattice.

Exact partition function zeros of the Wako-Saitô-Muñoz-Eaton β hairpin model

I compute exact partition function zeros of β hairpins, using both analytic and numerical methods, extending previous work [J. Lee, Phys. Rev. Lett. 110, 248101 (2013)] where only a restricted class of hairpins was considered. The zeros of β hairpins with an odd number of peptide bonds are computed and the difference of the distribution of zeros from those for an even number of peptide bonds is explained in terms of additional entropy of liberating the extra bond at the turn region. Upon the introduction of a hydrophobic core in the central region of the hairpin, the zeros are distributed uniformly on two concentric circles corresponding to the hydrophobic collapse and the transition to the fully folded conformation. One of the circles dissolves as the core moves toward the turn or the tip region, which is explained in terms of the similarity of the intermediate state with the folded or unfolded states. The exact partition function zeros for a hairpin with a more complex structure of native contacts, the 16 C-terminal residues of streptococcal protein G B1, are numerically computed and their loci are closely approximated by concentric circles.

Partition function zeros and phase transitions for a square-well polymer chain

The zeros of the canonical partition functions for flexible square-well polymer chains have been approximately computed for chains up to length 256 for a range of square-well diameters. We have previously shown that such chain molecules can undergo a coil-globule and globule-crystal transition as well as a direct coil-crystal transition. Here we show that each of these transitions has a well-defined signature in the complex-plane map of the partition function zeros. The freezing transitions are characterized by nearly circular rings of uniformly spaced roots, indicative of a discontinuous transition. The collapse transition is signaled by the appearance of an elliptical horseshoe segment of roots that pinches down towards the positive real axis and defines a boundary to a root-free region of the complex plane. With increasing chain length, the root density on the circular ring and in the space adjacent to the elliptical boundary increases and the leading roots move towards the positive real axis. For finite-length chains, transition temperatures can be obtained by locating the intersection of the ellipse and/or circle of roots with the positive real axis. A finite-size scaling analysis is used to obtain transition temperatures in the long-chain (thermodynamic) limit. The collapse transition is characterized by crossover and specific-heat exponents of φ≈0.76(2) and α≈0.66(2), respectively, consistent with a second-order phase transition.

Exact partition function zeros of the Wako-Saitô-Muñoz-Eaton protein model

I compute exact partition function zeros of the Wako-Saitô-Muñoz-Eaton model for various secondary structural elements and for two proteins, 1BBL and 1I6C, by using both analytic and numerical methods. Two-state and barrierless downhill folding transitions can be distinguished by a gap in the distribution of zeros at the positive real axis.

Partition function zeros of a square-lattice homopolymer with nearest- and next-nearest-neighbor interactions

We study distributions of the partition function zeros in the complex temperature plane for a square-lattice homopolymer with nearest-neighbor (NN) and next-nearest-neighbor (NNN) interactions. The dependence of distributions on the ratio of NN and NNN interaction strengths R is examined. The finite-size scaling of the zeros is performed to obtain the crossover exponent, which is shown to be independent of R within error bars, suggesting that all of these models belong to the same universality class. The transition temperatures are also computed by the zeros to obtain the phase diagram, and the results confirm that the model with stronger NNN interaction exhibits stronger effects of cooperativity.