Coarse-grained molecular dynamics simulations of the sphere to rod transition in surfactant micelles

Surfactant molecules self-assemble in aqueous solutions to form various micellar structures such as spheres, rods, or lamellae. Although phase transitions in surfactant solutions have been studied experimentally, their molecular mechanisms are still not well understood. In this work, we show that molecular dynamics (MD) simulations using the coarse-grained (CG) MARTINI force field and explicit CG solvent, validated against atomistic MD studies, can accurately represent micellar assemblies of cetyltrimethylammonium chloride (CTAC). The effect of salt on micellar structures is studied for aromatic anionic salts, e.g., sodium salicylate, and simple inorganic salts, e.g., sodium chloride. Above a threshold concentration, sodium salicylate induces a sphere to rod transition in the micelle. CG MD simulations are shown to capture the dynamics of this shape transition and support a mechanism based on the reduction in the micelle-water interfacial tension induced by the adsorption of the amphiphilic salicylate ions. At the threshold salt concentration, the interface is nearly saturated with adsorbed salicylate ions. Predictions of the effect of salt on the micelle structure in different CG solvent models, namely, single-site standard water and three-site polarizable water, show qualitative agreement. This suggests that phase transitions in aqueous micelle solutions could be investigated by using standard CG water models which allow for 3 orders of magnitude reduction in the computational time as compared to that required for atomistic MD simulations.     

Surface segregation of dissolved salt ions

Surface segregation of iodide, but not of fluoride or cesium ions, is observed by a combination of metastable impact electron spectroscopy (MIES) and ultraviolet photoelectron spectroscopy (UPS(HeI)) of amorphous solid water exposed to CsI or CsF vapor. The same surface ionic behavior is also derived from molecular dynamics (MD) simulations of the corresponding aqueous salt solutions. The MIES results show the propensity of iodide, but not fluoride, for the surface of the amorphous solid water film, providing thus strong evidence for the suggested presence of heavier halides (iodide, bromide, and to a lesser extent chloride) at the topmost layer of aqueous surfaces. In contrast, no appreciable surface segregation of ions is observed in methanol, neither in the experiment nor in the simulation. Furthermore, the present results indicate that, as far as the thermodynamic aspects of solvation of alkali halides are concerned, amorphous solid water and methanol surfaces behave similarly as surfaces of the corresponding liquids.     

Air-liquid interfaces of aqueous solutions containing ammonium and sulfate: spectroscopic and molecular dynamics studies

Investigations of the air-liquid interface of aqueous salt solutions containing ammonium (NH(4)(+)) and sulfate (SO(4)(2-)) ions were carried out using molecular dynamics simulations and vibrational sum frequency generation spectroscopy. The molecular dynamics simulations show that the predominant effect of SO(4)(2-) ions, which are strongly repelled from the surface, is to increase the thickness of the interfacial region. The vibrational spectra reported are in the O-H stretching region of liquid water. Isotropic Raman and ATR-FTIR (attenuated total reflection Fourier transform infrared) spectroscopies were used to study the effect of ammonium and sulfate ions on the bulk structure of water, whereas surface sum frequency generation spectroscopy was used to study the effect of these ions on the interfacial structure of water. Analysis of the interfacial and bulk vibrational spectra reveal that aqueous solutions containing SO(4)(2-) perturb the interfacial water structure differently than the bulk and, consistent with the molecular dynamics simulations, reveal an increase in the thickness of the interfacial region.     

Autoionization at the surface of neat water: is the top layer pH neutral, basic, or acidic?

Autoionization of water which gives rise to its pH is one of the key properties of aqueous systems. Surfaces of water and aqueous electrolyte solutions are traditionally viewed as devoid of inorganic ions; however, recent molecular simulations and spectroscopic experiments show the presence of certain ions including hydronium in the topmost layer. This raises the question of what is the pH (defined using proton concentration in the topmost layer) of the surface of neat water. Microscopic simulations and measurements with atomistic resolution show that the water surface is acidic due to a strong propensity of hydronium (but not of hydroxide) for the surface. In contrast, macroscopic experiments, such as zeta potential and titration measurements, indicate a negatively charged water surface interpreted in terms of preferential adsorption of OH(-). Here we review recent simulations and experiments characterizing autoionization at the surface of liquid water and ice crystals in an attempt to present and discuss in detail, if not fully resolve, this controversy.     

Nanophase segregation and water dynamics in the dendrion diblock copolymer formed from the Fréchet polyaryl ethereal dendrimer and linear PTFE

We propose a new material consisting of a dendrion copolymer formed from (a) a water-soluble dendritic polymer and (b) a hydrophobic backbone. Using molecular dynamics simulations techniques, we determine the structure and dynamics of the dendrion formed by second-generation Fréchet polyaryl ethereal dendrimer as the hydrophilic component and linear polytetrafluoroethylene (PTFE) as the hydrophobic polymer, with 5 and 10 wt % of water. We find that this material produces a well-developed nanoscale structure in which water forms a continuous nanophase, making this new family of compounds promising candidates for applications in fuel cell membranes. We find that the water molecules are incorporated into the dendrimer block of the copolymer to form a nanophase-segregated structure. The well-developed nanophase-segregated structures rendered by this material have characteristic dimensions of segregation ( approximately 30 Angstrom) and dendrimer conformational properties that are independent of water content. Calculations of water dynamics and proton transport in these nanophase-segregated structures indicate that the dendrion copolymer membrane with 10 wt % of water content has a water structure and transport properties equivalent to that of the hydrated Nafion membrane with 20 wt % of water content.     

Kinetics of oxygen absorption by aqueous electrolyte solutions in the presence of microencapsulated quartz particles activating the mass transfer in the liquid phase

Effect of the ion composition of aqueous solutions on the oxygen absorption kinetics in a system constituted by a gas (air) and a liquid (aqueous solution) in the presence of microencapsulated quartz particles activating the mass transfer in the liquid phase was studied. It was found that ions with positive hydration cause a substantial decrease in the O₂ mass-transfer enhancement factor, whereas ions with negative hydration lead to its increase under the same conditions. It is shown that the effect of ions on the rate of oxygen absorption by aqueous electrolyte solutions can be prognosticated on the basis of data on the influence of these ions on the structure and viscosity of water. The results of the study can serve as a basis for varying the rate of heterogeneous reactions in gas-liquid systems, whose rate is limited by the mass transfer of oxygen into aqueous media, by purposeful control over their ion composition.     

Structural and dynamic properties of concentrated alkali halide solutions: a molecular dynamics simulation study

The physicochemical properties of alkali halide solutions have long been attributed to the collective interactions between ions and water molecules in the solution, yet the structure of water in these systems and its effect on the equilibrium and dynamic properties of these systems are not clearly understood. Here, we present a systematic view of water structure in concentrated alkali halide solutions using molecular dynamics simulations. The results of the simulations show that the size of univalent ions in the solution has a significant effect on the dynamics of ions and other transport properties such as the viscosity that are correlated with the structural properties of water in aqueous ionic solution. Small cations (e.g., Li+) form electrostatically stabilized hydrophilic hydration shells that are different from the hydration shells of large ions (e.g., Cs+) which behave more like neutral hydrophobic particles, encapsulated by hydrogen-bonded hydration cages. The properties of solutions with different types of ion solvation change in different ways as the ion concentration increases. Examples of this are the diffusion coefficients of the ions and the viscosities of solutions. In this paper we use molecular dynamics (MD) simulations to study the changes in the equilibrium and transport properties of LiCl, RbCl, and CsI solutions at concentrations from 0.22 to 3.97 M.     

Proton transfer and the mobilities of the H+ and OH- ions from studies of a dissociating model for water

Hydrogen (H(+)) and hydroxide (OH(-)) ions in aqueous solution have anomalously large diffusion coefficients, and the mobility of the H(+) ion is nearly twice that of the OH(-) ion. We describe molecular dynamics simulations of a dissociating model for liquid water based on scaling the interatomic potential for water developed by Ojama?e-Shavitt-Singer from ab initio studies at the MP2 level. We use the scaled model to study proton transfer that occurs in the transport of hydrogen and hydroxide ions in acidic and basic solutions containing 215 water molecules. The model supports the Eigen-Zundel-Eigen mechanism of proton transfer in acidic solutions and the transient hyper-coordination of the hydroxide ion in weakly basic solutions at room temperature. The free energy barriers for proton transport are low indicating significant proton delocalization accompanying proton transfer in acidic and basic solutions. The reorientation dynamics of the hydroxide ion suggests changes in the proportions of hyper-coordinated species with temperature. The mobilities of the hydrogen and hydroxide ions and their temperature dependence between 0 and 50 °C are in excellent agreement with experiment and the reasons for the large difference in the mobilities of the two ions are discussed. The model and methods described provide a novel approach to studies of liquid water, proton transfer, and acid-base reactions in aqueous solutions, channels, and interfaces.     

Insights into Tris‐(2‐Hydroxylethyl)methylammonium Methylsulfate Aqueous Solutions

We report herein a combined experimental–computational study on tris‐(2‐hydroxylethyl)methylammonium methylsulfate in water solutions, as a representative ionic liquid of the aqueous‐solution behavior of hydroxylammonium‐based ionic liquids. Relevant thermophysical properties were measured as a function of mixture composition and temperature. Classical molecular dynamics simulations were performed to infer microscopic structural features. The reported results for ionic liquid in water‐rich solutions show that it behaves as isolated non‐interacting ions solvated by water molecules, through well‐defined solvation shells, exerting a disrupting effect on the water hydrogen bonding network. Nevertheless, as ionic liquid concentration increase, interionic association increases, even for diluted water solutions, evolving from the typical behavior of strong electrolytes in solution toward large interacting structures. For ionic‐liquid‐rich mixtures, water exerts a minor disrupting effect on the fluid’s structuring because it occupies regions around each ion (developing water–ion hydrogen bonds) but without significantly weakening anion–cation interactions.     

Electron binding energies of aqueous alkali and halide ions: EUV photoelectron spectroscopy of liquid solutions and combined ab initio and molecular dynamics calculations

Photoelectron spectroscopy combined with the liquid microjet technique enables the direct probing of the electronic structure of aqueous solutions. We report measured and calculated lowest vertical electron binding energies of aqueous alkali cations and halide anions. In some cases, ejection from deeper electronic levels of the solute could be observed. Electron binding energies of a given aqueous ion are found to be independent of the counterion and the salt concentration. The experimental results are complemented by ab initio calculations, at the MP2 and CCSD(T) level, of the ionization energies of these prototype ions in the aqueous phase. The solvent effect was accounted for in the electronic structure calculations in two ways. An explicit inclusion of discrete water molecules using a set of snapshots from an equilibrium classical molecular dynamics simulations and a fractional charge representation of solvent molecules give good results for halide ions. The electron binding energies of alkali cations computed with this approach tend to be overestimated. On the other hand, the polarizable continuum model, which strictly provides adiabatic binding energies, performs well for the alkali cations but fails for the halides. Photon energies in the experiment were in the EUV region (typically 100 eV) for which the technique is probing the top layers of the liquid sample. Hence, the reported energies of aqueous ions are closely connected with both structures and chemical reactivity at the liquid interface, for example, in atmospheric aerosol particles, as well as fundamental bulk solvation properties.     

Molecular dynamics simulations of aqueous solutions of ethanolamines

We report on molecular dynamics simulations performed at constant temperature and pressure to study ethanolamines as pure components and in aqueous solutions. A new geometric integration algorithm that preserves the correct phase space volume is employed to study molecules having up to three ethanol chains. The most stable geometry, rotational barriers, and atomic charges were obtained by ab initio calculations in the gas phase. The calculated dipole moments agree well with available experimental data. The most stable conformation, due to intramolecular hydrogen bonding interactions, has a ringlike structure in one of the ethanol chains, leading to high molecular stability. All molecular dynamics simulations were performed in the liquid phase. The interaction parameters are the same for the atoms in the ethanol chains, reducing the number of variables in the potential model. Intermolecular hydrogen bonding is also analyzed, and it is shown that water associates at low water mole fractions. The force field reproduced (within 1%) the experimental liquid densities at different temperatures of pure components and aqueous solutions at 313 K. The excess and partial molar volumes are analyzed as a function of ethanolamine concentration.     

On the atomistic mechanisms of alkane (methane-pentane) separation by distillation: a molecular dynamics study

Insights into the liquid-vapor transformation of methane-pentane mixtures were obtained from transition path sampling molecular dynamics simulations. This case study of the boiling of non-azeotropic mixtures demonstrates an unprejudiced identification of the atomistic mechanisms of phase separation in the course of vaporization which form the basis of distillation processes. From our simulations we observe spontaneous segregation events in the liquid mixture to trigger vapor nucleation. The formation of vapor domains stabilizes and further promotes the separation process by preferential evaporation of methane molecules. While this discrimination holds for all mixtures of different composition studied, a full account of the boiling process requires a more complex picture. At low methane concentration the nucleation of the vapor domains includes both methane and pentane molecules. The pentane molecules, however, tend to form small aggregates and undergo rapid re-condensation within picoseconds to nanoseconds scales. Accordingly, two aspects of selectivity accounting for methane-pentane separation in the course of liquid-vapor transformations were made accessible to molecular dynamics simulations: spontaneous segregation in the liquid phase leading to selective vapor nucleation and growth favoring methane vaporization and selective re-condensation of pentane molecules.     

Phase equilibria and critical phenomena in the cesium nitrate–water–pyridine ternary system

The solubilities of components, phase equilibria, and critical phenomena in the cesium nitrate–water–pyridine ternary system are studied in the 5–100°C temperature range by the visual–polythermal method. Cesium nitrate is found to exhibit a salting-out effect at temperatures above 79.9°C causing phase separation in homogeneous water–pyridine solutions. The temperature of formation of the critical monotectic tie line (79.9°C) and the compositions of solutions corresponding to the liquid–liquid critical points at three temperatures are determined. The pyridine distribution coefficients between the aqueous and organic phases of the monotectic state at 85.0, 90.0, and 100.0°C are calculated. Their values demonstrate that salting-out of pyridine from aqueous solutions by cesium nitrate increases at higher temperatures. The plotted isotherms of phase diagrams confirm the fragment of the scheme of topological transformation of the phase diagrams of salt–binary solvent ternary systems with salting-in and salting-out phenomena.     

The structures of ozone and HOx radicals in aqueous solution from combined quantum/classical molecular dynamics simulations

Ozone in aqueous solution decomposes through a complex mechanism that involves initial reaction with a hydroxide ion followed by formation of a variety of oxidizing species such as HO, HO(2), and HO(3) radicals. Though a number of hydrogen-bonded complexes have been described in the gas phase, both theoretically and experimentally, the structures of ozone and HO(x) in liquid water remain uncertain. In this work, combined quantum/classical computer simulations of aqueous solutions of these species have been reported. The results show that ozone undergoes noticeable electron polarization but it does not participate in hydrogen bonds with liquid water. The main contribution of the solvation energy comes from dispersion forces. In contrast, HO(x) radicals form strong hydrogen bonds. They are better proton donors but weaker proton acceptors than water. Their electronic and geometrical structures are significantly modified by the solvent, especially in the case of HO(3). In all cases, fluctuations in amplitudes of electronic properties are considerable, suggesting that solvent effects might play a crucial role on oxidation mechanisms initiated by ozone in liquid water. These mechanisms are important in a broad range of domains, such as atmospheric processes, plant response to ambient ozone, and medical and industrial applications.     

The [BMI][Tf2N] ionic liquid/water binary system: a molecular dynamics study of phase separation and of the liquid-liquid interface

We report molecular dynamics (MD) simulations of the aqueous interface of the hydrophobic [BMI][Tf2N] ionic liquid (IL), composed of 1-butyl-3-methylimidazolium cations (BMI+) and bis(trifluoromethylsulfonyl)imide anions (Tf2N-). The questions of water/IL phase separation and properties of the neat interface are addressed, comparing different liquid models (TIP3P vs TIP5P water and +1.0/-1.0 vs +0.9/-0.9 charged IL ions), the Ewald vs the reaction field treatments of the long range electrostatics, and different starting conditions. With the different models, the "randomly" mixed liquids separate much more slowly (in 20 to 40 ns) than classical water-oil mixtures do (typically, in less than 1 ns), finally leading to distinct nanoscopic phases separated by an interface, as in simulations which started with a preformed interface, but the IL phase is more humid. The final state of water in the IL thus depends on the protocol and relates to IL heterogeneities and viscosity. Water mainly fluctuates in hydrophilic basins (rich in O(Tf2N) and aromatic CH(BMI) groups), separated by more hydrophobic domains (rich in CF3(Tf2N) and alkyl(BMI) groups), in the form of monomers and dimers in the weakly humid IL phase, and as higher aggregates when the IL phase is more humid. There is more water in the IL than IL in water, to different extents, depending on the model. The interface is sharper and narrower (approximately 10 A) than with the less hydrophobic [BMI][PF6] IL and is overall neutral, with isotropically oriented molecules, as in the bulk phases. The results allow us to better understand the analogies and differences of aqueous interfaces with hydrophobic (but hygroscopic) ILs, compared to classical organic liquids.     

Photochemical Method of Controlling Dispersity of Nanostructures of Transition Metals

Original methods of photochemical preparation of stable mono- (Cu, Ag, Au, Ni, Pd, Pt) and bimetallic nanoparticles in the form of optically transparent compact films on quartz surface and of volume dispersions in porous inorganic (silica glasses) and organic (MF-4SK fluorocarbon) materials, solid polymers (polyvinyl alcohol, polyethylene glycol, gelatin, latexes), and in liquid media (glycerol) are reviewed. The results of studies of spectral and structural characteristics of nanophase films under various experimental conditions are presented. Experimental mechanistic models and ways of controlling disperse composition of metal colloids, which can be used for photochemical synthesis of nanophase systems, are proposed.     

Structural and Electrosurface Characteristics of Sulfonated Cation-Exchange Perfluoropolymer Membranes in 1 : 1 Electrolyte Solutions

The equilibrium (the exchange capacity, the structural resistance coefficient, and contact angles) and transport (the conductivity) characteristics of differently obtained sulfonated cation-exchange perfluoropolymeric membranes in 1 : 1 electrolyte solutions were investigated. It was shown that the transformation of membranes from the Na⁺form to the K⁺form sharply decreases their moisture content, which is accompanied by an increase in the structural resistance coefficient and the counterion concentration in membranes. Experimental data were used for calculating the electrochemical characteristics of membranes: efficiency coefficients, ion mobilities and transfer numbers of ions in the intramembrane liquid, as well as Donnan potentials. Measurements of the wettability of fluoroplastic ion-exchange membranes with water and electrolyte solutions showed that the presence of strongly acidic ionogenic groups significantly decreases contact angles as compared with that of the polytetrafluoroethylene surface. It was also established that, for the investigated systems, the contact angle is virtually independent of the composition of the liquid phase.     

Aqueous solutions of divalent chlorides: ions hydration shell and water structure

By combining neutron diffraction and Monte Carlo simulations, we have determined the microscopic structure of the hydration ions shell in aqueous solutions of MgCl(2) and CaCl(2), along with the radial distribution functions of the solvent. In particular the hydration shell of the cations, show cation specific symmetry, due to the strong and directional interaction of ions and water oxygens. The ions and their hydration shells likely form molecular moieties and bring clear signatures in the water-water radial distribution functions. Apart from these signatures, the influence of divalent salts on the microscopic structure of water is similar to that of previously investigated monovalent solutes, and it is visible as a shift of the second peak of the oxygen-oxygen radial distribution function, caused by distortion of the hydrogen bond network of water.     

Quantitative characterization of ion pairing and cluster formation in strong 1:1 electrolytes

Aqueous solutions of 1:1 strong electrolytes are considered to be the prototype for complete ionic dissociation. Nonetheless, clustering of strong 1:1 electrolytes has been widely reported in all atom molecular dynamics simulations, and their presence is indirectly implicated in a diverse range of experimental results. Is there a physical basis for nonidealities such as ion pairing and cluster formation in aqueous solutions of strong 1:1 electrolytes? We attempt to answer this question by direct comparison of results from detailed molecular dynamics simulations to experimentally observed properties of 1:1 electrolytes. We report the analysis of a series of lengthy molecular dynamics simulations of alkali-halide solutions carried out over a wide range of physiologically relevant concentrations using explicit representations of water molecules. We find evidence for pronounced nonideal behavior of ions at all concentrations in the form of ion pairs and clusters which are in rapid equilibrium with dissociated ions. The phenomenology for ion pairing seen in these simulations is congruent with the multistep scheme proposed by Eigen and Tamm based on data from ultrasonic absorption experiments. For a given electrolyte, we show that the dependence of cluster populations on concentration can be described through a single set of equilibrium constants. We assess the accuracy of calculated ion pairing constants by favorable comparison to estimates obtained by Fuoss and co-workers and based on conductometric experiments. Ion pairs and clusters form on length scales where the size of individual water molecules is as important as the hard core radius of ions. Ion pairing results as a balance between the favorable Coulomb interactions and the unfavorable partial desolvation of ions needed to form a pair.