Revisiting the birth of NaCl crystals using molecular dynamics simulation.

Molecular dynamics simulations were used to investigate the initial stage of phase separation mechanisms for an oversaturated electrolytic solution. We developed a low computational cost methodology to determine the simulation frames where the first ionic clusters are formed. By discretizing the simulation box, we obtain a density profile in the moments preceding and succeeding the nuclei's formation. The growth of the clusters identified with our methodology was analyzed until the end of the simulation. Calculation of the Steinhardt parameter showed symmetry of the solid, giving indications that the classical nucleation theory explains the mechanism of the solid formation. The methodology developed was useful for identifying phase separation mechanisms in the nucleation process. At lower concentrations, there was no formation of stable clusters. At intermediate concentrations, the analyses indicate a transition of phases in one stage, from a oversaturate electrolytic solution to a crystalline solid. At high concentration, a transition of phases in two stages, initially, is the formation of a dense liquid, and only after that, crystalline solid formed inside the dense liquid. The change in phase separation mechanism due to increasing oversaturation underscores the importance of precise determination of the driving force for phase separation and concentration limits for each mechanism.

Transitional phase inversion of emulsions monitored by in situ near-infrared spectroscopy.

Water-heptane/toluene model emulsions were prepared to study emulsion transitional phase inversion by in situ near-infrared spectroscopy (NIR). The first emulsion contained a small amount of ionic surfactant (0.27 wt % of sodium dodecyl sulfate) and n-pentanol as a cosurfactant. In this emulsion, the study was guided by an inversion coordinate route based on a phase behavior study previously performed. The morphology changes were induced by rising aqueous phase salinity in a "steady-state" inversion protocol. The second emulsion contained a nonionic surfactant (ethoxylated nonylphenol) at a concentration of 3 wt %. A continuous temperature change induced two distinct transitional phase inversions: one occurred during the heating of the system and another during the cooling. NIR spectroscopy was able to detect phase inversion in these emulsions due to differences between light scattered/absorbed by water in oil (W/O) and oil in water (O/W) morphologies. It was observed that the two model emulsions exhibit different inversion mechanisms closely related to different quantities of the middle phases formed during the three-phase behavior of Winsor type III.

Anion-specific partitioning in two-phase finite volume systems: possible implications for mechanisms of ion pumps.

In two-phase finite volume systems of electroneutral phospholipids, the electrolyte concentration is different in the two phases. The partitioning is highly anion-specific, a phenomenon not accounted for by classical electrolyte theories. It is explained if ionic dispersion forces that lead to specific ion binding are taken into account. The mechanism provides a contribution to active ion pumps not previously considered.

The influence of ion binding and ion specific potentials on the double layer pressure between charged bilayers at low salt concentrations.

Measurements of surface forces between double-chained cationic bilayers adsorbed onto molecularly smooth mica surfaces across different millimolar salt solutions have revealed a large degree of ion specificity [Pashley et al., J. Phys. Chem. 90, 1637 (1986)]. This has been interpreted in terms of highly specific anion binding to the adsorbed bilayers. We show here that inclusion in the double layer theory of nonspecific ion binding and ion specific nonelectrostatic potentials acting between ions and the two surfaces can account for the phenomenon. It also gives the right Hofmeister series for the double layer pressure.

Co-ion and ion competition effects: ion distributions close to a hydrophobic solid surface in mixed electrolyte solutions.

We consider within a modified Poisson-Boltzmann theory an electrolyte, with different mixtures of NaCl and NaI, near uncharged and charged solid hydrophobic surfaces. The parametrized potentials of mean force acting on Na+, Cl-, and I- near an uncharged self-assembled monolayer were deduced from molecular simulations with polarizable force fields. We study what happens when the surface presents negative charges. At moderately charged surfaces, we observe strong co-ion adsorption and clear specific ion effects at biological concentrations. At high surface charge densities, the co-ions are pushed away from the interface. We predict that Cl- ions can also be excluded from the surface by increasing the concentration of NaI. This ion competition effect (I- versus Cl-) may be relevant for ion-specific partitioning in multiphase systems where polarizable ions accumulate in phases with large surface areas.

Specific ion adsorption and surface forces in colloid science.

Mean-field theories that include nonelectrostatic interactions acting on ions near interfaces have been found to accommodate many experimentally observed ion specific effects. However, it is clear that this approach does not fully account for the liquid molecular structure and hydration effects. This is now improved by using parametrized ionic potentials deduced from recent nonprimitive model molecular dynamics (MD) simulations in a generalized Poisson-Boltzmann equation. We investigate how ion distributions and double layer forces depend on the choice of background salt. There is a strong ion specific double layer force set up due to unequal ion specific short-range potentials acting between ions and surfaces.

Ion-specific forces between a colloidal nanoprobe and a charged surface.

We investigate the effect of ion-specific potentials on the force between a nanoprobe attached to a cantilever tip, and a charged surface. The probe is treated as a spherical nanoparticle with constant charge. A modified Poisson-Boltzmann equation in bispherical coordinates is used to address this problem in a more quantitative way. We predict that the ion-specific series of measured forces depend on the sign and magnitude of surface charge densities.

Specific ion effects in solutions of globular proteins: comparison between analytical models and simulation.

Monte Carlo simulations have been performed for ion distributions outside a single globular macroion and for a pair of macroions, in different salt solutions. The model that we use includes both electrostatic and van der Waals interactions between ions and between ions and macroions. Simulation results are compared with the predictions of the Ornstein-Zernike equation with the hypernetted chain closure approximation and the nonlinear Poisson-Boltzmann equation, both augmented by pertinent van der Waals terms. Ion distributions from analytical approximations are generally very close to the simulation results. This demonstrates that properties that are related to ion distributions in the double layer outside a single interface can to a good approximation be obtained from the Poisson-Boltzmann equation. We also present simulation and integral equation results for the mean force between two globular macroions (with properties corresponding to those of hen-egg-white lysozyme protein at pH 4.3) in different salt solutions. The mean force and potential of mean force between the macroions become more attractive upon increasing the polarizability of the counterions (anions), in qualitative agreement with experiments. We finally show that the deduced second virial coefficients agree quite well with experimental results.

Why forces between proteins follow different Hofmeister series for pH above and below pI.

The relative effectiveness of different anions in crystallizing proteins follows a reversed Hofmeister sequence for pHpI. The phenomenon has been known almost since Hofmeister's original work but it has not been understood. It is here given a theoretical explanation. Classical electrolyte and double layer theory deals only with electrostatic forces acting between ions and proteins. Hydration and hydration interactions are dealt with usually only in terms of assumed hard core models. But there are, at and above biological salt concentrations, other non-electrostatic (NES) ion-specific forces acting that are ignored in such modeling. Such electrodynamic fluctuation forces are also responsible for ion-specific hydration. These missing forces are variously comprehended under familiar but generally unquantified terms, typically, hydration, hydrogen bonding, pi-electron-cation interactions, dipole-dipole, dipole-induced dipole and induced dipole-induced dipole forces and so on. The many important body electrodynamic fluctuation force contributions are accessible from extensions of Lifshitz theory from which, with relevant dielectric susceptibility data on solutions as a function of frequency, the forces can be extracted quantitatively, at least in principle. The classical theories of colloid science that miss such contributions do not account for a whole variety of ion-specific phenomena. Numerical results that include these non-electrostatic forces are given here for model calculations of the force between two model charge-regulated hen-egg-white protein surfaces. The surfaces are chosen to carry the same charge groups and charge density as the protein. What emerges is that for pHpI (where anions are co-ions) the forces increase in the order NaCl

Centrifugation equilibrium for spheres and spherocylinders.

The centrifugation equilibrium problem is formulated and solved using a new procedure in which the specified variables are the temperature, system volume, particle dimensions and concentrations, angular speed, cell length, and cell distance from the rotation axis. As a result, we obtain the concentration profiles for all types of particles present in the system, which are considered to be immersed in a fluid. The particles are modeled as hard nonattractive spherocylinders using an equation of state, but the procedure is not restricted to any geometrical shape, and can be used with any equation of state available. The fluid is treated as a continuous medium, responsible for centrifugal buoyancy. We make calculations for colloidal suspensions of silica, often used for separations in biotechnology. Results are in good agreement with experiments and show excellent agreement in comparison with Monte Carlo simulations. Our calculations also predict focusing and shifting phenomena that have been experimentally observed in separations of fine particles.

Phase behavior of aqueous solutions containing dipolar proteins from second-order perturbation theory.

Due to the interplay of Coulombic repulsion and attractive dipolar and van der Waals interactions, solutions of globular proteins display a rich variety of phase behavior featuring fluid-fluid and fluid-solid transitions that strongly depend on solution pH and salt concentration. Using a simple model for charge, dispersion and dipole-related contributions to the interprotein potential, we calculate phase diagrams for protein solutions within the framework of second-order perturbation theory. For each phase, we determine the Helmholtz energy as the sum of a hard-sphere reference term and a perturbation term that reflects both the electrostatic and dispersion interactions. Dipolar effects can induce fluid-fluid phase separation or crystallization even in the absence of any significant dispersion attraction. Because dissolved electrolytes screen the charge-charge repulsion more strongly than the dipolar attraction, the ionic strength dependence of the potential of mean force can feature a minimum at intermediate ionic strengths offering an explanation for the observed nonmonotonic dependence of the phase behavior on salt concentration. Inclusion of correlations between charge-dipole and dipole-dipole interactions is essential for a reliable calculation of phase diagrams for systems containing charged dipolar proteins and colloids.