Vibrational self-consistent field approach to anharmonic spectroscopy of molecules in solids: Application to iodine in argon matrix

On the synergy of matrix-isolation infrared spectroscopy and vibrational configuration interaction computations

The key feature of matrix-isolation infrared (MI-IR) spectroscopy is the isolation of single guest molecules in a host system at cryogenic conditions. The matrix mostly hinders rotation of the guest molecule, providing access to pure vibrational features. Vibrational self-consistent field (VSCF) and configuration interaction computations (VCI) on ab initio multimode potential energy surfaces (PES) give rise to anharmonic vibrational spectra. In a single-sourced combination of these experimental and computational approaches, we have established an iterative spectroscopic characterization procedure. The present article reviews the scope of this procedure by highlighting the strengths and limitations based on the examples of water, carbon dioxide, methane, methanol, and fluoroethane. An assessment of setups for the construction of the multimode PES on the example of methanol demonstrates that CCSD(T)-F12 level of theory is preferable to compute (a) accurate vibrational frequencies and (b) equilibrium or vibrationally averaged structural parameters. Our procedure has allowed us to uniquely assign unknown or disputed bands and enabled us to clarify problematic spectral regions that are crowded with combination bands and overtones. Besides spectroscopic assignment, the excellent agreement between theory and experiment paves the way to tackle questions of rather fundamental nature as to whether or not matrix effects are systematic, and it shows the limits of conventional notations used by spectroscopists.

Decomposing anharmonicity and mode-coupling from matrix effects in the IR spectra of matrix-isolated carbon dioxide and methane

A combined experimental and computational approach revealed similarities and differences in the vibrational signature of matrix-isolated carbon dioxide and methane.

Understanding of matrix embedding: a theoretical spectroscopic study of CO interacting with Ar clusters, surfaces and matrices

Through benchmark studies, we explore the performance of PBE density functional theory, with and without Grimme's dispersion correction (DFT-D3), in predicting spectroscopic properties for molecules interacting with rare gas matrices.

Hybrid MP2/MP4 potential surfaces in VSCF calculations of IR spectra: applications for organic molecules

This study introduces an improved hybrid MP2/MP4 ab initio potential for vibrational spectroscopy calculations which is very accurate, yet without high computational demands. The method uses harmonic vibrational calculations with the MP4(SDQ) potential to construct an improved MP2 potential by coordinate scaling. This improved MP2 potential is used for the anharmonic VSCF calculation. The method was tested spectroscopically for four molecules: butane, acetone, ethylene and glycine. Very good agreement with experiment was found. For most of the systems, the more accurate harmonic treatment considerably improved the MP2 anharmonic results.

Raman spectra of long chain hydrocarbons: anharmonic calculations, experiment and implications for imaging of biomembranes

First-principles anharmonic vibrational calculations are carried out for the Raman spectrum of the C-H stretching bands in dodecane, and for the C-D bands in the deuterated molecule. The calculations use the Vibrational Self-Consistent Field (VSCF) algorithm. The results are compared with liquid-state experiments, after smoothing the isolated-molecule sharp-line computed spectra. Very good agreement between the computed and experimental results is found for the two systems. The combined theoretical and experimental results provide insights into the spectrum, elucidating the roles of symmetric and asymmetric CH(3) and CH(2) hydrogenic stretches. This is expected to be very useful for the interpretation of spectra of long-chain hydrocarbons. The results show that anharmonic effects on the spectrum are large. On the other hand, vibrational degeneracy effects seem to be rather modest at the resolution of the experiments. The degeneracy effects may have more pronounced manifestations in higher-resolution experiments. The results show that first-principles anharmonic vibrational calculations for hydrocarbons are feasible, in good agreement with experiment, opening the way for applications to many similar systems. The results may be useful for the analysis of CARS imaging of lipids, for which dodecane is a representative molecule. It is suggested that first-principles vibrational calculations may be useful also for CARS imaging of other systems.

Applications of terahertz spectroscopy in biosystems

Terahertz (THz) spectroscopic investigations of condensed-phase biological samples are reviewed ranging from the simple crystalline forms of amino acids, carbohydrates and polypeptides to the more complex aqueous forms of small proteins, DNA and RNA. Vibrationally resolved studies of crystalline samples have revealed the exquisite sensitivity of THz modes to crystalline order, temperature, conformational form, peptide sequence and local solvate environment and have given unprecedented measures of the binding force constants and anharmonic character of the force fields, properties necessary to improve predictability but not readily obtainable using any other method. These studies have provided benchmark vibrational data on extended periodic structures for direct comparisons with classical (CHARMm) and quantum chemical (density functional theory) theories. For the larger amorphous and/or aqueous phase samples, the THz modes form a continuum-like absorption that arises because of the full accessibility to conformational space and/or the rapid time scale for inter-conversion in these environments. Despite severe absorption by liquid water, detailed investigations have uncovered the photo- and hydration-induced conformational flexibility of proteins, the solvent shell depth of the water/biomolecule boundary layers and the solvent reorientation dynamics occurring in these interfacial layers that occur on sub-picosecond time scales. As such, THz spectroscopy has enhanced and extended the accessibility to intermolecular forces, length- and timescales important in biological structure and activity.

Molecular dynamics simulations of matrix deposition. III. Site structure analysis for porphycene in argon and xenon

Porphycene (1) and porphyrin (2), two constitutional isomers, reveal completely different electronic spectral patterns in argon and xenon matrices. For the former the spectra recorded in the two solidified gases resemble each other, whereas for the latter they are completely different. This difference can be rationalized by molecular-dynamics simulations of the structure of the microenvironment carried out for the two chromophores embedded in argon and xenon hosts. For 1, the structure of the main substitutional site is the same for Ar and Xe and consists of a hexagonal cavity obtained by removing seven host atoms from the [111] crystallographic plane. An analogous structure is obtained for 2 in xenon. However, in argon the porphyrin chromophore environment is shared between several different sites, with the number of replaced host atoms ranging from seven to ten. These results demonstrate that a relatively minor structural alternation may lead to major changes in the spectral pattern of molecules embedded in rare-gas cryogenic matrices.

Femtosecond coherent anti-Stokes Raman-scattering polarization beat spectroscopy of I2–Xe complex in solid krypton

Time-resolved coherent anti-Stokes Raman-scattering (CARS) measurements are carried out to study the interaction between xenon atom and iodine molecule in a solid krypton matrix. Interference between the CARS polarizations of the "free" and complexed iodine molecules is observed, while the quantum beats of the complex are not detected due to low concentration. Vibrational analysis based on the polarization beats yields accurate molecular constants for the I2-Xe complex. The harmonic frequency of the I2-Xe complex is found to be redshifted by 0.90 cm-1 when compared to the free I2, whereas the anharmonicity is approximately the same. The dephasing rate of the complex is found to be somewhat higher than that of the free iodine molecule in solid Kr, showing that the complexation affects dephasing, although not dramatically. Molecular dynamics simulations are carried out to find the conformation of the complex, and wave packet simulations are used to reproduce the CARS signal to confirm the assignments of the observed beatings as quantum and polarization beats. The results show that the polarization beats are a useful tool for investigating weak interactions in condensed phase.

Formation of novel rare-gas molecules in low-temperature matrices

Progress in the study of a new class of chemically bound compounds of noble-gas atoms is reviewed. The focus is on rare-gas molecules of the form HNgY, where Ng is a noble-gas atom and Y is an electronegative group, prepared by photolysis of HY in the rare-gas matrix. Other related types of new molecules of noble-gas atoms are discussed as well. Topics discussed in this review include: (a) The nature of bonding and the energetic stability of the compounds. (b) The vibrational spectroscopy of the molecules, and its role in identification of the species. (c) The mechanism and dynamics of photochemical formation of HNgY in the matrix, and the pathways for thermal and infrared (IR)-induced decomposition. Specifically, attention is given to the issue of "direct" formation following photolysis of HY versus "delayed" formation involving H atom diffusion. (d) Molecules of the lighter rare gases Ar, Ne, and He, focusing on the experimentally prepared HArF and on theoretical predictions suggesting the existence of other molecules. (e) The most-recently discovered photochemically induced insertion compounds of Ng into hydrocarbons, such as HXeCCH. (f) Clusters of HNgY with other molecules. The possible existence of neat aggregates and crystals of HNgY. The reviewed state-of-the-art suggests this field is at an early stage of development with major open questions bearing on the surprising properties of the molecules and on the formation mechanisms. These are part of the challenge for the future.

Dynamics and the breaking of a driven cage: I2 in solid Ar

Pump-probe measurements of I2 in solid Ar are reported and analyzed to extract a description of cage response to impulsive excitation, from the gentle kick, up to the breaking point. The most informative data are obtained through wavepacket motion on cage-bound, but otherwise dissociative, potentials where the chromophore acts as a transducer to drive the cage and to report on the local dynamics. This general class of dynamics is identified and analyzed as a function of energy in Ar, Kr, and Xe. The overdriven cage rebounds with a characteristic period of 1.2 ps that shows little dependence on excitation amplitude in all hosts. After rebound, the cage rings as a local resonant mode in Ar, with a period of 1 ps and dephasing time of 3 ps. This mode remains at the Debye edge in Kr and Xe, with periods of 630 and 800 fs, and dephasing times of 8 and 6 ps, respectively. In the bound B-state, the cage fluctuates toward its dilated equilibrium structure on a time scale of 3 ps, which is extracted from the down-chirp in the molecular vibrational frequency. When kicked with excess energy of 4 eV, the Ar cage breaks with 50% probability, and the molecule dissociates. The kinetics of polarization selective, multiphoton dissociation with Gaussian laser intensity profiles is delineated and the ballistics of cage breakout is described: The photodissociation proceeds by destruction of the local lattice, by creating interstitials and vacancies. During large amplitude motion on cage-bound potentials, sudden, nonadiabatic spin-flip transitions can be observed and quantified in space and time. The spin-flip occurs with unit probability in Ar when the I*-I bond is stretched beyond 6 A.

Vibrational coherence of I2 in solid Kr

Time-resolved coherent anti-Stokes Raman scattering, with a resolution of 20 fs, is used to prepare a broadband vibrational superposition on the ground electronic state of I2 isolated in solid Kr. The coherent evolution of a packet consisting of nu=1-6 is monitored for as many as 1000 periods, allowing a precise analysis of the material response and radiation coherence. The molecular vibrations are characterized by omega(e)=211.330(2) cm(-1), omega(e)x(e)=0.6523(6) cm(-1), omega(e)y(e)=2.9(1) x 10(-3) cm(-1); the dephasing rates at 32 K range from 110 ps for nu=1 to 34 ps for nu=6, with nu dependence: gamma(nu)=8.5 x 10(-3)+4.9 x 10(-4)nu2+2.1 x 10(-6)nu4 ps(-1). The signal amplitude is also modulated at omega(q)=41.56(3) cm(-1); which can be interpreted as coupling between the molecule and a local mode. The surprising implication is that this resonant local mode is decoupled from the lattice phonons, a finding that cannot be rationalized based on a normal-mode analysis.

Nonperturbative vibrational energy relaxation effects on vibrational line shapes

A general formulation of nonperturbative quantum dynamics of solutes in a condensed phase is proposed to calculate linear and nonlinear vibrational line shapes. In the weak solute-solvent interaction limit, the temporal absorption profile can be approximately factorized into the population relaxation profile from the off-diagonal coupling and the pure-dephasing profile from the diagonal coupling. The strength of dissipation and the anharmonicity-induced dephasing rate are derived in Appendix A. The vibrational energy relaxation (VER) rate is negligible for slow solvent fluctuations, yet it does not justify the Markovian treatment of off-diagonal contributions to vibrational line shapes. Non-Markovian VER effects are manifested as asymmetric envelops in the temporal absorption profile, or equivalently as side bands in the frequency domain absorption spectrum. The side bands are solvent-induced multiple-photon effects which are absent in the Markovian VER treatment. Exact path integral calculations yield non-Lorentzian central peaks in absorption spectrum resulting from couplings between population relaxations of different vibrational states. These predictions cannot be reproduced by the perturbative or the Markovian approximations. For anharmonic potentials, the absorption spectrum shows asymmetric central peaks and the asymmetry increases with anharmonicity. At large anharmonicities, all the approximation schemes break down and a full nonperturbative path integral calculation that explicitly accounts for the exact VER effects is needed. A numerical analysis of the O-H stretch of HOD in D(2)O solvent reveals that the non-Markovian VER effects generate a small recurrence of the echo peak shift around 200 fs, which cannot be reproduced with a Markovian VER rate. In general, the nonperturbative and non-Markovian VER contributions have a stronger effect on nonlinear vibrational line shapes than on linear absorption.