General Differential Contact Identities for Macromolecules.

We discuss general Maxwell identities relating a macromolecule's charge, the forces acting at its surface, and the osmotic pressure of the solution in which it sits. The identities are closely related to the contact value relations that hold for certain special geometries, but are more general. In particular, the Maxwell identities can be applied to any macromolecule geometry, and they hold both within and outside of mean-field theory. Examples illustrate that combining the identities with approximate treatments of screening can often return simple, accurate osmotic pressure estimates.

Effect of charge inhomogeneity and mobility on colloid aggregation.

The aggregation of inhomogeneously charged colloids with the same average charge is analyzed using Monte Carlo simulations. We find aggregation of colloids for sizes in the range 10-200 nm, which is similar to the range in which aggregation is observed in several experiments. The attraction arises from the strongly correlated electrostatic interactions associated with the increase in the counterion density in the region between the particles; this effect is enhanced by the discreteness and mobility of the surface charges. Larger colloids attract more strongly when their surface charges are discrete. We study the aggregation as functions of the surface charge density, counterion valence, and volume fraction.

Long-range interaction between heterogeneously charged membranes.

Despite their neutrality, surfaces or membranes with equal amounts of positive and negative charge can exhibit long-range electrostatic interactions if the surface charge is heterogeneous; this can happen when the surface charges form finite-size domain structures. These domains can be formed in lipid membranes where the balance of the different ranges of strong but short-ranged hydrophobic interactions and longer-ranged electrostatic repulsion result in a finite, stable domain size. If the domain size is large enough, oppositely charged domains in two opposing surfaces or membranes can be strongly correlated by the electrostatic interactions; these correlations give rise to an attractive interaction of the two membranes or surfaces over separations on the order of the domain size. We use numerical simulations to demonstrate the existence of strong attractions at separations of tens of nanometers. Large line tensions result in larger domains but also increase the charge density within the domain. This promotes correlations and, as a result, increases the intermembrane attraction. On the other hand, increasing the salt concentration increases both the domain size and degree of domain anticorrelation, but the interactions are ultimately reduced due to increased screening. The result is a decrease in the net attraction as salt concentration is increased.

Monte carlo simulations of tau proteins: effect of phosphorylation.

We perform Monte Carlo simulations of tau proteins bound to a cylinder that mimics a microtubule (MT), and then study them in solution. Tau protein binds to a highly anionic MT surface to stabilize the cylindrical structure of MT. The negatively charged tail domain floats away from the anionic MT surface while positively charged tau segments localize near the MT surface. Monte Carlo simulations demonstrate that, in 3RS tau isoform (which has three imperfect repeats (R) short (S) isoform), amino acids are more condensed near a highly charged interface compared to 4RL isoform (which has four imperfect repeats (R) long (L) isoform). In 4RL isoform, amino acids in tail domain stay mostly apart from the MT surface. In the bulk solution, dephosphorylated taus are separated due to Coulomb repulsion between similarly charged isoforms. Moderate phosphorylation of 3RS isoform decreases average intermolecular distance between dephosphorylated and phosphorylated taus and lead to their overlap. Further phosphorylation does not change noticeably the intermolecular distances.

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.

Hybrid lipids as a biological surface-active component.

Cell membranes contain small domains (on the order of nanometers in size, sometimes called rafts) of lipids whose hydrocarbon chains are more ordered than those of the surrounding bulk-phase lipids. Whether these domains are fluctuations, metastable, or thermodynamically stable, is still unclear. Here, we show theoretically how a lipid with one saturated hydrocarbon chain that prefers the ordered environment and one partially unsaturated chain that prefers the less ordered phase, can act as a line-active component. We present a unified model that treats the lipids in both the bulk and at the interface and show how they lower the line tension between domains, eventually driving it to zero at sufficiently large interaction strengths or at sufficiently low temperatures. In this limit, finite-sized domains stabilized by the packing of these hybrid lipids can form as equilibrium structures.

A numerical study of the electrostatic properties of two finite-width charged dielectric slabs in water.

A numerical algorithm based on the image charge method is introduced to calculate the electrostatic potential, energy, and forces present in systems involving multiple point charges embedded in an inhomogeneous dielectric environment composed of five parallel dielectric slabs. The methodology is implemented within Monte Carlo simulations to calculate the thermal properties of two charged dielectric plates of finite thickness immersed in water.

Corrections to the Saffman-Delbruck mobility for membrane bound proteins.

Recent experiments (Gambin, Y., R. Lopez-Esparza, M. Reffay, E. Sierecki, N. S. Gov, M. Genest, R. S. Hodes, and W. Urbach. 2006. Proc. Natl. Acad. Sci. USA. 103:2098-2102) have called into question the applicability of the Saffman-Delbrück diffusivity for proteins embedded in the lipid bilayers. We present a simple argument to account for this observation that should be generically valid for a large class of transmembrane and membrane bound proteins. Whenever the protein-lipid interactions locally deform the membrane, that deformation generates new hydrodynamic stresses on the protein-membrane complex leading to a suppression of its mobility. We show that this suppression depends on the protein size in a manner consistent with the work of Gambin et al.

Effects of dielectric discontinuities on two charged plates.

Counterions in a biological system are charged in water and interact with charged macroions, which are generally made up of hydrocarbons. The dielectric difference between water and the hydrocarbon substrates occurs naturally, and may greatly affect the electrostatic properties of biological systems. Particularly for a slab geometry, bulk counterions that are dissolved in water are driven to the midplane of the slab because of their repulsive interaction with their image charges. The pressure between two charged plates becomes less repulsive since the low dielectric constant of the hydrocarbon substrate creates stronger association between counterions and surface charges as compared to the case of no dielectric discontinuity.

Charge Renormalization for Ellipsoidal Macroions.

We study the problem of counterion condensation for ellipsoidal macroions, a geometry well-suited for modeling liquid crystals, anisotropic vesicles, and polymers. We find that the ions within an ellipsoid's condensation layer are relatively unrestricted in their motions, and consequently work to establish a quasi-equipotential at its surface. This simplifies the application of Alexander et al.'s procedure, enabling us to obtain accurate analytic estimates for the critical valence of a general ellipsoid in the weak screening limit. Interestingly, we find that the critical valence of an eccentric ellipsoid is always larger than that of the sphere of equal volume, implying that counterion condensation provides a force resisting the deformation of spherical macroions. This contrasts with a recent study of flexible spherical macroions, which observed a preference for deformation into flattened shapes when considering only linear effects. Our work suggests that the balance of these competing forces might alter the nature of the transition.

The role of multipoles in counterion-mediated interactions between charged surfaces: strong and weak coupling.

We present general arguments for the importance, or lack thereof, of structure in the charge distribution of counterions for counterion-mediated interactions between bounding symmetrically charged surfaces. We show that on the mean field or weak coupling level, the charge quadrupole contributes the lowest order modification to the contact value theorem and thus to the intersurface electrostatic interactions. The image effects are non-existent on the mean field level even with multipoles. On the strong coupling level the quadrupoles and higher order multipoles contribute additional terms to the interaction free energy only in the presence of dielectric inhomogeneities. Without them, the monopole is the only multipole that contributes to the strong coupling electrostatics. We explore the consequences of these statements in all their generality.

Attractions between like-charged surfaces with dumbbell-shaped counterions.

We study the effect of dumbbell-like counterions on the interactions between similarly charged surfaces. Via a systematic study using Monte Carlo simulations and field theory, we fully consider electrostatic correlations and ion structure and find that their intricate coupling determines the equilibrium phase behaviors. In particular, an energetic bridging mechanism is revealed to cause surface attractions for a finite range of surface separations, even in the Poisson-Boltzmann limit.

Self assembly modulated by interactions of two heterogeneously charged surfaces.

Recent experiments have measured attractive interactions between two surfaces that each bear two molecular species with opposite charge. Such surfaces form charged domains of finite size. We present a theoretical model that predicts the dependence of the domain size, phase behavior and the interlayer forces as a function of spacing and salt concentration for two such interacting surfaces. A strong correlation between two length scales, the screening length and the surface separation, at the spinodal is shown. Remarkably, the first-order phase transition to infinite sized domains depends logarithmically on the ratio of the domain size to the molecular size. Finally, we fit the predicted pressure with experiments.

Strong-coupling electrostatics in the presence of dielectric inhomogeneities.

We study the strong-coupling (SC) interaction between two like-charged membranes of finite thickness embedded in a medium of higher dielectric constant. A generalized SC theory is applied along with extensive Monte Carlo simulations to study the image charge effects induced by multiple dielectric discontinuities in this system. These effects lead to strong counterion crowding in the central region of the intersurface space upon increasing the solvent-membrane dielectric mismatch and change the membrane interactions from attractive to repulsive at small separations. These features agree quantitatively with the SC theory at elevated couplings or dielectric mismatch where the correlation hole around counterions is larger than the thickness of the central counterion layer.

Equilibrium domains on heterogeneously charged surfaces.

Recent experiments have shown that salt solutions containing surfaces with two oppositely charged species show stable, possibly equilibrium, structures with finite domain sizes. The short-range interactions between the two species would normally result in phase separation that is driven by the line tension with macroscopically large domains of each species. In this paper, we show that, when at least one of the charged species is mobile, finite domains can occur in equilibrium. The domain size is determined by a competition of the electrostatic free energy that promotes charge mixing and small domains, with the line tension that promotes macroscopic phase separation. We calculate the equilibrium patch size as a function of the surface charge and the concentration of dissolved monovalent salts in the bulk phase. An important finding is the prediction of a first-order transition from finite patches to macroscopic phase separation of the two charge species as the salt concentration is increased.

Optical propagation in laboratory-generated turbulence.

A versatile and useful facility for simulating the effects of atmospheric turbulence on optical propagation is described, and the relevant system parameters are characterized. The scattering medium is a turbulent liquid (ethanol) with the turbulence created by unstable convection generated by a strong vertical thermal gradient. Measurements of the structure function and the spatial spectrum of the resulting refractive index fluctuations are presented and compared with the theoretically predicted forms. The effects of this scattering medium on laser beam propagation are determined and compared to the first-order Rytov theory. In particular, probability density functions, moments, and spatial covariance functions of the irradiance resulting from propagation through the system with a variety of turbulence levels and path lengths are presented.

Time delayed covariance of the received intensity of a monochromatic speckle pattern in the turbulent atmosphere.

An approximate numerical evaluation of the previously formulated time delayed covariance of the received intensity of a laser speckle pattern is described and the results are compared with experimental data.

Covariance of the received intensity of a partially coherent laser speckle pattern in the turbulent atmosphere.

A theoretical expression for the covariance of the received intensity of a partially coherent laser speckle pattern after propagation through the turbulent atmosphere is developed. It is shown that the atmospheric perturbation on a partially coherent speckle pattern can be decomposed into a coherent term and an incoherent term. The dependence of contributions of these components on the level of turbulence, vacuum speckle-contrast ratio, and detector spacing is studied in detail and the results are compared with the available experimental data.