Thiocyanate as a Local Probe of Ultrafast Structure and Dynamics in Imidazolium-Based Ionic Liquids: Water-Induced Heterogeneity and Cation-Induced Ion Pairing

Uncovering the binding nature of thiocyanate in contact ion pairs with lithium ions

Ion pair formation is a fundamental molecular process that occurs in a wide variety of systems, including electrolytes, biological systems, and materials. In solution, the thiocyanate (SCN-) anion interacts with cations to form contact ion pairs (CIPs). Due to its ambidentate nature, thiocyanate can bind through either its sulfur or nitrogen atoms, depending on the solvent. This study focuses on the binding nature of thiocyanate with lithium ions as a function of the solvents using FTIR, 2D infrared spectroscopy (2DIR) spectroscopies, and theoretical calculations. The study reveals that the SCN- binding mode (S or N end) in CIPs can be identified through 2DIR spectroscopy but not by linear IR spectroscopy. Linear IR spectroscopy shows that the CN stretch frequencies are too close to one another to separate N- and S-bound CIPs. Moreover, the IR spectrum shows that the S-C stretch presents different frequencies for the salt in different solvents, but it is related to the anion speciation rather than to its binding mode. A similar trend is observed for the anion bend. 2DIR spectra show different dynamics for N-bound and S-bound thiocyanate. In particular, the frequency-frequency correlation function (FFCF) dynamics extracted from the 2DIR spectra have a single picosecond exponential decay for N-bound thiocyanate and a biexponential decay for S-bound thiocyanate, consistent with the binding mode of the anion. Finally, it is also observed that the binding mode also affects the line shape parameters, probably due to the different molecular mechanisms of the FFCF for N- and S-bound CIPs.

Transitioning from Regular Electrolytes to Solvate Ionic Liquids to High-Concentration Electrolytes: Changes in Transport Properties and Ionic Speciation

Glyme-based lithium-ion electrolytes have received considerable attention from the scientific community due to their improved safety, as well as electrochemical and thermal stability over carbonate-based electrolytes. However, these electrolytes suffer from major drawbacks such as high viscosities. To overcome the challenges that hinder their full potential, the molecular description of glyme-based lithium electrolytes in the high-concentration regime, particularly in the solvate ionic liquid (SIL) and high-concentration electrolyte (HCE) regimes, is needed. In this study, model glyme-based electrolytes based on a lithium thiocyanate and either tetraglyme (G4) or a mixture of monoglyme (G1) and diglyme (G2) were investigated as a function of the solvent-to-lithium ratio using linear and nonlinear IR spectroscopies, in combination with ab initio computations as well as electrochemical methods . The transport properties reveal enhanced ionicities in the HCE and SIL regimes ([O]/[Li] ≤ 5) compared to the regular electrolytes (REs, with [O]/[Li] > 5) in both pure (G4) and mixed (G1:G2) glymes. IR and ab initio computations relate these larger ionicities to the higher concentration of charged aggregates in the HCE and SIL electrolytes ([O]/[Li] ≤ 5). Moreover, it was observed that the use of mixed glymes appears to have a minimal effect on the transport properties of REs but exhibits deleterious effects on SILs. Overall, the results provide a molecular framework for describing the local structure of lithium glyme-based electrolytes and demonstrate the key role that the nature of glyme solvation plays in the molecular structure and consequently the macroscopic properties of the Li-glyme SILs, HCEs, and REs.

Chain-like Structures Facilitate Li+ Transport in Concentrated Aqueous Electrolytes: Insights from Ultrafast Infrared Spectroscopy and Molecular Dynamics Simulations

Highly concentrated aqueous electrolytes have attracted attention due to their unique applications in lithium ion batteries (LIBs). However, the solvation structure and transport mechanism of Li+ cations at concentrated concentrations remain largely unexplored. To address this gap in knowledge, we employ ultrafast infrared spectroscopy and molecular dynamics (MD) simulations to reveal the dynamic and spatial structural heterogeneity in aqueous lithium chloride (LiCl) solutions. The coupling between the reorientation dynamics of the extrinsic probe and the macroscopic viscosity in aqueous LiCl solutions was analyzed using the Stokes-Einstein-Debye (SED) equations. MD simulations reveal that the Cl- and Li+ form chain-like structures through electrostatic interactions, supporting the vehicular migration of Li+ through the chain-like structure. The concentration dependent conductivity of the LiCl solution is well reproduced, where Li(H2O)2+ and Li(H2O)3+ are the dominant species that contribute to the conduction of Li+. This study is expected to establish correlations between ion pair structures and macroscopic properties.

Dicyanamide Anion Reports on Water Induced Local Structural and Dynamic Heterogeneity in Ionic Liquid Mixtures

The effects of limited amounts (under 21.6% χ_Water) of water on 1-butyl-3-methylimidazolium tetrafluoroborate (BmimBF₄) and 1-butyl-3-methylimidazolium dicyanamide (BmimDCA) room-temperature ionic liquid (RTIL) mixtures were characterized by tracking changes in the linear and two-dimensional infrared (2D IR) vibrational features of the dicyanamide anion (DCA). Peak shifts with increasing water suggest the formation of water-associated and nonwater-associated DCA populations. Further results showed clear differences in the dynamic behavior of these different populations of DCA at low (defined here as below 2.5% χ_Water), mid (defined here as between 2.5% χ_Water and 9.6% χ_Water), and high (defined here as between 11.6% χ_Water and 21.6% χ_Water) range water concentrations. Vibrational relaxation is accelerated with increasing water content for water-associated populations of DCA, indicating water facilitates population relaxation, possibly through the provision of additional bath modes. Conversely, spectral diffusion of water-associated populations slowed dramatically with increasing water, suggesting that water drives the formation of distinct and noninterchangeable or very slowly interchangeable local solvent environments.

Molten salt electro-preparation of graphitic carbons

Graphite has been used in a wide range of applications since the discovery due to its great chemical stability, excellent electrical conductivity, availability, and ease of processing. However, the synthesis of graphite materials still remains energy-intensive as they are usually produced through a high-temperature treatment (>3000°C). Herein, we introduce a molten salt electrochemical approach utilizing carbon dioxide (CO2) or amorphous carbons as raw precursors for graphite synthesis. With the assistance of molten salts, the processes can be conducted at moderate temperatures (700-850°C). The mechanisms of the electrochemical conversion of CO2 and amorphous carbons into graphitic materials are presented. Furthermore, the factors that affect the graphitization degree of the prepared graphitic products, such as molten salt composition, working temperature, cell voltage, additives, and electrodes, are discussed. The energy storage applications of these graphitic carbons in batteries and supercapacitors are also summarized. Moreover, the energy consumption and cost estimation of the processes are reviewed, which provides perspectives on the large-scale synthesis of graphitic carbons using this molten salt electrochemical strategy.

Direct Correlation between Short-Range Vibrational Spectral Diffusion and Localized Ion-Cage Dynamics of Water-in-Salt Electrolytes

The molecular dynamics simulations of a "water-in-salt" electrolyte, lithium bis(trifluoromethyl sulfonyl) imide (LiNTf₂), with a varying concentration range of 3 to 20 m were performed to establish a direct connection between a dynamic property like the ion-cage lifetime with the short-range vibrational stretching frequency shift of the used probe, HOD. The properties reported here are compared to that obtained from experiments performed at the same concentrations. The time-series wavelet transform was adopted as a preferable mathematical tool for calculating the instantaneous fluctuating frequencies of the probe O-D stretch mode and the concentration-dependent vibrational stretch spectral signature based on the variable functions associated with a particular chemical bond derived from classical molecular dynamics trajectories. The decay time constants of frequency fluctuations and the lifetime of the ion cage (τ_IC) were estimated as a function of salt concentration. Herein, we emphasize the correlation between the slowest time constant (τ₃) of the decay of O-D stretch frequency fluctuations and the timescales associated with the lifetime of ion cages (τ_IC). The results exhibit that the existing relationships were also concentration-dependent. Therefore, this study highlights the connection between the ionic motions that regulate the overall system dynamics with the short-range vibrational frequency shift of the used probe, which was used similar to experiments. It also provides an understanding of the interionic interactions and the dynamical and spectral properties of the electrolytic mixtures. We establish a direct correlation between short-range frequency profile and localized ion-cage lifetime, which can fill the gap of understanding between viscosity, vibrational frequency, and ion-cage dynamics of electrolytes.

Molecular Simulation-Guided Spectroscopy of Imidazolium-Based Ionic Liquids and Effects of Methylation on Ion-Cage and -Pair Dynamics

Classical molecular dynamics simulations were performed to assess an atomistic interpretation of the ion-probe structural interactions in two typical ionic liquids (ILs), 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide [BMIm][NTf₂] and 1-butyl-2,3-dimethylimidazolium bis(trifluoromethylsulfonyl)imide [BDimIm][NTf₂] through computational ultrafast spectroscopy. The nitrile stretching vibrations of the thiocyanate anion, [SCN]^-, serve as the local mode of the ultrafast system dynamics within the imidazolium-based ionic liquid environment. The wavelet transform of classical trajectories determines the time-varying fluctuating frequencies and the stretch spectral signatures of SCN^- in the normalized distribution. However, computational modeling of the two-dimensional (2D) spectra from the wavelet-derived vibrational frequencies yields time evolution of the local molecular structure along with the varied time-dependent dynamics of the spectral diffusion process. We calculated the frequency-frequency correlation functions (FFCFs), time correlations associated with the ion-pair and -cage dynamics, and mean square displacements as a function of time, depicting diffusive dynamics. The calculated results based on the pair correlation functions and the distribution of atomic density suggest that the hydrogen and methylated carbon at the two-position of the imidazolium ring of [BMIm] and [BDimIm] cations, respectively, strongly interact with the probe through the N of the thiocyanate anion rather than the S atom. The center-of-mass center-of-mass (COM-COM) cation-probe radial distribution functions (RDFs) in conjunction with the site-specific structural analysis further reveal well-structured interactions of the thiocyanate ion and [BMIm]⁺ cation rather than the [BDimIm] cation. In contrast, the anion-probe COM-COM RDFs depict weak interactive associations within the vibrational probe [SCN]^- and [NTf₂]^- ions. Methylation at the two-position of the imidazolium ring predicts slower structural reorganization and breaking and reformation dynamics of the ion pairs and cages within the ionic liquid framework.

Association-Dissociation Dynamics of Ionic Electrolytes in Low Dielectric Medium

Ionic electrolytes are known to form various complexes which exist in dynamic equilibrium in a low dielectric medium. However, structural characterization of these complexes has always posed a great challenge to the scientific community. An additional challenge is the estimation of the dynamic association-dissociation time scales (lifetime of the complexes), which are key to the fundamental understanding of ion transport. In this work, we have used a combination of infrared absorption spectroscopy, two-dimensional infrared spectroscopy, molecular dynamics simulations, and density functional theory calculations to characterize the various ion complexes formed by the thiocyanate-based ionic electrolytes as a function of different cations in a low dielectric medium. Our results demonstrate that thiocyanate is an excellent vibrational reporter of the heterogeneous ion complexes undergoing association-dissociation dynamics. We find that the ionic electrolytes exist as contact ion pairs, dimers, and clusters in a low dielectric medium. The relative ratios of the various ion complexes are sensitive to the cations. In addition to the interactions between the thiocyanate anion and the countercation, the solute-solvent interactions drive the dynamic equilibrium. We have estimated the association-dissociation dynamics time scales from two-dimensional infrared spectroscopy. The exchange time scale involving the cluster is faster than that between a dimer and an ion pair. Moreover, we find that the dynamic equilibrium between the cluster and another ion complex is correlated to the solvent fluctuations.

Least-Squares Fitting of Multidimensional Spectra to Kubo Line-Shape Models

We report a comprehensive study of the efficacy of least-squares fitting of multidimensional spectra to generalized Kubo line-shape models and introduce a novel least-squares fitting metric, termed the scale invariant gradient norm (SIGN), that enables a highly reliable and versatile algorithm. The precision of dephasing parameters is between 8× and 50× better for nonlinear model fitting compared to that for the centerline-slope (CLS) method, which effectively increases data acquisition efficiency by 1–2 orders of magnitude. Whereas the CLS method requires sequential fitting of both the nonlinear and linear spectra, our model fitting algorithm only requires nonlinear spectra but accurately predicts the linear spectrum. We show an experimental example in which the CLS time constants differ by 60% for independent measurements of the same system, while the Kubo time constants differ by only 10% for model fitting. This suggests that model fitting is a far more robust method of measuring spectral diffusion than the CLS method, which is more susceptible to structured residual signals that are not removable by pure solvent subtraction. Statistical analysis of the CLS method reveals a fundamental oversight in accounting for the propagation of uncertainty by Kubo time constants in the process of fitting to the linear absorption spectrum. A standalone desktop app and source code for the least-squares fitting algorithm are freely available, with example line-shape models and data. We have written the MATLAB source code in a generic framework where users may supply custom line-shape models. Using this application, a standard desktop fits a 12-parameter generalized Kubo model to a 10⁶ data-point spectrum in a few minutes.

Structural and dynamical changes observed when transitioning from an ionic liquid to a deep eutectic solvent

The microscopic molecular structure and dynamics of a new deep eutectic solvent (DES) composed of an ionic liquid (1-hexyl-3-methylimidazolium chloride) and an amide (trifluoroacetamide) at various molar ratios were investigated using linear and non-linear infrared spectroscopy with a vibrational probe. The use of the ionic liquid allows us to investigate the changes that the system undergoes with the addition of the amide or, equivalently, the changes from an ionic liquid to a DES. Our studies revealed that the vibrational probe in the DES senses a very similar local environment irrespective of the cation chemical structure. In addition, the amide also appears to perceive the same molecular environment. The concentration dependence studies also showed that the amide changes from being isolated from other amides in the ionic liquid environment to an environment where the amide-amide interactions are favored. In the case of the vibrational probe, the addition of the amide produced significant changes in the slow dynamics associated with the making and breaking of the ionic cages but did not affect the rattling-in-cage motions perceived by it. Furthermore, the concentration dependence of slow dynamics showed two regimes which are linked to the changes in the overall structure of the solution. These observations are interpreted in the context of a nanoscopic heterogeneous environment in the DES which, according to the observed dynamical regimes, appears at very large concentrations of the amide (molar ratio of greater than 1:1) since for lower amide molar ratios, the amide appears to be not segregated from the ionic liquid. This proposed molecular picture is supported by small angle x-ray scattering experiments.

Dynamics of Ionic Liquid through Intrinsic Vibrational Probes Using the Dispersion-Corrected DFT Functionals

First principles molecular dynamics simulations have been utilized to study the spectral properties of the protic ionic liquid, methylammonium formate (MAF). All simulations were performed using density functional theory (DFT) and various van der Waals-corrected exchange-correlation functionals. We calculated the vibrational stretch frequency distributions, determined the time-frequency correlations of the intrinsic vibrational probes, the N-H and C-O modes in MAF, and the frequency-structure correlations. We also estimated the average hydrogen-bond lifetimes and orientation dynamics to capture the ultrafast spectral response. The spectroscopic signature of the N-H stretching vibrations using the Becke-Lee-Yang-Parr (BLYP) and Perdew-Burke-Ernzerhof (PBE) functionals displays a spectral shift in the lower frequency side, suggesting stronger hydrogen-bonding interactions represented by the gradient approximation functionals than the van der Waals (vdW)-corrected simulations. The carboxylate frequency profiles with the dispersion-corrected representations are almost similar without a significant difference in the normalized distributions. Besides, the COO stretching frequencies at the peak maxima positions of the PBE functionals exhibit a lesser deviation from the experimental data. Spectral diffusion dynamics of the intrinsic vibrational probes on the cationic and anionic sites of the ionic liquid proceed through a short time relaxation of the intact hydrogen bonds followed by an intermediate time constant and a longer time decay indicating the switchover of hydrogen bonds. Dispersion-corrected atom-centered one-electron potential (DCACP) correction added to the BLYP system slows down the picosecond time scales of frequency correlation and the time constants of rotational motion, lengthening the overall system dynamics. The observed trends in the time-dependent decays of frequency fluctuations and the orientation autocorrelation functions correlate with the structural interactions in liquid MAF and hydrogen-bond dynamics. In this study, we examine the predictions made by different density functional treatments comparing the results of the uncorrected BLYP and PBE representations with the semiempirical vdW methods of Grimme and matching our calculated data with the experimental observations.

An ultrafast vibrational study of dynamical heterogeneity in the protic ionic liquid ethyl-ammonium nitrate. I. Room temperature dynamics

Using ultrafast two-dimensional infrared spectroscopy (2D-IR), a vibrational probe (thiocyanate, SCN^-) was used to investigate the hydrogen bonding network of the protic ionic liquid ethyl-ammonium nitrate (EAN) in comparison to H₂O. The 2D-IR experiments were performed in both parallel (⟨ZZZZ⟩) and perpendicular (⟨ZZXX⟩) polarizations at room temperature. In EAN, the non-Gaussian lineshape in the FTIR spectrum of SCN^- suggests two sub-ensembles. Vibrational relaxation rates extracted from the 2D-IR spectra provide evidence of the dynamical differences between the two sub-ensembles. We support the interpretation of two sub-ensembles with response function simulations of two overlapping bands with different vibrational relaxation rates and, otherwise, similar dynamics. The measured rates for spectral diffusion depend on polarization, indicating reorientation-induced spectral diffusion (RISD). A model of restricted molecular rotation (wobbling in a cone) fully describes the observed spectral diffusion in EAN. In H₂O, both RISD and structural spectral diffusion contribute with similar timescales. This complete characterization of the dynamics at room temperature provides the basis for the temperature-dependent measurements in Paper II of this series.

Orientational Pair Correlations and Local Structure of Benzonitrile from Molecular Dynamics Simulations with Comparisons to Experiments

We present an experimentally parametrized molecular dynamics study of single-molecule and collective orientational relaxation in neat benzonitrile through the analysis of the reorientational anisotropy and polarizability anisotropy time correlation function (PA-TCF). The simulations show that the PA-TCF is dominated by collective reorientation after 20 ps. Collective reorientation is found to be slower than single-molecule reorientation by a factor of 1.67, consistent with recent experiments. The simulations provide direct evidence of local antiparallel benzonitrile configurations. These structures, which have been the center of some debate, are responsible for the slower rate of collective versus single-molecule reorientation in the liquid. Further structural analysis indicates that significant Coulombic interactions between the nitrile group and hydrogen atoms on adjacent molecules play a role in the formation of the antiparallel structures. The single-molecule dynamics reflected in the anisotropy are complex and consist of a ballistic regime, restricted angular diffusion, and spatially anisotropic free diffusion. The principal components of the rotational diffusion tensor are independently obtained and shown to reproduce the free diffusion regime of the anisotropy for each principal axis according to the predictions of a previous theory.

Specific Host-Guest Interactions in the Crown Ether Complexes with K+ and NH4+ Revealed from the Vibrational Relaxation Dynamics of the Counteranion

The specific host-guest interactions in the corresponding complexes of K⁺ and NH₄⁺ with typical crown ethers were investigated by using FTIR and ultrafast IR spectroscopies. The counteranions, i.e., SCN^-, were employed as a local vibrational probe to report the structural dynamics of the complexation. It was found that the vibrational relaxation dynamics of the SCN^- was strongly affected by the cations confined in the cavities of the crown ethers. The time constant of the vibrational population decay of SCN^- in the complex of NH₄⁺ with the 18-crown-6 was determined to be 6 ± 2 ps, which is ∼30 times faster than that in the complex of K⁺ with the crown ethers. Control experiments showed that the vibrational population decay of SCN^- depended on the size of the cavities of the crown ethers. A theoretical calculation further indicated that the nitrogen atom of SCN^- showed preferential coordination to the K⁺ ions hosted by the crown ethers, while the NH₄⁺ can form hydrogen bonds with the oxygen atoms in the studied crown ethers. The geometric constraints formed in the complex of crown ethers can cause a specific interaction between the NH₄⁺ and SCN^-, which can facilitate the intermolecular vibrational energy redistribution of the SCN^-.

Strong Propensity of Ionic Liquids in Their Aqueous Solutions for an Organic-Modified Metal Surface

Understanding ionic structure and electrostatic environments near a surface has both fundamental and practical value. In electrochemistry, especially when room temperature ionic liquids (ILs) are involved, the complex ionic structure near the interface is expected to crucially influence reactions. Here we report evidence that even in dilute aqueous solutions of several ILs, the ions aggregate near the surface in ways that are qualitatively different from simple electrolytes. We have used a vibrational probe molecule, 4-mercaptobenzonitrile (MBN), tethered to a metal surface to monitor the behavior of the ionic layers. The characteristic nitrile vibrational frequency of this molecule has distinct values in the presence of pure water (∼2232 cm^-1) and pure IL (for example, ∼2226 cm^-1 for ethylmethylimidazolium tetrafluoroborate, [EMIM][BF₄]). This difference reflects the local electrostatic field and the hydrogen-bonding variations between these two limiting cases. We tracked this frequency shift as a function of IL concentration in water all the way from pure water to pure IL. We report two important findings. First, only one nitrile peak is observed for the entire concentration range, indicating that at least on the length scale of the probe molecule water and ILs do not phase separate within the interface, and no heterogeneously distinct electrostatic environments are formed. Second, and more importantly, we find that even up to a significant mole fraction of bulk water (x ∼ 0.95), the nitrile frequency does not change from that indicative of a pure IL for [EMIM][BF₄], indicating preferential aggregation of the ions near the surface. Because this behavior is very similar to surfactants, we chose an imidazolium cation with a longer side chain which resulted in behavior expected from a surfactant, with a preferential layer of the ions on the surface even in dilute water solutions (x ∼ 0.995). This observation indicates that even those ILs that are not nominally categorized as surfactants have a strong tendency to aggregate at the surface. Because ILs serve as electrolytes in a range of electrochemical reactions, including those requiring water, our results are likely useful for mechanistic understanding and tuning of such reactions.

Microstructural and Dynamical Heterogeneities in Ionic Liquids

Ionic liquids (ILs) are a special category of molten salts solely composed of ions with varied molecular symmetry and charge delocalization. The versatility in combining varied cation-anion moieties and in functionalizing ions with different atoms and molecular groups contributes to their peculiar interactions ranging from weak isotropic associations to strong, specific, and anisotropic forces. A delicate interplay among intra- and intermolecular interactions facilitates the formation of heterogeneous microstructures and liquid morphologies, which further contributes to their striking dynamical properties. Microstructural and dynamical heterogeneities of ILs lead to their multifaceted properties described by an inherent designer feature, which makes ILs important candidates for novel solvents, electrolytes, and functional materials in academia and industrial applications. Due to a massive number of combinations of ion pairs with ion species having distinct molecular structures and IL mixtures containing varied molecular solvents, a comprehensive understanding of their hierarchical structural and dynamical quantities is of great significance for a rational selection of ILs with appropriate properties and thereafter advancing their macroscopic functionalities in applications. In this review, we comprehensively trace recent advances in understanding delicate interplay of strong and weak interactions that underpin their complex phase behaviors with a particular emphasis on understanding heterogeneous microstructures and dynamics of ILs in bulk liquids, in mixtures with cosolvents, and in interfacial regions.

Assessing the Location of Ionic and Molecular Solutes in a Molecularly Heterogeneous and Nonionic Deep Eutectic Solvent

Deep eutectic solvents (DES) are emerging sustainable designer solvents viewed as greener and better alternatives to ionic liquids. Nonionic DESs possess unique properties such as viscosity and hydrophobicity that make them desirable in microextraction applications such as oil-spill remediation. This work builds upon a nonionic DES, NMA–LA DES, previously designed by our group. The NMA–LA DES presents a rich nanoscopic morphology that could be used to allocate solutes of different polarities. In this work, the possibility of solvating different solutes within the nanoscopically heterogeneous molecular structure of the NMA–LA DES is investigated using ionic and molecular solutes. In particular, the localized vibrational transitions in these solutes are used as reporters of the DES molecular structure via vibrational spectroscopy. The FTIR and 2DIR data suggest that the ionic solute is confined in a polar and continuous domain formed by NMA, clearly sensing the direct effect of the change in NMA concentration. In the case of the molecular nonionic and polar solute, the data indicates that the solute resides in the interface between the polar and nonpolar domains. Finally, the results for the nonpolar and nonionic solute (W(CO)₆) are unexpected and less conclusive. Contrary to its polarity, the data suggest that the W(CO)₆ resides within the NMA polar domain of the DES, probably by inducing a domain restructuring in the solvent. However, the data are not conclusive enough to discard the possibility that the restructuring comprises not only the polar domain but also the interface. Overall, our results demonstrate that the NMA–LA DES has nanoscopic domains with affinity to particular molecular properties, such as polarity. Thus, the presented results have a direct implication to separation science.

Experimental observation of nanophase segregation in aqueous salt solutions around the predicted liquid-liquid transition in water

The liquid-liquid transition in supercooled liquid water, predicted to occur around 220 K, is controversial due to the difficulty of studying it caused by competition from ice crystallization (the so-called "no man's land"). In aqueous solutions, it has been predicted to give rise to phase separation on a nanometer scale between a solute-rich high-density phase and a water-rich low-density phase. Here we report direct experimental evidence for the formation of a nanosegregated phase in eutectic aqueous solutions of LiCl and LiSCN where the presence of crystalline water can be experimentally excluded. Femtosecond infrared and Raman spectroscopies are used to determine the temperature-dependent structuring of water, the solvation of the SCN^- anion, and the size of the phase segregated domains.

Ultrafast vibrational dynamics of a trigonal planar anionic probe in ionic liquids (ILs): A two-dimensional infrared (2DIR) spectroscopic investigation

A major impediment limiting the widespread application of ionic liquids (ILs) is their high shear viscosity. Incorporation of a tricyanomethanide (TCM^-) anion in ILs leads to low shear viscosity and improvement of several characteristics suitable for large scale applications. However, properties including interactions of TCM^- with the local environment and dynamics of TCM^- have not been thoroughly investigated. Herein, we have studied the ultrafast dynamics of TCM^- in several imidazolium ILs using linear IR and two-dimensional infrared spectroscopy techniques. The spectral diffusion dynamics of the CN stretching modes of TCM^- in all ILs exhibit a nonexponential behavior with a short time component of ∼2 ps and a long time component spanning ∼9 ps to 14 ps. The TCM^- vibrational probe reports a significantly faster relaxation of ILs compared to those observed previously using linear vibrational probes, such as thiocyanate and selenocyanate. Our results indicate a rapid relaxation of the local ion-cage structure embedding the vibrational probe in the ILs. The faster relaxation suggests that the lifetime of the local ion-cage structure decreases in the presence of TCM^- in the ILs. Linear IR spectroscopic results show that the hydrogen-bonding interaction between TCM^- and imidazolium cations in ILs is much weaker. Shorter ion-cage lifetimes together with weaker hydrogen-bonding interactions account for the low shear viscosity of TCM^- based ILs compared to commonly used ILs. In addition, this study demonstrates that TCM^- can be used as a potential vibrational reporter to study the structure and dynamics of ILs and other molecular systems.

Vibrational Relaxation Dynamics of a Semiconductor Copper(I) Thiocyanate (CuSCN) Film as a Hole-Transporting Layer

The semiconductor CuSCN film, which is typically used as the hole-transporting layer (HTL) in solar cell studies, has been investigated by Fourier transform infrared (FTIR) spectroscopy and ultrafast transient infrared (IR) spectroscopy. A sharp peak at 2175 cm^-1 corresponding to the CN vibrational stretching mode in CuSCN was observed, and the peak frequency remained unchanged by varying the thickness of the CuSCN thin film. Vibrational relaxation measurements showed that the 0-1 and 1-2 transitions of CN stretching can be observed at 2175 and 2140 cm^-1, respectively. The heat-induced absorption and bleaching peaks (2167 and 2175 cm^-1) can be clearly seen at a waiting time of 40 ps. The vibrational relaxation of the CN stretching mode determined from the 1-2 transition exhibited a biexponential decay with time constants of 7.4 ± 0.5 (90%) and 158 ± 50 ps (10%). Importantly, the abnormal anisotropy decay of the CN stretching mode in the CuSCN thin film was also observed for the first time. A detailed analysis showed that the distinct anisotropy decay curve could be described using a triexponential decay function, which was explained by three different processes: resonance energy transfer (∼8 ps), a thermalization process (∼40 ps), and molecular rotation (∼150 ps). The time scale of the thermalization process caused by the vibrational relaxation in CuSCN is at a time scale of 40 ps, which is important for us to understand the thermally activated charge-transport property of the CuSCN film employed as the HTL. Further UV pump-IR probe measurement revealed that the carrier scattering and relaxation processes in the CuSCN film are strongly associated with the vibrational excitation and relaxation dynamics of the CN stretching mode. It is expected that the fundamental understanding of the vibrational relaxation dynamics of the CuSCN thin film should provide helpful insight to elucidate its role as the HTL in solar cell studies at the molecular level.

Dynamics and Vibrational Spectroscopy of Alcohols in Ionic Liquids: Methanol and Ethanol

The structure, dynamics, and vibrational spectroscopy of dilute HOD, methanol, and ethanol in the 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide, [emim][NTf₂], ionic liquid (IL) are investigated with molecular dynamics (MD) simulations. The structure of the ILs around the solutes is qualitatively similar, where the OD bond of the deuterated alcohols donates an interaction to an [NTf₂] anion and the [emim] cations interact with the oxygen atom of the OD group. The slowest time scale for the reorientational dynamics of the OD bond varied considerably for HOD, methanol, and ethanol (27, 71, and 87 ps, respectively). In contrast, the slowest time scales for spectral diffusion of the OD vibrational frequency were 11 ps for each of the three solutes, which indicates that the dynamics of the IL is relatively unchanged by the presence of the alcohols at dilute concentration. The theoretical results for the reorientational and spectral diffusion dynamics compare favorably with prior two-dimensional infrared (2D IR) spectroscopic measurements.

C-D Vibration at C2 Position of Imidazolium Cation as a Probe of the Ionic Liquid Microenvironment

Unlike molecular solvents, imidazolium-based ionic liquids are entirely made of ions with spatial heterogeneity. There is a need for spectroscopic probes that can assess the microenvironment near the cations of these complex liquids. In this manuscript, we describe simple chemical procedures to label the C2 position of imidazolium cation with a C-D vibrational probe and show, through linear and nonlinear vibrational spectroscopies, that this C-D stretching mode can be a useful analytical tool to assess both the solvent microenvironment and solute-solvent interactions in imidazolium-based ionic liquids from the cation point of view. It is expected that this C-D vibration probe on the cation will lead to the development of innovative experimental strategies that can provide a better understanding of such ionic liquids.

Ultrafast Hydrogen Bond Exchanging between Water and Anions in Concentrated Ionic Liquid Aqueous Solutions

The mixtures of 1-butyl-3-methylimidazolium tetrafluoroborate ([BMIM]BF₄) ionic liquids (ILs) and water as a function of IL concentrations have been investigated by Fourier transform infrared (FTIR) spectroscopy and ultrafast two-dimensional IR (2D IR) spectroscopy. FTIR spectra of the mixtures resolve two different types of water species, one interacting with the BF₄^- anions and the other associated with bulklike water molecules. These two water species are in a dynamic equilibrium through forming different hydrogen bonding configurations which are separated by more than 100 cm^-1 in the IR spectra. The structural dynamics of the IL mixtures are further revealed by monitoring the vibrational relaxation dynamics of the OD stretching group of interfacial water molecules hydrogen bonded to BF₄^- anions. With the increase of the IL bulk concentration, vibrational population and rotational dynamics of the interfacial water molecules can be described by a biexponential decay function and are strongly dependent on the IL concentrations. Furthermore, the ultrafast hydrogen bond exchanging between water and BF₄^- anions in the ILs are also measured using 2D IR spectroscopy. The average hydrogen bond exchanging rate is determined to be 19 ± 4 ps, which is around 3 times slower than that in the NaBF₄ electrolyte aqueous solution. The much slower hydrogen bond exchanging rate indicates that the local structure of ILs and water molecules are strongly mediated by the steric effect of the cationic group in the ILs, which is proposed to be responsible for the formation of the heterogeneous structure in the IL mixtures. By using SCN^- as the anionic probe, the structural inhomogeneity in the IL solutions can be confirmed from the distinct rotational dynamics of the SCN^-, which is segregated from the rotational dynamics of water molecules in the IL mixtures.

Modeling Carbon Dioxide Vibrational Frequencies in Ionic Liquids: IV. Temperature Dependence

In previous papers in the series, the vibrational spectroscopy of CO₂ in ionic liquids (ILs) was investigated at ambient conditions. Here, we extend these studies to understand the temperature dependence of the structure, dynamics, and thermodynamics of CO₂ in the 1-butyl-3-methylimidazolium hexafluorophosphate, [bmim][PF₆], IL. Using spectroscopic mapping techniques, the infrared absorption spectrum of the CO₂ asymmetric stretch mode is simulated at a number of temperatures, and the results are found to be consistent with similar experimental studies. Structural correlation functions are used to reveal the thermodynamics of complete CO₂ solvent cage breakdown. The enthalpy and entropy of activation for solvent cage reorganization are found to be 6.9 and 7.6 (kcal/mol)/K, respectively, and these values are similar to the those for spectral, orientational, and translational diffusion. Caging times for CO₂ are calculated, and it is shown that the short time dynamics of CO₂ are unaffected by temperature, even though the long-time dynamics are highly sensitive to temperature.

Dynamics of Dicyanamide in Ionic Liquids is Dominated by Local Interactions

The dynamics of probe molecules is commonly used to investigate the structural dynamics of room-temperature ionic liquids; however, the extent to which this dynamics reflects the dynamics of the ionic liquids or is probe specific has remained debated. Here, we explore to what extent the vibrational and rotational dynamics of the dicyanamide anion, a common ionic liquid anion, correlates with the structural relaxation of ionic liquids. We use polarization-resolved, ultrafast infrared spectroscopy to probe the temperature- and probe-concentration-dependent dynamics of samples with small amounts of 1-ethyl-3-methylimidazolium ([emim]⁺) dicyanamide ([DCA]⁻) dissolved in four [emim]⁺-based ionic liquids with tetrafluoroborate ([BF₄]⁻), bis(trifluoromethylsulfonyl)imide ([NTf₂]⁻), ethylsufate ([EtSO₄]⁻), and triflate ([OTf]⁻) as anions. The transient spectra after broad-band excitation at 2000–2300 cm^–1, resonant with the symmetric and antisymmetric C≡N stretching vibrations, initially contain oscillatory signatures due to the vibrational coherence between both modes. Vibrational population relaxation occurs on two distinct time scales, ∼6–7 and ∼15–20 ps. The vibrational dynamics is rather insensitive to the details of the ionic liquid anion and temperature, except for the slow vibrational relaxation component. The decay of the excitation anisotropy, a measure of the rotational dynamics of [DCA]⁻, markedly depends on temperature, and the obtained decay time exhibits an activation energy of ∼15–21 kJ/mol. Remarkably, neither the rotation time nor the activation energy can be simply explained by the variation of the macroscopic viscosity. Hence, our results suggest that the dynamics of dicyanamide is only in part representative of the ionic liquid structural dynamics. Rather, the dynamics of the probe anion seems to be determined by the specific interaction of [DCA]⁻ with the ionic liquid’s ions for the class of [emim]⁺-based ionic liquids studied here.

Counterion Effect on Vibrational Relaxation and the Rotational Dynamics of Interfacial Water and an Anionic Vibrational Probe in the Confined Reverse Micelles Environment

Vibrational relaxation and the rotational dynamics of water molecules encapsulated in reverse micelles (RMs) have been investigated by ultrafast infrared (IR) spectroscopy and two-dimensional IR (2D IR) spectroscopy. By changing the counterion of the hydrophilic headgroup in the RMs formed by Aerosol-OT (AOT) from Na⁺ to K⁺, Cs⁺ and Ca²⁺, we could determine the specific counterion effects on the rotational dynamics of water molecules. The orientational relaxation time constant of water decreases in the order Ca²⁺ > Na⁺ > K⁺ > Cs⁺. The SCN^- anionic probe and counterion can form ion pairs at the interfacial region of the RMs. The rotational dynamics of SCN^- anion significantly decreases because of the synergistic effects of confinement and the surface interactions in the interfacial region of the RMs. The results can provide a new understanding of the cationic Hofmeister effect at the molecular level observed in biological studies.

Modeling Carbon Dioxide Vibrational Frequencies in Ionic Liquids: III. Dynamics and Spectroscopy

In recent years, interest in carbon capture and sequestration has led to numerous investigations of the ability of ionic liquids to act as recyclable CO₂-sorbent materials. Herein, we investigate the structure and dynamics of a model physisorbing ionic liquid, 1-butyl-3-methylimidazolium hexafluorophosphate ([C₄C₁Im][PF₆]), from the perspective of CO₂ using two-dimensional (2D) IR spectroscopy and molecular dynamics simulations. A direct comparison of experimentally measured and calculated 2D IR line shapes confirms the validity of the simulations and spectroscopic calculations. Taken together, the simulations and experiments reveal new insights into the interactions of a CO₂ solute with the surrounding ionic liquid and how these interactions manifest in the 2D IR spectra. In particular, higher CO₂ asymmetric stretch vibrational frequencies are associated with softer, less populated solvent cages and lower frequencies are associated with tighter, more highly populated solvent cages. The CO₂ interacts most strongly with the anions, and these interactions persist for more than 1 ns. The second strongest interactions are with the imidazolium cation ring that last 100 ps, and the weakest interactions are with the cation butyl tail that persist for 10 ps. The principal contributors to spectral diffusion of the CO₂ asymmetric stretch vibrational frequency due to the dynamical evolution of the solvent are through Lennard-Jones interactions at short times and electrostatics at long times.

Dynamically Disordered Lattice in a Layered Pb-I-SCN Perovskite Thin Film Probed by Two-Dimensional Infrared Spectroscopy

The dynamically flexible lattices in lead halide perovskites may play important roles in extending carrier recombination lifetime in 3D perovskite solar-cell absorbers and in exciton self-trapping in 2D perovskite white-light phosphors. Two-dimensional infrared (2D IR) spectroscopy was applied to study a recently reported Pb-I-SCN layered perovskite. The Pb-I-SCN perovskite was spin-coated on a SiO₂ surface as a thin film, with a thickness of ∼100 nm, where the S¹²CN^- anions were isotopically diluted with the ratio of S¹²CN:S¹³CN = 5:95 to avoid vibrational coupling and excitation transfer between adjacent SCN^- anions. The ¹²CN stretch mode of the minor S¹²CN^- component was the principal vibrational probe that reported on the structural evolution through 2D IR spectroscopy. Spectral diffusion was observed with a time constant of 4.1 ± 0.3 ps. Spectral diffusion arises from small structural changes that result in sampling of frequencies within the distribution of frequencies comprising the inhomogeneously broadened infrared absorption band. These transitions among discrete local structures are distinct from oscillatory phonon motions of the lattice. To accurately evaluate the structural dynamics through measurement of spectral diffusion, the vibrational coupling between adjacent SCN^- anions had to be carefully treated. Although the inorganic layers of typical 2D perovskites are structurally isolated from each other, the 2D IR data demonstrated that the layers of the Pb-I-SCN perovskite are vibrationally coupled. When both S¹²CN^- and S¹³CN^- were pumped simultaneously, cross-peaks between S¹²CN and S¹³CN vibrations and an oscillating 2D band shape of the S¹²CN^- vibration were observed. Both observables demonstrate vibrational coupling between the closest SCN^- anions, which reside in different inorganic layers. The thin films and the isotopic dilution produced exceedingly small vibrational echo signal fields; measurements were made possible using the near-Brewster's angle reflection pump-probe geometry.

Enthalpic Driving Force for the Selective Absorption of CO2 by an Ionic Liquid

Molecular dynamics (MD) simulations validated against two-dimensional infrared (2D-IR) measurements of CO₂ in an imidazolium-based ionic liquid have revealed new insights into the mechanism of CO₂ solvation. The first solvation shell around CO₂ has a distinctly quadrupolar structure, with strong negative charge density around the CO₂ carbon atom and positive charge density near the CO₂ oxygen atoms. When CO₂ is modeled without atomic charges (thus removing its strong quadrupole moment), its solvation shell weakens and changes significantly into a structure that is similar to that of N₂ in the same liquid. The solvation shell of CO₂ evolves more quickly when its quadrupole is removed, and we find evidence that solvent cage dynamics is measured by 2D-IR spectroscopy. We also find that the solvent cage evolution of N₂ is similar to that of CO₂ with no atomic charges, implying that the weaker quadrupole of N₂ is responsible for its higher diffusion and lower absorption in ionic liquids.

Micellar and bicontinuous microemulsion structures show different solute–solvent interactions: a case study using ultrafast nonlinear infrared spectroscopy

Using aqueous and organic probes to simultaneously explore the structural dynamics of reverse micellar and bicontinuous microemulsion structures.

Ultrafast dynamics of ionic liquids in colloidal dispersion

Ionic liquid (IL)-surfactant complexes have significance both in applications and fundamental research, but their underlying dynamics are not well understood. We apply polarization-controlled two-dimensional infrared spectroscopy (2D-IR) to study the dynamics of [BMIM][SCN]/surfactant/solvent model systems. We examine the effect of the choice of surfactants and solvent, and the IL-to-surfactant ratio (W-value), with a detailed analysis of the orientation and structural dynamics of each system. Different surfactants create very different environments for the entrapped ILs, ranging from a semi-static micro-environment to a fluxional environment that evolves even faster than the bulk IL. The oil-phase also clearly affects the microscopic dynamics. The anisotropy decay for entrapped ILs completes within 10 ps, which is similar to free thiocyanate ion in water, while a significant reorientation-induced spectral diffusion (RISD) effect is observed. The entrapped ionic liquid are highly dynamic for all W-values, and no core-shell structure is observed. We hypothesize that, instead of an ionic liquid-reverse micelle (IL-RM), the microscopic structure of this system is small colloidal dispersions or pairs of IL and surfactants. A detailed analysis of the polarization-controlled 2D-IR spectra of AOT system reveals a potential ion-exchange mechanism.

Identification of Ion-Pair Structures in Solution by Vibrational Stark Effects

Ion pairing is a fundamental consideration in many areas of chemistry and has implications in a wide range of sciences and technologies that include batteries and organic photovoltaics. Ions in solution are known to inhabit multiple possible states, including free ions (FI), contact ion pairs (CIP), and solvent-separated ion pairs (SSIP). However, in solutions of organic radicals and nonmetal electrolytes, it is often difficult to distinguish between these states. In the first part of this work, we report evidence for the formation of SSIPs in low-polarity solvents and distinct measurements of CIP, SSIP, and FI, by using the ν(C≡N) infrared (IR) band of a nitrile-substituted fluorene radical anion. Use of time-resolved IR detection following pulse radiolysis allowed us to unambiguously assign the peak of the FI. In the presence of nonmetal electrolytes, two distinct red-shifted peaks were observed and assigned to the CIP and SSIP. The assignments are interpreted in the framework of the vibrational Stark effect (VSE) and are supported by (1) the solvent dependence of ion-pair populations, (2) the observation of a cryptand-separated sodium ion pair that mimics the formation of SSIPs, and (3) electronic structure calculations. In the second part of this work, we show that a blue-shift of the ν(C≡N) IR band due to the VSE can be induced in a nitrile-substituted fluorene radical anion by covalently tethering it to a metal-chelating ligand that forms an intramolecular ion pair upon reduction and complexation with sodium ion. This adds support to the conclusion that the shift in IR absorptions by ion pairing originates from the VSE. These results combined show that we can identify ion-pair structures by using the VSE, including the existence of SSIPs in a low-polarity solvent.

Molten salts and energy related materials

Molten salts have been known for centuries and have been used for the extraction of aluminium for over one hundred years and as high temperature fluxes in metal processing. This and other molten salt routes have gradually become more energy efficient and less polluting, but there have been few major breakthroughs. This paper will explore some recent innovations that could lead to substantial reductions in the energy consumed in metal production and in carbon dioxide production. Another way that molten salts can contribute to an energy efficient world is by creating better high temperature fuel cells and novel high temperature batteries, or by acting as the medium that can create novel materials that can find applications in high energy batteries and other energy saving devices, such as capacitors. Carbonate melts can be used to absorb carbon dioxide, which can be converted into C, CO and carbon nanoparticles. Molten salts can also be used to create black silicon that can absorb more sunlight over a wider range of wavelengths. Overall, there are many opportunities to explore for molten salts to play in an efficient, low carbon world.

Two-dimensional ultrafast vibrational spectroscopy of azides in ionic liquids reveals solute-specific solvation

The stereochemistry and the reaction rates of bimolecular nucleophilic substitution reactions involving azides in ionic liquids are governed by solute-solvent interactions. Two-dimensional ultrafast vibrational spectroscopy (2D-IR) shows that the picosecond dynamics of inorganic azides are substantially slower than organic azides in a series of homologous imidazolium ionic liquids. In water, both organic and inorganic azides spectrally diffuse with a ∼2 ps time constant. In the aprotic solvent tetrahydrofuran, both kinds of azides spectrally diffuse on a timescale >5 ps. In ionic liquids, like 1-butyl-3-methylimidazolium tetrafluoroborate ([BMIM][BF4]), organic azides spectrally diffuse with a 2-4 ps time constant, and inorganic azides spectrally diffuse with a >40 ps time constant. Such a striking difference suggests that neutral (organic) and charged (inorganic) azides are incorporated in the ionic liquids with different solvation structures.