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

Low-Frequency Spectra of Hydrated Ionic Liquids with Kosmotropic and Chaotropic Anions

In this study, we investigated the water concentration dependence of the intermolecular vibrations of two hydrated ionic liquids (ILs), cholinium dihydrogen phosphate ([ch][dhp]) and cholinium bromide ([ch]Br), using femtosecond Raman-induced Kerr effect spectroscopy (fs-RIKES). The anions of the former and latter hydrated ILs are kosmotropic and chaotropic, respectively. We found that the spectral peak of ∼50 cm-1 shifted to the low-frequency side in hydrated [ch][dhp], indicating the weakening of its intermolecular interactions. In contrast, no change in the peak frequency of the low-frequency band at ∼50 cm-1 was observed with increasing water concentration in hydrated [ch]Br. The vibrational density of states (VDOS) spectra generated from molecular dynamics (MD) simulations were in qualitative agreement with the experimental results. Decomposition analysis of the VDOS spectra for each component revealed that the red shift of the low-frequency band in the hydrated [ch][dhp] upon water addition was essentially due to the contributions of anions and water rather than that of the cholinium cation. We also found from the low-frequency spectra of the two hydrated ILs that they differed in the concentration dependence of the 180 cm-1 band, which is assigned as a hindered translational motion of water molecules combined to form O···O stretching motions. From the relationship between the peak frequency of the low-frequency band and the bulk parameter, which is the square root of the surface tension divided by the density, we found that the peak frequency in the hydrated IL with kosmotropic [dhp]- depends on the bulk parameter, similar to the case for an aqueous solution of the typical deep eutectic solvent reline. However, the peak frequency of the hydrated IL with chaotropic Br- is constant with the bulk parameter.

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

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

The polarizability response of a glass-forming liquid reveals intrabasin motion and interbasin transitions on a potential energy landscape

Potential energy landscape (PEL) concepts have been useful in conceptualizing the effects of intermolecular interactions on dynamic and thermodynamic properties of liquids and glasses. "Basins", or regions of reduced potential energy associated with locally preferred molecular packing are important PEL features. The molecular configurations at the bottom of these basins are referred to as inherent structures (ISs). Experimental methods for directly characterizing PEL features such as these are rare, largely relegating PEL concepts to theory and simulation studies, and impeding their exploration in real systems. Recently, we showed that quasielastic neutron scattering (QENS) data from propylene carbonate (PC) exhibit signatures of picosecond timescale motion that are consistent with intrabasin motion and interbasin transitions [Cicerone et al., J. Chem. Phys., 2017, 146, 054502]. Here we present optically-heterodyne-detected optical Kerr effect (OHD-OKE) spectroscopy studies on PC. The data exhibit signatures of motion within and transitions between basins that agree quantitatively with and extend the QENS results. We show that the librational component of the OKE response corresponds to intrabasin dynamics, and the enigmatic intermediate OKE response corresponds to interbasin transition events. The OKE data extend the measurement range of these parameters and reveal their utility in characterizing PEL features of real systems.

Empirical Analysis of Optical Kerr Effect Spectra: A Case for Constraint

Ultrafast optical Kerr effect (OKE) spectroscopy is a widely used method for studying the depolarized, Raman-active intermolecular dynamics of liquids. Through appropriate manipulation of OKE data, it is possible to determine the reduced spectral density (RSD), which is the Bose-Einstein-corrected, low-frequency Raman spectrum with the contribution of diffusive reorientation removed. OKE RSDs for van der Waals liquids can often be fit well to an empirical function that is the sum of a Bucaro-Litovitz function and an antisymmetrized Gaussian (AG). Although these functions are not directly representative of specific intermolecular dynamics, the AG fit parameters can provide useful insights into the microscopic properties of liquids. Here we show that fits using the AG function are typically not well-determined, and that equally good results can be obtained with a wide range of fitting parameters. We propose the use of a physically motivated constraint on the amplitude of the AG function, and demonstrate that this constraint leads to more intuitive trends in the fit parameters for temperature-dependent RSDs in 1,3,5-trifluorobenzene and hexafluorobenzene.

Comparative study of the intermolecular dynamics of imidazolium-based ionic liquids with linear and branched alkyl chains: OHD-RIKES measurements

This article describes a comparative study of the low-frequency (0-450 cm^-1) Kerr spectra of the branched 1-(iso-alkyl)-3-methylimidazolium bis(trifluoromethylsulfonyl)amide ([(N - 2)mC_N-1C₁im][NTf₂] with N = 3-7) ionic liquids (ILs) and that of the linear 1-(n-alkyl)-3-methylimidazolium bis(trifluoromethylsulfonyl)amide ([C_NC₁im][NTf₂] with N = 2-7) ILs. The spectra were obtained by use of femtosecond optical heterodyne-detected Raman-induced Kerr effect spectroscopy (OHD-RIKES). The intermolecular spectrum of a branched IL is similar to that of a linear IL that is of the same alkyl chain length rather than of the same number of carbon atoms in the alkyl chain. This similarity and the lack of a correlation of the first spectral moments and widths of the intermolecular spectra with chain length is mainly attributed to the increase in the dispersion contribution to the total molar cohesive energy being compensated by stretching of the ionic network due to the increasing size of the nonpolar domains, which is dependent only on the length of the alkyl chain.

The importance of polarizability: comparison of models of carbon disulphide in the ionic liquids [C1C1im][NTf2] and [C4C1im][NTf2]

Three dimensional distribution of CS₂ around a [C₁C₁im]⁺ ion showing the difference in behaviour of polarizable (red) and unpolarizable (blue) models of CS₂.

Assessing the role of moment of inertia in optical Kerr effect spectroscopy

Optical Kerr effect (OKE) spectroscopy allows for the acquisition of high-quality, Bose-Einstein-corrected, low-frequency Raman spectra in liquids. However, the assignment of a molecular interpretation to these spectra remains an open problem. To address this issue, here we present an OKE study of benzene and four of its isotopologues. Our results indicate that hindered rotations are the major contributor to the OKE reduced spectral density (RSD) over the entire intermolecular spectral region (0-250 cm(-1)). We also have found that isotopologues with six (13)C atoms have RSDs that are enhanced at frequencies below 30 cm(-1). We further demonstrate that the collective orientational correlation time of these liquids scales with the inverse square root of the tumbling moment of inertia, indicating that there is strong translation-rotation coupling in benzene.

Nanostructural Organization and Anion Effects on the Temperature Dependence of the Optical Kerr Effect Spectra of Ionic Liquids

The intermolecular spectra of three imidazolium ionic liquids were studied as a function of temperature by the use of optical heterodyne-detected Raman-induced Kerr effect spectroscopy. The ionic liquids comprise the 1,3-pentylmethylimidazolium cation ([C(5)mim]+), and the anions, bromide (Br-), hexafluorophosphate (PF(6)-), and bis(trifluoromethanesulfonyl)imide (NTf(2)-). Whereas the optical Kerr effect (OKE) spectrum of [C(5)mim][NTf(2)] is temperature-dependent, the OKE spectra of [C(5)mim]Br and [C(5)mim][PF6] are temperature-independent. These results are surprising in light of the fact that the bulk densities of these room temperature ionic liquids (RTILs) are temperature-dependent. The temperature independence of the OKE spectra and the temperature dependence of the bulk density in [C(5)mim]Br and [C(5)mim][PF(6)] suggest that there are inhomogeneities in the densities of these liquids. The existence of density inhomogeneities is consistent with recent molecular dynamics simulations that show RTILs to be nanostructurally organized with nonpolar regions arising from clustering of the alkyl chains and ionic networks arising from charge ordering of the anions and imidazolium rings of the cations. Differences in the temperature dependences of the OKE spectra are rationalized on the basis of the degree of charge ordering in the polar regions of the RTILs.

Ultrafast dynamics in complex fluids observed through the ultrafast optically-heterodyne-detected optical-Kerr-effect (OHD-OKE)

The ultrafast molecular dynamics of complex fluids have been recorded using the optically-heterodyne-detected optical-Kerr-effect (OHD-OKE). The OHD-OKE method is reviewed and some recent refinements to the method are described. Applications to a range of complex fluids, including microemulsions, polymer melts and solutions, liquid crystal and ionic liquids are surveyed. The level of detail attainable with the OHD-OKE method in these complex fluids is discussed. The prospects for future experiments are discussed.