Persistence-length renormalization of polymers in a crowded environment of hard disks

Bistability of a helical filament confined on a cylinder

The natural configuration of an intrinsically curved and twisted filament is uniquely a helix so that it can be referred to as a helical filament. We find that confining a helical filament on a cylinder can create a bistable state. When c_{0}R=0.5, where c_{0} is the intrinsic curvature of filament and R is the radius of cylinder, the phase diagram for the stability of a helix contains three regimes. Regime I has a small intrinsic twisting rate (ITR) and exhibits a bistable state which consists of two isoenergic helices. In regime II, the filament has a moderate ITR and the bistable state consists of a metastable low-pitch helix and a stable nonhelix. In regime III, the helix is unstable, owing to a large ITR. A similar phenomenon occurs when c_{0}R∼0.5. Monte Carlo simulation confirms these conclusions and indicates further that there are bistable nonhelices in regime III. This bistable system offers a prospective green material since the wide range of parameters and distinctive configurations for bistable states favor its realization and application.

Elastic behavior of a semiflexible polymer in 3D subject to compression and stretching forces

We elucidate the elastic behavior of a wormlike chain in 3D under compression and provide exact solutions for the experimentally accessible force-extension relation in terms of generalized spheroidal wave functions. In striking contrast to the classical Euler buckling instability, the force-extension relation of a clamped semiflexible polymer exhibits a smooth crossover from an almost stretched to a buckled configuration. In particular, the associated susceptibility, which measures the strength of the response of the polymer to the applied force, displays a prominent peak in the vicinity of the critical Euler buckling force. For increasing persistence length, the force-extension relation and the susceptibility of semiflexible polymers approach the behavior of a classical rod, whereas thermal fluctuations permit more flexible polymers to resist the applied force. Furthermore, we find that semiflexible polymers confined to 2D can oppose the applied force more strongly than in 3D.

Bimodal probability density characterizes the elastic behavior of a semiflexible polymer in 2D under compression

We explore the elastic behavior of a wormlike chain under compression in terms of exact solutions for the associated probability densities. Strikingly, the probability density for the end-to-end distance projected along the applied force exhibits a bimodal shape in the vicinity of the critical Euler buckling force of an elastic rod, reminiscent of the smeared discontinuous phase transition of a finite system. These two modes reflect the almost stretched and the S-shaped configuration of a clamped polymer induced by the compression. Moreover, we find a bimodal shape of the probability density for the transverse fluctuations of the free end of a cantilevered polymer as fingerprint of its semiflexibility. In contrast to clamped polymers, free polymers display a circularly symmetric probability density and their distributions are identical for compression and stretching forces.

Exact solution for the force-extension relation of a semiflexible polymer under compression

Exact solutions for the elastic and thermodynamic properties for the wormlike chain model are elaborated in terms of Mathieu functions. The smearing of the classical Euler buckling instability for clamped polymers is analyzed for the force-extension relation. Interestingly, at strong compression forces the thermal fluctuations lead to larger elongations than for the elastic rod. The susceptibility defined as the derivative of the force-extension relation displays a prominent maximum at a force that approaches the critical Euler buckling force as the persistence length is increased. We also evaluate the excess entropy and heat capacity induced by the compression and find that they vary nonmonotonically with the load. These findings are corroborated by pseudo-Brownian simulations.

Active Brownian particles moving in a random Lorentz gas

Biological microswimmers often inhabit a porous or crowded environment such as soil. In order to understand how such a complex environment influences their spreading, we numerically study non-interacting active Brownian particles (ABPs) in a two-dimensional random Lorentz gas. Close to the percolation transition in the Lorentz gas, they perform the same subdiffusive motion as ballistic and diffusive particles. However, due to their persistent motion they reach their long-time dynamics faster than passive particles and also show superdiffusive motion at intermediate times. While above the critical obstacle density [Formula: see text] the ABPs are trapped, their long-time diffusion below [Formula: see text] is strongly influenced by the propulsion speed v₀. With increasing v₀, ABPs are stuck at the obstacles for longer times. Thus, for large propulsion speed, the long-time diffusion constant decreases more strongly in a denser obstacle environment than for passive particles. This agrees with the behavior of an effective swimming velocity and persistence time, which we extract from the velocity autocorrelation function.

Microstructure of Sheared Entangled Solutions of Semiflexible Polymers

We study the influence of finite shear deformations on the microstructure and rheology of solutions of entangled semiflexible polymers theoretically and by numerical simulations and experiments with filamentous actin. Based on the tube model of semiflexible polymers, we predict that large finite shear deformations strongly affect the average tube width and curvature, thereby exciting considerable restoring stresses. In contrast, the associated shear alignment is moderate, with little impact on the average tube parameters, and thus expected to be long-lived and detectable after cessation of shear. Similarly, topologically preserved hairpin configurations are predicted to leave a long-lived fingerprint in the shape of the distributions of tube widths and curvatures. Our numerical and experimental data support the theory.

Thermodynamics and structure of macromolecules from flat-histogram Monte Carlo simulations

Over the last decade flat-histogram Monte Carlo simulations, especially multi-canonical and Wang-Landau simulations, have emerged as a strong tool to study the statistical mechanics of polymer chains. These investigations have focused on coarse-grained models of polymers on the lattice and in the continuum. Phase diagrams of chains in bulk as well as chains attached to surfaces were studied, for homopolymers as well as for protein-like models. Also, aggregation behavior in solution of these models has been investigated. We will present here the theoretical background for these simulations, explain the algorithms used and discuss their performance and give an overview over the systems studied with these methods in the literature, where we will limit ourselves to studies of coarse-grained model systems. Implementations of these algorithms on parallel computers will be also briefly described. In parallel to the development of these simulation methods, the power of a micro-canonical analysis of such simulations has been recognized, and we present the current state of the art in applying the micro-canonical analysis to phase transitions in nanoscopic polymer systems.