A molecular-dynamics study of a model S(N)1 dissociation reaction at the water liquid/vapor interface

The thermodynamics and dynamics of a model S(N)1 reaction: t-BuCl --> t-Bu+ + Cl- is studied at the water liquid/vapor interface using molecular-dynamics computer simulations. The empirical valence bond approach is used to couple two diabatic states, covalent and ionic, in the electronically adiabatic limit. Umbrella sampling calculations are used to calculate the potential of mean force along the reaction coordinate (defined as the t-Bu to Cl distance) in bulk water and in several locations at the interface. We find a significant increase of the dissociation barrier height and of the reaction free energy at the interface relative to the bulk. This is shown to be due to the reduced polarity of the interface. Reactive flux correlation function calculations show significant deviation of the rate constant from the transition-state theory: The transmission coefficients range from 0.49 in the bulk to 0.05 above the Gibbs surface. The low transmission coefficient at the interface despite the lower friction is shown to be due to slow vibrational relaxation.     

Ion-transfer voltammetry of local anesthetics at an organic solvent/water interface and pharmacological activity vs. ion partition coefficient relationship.

The ion-transfer reaction of local anesthetics at an organic solvent/water interface has been studied using cyclic voltammetry (CV) with a stationary nitrobenzene (NB)/water (W) interface. Procaine and seven other local anesthetics gave reversible or quasi-reversible voltammograms at the NB/W interface in the pH range between 0.9 and 9.6. These drugs are present in aqueous solution in either neutral or ionic form, or both forms. The half-wave potential, as determined by the midpoint potential in CV, vs. pH curves, were determined and analyzed to determine the partition coefficients of both neutral and ionic forms of the drugs between NB and W. The partition coefficients of the ionic forms were derived from their formal potential of transfer at an NB/W interface. The dissociation constants of ionic forms of the drugs in NB were also deduced. A high correlation between the pharmacological activity and the partition coefficient of the ionic form of amide-linked local anesthetics has been shown.     

On the ion‐pair dissociation mechanisms in the small NaCl·(H2O)6 cluster: A perspective from reaction path search calculations

We performed reaction path search calculations for the NaCl·(H₂O)₆ cluster using the global reaction route mapping (GRRM) code to understand the atomic‐level mechanisms of the NaCl → Na⁺ + Cl⁻ ionic dissociation induced by water solvents. Low‐lying minima, transition states connecting two local minima and corresponding intrinsic reaction coordinates on the potential energy surface are explored. We found that the Na Cl distances at the transitions states for the dissociation pathways were distributed in a relatively wide range of 2.7–3.7 Å and that the Na Cl distance at the transition state did not correlate with the commonly used solvation coordinates. This suggests that the definition of the transition states with specific structures as well as good reaction coordinate is very difficult for the ionic dissociation process even in a small water cluster. © 2018 Wiley Periodicals, Inc.     

Vinyl oxidative coupling as a synthetic route to catalytically active monovalent chromium

Reaction of the deprotonated form of cis-{(t-Bu)N(H)P[μ-N(t-Bu)](2)PN(H)(t-Bu)} with CrCl(3)(THF)(3) afforded the trivalent cis-{(t-Bu)NP[μ-N(t-Bu)](2)PN(t-Bu)}[Li (THF)])CrCl(2) (1). Subsequent reaction with 2 equiv of vinyl Grignard (CH(2)=CH)Mg Cl gave the butadiene derivative (cis-{(t-Bu)NP[μ-N(t-Bu)](2)PN(t-Bu)}[Li(THF)])Cr(cis-η(4)-butadiene) (3) formally containing the metal in its monovalent state. The presence of the monovalent state was thereafter confirmed by DFT calculations. The coordination of the butadiene unit appears to be rather robust since reaction with Me(3)P afforded cleavage of the dimeric ligand core but not its displacement. The reaction formed the new butadiene complex [(t-Bu)N-P-N(t-Bu)]Cr(cis-η(4)-butadiene)PMe(3) (4) containing a regular NPN monoanion. In agreement with the presence of monovalent chromium, complexes 3 and 4 act as single-component self-activating catalysts for selective ethylene trimerization and dimerization, respectively.     

Solvation of uranyl(II) and europium(III) cations and their chloro complexes in a room-temperature ionic liquid. A theoretical study of the effect of solvent "humidity"

We report a molecular dynamics study of the solvation of the UO2(2+) and Eu3+ cations and their chloro complexes in the [BMI][PF6][H2O] "humid" room-temperature ionic liquid (IL) composed of 1-butyl-3-methylimidazolium+ and PF6- ions and H2O in a 1:1:1 ratio. When compared to the results obtained in dry [BMI][PF6], the present results reveal the importance of water. The "naked" cations form UO2(H2O)5(2+) and Eu(H2O)9(3+) complexes, embedded in a shell of 7 and 8 PF6- anions, respectively. All studied UO2Cln(2-n) and EuCln(3-n) chloro complexes remain stable during the dynamics and coordinate additional H2O molecules in their first shell. UO2Cl4(2-) and EuCl6(3-) are surrounded by an "unsaturated" water shell, followed by a shell of BMI+ cations. According to an energy component analysis, the UO2Cl4(2-) and EuCl6(3-) species, intrinsically unstable toward dissociation, are more stable than their less halogenated analogues in the IL solution, due to the solvation forces. The different chloro species also interact better with the humid than with the dry IL, which hints at the importance of solvent humidity to improve their solubility. Humidity markedly modifies the local ion environment, with major consequences as far as their spectroscopic properties are concerned. We finally compare the aqueous interface of [BMI][PF6] and [OMI][PF6] ionic liquids, demonstrating the importance of imidazolium substituents (N-butyl versus N-octyl) to the nature of the interface and miscibility with water.     

SECM study of solute partitioning and electron transfer at the ionic liquid/water interface

Molecular partitioning and electron-transfer kinetics have been studied at the ionic liquid/water (IL/water) interface by scanning electrochemical microscopy (SECM). The ionic liquid C8mimC1C1N is immiscible with water and forms a nonpolarizable interface when in contact with it. Partitioning of ferrocene (Fc) across the IL/water interface was studied by SECM and found to be kinetically fast with a partition coefficient CIL/CW of 2400:1. The partition coefficient value was measured by SECM under quasi-steady-state conditions without waiting for complete solute equilibration. To investigate the kinetics of the electron transfer (ET) between aqueous ferricyanide and Fc dissolved in IL, a new approach to the analysis of the SECM current-distance curves was developed to separate the contributions of Fc partitioning and the ET reaction to the tip current. Several combinations of different aqueous and nonaqueous redox species were investigated; however, only the Fc/Fe(CN)63- system behaved according to the Butler-Volmer formalism over the entire accessible potential range.     

Kinetics of ion transfer at the ionic liquid/water nanointerface

Ion transfer (IT) processes in ionic liquids (ILs) are essential for their applications in electrochemical systems and chemical separations. In this Article, the first measurements of IT kinetics at the IL/water interface are reported. Steady-state voltammetry was performed at the nanometer-sized polarizable interface between water and ionic liquid, [THTDP(+)][C(4)C(4)N(-)], immiscible with it that was formed at the tip of a nanopipet. Kinetic measurements at such interfaces are extremely challenging because of slow mass-transfer rates in IL, which is ～700 times more viscous than water. The recently developed new mode of nanopipet voltammetry, common ion voltammetry, was used to overcome technical difficulties and ensure the reliability of the extracted kinetic parameters of IT. The results suggest that the rate of interfacial IT depends strongly on solution viscosity. Voltammetric responses of nanopipets of different radii were analyzed to evaluate the effect of the electrical double layer at the liquid/liquid interface on IT kinetics. The possibility of the influence of the charged pipet wall on ion transport was investigated by comparing currents produced by cationic and anionic species. Possible effects of relaxation phenomena at the IL/water interface on IT voltammograms have also been explored.     

Probing the neutral graphene-ionic liquid interface: insights from molecular dynamics simulations

We study basic mechanisms of the interfacial layer formation at the neutral graphite monolayer (graphene)-ionic liquid (1,3-dimethylimidazolium chloride, [dmim][Cl]) interface by fully atomistic molecular dynamics simulations. We probe the interface area by a spherical probe varying the charge (-1e, 0, +1e) as well as the size of the probe (diameter 0.50 nm and 0.38 nm). The molecular modelling results suggest that: there is a significant enrichment of ionic liquid cations at the surface. This cationic layer attracts Cl(-) anions that leads to the formation of several distinct ionic liquid layers at the surface. There is strong asymmetry in cationic/anionic probe interactions with the graphene wall due to the preferential adsorption of the ionic liquid cations at the graphene surface. The high density of ionic liquid cations at the interface adds an additional high energy barrier for the cationic probe to come to the wall compared to the anionic probe. Qualitatively the results from probes with diameter 0.50 nm and 0.38 nm are similar although the smaller probe can approach closer to the wall. We discuss the simulation results in light of available experimental data on the interfacial structure in ionic liquids.     

Hydrogen bond dynamics at water/organic liquid interfaces

Hydrogen bond dynamics at the neat interface between water and a series of organic liquids are studied with molecular dynamics computer simulation. The organic liquids are nonpolar (carbon tetrachloride), weakly polar (1,2-dichloroethane), and polar (nitrobenzene). The effect of surface polarity and surface roughness is examined. The dynamics are expressed in terms of the hydrogen bond population autocorrelation functions and are found to be nonexponential and strongly dependent on the nature of the organic phase. In particular, at all interfaces, the dynamics are slower at the interface than in the bulk and sensitive to the location of the water molecules along the interface normal.     

Femtosecond time-resolved electronic sum-frequency generation spectroscopy: a new method to investigate ultrafast dynamics at liquid interfaces

We developed a new surface-selective time-resolved nonlinear spectroscopy, femtosecond time-resolved electronic sum-frequency generation (TR-ESFG) spectroscopy, to investigate ultrafast dynamics of molecules at liquid interfaces. Its advantage over conventional time-resolved second harmonic generation spectroscopy is that it can provide spectral information, which is realized by the multiplex detection of the transient electronic sum-frequency signal using a broadband white light continuum and a multichannel detector. We studied the photochemical dynamics of rhodamine 800 (R800) at the air/water interface with the TR-ESFG spectroscopy, and discussed the ultrafast dynamics of the molecule as thoroughly as we do for the bulk molecules with conventional transient absorption spectroscopy. We found that the relaxation dynamics of photoexcited R800 at the air/water interface exhibited three characteristic time constants of 0.32 ps, 6.4 ps, and 0.85 ns. The 0.32 ps time constant was ascribed to the lifetime of dimeric R800 in the lowest excited singlet (S(1)) state (S(1) dimer) that is directly generated by photoexcitation. The S(1) dimer dissociates to a monomer in the S(1) state (S(1) monomer) and a monomer in the ground state with this time constant. This lifetime of the S(1) dimer was ten times shorter than the corresponding lifetime in a bulk aqueous solution. The 6.4 ps and 0.85 ns components were ascribed to the decay of the S(1) monomer (as well as the recovery of the dimer in the ground state). For the 6.4 ps time constant, there is no corresponding component in the dynamics in bulk water, and it is ascribed to an interface-specific deactivation process. The 0.85 ns time constant was ascribed to the intrinsic lifetime of the S(1) monomer at the air/water interface, which is almost the same as the lifetime in bulk water. The present study clearly shows the feasibility and high potential of the TR-ESFG spectroscopy to investigate ultrafast dynamics at the interface.     

Lactone cleavage reaction kinetics of rhodamine dye at liquid/liquid interfaces studied by micro-two-phase sheath flow/two-photon excitation fluorescence microscopy

Two-photon excitation fluorescence microscopy was combined with the two-phase microflow system in order to measure the fast interfacial reaction rate at liquid/liquid interfaces. The lactone cleavage kinetics of octadecylrhodamine B (C(18)RB) at the toluene/water and heptane/water interfaces was studied by this new method. The organic solution containing the nonfluorescent lactone of C(18)RB was made to flow as an inner flow with an aqueous outer sheath flow. The diameter of the inner flow was <20 microm. A focused fundamental beam of a Ti:sapphire pulse laser of 780 nm was irradiated to the interface, and emitted fluorescence from the fluorescent product was detected by a charge-coupled device (CCD) camera or a streakscope. The increase in the concentration of the fluorescent form of C(18)RB was measured along the interface of the inner flow of the toluene/water and heptane/water systems for 80 micros just after the contact of two phases. The analysis made by the time-dependent Langmuir adsorption model with the aid of the digital simulation method gave the cleavage reaction rate constants of the lactone form of C(18)RB at the liquid/liquid interfaces.     

Heterogeneous electron transfer kinetics at the ionic liquid/metal interface studied using cyclic voltammetry and scanning electrochemical microscopy

The electrochemical behavior of a redox-active, ferrocene-modified ionic liquid (1-ferrocenylmethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide) in acetonitrile and in an ionic liquid electrolyte (1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide) is reported. Reversible electrochemical behavior was observed in each electrolyte with responses typical of those for unmodified ferrocene observed in each medium. In the ionic liquid electrolyte, the diffusion coefficient of the redox-active ionic liquid increased by a factor of 5 upon increasing the temperature from 27 to 90 degrees C. The kinetics of electron transfer across the ionic liquid/electrode interface were studied using cyclic voltammetry, and the standard heterogeneous electron transfer rate constant, k (0) was determined to be 4.25 x 10 (-3) cm s (-1). Scanning electrochemical microscopy was then also used to probe the heterogeneous kinetics at the interface between the ionic liquid and the solid electrode and conventional kinetic SECM theory was used to determine k (0). The k (0) value obtained using SECM was higher than that determined using cyclic voltammetry. These results indicate that SECM is a very useful technique for studying electron transfer dynamics in ionic liquids.     

Solvation of sodium chloride in the 1-butyl-3-methyl-imidazolium bis(trifluoromethylsulfonyl)imide ionic liquid: a molecular dynamics study

We report molecular dynamics studies on the solvation of sodium chloride in the 1-butyl-3-methyl-imidazolium bis(trifluoromethylsulfonyl)imide ionic liquid ([BMI][Tf2N] IL). We first consider the potential of mean force for dissociating a single Na+Cl- ion pair, showing that the latter prefers to be undissociated rather than dissociated (by ca. 9 kcal/mol), with a free energy barrier of ca. 5 kcal/mol (at d approximately 5.2 A) for the association process. The preference for Na+Cl- association is also observed from a 100 ns molecular dynamics simulation of a concentrated solution, where the Na+Cl- ions tend to form oligomers and microcrystals in the IL. Conversely, the simulation of Na13Cl14- and Na14Cl13+ cubic microcrystals (with, respectively, Cl- and Na+ at the vertices) does not lead to dissolution in the IL. Among these, Na14Cl13+ is found to be better solvated than Na13Cl14-, mainly due to the stronger Na+...Tf2N- interactions as compared to the Cl-...BMI+ interactions at the vertices of the cube. We finally consider the solid/liquid interface between the 100 face of NaCl and the IL, revealing that, in spite of its polar nature, the crystal surface is solvated by the less polar IL components (CF3(Tf2N) and butyl(BMI) groups) rather than by the polar ones (O(Tf2N) and imidazolium(BMI) ring). Specific ordering at the interface is described for both Tf2N- anions and BMI+ cations. In the first IL layer, the ions are rather parallel to the surface, whereas in the second "layer" they are more perpendicular. A similar IL structure is found at the surface of the all-neutral Na0Cl0 solid analogue, confirming that the solvation of the crystal is rather "apolar", due to the mismatch between the IL and the crystal ions. Several comparisons with water, methanol, or different BMI+-based ILs as solvents are presented, allowing us to better understand the specificity of the ionic liquid-NaCl interactions.     

Theoretical studies of solute vibrational energy relaxation at liquid interfaces

Recent advances in the theoretical understanding of solute vibrational energy relaxation at liquid interfaces and surfaces are described. Non-equilibrium molecular dynamics simulations of the relaxation of an initially excited solute molecule are combined with equilibrium force autocorrelation calculations to gain insight into the factors that influence the vibrational relaxation rate. Diatomic and triatomic nonpolar, polar, and ionic solute molecules adsorbed at the liquid/vapor interface of several liquids as well as at the water/CCl(4) liquid/liquid interface are considered. In general, the vibrational relaxation rate is significantly slower (a factor of 3 to 4) at the liquid/vapor and liquid/liquid interface than in the bulk due to the reduced density, which gives rise to a reduced contribution of the repulsive solvent-solute forces on the vibrational mode. The surface effects on the ionic solutes are much smaller (50% or less slower relaxation relative to the bulk). This is due to the fact that ionic solutes at the interface are able to keep part of their solvation shell to a degree that depends on their size. Thus, a significant portion of the repulsive forces is maintained. A high degree of correlation is found between the peak height of the solvent-solute radial distribution function and the vibrational relaxation rate. The relaxation rate at the liquid/liquid interface strongly depends on the location of the solute across the interface and correlates with the change in the density and polarity profile of the interface.     

Rotational dynamics of water and benzene controlled by anion field in ionic liquids: 1-butyl-3-methylimidazolium chloride and hexafluorophosphate

The rotational correlation time (tau(2R)) is determined for D(2)O (polar) and C(6)D(6) (apolar) in 1-butyl-3-methylimidazolium chloride ([bmim][Cl]) and hexafluorophosphate ([bmim][PF(6)]) by measuring (2)H (D) nuclear magnetic resonance spin-lattice relaxation time (T(1)) in the temperature range from -20 to 110 degrees C. The tau(2R) ratio of water to benzene (tau(WB)) was used as a measure of solute-solvent attraction. tau(WB) is 0.73 and 0.52 in [bmim][Cl] and [bmim][PF(6)], respectively, whereas the molecular volume ratio is as small as 0.11. The slowdown of the water dynamics compared to the benzene dynamics in ionic liquids is interpreted by the Coulombic attractive interaction between the polar water molecule and the anion. As for the anion effect, the rotational dynamics of water solvated by Cl(-) is slower than that solvated by PF(6) (-), whereas the rotational dynamics of benzene is similar in the two ionic liquids. This is interpreted as an indication of the stronger solvation by the anion with a larger surface charge density. The slowdown of the water dynamics via Coulombic solvation is actually significant only at water concentrations lower than approximately 9 mol dm(-3) at room temperature, and it is indistinguishable at temperatures above approximately 100 degrees C. The quadrupolar coupling constants determined for D(2)O and C(6)D(6) in the ionic liquids were smaller by a factor of 2-3 than those in the pure liquid state.     

Structure of water at ionic liquid/Ag interface probed by surface enhanced Raman spectroscopy

The potential-dependent adsorption behavior of water and ionic liquid was probed by surface-enhanced Raman spectroscopy (SERS) at the Ag electrode surface in the ionic liquids containing water with different concentrations.The configuration of water at the ionic liquid/electrode interface and the relationship between the potential of zero charge (pzc) and the molar fraction of water were deduced through the changes in the vibrational frequency of OH stretching mode.A small Stark effect value was determined ...     

Excess electron relaxation dynamics at water/air interfaces

We have performed mixed quantum-classical molecular dynamics simulations of the relaxation of a ground state excess electron at interfaces of different phases of water with air. The investigated systems included ambient water/air, supercooled water/air, Ih ice/air, and amorphous solid water/air interfaces. The present work explores the possible connections of the examined interfacial systems to finite size cluster anions and the three-dimensional infinite, fully hydrated electron. Localization site analyses indicate that in the absence of nuclear relaxation the electron localizes in a shallow potential trap on the interface in all examined systems in a diffuse, surface-bound (SB) state. With relaxation, the weakly bound electron undergoes an ultrafast localization and stabilization on the surface with the concomitant collapse of its radius. In the case of the ambient liquid interface the electron slowly (on the 10 ps time scale) diffuses into the bulk to form an interior-bound state. In each other case, the excess electron persists on the interface in SB states. The relaxation dynamics occur through distinct SB structures which are easily distinguishable by their energetics, geometries, and interactions with the surrounding water bath. The systems exhibiting the most stable SB excess electron states (supercooled water/air and Ih ice/air interfaces) are identified by their characteristic hydrogen-bonding motifs which are found to contain double acceptor-type water molecules in the close vicinity of the electron. These surface states correlate reasonably with those extrapolated to infinite size from simulated water cluster anions.     

Photodissociation of ICN at the liquid/vapor interface of water

The photodissociation of ICN adsorbed at the liquid/vapor interface of water is studied using classical molecular dynamics with nonadiabatic surface hopping. The cage escape, geminate recombination to form ICN and INC and the subsequent vibrational relaxation of these two molecules (on their ground electronic states) is compared with the same process in bulk water and with previous photodissociation studies at liquid interfaces. We find that the reduced surface density and weaker solvent-solute interactions give rise to reduced rate of nonadiabatic transitions and that the probability for cage escape at the interface is significantly enhanced due to the possibility that one or both of the photodissociation fragments desorb into the gas phase. The overall desorption probability varies from 75% to 92% for ICN initially located just below the Gibbs surface (50% bulk density) to ICN located just above the Gibbs surface, respectively. The corresponding geminate recombination probabilities are 18% and 9%, respectively. The vibrational relaxation rate of the recombined ICN is slower than in the bulk by a factor of 2.3.     

Electrochemical behaviour of EOP at the liquid/liquid interface and its determination

A new approach to the determination of emulsifier OP (EOP, a nonionic surfactant) is proposed, based on the voltammetry of the EOP facilitating Ca(2+) transfer across the water/nitrobenzene (w/nb) interface. The peaks obtained with cyclic voltammetry at the w/nb interface correspond to the transfer of the coordination complex ion [Ca(EOP)](2+) from the aqueous to the nonaqueous phase. The calibration graph is linear in the range of 1.5x10(-5)-5x10(-4) mol l(-1). The ionic surfactants have a very minor effect on detection.