Polymer collapse in the presence of hydrodynamic interactions

Kinetically driven helix formation during the homopolymer collapse process

Using Langevin simulations, we find that simple "generic" bead-and-spring homopolymer chains in a sufficiently bad solvent spontaneously develop helical order during the process of collapsing from an initially stretched conformation. The helix formation is initiated by the unstable modes of the straight chain, which drive the system towards a long-lived metastable transient state. The effect is most pronounced if hydrodynamic interactions are screened.

Analysis of microtubule curvature

The microtubule cytoskeleton in living cells generate and resist mechanical forces to mediate fundamental cell processes, including cell division and migration. Recent advances in digital fluorescence microscopy have enabled the direct observation of bending of individual microtubules in living cells, which has enabled quantitative estimation of the mechanical state of the microtubule array. Although a variety of mechanisms have been proposed, the precise origins of microtubule deformation in living cells remain largely obscure. To investigate these mechanisms and their relative importance in cellular processes, a method is needed to accurately quantify microtubule bending within living cells. Here we describe a method for quantification of bending, using digital fluorescence microscope images to estimate the distribution of curvature in the microtubule. Digital images of individual microtubules can be used to obtain a set of discrete x-y coordinates along the microtubule contour, which is then used to estimate the curvature distribution. Due to system noise and digitization error, the estimate will be inaccurate to some degree. To quantify the inaccuracy, a computational model is used to simulate both the bending of thermally driven microtubules and their observation by digital fluorescence microscopy. This allows for direct comparison between experimental and simulated images, a method which we call model convolution microscopy. We assess the accuracy of various methods and present a suitable method for estimating the curvature distribution for thermally driven semiflexible polymers. Finally, we discuss extensions of the method to quantify microtubule curvature in living cells.

Dynamic correlations in stochastic rotation dynamics

The dynamic structure factor, vorticity and entropy density dynamic correlation functions are measured for stochastic rotation dynamics (SRD), a particle based algorithm for fluctuating fluids. This allows us to obtain unbiased values for the longitudinal transport coefficients such as thermal diffusivity and bulk viscosity. The results are in good agreement with earlier numerical and theoretical results, and it is shown for the first time that the bulk viscosity is indeed zero for this algorithm. In addition, corrections to the self-diffusion coefficient and shear viscosity arising from the breakdown of the molecular chaos approximation at small mean free paths are analyzed. In addition to deriving the form of the leading correlation corrections to these transport coefficients, the probabilities that two and three particles remain collision partners for consecutive time steps are derived analytically in the limit of small mean free path. The results of this paper verify that we have an excellent understanding of the SRD algorithm at the kinetic level and that analytic expressions for the transport coefficients derived elsewhere do indeed provide a very accurate description of the SRD fluid.

Mesoscopic description of solvent effects on polymer dynamics

Solvent effects on polymer dynamics and structure are investigated using a mesoscopic solvent model that accounts for hydrodynamic interactions among the polymer beads. The simulation method combines molecular dynamics of the polymer chain, interacting with the solvent molecules through intermolecular forces, with mesoscopic multiparticle collision dynamics for the solvent molecules. Changes in the intermolecular forces between the polymer beads and mesoscopic solvent molecules are used to vary the solvent conditions from those for good to poor solvents. Polymer collapse and expansion dynamics following changes in solvent conditions are studied for homopolymer and block copolymer solutions. The frictional properties of polymers are also investigated.

Dynamics of polymers in a particle-based mesoscopic solvent

We study the dynamics of flexible polymer chains in solution by combining multiparticle-collision dynamics (MPCD), a mesoscale simulation method, and molecular-dynamics simulations. Polymers with and without excluded-volume interactions are considered. With an appropriate choice of the collision time step for the MPCD solvent, hydrodynamic interactions build up properly. For the center-of-mass diffusion coefficient, scaling with respect to polymer length is found to hold already for rather short chains. The center-of-mass velocity autocorrelation function displays a long-time tail which decays algebraically as (Dt)(-3/2) as a function of time t, where D is the diffusion coefficient. The analysis of the intramolecular dynamics in terms of Rouse modes yields excellent agreement between simulation data and results of the Zimm model for the mode-number dependence of the mode-amplitude correlation functions.

Equilibrium calculation of transport coefficients for a fluid-particle model

A recently introduced particle-based model for fluid flow, called stochastic rotation dynamics, can be made Galilean invariant by introducing a random shift of the computational grid before collisions. In this paper, it is shown how the Green-Kubo relations derived previously can be resummed to obtain exact expressions for the collisional contributions to the transport coefficients. It is also shown that the collisional contribution to the microscopic stress tensor is not symmetric, and that this leads to an additional viscosity. The resulting identification of the transport coefficients for the hydrodynamic modes is discussed in detail, and it is shown that this does not impose restrictions on the applicability of the model. The collisional contribution to the thermal conductivity, which becomes important for small mean free path and small average particle number per cell, is also derived.

Kinetics of the polymer collapse transition: the role of hydrodynamics

We investigate numerically the dynamical behavior of a polymer chain collapsing in a dilute solution. The rate of collapse is measured with and without the presence of hydrodynamic interactions. We find that hydrodynamic interactions accelerate polymer collapse. We present a scaling theory describing the physical process responsible for the collapse kinetics. Predicted collapse times in a hydrodynamic (tauH approximately N(4/3)) and a Brownian heat bath (tauB approximately N2) agree well with the numerical results (tauH approximately N(1.40+/-0.08) and tauB approximately N(1.89+/-0.09)) where N denotes chain length. The folding kinetics of Go models of proteins is also examined. We show that for these systems, where many free energy minima compete, hydrodynamics has little effect on the kinetics.

Two-particle friction in a mesoscopic solvent

The effects of hydrodynamic interactions on the friction tensors for two particles in solution are studied. The particles have linear dimensions on nanometer scales and are either simple spherical particles interacting with the solvent through repulsive Lennard-Jones forces or are composite cluster particles whose atomic components interact with the solvent through repulsive Lennard-Jones forces. The solvent dynamics is modeled at a mesoscopic level through multiparticle collisions that conserve mass, momentum, and energy. The dependence of the two-particle relative friction tensors on the interparticle separation indicates the importance of hydrodynamic interactions for these nanoparticles.

Combining molecular dynamics with Lattice Boltzmann: A hybrid method for the simulation of (charged) colloidal systems

We present a hybrid method for the simulation of colloidal systems that combines molecular dynamics (MD) with the Lattice Boltzmann (LB) scheme. The LB method is used as a model for the solvent in order to take into account the hydrodynamic mass and momentum transport through the solvent. The colloidal particles are propagated via MD and they are coupled to the LB fluid by viscous forces. With respect to the LB fluid, the colloids are represented by uniformly distributed points on a sphere. Each such point [with a velocity V(r) at any off-lattice position r] is interacting with the neighboring eight LB nodes by a frictional force F = xi0(V(r)-u(r)), with xi0 being a friction coefficient and u(r) being the velocity of the fluid at the position r. Thermal fluctuations are introduced in the framework of fluctuating hydrodynamics. This coupling scheme has been proposed recently for polymer systems by Ahlrichs and Dunweg [J. Chem. Phys. 111, 8225 (1999)]. We investigate several properties of a single colloidal particle in a LB fluid, namely, the effective Stokes friction and long-time tails in the autocorrelation functions for the translational and rotational velocity. Moreover, a charged colloidal system is considered consisting of a macroion, counterions, and coions that are coupled to a LB fluid. We study the behavior of the ions in a constant electric field. In particular, an estimate of the effective charge of the macroion is yielded from the number of counterions that move with the macroion in the direction of the electric field.

Kinetic Theory Derivation of the Transport Coefficients of Stochastic Rotation Dynamics

We use a kinetic theory approach to derive the continuum Navier-Stokes and heat conduction equations for stochastic rotation dynamics, a particle based algorithm for simulating a fluid. Hence we obtain expressions for the viscosity and thermal conductivity in two and three dimensions. The predictions are tested numerically and good agreement is found.

Influence of hydrodynamics on many-particle diffusion in 2D colloidal suspensions

We study many-particle diffusion in 2D colloidal suspensions with full hydrodynamic interactions through a novel mesoscopic simulation technique. We focus on the behaviour of the effective scaled tracer and collective-diffusion coefficients DT(rho)/D0 and DC(rho)/D0 respectively, where D0 is the single-particle diffusion coefficient, as a function of the density of the colloids rho. At low Schmidt numbers Sc - 1, we find that hydrodynamics has essentially no effect on the behaviour of DT (rho)/D0. At larger Sc, DT (rho)/D0 seems to be enhanced at all densities, although the differences compared to the case without hydrodynamics are rather minor. The collective-diffusion coefficient, on the other hand, is much more strongly coupled to hydrodynamical conservation laws and is distinctly different from the purely dissipative case without hydrodynamic interactions.