Hofmeister effects in surface tension of aqueous electrolyte solution

The surface tension of electrolyte solutions shows marked specific ion effects. We here show an important role for both ionic solvation energies and ionic dispersion potentials in determining this ion specific surface tension of salt solutions. The ion self-free energy changes when an ion moves from bulk solution into the interfacial region, with its decreasing water density profile. We will show that the solvation energies of different ions correlate very well with the surface tension of salt solutions. Inclusion of this distance-dependent self-free energy contribution brings qualitative agreement with experiments and the right Hofmeister series. This is so not only for surface tension changes but also for measured surface potentials. The inclusion of ionic dispersion interaction potentials further improves the agreement with experiments. We discuss how further progress in the theory of the surface tension of salts can be achieved.     

Ion hydration near air/water interfaces and the structure of liquid surfaces

Negative adsorption of ions, commonly observed at air-water interfaces, is examined in terms of models of restricted polarization of the solvent by ions at the interface and the structure of the liquid interface. The Born and other models of ionic hydration are applied to evaluate the self-energy of the ion arising in the region of solvent near its interface and in the vacuum or vapour beyond. The adsorption energy of an ion varies substantially with distance from the liquid interface so that a distribution of ions arises as a function of distance from the interface. Integration of this distribution gives an expression, and results, for the ionic surface excess. The diffuse-layer potential, which an unequal distribution of cations and anions give rise to, gives a contribution to the surface potential of the electrolyte solution at finite concentrations.Structural aspects of the liquid interface at which ions are negatively adsorbed are discussed in terms of Stefan's ratio and the superficial excess entropies of various liquid surfaces. These entropies are related to the cohesive energy densities of the bulk liquids. Ion solvent-structure co-sphere interactions with structured interfaces will lead to specificity of negative adsorption of ions.     

A classical point charge model study of system size dependence of oxidation and reorganization free energies in aqueous solution

The response of water to a change of charge of a solvated ion is, to a good approximation, linear for the type of iron-like ions frequently used as a model system in classical force field studies of electron transfer. Free energies for such systems can be directly calculated from average vertical energy gaps. Exploiting this feature, we have computed the free energy and the reorganization energy of the M2+/M3+ and M1+/M2+ oxidations in a series of model systems all containing a single Mn+ ion and an increasing number of simple point charge water molecules. Long-range interactions are taken into account by Ewald summation methods. Our calculations confirm the observation made by Hummer, Pratt, and Garcia (J. Phys. Chem. 1996, 100, 1206) that the finite size correction to the estimate of solvation energy (and hence oxidation free energy) in such a setup is effectively proportional to the inverse third power (1/L3) of the length L of the periodic cell. The finite size correction to the reorganization energy is found to scale with 1/L. These simulation results are analyzed using a periodic generalization of the Born cavity model for solvation, yielding three different estimates of the cavity radius, namely, from the infinite system size extrapolation of oxidation free energy and reorganization energy, and from the slope of the linear dependence of oxidation free energy on 1/L3. The cavity radius for the reorganization energy is found to be significantly larger compared to the radius for the oxidation (solvation) free energy. The radius controlling the 1/L3 dependence of oxidation free energy is found to be comparable to the radius for reorganization. The implication of these results for density functional theory-based ab initio molecular dynamics calculation of redox potentials is discussed.     

Calculation of electron affinities of polycyclic aromatic hydrocarbons and solvation energies of their radical anion

Electron affinities (EAs) and free energies for electron attachment (DeltaGo(a,298K)) have been directly calculated for 45 polynuclear aromatic hydrocarbons (PAHs) and related molecules by a variety of theoretical methods, with standard regression errors of about 0.07 eV (mean unsigned error = 0.05 eV) at the B3LYP/6-31 + G(d,p) level and larger errors with HF or MP2 methods or using Koopmans' Theorem. Comparison of gas-phase free energies with solution-phase reduction potentials provides a measure of solvation energy differences between the radical anion and neutral PAH. A simple Born-charging model approximates the solvation effects on the radical anions, leading to a good correlation with experimental solvation energy differences. This is used to estimate unknown or questionable EAs from reduction potentials. Two independent methods are used to predict DeltaGo(a,298K) values: (1) based upon DFT methods, or (2) based upon reduction potentials and the Born model. They suggest reassignments or a resolution of conflicting experimental EAs for nearly one-half (17 of 38) of the PAH molecules for which experimental EAs have been reported. For the antiaromatic molecules, 1,3,5-tri-tert-butylpentalene and the dithia-substituted cyclobutadiene 1, the reduction potentials lead to estimated EAs close to those expected from DFT calculations and provide a basis for the prediction of the EAs and reduction potentials of pentalene and cyclobutadiene. The Born model has been used to relate the electrostatic solvation energies of PAH and hydrocarbon radical anions, and spherical halide anions, alkali metal cations, and ammonium ions to effective ionic radii from DFT electron-density envelopes. The Born model used for PAHs has been successfully extended here to quantitatively explain the solvation energy of the C60 radical anion.     

Computer simulation of dissociative equilibrium in aqueous NaCl electrolyte with account for polarization and ion recharging. Model of interactions

Comparative analysis of different models applied in theoretical studies of electrolytes points to the considerable role of polarization interactions. In order to study aqueous electrolytes on the molecular level, a detailed model is presented of intermolecular interactions that accounts in explicit form, apart from Coulombic, exchange, and dispersion interactions, also many-body interactions caused by polarization of solvent molecules in the field of ions and ion polarization in the field of solvent molecules, and also many-body covalent interactions, effects of excess charge transfer and effects of partial counterion recharging. Numeric values of parameters for an aqueous NaCl solution are obtained by correlation of the calculated values of free energies and entropy of the first reactions of vapor molecule addition to hydrate ion shells with the corresponding experimental values and also from the experimental data on the dissociation energy and IR vibration frequencies of the ion pair and quantum-chemical calculation of the energy of stable configurations of a hydrated ion pair. A special form of describing many-body interactions allows by more than an order of magnitudes reducing the extent of the required calculations and makes it possible to apply the developed model for computer simulation of aqueous electrolytes at the room temperatures.     

Real energies of ferricinium cation transfer from water to organo-aqueous solvents

The redox potentials of ferrocene-ferricinium couple are determined for five organo-aqueous systems with respect to a reversible chloride electrode in the same solvent. Combined with the literature data on real (and some chemical) transfer energies of chloride ion, the effective real energies of ferricinium transfer from water into organo-aqueous mixtures are determined. These values do not obey the Born equation. Most probably, this is due to the selective solvation of cation by water and of neutral ferrocene molecule by an organic component.     

Free energy of solvation of simple ions: molecular-dynamics study of solvation of Cl- and Na+ in the ice/water interface

Molecular-dynamics simulations of Cl(-) and Na(+) ions are performed to calculate ionic solvation free energies in both bulk simple point-charge/extended water and ice 1 h at several different temperatures, and at the basal ice 1 h/water interface. For the interface we calculate the free energy of "transfer" of the ions across the ice/water interface. For the ions in bulk water in the NPT ensemble at 298 K and 1 atm, results are found to be in good agreement with experiments, and with other simulation results. Simulations performed in the NVT ensemble are shown to give equivalent solvation free energies, and this ensemble is used for the interfacial simulations. Solvation free energies of Cl(-) and Na(+) ions in ice at 150 K are found to be approximately 30 and approximately 20 kcal mol(-1), respectively, less favorable than for water at room temperature. Near the melting point of the model the solvation of the ions in water is the same (within statistical error) as that measured at room temperature, and in the ice is equivalent and approximately 10 kcal mol(-1) less favorable than the liquid. The free energy of transfer for each ion across ice/water interface is calculated and is in good agreement with the bulk observations for the Cl(-) ion. However, for the model of Na(+) the long-range electrostatic contribution to the free energy was more negative in the ice than the liquid, in contrast with the results observed in the bulk calculations.     

Anisotropic interfacial free energies of the hard-sphere crystal-melt interfaces

We present a reliable method to define the interfacial particles for determining the crystal-melt interface position, which is the key step for the crystal-melt interfacial free energy calculations using capillary wave approach. Using this method, we have calculated the free energies gamma of the fcc crystal-melt interfaces for the hard-sphere system as a function of crystal orientations by examining the height fluctuations of the interface using Monte Carlo simulations. We find that the average interfacial free energy gamma(0) = 0.62 +/- 0.02k(B)T/sigma(2) and the anisotropy of the interfacial free energies are weak, gamma(100) = 0.64 +/- 0.02, gamma(110) = 0.62 +/- 0.02, gamma(111) = 0.61 +/- 0.02k(B)T/sigma(2). The results are in good agreement with previous simulation results based on the calculations of the reversible work required to create the interfaces (Davidchack and Laird, Phys. Rev. Lett. 2000, 85, 4571). In addition, our results indicate gamma(100) > gamma(110) > gamma(111) for the hard-sphere system, similar to the results of the Lennard-Jones system.     

Chemical oscillation with periodic adsorption and desorption of surfactant ions at a water/nitrobenzene interface.

Chemical oscillations with periodic adsorption and desorption of surfactant ions, alkyl sulfate ions, at a water/nitrobenzene interface have been investigated. The interfacial tension was measured with a quasi elastic laser scattering (QELS) method and the interfacial electrical potential was obtained. We found that this oscillation consists of a series of abrupt adsorptions of ions, followed by a gradual desorption. In addition, we observed that each abrupt adsorption was always accompanied by a small waving motion of the liquid interface. From the analysis of the video images of the liquid interface or bulk phase, we could conclude that each abrupt adsorption is caused by nonlinear amplification of mass transfer of ions from the bulk phase to the liquid interface by a Marangoni convection, which was generated due to local adsorption of the surfactant ions at the liquid interface that resulted in the heterogeneity of the interfacial tension. In the present paper, we describe the mechanism of the chemical oscillation in terms of the hydrodynamic effect on the ion adsorption processes, and we also show the interfacial chemical reaction with ion exchange during the ion desorption process.     

Effects of ion solvation on phase equilibrium and interfacial tension of liquid mixtures

We study the bulk thermodynamics and interfacial properties of electrolyte solution mixtures by accounting for electrostatic interaction, ion solvation, and inhomogeneity in the dielectric medium in the mean-field framework. Difference in the solvation energy between the cations and anions is shown to give rise to local charge separation near the interface, and a finite Galvani potential between two coexisting solutions. The ion solvation affects the phase equilibrium of the solvent mixture, depending on the dielectric constants of the solvents, reflecting the competition between the solvation energy and translation entropy of the ions. Miscibility is decreased if both solvents have low dielectric constants and is enhanced if both solvents have high dielectric constant. At the mean-field level, the ion distribution near the interface is determined by two competing effects: accumulation in the electrostatic double layer and depletion in a diffuse interface. The interfacial tension shows a nonmonotonic dependence on the salt concentration: it increases linearly with the salt concentration at higher concentrations and decreases approximately as the square root of the salt concentration for dilute solutions, reaching a minimum near 1 mM. We also find that, for a fixed cation type, the interfacial tension decreases as the size of anion increases. These results offer qualitative explanations within one unified framework for the long-known concentration and ion size effects on the interfacial tension of electrolyte solutions.     

Air-liquid interfaces of imidazolium-based [TF2N-] ionic liquids: insight from molecular dynamics simulations

We present molecular dynamics simulations of the air-liquid interface for three room temperature ionic liquids with a common anion: bis(trifluoromethylsulfonyl) imide ([Tf(2)N]), and imidazolium-based cations that differ in the alkyl tail length: 1-butyl-3-methylimidazolium ([C(4)mim]), 1-hexyl-3-methylimidazolium ([C(6)mim]), and 1-octyl-3-methylimidazolium ([C(8)mim]). The CHARMM type force field is used with the partial charges based on quantum calculations for isolated ion pairs. The total charge on cations and anions is around 0.9e and -0.9e, respectively, which somewhat mimics the anion to cation charge transfer and many-body effects. The surface tension at 300 K is computed using the mechanical route and its value slightly overpredicts experimental values. The air-liquid interface is analyzed using the intrinsic method of Identification of the Truly Interfacial Molecules. Structural and dynamic properties of the interfacial, sub-interfacial and central layers are determined. To describe the structure of the interface, we compute the surface roughness, number density and charge density profiles, and orientation ordering of the ions. We further determine the survival probability, normal and lateral self-diffusion coefficients, and re-orientation correlation functions to characterize the dynamics of the cations and anions in the layers. We found a significant enhancement of the cation density and preferential orientation ordering of both the cations and anions at the interface. Overall, the surface of the interfacial layer is smoother than the surface of the sub-interfacial layer and the roughness of both the interfacial and sub-interfacial layers increases with the increase of the length of the cation alkyl tail. Finally, the ions stay considerably longer in the interfacial layer than in the sub-interfacial layer and dynamics of exchange of the ions between the consecutive layers is related to the distinct diffusion and re-orientation dynamics behavior of the ions within the layers.     

Grid positioning independence and the reduction of self-energy in the solution of the Poisson—Boltzmann equation

A common problem in the solution of the Poisson–Boltzmann equation using finite difference methods is the self-energy of the system, also known as the grid energy. Because atoms are typically modeled as a point charge, the infinite self-energy of a point charge is likewise modeled. In this article, a simple, alternate treatment of atomic charge is described where each atom is represented as a sphere of uniform charge. Unlike the point charge model, this method converges as the grid spacing is reduced. The uniform charge model generates the same electrostatic field outside the atoms. In addition, the use of fine grids reduces the variations in the potential due to variations in the position of atoms relative to the grid. Calculations of Born ion solvation energies, small-molecule solvation energies, and the electrostatic field of superoxide dismutase are used to demonstrate that this method yields the same results as the point charge model. © John Wiley & Sons, Inc.     

An electrochemical method for determination of the standard Gibbs energy of anion transfer between water and n-octanol

The transfer of the ions Cl⁻, Br⁻, I⁻, ClO₄⁻, SCN⁻, NO₃⁻, BF₄⁻, and (C₆H₅)₄B⁻ across the water|n-octanol (W|OC) liquid interface was studied and the standard Gibbs energies of ion transfer were determined. The ion transfer was achieved by oxidation of decamethylferrocene dissolved in a droplet of n-octanol that was attached to a graphite electrode immersed in the aqueous solutions of the respective alkali salts of the anions. The electrode reaction can be described by the equation: dmfc_(OC)+X_(W)⁻⇄dmfc⁺_(OC)+X⁻_(OC)+e⁻, where X⁻ is the transferred anion. Square-wave voltammetry at this three-phase arrangement was utilised to determine the formal potential of the decamethylferrocene/decamethylferrocenium (dmfc/dmfc⁺) couple under the condition of ion transfer across the water|n-octanol interface. For calibration the standard Gibbs energies of ion transfer have been extrapolated to octanol from the series of known data for methanol, ethanol, n-propanol, and n-butanol. All these data are consistent and the experimental dependence of the formal potentials on the standard Gibbs energies is as predicted by theory. The validity of data is further supported by calculations of Gibbs energies of ion transfer using the Born theory. Until now it was not possible to perform electrochemical measurements at the water|n-octanol interface because in the conventional four-electrode cells this interface cannot be polarised. With the new method it is now for the first time possible to determine the Gibbs energies of transfer of ions across the water|n-octanol interface. These values are of very wide use for assessing the lipophilicity of compounds in chemistry, medicine, and pharmacology.     

Unified molecular picture of the surfaces of aqueous acid, base, and salt solutions

The molecular structure of the interfacial regions of aqueous electrolytes is poorly understood, despite its crucial importance in many biological, technological, and atmospheric processes. A long-term controversy pertains between the standard picture of an ion-free surface layer and the strongly ion specific behavior indicating in many cases significant propensities of simple inorganic ions for the interface. Here, we present a unified and consistent view of the structure of the air/solution interface of aqueous electrolytes containing monovalent inorganic ions. Molecular dynamics calculations show that in salt solutions and bases the positively charged ions, such as alkali cations, are repelled from the interface, whereas the anions, such as halides or hydroxide, exhibit a varying surface propensity, correlated primarily with the ion polarizability and size. The behavior of acids is different due to a significant propensity of hydronium cations for the air/solution interface. Therefore, both cations and anions exhibit enhanced concentrations at the surface and, consequently, these acids (unlike bases and salts) reduce the surface tension of water. The results of the simulations are supported by surface selective nonlinear vibrational spectroscopy, which reveals among other things that the hydronium cations are present at the air/solution interface. The ion specific propensities for the air/solution interface have important implications for a whole range of heterogeneous physical and chemical processes, including atmospheric chemistry of aerosols, corrosion processes, and bubble coalescence.     

Enhancement of ion transmission at low collision energies via modifications to the interface region of a spectrometer

The transmission efficiency of precursor and product ions decreases significantly at lower collision energies in a four-sector tandem mass spectrometer. In an effort to improve the overall ion transmission in this energy regime three modifications were made in the interface region between the two stages of mass analysis. An einzel lens was inserted prior to the deceleration lens of the collision cell block to reduce the precursor ion beam diameter. The collision cel1 block was reduced in thickness while maintaining the collision path length, thus increasing the number of ions which entered and exited the gas chamber, while removing any stray electrical fields. Finally, a second active focusing element was incorporated after the collision cell block to enhance the collection efficiency of the product ions. A tandem mass spectrum of angiotensin I obtained with this interface, at a collision cell block potential of 9200 volts, exhibited classical high energy collision-induced dissociation (CID) fragmentation patterns, a precursor ion transmission of 92% and an overall CID efficiency of approximately 7.5%. These improvements have resulted in a dramatically higher overall ion transmission at high collision cell potentials as well as sufficient sensitivity in acquiring good quality CID spectra in the lower collision energy regime (i.e., 60 eV). (460-469)     

Vibrational–rotational energies of all H2 isotopomers using Monte Carlo methods

Using variational Monte Carlo techniques, we have computed several of the lowest rotational–vibrational energies of all the hydrogen molecule isotopomers (H₂, HD, HT, D₂, DT, and T₂). These calculations do not require the excited states to be explicitly orthogonalized. We have examined both the usual Gaussian wave function form as well as a rapidly convergent Padé form. The high‐quality potential energy surfaces used in these calculations are taken from our earlier work and include the Born–Oppenheimer energy, the diagonal correction to the Born–Oppenheimer approximation, and the lowest‐order relativistic corrections at 24 internuclear points. Our energies are in good agreement with those determined by other methods. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

A bulk adjusted linear combination of atomic orbitals (BA-LCAO) approach for nanoparticles

We describe a bulk adjusted linear combination of atomic orbitals (BA-LCAO) approach for nanoparticles. In this method, we apply a many-body scaling function (in similar manner as in the environment-modified total energy based tight-binding method) to the DFT-derived diatomic AO interaction potentials (like in the conventional orbital-based density-functional tight binding approach) strictly according to atomic valences acquired naturally in a bulk structure. This modification, (a) facilitates all atom orbital-based electronic structure calculations of charge carrier dynamics in nanoscale structures with a molecular acceptor, and (b) allows to closely match high-level density functional calculation data (previously adjusted to the available experimental findings) for bulk structures. To advance practical application of the BA-LCAO approach we parameterize the Hamiltonian of wurtzite CdSe by fitting its band structure to a high-level DFT reference, corrected for experimentally measured band edges. Here, unlike in conventional DFTB approach, we: (1) use hydrogen-like AOs for the basis as exact atomic eigenfunctions, while orbital energies of which are taken from experimentally measured ionization potentials, and (2) parameterize the many-body scaling functions rather than the atomic wavefunctions. Development of this approach and parameters is guided by our goals to devise a method capable of simultaneously treating the problems of (i) interfacial electron/hole transfer between finite, variable size nanoparticles and electron scavenging molecules, and (ii) high-energy electronic transitions (Auger transitions) that mediate multi-exciton decay in quantum dots. Electronic structure results are described for CdSe quantum dots of various sizes. © 2018 Wiley Periodicals, Inc.     

Ion size effects on the dielectric and electrokinetic properties in aqueous colloidal suspensions

One of the main assumptions of the classical theory most widely used to characterize electrokinetic phenomena is that ions behave as point-like entities. While the realization of the importance of the finite ion size goes back to Stern, 1924, it was Bikerman who presented in 1942 the first expression for the steric interactions among ions. Even now, this is the most often used expression, mainly due to its analytic simplicity. However, once ions are considered to have a finite size, other consequences besides the steric interactions have to be considered. For example, the finite closest approach distance of ions to the interface, the dielectrophoretic force acting on ions in a non-uniform electric field, the variation of the electrolyte solution permittivity with the local ion concentration, and the corresponding Born force acting on the ions, have to be taken into account. In this work, we examine these items in detail and discuss the main contributions made in this field. They show that even for the relatively low surface charge and electrolyte concentration values encountered in colloidal suspension studies, corrections to the classical theory due to ion size effects are far from negligible.