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

Distinctive signatures and ultrafast dynamics of Brønsted sites, silanol nests and adsorbed water in zeolites revealed by 2D-IR spectroscopy

Characterising hydroxyl groups in zeolites and other amorphous solids often relies on methods such as IR and NMR spectroscopy. Their power to distinguish different types of hydroxyl groups diminishes when band broadening from hydrogen bonding and structural heterogeneity occurs. In support of this problem, we report in situ femtosecond 2D-IR spectroscopy of some of the different types of hydroxyl groups present in zeolites. Despite the samples studied being optically scattering pellets, we show that their structural and rotational dynamics can be determined. We show that the hydroxyl groups of Brønsted acid sites, silanol defects and water of hydration display distinct features in their 2D-IR spectra. Brønsted site hydroxyl group structural distributions have characteristic inhomogeneously broadened 2D-IR bandshapes. Water of hydration and partially hydrogen bonded silanol groups give unique 2D-IR cross peak signatures off-diagonal. Hydrogen bonded silanol groups arising from vacancy defects (silanol nests) show a distinctive 2D-IR signature with unique ultrafast dynamics observed to be identical between ZSM-5 and silicalite-1. 2D-IR spectroscopy makes IR measurements quantitative, and we use this property to estimate the concentration of ZSM-5 silanol nest hydroxyl groups relative to the number of Brønsted sites. Overlapping silanol nest spectral features are revealed by frequency dependence of their vibrational lifetime. In contrast to other framework hydroxyls, the silanol nest band shows picosecond 2D-IR anisotropy decay and spectral diffusion. The signatures of nest structural mobility revealed here presents new opportunities to understand these hitherto elusive structural defects.

Evolution of Water State in the "Ionic Liquid EAN/Al(NO3)3/H2O" System Depending on Temperature and Aluminum Salt Concentration

Although lithium salts are the most successful materials for electrochemical devices so far, there is growing interest in electrolytes containing divalent and trivalent salts, which are characterized by a much higher energy density than lithium systems. For this purpose, the physical and chemical properties of solutions of aluminum nitrate in ethylammonium nitrate (EAN) were studied in a number of works. However, it was demonstrated by the 1H NMR method that solutions contain water (from aluminum crystallohydrates) in the first solvated shell of the Al3+ cation even after usual drying procedures, and this solvated water cannot be directly detected by the standard Karl Fischer analysis. The goal of this work is to check the possibility/impossibility of removing water from the "Ionic Liquid EAN/Al(NO3)3/H2O" system and to investigate the evolution of water state in the system depending on temperature and aluminum salt concentration. To solve the problem, the NMR spectroscopy of 1H, 27Al, and 14N nuclei was used. To interpret the NMR spectra, quantum-chemical calculations of electron shielding constants for the 1H and 27Al nuclei in mixed cation-anion complexes were involved.

Fast Collective Hydrogen-Bond Dynamics in Hexafluoroisopropanol Related to its Chemical Activity

Using fluorinated mono-alcohols, in particular hexafluoro-isopropanol (HFIP), as a solvent can enhance chemical reaction rates in a spectacular manner. Previous work has shown evidence that this enhancement is related to the hydrogen-bond structure of these liquids. Here, we investigate the hydrogen-bond dynamics of HFIP and compare it to that of its non-fluorinated analog, isopropanol. Ultrafast infrared spectroscopy experiments show that the dynamics of individual hydrogen-bonds is about twice as slow in HFIP as in isopropanol. Surprisingly, from dielectric spectroscopy we find the opposite behavior for the dynamics of hydrogen-bonded clusters: collective rearrangements are 3 times faster in HFIP than in isopropanol. This difference indicates that the hydrogen-bonded clusters in HFIP are smaller than in isopropanol. The differences in cluster size can be traced to changes in the hydrogen-bond donor and acceptor strengths upon fluorination. The smaller cluster size can boost reaction rates in HFIP by increasing the concentration of reactive, terminal OH-groups of the clusters, whereas the fast collective dynamics can increase the rate of formation of hydrogen-bonds with the reactants. The longer lifetime of the individual hydrogen-bonds in HFIP can enhance the stability of the hydrogen-bonded clusters, and so increase the probability of reactant-solvent hydrogen-bonding.

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.

Cutting through the Noise: Extracting Dynamics from Ultrafast Spectra Using Dynamic Mode Decomposition

Coherent multidimensional spectroscopy provides experimental access to molecular structure and subpicosecond dynamics in solution. Dynamics are typically inferred from the evolution of lineshapes over a function of waiting time. Numerous spectral analysis methods, such as center/nodal line slope, have been developed to extract these dynamics. However, the extracted dynamics can depend heavily on subjective choices, such as the region selected for CLS analysis or the chosen models. In this study, we introduce a novel approach to extracting dynamics from ultrafast two-dimensional infrared (2D IR) spectra by using dynamic mode decomposition (DMD). As a data-driven method, DMD directly extracts spatiotemporal structures from the complex 2D IR spectra. We evaluated the performance of DMD in simulated and experimental spectra containing overlapped peaks. We show that DMD can retrieve the dynamics of overlapped transitions and cross peaks that are typically challenging to extract with traditional methods. In addition, we demonstrate that combining conditional generative adversarial neural networks with DMD can recover dynamics even at low signal-to-noise ratios. DMD methods do not require preliminary assumptions and can be readily extended to other multidimensional spectroscopies.

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.

Spectrophotometric Concentration Analysis Without Molar Absorption Coefficients by Two-Dimensional-Infrared and Fourier Transform Infrared Spectroscopy

A spectrophotometric method for determining relative concentrations of infrared (IR)-active analytes with unknown concentration and unknown molar absorption coefficient is explored. This type of method may be useful for the characterization of complex/heterogeneous liquids or solids, the study of transient species, and for other scenarios where it might be difficult to gain concentration information by other means. Concentration ratios of two species are obtained from their IR absorption and two-dimensional (2D)-IR diagonal bleach signals using simple ratiometric calculations. A simple calculation framework for deriving concentration ratios from spectral data is developed, extended to IR-pump-probe signals, and applied to the calculation of transition dipole ratios. Corrections to account for the attenuation of the 2D-IR signal caused by population relaxation, spectral overlap, wavelength-dependent pump absorption, inhomogeneous broadening, and laser intensity variations are described. A simple formula for calculating the attenuation of the 2D-IR signal due to sample absorption is deduced and by comparison with 2D-IR signals at varying total sample absorbance found to be quantitatively accurate. 2D-IR and Fourier transform infrared spectroscopy of two carbonyl containing species acetone and N-methyl-acetamide dissolved in D2O are used to experimentally confirm the validity of the ratiometric calculations. Finally, to address ambiguities over units and scaling of 2D-IR signals, a physical unit of 2D-IR spectral amplitude in mOD/c m - 1 is proposed.

Does Viscosity Drive the Dynamics in an Alcohol-Based Deep Eutectic Solvent?

Deep eutectic solvents, consisting of heterogeneous nanodomains of hydrogen-bonded networks, have evolved into a range of viscous fluids that find applications in several fields. As viscosity is known to influence the dynamics of other neoteric solvents like ionic liquids, understanding the effect of viscosity on deep eutectic solvents is crucial to realize their full potential. Herein, we combine polarization-selective pump-probe spectroscopy, two-dimensional infrared spectroscopy, and molecular dynamics simulations to elucidate the impact of viscosity on the dynamics of an alcohol-based deep eutectic solvent, ethaline. We compare the solvent fluctuation and solute reorientation time scales in ethaline with those in ethylene glycol, an ethaline constituent. Interestingly, we find that the solute's reorientation apparently scales the bulk viscosity of the solvent, but the same is not valid for the overall solvation dynamics. Using the variations in the estimated intercomponent hydrogen bond lifetimes, we show that a dissolved solute does not sense the bulk viscosity of the deep eutectic solvent; instead, it senses domain-specific local viscosity determined by the making and breaking of the hydrogen bond network. Our results indicate that understanding the domain-specific local environment experienced by the dissolved solute is of utmost importance in deep eutectic solvents.

Multiple Ensembles of the Hydrogen-bonded Network in Ethylammonium Nitrate versus Water from Vibrational Spectral Dynamics of SCN- Probe

We performed classical molecular dynamics simulations to monitor the structural interactions and ultrafast dynamical and spectral response in the protic ionic liquid, ethylammonium nitrate (EAN) and water using the nitrile stretching mode of thiocyanate ion (SCN-) as the vibrational probe. The normalized stretch frequency distribution of nitrile stretch of SCN- attains an asymmetric shape in EAN, indicating the existence of more than one hydrogen-bonding environment in EAN. We computed the 2D IR spectrum from classical trajectories, applying the response function formalism. Spectral diffusion dynamics in EAN undergo an initial rattling of the SCN - inside the local ion-cage occurring at a timescale of 0.10 ps, followed by the breakup of the ion-cage activating molecular diffusion at 7.86 ps timescale. In contrast, the dynamics of structural reorganization occur at a timescale of 0.58 ps in H 2 O. Hence, the time dependence of the frequency-frequency correlation function decay hints at the local molecular structure and ultrafast ion dynamics of the SCN - probe. The loss of frequency correlation read from the peak shape changes in the 2D correlation spectrum as a function of waiting time is faster in H 2 O than in EAN due to the enhanced structural ordering and higher viscosity of the latter. We provide an atomic-level interpretation of the solvation environment around SCN - in EAN and water, which indicates the multiple ensembles of the hydrogen bond network in EAN.

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.

Distinguishing steric and electrostatic molecular probe orientational ordering via their effects on reorientation-induced spectral diffusion

The theoretical framework for reorientation-induced spectral diffusion (RISD) describes the polarization dependence of spectral diffusion dynamics as measured with two-dimensional (2D) correlation spectroscopy and related techniques. Generally, RISD relates to the orientational dynamics of the molecular chromophore relative to local electric fields of the medium. The predictions of RISD have been shown to be very sensitive to both restricted orientational dynamics (generally arising from steric hindrance) and the distribution of local electric fields relative to the probe (electrostatic ordering). Here, a theory that combines the two effects is developed analytically and supported with numerical calculations. The combined effects can smoothly vary the polarization dependence of spectral diffusion from the purely steric case (least polarization dependence) to the purely electrostatic case (greatest polarization dependence). Analytic approximations of the modified RISD equations were also developed using the orientational dynamics of the molecular probe and two order parameters describing the degree of electrostatic ordering. It was found that frequency-dependent orientational dynamics are a possible consequence of the combined electrostatic and steric effects, providing a test for the applicability of this model to experimental systems. The modified RISD equations were then used to successfully describe the anomalous polarization-dependent spectral diffusion seen in 2D infrared spectroscopy in a polystyrene oligomer system that exhibits frequency-dependent orientational dynamics. The degree of polarization-dependent spectral diffusion enables the extent of electrostatic ordering in a chemical system to be quantified and distinguished from steric ordering.