Intermolecular vibrational energy exchange directly probed with ultrafast two dimensional infrared spectroscopy

Blue-Shifted and Broadened Fluorescence Enhancement by Visible and Mode-Selective Infrared Double Excitations

Mode-selective vibrational excitations to modify the electronic states of fluorescein dianion in methanol solutions are carried out with a femtosecond visible pulse synchronized with a tunable high-power, narrow-band picosecond infrared (IR) pulse. In this work, simultaneous intensity enhancement, peak blueshift, and line width broadening of fluorescence are observed in the visible/IR double resonance experiments. Comprehensive investigations on the modulation mechanism with scanning the vibrational excitation frequencies, tuning the time delay between the two excitation pulses, theoretical calculations, and nonlinear and linear spectroscopic measurements suggest that the fluorescence intensity enhancement is caused by the increase of the Franck-Condon factor induced by the vibrational excitations at the electronic ground state. Various enhancement effects are observed as vibrations initially excited by the IR photons relax to populate the vibrational modes of lower frequencies. The peak blueshift and line width broadening are caused by both increasing the Franck-Condon factors among different subensembles because of IR pre-excitation and the long-lived processes induced by the initial IR excitation. The results demonstrate that the fluorescence from the visible/IR double resonance experiments is not a simple sum frequency effect, and vibrational relaxations can produce profound effects modifying luminescence.

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.

Vibrational relaxation by methylated xanthines in solution: Insights from 2D IR spectroscopy and calculations

Two-dimensional infrared (2D IR) spectroscopy, infrared pump-infrared probe spectroscopy, and density functional theory calculations were used to study vibrational relaxation by ring and carbonyl stretching modes in a series of methylated xanthine derivatives in acetonitrile and deuterium oxide (heavy water). Isotropic signals from the excited symmetric and asymmetric carbonyl stretch modes decay biexponentially in both solvents. Coherent energy transfer between the symmetric and asymmetric carbonyl stretching modes gives rise to a quantum beat in the time-dependent anisotropy signals. The damping time of the coherent oscillation agrees with the fast decay component of the carbonyl bleach recovery signals, indicating that this time constant reflects intramolecular vibrational redistribution (IVR) to other solute modes. Despite their similar frequencies, the excited ring modes decay monoexponentially with a time constant that matches the slow decay component of the carbonyl modes. The slow decay times, which are faster in heavy water than in acetonitrile, approximately match the ones observed in previous UV pump-IR probe measurements on the same compounds. The slow component is assigned to intermolecular energy transfer to solvent bath modes from low-frequency solute modes, which are populated by IVR and are anharmonically coupled to the carbonyl and ring stretch modes. 2D IR measurements indicate that the carbonyl stretching modes are weakly coupled to the delocalized ring modes, resulting in slow exchange that cannot explain the common solvent-dependence. IVR is suggested to occur at different rates for the carbonyl vs ring modes due to differences in mode-specific couplings and not to differences in the density of accessible states.

Rethinking Vibrational Stark Spectroscopy: Peak Shifts, Line Widths, and the Role of Non-Stark Solvent Coupling

A vibration's transition frequency is partly determined by the first-order Stark effect, which accounts for the electric field experienced by the mode. Using ultrafast infrared pump-probe and FT-IR spectroscopies, we characterized both the 0 → 1 and 1 → 2 vibrational transitions' field-dependent peak positions and line widths of the CN stretching mode of benzonitrile (BZN) and phenyl selenocyanate (PhSeCN) in ten solvents. We present a theoretical model that decomposes the observed line width into a field-dependent Stark contribution and a field-independent non-Stark solvent coupling contribution (NSC). The model demonstrates that the field-dependent peak position is independent of the line width, even when the NSC dominates the latter. Experiments show that when the Stark tuning rate is large compared to the NSC (PhSeCN), the line width has a field dependence, albeit with major NSC-induced excursions from linearity. When the Stark tuning rate is small relative to the NSC (BZN), the line width is field-independent. BZN's line widths are substantially larger for the 1 → 2 transition, indicating a 1 → 2 transition enhancement of the NSC. Additionally, we examine, theoretically and experimentally, the difference in the 0 → 1 and 1 → 2 transitions' Stark tuning rates. Second-order perturbation theory combined with density functional theory explain the difference and show that the 1 → 2 transition's Stark tuning rate is ∼10% larger. The Stark tuning rate of PhSeCN is larger than BZN's for both transitions, consistent with the theoretical calculations. This study provides new insights into vibrational line shape components and a more general understanding of the vibrational response to external electric fields.

Time-resolved spectroscopic mapping of vibrational energy flow in proteins: Understanding thermal diffusion at the nanoscale

Vibrational energy exchange between various degrees of freedom is critical to barrier-crossing processes in proteins. Hemeproteins are well suited for studying vibrational energy exchange in proteins because the heme group is an efficient photothermal converter. The released energy by heme following photoexcitation shows migration in a protein moiety on a picosecond timescale, which is observed using time-resolved ultraviolet resonance Raman spectroscopy. The anti-Stokes ultraviolet resonance Raman intensity of a tryptophan residue is an excellent probe for the vibrational energy in proteins, allowing the mapping of energy flow with the spatial resolution of a single amino acid residue. This Perspective provides an overview of studies on vibrational energy flow in proteins, including future perspectives for both methodologies and applications.

Tuning Molecular Vibrational Energy Flow within an Aromatic Scaffold via Anharmonic Coupling

From guiding chemical reactivity in synthesis or protein folding to the design of energy diodes, intramolecular vibrational energy redistribution harnesses the power to influence the underlying fundamental principles of chemistry. To evaluate the ability to steer these processes, the mechanism and time scales of intramolecular vibrational energy redistribution through aromatic molecular scaffolds have been assessed by utilizing two-dimensional infrared (2D IR) spectroscopy. 2D IR cross peaks reveal energy relaxation through an aromatic scaffold from the azido- to the cyano-vibrational reporters in para-azidobenzonitrile (PAB) and para-(azidomethyl)benzonitrile (PAMB) prior to energy relaxation into the solvent. The rates of energy transfer are modulated by Fermi resonances, which are apparent by the coupling cross peaks identified within the 2D IR spectrum. Theoretical vibrational mode analysis allowed the determination of the origins of the energy flow, the transfer pathway, and a direct comparison of the associated transfer rates, which were in good agreement with the experimental results. Large variations in energy-transfer rates, approximately 1.9 ps for PAB and 23 ps for PAMB, illustrate the importance of strong anharmonic coupling, i.e., Fermi resonance, on the transfer pathways. In particular, vibrational energy rectification is altered by Fermi resonances of the cyano- and azido-modes allowing control of the propensity for energy flow.

FeCl3-Assisted Niobium-Catalyzed Cycloaddition of Nitriles and Alkynes: Synthesis of Alkyl- and Arylpyrimidines Based on Independent Functions of NbCl5 and FeCl3 Lewis Acids

NbCl₅-catalyzed [2 + 2 + 2] cycloaddition of nitriles with alkynes was used to synthesize pyrimidine derivatives. In this reaction, the use of individual Lewis acids, namely NbCl₅ and FeCl₃, is a key strategy for achieving the reaction using a catalytic amount of NbCl₅. The roles of the two Lewis acids were investigated using FT-IR spectroscopy. The results showed that NbCl₅ served as an efficient Lewis acid catalyst for nitrile activation, whereas FeCl₃ showed stronger Lewis acidity toward pyrimidines, releasing NbCl₅ into the catalytic cycle.

Ultrafast fluxional exchange dynamics in electrolyte solvation sheath of lithium ion battery

Lithium cation is the charge carrier in lithium-ion battery. Electrolyte solution in lithium-ion battery is usually based on mixed solvents consisting of polar carbonates with different aliphatic chains. Despite various experimental evidences indicating that lithium ion forms a rigid and stable solvation sheath through electrostatic interactions with polar carbonates, both the lithium solvation structure and more importantly fluctuation dynamics and functional role of carbonate solvent molecules have not been fully elucidated yet with femtosecond vibrational spectroscopic methods. Here we investigate the ultrafast carbonate solvent exchange dynamics around lithium ions in electrolyte solutions with coherent two-dimensional infrared spectroscopy and find that the time constants of the formation and dissociation of lithium-ion···carbonate complex in solvation sheaths are on a picosecond timescale. We anticipate that such ultrafast microscopic fluxional processes in lithium-solvent complexes could provide an important clue to understanding macroscopic mobility of lithium cation in lithium-ion battery on a molecular level.

Dynamics of rhenium photocatalysts revealed through ultrafast multidimensional spectroscopy

Rhenium catalysts have shown promise to promote carbon neutrality by reducing a prominent greenhouse gas, CO2, to CO and other starting materials. Much research has focused on identifying intermediates in the photocatalysis mechanism as well as time scales of relevant ultrafast processes. Recent studies have implemented multidimensional spectroscopies to characterize the catalyst's ultrafast dynamics as it undergoes the many steps of its photocycle. Two-dimensional infrared (2D-IR) spectroscopy is a powerful method to obtain molecular structure information while extracting time scales of dynamical processes with ultrafast resolution. Many observables result from 2D-IR experiments including vibrational lifetimes, intramolecular redistribution time scales, and, unique to 2D-IR, spectral diffusion, which is highly sensitive to solute-solvent interactions and motional dynamics. Spectral diffusion, a measure of how long a vibrational mode takes to sample its frequency space due to multiple solvent configurations, has various contributing factors. Properties of the solvent, the solute's structural flexibility, and electronic properties, as well as interactions between the solvent and solute, complicate identifying the origin of the spectral diffusion. With carefully chosen experiments, however, the source of the spectral diffusion can be unveiled. Within the context of a considerable body of previous work, here we discuss the spectral diffusion of several rhenium catalysts at multiple stages in the catalysis. These studies were performed in multiple polar liquids to aid in discovering the contributions of the solvent. We also performed electronic ground state 2D-IR and electronic excited state transient-2D-IR experiments to observe how spectral diffusion changes upon electronic excitation. Our results indicate that with the original Lehn catalyst in THF, relative to the ground state, the spectral diffusion slows by a factor of 3 in the equilibrated triplet metal-to-ligand charge transfer state. We attribute this slowdown to a decrease in dielectric friction as well as an increase in molecular flexibility. It is possible to partially simulate the charge transfer by altering the electron density moderately by adding electron donating or withdrawing substituents symmetrically to the bipyridine ligand. We find that unlike the significant electronic structure change induced by MLCT, such small substituent effects do not influence the spectral diffusion. A solvent study in THF, DMSO, and CH3CN found there to be an explicit solvent dependence that we can correlate to the solvent donicity, which is a measure of its nucleophilicity. Future studies focused on the solvent effects on spectral diffusion in the crucial photoinitiated state can illuminate the role the solvent plays in the catalysis.

Second hyperpolarizability of C-H, C-D, and C≡N stretch vibrations determined from computational Raman activities and a comparison with experiments

We have recently demonstrated that the second hyperpolarizability γ of a selected vibrational mode of a molecule can be determined by using the computational Raman activity against an internal standard with a known Raman γ value. This approach provides a convenient way for prediction of the γ magnitude of DOVE four wave mixing spectroscopy, an optical analogue to two-dimensional (2D) NMR. Here, by using the Hartree-Fock (HF) method, the density functional theory (DFT) method, and the second-order Møller-Plesset perturbation theory (MP2) method, we extend our early work from the less anharmonic region <2000 cm(-1) into the more anharmonic region >2000 cm(-1) covering C-H, C-D, and C≡N stretching modes of benzene, deuterated benzene, acetonitrile, deuterated acetonitrile, and tetrahydrofuran. The computed Raman γ values of these vibrational modes have been determined by using either the 992 cm(-1) Raman band of benzene or the compound's own Raman band (C-C stretch) around 800-1000 cm(-1) as an internal standard. In this more anharmonic region, the HF method with a larger basis set provides the best outputs and the predicted Raman γ values agree well with experimental values for most of the vibrational modes studied. By choosing a suitable method and basis set, this facile approach could be applied to a broader spectral range for Raman γ estimation of various materials.

Raman second hyperpolarizability determination using computational Raman activities and a comparison with experiments

Doubly vibrationally enhanced (DOVE) four-wave mixing spectroscopy, an optical analogue to 2D NMR, involves two infrared transitions and a Raman transition. The magnitude of the DOVE second hyperpolarizability γ (or third-order susceptibility χ((3))) can be theoretically estimated if the values of the dipolar moments of the two infrared transitions and the γ of the Raman transition are known. The Raman γ can be measured by using the four-wave mixing interferometric method or conventional Raman spectroscopy in the presence of an internal standard. In this work, we examine if one can use the Raman activity computed from density functional theory calculation to determine the Raman γ of selected vibrational modes of several samples including deuterated benzene, acetonitrile, tetrahydrofuran, and sodium benzoate aqueous solution. The 992 cm(-1) Raman band of benzene serves as an internal standard for organic solvents, and the 880 cm(-1) Raman band of hydrogen peroxide is for the aqueous solution sample with known γ values. We have found that the predicted Raman γ values from the computational Raman activities match experimental data reasonably well, suggesting a facile approach to predict the Raman γ of interested systems.

Probing ion/molecule interactions in aqueous solutions with vibrational energy transfer

Interactions between model molecules representing building blocks of proteins and the thiocyanate anion, a strong protein denaturant agent, were investigated in aqueous solutions with intermolecular vibrational energy exchange methods. It was found that thiocyanate anions are able to bind to the charged ammonium groups of amino acids in aqueous solutions. The interactions between thiocyanate anions and the amide groups were also observed. The binding affinity between the thiocyanate anion and the charged amino acid residues is about 20 times larger than that between water molecules and the amino acids and about 5-10 times larger than that between the thiocyanate anion and the neutral backbone amide groups. The series of experiments also demonstrates that the chemical nature, rather than the macroscopic dielectric constant, of the ions and molecules plays a critical role in ion/molecule interactions in aqueous solutions.

Measurement of Raman χ(3) and theoretical estimation of DOVE four wave mixing of hydrogen peroxide

Hydrogen peroxide has strong infrared (IR) transitions ν(6) and its combination band ν(2)+ν(6), which may provide a unique opportunity to implement doubly vibrationally enhanced (DOVE) four wave mixing (FWM) for directly measuring hydrogen peroxide in spectrally overcrowded mixtures. In this work, the magnitude of the DOVE third-order susceptibility χ(3) was theoretically estimated. By using a FWM interferometric method, one of the strongest Raman bands, O-O stretch ν(3) Raman χ(3) of 30 wt % H(2)O(2), was first measured to be 1.2 × 10(-14) esu. The Raman χ(3) of ν(2) was then determined to be 5.3 × 10(-15) esu based on their relative Raman intensities. The resulting Raman χ(3) of ν(2) was used to calculate the DOVE χ(3) of (ν(6), ν(2)+ν(6)), together with the dipolar moments of the two IR transitions determined from IR absorption measurement. The calculated value of DOVE-IR χ(3) was 1.1 × 10(-13) esu for pure H(2)O(2), about 1.5 times larger than that of the strong ring breathing Raman band of benzene. The large DOVE χ(3) suggests the feasibility of direct measurement of hydrogen peroxide in an aqueous environment using DOVE four wave mixing.