Infrared and Raman spectra of silica polymorphs from an ab initio parametrized polarizable force field

The general aim of this study is to test the reliability of polarizable model potentials for the prediction of vibrational (infrared and Raman) spectra in highly anharmonic systems such as high temperature crystalline phases. By using an ab initio parametrized interatomic potential for SiO2 and molecular dynamics simulations, we calculate the infrared and Raman spectra for quartz, cristobalite, and stishovite at various thermodynamic conditions. The model is found to perform very well in the prediction of infrared spectra. Raman peak positions are also reproduced very well by the model; however, Raman intensities calculated by explicitly taking the derivative of the polarizability with respect to the atomic displacements are found to be in poorer agreement than intensities calculated using a parametrized "bond polarizability" model. Calculated spectra for the high temperature beta phases, where the role of dynamical disorder and anharmonicities is predominant, are found to be in excellent agreement with experiments. For the octahedral phases, our simulations are able to reproduce changes in the Raman spectra across the rutile-to-CaCl2 transition around 50 GPa, including the observed phonon softening.     

Electronic polarization effect on low-frequency infrared and Raman spectra of aprotic solvent: molecular dynamics simulation study with charge response kernel by second order Møller-Plesset perturbation method

Low-frequency infrared (IR) and depolarized Raman scattering (DRS) spectra of acetonitrile, methylene chloride, and acetone liquids are simulated via molecular dynamics calculations with the charge response kernel (CRK) model obtained at the second order M?ller-Plesset perturbation (MP2) level. For this purpose, the analytical second derivative technique for the MP2 energy is employed to evaluate the CRK matrices. The calculated IR spectra reasonably agree with the experiments. In particular, the agreement is excellent for acetone because the present CRK model well reproduces the experimental polarizability in the gas phase. The importance of interaction induced dipole moments in characterizing the spectral shapes is stressed. The DRS spectrum of acetone is mainly discussed because the experimental spectrum is available only for this molecule. The calculated spectrum is close to the experiment. The comparison of the present results with those by the multiple random telegraph model is also made. By decomposing the polarizability anisotropy time correlation function to the contributions from the permanent, induced polarizability and their cross term, a discrepancy from the previous calculations is observed in the sign of permanent-induce cross term contribution. The origin of this discrepancy is discussed by analyzing the correlation functions for acetonitrile.     

Two-dimensional Raman spectra of atomic solids and liquids

We calculate third- and fifth-order Raman spectra of simple atoms interacting through a soft-core potential by means of molecular-dynamics (MD) simulations. The total polarizability of molecules is treated by the dipole-induced dipole model. Two- and three-body correlation functions of the polarizability at various temperatures are evaluated from equilibrium MD simulations based on a stability matrix formulation. To analyze the processes involved in the spectroscopic measurements, we divide the fifth-order response functions into symmetric and antisymmetric integrated response functions; the symmetric one is written as a simple three-body correlation function, while the antisymmetric one depends on a stability matrix. This analysis leads to a better understanding of the time scales and molecular motions that govern the two-dimensional (2D) signal. The 2D Raman spectra show novel differences between the solid and liquid phases, which are associated with the decay rates of coherent motions. On the other hand, these differences are not observed in the linear Raman spectra.     

Molecular dynamics study of the temperature-dependent Optical Kerr effect spectra and intermolecular dynamics of room temperature ionic liquid 1-methoxyethylpyridinium dicyanoamide

We have performed classical molecular dynamics simulations to calculate the Optical Kerr effect (OKE) spectra of 1-methoxyethylpyridinium dicyanoamide, a room-temperature ionic liquid (IL) which has been recently studied by Shirota and Castner (Shirota, H. ; Castner, E. J. Phys. Chem. A 2005, 109, 9388-9392) in comparison to its neutral isoelectronic solvent mixture. Our theoretical and computational studies show that the decay of the collective polarizability anisotropy correlation exhibits several different time scales originating from inter- and intramolecular dynamics, in good agreement with experiments. What's more, we find that the portion of the collective anisotropic polarizability relaxation due to "interaction-induced" phenomena is important at times much longer than those observed in normal solvents when these are far from their glass transition temperature. From our long (60 ns) molecular dynamics simulations, we are able to determine the appropriate time scales for orientational relaxation and interaction-induced processes occurring in the liquid. We find that the cationic contribution to the OKE signal is predominant. Because of the slow nature of relaxation processes in ILs, these calculations are very time, memory, and storage intensive. In the context of this research, we have developed a polarizable force field for this system and also theoretical methodology to generate molecular polarizabilities for arbitrarily shaped molecules and ions from corresponding atomic polarizabilities. We expect this methodology to have an important impact on the speed of molecular dynamics simulations of polarizable systems in the future.     

Infrared and Raman line shapes of dilute HOD in liquid H2O and D2O from 10 to 90 degrees C

A combined electronic structure/molecular dynamics approach was used to calculate infrared and isotropic Raman spectra for the OH or OD stretches of dilute HOD in D2O or H2O, respectively. The quantities needed to compute the infrared and Raman spectra were obtained from density functional theory calculations performed on clusters, generated from liquid-state configurations, containing an HOD molecule along with 4-9 solvent water molecules. The frequency, transition dipole, and isotropic transition polarizability were each empirically related to the electric field due to the solvent along the OH (or OD) bond, calculated on the H (or D) atom of interest. The frequency and transition dipole moment of the OH (or OD) stretch of the HOD molecule were found to be very sensitive to its instantaneous solvent environment, as opposed to the isotropic transition polarizability, which was found to be relatively insensitive to environment. Infrared and isotropic Raman spectra were computed within a molecular dynamics simulation by using the empirical relationships and semiclassical expressions for the line shapes. The line shapes agree well with experiment over a temperature range from 10 to 90 degrees C.     

Theoretical study of the local structure and Raman spectra of CaO-SiO2 binary melts

A procedure for the Raman spectra calculation of vitreous and molten silicates was presented in this paper. It includes molecular dynamics MD simulation for the generation of equilibrium configurations, Wilson's GF matrix method for the calculations of eigenfrequencies and corresponding vectors, electro-optical parameters method (EOPM) for the Raman intensity calculations, and the bond polarizability model (BPM) for the determination of polarizability and polarizability derivative. One of the most important characteristics of this procedure is the achievement of the partial Raman spectra of five tetrahedral units, as well as the total spectral envelope. In this paper, the calculation was carried out for the vitreous and molten calcium silicates with different compositions and at various temperatures. It is worthwhile to note that the calculation is based on statistical configurations distribution in the space and so it is not needed to artificially adjust the full width at half maximum (FWHM) of spectra. It was also tested through the good agreement of the calculated spectra with the experimental, including some regularity of spectral properties. According to the calculation, the symmetrical stretching of whole tetrahedral units, to which the stretching of Si-O(nb) bond gives the main contribution to intensity, is proven to be the dominance in the high-frequency range (800-1200 cm(-1)) and the symmetrical bending of Si-O(b)-Si, to which the stretching of Si-O(b) bond exhibits the main contribution, is the dominance in the medium-frequency range (400-700 cm(-1)). As the first theoretical results, the Raman scattering coefficient of each Q(i) was found little change along with the variation of composition and temperature.     

Raman spectra of ionic liquids: interpretation via computer simulation

Theoretical Raman spectra of the complex-forming ionic liquids LaCl3 and ScCl3, derived from molecular dynamics computer simulations, are presented. These simulations, which use polarizable ion interaction models, have previously been shown to predict structural properties in excellent agreement with diffraction experiments. The dependence of the polarizability of the melt on the ionic positions, which determines the Raman spectrum through the time dependence of the polarizability correlation function, is modeled on the basis of ab initio electronic structure calculations carried out on alkali chlorides. New simulation techniques are introduced in order to allow the spectrum to be calculated with acceptable statistics. The calculated spectra are in semiquantitative agreement with experimental data. The distinctive bands which appear in the spectra of such complex melts are linked to the vibrations of the transient coordination complexes which form in these melts and new interpretations for the origin of several well-known features are proposed. The simulations thus enable a link between the structure of a melt as perceived through Raman spectroscopy and through diffraction experiments to be made.     

Probing excited-state dynamics and intramolecular proton transfer in 1-acylaminoanthraquinones via the intermolecular solvent response

The dynamics of a series of 1-acylaminoanthraquinones with varying degrees of excited-state intramolecular proton transfer are studied in acetonitrile and dichloromethane. Events are followed via changes in the third-order intermolecular Raman response as a function of time after resonant excitation of the chromophore. Compared to electronically resonant probes of the solute, measuring the ultrafast dynamics using the nonresonant solvent response offers a new and complementary perspective on the events that accompany excitation and proton transfer. Experimentally observed changes in the nuclear polarizability of the solvent follow dynamic changes in the solvent-solute interactions. Reorganization of the solvent in response to the significant changes in the intermolecular interactions upon proton transfer is found to play an important role in the reaction dynamics. With transfer of the proton taking place rapidly, the solvent controls the dynamics via the time-dependent evolution of the free energy surface, even on subpicosecond time scales. In addition, the solvent response probes the effects of intermolecular energy transfer as energy released during the reactive event is rapidly transferred to the local solvent environment and then dissipates to the bulk solvent on about a 10 ps time scale. A brief initial account of a portion of this work has appeared previously, J. Am. Chem. Soc. 2004, 126, 8620-8621.     

Shielding of ionic interactions by sulfur dioxide in an ionic liquid

The effect of adding SO2 on the structure and dynamics of 1-butyl-3-methylimidazolium bromide (BMIBr) was investigated by low-frequency Raman spectroscopy and molecular dynamics (MD) simulations. The MD simulations indicate that the long-range structure of neat BMIBr is disrupted resulting in a liquid with relatively low viscosity and high conductivity, but strong correlation of ionic motion persists in the BMIBr-SO2 mixture due to ionic pairing. Raman spectra within the 5相似文献   

Correlation between quasielastic Raman scattering and configurational entropy in an ionic liquid

Low-frequency (5-200 cm(-1)) Raman spectra are reported for the ionic liquid 1-butyl-3-methylimidazolium hexafluorophosphate, [bmim]PF(6), in glassy, supercooled liquid, and normal liquid phases (77-330 K). Raman spectra of [bmim]PF(6) agree with previous results obtained by optical Kerr effect spectroscopy and molecular dynamics simulation. Both the superposition model and the coupling model give reasonable fit to low-frequency Raman spectra of [bmim]PF(6). The configurational entropy of [bmim]PF(6) has been evaluated as a function of temperature using recently reported data of heat capacity. The calculated configurational entropy is inserted in the Adam-Gibbs theory for supercooled liquids, giving a good fit to non-Arrhenius behavior of viscosity and diffusive process, with the latter revealed by a recent neutron scattering investigation of [bmim]PF(6). There is a remarkable linear dependence between intensity of quasielastic Raman scattering and configurational entropy from 77 K up to the melting point of [bmim]PF(6). This correlation offers insight into the nature of dynamical processes probed by low-frequency Raman spectra of ionic liquids.     

Intermolecular polarizability dynamics of aqueous formamide liquid mixtures studied by molecular dynamics simulations

A molecular dynamics simulation study is presented for the relaxation of the polarizability anisotropy in liquid mixtures of formamide and water, using a dipolar induction scheme that involves the intrinsic polarizability and first hyperpolarizability tensors of the molecules, and the dipole-quadrupole polarizability of water species. The long time diffusive decay of the collective polarizability anisotropy correlations exhibits a substantial slowing down as the formamide mole fraction increases in the mixture. The diffusive times for the polarizability relaxation obtained from the authors' simulations are in good agreement with optical Kerr effect experimental data, and they are found to correlate nearly linearly with the estimated mean lifetimes of the hydrogen bonds within the mixture, suggesting that the relaxation of the hydrogen bond network is responsible to some extent for the collective relaxation of the polarizability anisotropy of the mixture. The short time behavior of the polarizability anisotropy relaxation was investigated by computing the nuclear response function, R(t), which is very rapidly dominated by the formamide contribution as it is added to water, due to the much larger polarizability anisotropy of formamide molecules compared to that of water. Several contributions to the Raman spectrum were also analyzed as a function of composition, and the dynamical origin of the different bands was determined.     

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.     

Multidimensional infrared spectroscopy of water. I. Vibrational dynamics in two-dimensional IR line shapes

In this and the following paper, we describe the ultrafast structural fluctuations and rearrangements of the hydrogen bonding network of water using two-dimensional (2D) infrared spectroscopy. 2D IR spectra covering all the relevant time scales of molecular dynamics of the hydrogen bonding network of water were studied for the OH stretching absorption of HOD in D2O. Time-dependent evolution of the 2D IR line shape serves as a spectroscopic observable that tracks how different hydrogen bonding environments interconvert while changes in spectral intensity result from vibrational relaxation and molecular reorientation of the OH dipole. For waiting times up to the vibrational lifetime of 700 fs, changes in the 2D line shape reflect the spectral evolution of OH oscillators induced by hydrogen bond dynamics. These dynamics, characterized through a set of 2D line shape analysis metrics, show a rapid 60 fs decay, an underdamped oscillation on a 130 fs time scale induced by hydrogen bond stretching, and a long time decay constant of 1.4 ps. 2D surfaces for waiting times larger than 700 fs are dominated by the effects of vibrational relaxation and the thermalization of this excess energy by the solvent bath. Our modeling based on fluctuations with Gaussian statistics is able to reproduce the changes in dispersed pump-probe and 2D IR spectra induced by these relaxation processes, but misses the asymmetry resulting from frequency-dependent spectral diffusion. The dynamical origin of this asymmetry is discussed in the companion paper.     

Selective excitation of molecular modes in a mixture by optimal control of electronically nonresonant femtosecond four-wave mixing spectroscopy

Femtosecond time-resolved coherent anti-Stokes Raman scattering (fs-CARS) gives access to ultrafast molecular dynamics. Due to the spectrally broad laser pulses, usually poorly resolved spectra result from this spectroscopy. However, it can be demonstrated that by shaping the femtosecond pulses a selective excitation of specific vibrational modes is possible. We demonstrate that using a feedback-controlled optimization technique, molecule-specific CARS spectra can be obtained from a mixture of different substances. A careful analysis of the experimental results points to a nontrivial control of the vibrational mode dynamics in the electronic ground state of the molecules as underlying mechanism.     

Induced contributions in the rayleigh spectra of water: A molecular dynamics simulation

A molecular dynamics simulation of TIP4P water at T=280 K has been performed to calculate the depolarized and isotropic Raman spectra in the translational and librational region. Comparison of the isotropic spectra with experiment confirms the assignment of the low-frequency scattered light to collision-induced phenomena and allows identification of the spectral contributions arising from polarizability modulation due to short-range interactions.     

Probing structural evolution along multidimensional reaction coordinates with femtosecond stimulated Raman spectroscopy

Mapping out multidimensional potential energy surfaces has been a goal of physical chemistry for decades in the quest to both predict and control chemical reactivity. Recently a new spectroscopic approach called Femtosecond Stimulated Raman Spectroscopy or FSRS was introduced that can structurally interrogate multiple dimensions of a reactive potential energy surface. FSRS is an ultrafast laser technique which provides complete time-resolved, background-free Raman spectra in a few laser shots. The FSRS technique provides simultaneous ultrafast time (~50 fs) and spectral (~8 cm(-1)) resolution, thus enabling one to follow reactive structural evolutions as they occur. In this perspective we summarize how FSRS has been used to follow structural dynamics and provide mechanistic detail on three classical chemical reactions: a structural isomerization, an electron transfer reaction, and a proton transfer reaction.     

Structural dynamics of surfaces by ultrafast electron crystallography: experimental and multiple scattering theory

Recent studies in ultrafast electron crystallography (UEC) using a reflection diffraction geometry have enabled the investigation of a wide range of phenomena on the femtosecond and picosecond time scales. In all these studies, the analysis of the diffraction patterns and their temporal change after excitation was performed within the kinematical scattering theory. In this contribution, we address the question, to what extent dynamical scattering effects have to be included in order to obtain quantitative information about structural dynamics. We discuss different scattering regimes and provide diffraction maps that describe all essential features of scatterings and observables. The effects are quantified by dynamical scattering simulations and examined by direct comparison to the results of ultrafast electron diffraction experiments on an in situ prepared Ni(100) surface, for which structural dynamics can be well described by a two-temperature model. We also report calculations for graphite surfaces. The theoretical framework provided here allows for further UEC studies of surfaces especially at larger penetration depths and for those of heavy-atom materials.     

Comparative study of normal and branched alkane monolayer films adsorbed on a solid surface. II. Dynamics

The dynamics of monolayer films of the n-alkane tetracosane (n-C24H52) and the branched alkane squalane (C30H62) adsorbed on graphite have been studied by quasielastic and inelastic neutron scattering and molecular dynamics (MD) simulations. Both molecules have 24 carbon atoms along their carbon backbone, and squalane has an additional six methyl side groups symmetrically placed along its length. The authors' principal objective has been to determine the influence of the side groups on the dynamics of the squalane monolayer and thereby assess its potential as a nanoscale lubricant. To investigate the dynamics of these monolayers they used both the disk chopper spectrometer (DCS) and the high flux backscattering spectrometer (HFBS) at the National Institute of Standards and Technology. These instruments made it possible to study dynamical processes such as molecular diffusive motions and vibrations on very different time scales: 1-40 ps (DCS) and 0.1-4 ns (HFBS). The MD simulations were done on corresponding time scales and were used to interpret the neutron spectra. The authors found that the dynamics of the two monolayers are qualitatively similar on the respective time scales and that there are only small quantitative differences that can be understood in terms of the different masses and moments of inertia of the two molecules. In the course of this study, the authors developed a procedure to separate out the low-frequency vibrational modes in the spectra, thereby facilitating an analysis of the quasielastic scattering. They conclude that there are no major differences in the monolayer dynamics caused by intramolecular branching. It remains to be seen whether this similarity in monolayer dynamics also holds for the lubricating properties of these molecules in confined geometries.