Ab initio vibrational calculations for H2SO4 and H2SO4 x H2O: spectroscopy and the nature of the anharmonic couplings

Vibrational frequencies for fundamental, overtone, and combination excitations of sulfuric acid (H2SO4) and of sulfuric acid monohydrate cluster (H2SO4 x H2O) are computed directly from ab initio MP2/TZP potential surface points using the correlation-corrected vibrational self-consistent field (CC-VSCF) method, which includes anharmonic effects. The results are compared with experiment. The computed transitions show in nearly all cases good agreement with experimental data and consistent improvement over the harmonic approximation. The CC-VSCF improvements over the harmonic approximation are largest for the overtone and combination excitations and for the OH stretching fundamental. The agreement between the calculations and experiment also supports the validity of the MP2/TZP potential surfaces. Anharmonic coupling between different vibrational modes is found to significantly affect the vibrational frequencies. Analysis of the mean magnitude of the anharmonic coupling interactions between different pairs of normal modes is carried out. The results suggest possible mechanisms for the internal flow of vibrational energy in H2SO4 and H2SO4 x H2O.     

Vibrational spectroscopy of protonated imidazole and its complexes with water molecules: ab initio anharmonic calculations and experiments

The results of anharmonic frequency calculations on neutral imidazole (C3N2H4, Im), protonated imidazole (ImH+), and its complexes with water (ImH+)(H2O)n, are presented and compared to gas phase infrared photodissociation spectroscopy (IRPD) data. Anharmonic frequencies are obtained via ab initio vibrational self-consistent field (VSCF) calculations taking into account pairwise interactions between the normal modes. The key results are: (1) Prediction of anharmonic vibrational frequencies on an MP2 ab initio potential energy surface show excellent agreement with experiment and outstanding improvement over the harmonic frequencies. For example, the ab initio calculated anharmonic frequency for (ImH+)(H2O)N2 exhibits an overall average percentage error of 0.6% from experiment. (2) Anharmonic vibrational frequencies calculated on a semiempirical potential energy surface fitted to ab initio harmonic data represents spectroscopy well, particularly for water complexes. As an example, anharmonic frequencies for (ImH+)H2O and (ImH+)(H2O)2 show an overall average deviation of 1.02% and 1.05% from experiment, respectively. This agreement between theory and experiment also supports the validity and use of the pairwise approximation used in the calculations. (3) Anharmonic coupling due to hydration effects is found to significantly reduce the vibrational frequencies for the NH stretch modes. The frequency of the NH stretch is observed to increase with the removal of a water molecule or replacement of water with N2. This result also indicates the ability of the VSCF method to predict accurate frequencies in a matrix environment. The calculation provides insights into the nature of anharmonic effects in the potential surface. Analysis of percentage anharmoncity in neutral Im and ImH+ shows a higher percentage anharmonicity in the NH and CH stretch modes of neutral Im. Also, we observe that anharmonicity in the NH stretch modes of ImH+ have some contribution from coupling effects, while that of neutral Im has no contribution whatsoever from mode-mode coupling. It is concluded that the incorporation of anharmonic effects in the calculation brings theory and experiment into much closer agreement for these systems.     

DFT calculations, DRIFT spectroscopy and kinetic studies on the hydrolysis of isocyanic acid on the TiO2-anatase (1 0 1) surface

The co-adsorption of isocyanic acid (HNCO) and water (H₂O) and their reaction to ammonia and carbon dioxide on the anatase phase of TiO₂ were studied with ab initio density functional theory (DFT) calculations using a cluster model as well as with in situ DRIFTS investigations and kinetic experiments. We found that isocyanic acid can in principle adsorb both molecularly and dissociatively on the TiO₂(1 0 1) surface, but the moment at which water gets involved in the process, is vital for determining the further course of the surface reaction. In the absence of water, it was found that HNCO can adsorb in molecular form on the TiO₂ surface. Assuming this case to be the first step of the HNCO hydrolysis, the surface HNCO rearranges into an intermediate complex with a modified NCO skeleton. After decarboxylation water attacks the complex from the gas phase and ammonia is finally formed.

However, when water is present at the beginning of the hydrolysis reaction, it immediately attacks the NCO group present at the surface, yielding a carbamic acid complex, which is further transformed into a carbamate complex. After decarboxylation an NH₂ group remains at the surface. Finally, NH₃ is formed by hydrogen transfer from molecularly adsorbed water at a neighboring titanium center and the hydrolysis reaction is completed.

Since water is always present in diesel exhaust gas, only the second mechanism is relevant under practical conditions. Moreover, the calculated energy barrier is lower for the second mechanism compared to the first reaction pathway. The comparison between the sum of the theoretical vibrational spectra of the reaction intermediates with the in situ DRIFT spectra also strongly supports the accuracy of the second reaction pathway. The experimental investigation of the kinetics of the HNCO hydrolysis on TiO₂-anatase revealed a second order reaction—first order with respect to HNCO and first order with respect to water, which can only be reconciled with the second mechanism.     

Vibrational spectroscopy of (SO4(2-)).(H2O)n clusters, n=1-5: harmonic and anharmonic calculations and experiment

The vibrational spectroscopy of (SO4(2-)).(H2O)n is studied by theoretical calculations for n=1-5, and the results are compared with experiments for n=3-5. The calculations use both ab initio MP2 and DFT/B3LYP potential energy surfaces. Both harmonic and anharmonic calculations are reported, the latter with the CC-VSCF method. The main findings are the following: (1) With one exception (H2O bending mode), the anharmonicity of the observed transitions, all in the experimental window of 540-1850 cm(-1), is negligible. The computed anharmonic coupling suggests that intramolecular vibrational redistribution does not play any role for the observed linewidths. (2) Comparison with experiment at the harmonic level of computed fundamental frequencies indicates that MP2 is significantly more accurate than DFT/B3LYP for these systems. (3) Strong anharmonic effects are, however, calculated for numerous transitions of these systems, which are outside the present observation window. These include fundamentals as well as combination modes. (4) Combination modes for the n=1 and n=2 clusters are computed. Several relatively strong combination transitions are predicted. These show strong anharmonic effects. (5) An interesting effect of the zero point energy (ZPE) on structure is found for (SO4(2-)).(H2O)(5): The global minimum of the potential energy corresponds to a C(s) structure, but with incorporation of ZPE the lowest energy structure is C2v, in accordance with experiment. (6) No stable structures were found for (OH-).(HSO4-).(H2O)n, for n相似文献   

Molecular orientation of CO adsorbed on clean and Al-modified Ru(0 0 0 1) surfaces

Chemisorption of CO on the bare and Al precovered Ru(0 0 0 1) surface as well as on the Al₂Ru surface alloy, was studied by means of angle-resolved photoelectron spectroscopy using unpolarised He II radiation. Whereas CO molecules stand upright on the clean Ru(0 0 0 1) surface, as deduced by applying the symmetry selection rules of photoemission, they adsorb in a tilted geometry on the Al-precovered Ru(0 0 0 1) surface. On an Al₂Ru surface alloy, produced by annealing the Al-covered Ru(0 0 0 1) surface up to 1100 K, the CO molecules are again found to stand upright, as on the clean Ru(0 0 0 1) surface. The Al atoms, incorporated into the topmost Ru layer, do not have a significant influence on the orientation of the axis of the CO molecules.     

Ab initio calculations of anharmonic vibrational spectroscopy for hydrogen fluoride (HF)n (n = 3, 4) and mixed hydrogen fluoride/water (HF)n(H2O)n (n = 1, 2, 4) clusters

Anharmonic vibrational frequencies and intensities are computed for hydrogen fluoride clusters (HF)n, with n = 3, 4 and mixed clusters of hydrogen fluoride with water (HF)n(H2O)n where n = 1, 2. For the (HF)4(H2O)4 complex, the vibrational spectra are calculated at the harmonic level, and anharmonic effects are estimated. Potential energy surfaces for these systems are obtained at the MP2/TZP level of electronic structure theory. Vibrational states are calculated from the potential surface points using the correlation-corrected vibrational self-consistent field method. The method accounts for the anharmonicities and couplings between all vibrational modes and provides fairly accurate anharmonic vibrational spectra that can be directly compared with experimental results without a need for empirical scaling. For (HF)n, good agreement is found with experimental data. This agreement shows that the M?ller-Plesset (MP2) potential surfaces for these systems are reasonably reliable. The accuracy is best for the stiff intramolecular modes, which indicates the validity of MP2 in describing coupling between intramolecular and intermolecular degrees of freedom. For (HF)n(H2O)n experimental results are unavailable. The computed intramolecular frequencies show a strong dependence on cluster size. Intensity features are predicted for future experiments.     

Addition reactions of nitrones on the reconstructed C(1 0 0)-2 × 1 surfaces

In this paper, the reactions of nitrone, N-methyl nitrone, N-phenyl nitrone and their hydroxylamine tautomers (vinyl-hydroxylamine, N-methyl-vinyl-hydroxylamine and N-phenyl-vinyl-hydroxylamine) on the reconstructed C(1 0 0)-2 × 1 surface have been investigated using hybrid density functional theory (B3LYP), Møller–Plesset second-order perturbation (MP2) and multi-configuration complete-active-space self-consistent-field (CASSCF) methods. The calculations showed that all the nitrones can react with the surface “dimer” via facile 1,3-dipolar cycloaddition with small activation barriers (less than 12.0 kJ/mol at B3LYP/6-31g(d) level). The [2+2] cycloaddition of hydroxylamine tautomers on the C(1 0 0) surface follows a diradical mechanism. Hydroxylamine tautomers first form diradical intermediates with the reconstructed C(1 0 0)-2 × 1 surface by overcoming a large activation barrier of 50–60 kJ/mol (B3LYP), then generate [2+2] cycloaddition products via diradical transition states with negligible activation barriers. The surface reactions result in hydroxyl or amino-terminated diamond surfaces, which offers new opportunity for further modifications.     

The emergence of collective vibrations in cluster models: quantum chemical study of the methyl-terminated Si(111) surface

In this paper we present structures and harmonic vibrational frequencies for the methylated silicon (111) surface from quantum chemical calculations using both cluster models and periodic boundary conditions. The results from both calculations are in very good agreement with experimentally determined frequencies. We demonstrate that relatively small cluster models already show the emergence of collective vibrational modes and provide a general method for the assignment of vibrational frequencies for extended surfaces from cluster models. Finally, we discuss a vibrational mode that results from the coupling between near-surface phonons and the silicon-carbon bending modes.     

Prediction of accurate anharmonic experimental vibrational frequencies for water clusters, (H2O)n, n=2-5

Accurate anharmonic experimental vibrational frequencies for water clusters consisting of 2-5 water molecules have been predicted on the basis of comparing different methods with MP2/aug-cc-pVTZ calculated and experimental anharmonic frequencies. The combination of using HF/6-31G* scaled frequencies for intramolecular modes and anharmonic frequencies for intermolecular modes gives excellent agreement with experiment for the water dimer and trimer and are as good as the expensive anharmonic MP2 calculations. The water trimer, the cyclic Ci and S4 tetramers, and the cyclic pentamer all have unique peaks in the infrared spectrum between 500 and 800 cm-1 and between 3400 and 3700 cm-1. Under the right experimental conditions these different clusters can be uniquely identified using high-resolution IR spectroscopy.     

Collective vibrations in cluster models for semiconductor surfaces: vibrational spectra of acetylenyl and methylacetylenyl functionalized Si(111)

The geometries and harmonic vibrational frequencies of the acetylenyl and methylacetylenyl functionalized Si(111) surfaces are investigated using quantum chemical calculations. The vibrational spectra are computed using a previously introduced method whereby the collective vibrational modes that correspond to the vibrations of the infinite periodic system are derived from modest sized cluster models. Our predictions should be useful for the interpretation of the experimental spectra when they become available. The symmetry elements of the methylacetylenyl Si(111) surface that are derived from the space group of the optimized structure and a vibrational mode resulting from photon-adsorbate coupling are explored.     

Anharmonic overtone and combination states of glycine and two model peptides examined by vibrational self-consistent field theory

In this paper, the application of the vibrational self-consistent field (VSCF) and correction-corrected VSCF methods for calculating anharmonic parameters, including transition frequency, transition intensity and dipole, and vibrational anharmonicity of 3N-6 normal modes for formamide, glycine, N-methylacetamide and their deuterated derivatives are explored mainly at the level of density functional theory. The computed fundamental anharmonic frequencies are found to be in reasonable agreement with experimental results. Diagonal anharmonicities of the second overtone states were examined for multiple normal modes, whose magnitudes were found to correlate well with those of the first overtone states in the three small molecules. The results show that the VSCF-based approach can be utilized to predict anharmonic parameters of higher vibrational states that are essential to understanding multi-pulse infrared nonlinear experiments of peptides.     

2D calculation of anharmonic OH vibrations in a layered hydroxide crystal

Anharmonic vibrational frequencies for the Raman-active (A(1g)) and the IR-active (A(2u)) modes have been calculated for the LiOH crystal within a plane-wave density functional theory (DFT) framework. We find that a two-dimensional quantum-mechanical vibrational approach, allowing for anharmonic coupling between symmetric and antisymmetric OH stretching modes, produces OH frequencies--both absolute frequencies and gas-to-solid frequency shifts--in good agreement with experiment. Remaining errors in the absolute frequencies are largely a consequence of the DFT model chosen. A one-dimensional normal-mode following vibrational treatment, on the other hand, fails to reproduce both absolute anharmonic frequencies and gas-to-solid frequency shifts.     

Accurate calculation of anharmonic vibrational frequencies of medium sized molecules using local coupled cluster methods

Local coupled cluster methods were applied for the automated generation of accurate multidimensional potential energy surfaces for a set of test molecules ranging from six to nine atoms. Based on these surfaces anharmonic fundamental frequencies were computed using vibrational self-consistent field and configuration interaction methods. The computed vibrational frequencies are compared to those obtained from similar calculations using conventional coupled cluster methods and to experimental values. The results from local and conventional methods are found to be of similar accuracy and in close agreement with experimental values. In addition, an efficient parallelization of the fully automated surface generation code is presented.     

Anharmonic vibrational calculations modeling the raman spectra of intermediates in the photoactive yellow protein (PYP) photocycle

The role of anharmonic effects in the vibrational spectroscopy of the dark state and two major chromophore intermediates of the photoactive yellow protein (PYP) photocycle is examined via ab initio vibrational self-consistent field (VSCF) calculations and time-resolved resonance Raman spectroscopy. For the first time, anharmonicity is considered explicitly in calculating the vibrational spectra of an ensemble consisting of the PYP chromophore surrounded by model compounds used as mimics of the important active-site residues. Predictions of vibrational frequencies on an ab initio corrected semiempirical potential energy surface show remarkable agreement with experimental frequencies for all three states, thus shedding light on the potential along the reaction path. For example, calculated frequencies for vibrational modes of the red-shifted intermediate, PYPL, exhibit an overall average error of 0.82% from experiment. Upon analysis of anharmonicity patterns in the PYP modes we observe a decrease in anharmonicity in the C8-C9 stretching mode nu29 (trans-cis isomerization marker mode) with the onset of the cis configuration in PYPL. This can be attributed to the loss of the hydrogen-bonding character of the adjacent C9-O2 to the methylamine (Cys69 backbone). For several of the modes, the anharmonicity is mostly due to mode-mode coupling, while for others it is mostly intrinsic. This study shows the importance of the inclusion of anharmonicity in theoretical spectroscopic calculations, and the sensitivity of experiments to anharmonicity. The characterization of protein active-site residues by small molecular mimics provides an acceptable chemical structural representation for biomolecular spectroscopy calculations.     

Fundamental vibrational frequencies and dominant resonances in methylamine isotopologues by ab initio and density functional theory methods

Ab initio and density functional theory (DFT) calculations were performed for obtaining fundamental vibrational frequencies of methylamine, CH3NH2, and its deuterated variants CH3ND2, CD3NH2, and CD3ND2. The calculations were carried out using the CCSD(T) coupled cluster approximation with cc-pVTZ and cc-pVQZ basis sets, and by the DFT method with the semiempirical hybrid functional B97-1 with polarization consistent pc-2 and pc-3 basis sets. Reasonable performance of the DFT harmonic and ab initio harmonic calculations was found, which improved considerably upon combination of the harmonic fundamental frequencies with anharmonic corrections from the smaller, pc-2, basis. The computed anharmonic fundamental frequencies of methylamine isotopologues agree very well with the experimental values and represent a useful tool for assignment and analysis of the dominant resonances.     

Spectrum and dynamics of the CH chromophore in CD2HF. II. Ab initio calculations of the potential and dipole moment functions

We report ab initio (SCF-MP2 and CI) calculations on the three-dimensional anharmonic potential and dipole functions of the coupled vibrational modes of the CH chromopore in CD₂HF. Predictions of fundamental and overtone spectra are obtained from 3D solutions of the vibrational Schrödinger equation and are compared with experiment. The dominant Fermi and Darling-Dennison coupling constants in the effective Hamiltonian representation of experiment and theory agree well. Predicted overtone band strengths are satisfactory and very sensitive to the dipole function used.     

Intramolecular vibrational coupling contribution to temperature dependence of vibrational mode frequencies

High-frequency vibrational modes in molecules in solution are sensitive to temperature and shift either to lower or higher frequencies with the temperature increase. These frequency shifts are often attributed to specific interactions of the molecule and to the solvent polarization effect. We found that a substantial and often dominant contribution to sensitivity of vibrational high-frequency modes to temperature originates from anharmonic interactions with other modes in the molecule. The temperature dependencies were measured for several modes in ortho-, meta-, and para-isomers of acetylbenzonitrile in solution and in a solid matrix and compared to the theoretical predictions originated from the intramolecular vibrational coupling (IVC) evaluated using anharmonic density functional theory calculations. It is found that the IVC contribution is essential for temperature dependencies of all high-frequency vibrational modes and is dominant for many modes. As such, the IVC contribution alone permits predicting the main trend in the temperature dependencies, especially for vibrational modes with smaller transition dipoles. In addition, an Onsager reaction field theory was used to describe the solvent contribution to the temperature dependencies.     

Rapid anharmonic vibrational corrections derived from partial Hessian analysis

Vibrational analysis within a partial Hessian framework can successfully describe the vibrational properties of a variety of systems where the vibrational modes of interest are localized within a specific region of the system. We have developed a new approach to calculating anharmonic frequencies based on vibrational frequencies and normal modes obtained from a partial Hessian analysis using second-order vibrational perturbation theory and the transition optimized shifted Hermite method. This allows anharmonic frequencies for vibrational modes that are spatially localized to be determined at a significantly reduced computational cost. Several molecular systems are examined in order to demonstrate the effectiveness of this method including organic molecules adsorbed on the Si(100)-2×1 surface, model peptides in solution, and the C-H stretching region of polycyclic aromatic hydrocarbons. Overall, for a range of systems, anharmonic frequencies calculated using the partial Hessian approach are found to be in close agreement with the results obtained using full anharmonic calculations while providing a significant reduction in computational cost.     

Anharmonic vibrational spectroscopy and investigation of intramolecular mode couplings in adenine

Vibrational frequencies for the nucleobase adenine are calculated by the vibrational self-consistent field (VSCF) and correlation corrected vibrational self-consistent field (CC-VSCF) methods using Hartree-Fock (HF), density functional theory (DFT) and second order Møller-Plesset (MP2) theories. A large number of potential energy surface (PES) points were computed in the anharmonic calculations corresponding to each method. The quartic force field (QFF) approximation was used to generate the full grid of points for the VSCF solver. We have implemented our new procedure for computing the mode-mode coupling integrals in the 2-mode coupling representations of the quartic force field (2MR-QFF) for prediction of coupling magnitudes. Calculations were performed using the 6-31G(d,p) basis set. Comparison of the calculated ab initio anharmonic spectra with Ar matrix experimental data of adenine reported in the literature reveals that, the CC-VSCF (DFT) wavenumbers show the best agreement. The experimental geometric parameters of adenine are compared with the theoretically optimized molecular structural parameters. These are found to be in good agreement. Vibrational assignments are based on the calculated potential energy distribution (PED) values.     

Anharmonic Vibrational Motions In C<Subscript>60</Subscript>: A Potential Energy Surface Derived From Vibrational Self-Consistent Field Calculations

A potential energy surface (PES) is developed for C₆₀ designed to describe vibrational motions valid in the anharmonic limit. The PES is based on a previously existing one that is fit to the harmonic fundamentals and is then modified to generate anharmonicity of all orders and in all terms, but without additional fitted parameters. The resulting Cartesian vibrational motions are decomposed into normal modes, and the anharmonic expansion coefficients are calculated including 2-mode couplings and up to 4th order. The resulting PES is used in a vibrational self-consistent field (VSCF) algorithm to calculate the anharmonically corrected fundamental frequencies. The parameters are then fit to fundamental infrared and Raman frequencies. While it is not possible to assign combination and overtone transitions with sufficient experimental accuracy, conclusions about the effects of anharmonic vibrational coupling in C₆₀ are described.