The calculation of force constants and normal coordinates III: constrained force fields for the methyl halides

A method is discussed for imposing any desired constraint on the force field obtained in a force constant refinement calculation. The application of this method to force constant refinement calculations for the methyl halide molecules is reported. All available data on the vibration frequencies, Coriolis interaction constants and centrifugal stretching constants of CH3X and CD3X molecules were used in the refinements, but despite this apparent abundance of data it was found that constraints were necessary in order to obtain a unique solution to the force field. The results of unconstrained calculations, and of three different constrained calculations, are reported in this paper. The constrained models reported are a Urey—Bradley force field, a modified valence force field, and a constraint based on orbital-following bond-hybridization arguments developed in the following paper. The results are discussed, and compared with previous results for these molecules. The third of the above models is found to reproduce the observed data better than either of the first two, and additional reasons are given for preferring this solution to the force field for the methyl halide molecules.

The ν2, 2ν2, 3ν2, ν4 and ν2+ν4 bands of 15NH3

The ir absorption of gaseous 15NH3 between 510 and 3040 cm−1 was recorded with a resolution of 0.06 cm−1. The ν2, 2ν2, 3ν2, ν4, and ν2 + ν4 bands were measured and analyzed on the basis of the vibration-rotation Hamiltonian developed by V. Špirko, J. M. R. Stone, and D. Papoušek (J. Mol. Spectrosc. 60, 159–178 (1976)). A set of effective molecular parameters for the ν2 = 1, 2, 3 states was derived, which reproduced the transition frequencies within the accuracy of the experimental measurements. For ν4 and ν2 + ν4 bands the standard deviation of the calculated spectrum is about four times larger than the measurements accuracy: a similar result was found for ν4 in 14NH3 by Š. Urban et al. (J. Mol. Spectrosc. 79, 455–495 (1980)). This result suggests that the present treatment takes into account only the most significant part of the rovibration interaction in the doubly degenerate vibrational states of ammonia.

Nitrous acid: vibrational frequency shifts due to 15N substitution, and the harmonic force field

Vibration rotation spectra of HO15 NO and DO15 NO have been measured at a resolution of 0•04 cm-1 to determine the isotopic shifts in the vibrational band origins. These have been used together with recently determined data on the vibrational band origins, Coriolis constants, and centrifugal distorition constants, to determine the harmonic force field of both cis and trans nitrous acid in least squares refinement calculations. The results are discussed in relation to recent ab initio calculations, the inertia defects, and the torsional potential function.

Second order Coriolis resonance in the v2, v5 Raman bands of CH3F

The vibration-rotation Raman spectrum of the ν2 and ν5 fundamentals of CH3F is reported, from 1320 to 1640 cm−1, with a resolution of about 0.3 cm−1. The Coriolis resonance between the two bands leads to many perturbation-allowed transitions. Where the resonance is still sufficiently weak that the quantum number K′ retains its meaning, perturbation-allowed transitions are observed for all values of ΔK from +4 to −4; in regions of strong resonance, however, we can only say that the observed transitions obey the selection rule Δ(k−l) = 0 or ±3. The spectrum has been analyzed by band contour simulation using a computer program based on exact diagonalization of the Hamiltonian within the ν2, ν5 vibrational levels, and improved vibration-rotation constants for these bands are reported. The relative magnitudes and relative sings of polarizability derivatives involved in these vibrations are also reported.

Perturbations in the infrared spectrum of monofluoroacetylene, and their relationship to intramolecular vibrational energy redistribution

Several high-order vibration-rotation perturbations in the high-resolution infrared spectrum of monofluoroacetylene, HCCF, are assigned and analyzed in detail. They result in avoided crossings in the rotational structure of several bands, and precise values for the effective high-order terms in the Hamiltonian have been determined. The significance of these results for intramolecular vibrational redistribution is discussed.

Intensities in infrared bands V: Bond polar properties in dimethylacetylene

The absolute intensities of all except one of the infra-red fundamental vibration bands of dimethyl acetylene have been determined, and the results have been used to compute polar properties of the C—H and C—C bonds. It has been found that if the very probable assumption is made that the acetylenic carbon atoms carry a residual negative charge, the hydrogen atoms in the C—H bonds must carry a residual positive charge. The probable value of the C—H dipole is about 04 Debye, and that of the C—C bond about 1 Debye. Comparisons have been made with the results of similar work with related molecules.

Microwave spectrum of CF3CCH

Attempts to observe ΔK = ±1 transitions in the rotational spectrum of CF3CCH and CF3H in the first excited state of a degenerate vibration, by direct absorption in the Y band and K band regions of the microwave spectrum, have not been successful. In the course of this work the J = 3-2 and 4-3 rotational spectrum of CF3CCH has been observed with higher sensitivity than previously, and from the positions of the vibrational satellites several new rB values have been determined.

The microwave spectrum of thietan-2,2,4,4-d4

The microwave spectrum for thietan-2,2,4,4-d4 is analysed in six of its lowest puckering states and up to J = 25. The close lying pairs of states of vp = 0/1, 2/3 and 4/5 are treated with a vibration-rotation hamiltonian which includes an off-diagonal coupling term in vp. Additional corrections to this coupling term in higher powers of the angular momentum operator are derived and their importance for improving the fit of calculated to observed data is tested. The variation of the centrifugal distortion constants with vp follows the model of Creswell and Mills (1974, J. molec. Spectrosc., 52, 392). A value is determined for the derivative with respect to the puckering coordinate of the ac-component of the inverse moment of inertia tensor.

Signs of l-doubling constants in NF3

Computed infrared band contours are presented for the two degenerate fundamentals of the NF3 molecule, using the l-resonance band contour program described by Cartwright and Mills, with values of the vibration-rotation constants determined from the microwave spectrum by Morino and co-workers. Computed contours are presented for both possible signs of the l-doubling constant, and comparison with the observed band contours leads to the conclusion that q3 = −121.4 MHz and q4 = +51.4 MHz.

Carbon suboxide: the vibrational dependence of the v7 bending potential function

Data on the vibrational energy levels and rotational constants of carbon suboxide for the low-wavenumber bending mode ν7 are reviewed, in the ground-state manifold, and in the ν2-, ν3-, ν4-, and ν2 + ν4-state manifolds. Following the procedure developed by Duckett, Mills, and Robiette [J. Mol. Spectrosc. 63, 249 (1976)] the data have been inverted to give the effective bending potential in ν7 for each of these five states. Values are obtained for various other parameters in the effective vibration-rotation Hamiltonian. The potential and rotational constants in ν2 + ν4 are given to a close approximation by linear extrapolation from the ground state through the ν2 and ν4 states.

The microwave spectrum of 1-pyrazoline

The microwave spectrum of 1-pyrazoline has been observed from 18 to 40 GHz in the six lowest states of the ring-puckering vibration. It is an a-type spectrum of a near oblate asymmetric top. Each vibrational state has been fitted to a separate effective Hamiltonian, and the vibrational dependence of both the rotational constants and the quartic centrifugal distortion constants has been observed and analyzed. The v = 0 and 1 states have also been analyzed using a coupled Hamiltonian; this gives consistent results, with an improved fit to the high J data. The preferred choice of Durig et al. [J. Chem. Phys. 52, 6096 (1970)] for the ring-puckering potential is confirmed as essentially correct, but the A and B inertial axes are shown to be interchanged from those assumed by Durig et al. in their analysis of the mid-infrared spectrum.

ASYM20: a program for force constant and normal coordinate calculations, with a critical review of the theory involved

The theory of harmonic force constant refinement calculations is reviewed, and a general-purpose program for force constant and normal coordinate calculations is described. The program, called ASYM20. is available through Quantum Chemistry Program Exchange. It will work on molecules of any symmetry containing up to 20 atoms and will produce results on a series of isotopomers as desired. The vibrational secular equations are solved in either nonredundant valence internal coordinates or symmetry coordinates. As well as calculating the (harmonic) vibrational wavenumbers and normal coordinates, the program will calculate centrifugal distortion constants, Coriolis zeta constants, harmonic contributions to the α′s. root-mean-square amplitudes of vibration, and other quantities related to gas electron-diffraction studies and thermodynamic properties. The program will work in either a predict mode, in which it calculates results from an input force field, or in a refine mode, in which it refines an input force field by least squares to fit observed data on the quantities mentioned above. Predicate values of the force constants may be included in the data set for a least-squares refinement. The program is written in FORTRAN for use on a PC or a mainframe computer. Operation is mainly controlled by steering indices in the input data file, but some interactive control is also implemented.

The calculation of force constants and normal co-ordinates--VII Normal co-ordinates in the E species of allene

Force constant and normal co-ordinate calculations are reported for the E species vibrations of the allene molecule. Data on the fundamental vibration frequencies of allene-h4, allene-d4 and allene-1.1-d2 and on the five experimentally determined Coriolis zeta constants of C3H4 and C3D4, were used in a force constant refinement procedure. Allowing for product and sum rules this gives 21 independent data which were used to refine to the most general harmonic force field (10 parameters) with one constraint (in the absence of any constraints the refinement was not satisfactory). The results have been used to calculate the complete ζz Coriolis interaction matrix for the allene-1.1-d2 molecule, and hence to calculate the expected rotational structure of the perpendicular bending vibrations of this molecule; the good agreement obtained with the observed spectra is a check on our results.

Potential energy functions for ground state surfaces of HCO and HNO

Analytical potential functions are reported for the ground state surfaces of HCO and HNO, the functions being derived from spectroscopic and ab initio data. Harmonized force fields have been deduced for the stable configurations of both molecules and vibration frequencies predicted for the metastable species COH and NOH.

Rovibrational energy levels of hydrogen peroxide, studied by MULTIMODE with a reaction path Hamiltonian

Recently. Carter and Handy [J. Chem. Phys. 113 (2000) 987] have introduced the theory of the reaction path Hamiltonian (RPH) [J. Chem. Phys. 72 (1980) 99] into the variational scheme MULTIMODE, for the calculation of the J = 0 vibrational levels of polyatomic molecules, which have a single large-amplitude motion. In this theory the reaction path coordinate s becomes the fourth dimension of the moment-of-inertia tensor, and must be treated separately from the remaining 3N - 7 normal coordinates in the MULTIMODE program. In the modified program, complete integration is performed over s, and the M-mode MULTIMODE coupling approximation for the evaluation of the matrix elements applies only to the 3N - 7 normal coordinates. In this paper the new algorithm is extended to the calculation of rotational-vibration energy levels (i.e. J > 0) with the RPH, following from our analogous implementation for rigid molecules [Theoret. Chem. Acc. 100 (1998) 191]. The full theory is given, and all extra terms have been included to give the exact kinetic energy operator. In order to validate the new code, we report studies on hydrogen peroxide (H2O2), where the reaction path is equivalent to torsional motion. H2O2 has previously been studied variationally using a valence coordinate Hamiltonian; complete agreement for calculated rovibrational levels is obtained between the previous results and those from the new code, using the identical potential surface. MULTIMODE is now able to calculate rovibrational levels for polyatomic molecules which have one large-amplitude motion. (C) 2003 Elsevier B.V. All rights reserved.