Interaction potentials, structural ordering and effective charges in dispersions of charged colloidal particles

Theory of Charge Regulation of Colloidal Particles in Electrolyte Solutions

We present a theory that enables us to (i) calculate the effective surface charge of colloidal particles and (ii) efficiently obtain titration curves for different salt concentrations. The theory accounts for the shift of pH of solution due to the presence of 1:1 electrolyte. It also accounts self-consistently for the electrostatic potential produced by the deprotonated surface groups. To examine the accuracy of the theory, we have performed extensive reactive Monte Carlo simulations, which show excellent agreement between theory and simulations without any adjustable parameters.

Brownian dynamics simulations of shear-induced aggregation of charged colloidal particles in the presence of hydrodynamic interactions

HYPOTHESIS

In spite of the abundant literature on Brownian simulations of the aggregation behavior of colloidal suspensions both under quiescent conditions and in the presence of shear, few works performed simulations including the effect of hydrodynamic interactions. Even fewer works have investigated the effects of shear on the aggregation of electrostatically-stabilized colloidal suspensions.

SIMULATIONS

In this work, we employed Brownian dynamics simulations implementing the Rotne-Prager-Yamakawa approximation to account for hydrodynamic interactions and investigated the aggregation kinetics of electrostatically-stabilized colloidal suspensions exposed to simple shear, for various Péclet number values, particle volume fractions and surface potential values.

RESULTS

The increase in Péclet number (i.e., in the shear rate), leads to an overall increase in the aggregation rate and the formation of large aggregates that, for sufficiently high volume fractions, rapidly grow, leading to either breakup and restructuring phenomena or percolation of the system. In some cases, a bimodal distribution of the cluster population was observed. Our simulations further indicate that at the highest Péclet, the aggregation dynamics is independent of the energy barrier and entirely controlled by shear. A comparison with a simple BD method reveals that neglecting long-range hydrodynamic interactions leads to a substantial underestimation of the aggregation rate.

Aggregation behavior of nanoparticles: Revisiting the phase diagram of colloids

Surface functionalization of metal nanoparticles (NPs), e.g., using peptides and proteins, has recently attracted a considerable attention in the field of design of therapeutics and diagnostics. The possibility of diverse functionalization allows them to selectively interact with proteins, while the metal core ensures solubility, making them tunable therapeutic agents against diseases due to mis-folding or aggregation. On the other hand, their action is limited by possible self-aggregation, which could be, however, prevented based on the full understanding of their phase diagram as a function of the environmental variables (temperature, ionic strength of the solution, concentration) and intrinsic characteristics (size, charge, amount, and type of functional groups). A common modeling strategy to study the phase behavior is to represent the NPs as spheres interacting via effective potentials implicitly accounting for the solvation effects. Their size put the NPs into the class of colloids, albeit with particularly complex interactions including both attractive and repulsive features, and a consequently complex phase diagram. In this work, we review the studies exploring the phases of these systems starting from those with only attractive or repulsive interactions, displaying a simpler disperse-clustered-aggregated transitions. The phase diagram is here interpreted focusing on the universal aspects, i.e., those dependent on the general feature of the potentials, and available data are organized in a parametric phase diagram. We then consider the potentials with competing attractive short range well and average-long-range repulsive tail, better representing the NPs. Through the proper combination of the attractive only and repulsive only potentials, we are able to interpret the appearance of novel phases, characterized by aggregates with different structural characteristics. We identify the essential parameters that stabilize the disperse phase potentially useful to optimize NP therapeutic activity and indicate how to tune the phase behavior by changing environmental conditions or the NP chemical–physical properties.

Ferrofluids and bio-ferrofluids: looking back and stepping forward

Ferrofluids investigated along for about five decades are ultrastable colloidal suspensions of magnetic nanoparticles, which manifest simultaneously fluid and magnetic properties. Their magnetically controllable and tunable feature proved to be from the beginning an extremely fertile ground for a wide range of engineering applications. More recently, biocompatible ferrofluids attracted huge interest and produced a considerable increase of the applicative potential in nanomedicine, biotechnology and environmental protection. This paper offers a brief overview of the most relevant early results and a comprehensive description of recent achievements in ferrofluid synthesis, advanced characterization, as well as the governing equations of ferrohydrodynamics, the most important interfacial phenomena and the flow properties. Finally, it provides an overview of recent advances in tunable and adaptive multifunctional materials derived from ferrofluids and a detailed presentation of the recent progress of applications in the field of sensors and actuators, ferrofluid-driven assembly and manipulation, droplet technology, including droplet generation and control, mechanical actuation, liquid computing and robotics.

Influence of Semifluorinated Alkane Surface Domains on Phase Behavior and Linear and Nonlinear Viscoelasticity of Phospholipid Monolayers

Semifluorinated alkanes self-assemble into 30-40 nm-large surface domains (hemimicelles) at the air/water interface. They have been drawing increasing attention to stabilize microbubbles coated with lipids, which are used for enhancing the contrast in sonographic imaging. Although previous studies suggested that semifluorinated alkanes increase the stability of phospholipid membranes, little is known about how semifluorinated alkanes influence phase behaviors and mechanical properties of lipid-coated microbubbles. As a well-defined model of microbubble surfaces, we prepared monolayers consisting of a mixture of phospholipids and semifluorinated alkanes at the air/water interface and investigated the influence of hemimicelles of semifluorinated alkanes on the phase behavior and interfacial viscoelastic properties of phospholipid monolayers. Hemimicelles are phase-separated from phospholipids and accumulate at the phase boundary, which strongly modulates the correlation between solid phospholipid domains. Intringuingly, we found that the mixed monolayer of semifluorinated alkanes and phospholipids possesses linear and nonlinear viscoelastic properties comparable to those of phospholipid monolayers. Since the mixing of semifluorinated alkanes and phospholipids enables one to overcome the intrinsically low stability of pure semifluorinated alkanes against the change in the surface area of microbubbles through the partial dissolution of gas into the aqueous phase, this is a promising strategy for the stable coating of microbubbles in ultrasound diagnosis.

Theoretical predictions of structures in dispersions containing charged colloidal particles and non-adsorbing polymers

We develop a theoretical model to describe structural effects on a specific system of charged colloidal polystyrene particles, upon the addition of non-adsorbing PEG polymers. This system has previously been investigated experimentally, by scattering methods, so we are able to quantitatively compare predicted structure factors with corresponding experimental data. Our aim is to construct a model that is coarse-grained enough to be computationally manageable, yet detailed enough to capture the important physics. To this end, we utilize classical polymer density functional theory, wherein all possible polymer configurations are accounted for, subject to a mean-field Boltzmann weight. We make efforts to counteract drawbacks with this mean-field approach, resulting in structural predictions that agree very well with computationally more demanding simulations. Electrostatic interactions are handled at the fully non-linear Poisson-Boltzmann level, and we demonstrate that a linearization leads to less accurate predictions. The particle charge is an experimentally unknown parameter. We define the surface charge such that the experimental and theoretical gel point at equal polymer concentration coincide. Assuming a fixed surface charge for a certain salt concentration, we find very good agreements between measured and predicted structure factors across a wide range of polymer concentrations. We also present predictions for other structural quantities, such as radial distribution functions, and cluster size distributions. Finally, we demonstrate that our model predicts the occurrence of equilibrium clusters at high polymer concentrations, but low particle volume fractions and salt levels.

Sterically stabilized colloids with tunable repulsions

When studying tunable electrostatic repulsions in aqueous suspensions of charged colloids, irreversible colloid aggregation or gelation may occur at high salt concentrations. For many commonly used synthetic colloids, such as polystyrene and silica particles, the reason for coagulation is the presence of unbalanced, strongly attractive, and short-ranged van der Waals (VDW) forces. Here, we present an aqueous polystyrene model colloid that is sterically stabilized against VDW attractions. We show that the synthesis procedure, based on a neutral initiator couple and a nonionic surfactant, introduces surface charges that can be further increased by the addition of charged comonomer methacrylic acid. Thus, the interactions between the polystyrene spheres can be conveniently tuned from hard-sphere-like to charge-stabilized with long-ranged electrostatic repulsions described by a Yukawa-type pair potential. The particle size, grafting density, core-shell structure, and surface charge are characterized by light and neutron scattering. Using X-ray and neutron scattering in combination with an accurate analytic integral equation scheme for the colloidal static structure factor, we deduce effective particle charges for colloid volume fractions ≥0.1 and salt concentrations in the range of 1.5 to 50 mM.

Weak and strong coupling theories for polarizable colloids and nanoparticles

A theory is presented which allows us to accurately calculate the density profile of monovalent and multivalent counterions in suspensions of polarizable colloids or nanoparticles. In the case of monovalent ions, we derive a weak-coupling theory that explicitly accounts for the ion-image interaction, leading to a modified Poisson-Boltzmann equation. For suspensions with multivalent counterions, a strong-coupling theory is used to calculate the density profile near the colloidal surface and a Poisson-Boltzmann equation with a renormalized boundary condition to account for the counterion distribution in the far field. All the results are compared with the Monte Carlo simulations, showing an excellent agreement between the theory and the simulations.

Regulation of adhesion behavior of murine macrophage using supported lipid membranes displaying tunable mannose domains

Highly uniform, strongly correlated domains of synthetically designed lipids can be incorporated into supported lipid membranes. The systematic characterization of membranes displaying a variety of domains revealed that the equilibrium size of domains significantly depends on the length of fluorocarbon chains, which can be quantitatively interpreted within the framework of an equivalent dipole model. A mono-dispersive, narrow size distribution of the domains enables us to treat the inter-domain correlations as two-dimensional colloidal crystallization and calculate the potentials of mean force. The obtained results demonstrated that both size and inter-domain correlation can precisely be controlled by the molecular structures. By coupling α-D-mannose to lipid head groups, we studied the adhesion behavior of the murine macrophage (J774A.1) on supported membranes. Specific adhesion and spreading of macrophages showed a clear dependence on the density of functional lipids. The obtained results suggest that such synthetic lipid domains can be used as a defined platform to study how cells sense the size and distribution of functional molecules during adhesion and spreading.

Macroscopic ordering of polystyrene carboxylate-modified nanospheres self-assembled at the water-air interface

We present results from an experimental study of ordering characteristics in monolayers of polystyrene nanospheres self-assembled at a water-air interface. We demonstrate that the interaction of spheres, governed by the dissemination of surface charge, leads to the formation of macroscopic close-packed ordered areas or "domains" with a well-defined orientation of the lattice axes over areas of 25 mm(2). It was found that by changing the surface chemistry of the spheres it is possible to modify the balance between the attractive and repulsive forces and thus to control the ordering characteristics. We implemented a model that simulates the process of self-assembly and examines the ordering characteristics for layers with different ratio between attractive and repulsive forces. A good qualitative agreement was found between the simulations and experiment. These studies are technologically relevant as a method of producing nanosphere templates for large area patterned materials.

Polyelectrolyte-induced aggregation of liposomes: a new cluster phase with interesting applications

Different charged colloidal particles have been shown to be able to self-assemble, when mixed in an aqueous solvent with oppositely charged linear polyelectrolytes, forming long-lived finite-size mesoscopic aggregates. On increasing the polyelectrolyte content, with the progressive reduction of the net charge of the primary polyelectrolyte-decorated particles, larger and larger clusters are observed. Close to the isoelectric point, where the charge of the adsorbed polyelectrolytes neutralizes the original charge of the particles' surface, the aggregates reach their maximum size, while beyond this point any further increase of the polyelectrolyte-particle charge ratio causes the formation of aggregates whose size is progressively reduced. This re-entrant condensation behavior is accompanied by a significant overcharging. Overcharging, or charge inversion, occurs when more polyelectrolyte chains adsorb on a particle than are needed to neutralize its original charge so that, eventually, the sign of the net charge of the polymer-decorated particle is inverted. The stability of the finite-size long-lived clusters that this aggregation process yields results from a fine balance between long-range repulsive and short-range attractive interactions, both of electrostatic nature. For the latter, besides the ubiquitous dispersion forces, whose supply becomes relevant only at high ionic strength, the main contribution appears due to the non-uniform correlated distribution of the charge on the surface of the polyelectrolyte-decorated particles ('charge-patch' attraction). The interesting phenomenology shown by these system has a high potential for biotechnological applications, particularly when the primary colloidal particles are bio-compatible lipid vesicles. Possible applications of these systems as multi-compartment vectors for the simultaneous intra-cellular delivery of different pharmacologically active substances will be briefly discussed.

A simple phenomenological fix for the dielectric constant within the reference interaction site model approach

The effective interactions among ions immersed in water are studied by means of the effective pair potentials (EPPs) [J. Chem. Phys. 2002, 117, 6133] obtained after contracting (integrating out) the degrees of freedom of the solvent molecules. The dressed interaction site theory (DIST) leads to a simple way of adjusting the effective dielectric constant of the model solvent to its experimental value at standard conditions. The molecular structure of the solvent is mirrored in the structure of the short-ranged component of the induced EPPs, with noticeable differences between the cases with trivial (ideal gas) and nontrivial (experimental) values of the dielectric constant. The shape of these EPPs remains almost invariant over the whole range of salt concentrations considered here. The asymptotic behavior of the EPP between two macroions obtained after contracting the supporting electrolyte (water molecules plus small ions) is also briefly discussed.

Smoluchowski equation and the colloidal charge reversal

Smoluchowski equation and the Monte Carlo simulations are used to study the conditions leading to the reversal of the electrophoretic mobility. Zeta (zeta) potential is identified with the diffuse potential at the shear plane which, we argue, must be placed at least one ionic diameter away from the colloidal surface. For sufficiently strongly charged colloids, zeta potential changes sign as a function of the multivalent electrolyte concentration, resulting in a reversal of the electrophoretic mobility. This behavior occurs even for very small ions of 4 A diameter as long as the surface charge density of the colloidal particles is sufficiently large and the concentration of 1:1 electrolyte is sufficiently low.

Renormalization in charged colloids: non-monotonic behaviour with the surface charge

The static structure factor S(q) is measured for a set of deionized latex dispersions with different numbers of ionizable surface groups per particle and similar diameters. For a given volume fraction, the height of the main peak of S(q), which is a direct measure of the spatial ordering of latex particles, does not increase monotonically with the number of ionizable groups. This behaviour cannot be described using the classical renormalization scheme based on the cell model. We analyse our experimental data using a renormalization model based on the jellium approximation, which predicts the weakening of the spatial order for moderate and large particle charges.

Measuring inter-DNA potentials in solution

Interactions between short strands of DNA can be tuned from repulsive to attractive by varying solution conditions and have been quantified using small angle x-ray scattering techniques. The effective DNA interaction charge was extracted by fitting the scattering profiles with the generalized one-component method and inter-DNA Yukawa pair potentials. A significant charge is measured at low to moderate monovalent counterion concentrations, resulting in strong inter-DNA repulsion. The charge and repulsion diminish rapidly upon the addition of divalent counterions. An intriguing short range attraction is observed at surprisingly low divalent cation concentrations, approximately 16 mM Mg2+. Quantitative measurements of inter-DNA potentials are essential for improving models of fundamental interactions in biological systems.

Ion-specific repulsive interactions in colloidal silica dispersions evidenced through osmotic compression measurements and implication in frontal ultrafiltration experiments

Comparative osmotic compression experiments were performed on colloidal silica dispersions in the presence of various chloride salts at the same 0.01 M concentration with different counterions and highlighted the influence of ionic specificity on the resistance to water removal. These results were complemented with frontal ultrafiltration measurements that demonstrate modulation of the permeate flux according to the salt used.

The liquidlike ordering of lipid A-diphosphate colloidal crystals: The influence of Ca2+, Mg2+, Na+, and K+ on the ordering of colloidal suspensions of lipid A-diphosphate in aqueous solutions

A comprehensive study was performed on electrostatically stabilized aqueous dispersion of lipid A-diphosphate in the presence of bound Ca2+, Mg2+, K+, and Na+ ions at low ionic strength (0.10-10.0-mM NaCl, 25 degrees C) over a range of volume fraction of 1.0 x 10(-4)< or =phi< or =4.95 x 10(-4). These suspensions were characterized by light scattering (LS), quasielastic light scattering, small-angle x-ray scattering, transmission electron microscopy, scanning electron microscopy, conductivity measurements, and acid-base titrations. LS and electron microscopy yielded similar values for particle sizes, particle size distributions, and polydispersity. The measured static structure factor, S(Q), of lipid A-diphosphate was seen to be heavily dependent on the nature and concentration of the counterions, e.g., Ca2+ at 5.0 nM, Mg2+ at 15.0 microM, and K+ at 100.0 microM (25 degrees C). The magnitude and position of the S(Q) peaks depend not only on the divalent ion concentration (Ca2+ and Mg2+) but also on the order of addition of the counterions to the lipid A-diphosphate suspension in the presence of 0.1-microM NaCl. Significant changes in the rms radii of gyration (R2G) 1/2 of the lipid A-diphosphate particles were observed in the presence of Ca2+ (24.8+/-0.8 nm), Mg2+ (28.5+/-0.7 nm), and K+ (25.2+/-0.6 nm), whereas the Na+ salt (29.1+/-0.8 nm) has a value similar to the one found for the de-ionized lipid A-diphosphate suspensions (29.2+/-0.8 nm). Effective particle charges were determined by fits of the integral equation calculations of the polydisperse static structure factor, S(Q), to the light-scattering data and they were found to be in the range of Z*=700-750 for the lipid A-diphosphate salts under investigation. The light-scattering data indicated that only a small fraction of the ionizable surface sites (phosphate) of the lipid A-diphosphate was partly dissociated (approximately 30%). It was also discovered that a given amount of Ca2+ (1.0-5.0 nM) or K+ (100 microM) influenced the structure much more than Na+ (0.1-10.0-mM NaCl) or Mg2+ (50 microM). By comparing the heights and positions of the structure factor peaks S(Q) for lipid A-diphosphate-Na+ and lipid A-diphosphate-Ca2+, it was concluded that the structure factor does not depend simply on ionic strength but more importantly on the internal structural arrangements of the lipid A-diphosphate assembly in the presence of the bound cations. The liquidlike interactions revealed a considerable degree of ordering in solution accounting for the primary S(Q) peak and also the secondary minimum at large particle separation. The ordering of lipid A-diphosphate-Ca2+ colloidal crystals in suspension showed six to seven discrete diffraction peaks and revealed a face-centered-cubic (fcc) lattice type (a=56.3 nm) at a volume fraction of 3.2 x 10(-4)< or =phi< or =3.9 x 10(-4). The K+ salt also exhibited a fcc lattice (a=55.92 nm) at the same volume fractions, but reveals a different peak intensity distribution, as seen for the lipid A-diphosphate-Ca2+ salt. However, the Mg2+ and the Na+ salts of lipid A-diphosphate showed body-centered-cubic (bcc) lattices with a=45.50 nm and a=41.50 nm, respectively (3.2 x 10(-4)< or =phi< or =3.9 x 10(-4)), displaying the same intensity distribution with the exception of the (220) diffraction peaks, which differ in intensity for both salts of lipid A-diphosphate.

Influence of Confinement on the Electrostatic Interaction between Charged Colloids: a (N,V,T) Monte Carlo Study within Hyperspherical Geometry

Monte Carlo simulations within closed hyperspherical geometry are used to analyze the ionic distribution around two confined charged colloids to determine the origin of the net attraction recently reported in the literature. A scaling procedure is used to compare our numerical results obtained with small ideal colloids with the conclusion of the measurements performed with large silica colloids. Although no electrostatic attraction is detected under confinement, our simulations exhibit a significant reduction of the electrostatic repulsion between charged colloids confined between two weakly charged walls. After rescaling to reproduce the electrostatic repulsion between large confined colloids, our numerical results are qualitatively consistent with the reported attraction because we reasonably expect a reduction of the electrostatic force between such confined colloids below the order of magnitude of their van der Waals attraction.

Effective charges of colloidal particles obtained from collective diffusion experiments

In this work, the collective diffusion coefficient of highly charged colloidal particles in dilute dispersions has been measured by means of dynamic light scattering. The possibility of obtaining valuable information about the particle charge from these data is looked into with the help of electrophoresis experiments. Our results suggest that this is possible in the case of slight or moderately interacting particles as long as experimental data are properly treated. For highly interacting colloids, however, such information could not be so reliable, presumably due to certain shortcomings of the experimental technique at low angle. The role of charge renormalization is also discussed in this work.

Overcharging in colloids: beyond the Poisson-Boltzmann approach

A broad range of manufactured products and biological fluids are colliods. The ability to understand and control the processes (of scientific, technological and industrial interest) in which such colloids are involved relies upon a precise knowledge of the electrical double layer. The traditional approach to describing this ion cloud around colloidal particles has been the Gouy-Chapman model developed on the basis of the Poisson-Boltzmann equation. Since the early 1980s, however, more sophisticated theoretical treatments have revealed both quantitative and qualitative deficiencies in the Poisson-Boltzmann theory, particularly at high ionic strengths and/or high surface charge densities. This review deals with these novel approaches, which are mostly computer simulations and approximate integral equation theories based on the so-called primitive model. Special attention is paid to phenomena that cannot be accounted for by the classic theory as a result of neglecting ion size correlations, such as overcharging, namely, the counterion concentration in the immediate neighborhood of the surface is so large that the particle surface is overcompensated. Other illustrative examples are the nonmonotonic behavior of the electrostatic potential and attractive interactions between equally charged surfaces. These predictions are certainly remarkable and, on paper, they can have an effect on experimentally measurable quantities (for instance, electrophoretic mobility). Even so, these new approaches have scarcely been applied in practice. Thus a critical survey on the relevance of ion size correlation in real systems is also included. Overcharging of macroions can also be brought about by adsorption of oppositely charged polyelectrolytes. Noteworthy examples and theoretical approaches for them are also briefly reviewed.