Showing papers in "arXiv: Soft Condensed Matter in 1999"
TL;DR: In this paper, a new efficient method for simulating polymer-solvent systems which combines a lattice Boltzmann approach for the fluid with a continuum molecular dynamics (MD) model for the polymer chain is presented.
Abstract: In this paper we establish a new efficient method for simulating polymer-solvent systems which combines a lattice Boltzmann approach for the fluid with a continuum molecular dynamics (MD) model for the polymer chain. The two parts are coupled by a simple dissipative force while the system is driven by stochastic forces added to both the fluid and the polymer. Extensive tests of the new method for the case of a single polymer chain in a solvent are performed. The dynamic and static scaling properties predicted by analytical theory are validated. In this context, the influence of the finite size of the simulation box is discussed. While usually the finite size corrections scale as L^{-1} (L denoting the linear dimension of the box), the decay rate of the Rouse modes is only subject to an L^{-3} finite size effect. Furthermore, the mapping to an existing MD simulation of the same system is done so that all physical input values for the new method can be derived from pure MD simulation. Both methods can thus be compared quantitatively, showing that the new method allows for much larger time steps. Comparison of the results for both methods indicates systematic deviations due to non-perfect match of the static chain conformations.
286 citations
TL;DR: In this paper, the role of ageing in the rheology of soft materials was studied theoretically and several generalized rheological response functions suited to ageing samples were defined, in which time translation invariance was lost.
Abstract: We study theoretically the role of ageing in the rheology of soft materials. We define several generalized rheological response functions suited to ageing samples (in which time translation invariance is lost). These are then used to study ageing effects within a simple scalar model (the "soft glassy rheology" or SGR model) whose constitutive equations relate shear stress to shear strain among a set of elastic elements, with distributed yield thresholds, undergoing activated dynamics governed by a "noise temperature", $x$. (Between yields, each element follows affinely the applied shear.) For $1
285 citations
TL;DR: In this article, the authors formulate the non-linear field theory for a fluctuating counter-ion distribution in the presence of a fixed, arbitrary charge distribution, and obtain the fluctuation corrections to the electrostatic potential and counterion density to one-loop order without further approximations.
Abstract: We formulate the non-linear field theory for a fluctuating counter-ion distribution in the presence of a fixed, arbitrary charge distribution. The Poisson-Boltzmann equation is obtained as the saddle-point, and the effects of fluctuations and correlations are included by a loop-wise expansion around this saddle point. We show that the Poisson equation is obeyed at each order in the loop expansion and explicitly give the expansion of the Gibbs potential up to two loops. We then apply our formalism to the case of an impenetrable, charged wall, and obtain the fluctuation corrections to the electrostatic potential and counter-ion density to one-loop order without further approximations. The relative importance of fluctuation corrections is controlled by a single parameter, which is proportional to the cube of the counter-ion valency and to the surface charge density. We also calculate effective interactions between charged particles, which reflect counter-ion correlation effects.
246 citations
TL;DR: In this article, the authors investigated the phenomenon of counterion condensation in a solution of highly charged rigid polyelectrolytes within the cell model and proposed a method based on the charge distribution function to identify both the fraction of condensed ions and the radial extension of the condensed layer.
Abstract: We investigate the phenomenon of counterion condensation in a solution of highly charged rigid polyelectrolytes within the cell model. A method is proposed which -- based on the charge distribution function -- identifies both the fraction of condensed ions and the radial extension of the condensed layer. Within salt-free Poisson-Boltzmann (PB) theory it reproduces the well known fraction 1-1/xi of condensed ions for a Manning parameter xi>1. Furthermore, it predicts a weak salt dependence of this fraction and a breakdown of the concept of counterion condensation in the high salt limit. We complement our theoretical investigations with molecular dynamics simulations of a cell-like model, which constantly yield a stronger condensation than predicted by PB theory. While the agreement between theory and simulation is excellent in the monovalent, weakly charged case, it deteriorates with increasing electrostatic interaction strength and, in particular, increasing valence. For instance, at a high concentration of divalent salt and large xi our computer simulations predict charge oscillations, which mean-field theory is unable to reproduce.
217 citations
TL;DR: In this article, an equilibrated model glass-forming liquid is studied by mapping successive configurations produced by molecular dynamics simulation onto a time series of inherent structures (local minima in the potential energy).
Abstract: An equilibrated model glass-forming liquid is studied by mapping successive configurations produced by molecular dynamics simulation onto a time series of inherent structures (local minima in the potential energy). Using this ``inherent dynamics'' approach we find direct numerical evidence for the long held view that below a crossover temperature, $T_x$, the liquid's dynamics can be separated into (i) vibrations around inherent structures and (ii) transitions between inherent structures (M. Goldstein, J. Chem. Phys. {\bf 51}, 3728 (1969)), i.e., the dynamics become ``dominated'' by the potential energy landscape. In agreement with previous proposals, we find that $T_x$ is within the vicinity of the mode-coupling critical temperature $T_c$. We further find that at the lowest temperature simulated (close to $T_x$), transitions between inherent structures involve cooperative, string like rearrangements of groups of particles moving distances substantially smaller than the average interparticle distance.
217 citations
TL;DR: A simple model of globular proteins which incorporates anisotropic attractions is proposed in this article, which is closely related to models used to model simple hydrogen-bonding molecules such as water.
Abstract: A simple model of globular proteins which incorporates anisotropic attractions is proposed. It is closely related to models used to model simple hydrogen-bonding molecules such as water. Theories for both the fluid and solid phases are presented, and phase diagrams calculated. The model protein exhibits a fluid-fluid transition which is metastable with respect to the fluid-solid transition for most values of the model parameters. This is behaviour often observed for globular proteins. The model offers an explanation of the difficulty observed in crystallising some globular proteins and suggests that some proteins may not have a solid phase at all under all but extreme conditions.
166 citations
TL;DR: In this article, the authors present results of molecular dynamics simulations of very long model polymer chains analyzed by various experimentally relevant techniques, showing that the segment motion of the chains is in very good agreement with the repatation model.
Abstract: We present results of molecular dynamics simulations of very long model polymer chains analyzed by various experimentally relevant techniques. The segment motion of the chains is found to be in very good agreement with the repatation model. We also calculated the plateau-modulus G_N. The predicitions of the entanglement length N_e from G_N and from the mean square displacements of the chains segments disagree by a factor of about 2.2(2), indicating an error in the prefactor in the standard formula for G_N. We show that recent neutron spin echo measurements were carried out for chain lengths which are too small for a correct determination of N_e.
161 citations
TL;DR: In this article, a reverse nonequilibrium molecular dynamics method was used to calculate Soret coefficients of various systems in a systematic fashion, and the authors observed that the magnitude of the Soret coefficient depends on all three ratios.
Abstract: The aim of this study is to understand deeper the thermal diffusion transport process (Ludwig-Soret effect) at the microscopic level. For that purpose, the recently developed reverse nonequilibrium molecular dynamics method was used to calculate Soret coefficients of various systems in a systematic fashion. We studied binary Lennard-Jones (LJ) fluids near the triple point (of one of the components) in which we separately changed the ratio of one of the LJ parameters mass, atomic diameter and interaction strength while keeping all other parameters fixed and identical. We observed that the magnitude of the Soret coefficient depends on all three ratios. Concerning its sign we found that heavier species, smaller species and species with higher interaction strengths tend to accumulate in the cold region whereas the other ones (lighter, bigger or weaker bound) migrate to the hot region of our simulation cell. Additionally, the superposition of the influence of the various parameters was investigated as well as more realistic mixtures. We found that in the experimentally relevant parameter range the contributions are nearly additive and that the mass ratio often is the dominating factor.
127 citations
TL;DR: In this article, the distance-resolved effective interaction between two star polymers in a good solvent is calculated by Molecular Dynamics computer simulations, and the results are compared with a pair potential proposed recently by Likos et al. which is exponentially decaying for large distances and crosses over, at the corona diameter of the star, to an ultrasoft logarithmic repulsion for small distances.
Abstract: The distance-resolved effective interaction between two star polymers in a good solvent is calculated by Molecular Dynamics computer simulations. The results are compared with a pair potential proposed recently by Likos et al. [Phys. Rev. Lett. 1998, 80, 4450] which is exponentially decaying for large distances and crosses over, at the corona diameter of the star, to an ultrasoft logarithmic repulsion for small distances. Excellent agreement is found in a broad range of star arm numbers.
107 citations
Posted Content•
TL;DR: The microscopic contact of a cell/semiconductor hybrid with nanostructured semiconductor is investigated with the aim of single channel recording of ion channels in cell membranes to overcome many limitations of the classical patch-clamp technique.
Abstract: We investigate the microscopic contact of a cell/semiconductor hybrid. The semiconductor is nanostructured with the aim of single channel recording of ion channels in cell membranes. This approach will overcome many limitations of the classical patch-clamp technique. The integration of silicon-based devices 'on-chip' promises novel types of experiments on single ion channels.
71 citations
TL;DR: In this paper, a method for computer simulations in the isothermal-isobaric (NPT) ensemble, based on Langevin-type equations of motion for the particle coordinates and the ''piston'' degree of freedom, is derived by straightforward application of the standard Kramers-Moyal formalism.
Abstract: A recently proposed method for computer simulations in the isothermal-isobaric (NPT) ensemble, based on Langevin-type equations of motion for the particle coordinates and the ``piston'' degree of freedom, is re-derived by straightforward application of the standard Kramers-Moyal formalism. An integration scheme is developed which reduces to a time-reversible symplectic integrator in the limit of vanishing friction. This algorithm is hence expected to be quite stable for small friction, allowing for a large time step. We discuss the optimal choice of parameters, and present some numerical test results.
Posted Content•
TL;DR: Excitations of standing and rotating quadrupolar and hexadecapolar modes were excited with high spatial and temporal resolution using the optical dipole force of a rapidly scanning laser beam, which should be useful for the study of rotating Bose-Einstein condensates and vortices.
Abstract: Surface modes in a Bose-Einstein condensate of sodium atoms have been studied. We observed excitations of standing and rotating quadrupolar and octopolar modes. The modes were excited with high spatial and temporal resolution using the optical dipole force of a rapidly scanning laser beam. This novel technique is very flexible and should be useful for the study of rotating Bose-Einstein condensates and vortices.
TL;DR: In this paper, it is shown that ions form strongly correlated liquid at the macroion surface and the cohesiveness of this liquid leads to strong additional attraction of counterions to the surface, which is taken into account by a new boundary condition for the Poisson-Boltzmann equation.
Abstract: Screening of a macroion by multivalent counterions is considered. It is shown that ions form strongly correlated liquid at the macroion surface. Cohesive energy of this liquid leads to strong additional attraction of counterions to the surface. Away from the surface this attraction is taken into account by a new boundary condition for the Poisson-Boltzmann equation. This equation is solved with the new boundary condition for a charged flat surface and a long cylinder. For a cylinder Onsager-Manning theory looses its universality so that apparent charge of the cylinder is smaller than their theory predicts and depends on its bare charge. It can also vanish or change sign.
Posted Content•
TL;DR: In this article, the authors present analytic solutions for steady flow of the Johnson-Segalman (JS) model with a diffusion term in various geometries and under controlled strain rate conditions, using matched asymptotic expansions.
Abstract: We present analytic solutions for steady flow of the Johnson-Segalman (JS) model with a diffusion term in various geometries and under controlled strain rate conditions, using matched asymptotic expansions. The diffusion term represents a singular perturbation that lifts the continuous degeneracy of stable, banded, steady states present in the absence of diffusion. We show that the stable steady flow solutions in Poiseuille and cylindrical Couette geometries always have two bands. For Couette flow and small curvature, two different banded solutions are possible, differing by the spatial sequence of the two bands.
TL;DR: In this paper, the possible phase diagrams for shear-induced phase transitions between two phases are collected, under conditions of both common strain rate and common stress in the two phases, and four fundamental shear stress vs. strain-rate curves and their concentration dependence are discussed.
Abstract: The possible ``phase diagrams'' for shear-induced phase transitions between two phases are collected. We consider shear-thickening and shear-thinning fluids, under conditions of both common strain rate and common stress in the two phases, and present the four fundamental shear stress vs. strain-rate curves and discuss their concentration dependence. We outline how to construct more complicated phase diagrams, discuss in which class various experimental systems fall, and sketch how to reconstruct the phase diagrams from rheological measurements.
Posted Content•
TL;DR: In this paper, the authors investigate the effective potential between two like-charged spheres in divalent counterion solution as a function of concentration and find a strong attractive interaction for high concentration and a global repulsive effective interaction for dilute systems.
Abstract: We simulate model systems of charged spherical particles in their counterion solution and measure the thermodynamic pressure and the pair distribution function from which we derive effective potentials of mean force. For a system with only electrostatic and hard core interactions, we investigate the effective potential between two like-charged spheres in divalent counterion solution as a function of concentration. We find a strong attractive interaction for high concentration and a global repulsive effective interaction for dilute systems. The results indicate a first order phase transition in sphere-counterion density as a function of global concentration and the effective sphere-sphere potentials in the dilute (solvated) regime suggest significant density fluctuations due to short range local minima in the effective energy surface. Our results arise from a minimal approach model of several recent experiments on polystyrene latex particles in monovalent counterion solution.
TL;DR: In this article, the average dipole moment of strongly polarized hydrated water molecules in the first solvation shell of the hexa-hydrate aluminum ion was found to be 5.02 Debye.
Abstract: The solvation of Al and its hydrolyzed species in water clusters has been studied by means of ab initio molecular dynamics simulations. The hexa-hydrate aluminum ion formed a stable complex in the finite temperature cluster simulation of one aluminum ion and 16 waters. The average dipole moment of strongly polarized hydrated water molecules in the first solvation shell of the hexa-hydrate aluminum ion was found to be 5.02 Debye. The deprotonated hexa-hydrate complex evolves into a tetra-coordinated aluminate ion with two water molecules in the second solvation shell forming hydrogen bonds to the hydroxyl groups in agreement with the observed coordination.
Posted Content•
TL;DR: In this article, the authors investigate the equilibrium charge distribution along a single annealed polyelectrolyte chain under different conditions and show that the charge inhomogeneity is localized close to the chain ends over a distance of order of the Debye length.
Abstract: We investigate the equilibrium charge distribution along a single annealed polyelectrolyte chain under different conditions. The coupling between the conformation of the chain and the local charge distribution is described for various solvent qualities and salt concentration. In salt free solution, we find a slight charge depletion in the central part of the chain: the charges accumulate at the ends. The effect is less important if salt is added to the solution since the charge inhomogeneity is localized close to the chain ends over a distance of order of the Debye length. In the case of poor solvent conditions we find a different charging of beads and strings in the framework of the necklace model. This inhomogeneity leads to a charge instability and a first order transition between spherical globules and elongated chains.
TL;DR: A lattice Boltzmann model for amphiphilic fluid dynamics is presented in this paper, which is a ternary model that conserves mass separately for each chemical species present (water, oil, amphiphile).
Abstract: A lattice Boltzmann model for amphiphilic fluid dynamics is presented. It is a ternary model, in that it conserves mass separately for each chemical species present (water, oil, amphiphile), and it maintains an orientational degree of freedom for the amphiphilic species. Moreover, it models fluid interactions at the microscopic level by introducing self-consistent forces between the particles, rather than by positing a Landau free energy functional. This combination of characteristics fills an important need in the hierarchy of models currently available for amphiphilic fluid dynamics, enabling efficient computer simulation and furnishing new theoretical insight. Several computational results obtained from this model are presented and compared to existing lattice-gas model results. In particular, it is noted that lamellar structures, which are precluded by the Peierls instability in two-dimensional systems with kinetic fluctuations, are not observed in lattice-gas models, but are easily found in the corresponding lattice Boltzmann models. This points out a striking difference in the phenomenology accessible to each type of model.
TL;DR: In this paper, the nonlinear stress-strain relation for crosslinked polymer networks is studied using molecular dynamics simulations, and new results for the stress versus strain for both dry and swollen networks are presented.
Abstract: The non-linear stress-strain relation for crosslinked polymer networks is studied using molecular dynamics simulations. Previously we demonstrated the importance of trapped entanglements in determining the elastic and relaxational properties of networks. Here we present new results for the stress versus strain for both dry and swollen networks. Models which limit the fluctuations of the network strands like the tube model are shown to describe the stress for both elongation and compression. For swollen networks, the total modulus is found to decrease like (V_0/V)^{2/3} and goes to the phantom model result only for short strand networks.
TL;DR: In this article, a Debye-Hueckel approach to the one-component plasma is proposed, which employs the exclusion of the charged background in the vicinity of the central ion, thereby yielding a thermodynamically stable free energy density.
Abstract: Based on a Debye-Hueckel approach to the one-component plasma we propose a new free energy for incorporating ionic correlations into Poisson-Boltzmann like theories. Its derivation employs the exclusion of the charged background in the vicinity of the central ion, thereby yielding a thermodynamically stable free energy density, applicable within a local density approximation. This is an improvement over the existing Debye-Hueckel plus hole theory, which in this situation suffers from a "structuring catastrophe". For the simple example of a strongly charged stiff rod surrounded by its counterions we demonstrate that the Poisson-Boltzmann free energy functional augmented by our new correction accounts for the correlations present in this system when compared to molecular dynamics simulations.
TL;DR: In this paper, the authors analyzed the effects of this caging process on the Rouse description of the polymer melt's dynamics and showed that the Rouse theory is applicable for length and time scales above the typical scales for the Caging process.
Abstract: We report results of molecular-dynamics simulations for a glassy polymer melt consisting of short, linear bead-spring chains. It was shown in previous work that this onset of the glassy slowing down is compatible with the predictions of the mode coupling theory. The physical process of `caging' of a monomer by its spatial neighbors leads to a distinct two step behavior in the particle mean square displacements. In this work we analyze the effects of this caging process on the Rouse description of the melt's dynamics. We show that the Rouse theory is applicable for length and time scales above the typical scales for the caging process. Futhermore, the monomer displacement is compared with simulation results for a binary Lennard-Jones mixture to point out the differences which are introduced by the connectivity of the particles.
Posted Content•
TL;DR: With the help of a simple 20‐letter lattice model of heteropolymers, the energy landscape in the space of designed good‐folder sequences is investigated, finding residues that play a key role in the foldability of the chain and in the stability of the native state are highly conserved.
Abstract: With the help of a simple 20 letters, lattice model of heteropolymers, we investigate the energy landscape in the space of designed good-folder sequences. Low-energy sequences form clusters, interconnected via neutral networks, in the space of sequences. Residues which play a key role in the foldability of the chain and in the stability of the native state are highly conserved, even among the chains belonging to different clusters. If, according to the interaction matrix, some strong attractive interactions are almost degenerate (i.e. they can be realized by more than one type of aminoacid contacts) sequence clusters group into a few super-clusters. Sequences belonging to different super-clusters are dissimilar, displaying very small ($\approx 10%$) similarity, and residues in key-sites are, as a rule, not conserved. Similar behavior is observed in the analysis of real protein sequences.
Posted Content•
TL;DR: In this article, a structural analysis of a stiff polyelectrolyte chain with intrachain attractions is presented, and the authors show that the mean-field solution corresponds to an Euler-like buckling instability.
Abstract: We present a systematic statistical mechanical analysis of the conformational properties of a stiff polyelectrolyte chain with intrachain attractions that are due to counterion correlations. We show that the mean-field solution corresponds to an Euler-like buckling instability. The effect of the conformational fluctuations on the buckling instability is investigated, first, qualitatively, within the harmonic (``semiclassical'') theory, then, systematically, within a 1/d-expansion, where d denotes the dimension of embedding space. Within the ``semiclassical'' approximation, we predict that the effect of fluctuations is to renormalize the effective persistence length to smaller values, but not to change the nature of the mean-field (i.e., buckling) behavior. Based on the 1/d-expansion we are, however, led to conclude that thermal fluctuations are responsible for a change of the buckling behavior which is turned into a polymer collapse. A phase diagram is constructed in which a sequence of collapse transitions terminates at a buckling instability that occurs at a place that varies with the magnitude of the bare persistence length of the polymer chain, as well as with the strength and range of the attractive potential.
Posted Content•
TL;DR: In this article, a field theory for particles which interact via Coulomb and hard core potentials is presented for one-component plasma (OCP) with hard cores, consisting of identical particles of fixed charge and diameter in a neutralizing background.
Abstract: A field theory is presented for particles which interact via Coulomb and hard-core potentials. We apply the method to the one-component plasma (OCP) with hard cores, consisting of identical particles of fixed charge and diameter in a neutralizing background, and the symmetric two-component plasma (TCP) with hard cores, consisting of equal numbers of positively and negatively charged particles of identical size. We obtain exactly the first few coefficients of a systematic low-density expansion of the free energy for both models. The OCP coefficients go over to the classical Abe result in the high-charge (or small hard-core diameter) limit. The TCP, on the other hand, exhibits diverging coefficients in the high-charge limit, which is due to the formation of strongly bound ion-pairs.
TL;DR: In this paper, the mean square displacement of a center monomer of a self-avoiding polymer embedded in regularly distributed obstacles with spacing a=20 and confined in two dimensions is studied numerically using the extended bond fluctuation method.
Abstract: Diffusion properties of a self-avoiding polymer embedded in regularly distributed obstacles with spacing a=20 and confined in two dimensions is studied numerically using the extended bond fluctuation method which we have developed recently. We have observed for the first time to our knowledge, that the mean square displacement of a center monomer $\phi_{M/2}(t)$ exhibits four dynamical regimes, i.e., $\phi_{M/2}(t) \sim t^{
u_m}$ with $
u_m\sim 0.6$, 3/8, 3/4, and 1 from the shortest to longest time regimes. The exponents in the second and third regimes are well described by segmental diffusion in the ``self-avoiding tube''. In the fourth (free diffusion) regime, we have numerically confirmed the relation between the reptation time $\tau_d$ and the number of segments $M, \tau_d\propto M^3$.
Posted Content•
TL;DR: In this article, the authors present new simulation results on a hard sphere system at high densities using three-time correlations, which can account for the anomalous diffusion, which results from a homogeneous back-dragging effect.
Abstract: We present new simulation results on a hard sphere system at high densities. Using three-time correlations, we can account for the anomalous diffusion, which results from a homogeneous back-dragging effect. Furthermore, we calculate the non-gaussian parameter and connect it to the existence of dynamic heterogeneities.
TL;DR: In this paper, the authors consider a porous solid covered with a water film and show that the dynamics of dewetting is strongly modified when the liquid is pumped in, either spontaneously (if the porous medium is hydrophilic) or by an external pump.
Abstract: We consider a porous solid covered with a water film (or with a drop) in situations where the liquid is pumped in, either spontaneously (if the porous medium is hydrophilic) or mechanically (by an external pump). The dynamics of dewetting is then strongly modified. We analyse a few major examples: a) horizontal films, which break at a certain critical thickness, b) the "modified Landau-Levich problem" where a porous plate moves up from a bath and carries a film: aspiration towards the plate limits the height H reached by the film, c) certain situation where the hysteresis of contact angles is important.
TL;DR: In this article, self-assembly techniques, x-ray reflectivity measurements, and computer simulations are used to measure the effective adhesion force acting on the rods in a model consisting of a planar array of uniformly positively charged, stiff rods and a negatively charged planar substrate.
Abstract: Using self-assembly techniques, x-ray reflectivity measurements, and computer simulations, we study the effective interaction between charged polymer rods and surfaces. Long-time Brownian dynamics simulations are used to measure the effective adhesion force acting on the rods in a model consisting of a planar array of uniformly positively charged, stiff rods and a negatively charged planar substrate in the presence of explicit monovalent counterions and added monovalent salt ions in a continuous, isotropic dielectric medium. This electrostatic model predicts an attractive polymer-surface adhesion force that is weakly dependent on the bulk salt concentration and that shows fair agreement with a Debye-Huckel approximation for the macroion interaction at salt concentrations near 0.1 M. Complementary x-ray reflectivity experiments on poly(diallyldimethyl ammonium) chloride (PDDA) monolayer films on the native oxide of silicon show that monolayer structure, electron density, and surface roughness are likewise independent of the bulk ionic strength of the solution.
TL;DR: In this article, the electron density maps of the ripple phase of the DMPC and POPC multibilayers at different temperatures and fixed relative humidity were calculated, and the existence of an average tilt of the hydrocarbon chains of the lipid molecules along the direction of ripple wave vector was established, which is responsible for the occurrence of asymmetric ripples in these systems.
Abstract: We have calculated the electron density maps of the ripple phase of dimyristoylphosphatidylcholine (DMPC) and palmitoyl-oleoyl phosphatidylcholine (POPC) multibilayers at different temperatures and fixed relative humidity. Our analysis establishes, for the first time, the existence of an average tilt of the hydrocarbon chains of the lipid molecules along the direction of the ripple wave vector, which we believe is responsible for the occurrence of asymmetric ripples in these systems.