Thermodynamics of small systems with conformational transitions: The case of two-state freely jointed chains with extensible units

Stretching multistate flexible chains and loops

Polymer loop structure commonly appears in biological phenomena, such as DNA looping and DNA denaturation. When a chain forms a loop, its elastic behavior differs from that of an open chain due to the loss of entropy. In the case of reversible loop formation, interesting behavior emerges related to the multistate nature of the conformations. In this study, we model a multistate reversible loop as a looping Gaussian chain, which can bind (close) reversibly at one or several points to form a loop, or a zipping Gaussian loop, which can zip reversibly to form a double-stranded chain. For each model, we calculate the force-extension relations in the fixed-extension (Helmholtz) and the fixed-force (Gibbs) statistical ensembles. Unlike the single Gaussian chain or loop, the multilevel systems demonstrate qualitatively distinct tensile elasticity and ensemble inequivalence. In addition, we investigate a Gaussian necklace consisting of reversible alternating blocks of the zipped chain and loop and obtain the force-temperature phase diagram. The phase diagram implies a force-induced phase transition from a completely looped (denatured) state to a mixed (bound) state.

Forced extension of a wormlike chain in the Gibbs and Helmholtz ensembles

A semiflexible polymer can be stretched by either applying a force to it or by fixing the positions of its endpoints. The two approaches generally yield different results and correspond to experiments performed in either the Gibbs or Helmholtz statistical ensembles. Here, we derive the Helmholtz force-extension relationship for the commonly used wormlike-chain model in the strongly stretched regime. By analyzing it in comparison with the Gibbs ensemble result, we show that equivalence between the two relationships is achieved only in the long-chain thermodynamic limit.

Elasticity of a Grafted Rod-like Filament with Fluctuating Bending Stiffness

Quite often polymers exhibit different elastic behavior depending on the statistical ensemble (Gibbs vs. Helmholtz). This is an effect of strong fluctuations. In particular, two-state polymers, which locally or globally fluctuate between two classes of microstates, can exhibit strong ensemble inequivalence with negative elastic moduli (extensibility or compressibility) in the Helmholtz ensemble. Two-state polymers consisting of flexible beads and springs have been studied extensively. Recently, similar behavior was predicted in a strongly stretched wormlike chain consisting of a sequence of reversible blocks, fluctuating between two values of the bending stiffness (the so called reversible wormlike chain, rWLC). In this article, we theoretically analyse the elasticity of a grafted rod-like semiflexible filament which fluctuates between two states of bending stiffness. We consider the response to a point force at the fluctuating tip in both the Gibbs and the Helmholtz ensemble. We also calculate the entropic force exerted by the filament on a confining wall. This is done in the Helmholtz ensemble and, under certain conditions, it yields negative compressibility. We consider a two-state homopolymer and a two-block copolymer with two-state blocks. Possible physical realizations of such a system would be grafted DNA or carbon nanorods undergoing hybridization, or grafted F-actin bundles undergoing collective reversible unbinding.

Temperature dependent model for the quasi-static stick-slip process on a soft substrate

The classical Prandtl-Tomlinson model is the most famous and efficient method to describe the stick-slip phenomenon and the resulting friction between a slider and a corrugated substrate. It is widely used in all studies of frictional physics and notably in nanotribology. However, it considers a rigid or undeformable substrate and therefore is hardly applicable for investigating the physics of soft matter and in particular biophysics. For this reason, we introduce here a modified model that is capable of taking into consideration a soft or deformable substrate. It is realized by a sequence of elastically bound quadratic energy wells, which represent the corrugated substrate. We study the quasi-static behavior of the system through the equilibrium statistical mechanics. We thus determine the static friction and the deformation of the substrate as a function of temperature and substrate stiffness. The results are of interest for the study of cell motion in biophysics and for haptic and tactile systems in microtechnology.

Thermal effects on fracture and the brittle-to-ductile transition

The fracture behavior of brittle and ductile materials can be strongly influenced by thermal fluctuations, especially in micro- and nanodevices as well as in rubberlike and biological materials. However, temperature effects, in particular on the brittle-to-ductile transition, still require a deeper theoretical investigation. As a step in this direction we propose a theory, based on equilibrium statistical mechanics, able to describe the temperature-dependent brittle fracture and brittle-to-ductile transition in prototypical discrete systems consisting in a lattice with breakable elements. Concerning the brittle behavior, we obtain closed form expressions for the temperature-dependent fracture stress and strain, representing a generalized Griffith criterion, ultimately describing the fracture as a genuine phase transition. With regard to the brittle-to-ductile transition, we obtain a complex critical scenario characterized by a threshold temperature between the two fracture regimes (brittle and ductile), an upper and a lower yield strength, and a critical temperature corresponding to the complete breakdown. To show the effectiveness of the proposed models in describing thermal fracture behaviors at small scales, we successfully compare our theoretical results with molecular dynamics simulations of Si and GaN nanowires.

Freely jointed chain models with extensible links

Analytical relations for the mechanical response of single polymer chains are valuable for modeling purposes, on both the molecular and the continuum scale. These relations can be obtained using statistical thermodynamics and an idealized single-chain model, such as the freely jointed chain model. To include bond stretching, the rigid links in the freely jointed chain model can be made extensible, but this almost always renders the model analytically intractable. Here, an asymptotically correct statistical thermodynamic theory is used to develop analytic approximations for the single-chain mechanical response of this model. The accuracy of these approximations is demonstrated using several link potential energy functions. This approach can be applied to other single-chain models, and to molecular stretching in general.

Rate-dependent force-extension models for single-molecule force spectroscopy experiments

Single-molecule force spectroscopy techniques allow for the measurement of several static and dynamic features of macromolecules of biological origin. In particular, atomic force microscopy, used with a variable pulling rate, provides valuable information on the folding/unfolding dynamics of proteins. We propose here two different models able to describe the out-of-equilibrium statistical mechanics of a chain composed of bistable units. These latter represent the protein domains, which can be either folded or unfolded. Both models are based on the Langevin approach and their implementation allows for investigating the effect of the pulling rate and of the device intrinsic elasticity on the chain unfolding response. The theoretical results (both analytical and numerical) have been compared with experimental data concerning the unfolding of the titin and filamin proteins, eventually obtaining a good agreement over a large range of the pulling rates.

The influence of device handles in single-molecule experiments

We deduce a fully analytical model to predict the artifacts of the device handles in single molecule force spectroscopy experiments. As we show, neglecting the handle stiffness can lead to crucial overestimation or underestimation of the stability properties and unfolding thresholds of multistable molecules.

Full Statistics of Conjugated Thermodynamic Ensembles in Chains of Bistable Units

The statistical mechanics and the thermodynamics of small systems are characterized by the non-equivalence of the statistical ensembles. When concerning a polymer chain or an arbitrary chain of independent units, this concept leads to different force-extension responses for the isotensional (Gibbs) and the isometric (Helmholtz) thermodynamic ensembles for a limited number of units (far from the thermodynamic limit). While the average force-extension response has been largely investigated in both Gibbs and Helmholtz ensembles, the full statistical characterization of this thermo-mechanical behavior has not been approached by evaluating the corresponding probability densities. Therefore, we elaborate in this paper a technique for obtaining the probability density of the extension when force is applied (Gibbs ensemble) and the probability density of the force when the extension is prescribed (Helmholtz ensemble). This methodology, here developed at thermodynamic equilibrium, is applied to a specific chain composed of units characterized by a bistable potential energy, which is able to mimic the folding and unfolding of several macromolecules of biological origin.