Distribution functions and dynamical properties of stiff macromolecules

An analytically tractable model for chain molecules with bending stiffness is presented and the dynamical properties of such chains are investigated. The partition function is derived via the maximum entropy principle taking into account the chain connectivity as well as the bending restrictions in form of constraints. We demonstrate that second moments agree exactly with those known from the Kratky-Porod wormlike chain. Moreover, various distribution functions are calculated. In particular, the static structure factor is shown to be proportional to 1/q at large scattering vectors q. The equations of motion for a chain in a melt as well as in dilute solution are presented. In the latter case the hydrodynamic interaction is taken into account via the Rotne-Prager tensor. The dynamical equations are solved by a normal mode analysis. In the limit of a flexible chain the model reproduces the well-known Rouse and Zimm dynamics, respectively, on large length scales, whereas in the rod limit the eigenfunctions correspond to bending motion only. In addition, the coherent and incoherent dynamic structure factor is discussed. For melts we show that at large scattering vectors the incoherent dynamic structure factor is a universal function of only the combination q^8/3tp^1/3, where 1/(2p) is the persistence length of the macromolecules. The comparison of the theoretical results with quasielastic neutron and light scattering experiments of various polymers in solution and melt exhibits good agreement. Our investigations show that local stiffness strongly influences the dynamics of macromolecules on small length scales even for long and flexible chains.     

Kinetics of interior loop formation in semiflexible chains

Loop formation between monomers in the interior of semiflexible chains describes elementary events in biomolecular folding and DNA bending. We calculate analytically the interior distance distribution function for semiflexible chains using a mean field approach. Using the potential of mean force derived from the distance distribution function we present a simple expression for the kinetics of interior looping by adopting Kramers theory. For the parameters, that are appropriate for DNA, the theoretical predictions in comparison with the case are in excellent agreement with explicit Brownian dynamics simulations of wormlike chain (WLC) model. The interior looping times (tauIC) can be greatly altered in the cases when the stiffness of the loop differs from that of the dangling ends. If the dangling end is stiffer than the loop then tauIC increases for the case of the WLC with uniform persistence length. In contrast, attachment of flexible dangling ends enhances rate of interior loop formation. The theory also shows that if the monomers are charged and interact via screened Coulomb potential then both the cyclization (tauc) and interior looping (tauIC) times greatly increase at low ionic concentration. Because both tauc and tauIC are determined essentially by the effective persistence length [lp(R)] we computed lp(R) by varying the range of the repulsive interaction between the monomers. For short range interactions lp(R) nearly coincides with the bare persistence length which is determined largely by the backbone chain connectivity. This finding rationalizes the efficacy of describing a number of experimental observations (response of biopolymers to force and cyclization kinetics) in biomolecules using WLC model with an effective persistence length.     

Statistics of time-limited ensembles of bent DNA conformations

The paper considers statistical properties of ensembles of chain conformations obtained by short-time Brownian dynamics (BD) of a coarse-grained DNA model in order to find out if the conditions necessary for accurate evaluation of the polymer elasticity are attainable in atom-level molecular dynamics (MD) simulations. To measure the bending persistence length (PL) with a 10% error using data accumulated in a single trajectory of a double helix of 15 base pairs, dynamics should be continued for a few microseconds. However, these estimates should be scaled down by about 2 orders of magnitude because the bending dynamics of short double helices in MD features much smaller relaxation times. As a result, good qualitative agreement with the worm-like chain (WLC) theory is reached in MD after tens of nanoseconds. The presently accessible durations of MD trajectories provide reasonably accurate evaluation of DNA elasticity and allow modeling of its mesoscopic properties. The surprisingly fast bending dynamics of short double helices in MD suggests that the microscopic mechanisms of DNA flexibility differ from a simple harmonic model.     

Diffusion and segmental dynamics of rodlike molecules by fluorescence correlation spectroscopy

The dynamics of weakly bending polymers is analyzed on the basis of a Gaussian semiflexible chain model and the fluorescence correlation spectroscopy (FCS) correlation function is determined. Particular attention is paid to the influence of the rotational motion on the decay of the FCS correlation function. An analytical expression for the correlation function is derived, from which the averaged segmental mean square displacement can be determined independent of any specific model for the polymer dynamics. The theoretical analysis exhibits a strong dependence of the correlation function on the rotational motion for semiflexible polymers with typical lengths and persistence lengths of actin filaments or fd viruses. Hence, FCS allows for a measurement of the rotational motion of such semiflexible polymers. The theoretical results agree well with experimental measurements on actin filaments and confirm the importance of large relaxation times.     

Interfaces in polymer blends

We investigate the structure and thermodynamics of interfaces in dense polymer blends using Monte Carlo (MC) simulations and self‐consistent field (SCF) calculations. For structurally symmetric blends we find quantitative agreement between the MC simulations and the SCF calculations for excess quantities of the interface (e.g., interfacial tension or enrichment of copolymers at the interface). However, a quantitative comparison between profiles across the interface in the MC simulations and the SCF calculations has to take due account of capillary waves. While the profiles in the SCF calculations correspond to intrinsic profiles of a perfectly flat interface the local interfacial position fluctuates in the MC simulations. We test this concept by extensive Monte Carlo simulations and study the cross‐over between “intrinsic” fluctuations which build up the local profile and capillary waves on long (lateral) length scales. Properties of structurally asymmetric blends are exemplified by investigating polymers of different stiffness. At high incompatibilities the interfacial width is not much larger than the persistence length of the stiffer component. In this limit we find deviations from the predictions of the Gaussian chain model: while the Gaussian chain model yields an increase of the interfacial width upon increasing the persistence length, no such increase is found in the MC simulations. Using a partial enumeration technique, however, we can account for the details of the chain architecture on all length scales in the SCF calculations and achieve good agreement with the MC simulations. In blends containing diblock copolymers we investigate the enrichment of copolymers at the interface and the concomitant reduction of the interfacial tension. At weak segregation the addition of copolymers leads to compatibilization. At high incompatibilities, the homopolymer‐rich phase can accommodate only a small fraction of copolymer before the copolymer forms a lamellar phase. The analysis of interfacial fluctuations yields an estimate for the bending rigidity of the interface. The latter quantity is important for the formation of a polymeric microemulsion at intermediate segregation and the consequences for the phase diagram are discussed.     

Model theory of relaxation properties of macromolecules in the nematic LC state

The continuous and lattice dynamic models of the macromolecules in the LC state are considered. The relaxation properties of polymer chain in the nematic LC state may depend on both the variation of the conformation in the ordered state and on the activation barrier of the molecular field. The calculation and comparison of the longitudinal and transversal relaxation spectra for the continuous and discrete rotameric mechanism of the mobility are performed. For the simplest model of the heterogeneous polymer chain for continuous dynamic models the complicated relaxational behavior with two longitudinal and two transversal branches is manifested. The dynamic properties of the lattice dynamic model with spacer consisting of several bonds are considered.     

Heterogeneous thermal excitation and relaxation in supercooled liquids

We investigate a phenomenological model which rationalizes the effects of dielectric hole burning on the basis of heterogeneous dielectric and specific heat relaxation in supercooled liquids. The quantitative agreement between model predictions and dielectric hole-burning observations is lost if the assumption of correlated dielectric and thermal relaxation times is removed from the model. This suggests that dynamically distinct domains in real liquids are associated with a time constant which characterizes both the structural and thermal relaxation behaviors. The calculations demonstrate that the observed burn-induced modifications reflect the spectral selectivity and persistence time of the fictive temperatures within these domains, and that 100 or more cycles of the sinusoidal burn field can be required to saturate the heat accumulated in the slow degrees of freedom. It is also shown that the recovery of dielectric holes is entirely accounted for by the model, and that the persistence times do not provide direct insight into rate exchange processes. Additionally, the model predicts that the heating effects considered here are a significant source of nonlinear dielectric behavior, even in the absence of deliberate frequency selective hole burning.     

Influence of molecular stiffness on the dynamic structure factor

Characteristic features of the influence of molecular stiffness on the dynamic structure factor of macromolecules are briefly outlined. The relaxation times characterizing the internal dynamics of the macro‐molecules exhibit a crossover from Rouse‐Zimm to bending modes with increasing mode number. As a consequence the dynamic structure factor is strongly influenced by the molecular stiffness. In particular, a stretched exponential relaxation of the dynamic structure factor at scattering vectors larger than the inverse persistence length is predicted and confirmed by a comparison with experimental data. Moreover, the influence of polydispersity is discussed.     

Quantum two-state dynamics driven by stationary non-Markovian discrete noise: Exact results

We consider the problem of stochastic averaging of a quantum two-state dynamics driven by non-Markovian, discrete noises of the continuous time random walk type (multistate renewal processes). The emphasis is put on the proper averaging over the stationary noise realizations corresponding, e.g., to a stationary environment. A two-state non-Markovian process with an arbitrary non-exponential distribution of residence times (RTDs) in its states with a finite mean residence time provides a paradigm. For the case of a two-state quantum relaxation caused by such a classical stochastic field we obtain the explicit exact, analytical expression for the averaged Laplace-transformed relaxation dynamics. In the limit of Markovian noise (implying an exponential RTD), all previously known results are recovered. We exemplify new more general results for the case of non-Markovian noise with a biexponential RTD. The averaged, real-time relaxation dynamics is obtained in this case by numerically exact solving of a resulting algebraic polynomial problem. Moreover, the case of manifest non-Markovian noise with an infinite range of temporal autocorrelation (which in principle is not accessible to any kind of perturbative treatment) is studied, both analytically (asymptotic long-time dynamics) and numerically (by a precise numerical inversion of the Laplace-transformed averaged quantum relaxation).     

Rotational correlation and dynamic heterogeneity in a kinetically constrained lattice gas

We study the dynamical heterogeneity and glassy dynamics in a kinetically constrained lattice-gas model which has both translational and rotational degrees of freedom. We find that the rotational relaxation time tracks the structural relaxation time as density is increased whereas the translational diffusion constant exhibits a strong decoupling. We investigate distributions of exchange and persistence times for both the rotational and translational degrees of freedom and compare our results on the distributions of rotational exchange times to recent single molecule studies.     

Modeling diffusional transport in the interphase cell nucleus

In this paper a lattice model for the diffusional transport of particles in the interphase cell nucleus is proposed. Dense networks of chromatin fibers are created by three different methods: Randomly distributed, noninterconnected obstacles, a random walk chain model, and a self-avoiding random walk chain model with persistence length. By comparing a discrete and a continuous version of the random walk chain model, we demonstrate that lattice discretization does not alter particle diffusion. The influence of the three dimensional geometry of the fiber network on the particle diffusion is investigated in detail while varying the occupation volume, chain length, persistence length, and walker size. It is shown that adjacency of the monomers, the excluded volume effect incorporated in the self-avoiding random walk model, and, to a lesser extent, the persistence length affect particle diffusion. It is demonstrated how the introduction of the effective chain occupancy, which is a convolution of the geometric chain volume with the walker size, eliminates the conformational effects of the network on the diffusion, i.e., when plotting the diffusion coefficient as a function of the effective chain volume, the data fall onto a master curve.     

Computer simulations of polymer chain structure and dynamics on a hypersphere in four-space

There is a rapidly growing interest in performing computer simulations in a closed space, avoiding periodic boundary conditions. To extend the range of potential systems to include also macromolecules, we describe an algorithm for computer simulations of polymer chain molecules on S3, a hypersphere in four dimensions. In particular, we show how to generate initial conformations with a bond angle distribution given by the persistence length of the chain and how to calculate the bending forces for a molecule moving on S3. Furthermore, we discuss how to describe the shape of a macromolecule on S3, by deriving the radius of gyration tensor in this non-Euclidean space. The results from both Monte Carlo and Brownian dynamics simulations in the infinite dilution limit show that the results on S3 and in R3 coincide, both with respect to the size and shape as well as for the diffusion coefficient. All data on S3 can also be described by master curves by suitable scaling by the corresponding values in R3. We thus show how to extend the use of spherical boundary conditions, which are most effective for calculating electrostatic forces, to polymer chain molecules, making it possible to perform simulations on S3 also for polyelectrolyte systems.     

The role of chromatin conformations in diffusional transport of chromatin-binding proteins: Cartesian lattice simulations

In this paper, a lattice model for the diffusional transport of chromatin-binding particles in the interphase cell nucleus is proposed. Sliding effects are studied in dense networks of chromatin fibers created by three different methods: Randomly distributed, noninterconnected obstacles, a random walk chain model with an attractive step potential, and a self-avoiding random walk chain model with a hard repulsive core and attractive surroundings. By comparing a discrete and continuous version of the random walk chain model, we demonstrate that lattice discretization does not alter the diffusion of chromatin-binding particles. The influence of conformational properties of the fiber network on the particle sliding is investigated in detail while varying occupation volume, sliding probability, chain length, and persistence length. It is observed that adjacency of the monomers, the excluded volume effect incorporated in the self-avoiding random walk model, and the persistence length affect the chromatin-binding particle diffusion. It is demonstrated that sliding particles sense local chain structures. When plotting the diffusion coefficient as a function of the accessible volume for diffusing particles, the data fall onto master curves depending on the persistence length. However, once intersegment transfer is involved, chromatin-binding proteins no longer perceive local chain structures.     

Universal properties of a single polymer chain in slit: Scaling versus molecular dynamics simulations

We revisit the classical problem of a polymer confined in a slit in both of its static and dynamic aspects. We confirm a number of well known scaling predictions and analyze their range of validity by means of comprehensive molecular dynamics simulations using a coarse-grained bead-spring model of a flexible polymer chain. The normal and parallel components of the average end-to-end distance, mean radius of gyration and their distributions, the density profile, the force exerted on the slit walls, and the local bond orientation characteristics are obtained in slits of width D=4/10 (in units of the bead diameter) and for chain lengths N=50/300. We demonstrate that a wide range of static chain properties in normal direction can be described quantitatively by analytic model-independent expressions in perfect agreement with computer experiment. In particular, the observed profile of confinement-induced bond orientation is shown to closely match theory predictions. The anisotropy of confinement is found to be manifested most dramatically in the dynamic behavior of the polymer chain. We examine the relation between characteristic times for translational diffusion and lateral relaxation. It is demonstrated that the scaling predictions for lateral and normal relaxation times are in good agreement with our observations. A novel feature is the observed coupling of normal and lateral modes with two vastly different relaxation times. We show that the impact of grafting on lateral relaxation is equivalent to doubling the chain length.     

From caging to Rouse dynamics in polymer melts with intramolecular barriers: a critical test of the mode coupling theory

By means of computer simulations and solution of the equations of the mode coupling theory (MCT), we investigate the role of the intramolecular barriers on several dynamic aspects of nonentangled polymers. The investigated dynamic range extends from the caging regime characteristic of glass-formers to the relaxation of the chain Rouse modes. We review our recent work on this question, provide new results, and critically discuss the limitations of the theory. Solutions of the MCT for the structural relaxation reproduce qualitative trends of simulations for weak and moderate barriers. However, a progressive discrepancy is revealed as the limit of stiff chains is approached. This disagreement does not seem related with dynamic heterogeneities, which indeed are not enhanced by increasing barrier strength. It is not connected either with the breakdown of the convolution approximation for three-point static correlations, which retains its validity for stiff chains. These findings suggest the need of an improvement of the MCT equations for polymer melts. Concerning the relaxation of the chain degrees of freedom, MCT provides a microscopic basis for time scales from chain reorientation down to the caging regime. It rationalizes, from first principles, the observed deviations from the Rouse model on increasing the barrier strength. These include anomalous scaling of relaxation times, long-time plateaux, and nonmonotonous wavelength dependence of the mode correlators.     

Rotational relaxation in simple chain models

The rotational dynamics of chemically similar systems based on freely jointed and freely rotating chains are studied. The second Legendre polynomial of vectors along chain backbones is used to investigate the rotational dynamics at different length scales. In a previous study, it was demonstrated that the additional bond-angle constraint in the freely rotating case noticeably perturbs the character of the translational relaxation away from that of the freely jointed system. Here, it is shown that differences are also apparent in the two systems' rotational dynamics. The relaxation of the end-to-end vector is found to display a long time, single-exponential tail and a stretched exponential region at intermediate times. The stretching exponents beta are found to be 0.75+/-0.02 for the freely jointed case and 0.68+/-0.02 for the freely rotating case. For both system types, time-packing-fraction superposition is seen to hold on the end-to-end length scale. In addition, for both systems, the rotational relaxation times are shown to be proportional to the translational relaxation times, demonstrating that the Debye-Stokes-Einstein law holds. The second Legendre polynomial of the bond vector is used to probe relaxation behavior at short length scales. For the freely rotating case, the end-to-end relaxation times scale differently than the bond relaxation times, implying that the behavior is non-Stokes-Einstein, and that time-packing-fraction superposition does not hold across length scales for this system. For the freely jointed case, end-to-end relaxation times do scale with bond relaxation times, and both Stokes-Einstein and time-packing-fraction-across-length-scales superposition are obeyed.     

Hyper-aging dynamics of nanoclay suspension

Aqueous suspension of nanoclay Laponite undergoes structural evolution as a function of time, which enhances its elasticity and relaxation time. In this work, we employ an effective time approach to investigate long-term relaxation dynamics by carrying out creep experiments. Typically, we observe that the monotonic evolution of elastic modulus shifts to lower aging times, while maxima in viscous moduli get progressively broader for experiments carried out on a later date after preparation (idle time) of the nanoclay suspension. Application of effective time theory produces a superposition of all the creep curves irrespective of their initial state. The resulting dependence of the relaxation time on aging time shows very strong hyper-aging dynamics at short idle times, which progressively weakens to demonstrate a linear dependence in the limit of very long idle times. Remarkably, this behavior of nanoclay suspensions is akin to that observed for polymeric glasses. Consideration of aging as a first-order process suggests that continued hyper-aging dynamics causes cessation of aging. The dependence of relaxation time on aging time, therefore, must attenuate eventually producing linear or weaker dependence on time in order to approach a progressively low-energy state in the limit of very long times as observed experimentally. We also develop a simple scaling model based on a concept of aging of an energy well, which qualitatively captures various experimental observations very well, leading to profound insight into the hyper-aging dynamics of nanoclay suspensions.     

Langevin dynamics simulations of polymer translocation through nanopores

We investigate the dynamics of polymer translocation through a nanopore using two-dimensional Langevin dynamics simulations. In the absence of an external driving force, we consider a polymer which is initially placed in the middle of the pore and study the escape time tau(e) required for the polymer to completely exit the pore on either side. The distribution of the escape times is wide and has a long tail. We find that tau(e) scales with the chain length N as tau(e) approximately N(1+2nu), where nu is the Flory exponent. For driven translocation, we concentrate on the influence of the friction coefficient xi, the driving force E, and the length of the chain N on the translocation time tau, which is defined as the time duration between the first monomer entering the pore and the last monomer leaving the pore. For strong driving forces, the distribution of translocation times is symmetric and narrow without a long tail and tau approximately E(-1). The influence of xi depends on the ratio between the driving and frictional forces. For intermediate xi, we find a crossover scaling for tau with N from tau approximately N(2nu) for relatively short chains to tau approximately N(1+nu) for longer chains. However, for higher xi, only tau approximately N(1+nu) is observed even for short chains, and there is no crossover behavior. This result can be explained by the fact that increasing xi increases the Rouse relaxation time of the chain, in which case even relatively short chains have no time to relax during translocation. Our results are in good agreement with previous simulations based on the fluctuating bond lattice model of polymers at intermediate friction values, but reveal additional features of dependency on friction.     

A physical model for the longitudinal polarizabilities of polymer chains

The aim of this work is to provide a physical model to relate the polarizability per unit cell of oligomers to that of their corresponding infinite polymer chains. For this we propose an extrapolation method for the polarizability per unit cell of oligomers by fitting them to a physical model describing the dielectric properties of polymer chains. This physical model is based on the concept of a dielectric needle in which we assume a polymer chain to be well described by a cylindrically shaped nonconducting rod with a radius much smaller than its length. With this model we study in which way the polarizability per unit cell approaches the limit of the infinite chain. We show that within this model the macroscopic contribution of the induced electric field to the macroscopic electric field vanishes in the limit of an infinite polymer chain, i.e., there is no macroscopic screening. The macroscopic electric field becomes equal to the external electric field in this limit. We show that this identification leads to a relation between the polarizability per unit cell and the electric susceptibility of the infinite polymer chain. We test our dielectric needle model on the polarizability per unit cell of oligomers of the hydrogen chain and polyacetylene obtained earlier using time-dependent current-density-functional theory in the adiabatic local-density approximation and with the Vignale-Kohn functional. We also perform calculations using the same theory on truly infinite polymer chains by employing periodic boundary conditions. We show that by extrapolating the oligomer results according to our dielectric needle model we get good agreement with our results from calculations on the corresponding infinite polymer chains.