Impact of Solvent Quality on the Hysteresis in the Coil–Stretch Transition of Flexible Polymers in Good Solvents

Conformations and Solvation of Synthetic Polymers in Water by Generalized Born Implicit-Solvent Molecular Dynamics Simulations: Stereoisomers of Poly(acrylic acid) and Poly(methacrylic acid)

We present the GB-OBC model as an approach for implicit-solvent MD simulations of a synthetic macromolecule in water. The model is tested and found to be successful in reproducing the chain dimensions and predicting the free energy of solvation of carboxylic acid vinyl polymers. The influence of stereochemistry and the hydrophobic nature of the polymer was investigated as a function of chain length (20 < N < 600) for poly(acrylic acid) (PAA) and poly(methacrylic acid) (PMA). The dimensionless parameters of the GB-OBC model were parameterized to be applicable to PAA and PMA. Scaling relations for chain dimensions obtained using implicit-solvent MD simulations in this study are in good agreement with those from experiments, theory of solvated chains in good solvents, and all-atom MD simulations in explicit water. Results show that ⟨R_g²⟩/NL² is greater for the atactic chain as compared to the isotactic chain, for PAA as well as PMA. ⟨R_g²⟩/NL² values of chains attain constancy in water for N = 200, with the values being greater for PMA. The PMA chain is conformationally more perturbed than the PAA chain, for both isotactic and atactic stereochemistry. The solvation free energy ΔG_hyd of PAA and PMA in water is negative for all chain lengths (N = 20-600) and becomes more favorable with an increase in molecular weight. The ΔG_hyd values for isotactic and atactic chains are identical at lower values of N but differ slightly for N > 300. Irrespective of the hydrophobic nature of the polymer, the atactic chain is thermodynamically more soluble in water as compared to the isotactic chain. The isotactic chain is less hydrophilic as compared to the atactic chain due to the closer proximity of the COOH groups along the backbone. This implicit solvent method is an effective way to accurately simulate the configurational properties and solvation of synthetic polymers in water.

Crooks Fluctuation Theorem for Single Polymer Dynamics in Time-Dependent Flows: Understanding Viscoelastic Hysteresis

Nonequilibrium work relations have fundamentally advanced our understanding of molecular processes. In recent years, fluctuation theorems have been extensively applied to understand transitions between equilibrium steady-states, commonly described by simple control parameters such as molecular extension of a protein or polymer chain stretched by an external force in a quiescent fluid. Despite recent progress, far less is understood regarding the application of fluctuation theorems to processes involving nonequilibrium steady-states such as those described by polymer stretching dynamics in nonequilibrium fluid flows. In this work, we apply the Crooks fluctuation theorem to understand the nonequilibrium thermodynamics of dilute polymer solutions in flow. We directly determine the nonequilibrium free energy for single polymer molecules in flow using a combination of single molecule experiments and Brownian dynamics simulations. We further develop a time-dependent extensional flow protocol that allows for probing viscoelastic hysteresis over a wide range of flow strengths. Using this framework, we define quantities that uniquely characterize the coil-stretch transition for polymer chains in flow. Overall, generalized fluctuation theorems provide a powerful framework to understand polymer dynamics under far-from-equilibrium conditions.

Characterization of hyaluronic acid and synovial fluid in stagnation point elongational flow

Hyaluronic acid (HA) is an important biomacromolecule, which fulfils a number of vital physiological functions (especially in the joint synovial fluid) and also has consumer and pharmaceutical applications. HA solution properties have already been quite thoroughly characterized in response to steady shear flows but are less well understood in highly deforming extensional flows. In this study, flow-induced birefringence measurements are made as a function of the strain rate in planar elongational flow at the stagnation point of a cross-slot device using HA solutions of a range of molecular weights (0.9×10(6) g mol(-1)≤Mw≤4.8×10(6) g mol(-1)) and at dilute concentrations. The results provide macromolecular relaxation times, molecular weight distributions and the extensional viscosities and Trouton ratios of the fluids. The HA relaxation time is found to vary as τ∼Mw1.8, which is consistent with a partially solvated, expanded coil. An intrinsic Trouton ratio is defined, which varies as [Tr]∼Mw2. The measurement of birefringence with strain rate is shown to be highly sensitive to the molecular weight distribution and can resolve subtle changes due to macromolecular degradation and the presence of fracture products. Mechanical degradation experiments in the cross-slots indicate midchain scission of HA macromolecules, strongly suggesting near full extension of the high-molecular weight fraction in the stagnation point extensional flow field. Taken together the results suggest a possible method for analysis of the HA in synovial fluid, and this concept is tested using synovial fluid obtained from porcine tarsal joint.

Coarse-grained model of conformation-dependent electrophoretic mobility and its influence on DNA dynamics

The electrophoretic mobility of molecules such as λ-DNA depends on the conformation of the molecule. It has been shown that electrohydrodynamic interactions between parts of the molecule lead to a mobility that depends on conformation and can explain some experimental observations. We have developed a new coarse-grained model that incorporates these changes of mobility into a bead-spring chain model. Brownian dynamics simulations have been performed using this model. The model reproduces the cross-stream migration that occurs in capillary electrophoresis when pressure-driven flow is applied parallel or antiparallel to the electric field. The model also reproduces the change of mobility when the molecule is stretched significantly in an extensional field. We find that the conformation-dependent mobility can lead to a new type of unraveling of the molecule in strong fields. This occurs when different parts of the molecule have different mobilities and the electric field is large.

Determining elasticity from single polymer dynamics

The ability to determine polymer elasticity and force-extension relations from polymer dynamics in flow has been challenging, mainly due to difficulties in relating equilibrium properties such as free energy to far-from-equilibrium processes. In this work, we determine polymer elasticity from the dynamic properties of polymer chains in fluid flow using recent advances in statistical mechanics. In this way, we obtain the force-extension relation for DNA from single molecule measurements of polymer dynamics in flow without the need for optical tweezers or bead tethers. We further employ simulations to demonstrate the practicality and applicability of this approach to the dynamics of complex fluids. We investigate the effects of flow type on this analysis method, and we develop scaling laws to relate the work relation to bulk polymer viscometric functions. Taken together, our results show that nonequilibrium work relations can play a key role in the analysis of soft material dynamics.

Fluctuations in the coil-stretch transition of flexible polymers in good solvents: a peak due to nonlinear force relation

Long flexible polymers undergo a coil to stretch transition (CST) in an elongational flow. Near the CST, a peak can be observed in the fluctuations of the size of a molecule (|R|). Solvent effects on the fluctuations are studied using Brownian dynamics simulations of a nonlinear spring force relation that can represent real molecules. Ignoring the influence of hydrodynamic interactions, a linear region in the spring force relation is known to cause the peak in |R| fluctuations. In contrast, we find that a peak in the fluctuations can be obtained even for the nonlinear spring force relation. We analyze the influence of hydrodynamic interactions on the fluctuations using a dumbbell model with a conformation-dependent drag coefficient.