Contributions of rotation and translation to polarizability anisotropy and solvation dynamics in acetonitrile

Collective Spectroscopy of Solvation Phenomena: Conflicts, Challenges, and Opportunities

Different spectroscopy types reveal different aspects of molecular processes in soft matter. In particular, collective observables can provide insights into intermolecular correlations invisible to the more popular single-particle methods. In this perspective we feature the dielectric relaxation spectroscopy (DRS) with an emphasis on the proper interpretation of this complex observable aided by computational spectroscopy. While we focus on the history and recent advances of DRS in the fields of biomolecular hydration and nanoconfinement, the discussion transcends this particular field and provides a guide for how collective spectroscopy types supported by computational decomposition can be employed to further our understanding of soft matter phenomena.

Altered relaxation dynamics of excited state reactions by confinement in reverse micelles probed by ultrafast fluorescence up-conversion

Chemical reactions in confined environments are important in areas as diverse as heterogenous catalysis, environmental chemistry and biochemistry, yet they are much less well understood than the equivalent reactions in either the gas phase or in free solution. The understanding of chemical reactions in solution was greatly enhanced by real time studies of model reactions, through ultrafast spectroscopy (especially when supported by molecular dynamics simulation). Here we review some of the efforts that have been made to adapt this approach to the investigation of reactions in confined media. Specifically, we review the application of ultrafast fluorescence spectroscopy to measure reaction dynamics in the nanoconfined water phase of reverse micelles, as a function of the droplet radius and the charge on the interface. Methods of measurement and modelling of the reactions are outlined. In all of the cases studied (which are focused on ultrafast intramolecular reactions) the effect of confinement was to suppress the reaction. Even in the largest micelles the result in the bulk aqueous phase was not usually recovered, suggesting an important role for specific interactions between reactant and environment, for example at the interface. There was no simple one-to-one correspondence with direct measures of the dynamics of the confined phase. Thus, understanding the effect of confinement on reaction rate appears to require not only knowledge of the dynamics of the reaction in solutions and the effect of confinement on the medium, but also of the interaction between reactant and confining medium.

Optical Kerr Effect of Liquid Acetonitrile Probed by Femtosecond Time-Resolved X-ray Liquidography

Optical Kerr effect (OKE) spectroscopy is a method that measures the time-dependent change of the birefringence induced by an optical laser pulse using another optical laser pulse and has been used often to study the ultrafast dynamics of molecular liquids. Here we demonstrate an alternative method, femtosecond time-resolved X-ray liquidography (fs-TRXL), where the microscopic structural motions related to the OKE response can be monitored using a different type of probe, i.e., X-ray solution scattering. By applying fs-TRXL to acetonitrile and a dye solution in acetonitrile, we demonstrate that different types of molecular motions around photoaligned molecules can be resolved selectively, even without any theoretical modeling, based on the anisotropy of two-dimensional scattering patterns and extra structural information contained in the q-space scattering data. Specifically, the dynamics of reorientational (libration and orientational diffusion) and translational (interaction-induced motion) motions are captured separately by anisotropic and isotropic scattering signals, respectively. Furthermore, the two different types of reorientational motions are distinguished from each other by their own characteristic scattering patterns and time scales. The measured time-resolved scattering signals are in excellent agreement with the simulated scattering signals based on a molecular dynamics simulation for plausible molecular configurations, providing the detailed structural description of the OKE response in liquid acetonitrile.

Low-Frequency Vibrational Motions of Polystyrene in Carbon Tetrachloride: Comparison with Model Monomer and Dependence on Concentration and Molecular Weight

In this study, the low-frequency vibrational dynamics of polystyrene (PS) in CCl₄ was investigated by femtosecond Raman-induced Kerr effect spectroscopy. Ethylbenzene (EBz) was also investigated as a model monomer of the polymer to elucidate the unique dynamical features of PS in solution. The broadened low-frequency spectrum of the PS/CCl₄ in the frequency region below 150 cm^-1 is significantly different from that of the EBz/CCl₄. Difference spectra between the PS or EBz solutions and neat CCl₄, normalized to an internal vibrational mode of CCl₄, clearly show a much lower spectral intensity for the PS/CCl₄ than the EBz/CCl₄ in the low-frequency region below ca. 20 cm^-1. This indicates that translational motions are suppressed in the PS/CCl₄ compared to the EBz/CCl₄. Moreover, the high-frequency motion at ca. 70 cm^-1, mainly due to phenyl ring librations, occurs at higher frequency in PS (78 cm^-1) than EBz (65 cm^-1). In addition, the results of concentration-dependent experiments show that the first moment (M₁) of the low-frequency difference spectra of both PS/CCl₄ and EBz/CCl₄ is almost independent of the concentration. The molecular weight dependence of the low-frequency spectrum in the PS/CCl₄ shows that the M₁ value of the low-frequency spectral band of PS shifts to higher frequencies when the molecular weight of PS increases up to M_w = ∼1000, which corresponds approximately to the decamer, and then remains constant upon further increasing the molecular weight.

Atomistic characterization of the active-site solvation dynamics of a model photocatalyst

The interactions between the reactive excited state of molecular photocatalysts and surrounding solvent dictate reaction mechanisms and pathways, but are not readily accessible to conventional optical spectroscopic techniques. Here we report an investigation of the structural and solvation dynamics following excitation of a model photocatalytic molecular system [Ir₂(dimen)₄]²⁺, where dimen is para-diisocyanomenthane. The time-dependent structural changes in this model photocatalyst, as well as the changes in the solvation shell structure, have been measured with ultrafast diffuse X-ray scattering and simulated with Born-Oppenheimer Molecular Dynamics. Both methods provide direct access to the solute–solvent pair distribution function, enabling the solvation dynamics around the catalytically active iridium sites to be robustly characterized. Our results provide evidence for the coordination of the iridium atoms by the acetonitrile solvent and demonstrate the viability of using diffuse X-ray scattering at free-electron laser sources for studying the dynamics of photocatalysis.

Interactions between reactive excited states of molecular photocatalysts and surrounding solvent can dictate reaction pathways, but are not readily accessible to conventional spectroscopic methods. Here the authors use diffuse X-ray scattering and theory to study the atomistic solvation dynamics of a photoexcited di-iridium complex in acetonitrile.

Local structural effects on orientational relaxation of OH-bond in liquid water over short to intermediate timescales

By simulating the rigid simple point charge extended model at temperature T = 300 K, the orientational relaxation of the OH-bond in water was investigated over short to intermediate timescales, within which molecules undergo inertial rotation and libration and then enter the rotational diffusion regime. According to the second-cumulant approximation, the orientational time correlation function (TCF) of each axis that is parallel or perpendicular to an OH-bond is related to an effective rotational density of states (DOS), which is determined using the power spectra of angular velocity autocorrelation functions (AVAFs) of the other two axes. In addition, the AVAF power spectrum of an axis was approximated as the rotational stable instantaneous normal mode (INM) spectrum of the axis. As described in a previous study [S. L. Chang, T. M. Wu, and C. Y. Mou, J. Chem. Phys. 121, 3605 (2004)], simulated molecules were classified into subensembles, according to either the local structures or the H-bond configurations of the molecules. For global molecules and the classified subensembles, the simulation results for the first- and second-rank orientational TCFs were compared with the second-cumulant predictions obtained using the effective rotational DOSs and the rotational stable-INM spectra. On short timescales, the OH-bond in water behaves similar to an inertial rotor and its anisotropy is lower than that of a water molecule. For molecules with three or more H-bonds, the OH-bond orientational TCFs are characterized by a recurrence, which is an indication for libration of the OH-bond. The recurrence can generally be described by the second-cumulant prediction obtained using the rotational stable-INM spectra; however, the orientational TCFs after the recurrence switch to a behavior similar to that predicted using the AVAF power spectra. By contrast, the OH-bond orientational TCFs of molecules initially connected with one or two H-bonds decay monotonically or exhibit a weak recurrence, indicating rapid relaxation into the rotational diffusion regime after the initial Gaussian decay. In addition to accurately describing the Gaussian decay, the second-cumulant predictions formulated using the rotational stable-INM spectra and the AVAF power spectra serve as the upper and lower limits, respectively, for the OH-bond orientational TCFs of these molecules after the Gaussian decay.

Hybrid quantum/classical simulations of the vibrational relaxation of the amide I mode of N-methylacetamide in D2O solution

Hybrid quantum/classical molecular dynamics (MD) is applied to simulate the vibrational relaxation (VR) of the amide I mode of deuterated N-methylacetamide (NMAD) in aqueous (D(2)O) solution. A novel version of the vibrational molecular dynamics with quantum transitions (MDQT) treatment is developed in which the amide I mode is treated quantum mechanically while the remaining degrees of freedom are treated classically. The instantaneous normal modes of the initially excited NMAD molecule (INM(0)) are used as internal coordinates since they provide a proper initial partition of the system in quantum and classical subsystems. The evolution in time of the energy stored in each individual normal mode is subsequently quantified using the hybrid quantum-classical instantaneous normal modes (INM(t)). The identities of both the INM(0)s and the INM(t)s are tracked using the equilibrium normal modes (ENMs) as templates. The results extracted from the hybrid MDQT simulations show that the quantum treatment of the amide I mode accelerates the whole VR process versus pure classical simulations and gives better agreement with experiments. The relaxation of the amide I mode is found to be essentially an intramolecular vibrational redistribution (IVR) process with little contribution from the solvent, in agreement with previous theoretical and experimental studies. Two well-defined relaxation mechanisms are identified. The faster one accounts for ≈40% of the total vibrational energy that flows through the NMAD molecule and involves the participation of the lowest frequency vibrations as short-life intermediate modes. The second and slower mechanism accounts for the remaining ≈60% of the energy released and is associated to the energy flow through specific mid-range and high-frequency modes.

Molecular Dynamics Study of Polarizability Anisotropy Relaxation in Aromatic Liquids and Its Connection with Local Structure

The collective polarizability anisotropy dynamics in a set of three aromatic liquids, benzene (Bz), hexafluorobenzene (HFB), and 1,3,5-trifluorobenzene (TFB), has been studied by molecular dynamics simulation. These liquids have very similar shapes, but different electrostatic interactions due to opposite polarities of C-H and C-F bonds, giving rise to different local intermolecular structures in the liquid phase. We have investigated how these structural arrangements affect polarizability anisotropy dynamics observed in optical Kerr-effect (OKE) spectroscopy. We have modeled the interaction-induced polarizability with the first-order dipole-induced dipole approximation, with the molecular polarizability distributed over the carbon sites. Local contributions to the librational OKE spectrum were computed separately for molecules participating in parallel or perpendicular relative orientations within the first coordination shell. We found that the relative locations of parallel and perpendicular librational bands of the OKE spectra are closely related to the corresponding pair energy distributions of the closest four neighbors of a given molecule, corresponding to a model of a harmonic oscillator in a cage of nearest neighbors. This model predicts higher librational frequencies for more attractive intermolecular interactions, which in all three liquids correspond to parallel local arrangements. On the diffusive orientational time scale, all three liquids exhibit slower relaxation of molecules in parallel arrangements, although the difference in relaxation rates is substantial only in TFB, which has the strongest tendency toward parallel stacking. The analysis of the collective polarizability relaxation was performed using two different approaches, the projection scheme (J. Chem. Phys. 1980, 72, 2801) and the theory developed by Steele (Mol. Phys. 1987, 61, 1031) for the second time derivatives applied to collective time correlations. Both approaches allow the decomposition of the OKE response into contributions from orientational relaxation and other dynamical processes. We find that they lead to different predictions on how the response depends on collective reorientation and processes arising from fluctuations in the interaction-induced polarizability. We discuss the reasons for these differences and the advantages and disadvantages of the two analysis schemes.

Solvation in supercritical water

Solvation in supercritical water under equilibrium and nonequilibrium conditions is studied via molecular dynamics simulations. The influence of solute charge distributions and solvent density on the solvation structures and dynamics is examined with a diatomic probe solute molecule. It is found that the solvation structure varies dramatically with the solute dipole moment, especially in low-density water, in accord with many previous studies on ion solvation. This electrostrictive effect has important consequences for solvation dynamics. In the case of a nonequilibrium solvent relaxation, if there are sufficiently many water molecules close to the solute at the outset of the relaxation, the solvent response measured as a dynamic Stokes shift is almost completely governed by inertial rotations of these water molecules. By contrast, in the opposite case of a low local solvent density near the solute, not only rotations but also translations of water molecules play an important role in solvent relaxation dynamics. The applicability of a linear response is found to be significantly restricted at low water densities.

Temperature- and solvation-dependent dynamics of liquid sulfur dioxide studied through the ultrafast optical Kerr effect

The ultrafast dynamics of liquid sulphur dioxide have been studied over a wide temperature range and in solution. The optically heterodyne-detected and spatially masked optical Kerr effect (OKE) has been used to record the anisotropic and isotropic third-order responses, respectively. Analysis of the anisotropic response reveals two components, an ultrafast nonexponential relaxation and a slower exponential relaxation. The slower component is well described by the Stokes-Einstein-Debye equation for diffusive orientational relaxation. The simple form of the temperature dependence and the agreement between collective (OKE) and single molecule (e.g., NMR) measurements of the orientational relaxation time suggests that orientational pair correlation is not significant in this liquid. The relative contributions of intermolecular interaction-induced and single-molecule orientational dynamics to the ultrafast part of the spectral density are discussed. Single-molecule librational-orientational dynamics appear to dominate the ultrafast OKE response of liquid SO2. The temperature-dependent OKE data are transformed to the frequency domain to yield the Raman spectral density for the low-frequency intermolecular modes. These are bimodal with the lowest-frequency component arising from diffusive orientational relaxation and a higher-frequency component connected with the ultrafast time-domain response. This component is characterized by a shift to higher frequency at lower temperature. This result is analyzed in terms of a harmonic librational oscillator model, which describes the data accurately. The observed spectral shifts with temperature are ascribed to increasing intermolecular interactions with increasing liquid density. Overall, the dynamics of liquid SO2 are found to be well described in terms of molecular orientational relaxation which is controlled over every relevant time range by intermolecular interactions.

Why Does the Intermolecular Dynamics of Liquid Biphenyl so Closely Resemble that of Liquid Benzene? Molecular Dynamics Simulation of the Optical-Kerr-Effect Spectra

The combination of optical-Kerr-effect (OKE) spectroscopy and molecular dynamics simulations has provided us with a newfound ability to delve into the librational dynamics of liquids, revealing, in the process, some surprising commonalities among aromatic liquids. Benzene and biphenyl, for example, have remarkably similar OKE spectra despite marked differences in their shapes, sizes, and moments of inertia--and even more chemically distinct aromatics tend to have noticeable similarities in their spectra. We explore this universality by using a molecular dynamics simulation to investigate the librational dynamics of molten biphenyl and to predict its OKE spectrum, comparing the results with our previous calculations for liquid benzene. We suggest that the impressive level of quantitative agreement between these two liquids is largely a reflection of the fact that librations in these and other aromatic liquids act as torsional oscillations with oscillator frequencies selected from the liquid's librational bands. Since these bands are centered about the librational Einstein frequencies, the quantitative similarities between the liquids are essentially reflections of the near identities of their Einstein frequencies. Why then are the Einstein frequencies themselves so insensitive to molecular details? We show that, for nearly planar molecules, mean-square torques and moments of inertia tend to scale with molecular dimensions in much the same way. We demonstrate that this near cancellation provides both a quantitative explanation of the close relationship between benzene and biphenyl and a more general perspective on the similarities seen in the ultrafast dynamics of aromatic liquids.

Polarizability response in polar solvents: Molecular-dynamics simulations of acetonitrile and chloroform

The relaxation of the many-body polarizability in liquid acetonitrile and chloroform at room temperature was studied by molecular-dynamics simulations. The collective polarizability induced by intermolecular interactions was included using first- and all-orders dipole-induced-dipole models and calculated considering both molecule-centered and distributed site polarizabilities. The anisotropic response was analyzed using a separation scheme that allows a decomposition of the total response in terms of orientational and collision-induced effects. We found the method effective in approximately separating the contributions of these relaxation mechanisms, although the orientational-collision-induced interference makes a non-negligible contribution to the total response. In both liquids the main contribution to the anisotropic response is due to orientational dynamics, but intermolecular collision-induced (or translational) effects are important, especially at short times. We found that higher-order interaction-induced effects were essentially negligible for both liquids. Larger differences were found between the center-center and site-site models, with the latter showing faster polarizability relaxation and better agreement with experiment. Isotropic and anisotropic spectra were computed from the corresponding time correlation functions. The lowest-frequency contributions are largely suppressed in the isotropic spectra and their overall shape is similar to the purely collision-induced contribution to the anisotropic spectra, but with an amplitude which is smaller by a factor of approximately 5 in acetonitrile and approximately 3 in chloroform.

Measuring the Change in the Intermolecular Raman Spectrum during Dipolar Solvation

We demonstrate a method to directly measure the change in the spectrum of intermolecular solvent fluctuations as a function of time after electronic excitation of a solute, and this method is applied to the dye Coumarin 102 (C102) in acetonitrile. The complete intermolecular response is captured following resonant excitation with time domain third-order Raman spectroscopy. In a previous report, we introduced this method and used it to probe one point in the intermolecular response as a function of time after solute excitation (Underwood, D. F., Blank, D. A. J. Phys. Chem. A 2003, 107 (7), 956). Here we extend this approach to recover the change in the entire intermolecular response as a function of time. To our knowledge the results provide the first direct measurement of the difference in the equilibrated intermolecular response after excitation of a solute and its evolution during a dipolar solvation event. Excitation of C102 results in a significant increase in the solvent-solute interaction due to a large increase in the dipole moment. The observed change in the intermolecular response is consistent with a rapid change in local solvent density, with intermolecular kinetic energy transfer changing the response on longer time scales. Evolution of the response exhibits a strong frequency dependence and suggests changes over longer distances at longer delay times. The measured change in the spectrum of solvent fluctuations represents a direct experimental confirmation of the breakdown of linear response and confirms predictions from molecular dynamics simulations.

Solvation Dynamics in Acetonitrile: A Study Incorporating Solute Electronic Response and Nuclear Relaxation

The solvent reorganization process after electronic excitation of a polar solute in a polar solvent such as acetonitrile is related mainly to the time evolution of the solute-solvent electrostatic interaction. Modern laser-based techniques have sufficient time resolution to follow this decay in real time, providing information to be confirmed and interpreted by theories and models. We present here a study aimed at the investigation of the different steps involved in the process taking place after a vertical S(0) --> S(1) excitation of a large size chromophore, coumarin 153 (C153), in acetonitrile, from both the solute and the solvent points of view. To do this, we use accurate quantum mechanical calculations for the solute properties within the polarizable continuum model (PCM) and classical molecular dynamics (MD) simulations, both equilibrium and nonequilibrium, for C153 in the presence of the solvent. The geometry of the solute is allowed to change in order to study the role of internal motions in the time-dependent solvation process. The solvent response function has been obtained from the simulation data and compared to experiment, while the comparison between equilibrium and nonequilibrium MD results for the solvation response confirms the validity of the linear response approximation in the C153-acetonitrile system. The MD trajectories have also been used to monitor the structure of the solvation shell and to determine its change in response to the change in the solute partial charges.