End-monomer Dynamics in Semiflexible Polymers

Spurred by an experimental controversy in the literature, we investigate the end-monomer dynamics of semiflexible polymers through Brownian hydrodynamic simulations and dynamic mean-field theory. Precise experimental observations over the last few years of end-monomer dynamics in the diffusion of double-stranded DNA have given conflicting results: one study indicated an unexpected Rouse-like scaling of the mean squared displacement (MSD) ?r(2)(t)? ~ t(1/2) at intermediate times, corresponding to fluctuations at length scales larger than the persistence length but smaller than the coil size; another study claimed the more conventional Zimm scaling ?r(2)(t)? ~ t(2/3) in the same time range. Using hydrodynamic simulations, analytical and scaling theories, we find a novel intermediate dynamical regime where the effective local exponent of the end-monomer MSD, α(t) = d log?r(2)(t)?/d log t, drops below the Zimm value of 2/3 for sufficiently long chains. The deviation from the Zimm prediction increases with chain length, though it does not reach the Rouse limit of 1/2. The qualitative features of this intermediate regime, found in simulations and in an improved mean-field theory for semiflexible polymers, in particular the variation of α(t) with chain and persistence lengths, can be reproduced through a heuristic scaling argument. Anomalously low values of the effective exponent α are explained by hydrodynamic effects related to the slow crossover from dynamics on length scales smaller than the persistence length to dynamics on larger length scales.     

Pair diffusion, hydrodynamic interactions, and available volume in dense fluids

We calculate the pair diffusion coefficient D(r) as a function of the distance r between two hard sphere particles in a dense monodisperse fluid. The distance-dependent pair diffusion coefficient describes the hydrodynamic interactions between particles in a fluid that are central to theories of polymer and colloid dynamics. We determine D(r) from the propagators (Green's functions) of particle pairs obtained from molecular dynamics simulations. At distances exceeding ~3 molecular diameters, the calculated pair diffusion coefficients are in excellent agreement with predictions from exact macroscopic hydrodynamic theory for large Brownian particles suspended in a solvent bath, as well as the Oseen approximation. However, the asymptotic 1/r distance dependence of D(r) associated with hydrodynamic effects emerges only after the pair distance dynamics has been followed for relatively long times, indicating non-negligible memory effects in the pair diffusion at short times. Deviations of the calculated D(r) from the hydrodynamic models at short distances r reflect the underlying many-body fluid structure, and are found to be correlated to differences in the local available volume. The procedure used here to determine the pair diffusion coefficients can also be used for single-particle diffusion in confinement with spherical symmetry.     

Critical polymer-polymer phase separation in ternary solutions

We study polymer-polymer phase separation in a common good solvent by means of Monte Carlo simulations of the bond-fluctuation model. Below a critical, chain-length-dependent concentration, no phase separation occurs. For higher concentrations, the critical demixing temperature scales nonlinearly with the total monomer concentration, with a power law relatively close to a renormalization-group prediction based on "blob" scaling arguments. We point out that earlier simulations and experiments have tested this power-law dependence at concentrations outside the validity regime of the scaling arguments. The critical amplitudes of the order parameter and the zero-angle scattering intensity also exhibit chain-length dependences that differ from the conventional predictions but are in excellent agreement with the renormalization-group results. In addition, we characterize the variation of the average coil shape upon phase separation.     

Structure, dynamics, and phase transitions of tethered membranes: a Monte Carlo simulation study

A coarse-grained model of a self-avoiding tethered membrane with hexagonal coordination, embedded in three-dimensional space, is studied by means of extensive Monte Carlo computer simulations. The simulations are performed at various temperatures for membranes with linear size 5< or =L< or =50. We find that the membrane undergoes several folding transitions from a high-temperature flat phase to multiple-folded structure as the temperature is steadily decreased. Using a suitable order parameter and finite size scaling analysis, these phase transitions are shown to be of first order. The equilibrium shape of the membranes is analyzed by calculating the eigenvalues lambda(max) (2)> or =lambda(med) (2)> or =lambda(min) (2) of the inertia tensor. We present a systematic finite size scaling analysis of the radius of gyration and the eigenvalues of the inertia tensor at different phases of the observed folding transitions. In the high-temperature flat phase, the radius of gyration R(g) grows with the linear size of the membrane L as R(g) proportional to L(nu), where the exponent nu is approximately equal to 1.0. The eigenvalues of the inertia tensor scale as lambda(max) proportional to lambda(med) proportional to L(nu) and lambda(min) proportional to L(nu(min) ), whereby the roughness exponent nu(min) is approximately equal to 0.7. We also find that the time tau(R) of a self-avoiding membrane to diffuse a distance R(g) scales as tau(R) proportional to L(2nu+2), which is in good agreement with the theoretical predictions.     

Scaling behavior of nonisothermal phase separation

The phase separation process in a critical mixture of polydimethylsiloxane and polyethylmethylsiloxane (PDMS/PEMS, a system with an upper critical solution temperature) was investigated by time-resolved light scattering during continuous quenches from the one-phase into the two-phase region. Continuous quenches were realized by cooling ramps with different cooling rates kappa. Phase separation kinetics is studied by means of the temporal evolution of the scattering vector qm and the intensity Im at the scattering peak. The curves qm(t) for different cooling rates can be shifted onto a single mastercurve. The curves Im(t) show similar behavior. As shift factors, a characteristic length Lc and a characteristic time tc are introduced. Both characteristic quantities depend on the cooling rate through power laws: Lc approximately kappa(-delta) and tc approximately kappa(-rho). Scaling behavior in isothermal critical demixing is well known. There the temporal evolutions of qm and Im for different quench depths DeltaT can be scaled with the correlation length xi and the interdiffusion coefficient D, both depending on DeltaT through critical power laws. We show in this paper that the cooling rate scaling in nonisothermal demixing is a consequence of the quench depth scaling in the isothermal case. The exponents delta and rho are related to the critical exponents nu and nu* of xi and D, respectively. The structure growth during nonisothermal demixing can be described with a semiempirical model based on the hydrodynamic coarsening mechanism well known in the isothermal case. In very late stages of nonisothermal phase separation a secondary scattering maximum appears. This is due to secondary demixing. We explain the onset of secondary demixing by a competition between interdiffusion and coarsening.     

Lattice-Boltzmann simulation of coalescence-driven island coarsening

A two-dimensional lattice-Boltzmann model (LBM) with fluid-fluid interactions was used to simulate first-order phase separation in a thin fluid film. The intermediate asymptotic time dependence of the mean island size, island number concentration, and polydispersity were determined and compared with the predictions of the distribution-kinetics model. The comparison revealed that the combined effects of growth, coalescence, and Ostwald ripening control the phase transition process in the LBM simulations. However, the overall process is dominated by coalescence, which is independent of island mass. As the phase transition advances, the mean island size increases, the number of islands decrease, and the polydispersity approaches unity, which conforms to the predictions of the distribution-kinetics model. The effects of the domain size on the intermediate asymptotic island size distribution, scaling form of the island size distribution, and the crossover to the long-term asymptotic behavior were elucidated.     

Phase separation of the Widom-Rowlinson mixture in confined geometry

The behavior of binary Widom-Rowlinson [J. Chem. Phys. 52, 1670 (1970)] mixture in confined geometry is investigated. The influence of confinement on phase separation is examined. Coexistence curves for the mixture in slitlike pores and pores of complex geometry are calculated in Monte Carlo simulations. Finite size scaling analysis is used to determine precisely the location of critical point for the bulk mixture.     

Gas-liquid phase separation in oppositely charged colloids: stability and interfacial tension

We study the phase behavior and the interfacial tension of the screened Coulomb (Yukawa) restricted primitive model (YRPM) of oppositely charged hard spheres with diameter sigma using Monte Carlo simulations. We determine the gas-liquid and gas-solid phase transitions using free energy calculations and grand-canonical Monte Carlo simulations for varying inverse Debye screening length kappa. We find that the gas-liquid phase separation is stable for kappasigma相似文献   

Temperature-induced gelation in dilute nanofluids

We report a temperature-induced gelation in dilute nanofluids containing surfactant capped iron oxide and alumina particles of average diameter ~10 nm. We observe a dramatic enhancement in the elastic modulus, viscous modulus, and viscosity, by 3-6 orders of magnitude for a volume fraction (φ) less than 0.035, above a critical shear rate ( ?γ(c)) and temperature (T(c)). The T(c) follows a weak power law scaling with φ as T(c) ~ φ(β), where the scaling exponent β is found to be -0.24. The observed gel-like transition at elevated temperature is attributed to strong van der Waals attractions on the kT energy scale due to poor solvent conditions, which is reminiscent of the phase behavior reported in polymer-coated colloids.     

Depletion induced isotropic-isotropic phase separation in suspensions of rod-like colloids

When non-adsorbing polymers are added to an isotropic suspension of rod-like colloids, the colloids effectively attract each other via depletion forces. We performed Monte Carlo simulations to study the phase diagram of such rod-polymer mixture. The colloidal rods were modeled as hard spherocylinders; the polymers were described as spheres of the same diameter as the rods. The polymers may overlap with no energy cost, while the overlap of polymers and rods is forbidden. Large amounts of depletant cause phase separation of the mixture. We estimated the phase boundaries of isotropic-isotropic coexistence both in the bulk and in confinement. To determine the phase boundaries we applied the grand canonical ensemble using successive umbrella sampling [J. Chem. Phys. 120, 10925 (2004)], and we performed a finite size scaling analysis to estimate the location of the critical point. The results are compared with predictions of the free volume theory developed by Lekkerkerker and Stroobants [Nuovo Cimento D 16, 949 (1994)]. We also give estimates for the interfacial tension between the coexisting isotropic phases and analyze its power-law behavior on the approach of the critical point.     

Hydrodynamic effects on spinodal decomposition kinetics in planar lipid bilayer membranes

The formation and dynamics of spatially extended compositional domains in multicomponent lipid membranes lie at the heart of many important biological and biophysical phenomena. While the thermodynamic basis for domain formation has been explored extensively in the past, domain growth in the presence of hydrodynamic interactions both within the (effectively) two-dimensional membrane and in the three-dimensional solvent in which the membrane is immersed has received little attention. In this work, we explore the role of hydrodynamic effects on spinodal decomposition kinetics via continuum simulations of a convective Cahn-Hilliard equation for membrane composition coupled to the Stokes equation. Our approach explicitly includes hydrodynamics both within the planar membrane and in the three-dimensional solvent in the viscously dominated flow regime. Numerical simulations reveal that dynamical scaling breaks down for critical lipid mixtures due to distinct coarsening mechanisms for elongated versus more isotropic compositional lipid domains. The breakdown in scaling should be readily observable in experiments on model membrane systems.     

Kinetics of phase separation and coarsening in dilute surfactant pentaethylene glycol monododecyl ether solutions

We investigated the phase separation phenomena in dilute surfactant pentaethylene glycol monodedecyl ether (C(12)E(5)) solutions focusing on the growth law of separated domains. The solutions confined between two glass plates were found to exhibit the phase inversion, characteristic of the viscoelastic phase separation; the majority phase (water-rich phase) nucleated as droplets and the minority phase (micelle-rich phase) formed a network temporarily, then they collapsed into an usual sea-island pattern where minority phase formed islands. We found from the real-space microscopic imaging that the dynamic scaling hypothesis did not hold throughout the coarsening process. The power law growth of the domains with the exponent close to 1/3 was observed even though the coarsening was induced mainly by hydrodynamic flow, which was explained by Darcy's law of laminar flow.     

The critical adsorption point of self-avoiding walks: a finite-size scaling approach

The critical adsorption of self-avoiding polymer chain in a simple cubic lattice onto a flat surface is studied with Monte Carlo simulations. The dependence of number of surface contacts M on chain length N and polymer-surface interaction epsilon is investigated by a finite-size scaling approach. We estimate the critical adsorption point epsilon(c)=0.291+/-0.002 and the exponent phi=0.54+/-0.01. The asymptotic behaviors M proportional variant N for epsilon>epsilon(c) and M proportional variant N(0) for epsilon相似文献   

Analytical phase diagrams for colloids and non-adsorbing polymer

We review the free-volume theory (FVT) of Lekkerkerker et al. [Europhys. Lett. 20 (1992) 559] for the phase behavior of colloids in the presence of non-adsorbing polymer and we extend this theory in several aspects: (i) We take the solvent into account as a separate component and show that the natural thermodynamic parameter for the polymer properties is the insertion work Pi(v), where Pi is the osmotic pressure of the (external) polymer solution and v the volume of a colloid particle. (ii) Curvature effects are included along the lines of Aarts et al. [J. Phys.: Condens. Matt. 14 (2002) 7551] but we find accurate simple power laws which simplify the mathematical procedure considerably. (iii) We find analytical forms for the first, second, and third derivatives of the grand potential, needed for the calculation of the colloid chemical potential, the pressure, gas-liquid critical points and the critical endpoint (cep), where the (stable) critical line ends and then coincides with the triple point. This cep determines the boundary condition for a stable liquid. We first apply these modifications to the so-called colloid limit, where the size ratio q(R)=R/a between the radius of gyration R of the polymer and the particle radius a is small. In this limit the binodal polymer concentrations are below overlap: the depletion thickness delta is nearly equal to R, and Pi can be approximated by the ideal (van't Hoff) law Pi=Pi(0)=phi/N, where phi is the polymer volume fraction and N the number of segments per chain. The results are close to those of the original Lekkerkerker theory. However, our analysis enables very simple analytical expressions for the polymer and colloid concentrations in the critical and triple points and along the binodals as a function of q(R). Also the position of the cep is found analytically. In order to make the model applicable to higher size ratio's q(R) (including the so-called protein limit where q(R)>1) further extensions are needed. We introduce the size ratio q=delta/a, where the depletion thickness delta is no longer of order R. In the protein limit the binodal concentrations are above overlap. In such semidilute solutions delta approximately xi, where the De Gennes blob size (correlation length) xi scales as xi approximately phi(-gamma), with gamma=0.77 for good solvents and gamma=1 for a theta solvent. In this limit Pi=Pi(sd) approximately phi(3gamma). We now apply the following additional modifications: With these latter two modifications we obtain again a fully analytical model with simple equations for critical and triple points as a function of q(R). In the protein limit the binodal polymer concentrations scale as q(R)(1/gamma), and phase diagrams phiq(R)(-1/gamma) versus the colloid concentration eta become universal (i.e., independent of the size ratio q(R)). The predictions of this generalized free-volume theory (GFVT) are in excellent agreement with experiment and with computer simulations, not only for the colloid limit but also for the protein limit (and the crossover between these limits). The q(R)(1/gamma) scaling is accurately reproduced by both simulations and other theoretical models. The liquid window is the region between phi(c) (critical point) and phi(t) (triple point). In terms of the ratio phi(t)/phi(c) the liquid window extends from 1 in the cep (here phi(t)-phi(c)=0) to 2.2 in the protein limit. Hence, the liquid window is narrow: it covers at most a factor 2.2 in (external) polymer concentration.     

Surface tension and vapor-liquid phase coexistence of confined square-well fluid

Phase equilibria of a square-well fluid in planar slit pores with varying slit width are investigated by applying the grand-canonical transition-matrix Monte Carlo (GC-TMMC) with the histogram-reweighting method. The wall-fluid interaction strength was varied from repulsive to attractive such that it is greater than the fluid-fluid interaction strength. The nature of the phase coexistence envelope is in agreement with that given in literature. The surface tension of the vapor-liquid interface is calculated via molecular dynamics simulations. GC-TMMC with finite size scaling is also used to calculate the surface tension. The results from molecular dynamics and GC-TMMC methods are in very good mutual agreement. The vapor-liquid surface tension, under confinement, was found to be lower than the bulk surface tension. However, with the increase of the slit width the surface tension increases. For the case of a square-well fluid in an attractive planar slit pore, the vapor-liquid surface tension exhibits a maximum with respect to wall-fluid interaction energy. We also report estimates of critical properties of confined fluids via the rectilinear diameter approach.     

Computer simulation of linkage of two ring chains

We performed off-lattice Monte Carlo simulations of links of two model ring chains with chain length N up to 32,768 in the theta solution or amorphous bulk state by using a random walk model (Model I), and molecular dynamics simulations of two model ring chains in solution with excluded volume interaction (Model II) to investigate topological effects on the geometry of link and ring conformation. In the case of Model I, the mean squared linking number, its distribution, and the size of two chains with fixed linking number are investigated. Our simulation results confirm the previous theoretical prediction that the mean squared linking number decays as pe(-qs(2)) with the distance of centers of chain mass s, where p and q are found to be chain length dependent and q asymptotically approaches to 0.75 as chain length increases. The linking number distribution of two chains has a universal form for long chains, but our simulation results clearly show that the distribution function deviates from the Gaussian distribution, a fact not predicted by any previous theoretical work. A scaling prediction is proposed to predict the link size, and is checked for our simulations for the Model II. The simulation results confirmed the scaling prediction of the blob picture that the link with linking number m occupies a compact volume of m blobs, and the size of the link is asymptotic to R(L) ≈ bN(ν)m(1/3-ν), where N is the chain length, and v is the Flory exponent of polymer in solutions.     

Behavior of a polymer chain immersed in a binary mixture of solvents

The behavior of a polymer chain immersed in a binary solvent mixture is investigated via a single-polymer simulation using an effective Hamiltonian, where the solvent effects are taken into account through a density-functional theory for polymer-solvent admixtures. The liquid-liquid phase separation of the binary solvent mixture is modeled as that of a Lennard-Jones binary fluid mixture with weakly attractive interactions between the different components. Two types of energetic preferences of the polymer chain for the better solvent-(A) no preferential solvophilicity and (B) strong preferential solvophilicity-are employed as polymer-solvent interaction models. The radius of gyration and the polymer-solvent radial distribution functions are determined from the simulations of various molar fractions along an isotherm slightly above the critical temperature of the liquid-liquid phase separation. These quantities near the critical point conspicuously depend on the strength of the preferential solvophilicity. In the case where the polymer exhibits a strong preferential solvophilicity, a remarkable expansion of the polymer chain is observed near the critical point. On the other hand, in the case where the polymer has no preferential solvophilicity, no characteristic variation of the polymer conformation is observed even near the critical point. These results indicate that the expansion of a polymer chain enhances the local phase separation around it, acting as a nucleus of demixing in the vicinity of the critical point. This phenomenon in binary solvents near the liquid-liquid critical point is similar to the expansion of the polymer chain in one-component supercritical solvents near the liquid-vapor critical point, which we have reported [T. Sumi and H. Sekino J. Chem. Phys. 122, 194910 (2005)].     

Capillary condensation in cylindrical pores: Monte Carlo study of the interplay of surface and finite size effects

When a fluid that undergoes a vapor to liquid transition in the bulk is confined to a long cylindrical pore, the phase transition is shifted (mostly due to surface effects at the walls of the pore) and rounded (due to finite size effects). The nature of the phase coexistence at the transition depends on the length of the pore: for very long pores, the system is axially homogeneous at low temperatures. At the chemical potential where the transition takes place, fluctuations occur between vapor- and liquidlike states of the cylinder as a whole. At somewhat higher temperatures (but still far below bulk criticality), the system at phase coexistence is in an axially inhomogeneous multidomain state, where long cylindrical liquid- and vaporlike domains alternate. Using Monte Carlo simulations for the Ising/lattice gas model and the Asakura-Oosawa model of colloid-polymer mixtures, the transition between these two different scenarios is characterized. It is shown that the density distribution changes gradually from a double-peak structure to a triple-peak shape, and the correlation length in the axial direction (measuring the equilibrium domain length) becomes much smaller than the cylinder length. The (rounded) transition to the disordered phase of the fluid occurs when the axial correlation length has decreased to a value comparable to the cylinder diameter. It is also suggested that adsorption hysteresis vanishes when the transition from the simple domain state to the multidomain state of the cylindrical pore occurs. We predict that the difference between the pore critical temperature and the hysteresis critical temperature should increase logarithmically with the length of the pore.     

Entropic wetting and the free isotropic-nematic interface of hard colloidal platelets

We study bulk and interfacial properties of a model suspension of hard colloidal platelets with continuous orientations and vanishing thickness using both density functional theory, based on either a second virial approach or fundamental measure theory (FMT), and Monte Carlo (MC) simulations. We calculate the bulk equation of state, bulk isotropic-nematic (IN) coexistence, and properties of the (planar) free IN interface and of adsorption at a planar hard wall, where we find complete wetting of the nematic phase at the isotropic-wall interface upon approaching bulk IN coexistence. We investigate in detail the asymptotic decay of correlations at large distances. In all cases, the results from FMT and MC agree quantitatively. Our findings are of direct relevance to understanding interfacial properties of dispersions of colloidal platelets.     

Structure, phase behavior, and inhomogeneous fluid properties of binary dendrimer mixtures

The effective pair potentials between different kinds of dendrimers in solution can be well approximated by appropriate Gaussian functions. We find that in binary dendrimer mixtures the range and strength of the effective interactions depend strongly upon the specific dendrimer architecture. We consider two different types of dendrimer mixtures, employing the Gaussian effective pair potentials, to determine the bulk fluid structure and phase behavior. Using a simple mean field density functional theory (DFT) we find good agreement between theory and simulation results for the bulk fluid structure. Depending on the mixture, we find bulk fluid-fluid phase separation (macrophase separation) or microphase separation, i.e., a transition to a state characterized by undamped periodic concentration fluctuations. We also determine the inhomogeneous fluid structure for confinement in spherical cavities. Again, we find good agreement between the DFT and simulation results. For the dendrimer mixture exhibiting microphase separation, we observe a rather striking pattern formation under confinement.