Dynamics of Water Molecules and Ions in Concentrated Lithium Chloride Solutions Probed with Ultrafast 2D IR Spectroscopy

Modulating vibrational energy redistribution in highly conjugated systems

Elucidating the nature of intramolecular vibrational energy redistribution (IVR) can guide the design of molecular wires. The ability to steer these processes through a mechanistic understanding of IVR is assessed by utilizing two-dimensional infrared (2D IR) spectroscopy. 2D IR spectroscopy allows for the direct investigation of timescales of energy transfer within three aromatic molecular scaffolds: 4'-azido-[1,1'-biphenyl]-4-carbonitrile (PAB), 2'-azido-[1,1'-biphenyl]-4-carbonitrile (OAB), and 4'-(azidomethyl)-[1,1'-biphenyl]-4-carbonitrile (PAMB). Energy transfer pathways between azido (N3)- and cyano (CN)-vibrational reporters uncover the importance of Fermi resonances, anharmonic coupling, and specific structural components in directing energy flow. Among these systems, PAB exhibits the fastest energy transfer (22 ps), facilitated by its co-planar biphenyl structure, enabling strong π-π stacking interactions to optimize vibrational coupling. In contrast, OAB demonstrates a moderate IVR timescale (38 ps) due to an orthogonal molecular plane and steric hindrance, which disrupts coupling pathways. PAMB, with a para-methylene group, introduces a structural bottleneck that significantly impedes energy flow, slowing down the energy transfer to 84 ps. The observed IVR rates align with computational predictions, highlighting intermediate ring modes in PAB as efficient energy transfer bridges, a mechanism that is less pronounced in OAB and PAMB. This study demonstrates that IVR is dictated not only by anharmonic coupling strengths but also by the extended alignment of vibrational modes across molecular planes and their delocalization within aromatic scaffolds. By modulating structural features, such as steric constraints and π-π interactions, we provide a framework for tailoring energy flow in conjugated molecular systems. These findings offer new insights into IVR dynamics for applications in molecular electronics.

Correlating Viscosity Trends and Ultrafast Structural Dynamics in TMSO-Water and SFL-Water Binary Mixtures

The hydrogen bonding dynamics in tetramethylene sulfoxide (TMSO)-water and sulfolane (SFL)-water binary mixtures were investigated by using FTIR spectroscopy, ultrafast IR spectroscopy, and molecular dynamics (MD) simulations. Despite TMSO and SFL sharing a similar cyclic backbone, markedly different viscosity behaviors were observed. By employing SCN- as a local probe, its distinct reorientational dynamics was observed that directly correlates with the contrasting viscosity trends in the two systems. In TMSO-water solutions, strong solute-water hydrogen bonding dominates water-water interactions, leading to a viscosity maximum at intermediate concentrations. Conversely, the dynamics of water molecules in SFL-water solutions is decoupled from bulk viscosity trends due to the weaker solute-water interactions. MD simulations further elucidate how the interplay between solute-water and water-water hydrogen bonding governs the viscosity trends. This work advances our understanding of hydrogen bonding in complex aqueous environments and provides a systematic approach for connecting molecular-level interactions to macroscopic fluid properties.

Revealing Lithium Ion Transport Mechanisms and Solvation Structures in Carbonate Electrolytes

Optimizing lithium-ion battery (LIB) electrolytes is essential for high-current applications such as electric vehicles, yet experimental techniques to characterize the complex structural dynamics responsible for the lithium transport within these electrolytes are limited. In this study, we used ultrafast infrared spectroscopy to measure chemical exchange, spectral diffusion, and solvation structures across a wide range of lithium concentrations in propylene carbonate-based LiTFSI (lithium bis(trifluoromethanesulfonimide) electrolytes, with the CN stretch of phenyl selenocyanate as the long-lived vibrational probe. Phenyl selenocyanate is shown to be an excellent dynamical surrogate for propylene carbonate in Li+ solvation clusters. A strong correlation between exchange times and ionic conductivity was observed. This correlation and other observations suggest structural diffusion as the primary transport mechanism rather than vehicular diffusion. Additionally, spectral diffusion observables measured by the probe were directly linked to the desolvation dynamics of the Li+ clusters, as supported by density functional theory and molecular dynamics simulations. These findings provide detailed molecular-level insights into LIB electrolytes' transport dynamics and solvation structures, offering rational design pathways to advanced electrolytes for next-generation LIBs.

Two-Dimensional Infrared Spectroscopy Resolves the Vibrational Landscape in Donor-Bridge-Acceptor Complexes with Site-Specific Isotopic Labeling

Donor-bridge-acceptor complexes (D-B-A) are important model systems for understanding of light-induced processes. Here, we apply two-color two-dimensional infrared (2D-IR) spectroscopy to D-B-A complexes with a trans-Pt(II) acetylide bridge (D-C≡C-Pt-C≡C-A) to uncover the mechanism of vibrational energy redistribution (IVR). Site-selective 13C isotopic labeling of the bridge is used to decouple the acetylide modes positioned on either side of the Pt-center. Decoupling of the D-acetylide- from the A-acetylide- enables site-specific investigation of vibrational energy transfer (VET) rates, dynamic anharmonicities, and spectral diffusion. Surprisingly, the asymmetrically labeled D-B-A still undergoes intramolecular IVR between acetylide groups even though they are decoupled and positioned across a heavy atom usually perceived as a "vibrational bottleneck". Further, the rate of population transfer from the bridge to the acceptor was both site-specific and distance dependent. We show that vibrational excitation of the acetylide modes is transferred to ligand-centered modes on a subpicosecond time scale, followed by VET to solvent modes on the time scale of a few picoseconds. We also show that isotopic substitution does not affect the rate of spectral diffusion, indicating that changes in the vibrational dynamics are not a result of differences in local environment around the acetylides. Oscillations imprinted on the decay of the vibrationally excited acceptor-localized carbonyl modes show they enter a coherent superposition of states after excitation that dephases over 1-2 ps, and thus cannot be treated as independent in the 2D-IR spectra. These findings elucidate the vibrational landscape governing IR-mediated electron transfer and illustrate the power of isotopic labeling combined with multidimensional IR spectroscopy to disentangle vibrational energy propagation pathways in complex systems.

Effects of Nanoconfinement on Dynamics in Concentrated Aqueous Magnesium Chloride Solutions

Water behavior in various natural and manufactured settings is influenced by confinement in organic or inorganic frameworks and the presence of solutes. Here, the effects on dynamics from both confinement and the addition of solutes are examined. Specifically, water and ion dynamics in concentrated (2.5-4.2 m) aqueous magnesium chloride solutions confined in mesoporous silica (2.8 nm pore diameter) were investigated using polarization selective pump-probe and 2D infrared spectroscopies. Fitting the rotational and spectral diffusion dynamics measured by the vibrational probe, selenocyanate, with a previously developed two-state model revealed distinct behaviors at the interior of the silica pores (core state) and near the wall of the confining framework (shell state). The shell dynamics are noticeably slower than the bulk, or core, dynamics. The concentration-dependent slowing of the dynamics aligns with behavior in the bulk solutions, but the spectrally separated water-associated and Mg2+-associated forms of the selenocyanate probe exhibit different responses to confinement. The disparity in the complete reorientation times is larger upon confinement, but the spectral diffusion dynamics become more similar near the silica surface. The length scales that characterize the transition from surface-influenced to bulk-like behavior for the salt solutions in the pores are discussed and compared to those of pure water and an organic solvent confined in the same pores. These comparisons offer insights into how confinement modulates the properties of different liquids.

Anionic Effects on Concentrated Aqueous Lithium Ion Dynamics

The dynamics, orientational anisotropy, diffusivity, viscosity, and density were measured for concentrated lithium salt solutions, including lithium chloride (LiCl), lithium bromide (LiBr), lithium nitrite (LiNO2), and lithium nitrate (LiNO3), with methyl thiocyanate as an infrared vibrational probe molecule, using two-dimensional infrared spectroscopy (2D IR), nuclear magnetic resonance (NMR) spectroscopy, and viscometry. The 2D IR, NMR, and viscosity results show that LiNO2 exhibits longer correlation times, lower diffusivity, and nearly 4 times greater viscosity compared to those of the other lithium salt solutions of the same concentration, suggesting that nitrite anions may strongly facilitate structure formation via strengthening water-ion network interactions, directly impacting bulk solution properties at sufficiently high concentrations. Additionally, the LiNO2 and LiNO3 solutions show significantly weakened chemical interactions between the lithium cations and the methyl thiocyanate when compared with those of the lithium halide salts.

Water-in-Salt: Fast Dynamics, Structure, Thermodynamics, and Bulk Properties

We present concentration-dependent dynamics of highly concentrated LiBr solutions and LiCl temperature-dependent dynamics for two high concentrations and compare the results to those of prior LiCl concentration-dependent data. The dynamical data are obtained using ultrafast optical heterodyne-detected optical Kerr effect (OHD-OKE). The OHD-OKE decays are composed of two pairs of biexponentials, i.e., tetra-exponentials. The fastest decay (t1) is the same as pure water's at all concentrations within error, while the second component (t2) slows slightly with concentration. The slower components (t3 and t4), not present in pure water, slow substantially, and their contributions to the decays increase significantly with increasing concentration, similar to LiCl solutions. Simulations of LiCl solutions from the literature show that the slow components arise from large ion/water clusters, while the fast components are from ion/water structures that are not part of large clusters. Temperature-dependent studies (15-95 °C) of two high LiCl concentrations show that decreasing the temperature is equivalent to increasing the room temperature concentration. The LiBr and LiCl concentration dependences and the two LiCl concentrations' temperature dependences all have bulk viscosities that are linearly dependent on τcslow, the correlation time of the slow dynamics (weighted averages of t3 and t4). Remarkably, all four viscosity vs 1/τCslow plots fall on the same line. Application of transition state theory to the temperature-dependent data yields the activation enthalpies and entropies for the dynamics of the large ion/water clusters, which underpin the bulk viscosity.

Sponge-phase Lipid Droplets as Synthetic Organelles: An Ultrafast Study of Hydrogen Bonding and Interfacial Environments

Bottom-up design of biomimetic organelles has gained recent attention as a route towards understanding the transition between non-living matter and life. Despite various artificial lipid membranes being developed, the specific relations between lipid structure, composition, interfacial properties, and morphology are not currently understood. Sponge-phase droplets contain dense, nonlamellar lipid bilayer networks that capture the complexities of the endoplasmic reticulum (ER), making them ideal artificial models of such organelles. Here, we combine ultrafast two-dimensional infrared (2D IR) spectroscopy and molecular dynamics simulations to investigate the interfacial H-bond networks in sponge-phase droplets composed of glycolipid and nonionic detergents. In the sponge phase, the interfacial environments are more hydrated and water molecules confined to the nanometer-scale aqueous channels in the sponge phase exhibit dynamics that are significantly slower compared to bulk water. Surfactant configurations and microscopic phase separation play a dominant role in determining membrane curvature and slow dynamics observed in the sponge phase. The studies suggest that H-bond networks within the nanometer-scale channels are disrupted not only by confinement but also by the interactions of surfactants, which extend 1-2 nm from the bilayer surface. The results provide a molecular-level description for controlling phase and morphology in the design of synthetic lipid organelles.

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.

Progress and Prospects in Optical Ultrafast Microscopy in the Visible Spectral Region: Transient Absorption and Two-Dimensional Microscopy

Ultrafast optical microscopy, generally employed by incorporating ultrafast laser pulses into microscopes, can provide spatially resolved mechanistic insight into scientific problems ranging from hot carrier dynamics to biological imaging. This Review discusses the progress in different ultrafast microscopy techniques, with a focus on transient absorption and two-dimensional microscopy. We review the underlying principles of these techniques and discuss their respective advantages and applicability to different scientific questions. We also examine in detail how instrument parameters such as sensitivity, laser power, and temporal and spatial resolution must be addressed. Finally, we comment on future developments and emerging opportunities in the field of ultrafast microscopy.

Inhibition of vibrational energy flow within an aromatic scaffold via heavy atom effect

The regulation of intramolecular vibrational energy redistribution (IVR) to influence energy flow within molecular scaffolds provides a way to steer fundamental processes of chemistry, such as chemical reactivity in proteins and design of molecular diodes. Using two-dimensional infrared (2D IR) spectroscopy, changes in the intensity of vibrational cross-peaks are often used to evaluate different energy transfer pathways present in small molecules. Previous 2D IR studies of para-azidobenzonitrile (PAB) demonstrated that several possible energy pathways from the N3 to the cyano-vibrational reporters were modulated by Fermi resonance, followed by energy relaxation into the solvent [Schmitz et al., J. Phys. Chem. A 123, 10571 (2019)]. In this work, the mechanisms of IVR were hindered via the introduction of a heavy atom, selenium, into the molecular scaffold. This effectively eliminated the energy transfer pathway and resulted in the dissipation of the energy into the bath and direct dipole-dipole coupling between the two vibrational reporters. Several structural variations of the aforementioned molecular scaffold were employed to assess how each interrupted the energy transfer pathways, and the evolution of 2D IR cross-peaks was measured to assess the changes in the energy flow. By eliminating the energy transfer pathways through isolation of specific vibrational transitions, through-space vibrational coupling between an azido (N3) and a selenocyanato (SeCN) probe is facilitated and observed for the first time. Thus, the rectification of this molecular circuitry is accomplished through the inhibition of energy flow using heavy atoms to suppress the anharmonic coupling and, instead, favor a vibrational coupling pathway.

Dynamics of Concentrated Aqueous Lithium Chloride Solutions Investigated with Optical Kerr Effect Experiments

We report the dynamics of concentrated lithium chloride aqueous solutions over a range of moderate to high concentrations. Concentrations (1-29 to 1-3.3 LiCl-water) were studied in which, at the highest concentrations, there are far too few water molecules to solvate the ions. The measurements were made with optically heterodyne-detected optical Kerr effect experiments, a non-resonant technique able to observe dynamics over a wide range of time scales and signal amplitudes. While the pure water decay is a biexponential, the LiCl-water decays are tetra-exponentials at all concentrations. The faster two decays arise from water dynamics, while the slower two decays reflect the dynamics of the ion-water network. The fastest decay (t₁) is the same as pure water at all concentrations. The second decay (t₂) is also the same as that of pure water at the lower concentrations, and then, it slows with increasing concentration. The slower dynamics (t₃ and t₄), which do not have counterparts in pure water, arise from ion-water complexes and, at the highest concentrations, an extended ion-water network. Comparisons are made between the concentration dependence of the observed dynamics and simulations of structural changes from the literature, which enable the assignment of dynamics to specific ion-water structures. The concentration dependences of the bulk viscosity and the ion-water network dynamics are directly correlated. The correlation provides an atomistic-level understanding of the viscosity.

Dynamics of Acrylamide Hydrogels, Polymers, and Monomers in Water Measured with Optical Heterodyne-Detected Optical Kerr Effect Spectroscopy

The ultrafast dynamics of acrylamide monomers (AAm), polyacrylamide (PAAm), and polyacrylamide hydrogels (PAAm-HG) in water were studied using optical heterodyne-detected optical Kerr effect (OHD-OKE) spectroscopy. Previous ultrafast infrared (IR) measurements of the water dynamics showed that at the same concentration of the acrylamide moiety, AAm, PAAm, and PAAm-HG exhibited identical water dynamics and that these dynamics slowed with increasing concentration. In contrast to the IR measurements, OHD-OKE experiments measure the dynamics of both the water and the acrylamide species, which occur on different time scales. In this study, the dynamics of all the acrylamide systems slowed with increasing concentration. We found that AAm exhibits tetraexponential decays, the longest component of which followed Debye-Stokes-Einstein behavior except for the highest concentration, 40% (w/v). Low concentrations of PAAm followed a single power law decay, while high concentrations of PAAm and all concentrations of PAAm-HG decayed with two power laws. The highest concentrations, 25% and 40%, of PAAm and PAAm-HG showed nearly identical dynamics. We interpreted this result as reflecting a similar extent of chain-chain interactions. At low concentrations, PAAm displays non-Markovian, single-chain dynamics (single power law), but PAAm displays entangled chain-chain interactions at high concentrations (two power laws). PAAm-HG has chain-chain interactions at all concentrations that arise from the cross-linking. At high concentrations, the dynamics of the entangled of PAAm become identical within error as those of the cross-linked PAAm-HG.

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.

2D spectroscopies from condensed phase dynamics: Accessing third-order response properties from equilibrium multi-time correlation functions

The third-order response lies at the heart of simulating and interpreting nonlinear spectroscopies ranging from two dimensional infrared (2D-IR) to 2D electronic (2D-ES), and 2D sum frequency generation (2D-SFG). The extra time and frequency dimensions in these spectroscopies provides access to rich information on the electronic and vibrational states present, the coupling between them, and the resulting rates at which they exchange energy that are obscured in linear spectroscopy, particularly for condensed phase systems that usually contain many overlapping features. While the exact quantum expression for the third-order response is well established it is incompatible with the methods that are practical for calculating the atomistic dynamics of large condensed phase systems. These methods, which include both classical mechanics and quantum dynamics methods that retain quantum statistical properties while obeying the symmetries of classical dynamics, such as LSC-IVR, Centroid Molecular Dynamics (CMD) and Ring Polymer Molecular Dynamics (RPMD) naturally provide short-time approximations to the multi-time symmetrized Kubo transformed correlation function. Here, we show how the third-order response can be formulated in terms of equilibrium symmetrized Kubo transformed correlation functions. We demonstrate the utility and accuracy of our approach by showing how it can be used to obtain the third-order response of a series of model systems using both classical dynamics and RPMD. In particular, we show that this approach captures features such as anharmonically induced vertical splittings and peak shifts while providing a physically transparent framework for understanding multidimensional spectroscopies.

Exchange-Mediated Transport in Battery Electrolytes: Ultrafast or Ultraslow?

Understanding the mechanisms of charge transport in batteries is important for the rational design of new electrolyte formulations. Persistent questions about ion transport mechanisms in battery electrolytes are often framed in terms of vehicular diffusion by persistent ion-solvent complexes versus structural diffusion through the breaking and reformation of ion-solvent contacts, i.e., solvent exchange events. Ultrafast two-dimensional (2D) IR spectroscopy can probe exchange processes directly via the evolution of the cross-peaks on picosecond time scales. However, vibrational energy transfer in the absence of solvent exchange gives rise to the same spectral signatures, hiding the desired processes. We employ 2D IR on solvent resonances of a mixture of acetonitrile isotopologues to differentiate chemical exchange and energy-transfer dynamics in a comprehensive series of Li⁺, Mg²⁺, Zn²⁺, Ca²⁺, and Ba²⁺ bis(trifluoromethylsulfonyl)imide electrolytes from the dilute to the superconcentrated regime. No exchange phenomena occur within at least 100 ps, regardless of the ion identity, salt concentration, and presence of water. All of the observed spectral dynamics originate from the intermolecular energy transfer. These results place the lower experimental boundary on the ion-solvent residence times to several hundred picoseconds, much slower than previously suggested. With the help of MD simulations and conductivity measurements on the Li⁺ and Zn²⁺ systems, we discuss these results as a continuum of vehicular and structural modalities that vary with concentration and emphasize the importance of collective electrolyte motions to ion transport. These results hold broadly applicable to many battery-relevant ions and solvents.

Water Dynamics and Structure of Highly Concentrated LiCl Solutions Investigated Using Ultrafast Infrared Spectroscopy

In highly concentrated salt solutions, the water hydrogen bond (H-bond) network is completely disrupted by the presence of ions. Water is forced to restructure as dictated by the water-ion and ion-ion interactions. Using ultrafast polarization-selective pump-probe (PSPP) spectroscopy measurements of the OD stretch of dilute HOD, we demonstrate that the limited water-water H-bonding present in concentrated lithium chloride solutions (up to four waters per ion pair) is, on average, stronger than that occurring in bulk water. Furthermore, information on the orientational dynamics and the angular restriction of water H-bonded to both water oxygens and chloride anions was obtained through analysis of the frequency-dependent anisotropy decays. It was found that, when the salt concentration increased, the water showed increasing restriction and slowing at frequencies correlated with strong H-bonding. The angular restriction of the water molecules and strengthening of water-water H-bonds are due to the formation of a water-ion network not present in bulk water and dilute salt solutions. The structural evolution of the ionic medium was also observed through spectral diffusion of the OD stretch using 2D IR spectroscopy. Compared to bulk water, there is significant slowing of the biexponential spectral diffusion dynamics. The slowest component of the spectral diffusion (13 ps) is virtually identical to the time for complete reorientation of HOD measured with the PSPP experiments. This result suggests that the slowest component of the spectral diffusion reflects rearrangement of water molecules in the water-ion network.

A polarization scheme that resolves cross-peaks with transient absorption and eliminates diagonal peaks in 2D spectroscopy

Two-dimensional (2D) optical spectroscopy contains cross-peaks that are helpful features for determining molecular structure and monitoring energy transfer, but they can be difficult to resolve from the much more intense diagonal peaks. Transient absorption (TA) spectra contain transitions similar to cross-peaks in 2D spectroscopy, but in most cases they are obscured by the bleach and stimulated emission peaks. We report a polarization scheme, <0°,0°,+θ₂(t₂),-θ₂(t₂)>, that can be easily implemented in the pump-probe beam geometry, used most frequently in 2D and TA spectroscopy. This scheme removes the diagonal peaks in 2D spectroscopies and the intense bleach/stimulated emission peaks in TA spectroscopies, thereby resolving the cross-peak features. At zero pump-probe delay, θ₂ = 60° destructively interferes two Feynman paths, eliminating all signals generated by field interactions with four parallel transition dipoles, and the intense diagonal and bleach/stimulated emission peaks. At later delay times, θ₂(t₂) is adjusted to compensate for anisotropy caused by rotational diffusion. When implemented with TA spectroscopy or microscopy, the pump-probe spectrum is dominated by the cross-peak features. The local oscillator is also attenuated, which enhances the signal two times. This overlooked polarization scheme reduces spectral congestion by eliminating diagonal peaks in 2D spectra and enables TA spectroscopy to measure similar information given by cross-peaks in 2D spectroscopy.

Characterization of Acetonitrile Isotopologues as Vibrational Probes of Electrolytes

Acetonitrile has emerged as a solvent candidate for novel electrolyte formulations in metal-ion batteries and supercapacitors. It features a bright local C≡N stretch vibrational mode whose infrared (IR) signature is sensitive to battery-relevant cations (Li⁺, Mg²⁺, Zn²⁺, Ca²⁺) both in pure form and in the presence of water admixture across a full possible range of concentrations from the dilute to the superconcentrated regime. Stationary and time-resolved IR spectroscopy thus emerges as a natural tool to study site-specific intermolecular interactions from the solvent perspective without introducing an extrinsic probe that perturbs solution morphology and may not represent the intrinsic dynamics in these electrolytes. The metal-coordinated acetonitrile, water-separated metal-acetonitrile pair, and free solvent each have a distinct vibrational signature that allows their unambiguous differentiation. The IR band frequency of the metal-coordinated acetonitrile depends on the ion charge density. To study the ion transport dynamics, it is necessary to differentiate energy-transfer processes from structural interconversions in these electrolytes. Isotope labeling the solvent is a necessary prerequisite to separate these processes. We discuss the design principles and choice of the CD₃¹³CN label and characterize its vibrational spectroscopy in these electrolytes. The Fermi resonance between ¹³C≡N and C-D stretches complicates the spectral response but does not prevent its effective utilization. Time-resolved two-dimensional (2D) IR spectroscopy can be performed on a mixture of acetonitrile isotopologues and much can be learned about the structural dynamics of various species in these formulations.

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.

Binary structure and dynamics of the hydrogen bonds in the hydration shells of ions

Ion-specific effects of cations (Li+, Na+, K+, Mg2+, Ca2+) and anions (F-, Cl-) on the hydrogen bond structure and dynamics of the coordination waters in the hydration shells have been studied using molecular dynamics simulations. Our simulations indicate that the hydrogen bonds between the first and second hydration shell waters show binary structural and dynamic properties. The hydrogen bond with a first shell water as the donor (HD) is strengthened, while those with a first shell water as the acceptor (HA) are weakened. For a hydrated anion, this binary effect reverses, but is less significant. This ion-specific binary effect correlates with the size and the valence of the ion, and is more significant for the strong kosmotropic ions of high charge density.

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.

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.

Cooperative CO2 absorption by amino acid-based ionic liquids with balanced dual sites

In this study, a variety of functionalized ILs with dual sites including amino acid group (AA) and basic anion (R) were synthesized to investigate the suppression and cooperation between the sites in CO₂ absorption. The basic anions selected in this study with different basicity include sulfonate (Su), carboxylate (Ac), imidazolium (Im), and indolium (Ind). These ILs ([P₆₆₆₁₄]₂[AA-R]) were applied to CO₂ absorption. The results present that CO₂ capacity increases first and then decreases later with the continuous increase in the activity of the anion site. Combined with CO₂ absorption experiments, IR and NMR spectroscopic analyses and DFT calculation demonstrate that the ability of one site to capture CO₂ would be suppressed when the activity of another site is much stronger. Thus, the cooperation of dual site-functionalized ILs and high CO₂ capacity might be achieved through balancing the two sites to be equivalent. Based on this point, [P₆₆₆₁₄]₂[5Am-iPA] was further synthesized by taking the advantage of the conjugated benzene ring. As expected, [P₆₆₆₁₄]₂[5Am-iPA] showed capacity as high as 2.38 mol CO₂ per mol IL at 30 °C and 1 bar without capacity decrease even after 10 times recycling performance of CO₂ absorption and desorption.