Understanding hydrogen scrambling and infrared spectrum of bare CH5+ based on ab initio simulations

A comprehensive study of the properties of protonated methane obtained from ab initio molecular dynamics simulations is presented. Comparing computed infrared spectra to the measured one gives further support to the high fluxionality of bare CH(5)(+). The computational trick to partially freezing out large-amplitude motion, in particular hydrogen scrambling and internal rotation of the H(2) moiety, leads to an understanding of the measured IR spectrum despite the underlying rapid hydrogen scrambling motion that interconverts dynamically structures of different symmetry and chemical bonding pattern. In particular, the fact that C-H stretching modes involving the carbon nucleus and those protons that form the H(2) moiety and the CH(3) tripod, respectively, result in distinct peaks is arguably experimental support for three-center two-electron bonding being operative at experimental conditions. It is proposed that hydrogen scrambling is associated with the softening of a mode that involves the bending of the H(2) moiety relative to the CH(3) tripod, which characterizes the C(s) ground-state structure. The potential energy surface that is mapped on to a two dimensional subspace of internal coordinates provides insight into the dynamical mechanism for exchange of hydrogens between CH(3) tripod and the three-center bonded H(2) moiety that eventually leads to full hydrogen scrambling.     

Analysis of the HO2 vibrational spectrum on an accurate ab initio potential energy surface

The complete vibrational spectrum of the HO2(X(2)A' ') radical, up to the H + O2 dissociation limit, has been determined quantum mechanically on an accurate potential energy surface (PES), based on approximately 15000 ab initio points at the icMRCI+Q/aug-cc-pVQZ level of theory. The vibrational states are found to be assignable at low energies but become more irregular as the energy approaches the dissociation limit. However, even at very high energies, regularity still exists, in sharp contrast to earlier results based on the double many-body expansion (DMBE) IV potential. Several Fermi resonances have been identified, and the spectrum is fit with a spectroscopic Hamiltonian. In addition, the vibrational dynamics is analyzed using a periodic orbit approach.     

The calculated infrared spectrum of Cl- H2O using a full dimensional ab initio potential surface and dipole moment surface

We report a full dimensional ab initio-based global potential energy surface (PES) and dipole moment surface (DMS) for Cl(-)H(2)O. Both surfaces are symmetric with respect to interchange of the H atoms. The PES is a fit to thousands of electronic energies calculated using the coupled-cluster method (CCSD(T)) with a moderately large basis (aug-cc-pVTZ). The infrared spectrum and vibrational dynamics are reported and compared to experiment. These results are analyzed by examination of wave function and the dipole surface.     

A three-dimensional ab initio potential energy surface and predicted infrared spectra for the He-N2O complex

We report a three-dimensional ab initio potential energy surface for He-N(2)O using a supermolecular method at the coupled-cluster singles and doubles with noniterative inclusion of connected triple level. Besides the intermolecular stretching and bending modes, we included the Q(3) normal mode for the nu(3) antisymmetric stretching vibration of N(2)O molecule in order to simulate the observed infrared spectra in the nu(3) region of N(2)O, especially to explain the frequency shift of the band origin in the infrared spectra. The harmonic oscillator approximation is used for the potential curve of the Q(3) mode of the isolate N(2)O molecule. The intermolecular potential energy surfaces are calculated for five potential-optimized discrete variable representation grid points of the Q(3) mode. The three-dimensional discrete variable representation method was employed to calculate the rovibrational states without separating the inter- and intramolecular nuclear motions. The calculated transition frequencies and line intensities of the rotational transitions in the nu(3) region of N(2)O for the van der Waals ground vibrational state are in good agreement with the observed infrared spectra. The calculated band shifts are found to be 0.1704 and 0.1551 cm(-1) for (4)He-N(2)O and (3)He-N(2)O, respectively, which agree well with the observed values of 0.2532 and 0.2170 cm(-1).     

The calculated infrared spectrum of Cl-H2O using a new full dimensional ab initio potential surface and dipole moment surface

We report a full dimensional, ab initio-based global potential energy surface (PES) and dipole moment surface for Cl-H2O. Both surfaces are symmetric with respect to interchange of the H atoms. The PES is a fit to thousands of electronic energies calculated using the coupled-cluster method [CCSD(T)] with a moderately large basis (aug-cc-pVTZ). Vibrational energies and wave functions are accurately obtained using MULTIMODE. The wave function and dipole moment surface are used to calculate and analyze the pure infrared spectrum at 0 K which is compared with experiment. Vibrational energies and the infrared spectra for DOD and HOD/DOH are also presented.     

Combined electronic structure/molecular dynamics approach for ultrafast infrared spectroscopy of dilute HOD in liquid H2O and D2O

We present a new approach that combines electronic structure methods and molecular dynamics simulations to investigate the infrared spectroscopy of condensed phase systems. This approach is applied to the OH stretch band of dilute HOD in liquid D2O and the OD stretch band of dilute HOD in liquid H2O for two commonly employed models of water, TIP4P and SPC/E. Ab initio OH and OD anharmonic transition frequencies are calculated for 100 HOD x (D2O)n and HOD x(H2O)n (n = 4-9) clusters randomly selected from liquid water simulations. A linear empirical relationship between the ab initio frequencies and the component of the electric field from the solvent along the bond of interest is developed. This relationship is used in a molecular dynamics simulation to compute frequency fluctuation time-correlation functions and infrared absorption line shapes. The normalized frequency fluctuation time-correlation functions are in good agreement with the results of previous theoretical approaches. Their long-time decay times are 0.5 ps for the TIP4P model and 0.9 ps for the SPC/E model, both of which appear to be somewhat too fast compared to recent experiments. The calculated line shapes are in good agreement with experiment, and improve upon the results of previous theoretical approaches. The methods presented are simple, and transferable to more complicated systems.     

An ab initio theoretical study of the optical spectrum of CF2

The ab initio calculation methods have been used to calculate the spectral and electronic characteristics of difluorocarbene in the ground electronic state (¹A₁), the lowest-lying singlet (¹B₁) and triplet (³B₁) states. The optimized equilibrium geometries, rotational constants, harmonic vibrational frequencies and energy gaps, electronic charges, dipole moments of these states have been computed with different basis sets. The calculated vibrational frequency of ³B₁ state (v₂=522 cm^?1) and the energy separation (2.26 eV) between ³B₁ and ¹A₁ states are in good agreement with the experimental results (519 cm^?1, 2.46 eV respectively). According to the calculations the previous assignment of vibrational symmetries of ¹B₁ state was corrected, and some experimentally undetermined vibrational frequencies were predicted.     

An ab initio study of the electronic spectrum of dichlorocarbene CCl2

The lower-lying electronic states of dichlorocarbene have been studied by ab initio methods at different levels of accuracy. For the ¹B₁ ← ¹A₁ transition, the calculated transition energy (1.95 eV) is in good agreement with the gas-phase value (2.10 eV) of Predmore, Murray and Harmony. The discrepancy with the previous CI study by Ha. Gremlich and Bühler is commented on. The vibrational frequencies of CCl₂ in both ¹A₁ and ¹B₁ states were calculated and compared with experiment. The first four triplet and four cationic states have also been examined.     

Effects of temperature and pressure on the near-infrared spectra of HOD in D2O

The near-infrared absorption spectra (9500 to 11000 cm^–1) of HOD, 20 mol% in D₂O were measured at temperatures between 4 and 55°C and pressures up to 500 MPa. From the analysis of the spectra, the following conclusions are drawn. (1) At temperatures below about 38°C, the ice I-like bulky structure is destroyed to form the dense structure which reflects the high-pressure ice-like structure as the pressure is increased. (2) At temperatures above about 38°C, the bulky structure hardly remains at atmospheric pressure and the formation of dense structure proceeds monotonically with increasing pressure. The results and conclusion obtained in the present paper agrees with those obtained for pure H₂O water in the previous investigation.     

Quantum calculations of vibrational energies of H3O2- on an ab initio potential

We report a full-dimensional potential energy surface for H3O2-, based on fitting 66,965 ab initio electronic energies. A major feature of this potential is a barrier of roughly 200 cm-1 to internal rotation of the two hydroxyl groups about a line connecting the two oxygen atoms and the bridging hydrogen atom. The potential is used in calculations of vibrational energies, performed with the "Reaction Path" version of the code "MULTIMODE". The results are compared to recent infrared messenger experiments and are used to propose interpretations of the experimental results.     

The vibrational Stokes shift of water (HOD in D2O)

The vibrational Stokes shift of the OH stretching transition nu(OH) of water is the shift between the ground-state absorption and the excited-state (v=1) emission. A recent measurement on HOD in D(2)O solvent [S. Woutersen and H. J. Bakker, Phys. Rev. Lett. 83, 2077 (1999)] of a 70 cm(-1) redshift, and a subsequent calculation of a 57 cm(-1) redshift using equilibrium molecular dynamics simulations [C. P. Lawrence and J. L. Skinner, J. Chem. Phys. 117, 8847 (2002)] were in good agreement. We now report extensive measurements of the vibrational Stokes shift in HOD/D(2)O using an ultrafast IR pump, Raman probe method. The vibrational Stokes shift is seen to depend on the pump pulse frequency and on time delay; by varying these parameters it can be made to range from 112 to -32 cm(-1) (negative values indicate a blueshift in the excited state). The equilibrium vibrational Stokes shift is actually a negative rather than a positive quantity. Possible reasons for the disagreement between experiment and theory are briefly discussed.     

OH-stretch vibrational relaxation of HOD in liquid to supercritical D2O

The population relaxation of the OH-stretching vibration of HOD diluted in D2O is studied by time-resolved infrared (IR) pump-probe spectroscopy for temperatures of up to 700 K in the density range 12 1 OH stretching transition with a 200 fs laser pulse centered at approximately 3500 cm(-1). Above 400 K these spectra show no indication of spectral diffusion after pump-probe delays of 0.3 ps. Over nearly the entire density range and for sufficiently high temperatures (T > 360 K), the vibrational relaxation rate constant, kr, is strictly proportional to the dielectric constant, epsilon, of water. Together with existing molecular dynamics simulations, this result suggests a simple linear dependence of kr on the number of hydrogen-bonded D2O molecules. It is shown that, for a given temperature, an isolated binary collision model is able to adequately describe the density dependence of vibrational energy relaxation even in hydrogen-bonded fluids. However, dynamic hydrogen bond breakage and formation is a source of spectral diffusion and affects the nature of the measured kr. For sufficiently high temperatures when spectral diffusion is much faster than energy transfer, the experimentally observed decays correspond to ensemble averaged population relaxation rates. In contrast, when spectral diffusion and vibrational relaxation occur on similar time scales, as is the case for ambient conditions, deviations from the linear kr(epsilon) relation occur because the long time decay of the v = 1 population is biased to slower relaxing HOD molecules that are only weakly connected to the hydrogen bond network.     

Full-dimensional vibrational calculations for H5O2+ using an ab initio potential energy surface

We report quantum diffusion Monte Carlo (DMC) and variational calculations in full dimensionality for selected vibrational states of H(5)O(2) (+) using a new ab initio potential energy surface [X. Huang, B. Braams, and J. M. Bowman, J. Chem. Phys. 122, 044308 (2005)]. The energy and properties of the zero-point state are focused on in the rigorous DMC calculations. OH-stretch fundamentals are also calculated using "fixed-node" DMC calculations and variationally using two versions of the code MULTIMODE. These results are compared with infrared multiphoton dissociation measurements of Yeh et al. [L. I. Yeh, M. Okumura, J. D. Myers, J. M. Price, and Y. T. Lee, J. Chem. Phys. 91, 7319 (1989)]. Some preliminary results for the energies of several modes of the shared hydrogen are also reported.     

Analytical potential from ab initio calculations for the Fe+-H2O and Feo-H2O systems

Analytical intermolecular potentials for the Fe⁺?H₂O and Fe^o?H₂O systems have been determined from ab initio calculations. Interaction energies for a lot of points along the two potential energy surfaces were calculated using Huzinga's MINI ?2 basis set. The results obtained were fitted to an analytical function containing 11 adjustable parameters that we have already used with success for the Fe²⁺?H₂O system. The goodness of the generated intermolecular potentials is discussed.     

Weakly bound complexes of N2O: an ab initio theoretical analysis toward the design of N2O receptors

Ab initio calculations at MP2/6-311++G(2d,2p) and MP2/6-311++G(3df,3pd) computational levels have been used to analyze the interactions between nitrous oxide and a series of small and large molecules that act simultaneously as hydrogen bond donors and electron donors. The basis set superposition error (BSSE) and zero point energy (ZPE) corrected binding energies of small N2O complexes (H2O, NH3, HOOH, HOO*, HONH2, HCO2H, H2CO, HCONH2, H2CNH, HC(NH)NH2, SH2, H2CS, HCSOH, HCSNH2) vary between -0.93 and -2.90 kcal/mol at MP2/6-311++G(3df,3pd) level, and for eight large complexes of N2O they vary between -2.98 and -3.37 kcal/mol at the MP2/6-311++G(2d,2p) level. The most strongly bound among small N2O complexes (HCSNH2-N2O) contains a NH..N bond, along with S-->N interactions, and the most unstable (H2S-N2O) contains just S-->N interactions. The electron density properties have been analyzed within the atoms in molecules (AIM) methodology. Results of the present study open a window into the nature of the interactions between N2O with other molecular moieties and open the possibility to design N2O abiotic receptors.     

New ab initio potential energy surface and the vibration-rotation-tunneling levels of (H2O)2 and (D2O)2

We report a new full-dimensional potential energy surface (PES) for the water dimer, based on fitting energies at roughly 30,000 configurations obtained with the coupled-cluster single and double, and perturbative treatment of triple excitations method using an augmented, correlation consistent, polarized triple zeta basis set. A global dipole moment surface based on Moller-Plesset perturbation theory results at these configurations is also reported. The PES is used in rigorous quantum calculations of intermolecular vibrational frequencies, tunneling splittings, and rotational constants for (H2O)2 and (D2O)2, using the rigid monomer approximation. Agreement with experiment is excellent and is at the highest level reported to date. The validity of this approximation is examined by comparing tunneling barriers within that model with those from fully relaxed calculations.     

Comment on gas phase infrared spectrum and ab initio calculations of phosphorus(III) thiocyanide, SPCN

We show that two of the three bands assigned to SPCN are assignable to cyanogen and that the bands assigned to SPCl are done so in error, with the 712 cm(-1) band in the spectrum of Allaf and Odeh assignable to HCN. There is no evidence for either ClPS or SPCN in the spectra shown. Finally low resolution IR spectroscopy by itself, whilst useful in assisting in the identification of pyrolysis products does not provide unambiguous identification and requires support by rigorous computation, band modelling and correct use of the literature. None of these are evidenced in the present paper of Allaf and Odeh nor indeed in the work that was the subject of our previous criticism and re-analysis.     

Interpretation of the microwave spectrum of 2-methoxy ethylamine using its ab initio structures

The 7₃ ← 6₃ μ_aK₋₁ doublet line series of 2-methoxy ethylamine, previously recorded with the radiofrequency microwave double resonance technique, was interpreted for the T and G conformations of this compound with a number of structural models including the 4-21G optimized geometries. A single model of G can reproduce the doublet series as the vibrational progression of the CO skeletal mode, starting with the υ_CO = 0 groundstate at 35 457 MHz, and stepping in approximately 95 MHz intervals up to the vibrational satellite of υ_CO = 9, with the υ_CO = 5 doublet missing. When model calculations are performed for T with the structural parameters of G which reproduce the 7₃ ← 6₃K₋₁ doublet series, except for rotating the NH₂ group by 120° to obtain T, the frequency of its groundstate doublet is found at lower frequencies than that of G. In contrast to this, when the empirically corrected 4-21G parameters of T are used in the analysis, the calculated groundstate of T coincides with the expected υ_CO = 6 doublet of G. The 4-21G geometries of G and T are in good agreement with the rotational constants and the conformational energy difference derived from the microwave spectra. Thus, this analysis clearly confirms the assignment of the observed 7₃ ← 6₃K₋₁ series as the progressions of the υ_CO = 0 to υ_CO = 4 states of G, and the υ_CO = 0 to υ_CO = 3 states of T.     

Full-dimensional quantum dynamics study of exchange processes for the D + H2O and D + HOD reactions

The exchange processes of D + H(2)O and D + HOD reactions are studied using initial state-selected time-dependent wave packet approach in full dimension. The total reaction probabilities for different partial waves, together with the integral cross sections, are obtained both by the centrifugal sudden (CS) approximation and exact coupled-channel (CC) calculations, for the H(2)O(HOD) reactant initially in the ground rovibrational state. In the CC calculations, small resonance peaks in the reaction probabilities and quick diminishing of the resonance peaks with the increase of total angular momenta J do not lead to clear step-like features just above the threshold in the cross sections for the title reactions, which are different in other isotopically substituted reactions where the hydrogen atom was included as the reactant instead of the deuterium atom [B. Fu, Y. Zhou, and D. H. Zhang, Chem. Sci. 3, 270 (2012); B. Fu and D. H. Zhang, J. Phys. Chem. A 116, 820 (2012)]. It is interesting that the shape resonance-induced features resulting from the reaction tunneling are significantly diminished accordingly in the reactions of the deuterium atom and H(2)O or HOD, owing to the weaker tunneling capability of the reagent deuterium atom in the title reactions than the reagent hydrogen atom in other reactions. In the CS calculations, the resonance peaks persist in many partial waves but cannot survive the partial-wave summations. The cross sections for the D(') + H(2)O → D(')OH + H and D(') + HOD → D(')OD + H reactions are substantially larger than those for the D(') + HOD → HOD(') + D reaction, indicating that the D(')/H exchange reactions are much more favored than the D(')/D exchange.