Competing interactions in two dimensional Coulomb systems: Surface charge heterogeneities in coassembled cationic-anionic incompatible mixtures

The long-range attraction between hydrophobic macroscopic surfaces

Direct measurements of the long-range strongly attractive force observed between macroscopic hydrophobic surfaces across aqueous solutions are reexamined in light of recent experiments and theoretical developments. The focus is on systems in the absence of submicroscopic bubbles (preexistent or induced) to avoid capillary bridging forces. Other possible interferences to the measurements are also eliminated. The force-distance profiles are obtained directly (no contributions from electrical double layer or hydrodynamics) between symmetric identical hydrophobic surfaces, overall charge-neutral, at the thermodynamic equilibrium and in a quenched state. Therefore in the well-defined geometry of crossed-cylinders, sphere-flat, or sphere-sphere, there is no additional interaction to be considered except the ever-present dispersion forces, negligible at large separations. For the three main categories of substrates rendered hydrophobic, namely surfaces obtained with surfactant monolayers physically adsorbed from solution to deposited ones, and substrates coated with a hydrophobizing agent bonded chemically onto the surface, the interaction energy scales as A exp (-2κD)/2κD at large separations, with measured decay lengths in accord with theoretical predictions, simply being half the Debye screening length, κ^-1/2, at least in non vanishing electrolyte. Taken together with the prefactor A scaling as the ionic strength, the interaction energy is demonstrated to have an electrostatic origin in all the systems. Thanks to our recent SFAX coupling force measurements with x-ray solution scattering under controlled nano-confinement, the microstructuration of the adsorbed film emerges as an essential feature in the molecular mechanism for explaining the observed attraction of larger magnitude than dispersion forces. The adsorption of pairs of positive and negative ions on small islands along the interface, the fluctuation of the surface charge density around a zero mean-value with desorption into or adsorption from the electrolyte solution, the correlations in the local surface ion concentrations along the surfaces, the redistribution of counterions upon intersurface variation, all contribute and are tuned finely by the inhomogeneities and defects present in the hydrophobic layers. It appears that the magnitude of the interacting energy can be described by a single master curve encompassing all the systems.

Direct Measurement of Lateral Correlations under Controlled Nanoconfinement

Lateral correlations along hydrophobic surfaces whose separation can be varied continuously are measured by x-ray scattering using a modified surface force apparatus coupled with synchrotron radiation, named SFAX. A weak isotropic diffuse scattering along the equatorial plane is revealed for mica surfaces rendered hydrophobic and charge neutral by immersion in cationic surfactant solutions at low concentrations. The peak corresponds to a lateral surface correlation length ξ≈12  nm, without long-range order. These findings are compatible with the atomic force microscopy imaging of a single surface, where adsorbed surfactant stripes appear surrounded by bare mica zones. Remarkably, the scattering patterns remain stable for gap widths D larger than the lateral period but change in intensity and shape (to a lesser extent) as soon as D<ξ. This evolution codes for a redistribution of counterions (counterion release from antagonistic patches) and the associated new x-ray labeling of the patterns. The redistribution of counterions is also the key mechanism to the long-range electrostatic attraction between similar, overall charge-neutral walls, reported earlier.

Electrostatics of patchy surfaces

In the study of colloidal, biological and electrochemical systems, it is customary to treat surfaces, macromolecules and electrodes as homogeneously charged. This simplified approach is proven successful in most cases, but fails to describe a wide range of heterogeneously charged surfaces commonly used in experiments. For example, recent experiments have revealed a long-range attraction between overall neutral surfaces, locally charged in a mosaic-like structure of positively and negatively charged domains ("patches"). Here, we review experimental and theoretical studies addressing the stability of heterogeneously charged surfaces, their effect on ionic profiles in solution, and the interaction between two such surfaces. We focus on electrostatics, and highlight the important new physical parameters appearing in the heterogeneous case, such as the largest patch size and inter-surface charge correlations.

Two-Dimensional Disorder in Black Phosphorus and Monochalcogenide Monolayers

Ridged, orthorhombic two-dimensional atomic crystals with a bulk Pnma structure such as black phosphorus and monochalcogenide monolayers are an exciting and novel material platform for a host of applications. Key to their crystallinity, monolayers of these materials have a 4-fold degenerate structural ground state, and a single energy scale EC (representing the elastic energy required to switch the longer lattice vector along the x- or y-direction) determines how disordered these monolayers are at finite temperature. Disorder arises when nearest neighboring atoms become gently reassigned as the system is thermally excited beyond a critical temperature Tc that is proportional to EC/kB. EC is tunable by chemical composition and it leads to a classification of these materials into two categories: (i) Those for which EC ≥ kBTm, and (ii) those having kBTm > EC ≥ 0, where Tm is a given material's melting temperature. Black phosphorus and SiS monolayers belong to category (i): these materials do not display an intermediate order-disorder transition and melt directly. All other monochalcogenide monolayers with EC > 0 belonging to class (ii) will undergo a two-dimensional transition prior to melting. EC/kB is slightly larger than room temperature for GeS and GeSe, and smaller than 300 K for SnS and SnSe monolayers, so that these materials transition near room temperature. The onset of this generic atomistic phenomena is captured by a planar Potts model up to the order-disorder transition. The order-disorder phase transition in two dimensions described here is at the origin of the Cmcm phase being discussed within the context of bulk layered SnSe.

Pickering emulsions stabilized by oppositely charged colloids: Stability and pattern formation

A binary mixture of oppositely charged colloids can be used to stabilize water-in-oil or oil-in-water emulsions. A Monte Carlo simulation study to address the effect of charge ratio of colloids on the stability of Pickering emulsions is presented. The colloidal particles at the interface are modeled as aligned dipolar hard spheres, with attractive interaction between unlike-charged and repulsive interaction between like-charged particles. The optimum composition (fraction of positively charged particles) required for the stabilization corresponds to a minimum in the interaction energy per particle. In addition, for each charge ratio, there is a range of compositions where emulsions can be stabilized. The structural arrangement of particles or the pattern formation at the emulsion interface is strongly influenced by the charge ratio. We find well-mixed isotropic, square, and hexagonal arrangements of particles on the emulsion surface for different compositions at a given charge ratio. The distribution of coordination numbers is calculated to characterize structural features. The simulation study is useful for the rational design of Pickering emulsifications wherein oppositely charged colloids are used, and for the control of pattern formation that can be useful for the synthesis of colloidosomes and porous shells derived thereof.

Theory of non-equilibrium critical phenomena in three-dimensional condensed systems of charged mobile nanoparticles

A study of 3d electrostatic self-assembly (SA) in systems of charged nanoparticles (NPs) is one of the most difficult theoretical problems. In particular, the limiting case of negligible or very low polar media (e.g. salt) concentration, where the long-range NP interactions cannot be reduced to commonly used effective short-range (Yukawa) potentials, remains unstudied. Moreover, the present study has demonstrated that unlike the Debye-Hückel theory, a complete screening of the charges in SA kinetics (dynamic SA) is not always possible. Generally speaking, one has to take into account implicitly how each NP interacts with all other NPs (the true long-range interactions). Traditional theoretical methods allow us to monitor such electrostatic 3d system kinetics only for very short times, which is far from sufficient for understanding the dynamic SA. In this paper, combining an integrated analytical approach (the non-linear integro-differential kinetic equation for correlation functions) and reverse Monte Carlo in the 3d case, we have obtained a self-consistent solution of this challenging problem. We demonstrate, in particular, the existence of critical points and critical phenomena in the non-equilibrium kinetics in a 3d system of oppositely charged mobile NPs.

Macroscopic lateral heterogeneity observed in a laterally mobile immiscible mixed polyelectrolyte-neutral polymer brush

We studied mixed poly(ethylene oxide) (PEO) and poly(2-(dimethylamino)ethyl methacrylate) (PDMAEMA) brushes. The question we attempted to answer was: when the chain grafting points are laterally mobile, how will this lateral mobility influence the structure and phase behavior of the mixed brush? Three different model mixed PEO/PDMAEMA brush systems were prepared: (1) a laterally mobile mixed brush by spreading onto the air-water interface a mixture of poly(ethylene oxide)-poly(n-butyl acrylate) (PEO-PnBA) and poly(2-(dimethylamino)ethyl methacrylate)-poly(n-butyl acrylate) (PDMAEMA-PnBA) diblock copolymers (the specific diblock copolymers used will be denoted as PEO113-PnBA100 and PDMAEMA118-PnBA100, where the subscripts refer to the number-average degrees of polymerization of the individual blocks), (2) a mobility-restricted (inseparable) version of the above mixed brush prepared using a PEO-PnBA-PDMAEMA triblock copolymer (denoted as PEO113-PnBA89-PDMAEMA120) having respective brush molecular weights matched with those of the diblock copolymers, and (3) a different laterally mobile mixed PEO and PDMAEMA brush prepared from a PEO113-PnBA100 and PDMAEMA200-PnBA103 diblock copolymer combination, which represents a further more height-mismatched mixed brush situation than described in (1). These three mixed brush systems were investigated by surface pressure-area isotherm and X-ray (XR) reflectivity measurements. These experimental data were analyzed within the theoretical framework of a continuum self-consistent field (SCF) polymer brush model. The combined experimental and theoretical results suggest that the mobile mixed brush derived using the PEO113-PnBA100 and PDMAEMA118-PnBA100 combination (i.e., mixed brush System #1) undergoes a lateral macroscopic phase separation at high chain grafting densities, whereas the more height-mismatched system (System #3) is only microscopically phase separated under comparable brush density conditions even though the lateral mobility of the grafted chains is unrestricted. The macroscopic phase separation observed in the laterally mobile mixed brush system is in contrast with the microphase separation behavior commonly observed in two-dimensional laterally mobile charged small molecule mixtures. Further study is needed to determine the detailed morphologies of the macro- and microphase-separated mixed PEO/PDMAEMA brushes.

Static and dynamic screening effects in the electrostatic self-assembly of nano-particles

3D structure of self-assembled and oppositely charged nanoparticles.

Electric Field Controlled Self-Assembly of Hierarchically Ordered Membranes

Self-assembly in the presence of external forces is an adaptive, directed organization of molecular components under nonequilibrium conditions. While forces may be generated as a result of spontaneous interactions among components of a system, intervention with external forces can significantly alter the final outcome of self-assembly. Superimposing these intrinsic and extrinsic forces provides greater degrees of freedom to control the structure and function of self-assembling materials. In this work we investigate the role of electric fields during the dynamic self-assembly of a negatively charged polyelectrolyte and a positively charged peptide amphiphile in water leading to the formation of an ordered membrane. In the absence of electric fields, contact between the two solutions of oppositely charged molecules triggers the growth of closed membranes with vertically oriented fibrils that encapsulate the polyelectrolyte solution. This process of self-assembly is intrinsically driven by excess osmotic pressure of counterions, and the electric field is found to modify the kinetics of membrane formation, and also its morphology and properties. Depending on the strength and orientation of the field we observe a significant increase or decrease of up to nearly 100% in membrane thickness, as well as the controlled rotation of nanofiber growth direction by 90 degrees, resulting in a significant increase in mechanical stiffness. These results suggest the possibility of using electric fields to control structure in self-assembly processes involving diffusion of oppositely charged molecules.

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.

Microscopic approach to the kinetics of pattern formation of charged molecules on surfaces

A microscopic formalism based on computing many-particle densities is applied to the analysis of the diffusion-controlled kinetics of pattern formation in oppositely charged molecules on surfaces or adsorbed at interfaces with competing long-range Coulomb and short-range Lennard-Jones interactions. Particular attention is paid to the proper molecular treatment of energetic interactions driving pattern formation in inhomogeneous systems. The reverse Monte Carlo method is used to visualize the spatial molecular distribution based on the calculated radial distribution functions (joint correlation functions). We show the formation of charge domains for certain combinations of temperature and dynamical interaction parameters. The charge segregation evolves into quasicrystalline clusters of charges, due to the competing long- and short-range interactions. The clusters initially co-exist with a gas phase of charges that eventually add to the clusters, generating "fingers" or line of charges of the same sign, very different than the nanopatterns expected by molecular dynamics in systems with competing interactions in two dimensions, such as strain or dipolar versus van der Waals interactions.

Morphologies of planar polyelectrolyte brushes in a poor solvent: molecular dynamics simulations and scaling analysis

Using molecular dynamics simulations and scaling analysis, we study the effect of the solvent quality for the polymer backbone, the strength of the electrostatic interactions, the chain degree of polymerization, and the brush grafting density on conformations of the planar polyelectrolyte brushes in salt-free solutions. Polyelectrolyte brush forms: (1) vertically oriented cylindrical aggregates (bundles of chains), (2) maze-like aggregate structures, or (3) thin polymeric layer covering a substrate. These different brush morphologies appear as a result of the fine interplay between electrostatic and short-range monomer-monomer interactions. The brush thickness shows nonmonotonic dependence on the value of the Bjerrum length. It first increases with the increasing value of the Bjerrum length, and then it begins to decrease. This behavior is a result of counterion condensation within a brush volume.

Surface patterning of low-dimensional systems: the chirality of charged fibres

Charged surfaces are interesting for their ability to have long-range correlations and their ability to be dynamically tuned. While the configurations of charged planar surfaces have been thoroughly mapped and studied, charged cylindrical surfaces show novel features. The surface patterning of cylindrically confined charges is discussed with emphasis on the role of chiral configurations. The origins of surface patterns due to competing interactions in charged monolayers are summarized along with their associated theoretical models. The electrostatically induced patterns described in this paper are important in many low-dimensional biological systems such as plasma membrane organization, filamentous virus capsid structure or microtubule interactions. A simple model effectively predicting some features of chiral patterns in biological systems is presented. We extend our model from helical lamellar patterns to elliptical patterns to consider asymmetrical patterns in assemblies of filamentous aggregates.

Faceting ionic shells into icosahedra via electrostatics

Shells of various viruses and other closed packed structures with spherical topology exhibit icosahedral symmetry because the surface of a sphere cannot be tiled without defects, and icosahedral symmetry yields the most symmetric configuration with the minimum number of defects. Icosahedral symmetry is different from icosahedral-shaped structures, which include some large viruses, cationic-anionic vesicles, and fullerenes. We present a faceting mechanism of ionic shells into icosahedral shapes that breaks icosahedral symmetry resulting from different arrangements of the charged components among the facets. These self-organized ionic structures may favor the formation of flat domains on curved surfaces. We show that icosahedral shapes without rotational symmetry can have lower energy than spheres with icosahedral symmetry caused by preferred bending directions in the planar ionic lattice. The ability to create icosahedral shapes without icosahedral symmetry may lead to the design of new functional materials. The electrostatically driven faceting mechanism we present here suggests that we can design faceted polyhedra with diverse symmetries by coassembling oppositely charged molecules of different stoichiometric ratios.

Spontaneous chirality via long-range electrostatic forces

We consider a model for periodic patterns of charges constrained over a cylindrical surface. In particular we focus on patterns of chiral helices, achiral rings, or vertical lamellae, with the constraint of global electroneutrality. We study the dependence of the patterns' size and pitch angle on the radius of the cylinder and salt concentration. We obtain a phase diagram by using numerical and analytic techniques. For pure Coulomb interactions, we find a ring phase for small radii and a chiral helical phase for large radii. At a critical salt concentration, the characteristic domain size diverges, resulting in an achiral macroscopic phase-segregated structure. We discuss possible consequences and generalizations of our model.

Charged particles on surfaces: coexistence of dilute phases and periodic structures at interfaces

We consider a mixture of two immiscible oppositely charged molecules strongly adsorbed to an interface, with a neutral nonselective molecular background. We determine the coexistence between a high density ionic periodic phase and a dilute isotropic ionic phase. We use a strong segregation approach for the periodic phase and determine the one-loop free energy for the dilute phase. Lamellar and hexagonal patterns are calculated for different charge stoichiometries of the mixture. Molecular dynamics simulations exhibit the predicted phase behavior. The periodic length scale of the solid phase is found to scale as epsilon/(lB psi3/2), where psi is the effective charge density, lB is the Bjerrum length, and epsilon is the cohesive energy.

Electrostatic attraction between cationic-anionic assemblies with surface compositional heterogeneities

Electrostatics play a key role in biomolecular assembly. Oppositely charged biomolecules, for instance, can be coassembled into functional units, such as DNA and histone proteins into nucleosomes and actin-binding protein complexes into cytoskeleton components, at appropriate ionic conditions. These cationic-anionic coassemblies often have surface charge heterogeneities that result from the delicate balance between electrostatics and packing constraints. Despite their importance, the precise role of surface charge heterogeneities in the organization of cationic-anionic coassemblies is not well understood. We show here that coassemblies with charge heterogeneities strongly interact through polarization of the domains. We find that this leads to symmetry breaking, which is important for functional capabilities, and structural changes, which is crucial in the organization of coassemblies. We determine the range and strength of the attraction as a function of the competition between the steric and hydrophobic constraints and electrostatic interactions.