Accurate Spectroscopic Models for Methane Polyads Derived from a Potential Energy Surface Using High-Order Contact Transformations

Analysis of experimental spectra of phosphine in the Tetradecad range near 2.3 μm using ab initio calculations

Due to its major interest for the chemistry of planetary atmospheres and exobiology, accurate spectroscopy data of phosphine are required for the search of signatures of this molecule in astronomical observations. In this work, high resolution infrared laboratory spectra of phosphine were analyzed for the first time in the full Tetradecad region (3769-4763 cm-1) involving 26 rotationally resolved bands. Overall, 3242 lines were assigned in spectra previously recorded by Fourier transform spectroscopy at temperatures 200 K and 296 K, using a combined theoretical model based on ab initio calculations. The total nuclear motion Hamiltonian of PH3 including ab initio potential energy surface, was reduced to an effective Hamiltonian using the high-order contact transformation method adapted to vibrational polyads of the AB3 symmetric top molecules, followed by empirical optimization of the parameters. At this step, the experimental line positions were reproduced with a standard deviation of 0.0026 cm-1 that provided unambiguous identification of observed transitions. The effective dipole transition moments of the bands were obtained by fitting to the intensities obtained from variational calculations using the ab initio dipole moment surface. The assigned lines were used to newly determine 1609 experimental vibration-rotational levels up to Jmax = 18 with energy in the range 3896-6037 cm-1 that represents significant extension towards higher energies compared to previous works. Transitions for all 26 sublevels of the Tetradecad were identified but with noticeably fewer transitions for fourfold excited bands because of their weaker intensity. At the final step, pressure-broadened half widths were attached to each transition and a composite line list adopting ab initio intensities and empirical line positions corrected to the accuracy of about 0.001 cm-1 for strong and medium transitions was validated against experimental spectra available in the literature.

Vibrational resonance analysis of linear molecules using resummation of divergent Rayleigh-Schrödinger perturbation theory series

The large order Rayleigh-Schrödinger perturbation theory (RSPT) was applied for calculating vibrational states of linear molecules. Two molecules (CO₂ and C₂H₂) were used as test cases with using of isomorphic Watson Hamiltonian and quartic force fields. For CO₂ the Sayvetz condition can remove all degeneracies for purely vibrational states and the non-degenerate perturbation theory can be applied. However, an existence of two degenerate modes in C₂H₂ requires using the upgraded degenerate version of RSPT that was employed in this context for the first time. The dominating divergent behavior of such series requires the resummation technique that mimics the multivalued nature of the underlying solutions, and the applied quartic Padé-Hermite approximants (QPHA) provided full solution of the problem. Moreover, some mathematical properties of QPHA proved to be an efficient tool for studying resonance effects through the Katz theorem that controls the singular points of the eigenvalues on the complex plane. In the case of C₂H₂, not only all earlier observed classical resonances were confirmed and quantified, but also subtle interpolyad resonances (K_2/55,K_3/4555), proposed recently by Herman (2011) were described as well. Following the analysis, we found several novel resonances, of which we proposed one independent interpolyad resonance K_2/4444. The complete analysis of such critical points provided the full resonance picture of all studied molecules.

Novel methodology for systematically constructing global effective models from ab initio-based surfaces: A new insight into high-resolution molecular spectra analysis

In this paper, a novel methodology is presented for the construction of ab initio effective rotation-vibration spectroscopic models from potential energy and dipole moment surfaces. Non-empirical effective Hamiltonians are obtained via the block-diagonalization of selected variationally computed eigenvector matrices. For the first time, the derivation of an effective dipole moment is carried out in a systematic way. This general approach can be implemented quite easily in most of the variational computer codes and turns out to be a clear alternative to the rather involved Van Vleck perturbation method. Symmetry is exploited at all stages to translate first-principles calculations into a set of spectroscopic parameters to be further refined on experiment. We demonstrate on H₂CO, PH₃, CH₄, C₂H₄, and SF₆ that the proposed effective model can provide crucial information to spectroscopists within a very short time compared to empirical spectroscopic models. This approach brings a new insight into high-resolution spectrum analysis of polyatomic molecules and will be also of great help in the modeling of hot atmospheres where completeness is important.

Derivation of ρ-dependent coordinate transformations for nonrigid molecules in the Hougen-Bunker-Johns formalism

In this paper, we report a series of transformations for the construction of a Hamiltonian model for nonrigid polyatomic molecules in the framework of the Hougen-Bunker-Johns formalism (HBJ). This model is expressed in normal mode coordinates for small vibrations and in a specific coordinate ρ to describe the large amplitude motion. For the first time, a general procedure linking the "true" curvilinear coordinates to ρ is proposed, allowing the expression of the potential energy part in the same coordinate representation as the kinetic energy operator, whatever the number of atoms. A Lie group-based method is also proposed for the derivation of the reference configuration in the internal axis system. This work opens new perspectives for future high-resolution spectroscopy studies of nonrigid, medium-sized molecules using HBJ-type Hamiltonians. Illustrative examples and computation of vibrational energy levels on semirigid and nonrigid molecules are given to validate this method.

Isotopic and symmetry breaking effects on phosphine spectra under H → D substitutions from ab initio variational calculations

Variationally computed infrared spectra in the range [0-5000] cm^-1 are reported for the deuterated PH₂D and PHD₂ molecules from accurate potential energy and dipole moment surfaces initially derived for the major isotopologue PH₃( C 3 v ). Energy level and line intensity calculations were performed by using a normal-mode model combined with isotopic and symmetry transformations for the H → D substitutions. Theoretical spectra were computed at 296 K up to J _max = 30 and will be made available through the TheoReTS information system (http://theorets.univ-reims.fr, http://theorets.tsu.ru). For the very first time, ab initio intensity predictions of PH₂D/PHD₂ are in good qualitative agreement with the literature. This work will be useful for spectral intensity analysis for which accurate spectral intensity data are still missing.

Highly excited vibrational levels of methane up to 10 300 cm-1: Comparative study of variational methods

In this work, we report calculated vibrational energy levels of the methane molecule up to 10 300 cm^-1. Two potential energy surfaces constructed in quite different coordinate systems with different analytical representations are employed in order to evaluate the uncertainty of vibrational predictions. To calculate methane energy levels, we used two independent techniques of the variational method. One method uses an exact kinetic energy operator in internal curvilinear coordinates. Another one uses an expansion of Eckart-Watson nuclear motion Hamiltonian in rectilinear normal coordinates. In the Icosad range (up to five vibrational quanta bands-below 7800 cm^-1), the RMS standard deviations between calculated and observed energy levels were 0.22 cm^-1 and 0.41 cm^-1 for these two quite different approaches. For experimentally well-known 3v₃ sub-levels, the calculation accuracy is estimated to be ∼1 cm^-1. In the Triacontad range (7660-9188 cm^-1), the average error of the calculation is about 0.5 cm^-1. The accuracy and convergence issues for higher energy ranges are discussed.

Understanding global infrared opacity and hot bands of greenhouse molecules with low vibrational modes from first-principles calculations: the case of CF4

Fluorine containing molecules have a particularly long atmospheric lifetime and their very big estimated global warming potentials are expected to rapidly increase in the future. This work is focused on the global theoretical prediction of infrared spectra of the tetrafluoromethane molecule that is considered as a potentially powerful greenhouse gas having the largest estimated lifetime of over 50 000 years in the atmosphere. The presence of relatively low vibrational frequencies makes the Boltzmann population of the excited levels important. Consequently, the "hot bands" corresponding to transitions among excited rovibrational states contribute significantly to the CF4 opacity in the infrared even at room temperature conditions but the existing laboratory data analyses are not sufficiently complete. In this work, we construct the first accurate and complete ab initio based line lists for CF4 in the range 0-4000 cm-1, containing rovibrational bands that are the most active in absorption. An efficient basis set compression method was applied to predict more than 700 new bands and subbands via variational nuclear motion calculations. We show that already at room temperature a quasi-continuum of overlapping weak lines appears in the CF4 infrared spectra due to the increasing density of bands and transitions. In order to converge the infrared opacity at room temperature, it was necessary to include a high rotational quantum number up to J = 80 resulting in 2 billion rovibrational transitions. In order to make the cross-section simulation faster, we have partitioned our data into two parts: (a) strong & medium line lists with lower energy levels for calculation of selective absorption features that can be used at various temperatures and (b) compressed "super-line" libraries of very weak transitions contributing to the quasi-continuum modelling. Comparisons with raw previously unassigned experimental spectra showed a very good accuracy for integrated absorbance in the entire range of the reported spectra predictions. The data obtained in this work will be made available through the TheoReTS information system (http://theorets.univ-reims.fr, http://theorets.tsu.ru) that contains ab initio born line lists and provides a user-friendly graphical interface for a fast simulation of the CF4 absorption cross-sections and radiance under various temperature conditions from 80 K to 400 K.

Anharmonic vibrational analysis of s-trans and s-cis conformers of acryloyl fluoride using numerical-analytic Van Vleck operator perturbation theory

A new gas-phase infrared (IR) spectrum of acryloyl fluoride (ACRF, CH₂CHCFO) with a resolution of 0.1cm^-1 in the range 4000-450cm^-1 was measured. Theoretical ab initio molecular structures, full quartic potential energy surfaces (PES), and cubic surfaces of dipole moments and polarizability tensor components (electro-optical properties, EOP) of the s-trans and s-cis conformers of the ACRF were calculated by the second-order Møller-Plesset electronic perturbation theory with a correlation consistent Dunning triple-ζ basis set. The numerical-analytic implementation of the second-order operator canonical Van Vleck perturbation theory was employed for predicting anharmonic IR and Raman scattering (RS) spectra of ACRF. To improve the anharmonic predictions, harmonic frequencies were replaced by their counterparts evaluated with the higher-level CCSD(T)/cc-pVTZ model, to form a "hybrid" PES. The original operator representation of the Hamiltonian is analytically reduced to a quasi-diagonal form, integrated in the harmonic oscillator basis and diagonalized to account for strong resonance couplings. Double canonical transformations of EOP expansions enabled prediction of integral intensities of both fundamental and multi-quanta transitions in IR/RS spectra. Enhanced band shape analysis reinforced the assignments. A thorough interpretation of the new IR experimental spectra and existing matrix-isolation literature data for the mixture of two conformers of ACRF was accomplished, and a number of assignments clarified.

Ab initio variational predictions for understanding highly congested spectra: rovibrational assignment of 108 new methane sub-bands in the icosad range (6280–7800 cm−1)

This work demonstrates for the first time how accurate first principles global calculations allow assigning complicated spectra of a molecule with more than 4 atoms.

[In Process Citation]

A novel model for the analysis of the rovibrational structure of a molecule based on anharmonic ladder operators associated with the vibrational degrees of freedom is presented. This is devised as an alternative method for the global spectral analysis of rovibrational data considering vibrational anharmonicities from the outset. The present method is thought up with an effective rovibrational Hamiltonian written in terms of angular momentum components and anharmonic Morse ladder operators, associated with rotational and vibrational degrees of freedom, respectively. The resulting Hamiltonian is diagonalized in a symmetry-adapted basis set expressed as a product of rotational states and individual 1D-Morse wave functions for each local vibrational degree of freedom. This approach has been successfully applied to the study of the vibrational structure (up to polyad 14) and the rovibrational structure (up to polyad 2 and J(max) = 20) of hydrogen sulfide. It was shown that this new global analysis formalism is able to reduce considerably the number of fitted parameters with respect to the spectral analysis carried out for separate polyads.

Does the "reef structure" at the ozone transition state towards the dissociation exist? New insight from calculations and ultrasensitive spectroscopy experiments

Since the discovery of anomalies in ozone isotope enrichment, several fundamental issues in the dynamics linked to the shape of the potential energy surface in the transition state region have been raised. The role of the reeflike structure on the minimum energy path is an intricate question previously discussed in the context of chemical experiments. In this Letter, we bring strong arguments in favor of the absence of a submerged barrier from ultrasensitive laser spectroscopy experiments combined with accurate predictions of highly excited vibrations up to nearly 95% of the dissociation threshold.