Computer Simulation of Cation-Exchange Membrane Structure: An Elementary Act of Hydrated Ion Transport

A fragment of the structure of a sulfo cation exchanger, which is the basis of most cation-exchange membranes, is calculated by an ab initio method. An analysis of interatomic bonds in the structure shows that, to detach a mobile ion from a fixed ion, it is necessary to break the hydrogen bond between hydration water molecules of the counterion in addition to overcoming the electrostatic attraction. As the hydrogen bond cleavage work for simple hydrated ions is ten times the electrostatic attraction energy, an elementary act of ion transport in a cation-exchange membrane is considered mostly as the hydrogen bond transfer reaction.     

On the structure and dynamics of the hydrated sulfite ion in aqueous solution--an ab initio QMCF MD simulation and large angle X-ray scattering study

Theoretical ab initio quantum mechanical charge field molecular dynamics (QMCF MD) formalism has been applied in conjunction to experimental large angle X-ray scattering to study the structure and dynamics of the hydrated sulfite ion in aqueous solution. The results show that there is a considerable effect of the lone electron-pair on sulfur concerning structure and dynamics in comparison with the sulfate ion with higher oxidation number and symmetry of the hydration shell. The S-O bond distance in the hydrated sulfite ion has been determined to 1.53(1) ? by both methods. The hydrogen bonds between the three water molecules bound to each sulfite oxygen are only slightly stronger than those in bulk water. The sulfite ion can therefore be regarded as a weak structure maker. The water exchange rate is somewhat slower for the sulfite ion than for the sulfate ion, τ(0.5) = 3.2 and 2.6 ps, respectively. An even more striking observation in the angular radial distribution (ARD) functions is that the for sulfite ion the water exchange takes place in close vicinity of the lone electron-pair directed at its sides, while in principle no water exchange did take place of the water molecules hydrogen bound to sulfite oxygens during the simulation time. This is also confirmed when detailed pathway analysis is conducted. The simulation showed that the water molecules hydrogen bound to the sulfite oxygens can move inside the hydration shell to the area outside the lone electron-pair and there be exchanged. On the other hand, for the hydrated sulfate ion in aqueous solution one can clearly see from the ARD that the distribution of exchange events is symmetrical around the entire hydration sphere.     

Structure and dynamics of N,N-diethyl-N-methylammonium triflate ionic liquid, neat and with water, from molecular dynamics simulations

We investigated by means of molecular dynamics simulations the properties (structure, thermodynamics, ion transport, and dynamics) of the protic ionic liquid N,N-diethyl-N-methylammonium triflate (dema:Tfl) and of selected aqueous mixtures of dema:Tfl. This ionic liquid, a good candidate for a water-free proton exchange membrane, is shown to exhibit high ion mobility and conductivity. The radial distribution functions reveal a significant long-range structural correlation. The ammonium cations [dema](+) are found to diffuse slightly faster than the triflate anions [Tfl](-), and both types of ions exhibit enhanced mobility at higher temperatures, leading to higher ionic conductivity. Analysis of the dynamics of ion pairing clearly points to the existence of long-lived contact ion pairs. We also examined the effects of water through characterization of properties of dema:Tfl-water mixtures. Water molecules replace counterions in the coordination shell of both ions, thus weakening their association. As water concentration increases, water molecules start to connect with each other and then form a large network that percolates through the system. Water influences ion dynamics in the mixtures. As the concentration of water increases, both translational and rotational motions of [dema](+) and [Tfl](-) are significantly enhanced. As a result, higher vehicular ionic conductivity is observed with increased hydration level.     

Mean force potential of interaction between Na<Superscript>+</Superscript> and Cl<Superscript>−</Superscript> ions in planar nanopores in contact with water under pressure

The mean force potential (MFP) of interaction between counterions Na⁺ and Cl^? in a planar nanopore with structureless hydrophobic walls is calculated via computer simulation under the condition that the nanopore is in contact with water at an external pressure that exceeds the saturation pressure but remains insufficient to fill the nanopore with water. For a nanopore with a liquid phase, the MFP dependence on the interionic distance indicates the dissociation of an ion pair into two hydrated ions in a nanopore that is not completely filled with water. Fluctuations in the number of water molecules drawn into the interionic space decisively influence the dissociation. The attraction between counterions, averaged over thermal fluctuations, depends largely on the pore width and grows as the shielding of the ions’ electric field by water molecules in a narrow pore diminishes. The contributions from energy and entropy to the free energy of hydration are analyzed.     

Equilibrium and Transport Characteristics of Perfluorinated Cation-Exchange Membranes Containing Carboxyl Groups

Equilibrium (the exchange capacity and moisture content) and transport (the ion transport numbers and the electrical conductance) properties of perfluorinated cation-exchange membranes containing carboxyl groups were studied as a function of pH and concentration of KCl solutions (10^–4–1 M). Dissociation constants of carboxyl groups, adsorption potentials of K⁺ ions, transport numbers and mobilities of counterions in the membranes, concentrations of fixed ions, co- and counterions, as well as the Donnan potentials were calculated.     

NMR and SANS studies of aggregation and microemulsion formation by phosphorus fluorosurfactants in liquid and supercritical carbon dioxide

(1)H NMR relaxation and diffusion studies were performed on water-in-CO(2) (W/C) microemulsion systems formed with phosphorus fluorosurfactants of bis[2-(F-hexyl)ethyl] phosphate salts (DiF(8)), having different counterions (Na(+), NH(4)(+), N(CH(3))(4)(+)) by means of high-pressure in situ NMR. Water has a low solubility in CO(2) and is mainly solubilized by the microemulsion droplets formed with surfactants added to CO(2) and water mixtures. There is rapid exchange of water between the bulk CO(2) and the microemulsion droplets; however, NMR relaxation measurements show that the entrapped water has restricted motion, and there is little "free" water in the core. Counterions entrapped by the droplets are mostly associated with the surfactant headgroups: diffusion measurements show that counterions and the surfactant molecules move together with a diffusion coefficient that is associated with the droplet. The outer shell of the microemulsion droplets consists of the surfactant tails with some associated CO(2). For W/C microemulsions formed with the phosphate-based surfactant having the ammonia counterion (A-DiF(8)), the (1)H NMR signal for NH(4)(+) shows a much larger diffusion coefficient than that of the surfactant tails. This apparent paradox is explained on the basis of proton exchange between water and the ammonium ion. The observed dependence of the relaxation time (T(2)) on W(0) (mole ratio of water to surfactant in the droplets) for water and NH(4)(+) can also be explained by this exchange model. The average hydrodynamic radius of A-DiF(8) microemulsion droplets estimated from NMR diffusion measurements (25 degrees C, 206 bar, W(0) = 5) was R(h) = 2.0 nm. Assuming the theoretical ratio of R(g)/R(h) = 0.775 for a solid sphere, where R(g) is the radius of gyration, the equivalent hydrodynamic radius from SANS is R(h) = 1.87 nm. The radii measured by the two techniques are in reasonable agreement, as the two techniques are weighted to measure somewhat different parts of the micelle structure.     

A first principles molecular dynamics study of lithium atom solvation in binary liquid mixture of water and ammonia: structural, electronic, and dynamical properties

The preferential solvation of solutes in mixed solvent systems is an interesting phenomenon that plays important roles in solubility and kinetics. In the present study, solvation of a lithium atom in aqueous ammonia solution has been investigated from first principles molecular dynamics simulations. Solvation of alkali metal atoms, like lithium, in aqueous and ammonia media is particularly interesting because the alkali metal atoms release their valence electrons in these media so as to produce solvated electrons and metal counterions. In the present work, first principles simulations are performed employing the Car-Parrinello molecular dynamics method. Spontaneous ionization of the Li atom is found to occur in the mixed solvent system. From the radial distribution functions, it is found that the Li(+) ion is preferentially solvated by water and the coordination number is mostly four in its first solvation shell and exchange of water molecules between the first and second solvation shells is essentially negligible in the time scale of our simulations. The Li(+) ion and the unbound electron are well separated and screened by the polar solvent molecules. Also the unbound electron is primarily captured by the hydrogens of water molecules. The diffusion rates of Li(+) ion and water molecules in its first solvation shell are found to be rather slow. In the bulk phase, the diffusion of water is found to be slower than that of ammonia molecules because of strong ammonia-water hydrogen bonds that participate in solvating ammonia molecules in the mixture. The ratio of first and second rank orientational correlation functions deviate from 3, which suggests a deviation from the ideal Debye-type orientational diffusion. It is found that the hydrogen bond lifetimes of ammonia-ammonia pairs is very short. However, ammonia-water H-bonds are found to be quite strong when ammonia acts as an acceptor and these hydrogen bonds are found to live longer than even water-water hydrogen bonds.     

Glass transitions of ethyl acrylate-based ionomers

The glass-transition temperatures of a series of copolymers of ethyl acrylate and acrylic acid neutralized with various cations were investigated. It was found that a plot of T_g as a function of ion content, for every type of ion investigated here, gives an unusual sigmoidal curve, which can be correlated with the onset of the failure of time–temperature superposition in viscoelastic studies, as will be shown in a future publication. Also, all of the T_g versus concentration curves for the various counterions can be superposed if the plots are made against cq/a, where c is the metal acrylate content, q the cation charge, and a the distance between centers of charge. Furthermore, in one region of water content, a linear relation is obtained between the glass transition and the water content (in weight-%) independent of the ion concentration over wide ranges of ion content. Finally, above an ion concentration of 12 mole-%, the rate of change in T_g per water molecule per ion pair at constant ion content, (?T_g/?n)_c is linear but with different slopes above and below two water molecules per ion pair.     

Is there a ‘matrix suppression effect’ under fast‐atom bombardment liquid secondary ion mass spectrometry of ionic surfactants in glycerol?

Some features of a ‘matrix suppression effect’ caused by ionic surface‐active compounds under fast‐atom bombardment (FAB) liquid secondary ion mass spectrometry (LSIMS) are being revised. It is shown that abundant transfer of the glycerol matrix molecules to the gas phase does occur under FAB‐LSIMS of ionic surfactants, contrary to popular belief. This process can be obscure because of the dependence of the charge state of the glycerol‐containing cluster ions on the type of ionic surfactant. It is revealed that, while glycerol matrix signals are really completely suppressed in the positive ion mass spectra of cationic surfactants (decamethoxinum, aethonium), abundant deprotonated glycerol and glycerol‐anion clusters are recorded in the negative ion mode. In the case of an anionic surfactant (sodium dodecyl sulfate), on the contrary, glycerol is completely suppressed in the negative ion mode, but is present in the protonated and cationized forms in the positive ion mass spectra. It is suggested that such patterns of positive and negative ion FAB‐LSIMS spectra of ionic surfactants solutions reflect the structure and composition of the electric double layer formed at the vacuum‐liquid interface by organic cations or anions and their counterions. Processes leading to the formation of the glycerol‐containing ions preferentially of positive or negative charge are discussed. The most obvious of them is efficient binding of glycerol to inorganic counterions of the salts Cl^? or Na⁺, which is confirmed by data from quantum chemical calculations. The high content of the counterions and relatively small content of glycerol in the sputtered zone may be responsible for the charge‐selective suppression of neat glycerol clusters of opposite charge to the counterions. In the case of a mixture of cationic and anionic surfactants the substitution of inorganic counterions by organic ones was observed. The dependence of the exchange rate in the surface layer is not a linear function of the bulk solution concentration, and an effect of abrupt recharging of the surface can be registered. No both positively or negatively charged pure glycerol and glycerol‐inorganic counterion clusters are recorded for the mixture. Correlations between the mass spectrometric observations and some phenomena of surface and colloid chemistry and physics are discussed. Copyright © 2007 John Wiley & Sons, Ltd.     

Local structure around chloride ion in anion exchange resin

Neutron diffraction measurements on 35Cl/37Cl isotopically substituted anion exchange resins were carried out in order to obtain direct information on the local structure around the chloride ion absorbed in the resin. Structural parameters concerning the first nearest-neighbor interaction of chloride ions were determined through a least-squares fitting procedure of the observed first-order difference function, DeltaCl(Q). It has been revealed that the chloride ion is neighboring an ion exchange group (-CH2(CH3)3N+) with a Cl-...N distance of 3.10(3) A, and simultaneously bonded with 2.4(1) D2O molecules with a Cl-...D distance of 2.25(2) A. The second and third nearest water molecules around Cl- have also been observed. These results indicate that the direct ionic interaction between Cl- and -CH2(CH3)3N+ drastically reduces the number of first-neighbor water molecules around Cl- but enhances the long-distance structuring of the remaining water molecules in the environment surrounded by a hydrophobic polymer matrix.     

Effect of nonelectrostatic interactions on the macromolecular behavior of polyelectrolyes. II. Influence of urea and dyes on the viscometric and optical behavior of PLL

The effect of urea and of dyes on the viscometric behavior of aqueous solutions of PLL has been investigated. The compact structure, which is characteristic of PLL in pure water, is found to disappear progressively as the amount of urea is increased. This is interpreted as due to the destruction in such media of nonelectrostatic interactions responsible for the stabilization of the compact structure of PLL in pure water. In 8M urea, the macromolecular behavior of PLL during the ionization is dependent only on the repulsive interactions between the charged groups. The extension of PLL is then practically independent of the nature of counterions, contrary to what was observed in pure water. In the presence of dye (acridine orange), the compact form of PLL is stabilized, and the dye is found to bind for the two structures of PLL. The analogy between the behavior of Bu₄N⁺ and that of dye is in agreement with the fact the Bu₄N⁺ leads to a stabilization of the compact structure by way of specific interactions between N-tetraalkylammonium counterions and the PLL chain through the structure of water. The hydrodynamic behavior of PLL is in good agreement with data obtained from potentiometric and optical activity measurements.     

Effect of nonelectrostatic interactions on the macromolecular behavior of polyelectrolytes. I. Potentiometric titrations of PLL in water and aqueous–organic mixtures

The potentiometric and viscosity behavior of salt-free solutions of the polycondensate between L-lysine and 1,3-benzenedisulfonyl chloride (PLL) and related polyelectrolytes, have been investigated in water and in aqueous–organic mixtures (acetone–water). The effects of counterions and solvent on the pH-induced conformational transition is related to hydrophobic interactions. The transition from a compact state to a more extended one shifts towards high values of neutralization with increasing length of N-tetraalkylammonium counterions. Moreover, the above transition gradually disappears as the percentage of acetone in the mixture is increased, which is interpreted as reflecting the breakdown of the compact structure due to hydrophobic interactions stabilized by bulky counterions. A linear decrease in the free energy variation ΔG_t° is observed with increasing percentage of acetone up to about 40% of acetone, and a value of 1200 cal/unit is obtained for PLL in water, which is a reasonable value compared to other data for hydrophobic polyelectrolytes.     

Explicit water molecular dynamics study of the mobility of halide ions in presence of ionene oligocations

We present explicit water molecular dynamics simulations of solutions of aliphatic 3,3- and 6,6-ionene oligocations neutralized with (i) fluoride, chloride, bromide, or iodide counterions, respectively, or (ii) with a 1:1 mixture of chloride and bromide anions in presence of a low molecular weight salt at 298 K. The SPC/E model was used to describe water molecules. Results of the simulation are presented in form of the pair distribution functions between various atoms on the ionene oligoion and counterions in solution. In addition, we were interested in the dynamics of counterions around model ionenes. We showed that counterions residing in the vicinity of the oligoion exchange rapidly with those in the bulk solution, with the frequency depending on the nature of the counterion and on the charge density of the oligoion. We calculated the average residence times of the various counterion species to the oligoions and proposed the model which divides the counterions into "free" and "bound" and calculated the fraction of "free" counterions. In the second part of the study, we investigated interaction of the sodium chloride and sodium bromide, being simultaneously present in the solution, with differently charged ionenes in water. The selectivity effect was clearly observed: bromide ions tend to replace chloride ions in the immediate vicinity of the ionene oligoions. Simulation results are discussed in light of our recent measurements of thermodynamic and transport properties of aqueous ionene solutions.     

Solvated calcium ions in charged silica nanopores

Hydroxyl surface density in porous silica drops down to nearly zero when the pH of the confined aqueous solution is greater than 10.5. To study such extreme conditions, we developed a model of slit silica nanopores where all the hydrogen atoms of the hydroxylated surface are removed and the negative charge of the resulting oxygen dangling bonds is compensated by Ca(2+) counterions. We employed grand canonical Monte Carlo and molecular dynamics simulations to address how the Ca(2+) counterions affect the thermodynamics, structure, and dynamics of confined water. While most of the Ca(2+) counterions arrange themselves according to the so-called "Stern layer," no diffuse layer is observed. The presence of Ca(2+) counterions affects the pore filling for strong confinement where the surface effects are large. At full loading, no significant changes are observed in the layering of the first two adsorbed water layers compared to nanopores with fully hydroxylated surfaces. However, the water structure and water orientational ordering with respect to the surface is much more disturbed. Due to the super hydrophilicity of the Ca(2+)-silica nanopores, water dynamics is slowed down and vicinal water molecules stick to the pore surface over longer times than in the case of hydroxylated silica surfaces. These findings, which suggest the breakdown of the linear Poisson-Boltzmann theory, provide important information about the properties of nanoconfined electrolytes upon extreme conditions where the surface charge and ion concentration are large.     

A reversed phase/hydrophilic interaction/ion exchange mixed-mode stationary phase for liquid chromatography

A new stationary phase demonstrated effective separation towards polar analytes or their counterions within a single run.     

New concepts of the mechanism of hydrogen exchange between organic molecules and strong acidic centers

Ab initio calculations of the activation energy of the reaction of hydrogen exchange between the methane molecule and the H₃O⁺ ion in the gaseous phase have been carried out by using Hartree-Fock methods and Moeller-Plesset second order perturbation theory (MP2) methods. The structure of the transition state of this reaction has been found. The interaction of the H₃O⁺ ion with molecules of aliphatic hydrocarbons and amino acids has been studied by the Serniernpirical AMI method.Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 7, pp. 1852–1854, July, 1996.     

H+-Cu2+ ion exchange on a Cu0-KU-23 sulfocation exchanger nanocomposite in solutions with various pH values

The influence of solution pH on H⁺-Cu²⁺ exchange equilibrium was studied for the disperse copper- KU-23 (15/100 S) macroporous strongly acid sulfocation exchanger nanocomposite. In dilute solutions of copper(II) sulfate in weakly acid media (pH 3–6), ion exchange limitation is related to the consumption of hydrogen counterions and the formation of copper counterions in a side reaction between copper particles and impurity oxygen. The concentration of hydrogen counterions in the nanocomposite is replenished as pH decreases because of their sorption from solution together with copper ions. Original Russian Text ? E.V. Zolotukhina, T.A. Kravchenko, 2009, published in Zhurnal Fizicheskoi Khimii, 2009, Vol. 83, No. 5, pp. 934–938.     

Fourier transform infrared spectroscopy of azide and cyanate ion pairs in AOT reverse micelles

Evidence for ion pair formation in aqueous bis(2-ethylhexyl) sulfosuccinate (AOT) reverse micelles (RMs) was obtained from infrared spectra of azide and cyanate with Li(+), Na(+), K(+), and NH(4)(+) counterions. The anions' antisymmetric stretching bands near 2000 cm(-1) are shifted to higher frequency (blueshifted) in LiAOT and to a lesser extent in NaAOT, but they are very similar to those in bulk water with K(+) and NH(4)(+) as the counterions. The shifts are largest for low values of w(o) = [water]/[AOT] and approach the bulk value with increasing w(o). The blueshifts are attributed to ion pairing between the anions and the counterions. This interpretation is reinforced by the similar trend (Li(+)>Na(+)>K(+)) for producing contact ion pairs with the metal cations in bulk dimethyl sulfoxide (DMSO) solutions. We find no evidence of ion pairs being formed in NH(4)AOT RMs, whereas ammonium does form ion pairs with azide and cyanate in bulk DMSO. Studies are also reported for the anions in formamide-containing AOT RMs, in which blueshifts and ion pair formation are observed more than in the aqueous RMs. Ion pairs are preferentially formed in confined RM systems, consistent with the well established ideas that RMs exhibit reduced polarity and a disrupted hydrogen bonding network compared to bulk water and that ion-specific effects are involved in mediating the structure of species at interfaces.     

Interpretation of protein/ligand crystal structure using QM/MM calculations: case of HIV-1 protease/metallacarborane complex

Deltahedral metallacarborane compounds have recently been discovered as potent, specific, stable, and nontoxic inhibitors of HIV-1 protease (PR), the major target for AIDS therapy. The 2.15 A-resolution X-ray structure has exhibited a nonsymmetrical binding of the parental compound [Co(3+)-(C2B9H11)2](-) (GB-18) into PR dimer and a symmetrical arrangement in the crystal of two PR dimer complexes into a tetramer. In order to explore structural and energetic details of the inhibitor binding, quantum mechanics coupled with molecular mechanics approach was utilized. Realizing the close positioning of anionic inhibitors in the active site cavity, the possibility of an exchange of structural water molecules Wat50 and Wat128 by Na+ counterions was studied. The energy profiles for the rotation of the GB-18 molecules along their longitudinal axes in complex with PR were calculated. The results show that two Na+ counterions are present in the active site cavity and provide energetically favorable and unfavorable positions for carbon atoms within the carborane cages. Eighty-one rotamer combinations of four molecules of GB-18 bound to PR out of 4 x 10(5) are predicted to be highly populated. These results lay ground for further calculations of interaction energies between GB-18 and amino acids of PR active site and will make it possible to interpret computationally the binding of similar metallacarborane molecules to PR as well as to resistant PR variants. Moreover, this computational tool will allow the design of new, more potent metallacarborane-based HIV-1 protease inhibitors.