Equilibrium Solvation, Electron-Transfer Reactions, and Stokes-Shift Dynamics in Ionic Liquids

Structural Properties of [N1888][TFSI] Ionic Liquid: A Small Angle Neutron Scattering and Polarizable Molecular Dynamics Study

In this study, we investigate the quaternary ammonium-based ionic liquid (QAIL), methyltrioctylammonium bis(trifluoromethylsulfonyl)imide, [N1888][TFSI], utilizing small angle neutron scattering (SANS) measurements and polarizable molecular dynamics (MD) simulations to characterize the short- and long-range liquid structure. Scattering structure factors show signatures of three length scales in reciprocal space indicative of alternating polarity (k ∼ 0.44 Å-1), charge (k ∼ 0.75 Å-1), and neighboring or adjacent (k ∼ 1.46 Å-1) domains. Excellent agreement between simulation and experimental scattering structure factors validates various simulation analyses that provide detailed atomistic characterization of the different length scale correlations. The first solvation shell structure is illustrated by obtaining radial, angular, dihedral, and combined distribution functions, where two dominant spatial motifs, N+···N- and N+···O-, compete for optimal packing around the polar head of the [N1888]+ cation. Intermediate and long-range structures are governed by the balance between local electroneutrality and octyl chain networking, respectively. By computing the charge-correlation structure factor, SZZ, and the spatial extent of the octyl chain network using graph theory, the bulk-phase structure of [N1888][TFSI] is characterized in terms of electrostatic screening and apolar domain formation length scales.

Redox potentials in ionic liquids: Anomalous behavior?

Redox potentials depend on the nature of the solvent/electrolyte through the solvation energies of the ionic solute species. For concentrated electrolytes, ion solvation may deviate significantly from the Born model predictions due to ion pairing and correlation effects. Recently, Ghorai and Matyushov [J. Phys. Chem. B 124, 3754-3769 (2020)] predicted, on the basis of linear response theory, an anomalous trend in the solvation energies of room temperature ionic liquids, with deviations of hundreds of kJ/mol from the Born model for certain size solutes/ions. In this work, we computationally evaluate ionic solvation energies in the prototypical ionic liquid, 1-butyl-3-methylimidazolium tetrafluoroborate (BMIM/BF4), to further explore this behavior and benchmark several of the approximations utilized in the solvation energy predictions. For comparison, we additionally compute solvation energies within acetonitrile and molten NaCl salt to illustrate the limiting behavior of purely dipolar and ionic solvents. We find that the overscreening effect, which results from the inherent charge oscillations of the ionic liquid, is substantially reduced in magnitude due to screening from the dipoles of the molecular ions. Therefore, for the molten NaCl salt, for which the ions do not have permanent dipoles, modulation of ionic solvation energies from the overscreening effect is most significant. The conclusion is that ionic liquids do indeed exhibit unique solvation behavior due to peak(s) in the electrical susceptibility caused by the ion shell structure; redox potential shifts for BMIM/BF4 are of more modest order ∼0.1 V, but may be larger for other ionic liquids that approach molten salt behavior.

A Gaussian field approach to the solvation of spherical ions in electrolyte solutions

In this work, the electrostatic response of an electrolyte solution to a spherical ion is studied with a Gaussian field theory. In order to capture the ionic correlation effect in concentrated solutions, the bulk dielectric response function is described by a two-Yukawa response function. The modified response function of the solution is solved analytically in the spherical geometry, from which the induced charge density and the electrostatic energy are also derived analytically. Comparisons with results for small ions in electrolyte solutions from the hyper-netted chain theory demonstrate the validity of the Gaussian field theory.

A cavity formation energy formula for hard spheres in simple electrolyte solutions

A formula for cavity formation energy of a hard sphere in restricted primitive electrolyte solutions is derived based on the integral equation theory. Specifically, the contact values of radial distribution functions between the hard sphere and the ionic species, determined analytically from the first-order mean spherical approximation theory, are used to evaluate the cavity formation energy. In the large solute-size limit, the scaling relation of the cavity formation energy further leads to an analytical expression for the surface tension of the electrolyte solution near a curved interface. Our theory is applied to hard spheres immersed in restricted primitive electrolyte solutions, where the good agreement of the cavity formation energy with the hyper-netted chain theory demonstrates the accuracy of our theory.

Reorganization energy of electron transfer

The theory of electron transfer reactions establishes the conceptual foundation for redox solution chemistry, electrochemistry, and bioenergetics. Electron and proton transfer across the cellular membrane provide all energy of life gained through natural photosynthesis and mitochondrial respiration. Rates of biological charge transfer set kinetic bottlenecks for biological energy storage. The main system-specific parameter determining the activation barrier for a single electron-transfer hop is the reorganization energy of the medium. Both harvesting of light energy in natural and artificial photosynthesis and efficient electron transport in biological energy chains require reduction of the reorganization energy to allow fast transitions. This review article discusses mechanisms by which small values of the reorganization energy are achieved in protein electron transfer and how similar mechanisms can operate in other media, such as nonpolar and ionic liquids. One of the major mechanisms of reorganization energy reduction is through non-Gibbsian (nonergodic) sampling of the medium configurations on the reaction time. A number of alternative mechanisms, such as electrowetting of active sites of proteins, give rise to non-parabolic free energy surfaces of electron transfer. These mechanisms, and nonequilibrium population of donor-acceptor vibrations, lead to a universal phenomenology of separation between the Stokes shift and variance reorganization energies of electron transfer.

Reorganization Energy of Electron Transfer in Ionic Liquids

Bandshape analysis of charge-transfer optical bands in room-temperature ionic liquids (ILs) was performed to extract the reorganization energy of electron transfer. Remarkably, the reorganization energies in ILs are close to those in cyclohexane. This result runs against common wisdom in the field since conducting ILs, which are characterized by an infinite static dielectric constant, and nonpolar cyclohexane fall to the opposite ends of the polarity scale based on their dielectric constants. Theoretical calculations employing structure factors of ILs from molecular dynamics simulations support the low values of the reorganization energy. Standard dielectric arguments do not apply to solvation in ILs, and nonergodic reorganization energies are required for a quantitative analysis.

Electron-Induced Proton Transfer

The pathway of activationless proton transfer induced by an electron-transfer reaction is studied theoretically. Long-range electron transfer produces highly nonequilibrium medium polarization that can drive proton transfer through an activationless transition during the process of thermalization, dynamically altering the screening of the electron-proton Coulomb interaction by the medium. The cross electron-proton reorganization energy is the main energy parameter of the theory, which exceeds in magnitude the proton-transfer reorganization energy roughly by the ratio of the electron-transfer to proton-transfer distance. This parameter, which can be either positive or negative, is related to the difference in pK_a values in two electron-transfer states. The relaxation time of the medium is on the (sub)picosecond time scale, which establishes the characteristic time for activationless proton transfer. Microscopic calculations predict substantial retardation of the collective relaxation dynamics compared to the continuum estimates due to the phenomenology analogous to de Gennes narrowing. Nonequilibrium medium configuration promoting proton transfer can be induced by either thermal or photoinduced charge transfer.

The Solvent Effect on H2 O2 Generation at Room Temperature Ionic Liquid|Water Interface

H₂ O₂ is a versatile chemical and can be generated by the oxygen reduction reaction (ORR) in proton donor solution in molecular solvents or room temperature ionic liquids (IL). We investigated this reaction at interfaces formed by eleven hydrophobic ILs and acidic aqueous solution as a proton source with decamethylferrocene (DMFc) as an electron donor. H₂ O₂ is generated in colorimetrically detectable amounts in biphasic systems formed by alkyl imidazolium hexafluorophosphate or tetraalkylammonium bis(trifluoromethylsulfonyl)imide ionic liquids. H₂ O₂ fluxes were estimated close to liquid|liquid interface by scanning electrochemical microscopy (SECM). Contrary to the interfaces formed by hydrophobic electrolyte solution in a molecular solvent, H₂ O₂ generation is followed by cation expulsion to the aqueous phase. Weak correlation between the H₂ O₂ flux and the difference between DMFc/DMFc⁺ redox potential and 2 electron ORR standard potential indicates kinetic control of the reaction.

Understanding the Photoexcitation of Room Temperature Ionic Liquids

Photoexcitation of (neat) room temperature ionic liquids (RTILs) leads to the observation of transient species that are reminiscent of the composition of the RTILs themselves. In this minireview, we summarize state-of-the-art in the understanding of the underlying elementary processes. By varying the anion or cation, one aim is to generally predict radiation-induced chemistry and physics of RTILs. One major task is to address the fate of excess electrons (and holes) after photoexcitation, which implies an overview of various formation mechanisms considering structural and dynamical aspects. Therefore, transient studies on time scales from femtoseconds to microseconds can greatly help to elucidate the most relevant steps after photoexcitation. Sometimes, radiation may eventually result in destruction of the RTILs making photostability another important issue to be discussed. Finally, characteristic heterogeneities can be associated with specific physicochemical properties. Influencing these properties by adding conventional solvents, like water, can open a wide field of application, which is briefly summarized.

Photodetachment and Electron Dynamics in 1-Butyl-1-methyl-pyrrolidinium Dicyanamide

The ultrafast transient absorption spectrum of 1-butyl-1-methyl-pyrrolidinium dicyanamide, [Pyr_1,4⁺][DCA^-], was measured in the visible and near-infrared (IR) spectral regions. Excitation of the liquid at 4.6 eV created initially delocalized and highly reactive electrons that either geminately recombined (69%) or localized onto a cavity with a time constant of ∼300 fs. Electron localization was reflected in the evolution of the TA spectrum and the time-dependent loss of reactivity with a dichloromethane quencher. The delocalized initial state and spectrum of the free electrons were consistent with computational predictions by Xu and Margulis [ J. Phys. Chem. B, 2015, 119, 532-542] on excess electrons in [Pyr_1,4⁺][DCA^-]. The computational study considered two possible localization mechanisms for excess electrons, localization on ions, and localization on cavities. In the case of photogenerated electron-hole pairs, the results presented here demonstrate localization to cavities as the dominant channel. Following localization onto a cavity, the free electrons underwent solvation and loss of reactivity with the quencher with rates that slowed in time. The dynamics were similar to an analogous prior study on the related liquid [Pyr_1,x⁺][NTf₂^-]. One significant difference was the larger yield of free electrons from photoexcitation of [Pyr_1,4⁺][DCA^-]. This was found to primarily reflect more efficient localization onto cavities rather than a slower geminate recombination rate.