Lamellar phase coexistence induced by electrostatic interactions

Membranes containing highly charged biomolecules can have a minimal free-energy state at small separations that originates in the strongly correlated electrostatic interactions mediated by counterions. This phenomenon can lead to a condensed, lamellar phase of charged membranes that coexists in thermodynamic equilibrium with a very dilute membrane phase. Although the dilute phase is mostly water, entropy dictates that this phase must contain some membranes and counterions. Thus, electrostatics alone can give rise to the coexistence of a condensed and an unbound lamellar phase. We use numerical simulations to predict the nature of this coexistence when the charge density of the membrane is large, for the case of multivalent counterions and for a membrane charge that is characteristic of biomolecules. We also investigate the effects of counterion size and salt on the two coexisting phases. With increasing salt concentration, we predict that electrostatic screening by salt can destroy the phase separation.     

Electrophoresis of a charge-inverted macroion complex: Molecular-dynamics study

We have performed molecular-dynamics simulations to study the effect of an external electric field on a macroion in the solution of multivalent Z : 1 salt. To obtain plausible hydrodynamics of the medium, we explicitly make the simulation of many neutral particles along with ions. In a weak electric field, the macroion drifts together with the strongly adsorbed multivalent counterions along the electric field, in the direction proving inversion of the charge sign. The reversed mobility of the macroion is insensitive to the external field, and increases with salt ionic strength. The reversed mobility takes a maximal value at intermediate counterion valence. The motion of the macroion complex does not induce any flow of the neutral solvent away from the macroion, which reveals screening of hydrodynamic interactions at short distances in electrolyte solutions. A very large electric field, comparable to the macroion unscreened field, disrupts charge inversion by stripping the adsorbed counterions off the macroion. Received 5 December 2001 and Received in final form 10 April 2002     

Macroions in salty water with multivalent ions: giant inversion of charge

Screening of a strongly charged macroion by oppositely charged colloidal particles, micelles, or short polyelectrolytes is considered. Because of strong lateral repulsion such multivalent counterions form a strongly correlated liquid at the surface of the macroion. This liquid provides correlation-induced attraction of multivalent counterions to the macroion surface. As a result even a moderate concentration of multivalent counterions in the solution inverts the sign of the net macroion charge. We show that at high concentration of monovalent salt the absolute value of inverted charge can be larger than the bare one. This giant inversion of charge can be observed in electrophoresis.     

Simulations of stretching a flexible polyelectrolyte with varying charge separation

We calculated the force-extension curves for a flexible polyelectrolyte chain with varying charge separations by performing Monte Carlo simulations of a 5000 bead chain using a screened Coulomb interaction. At all charge separations, the force-extension curves exhibit a Pincus-like scaling regime at intermediate forces and a logarithmic regime at large forces. As the charge separation increases, the Pincus regime shifts to a larger range of forces and the logarithmic regime starts are larger forces. We also found that force-extension curve for the corresponding neutral chain has a logarithmic regime. Decreasing the diameter of bead in the neutral chain simulations removed the logarithmic regime, and the force-extension curve tends to the freely jointed chain limit. This result shows that only excluded volume is required for the high force logarithmic regime to occur.     

Measurement of long-range repulsive forces between charged particles at an oil-water interface

Using a laser tweezers method, we have determined the long-range repulsive force as a function of separation between two charged, spherical polystyrene particles (2.7 microm diameter) present at a nonpolar oil-water interface. At large separations (6 to 12 microm between particle centers) the force is found to decay with distance to the power -4 and is insensitive to the ionic strength of the aqueous phase. The results are consistent with a model in which the repulsion arises primarily from the presence of a very small residual electric charge at the particle-oil interface. This charge corresponds to a fractional dissociation of the total ionizable (sulfate) groups present at the particle-oil surface of approximately 3 x 10(-4).     

Variational theory for a single polyelectrolyte chain

Variational methods are applied to a single polyelectrolyte chain. The polymer is modeled as a Gaussian chain with screened electrostatic repulsion between all monomers. As a variational Hamiltonian, the most general Gaussian kernel, including the possibility of a classical or mean polymer path, is employed. The resulting self-consistent equations are systematically solved both for large and small monomer-monomer separations along the chain. In the absence of screening, the polymer is stretched on average. It is described by a straight classical path with Gaussian fluctuations around it. If the electrostatic repulsion is screened, the polymer is isotropically swollen for large separations, and for small separations the polymer correlation function is calculated as an analytic expansion in terms of the monomer-monomer separation along the chain. The electrostatic persistence length and the electrostatic blobsize are inferred from the crossover between distinct scaling ranges. We perform a global analysis of the scaling behavior as a function of the screening length and electrostatic interaction strength , where is the Bjerrum length and A is the distance of charges along the polymer chain. We find three different scaling regimes. i) A Gaussian-persistent regime with Gaussian behavior at small, persistent behavior at intermediate, and isotropically swollen behavior at large length scales. This regime occurs for weakly charged polymers and only for intermediate values of the screening length. The electrostatic persistence length is defined as the crossover length between the persistent and the asymptotically swollen behavior and is given by and thus disagrees with previous (restricted) variational treatments which predict a linear dependence on the screening length .ii) A Gaussian regime with Gaussian behavior at small and isotropically swollen behavior at large length scales. This regime occurs for weakly charged polymers and/or strong screening, and the electrostatic repulsion between monomers only leads to subfluent corrections to Gaussian scaling at small separations. The concept of a persistence length is without meaning in this regime. iii) A persistent regime , where the chain resembles a stretched rod on intermediate and small scales. Here the persistence length is given by the original Odijk prediction, , if the overstretching of the chain is avoided. We also investigate the effects of a finite polymer length and of an additional excluded-volume interaction, which modify the resultant scaling behavior. Applications to experiments and computer simulations are discussed. Received 24 December 1997     

Debye-Hückel theory for slab geometries

The electrostatic interaction of charged particles through or at a low-dielectric slab, such as a lipid bilayer immersed in water or a self-assembled monolayer (SAM) on a metal substrate, is considered theoretically in the presence of salt within the Gaussian approximation using a generalized Green's formalism. A number of separate situations are discussed: i) The presence of a low-dielectric slab leads to pronounced interactions of a single charge with the slab via the formation of polarization surface charges. For SAMs on metal substrates, there is an intricate crossover from image-charge attraction to the metal substrate (for large distances) to image-charge repulsion from the SAM (for small distances) with a stable minimum at a distance of roughly 20 times the thickness of the hydrophobic film. For bilayers in water, the interaction of a single charge is always repulsive. ii) The surface potential of a SAM is calculated for the case when the hydrophobic layer contains dipole moments, which might explain the recently observed long-ranged repulsion of hydrophobic scanning tips from PEG-terminated SAMs on gold. iii) The interaction between charged particles through the bilayer is weakened. Oppositely charged particles still attract each other through the membrane. The free-energy minimum occurs as a result of the competition between self-repulsion from the slab and interparticle attraction and is located at a separation from the membrane surface which equals 15 times the membrane thickness. iv) Surface charges on the two surfaces of a bilayer attract each other through the bilayer unless the surface charge densities are the same, even if the signs are the same. v) All these effects are strongly influenced by the presence of salt. Received 25 January 2000     

Calculating potential of mean force between like-charged nanoparticles: A comprehensive study on salt effects

Ions are critical to the structure and stability of polyelectrolytes such as nucleic acids. In this work, we systematically calculated the potentials of mean force between two like-charged nanoparticles in salt solutions by Monte Carlo simulations. The pseudo-spring method is employed to calculate the potential of mean force and compared systematically with the inversed-Boltzmann method. An effective attraction is predicted between two like-charged nanoparticles in divalent/trivalent salt solution and such attraction becomes weakened at very high salt concentration. Our analysis reveals that for the system, the configuration of ion-bridging nanoparticles is responsible for the attraction, and the invasion of anions into the inter-nanoparticles region at high salt concentration would induce attraction weakening rather than the charge inversion effect. The present method would be useful for calculating effective interactions during nucleic acid folding.     

Debye-Hückel theory for interfacial geometries

The Debye-Hückel theory for bulk electrolyte solutions is generalized to planar interfacial geometries, including screening effects due to mobile salt ions which are confined to the interface and solutions with in general different salt concentrations and dielectric constants on the two sides of the interface. We calculate the general Debye-Hückel interaction between fixed test charges, and analyze a number of relevant special cases as applicable to charged colloids and charged polymers. Salty interfaces, which are experimentally realized by monolayers or bilayers made of cationic and anionic surfactants or lipids, exert a strong attraction on charged particles of either sign at large separations from the interface; at short distances image-charge repulsion sets in. Likewise, the effective interactions between charged particles are strongly modified in the neighborhood of such a salty interface. On the other hand, charged particles which are immersed in a salt solution are repelled from the air (or a substrate) interface, and the interaction between two charges decays algebraically close to such an interface. These general results have experimentally measurable consequences for the adsorption of charged colloids or charged polymers at monolayers, solid substrates, and interfaces.     

Nonlinear effects in the electrophoresis of a spherical colloidal particle

The reduced electrophoretic mobility-reduced zeta potential relationship for a charged macroparticle is shown to be nonuniversal and to be highly nonlinear. In agreement with experimental results, a mobility reversal due to the macroion's charge inversion and a nonlinear dependence of the mobility on salt concentration is obtained.     

Simulation of charge reversal in salty environments: Giant overcharging?

We have performed MD simulations of a highly charged colloid in a solution of 3:1 and additional 1:1 salt. The dependency of the colloid's inverted charge on the concentration of the additional 1:1 salt has been studied. Most theories predict, that the inverted charge increases when the concentration of monovalent salt grows, up to what is called giant overcharging, while experiments and simulational studies observe the opposite. Our simulations agree with the experimental findings and shed light onto the weaknesses of the theories.     

Hierarchical model for the scale-dependent velocity of waves in random media

Elastic waves of short wavelength propagating through the upper layer of the Earth appear to move faster at large separations of source and receiver than at short separations. Existing perturbation theories predict a linear increase of the velocity shift with increasing separation and cannot describe the saturation of the velocity shift at large separations that is seen in computer simulations. We point out that this nonperturbative problem can be solved using a model developed originally for the study of directed polymers. The saturation velocity is found to scale with the four-thirds power of the root-mean-square amplitude of the velocity fluctuations, in good agreement with the computer simulations.     

Incommmensurability and unconventional superconductor to insulator transition in the hubbard model with bond-charge interaction

We determine the quantum phase diagram of the one-dimensional Hubbard model with bond-charge interaction X in addition to the usual Coulomb repulsion U>0 at half-filling. For large enough X相似文献   

Microscopic membrane elasticity and interactions among membrane inclusions: interplay between the shape, dilation, tilt and tilt-difference modes

A phenomenological Landau elasticity for the shape, dilation, and lipid-tilt of bilayer membranes is developed. The shape mode couples with the sum of the monolayers' tilt, while the dilation mode couples with the difference of the monolayers' tilts. Interactions among membrane inclusions within regular arrays are discussed. Inclusions modifying the membrane thickness and/or inducing a tilt-difference due to their convex or concave shape yield a dilation-induced attraction and a tilt-difference-induced repulsion. The resulting interaction can stabilize 2D crystal phases, with the possible coexistence of different lattice spacings when the dilation-tilt-difference coupling is large. Inclusions favoring crystals are those with either a long-convex or a short-concave hydrophobic core. Inclusions inducing a local membrane curvature due to their conical shape repel one another. At short inclusions separations, a tilt comparable with the inclusion's cone angle develops: it relaxes the membrane curvature and reduces the repulsion. At large separations the tilt vanishes, whatever the value of the shape-tilt coupling. Received 23 October 1998 and Received in final form 12 January 1999     

Role of quantum statistics in multi-particle decay dynamics

The role of quantum statistics in the decay dynamics of a multi-particle state, which is suddenly released from a confining potential, is investigated. For an initially confined double particle state, the exact dynamics is presented for both bosons and fermions. The time-evolution of the probability to measure two-particle is evaluated and some counterintuitive features are discussed. For instance, it is shown that although there is a higher chance of finding the two bosons (as oppose to fermions, and even distinguishable particles) at the initial trap region, there is a higher chance (higher than fermions) of finding them on two opposite sides of the trap as if the repulsion between bosons is higher than the repulsion between fermions. The results are demonstrated by numerical simulations and are calculated analytically in the short-time approximation. Furthermore, experimental validation is suggested.     

四臂星型聚电解质刷在外加三价盐离子和电场下的塌缩

本文采用Langevin动力学模拟的方法,研究存在三价盐离子时,接枝在带相反电荷的极板上的部分带电的四个臂星型链呈现出的塌缩现象. 结果表明,在电场作用下,接枝星型链的平均带电分数和盐浓度在带电单体和三价盐离子的竞争性吸附中起关键作用. 对于接枝在带相反电荷的极板上的带电分数较高的星型链,刷子会塌缩到接枝极板上,并会产生极板表面电荷的过度补偿现象. 当带电分数较低时,如果星型链所带电荷数与三价盐离子电荷数相同,即使在很低的盐浓度下,极板对三价盐离子的吸引能力也高于对星型链中的带电单体的吸引. 结果表明,星型链在带电分数较低的情况下,三价盐离子的加入不会导致接枝电极表面电荷的过度补偿. 此外,本文还研究了三价盐离子对电场作用下星型刷的拉伸的影响.     

Hexagonal to square lattice conversion in bilayer systems

We report the results of extensive molecular dynamics simulations of the reconstructive hexagonal to square lattice conversion in bilayer colloid systems. Two types of interparticle potential were used to represent the colloid-colloid interactions in the suspension. One potential, due to Marcus and Rice, is designed to describe the interaction of sterically stabilized colloid particles. This potential has a term that represents the attraction between colloid particles when there is incipient overlap between the stabilizing brushes on their surfaces, a (soft repulsion) term that represents the entropy cost associated with interpenetration of the stabilizing brushes, and a term that represents core-core repulsion. The other potential we used is an almost hard core repulsion with continuous derivatives. Our results clearly show that the character of the reconstructive hexagonal to square lattice conversion in bilayer colloid systems is potential dependent. For a system with colloid-colloid interactions of the Marcus-Rice type, the packing of particles in the square array exhibits a large interlayer lattice spacing, with the particles located at the minima of the attractive well. In this case the hexagonal to square lattice transition is first order. For a system with hard core colloid-colloid interactions there are two degenerate stable intermediate phases, linear and zigzag rhombic, that are separated from the square lattice by strong first order transitions, and from the hexagonal lattice by either weak first or second order transitions.     

Donor ordering and the statistics of quantum dot level spacings and widths from full 3D self-consistent electronic structure

Using full 3D self-consistent electronic structure calculations of small (electron numberN 100) lateral quantum dots formed on GaAs–AlGaAs HEMT devices we calculate the statistics of level spacings Δε_pand tunneling coefficients Γ_pbetween leads and confined states of the dot. We employ random and ordered donor layer charge distributions, the latter generated through Monte Carlo variable range hopping simulations, as well as a homogeneous (jellium) ionic charge distribution, and examine the effects on these statistics.It has recently been argued that the statistics of the level spacings and widths follow from random matrix theory when the Hamiltonian is described by the Gaussian orthogonal ensemble (GOE) for zero magnetic fieldB, and by the Gaussian unitary ensemble (GUE) forBsufficiently large to break time reversal symmetry. Specifically it is argued that when the dot wave functions are expanded in an arbitrary basis the expansion coefficients, according to the postulate of Porter and Thomas, are uniformly distributed in Hilbert space.In our calculation we obtain statistics of level spacings and widths by generating many configurations of disordered and ordered donor charge. This corresponds to the experimental situation of thermal cycling of the device. We find that a pronounced transition occurs in the level spacing statistics between the completely disordered donor layer ensemble, which seems to be well described by random matrix theory, and the ordered ensemble which is dominated by secular variations in the coefficients. In particular, a shell structure in the levels, which results from approximate parabolicity in the self-consistent confining potential, is observed. This, and the effects of symmetry under inversion and azimuthal symmetry, are speculated to undermine level repulsion and result in Poisson statistics for the levels here at the band edge.Finally we find that distortions in the dot shape are markedly less significant in varying the widths (and level spacings) than calculations based on a hard wall potential for the dot predict. This suggests that the notion of invariant atomic structure for sufficiently small dots is not invalidated by the randomness inherent in donor positions and shape distortion but, on the contrary, a systematic study of dot structure is possible.