Accuracy of the rotational constants

Accurate structures and rotational constants of bicyclic monoterpenes at DFT cost by means of the bond-corrected Pisa composite scheme (BPCS)

The structural, conformational, and spectroscopic properties in the gas phase of 20 bicyclic monoterpenes and monoterpenoids have been analyzed by a new accurate, reduced-cost computational strategy. In detail, the revDSD-PBEP86 double-hybrid functional in conjunction with the D3BJ empirical dispersion corrections and a suitable triple-zeta basis set provides accurate geometrical parameters, whence equilibrium rotational constants, which are further improved by proper account of core-valence correlation. Average deviations within 0.1% between computed and experimental rotational constants are reached when taking into account the vibrational corrections obtained by the B3LYP functional in conjunction with a double-zeta basis set in the framework of second-order vibrational perturbation theory. In addition to their intrinsic interest, the studied terpenes further extend the panel of systems for which the proposed strategy has provided accurate results at density functional theory cost. Therefore, a very accurate yet robust and user-friendly tool is now available for systematic investigations of the role of stereo-electronic effects on the properties of large systems of current technological and/or biological interest by experimentally oriented researchers.

Equilibrium Values for the Si-H Bond Length and Equilibrium Structures of Silyl Iodide and Halosilylenes

The equilibrium structures of silyl iodide, SiH3I, and silylene halides, SiHX (X = F, Cl, Br, I), were determined by using the mixed regression method, where approximate values of the rotational constants are supplemented by the structural parameters of a different origin. For this goal, it is shown that the r(Si-H) bond length can be determined by using the isolated SiH stretching frequency and that an accurate estimation of the bond angles is obtained by an MP2 calculation with a basis set of triple zeta quality. To check the accuracy of the experimental structures, they were also optimized by means of all electron CCSD(T) calculations using basis sets of quadruple zeta quality.

Unbiased Comparison between Theoretical and Experimental Molecular Structures and Properties: Toward an Accurate Reduced-Cost Evaluation of Vibrational Contributions

The tremendous development of hardware and software is constantly increasing the role of quantum chemical (QC) computations in the assignment and interpretation of experimental results. However, an unbiased comparison between theory and experiment requires the proper account of vibrational averaging effects. In particular, high-resolution spectra in the gas phase are now available for molecules containing up to about 50 atoms, which are too large for a brute-force approach with the available QC methods of sufficient accuracy. In the present paper, we introduce hybrid approaches, which allow the accurate evaluation of vibrational averaging effects for molecules of this size beyond the harmonic approximation, with special attention being devoted to rotational constants. After the validation of new tools for relatively small molecules, the β-estradiol hormone and a prototypical molecular motor have been considered to witness the feasibility of accurate computations for large molecules.

Accurate Structures and Rotational Constants of Steroid Hormones at DFT Cost: Androsterone, Testosterone, Estrone, β-Estradiol, and Estriol

A comprehensive analysis of the structural, conformational, and spectroscopic properties in the gas phase has been performed for five prototypical steroid hormones, namely, androsterone, testosterone, estrone, β-estradiol, and estriol. The revDSD-PBEP86 double-hybrid functional in conjunction with the D3BJ empirical dispersion and a suitable triple-ζ basis set provides accurate conformational energies and equilibrium molecular structures, with the latter being further improved by proper account of core-valence correlation. Average deviations within 0.1% between computed and experimental ground state rotational constants are reached when adding to those equilibrium values vibrational corrections obtained at the cost of standard harmonic frequencies thanks to the use of a new computational tool. Together with the intrinsic interest of the studied hormones, the accuracy of the results obtained at DFT cost for molecules containing about 50 atoms paves the way toward the accurate investigations of other flexible bricks of life.

Describing nuclear quantum effects in vibrational properties using molecular dynamics with Wigner sampling

In this work we discuss the generally applicable Wigner sampling and introduce a new, simplified Wigner sampling method, for computationally effective modeling of molecular properties containing nuclear quantum effects and vibrational anharmonicity. For various molecular systems test calculations of (a) vibrationally averaged rotational constants, (b) vibrational IR spectra and (c) photoelectron spectra have been performed. The performance of Wigner sampling has been assessed by comparing with experimental data and with results of other theoretical models, including harmonic and VPT2 approximations. The developed simplified Wigner sampling method shows advantages in application to large and flexible molecules.

Fourier Transform Microwave Spectrum of Propene-3-d1 (CH2═CHCH2D), Quadrupole Coupling Constants of Deuterium, and a Semiexperimental Equilibrium Structure of Propene

The ground-state rotational spectrum of propene-3-d₁, CH₂═CHCH₂D, was measured by Fourier transform microwave spectroscopy. Transitions were assigned for the two conformers, one with the D atom in the symmetry plane (S) and the other with the D atom out of the plane (A). The energy difference between the two conformers was calculated to be 6.5 cm^-1, the S conformer having lower energy. The quadrupole hyperfine structure due to deuterium was resolved and analyzed for both conformers. The experimental quadrupole coupling and the centrifugal distortion constants compared favorably to their ab initio counterparts. Ground-state rotational constants for the S conformer are 40582.157(9), 9067.024(1), and 7766.0165(12) MHz. Ground-state rotational constants for the A conformer are 43403.75(3), 8658.961(2), and 7718.247(2) MHz. For the A conformer, a small tunneling splitting (19 MHz) due to internal rotation was observed and analyzed. Using the new rotational constants of this work as well as those previously determined for the ¹³C species and for some deuterium-substituted species from the literature, a new semiexperimental equilibrium structure was determined and its high accuracy was confirmed. The difficulty in obtaining accurate coordinates for the out-of-plane hydrogen atom is discussed.

The equilibrium molecular structures of 2-deoxyribose and fructose by the semiexperimental mixed estimation method and coupled-cluster computations

The mixed estimation method yields very accurate equilibrium structures, as observed in deoxyribose and fructose.

Electron delocalization in polyenes: a semiexperimental equilibrium structure for (3E)-1,3,5-hexatriene and theoretical structures for (3Z,5Z)-, (3E,5E)-, and (3E,5Z)-1,3,5,7-octatetraene

Electronic structure theory reveals that π-electron delocalization increases with the chain length in polyenes. To analyze quantitatively this effect a semiexperimental equilibrium structure has been determined for trans-hexatriene by the mixed estimation method. For this fit rotational constants for a number of carbon and hydrogen isotopologues as well as a high-level ab initio structure have been used. The accuracy is 0.001 Å for bond lengths and 0.1° for bond angles. For the three isomers of octatetraene, high-level ab initio calculations have given a comparably accurate structure. These structures have been used in comparison with the structure of s-trans-butadiene to show that "C═C" bonds increase in length and "C-C" bonds decrease in length as the polyene chain lengthens. These structural effects of π-electron delocalization increase toward the center of polyenes. Most likely, π-π conjugation in the molecules studied plays a large part in their planarity that, in turn, forces the hydrogen atoms of cis fragments in bay regions to be in a close contact. Their distance is indeed shorter than the sum of their van der Waals radii, and they seem to participate in a six-membered ring.

Accurate determination of the deformation of the benzene ring upon substitution: equilibrium structures of benzonitrile and phenylacetylene

Accurate equilibrium, re, structures of the monosubstituted benzene molecules benzonitrile, C6H5CN, and phenylacetylene, C6H5CCH, have been determined using two different, to some extent complementary techniques. The semiexperimental, r(e)(SE), structural parameters are the result of a least-squares fit to equilibrium rotational constants derived from experimental effective ground-state rotational constants and rovibrational corrections based principally on an ab initio cubic force field. The composite ab initio Born-Oppenheimer, r(e)(BO), structural parameters are obtained from frozen-core and all-electron MP2 and the CCSD(T) geometry optimizations using Gaussian basis sets up to quintuple-zeta quality. The DFT(B3LYP) method, with two different Gaussian basis sets, 6-31G* and 6-311+G(3df,2pd), was used to calculate the cubic force field employed during the r(e)(SE) structure determination. With the 6-31G* basis set, the error of the rovibrational correction is to a large extent random, whereas with the 6-311+G(3df,2pd) basis set it is mainly systematic. As shown here, systematic errors do not have a significant effect on the accuracy of the derived structure; the quality of the structural fit, however, is sensitive to the true accuracy of the ground-state rotational constants. An even more important general conclusion of this study is that the addition of extra rotational constants from multisubstituted species does not seem to improve the accuracy of the r(e)(SE) structures, quite in contrast to the highly desirable availability of data corresponding to all singly substituted species.

Interplay of experiment and theory: determination of an accurate equilibrium structure of 1-methyluracil by the gas electron diffraction method and coupled-cluster computations

As far as fundamental knowledge is concerned, the methyl derivatives of uracil can be considered as the simplest objects for studying the structural effects due to the substitution in the pyrimidyne nucleobases. From this point of view, 1-methyluracil is of special importance in biochemistry because uracil attaches ribose in ribonucleic acid (RNA) just precisely at the N1 atom. The semi-experimental equilibrium structure (r(e)(se)) of 1-methyluracil has been determined for the first time by the gas electron diffraction (GED) method taking into account rovibrational corrections to the thermal-average internuclear distances calculated with harmonic and anharmonic (cubic) MP2/cc-pVTZ force constants with consideration of the methyl torsion as a large-amplitude motion. For the first time, the structure of the molecule has been optimized by the very time-consuming coupled-cluster method with single and double excitations and perturbative treatment of connected triples using the correlation-consistent polarized weighted core-valence triple-ζ basis set with all electrons being correlated (CCSD(T)(all)/cc-pwCVTZ) and extrapolated to the complete basis set (CBS) with the help of the MP2 calculations. Small differences between similar bond lengths of equilibrium configurations were assumed in the GED analysis at the CCSD(T)(all)/CBS values. A remarkable agreement between the semi-experimental and computed equilibrium structures points out the high accuracy of both the GED determination and the coupled-cluster computations. The effect of methylation on the structure of uracil has been analyzed.

Microwave spectra of the deuterium isotopologues of cis-hexatriene and a semiexperimental equilibrium structure

Microwave transitions and ground state rotational constants are reported for five newly synthesized deuterium isotopologues of cis-1,3,5-hexatriene (cHTE). These rotational constants along with those of the parent and the three (13)C species are used with vibration-rotation constants calculated from an MP2/cc-pVTZ model to derive an equilibrium structure. That structure is improved by the mixed estimation method. In this method, internal coordinates from good-quality quantum chemical calculations (with appropriate uncertainties) are fit simultaneously with moments of inertia of the full set of isotopologues. The new structure of cHTE is confirmed to be planar and is stabilized by an interaction between the hydrogen atoms H2 and H5, which form a bond and participate in a six-membered ring. cHTE shows larger structural effects of π-electron delocalization than does butadiene with the effects being magnified in the center of the molecule. Thus, strong structural evidence now exists for an increase in π-electron delocalization as the polyene chain lengthens.

Rotational spectroscopy meets theory

Rotational spectroscopy is known to be a technique that is widely used to infer information on molecular structure and dynamics. In the last few decades, its role in the field of atmospheric and astrophysical investigations has rapidly grown. However, several are the challenging aspects in rotational spectroscopy, since the detection and analysis of spectra as well as interpretation of obtained results are not at all straightforward. Quantum chemistry has reached such an accuracy that can be used to disentangle these challenging situations by guiding the experimental investigation, assisting in the determination of the spectroscopic parameters, and extracting information of chemical interest. This perspective provides an overview of the theoretical background and computational requirements needed for the accurate evaluation of the spectroscopic parameters of relevance to rotational spectroscopy. The role of theory in guiding and supporting experiment is detailed through a few examples, and the interplay of experiment and theory is discussed in terms of the information of physical and chemical interest that can be derived.