Equilibrium structures of selenium compounds: The torsionally flexible molecule of selenophenol.

The equilibrium structure of selenophenol has been investigated using rotational spectroscopy and high-level quantum mechanical calculations, offering electronic and structural insight into the scarcely studied selenium compounds. The jet-cooled broadband microwave spectrum was measured in the 2-8 GHz cm-wave region using broadband (chirped-pulse) fast-passage techniques. Additional measurements up to 18 GHz used narrow-band impulse excitation. Spectral signatures were obtained for six isotopic species of selenium (80Se, 78Se, 76Se, 82Se, 77Se, and 74Se), together with different monosubstituted 13C species. The (unsplit) rotational transitions associated with the non-inverting µa-dipole selection rules could be partially reproduced with a semirigid rotor model. However, the internal rotation barrier of the selenol group splits the vibrational ground state into two subtorsional levels, doubling the dipole-inverting µb transitions. The simulation of the double-minimum internal rotation gives a very low barrier height (B3PW91: 42 cm-1), much smaller than for thiophenol (277 cm-1). A monodimensional Hamiltonian then predicts a huge vibrational separation of 72.2 GHz, justifying the non-observation of µb transitions in our frequency range. The experimental rotational parameters were compared with different MP2 and density functional theory calculations. The equilibrium structure was determined using several high-level ab initio calculations. A final Born-Oppenheimer (reBO) structure was obtained at the coupled-cluster CCSD(T)_ae/cc-wCVTZ level of theory, including small corrections for the wCVTZ → wCVQZ basis set enlargement calculated at the MP2 level. The mass-dependent method with predicates was used to produce an alternative rm(2) structure. The comparison between the two methods confirms the high accuracy of the reBO structure and offers information on other chalcogen-containing molecules.

Determination of the semiexperimental equilibrium structure of 2-acetylthiophene in the presence of methyl internal rotation and substituent effects compared to thiophene.

The microwave spectra of thiophene and 2-acetylthiophene were recorded in the frequency range from 2 to 40 GHz using two molecular jet Fourier transform microwave spectrometers. For 2-acetylthiophene, two conformers with a syn and an anti orientation of the S1-C2 and C6O bonds (with respect to the C2-C6 bond) were identified, and the syn-conformer was more stable. The spectra of the 34S- and 13C-isotopologues of syn-2-acetylthiophene were also assigned, and the semiexperimental equilibrium structure could be determined. Compared to thiophene, at the substitution position, the S1-C2 and C2C3 bond lengths both increase by about 0.007 Å, and the bond angle S1-C2C3 decreases by 0.06°, noticeably larger than the experimental uncertainties. A-E torsional splittings were observed due to internal rotation of the methyl group hindered by a barrier height of 330.187(35) and 295.957(17) cm-1 for the syn-conformer and the anti-conformer, respectively. Geometry and internal rotation parameters are compared with those of related thiophene derivatives, as well as those of furan and 2-acetylthiophene to gain a better understanding of structure determination in the presence of methyl internal rotation.

Equilibrium Structure in the Presence of Methyl Internal Rotation: Microwave Spectroscopy and Quantum Chemistry Study of the Two Conformers of 2-Acetylfuran.

For 2-acetylfuran, quantum chemistry predicted and proton magnetic resonance study reported two conformers, anti and syn, differing in the position of the carbonyl group with respect to the O1-C2 bond of the furan ring. The microwave spectrum of the title molecule was recorded in the frequency range from 2 to 26.5 GHz using a molecular jet Fourier transform microwave spectrometer, confirming the presence of both conformers. Spectroscopic parameters such as the rotational and centrifugal distortion constants could be determined with high precision. The spectra of all 13C- and 18O-isotopologues of the energetically more favorable anti-conformer could be assigned, allowing the experimental determination of bond lengths and bond angles from the heavy atom substitution rs and the semi-experimental equilibrium reSE structures. Splittings arising from the internal rotation of the acetyl methyl group could be resolved for both conformers as well as for all assigned isotopologues, from which the barrier to methyl internal rotation was determined. The torsional barrier is largely invariant at around 319 cm-1 in the parent species of anti-2-acetylfuran and its isotopologues, showing that though isotopic substitution greatly influences the rotational properties of the molecule and causes a different microwave spectrum, its effect on the methyl torsion is negligible. On the other hand, conformational effects play a decisive role, as the torsional barrier of 239.780(13) cm-1 found for syn-2-acetylfuran differs significantly from the value for anti-2-acetylfuran. The results are compared and discussed with other methyl-substituted furan derivatives and acetyl group containing ketones for a better understanding of different effects influencing molecular geometry parameters and methyl internal rotations.

How accurate is the determination of equilibrium structures for van der Waals complexes? The dimer N2O⋯CO as an example.

Plausible methods for accurate determination of equilibrium structures of intermolecular clusters have been assessed for the van der Waals dimer N2O⋯CO. In order to assure a large initial dataset of rotational parameters, we first measured the microwave spectra of the 15N2O⋯12CO and 15N2O⋯13CO isotopologs, expanding previous measurements. Then, an anharmonic force field was calculated ab initio and a semi-experimental equilibrium structure was determined. The dimer structure was also calculated at the coupled-cluster level of theory using very large basis sets with diffuse functions and counterpoise correction. It was found that the contributions of the diffuse functions and the counterpoise correction are not additive and do not compensate each other although they have almost the same value but opposite signs. The semi-experimental and ab initio structures were found to be in fair agreement, with the equilibrium distance between the centers of mass of both monomers being 3.825(13) Å and the intermolecular bond length r(C⋯O) = 3.300(9) Å. In this case, the mass-dependent method did not permit us to determine reliable intermolecular parameters. The combination of experimental rotational constants and results of ab initio calculations thus proves to be very sensitive to examine the accuracy of structural determinations in intermolecular clusters, offering insight into other aggregates.

Accurate equilibrium structure of 3-aminophthalimide from gas electron diffraction and coupled-cluster computations and diverse structural effects due to electron density transfer.

For the first time, the molecular structure of 3-aminophthalimide has been determined by the gas electron diffraction (GED) method supported by a mass-spectrometric analysis of the gas phase and results of quantum-chemical computations up to coupled-cluster level of theory, CCSD(T). The semiexperimental equilibrium structure, rsee, has been derived from the GED data by taking into account harmonic and anharmonic vibrational corrections estimated from the quantum-chemical force field (up to cubic terms). High accuracy structures have been exploited for the observation of fine structural effects arising due the presence of the electron-donating amino group and the formation of a hydrogen bond. Natural bond orbital (NBO) analysis and quantum theory of atoms in molecules (QTAIM) have been applied to explain these effects.

The puzzling hyper-fine structure and an accurate equilibrium geometry of succinic anhydride.

An accurate semiexperimental equilibrium structure of succinic anhydride has been determined from a combination of experiment and theory. The cm-wave and mm-wave rotational spectra of succinic anhydride, 3,4-dihydrofuran-2,5-dione, were recorded in a pulsed supersonic jet using Fourier-transform microwave spectroscopy and in a free-jet using mm-wave absorption spectroscopy. Many lines in the cm-wave spectrum show fine structure and after eliminating all other possibilities the origin of this fine structure is determined to be from spin-spin interaction. Accurate rotational and quartic centrifugal distortion constants are determined. Assignments of 13C and 18O singly substituted isotopologues in natural abundance were used to obtain a substitution geometry for the heavy atoms of succinic anhydride. Theoretical approaches permitted the calculation of a Born-Oppenheimer ab initio structure and the determination of a semiexperimental equilibrium structure in which computed rovibrational corrections were utilized to convert vibrational ground state rotational constants into equilibrium constants. The agreement between the semiexperimental structure and the Born-Oppenheimer ab initio structure is excellent. Succinic anhydride has been shown to have a planar heavy atom equilibrium structure with the effects of a large amplitude vibration apparent in the resultant rotational constants.

The S-S Bridge: A Mixed Experimental-Computational Estimation of the Equilibrium Structure of Diphenyl Disulfide.

The disulfide bridge (-S-S-) is an important structural motif in organic and protein chemistry, but only a few accurate equilibrium structures are documented. We report the results of supersonic-jet microwave spectroscopy experiments on the rotational spectra of diphenyl disulfide, C6 H5 -S-S-C6 H5 (including all 13 C and 34 S monosubstituted isotopologues), and the determination of the equilibrium structure by the mixed estimation (ME) method. A single conformation of C2 symmetry was observed in the gas phase. This disulfide is a challenging target since its structure is determined by 34 independent parameters. Additionally, ab initio calculations revealed the presence of three low-frequency vibrations (<50 cm-1 ) associated to phenyl torsions which would prevent the calculation of an accurate force field. For this reason, instead of the semiexperimental method, we used the mass-dependent (rm ) method to fit the structural parameters concurrently to moments of inertia and predicate parameters, affected with appropriate uncertainties. The predicates were obtained by high-level quantum-chemical computations. A careful analysis of the results of different fits and a comparison with the ab initio optimizations confirms the validity of the used methods, providing detailed structural information on the title compound and the disulfide bridge.

How flexible is the disulfide linker? A combined rotational-computational investigation of diallyl disulfide.

The symmetrically substituted diallyl disulfide adopts a non-symmetric conformation in the gas-phase, as observed with supersonic-jet rotational spectroscopy. The determination of the equilibrium structure with a predicate mixed regression illustrates both the benefits of the mass-dependent method for moderately large molecules and the structural peculiarities of the disulfide bridge.

Experiment and theory at the convergence limit: accurate equilibrium structure of picolinic acid by gas-phase electron diffraction and coupled-cluster computations.

The accurate molecular structure of picolinic acid has been determined from experimental data and computed at the coupled cluster level of theory. Only one conformer with the O[double bond, length as m-dash]C-C-N and H-O-C[double bond, length as m-dash]O fragments in antiperiplanar (ap) positions, ap-ap, has been detected under conditions of the gas-phase electron diffraction (GED) experiment (Tnozzle = 375(3) K). The semiexperimental equilibrium structure, rsee, of this conformer has been derived from the GED data taking into account the anharmonic vibrational effects estimated from the ab initio force field. The equilibrium structures of the two lowest-energy conformers, ap-ap and ap-sp (with the synperiplanar H-O-C[double bond, length as m-dash]O fragment), have been fully optimized at the CCSD(T)_ae level of theory in conjunction with the triple-ζ basis set (cc-pwCVTZ). The quality of the optimized structures has been improved due to extrapolation to the quadruple-ζ basis set. The high accuracy of both GED determination and CCSD(T) computations has been disclosed by a correct comparison of structures having the same physical meaning. The ap-ap conformer has been found to be stabilized by the relatively strong NH-O hydrogen bond of 1.973(27) Š(GED) and predicted to be lower in energy by 16 kJ mol-1 with respect to the ap-sp conformer without a hydrogen bond. The influence of this bond on the structure of picolinic acid has been analyzed within the Natural Bond Orbital model. The possibility of the decarboxylation of picolinic acid has been considered in the GED analysis, but no significant amounts of pyridine and carbon dioxide could be detected. To reveal the structural changes reflecting the mesomeric and inductive effects due to the carboxylic substituent, the accurate structure of pyridine has been also computed at the CCSD(T)_ae level with basis sets from triple- to 5-ζ quality. The comprehensive structure computations for pyridine as well as for carbon dioxide have been used to examine the convergence with respect to the basis set size.

Semiexperimental and mass-dependent structures by the mixed regression method: Accurate equilibrium structure and failure of the Kraitchman method for ethynylcyclohexane.

The mixed regression method for determination of molecular structures is reviewed and applied to the investigation of ethynylcyclohexane, using both semiexperimental and mass-dependent methods. This methodology provides an efficient and computationally affordable route to obtain accurate molecular reference data, preventing ill-conditioning in the structural least-squares determinations from experimental rotational constants. New supersonic-jet microwave measurements are reported to obtain inertial data for the axial and equatorial species of ethynylcyclohexane, together with all 13C isotopologues of the equatorial form. The semiexperimental equilibrium (reSE) and mass-dependent (rm(2)) structures of the molecule are compared with high-level ab initio optimizations, showing that both methods deliver compatible structures with accuracies of about 0.002 Å for bond lengths and 0.2° for bond angles. We confirm that dependable predicates can be obtained for a large variety of bonds. Finally, we verify that the substitution method completely fails to determine a reliable structure for the title compound.

N-Methyl Inversion and Accurate Equilibrium Structures in Alkaloids: Pseudopelletierine.

A rotational spectroscopy investigation has resolved the conformational equilibrium and structural properties of the alkaloid pseudopelletierine. Two different conformers, which originate from inversion of the N-methyl group from an axial to an equatorial position, have been unambiguously identified in the gas phase, and nine independent isotopologues have been recorded by Fourier-transform microwave spectroscopy in a jet expansion. Both conformers share a chair-chair configuration of the two bridged six-membered rings. The conformational equilibrium is displaced towards the axial form, with a relative population in the supersonic jet of Naxial /Nequatorial ≈2/1. An accurate equilibrium structure has been determined by using the semiexperimental mixed-estimation method and alternatively computed by quantum-chemical methods up to the coupled-cluster level of theory. A comparison with the N-methyl inversion equilibria in related tropanes is also presented.

Axial-equatorial isomerism and semiexperimental equilibrium structures of fluorocyclohexane.

An experimental-computational methodology combining rotational data, high-level ab initio calculations and predicate least-squares fitting is applied to the axial-equatorial isomerism and semiexperimental equilibrium structure determination of fluorocyclohexane. New supersonic-jet microwave measurements of the rotational spectra of the two molecular conformations, together with all 13C isotopologues of both isomeric forms are reported. Equilibrium rotational constants are obtained from the ground-state rotational constants corrected for vibration-rotation interactions and electronic contributions. Equilibrium structures were determined by the mixed estimation (ME) method. Different computational methods were tested for the evaluation of predicate values of the structural parameters, and a computationally effective procedure for estimating reliable dihedral angles is proposed. Structural parameters were fitted concurrently to predicate parameters and moments of inertia, affected with appropriate uncertainties. The new structures of the title compound are regarded as accurate to 0.001 Å and 0.2°, illustrating the advantages of this methodology. Structural comparisons are offered with the cyclohexane prototype, revealing subtle substituent effects. For comparison purposes the equilibrium structures for the two fluorocyclohexane isomers and cyclohexanone are computed from high-level ab initio theory with inclusion of adjustments for basis set dependence and correlation of the core electrons.

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

Fructose and deoxyribose (24 and 19 atoms, respectively) are too large for determining accurate equilibrium structures, either by high-level ab initio methods or by experiments alone. We show in this work that the semiexperimental (SE) mixed estimation (ME) method offers a valuable alternative for equilibrium structure determinations in moderate-sized molecules such as these monosaccharides or other biochemical building blocks. The SE/ME method proceeds by fitting experimental rotational data for a number of isotopologues, which have been corrected with theoretical vibration-rotation interaction parameters (αi), and predicate observations for the structure. The derived SE constants are later supplemented by carefully chosen structural parameters from medium level ab initio calculations, including those for hydrogen atoms. The combined data are then used in a weighted least-squares fit to determine an equilibrium structure (r). We applied the ME method here to fructose and 2-deoxyribose and checked the accuracy of the calculations for 2-deoxyribose against the high level ab initio r structure fully optimized at the CCSD(T) level. We show that the ME method allows determining a complete and reliable equilibrium structure for relatively large molecules, even when experimental rotational information includes a limited number of isotopologues. With a moderate computational cost the ME method could be applied to larger molecules, thereby improving the structural evidence for subtle orbital interactions such as the anomeric effect.

Accurate Determination of Equilibrium Structure of 3-Aminophthalonitrile by Gas Electron Diffraction and Coupled-Cluster Computations: Structural Effects Due to Intramolecular Charge Transfer.

For the first time, the molecular structure of 3-aminophthalonitrile with unique electronic properties has been determined by the gas electron diffraction (GED) method supported by a mass spectrometric analysis of the gas phase. Moreover, it has been optimized at the high-level quantum-chemical coupled-cluster theory, CCSD(T), in conjunction with the triple-ζ basis set. The equilibrium structure has been determined from the GED data taking into account harmonic and anharmonic vibrational corrections estimated from the quantum-chemical force field (up to cubic terms). The computed CCSD(T) structure has been corrected for the core-core and core-valence electron-correlation effects estimated at the MP2 level and extrapolated to the basis set of quadruple-ζ quality. A remarkable agreement between the experimental and theoretical equilibrium structural parameters (bond lengths and angles) points to a high accuracy of both the molecular structure and applied theories. The high accuracy of structure computations and experimental determination allows the observation of structural changes due to the intramolecular charge transfer predicted by natural bond orbital (NBO) calculations. According to the NBO analysis, the amino group is the electron-donating substituent, whereas the nitrile groups are able to withdraw the π electrons from the benzene ring. Noticeable variations in the structural parameters are also explained by the interaction of σ → π* orbitals of the nearest C≡N and N-H bonds.

Interplay of experiment and theory: high resolution infrared spectrum and accurate equilibrium structure of BF2OH.

The high-resolution Fourier transform infrared (FTIR) spectrum of (11)BF2OH (difluoroboric acid) is analyzed taking into account numerous interactions. The ν1, ν2 and ν3 infrared bands are analyzed for the first time, whereas the parameters of the 6(1), 7(1), 8(1) and 9(1) states and for the 4(1) and 9(2) interacting states are redetermined. These results are used to check the quality of the ab initio force field. It is found that the ab initio rovibrational corrections are more accurate than the experimental ones. An earlier attempt to determine a semiexperimental structure did not allow us to obtain an accurate equilibrium structure. The reasons of this failure are investigated. This failure was mainly due to the lack of useful experimental information. Indeed, there is no isotopic substitution available for the fluorine atoms, and the boron atom is extremely close to the center of mass. Furthermore, the available isotopic substitutions (H → D and (16)O → (18)O) induce a large rotation of the principal axis system which amplifies the errors. However, the mixed estimation method has allowed us to determine a complete and reliable equilibrium structure. Thanks to this method, it is possible to determine an accurate structure, even in extremely difficult cases. An extensive analysis of the quality of structure calculations at the CCSD(T) level is also performed using basis sets up to five ζ quality. It was found that, at the convergence limit, the effects of the diffuse functions are practically disappearing, whereas the core-core and core-valence electron correlation effects are quite important for the bond lengths.

Ab Initio Anharmonic Analysis of Vibrational Spectra of Uracil Using the Numerical-Analytic Implementation of Operator Van Vleck Perturbation Theory.

The numerical-analytic implementation of the operator version of the canonical Van Vleck second-order vibrational perturbation theory (CVPT2) is employed for a purely ab initio prediction and interpretation of the infrared (IR) and Raman anharmonic spectra of a medium-size molecule of the diketo tautomer of uracil (2,4(1H,3H)-pyrimidinedione), which has high biological importance as one of the four RNA nucleobases. A nonempirical, semidiagonal quartic potential energy surface (PES) expressed in normal coordinates was evaluated at the MP2/cc-pVTZ level of theory. The quality of the PES was improved by replacing the harmonic frequencies with the "best" estimated CCSD(T)-based values taken from the literature. The theoretical method is enhanced by an accurate treatment of multiple Fermi and Darling-Dennison resonances with evaluation of the corresponding resonance constants W and K (CVPT2+WK method). A prediction of the anharmonic frequencies as well as IR and Raman intensities was used for a detailed interpretation of the experimental spectra of uracil. Very good agreement between predicted and observed vibrational frequencies has been achieved (RMSD ∼4.5 cm(-1)). The model employed gave a theoretically robust treatment of the multiple resonances in the 1680-1790 cm(-1) region. Our new analysis gives the most reliable reassignments of IR and Raman spectra of uracil available to date.

From the determination of the accurate equilibrium structure of 1-methylthymine by gas electron diffraction and coupled cluster computations to the observation of methylation and flexibility effects in pyrimidine nucleobases.

1-Methylthymine is of particular interest because it can be considered as a simple model of thymidine, in which deoxyribose attaches to thymine precisely at the N1 atom. The structure of this molecule is still unknown and so far has been experimentally studied for the first time in this work. The semiexperimental equilibrium structural parameters (r(e)(se)/∠(e)(se)) of 1-methylthymine have been determined by the gas electron diffraction (GED) method, taking into account vibrational corrections calculated with the use of the anharmonic (cubic) MP2/cc-pVTZ force constants. The methyl torsion around the C­N bond has been treated as a large-amplitude motion. For the first time, the structure of this molecule has been optimized by the very time-consuming coupled-cluster method (CCSD(T)(ae)) with the triple-ζ (cc-pwCVTZ) basis set. The obtained results have been extrapolated to the quadruple-ζ basis set at the MP2 level. It has been revealed that the methylation of uracil, especially at the nitrogen atom, leads to an increase in the flexibility of the nucleobase as well as a noticeable deformation of the pyrimidine ring.

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.

Why it is sometimes difficult to determine the accurate position of a hydrogen atom by the semiexperimental method: structure of molecules containing the OH or the CH3 group.

The semiexperimental (SE) technique, whereby equilibrium rotational constants are derived from experimental ground-state rotational constants and corrections based on an ab initio cubic force field, has the reputation to be one of the most accurate methods to determine an equilibrium structure ( reSE). However, in some cases, it cannot determine accurately the position of the hydrogen. To investigate the origins of this difficulty, the SE structures of several molecules containing either the OH or the CH3 group are determined and compared to their best ab initio counterparts. It appears that an important factor is the accuracy of the geometry used to calculate the force field, in particular when the least-squares system is not well conditioned. In this case, the mixed regression method is often an easy way to circumvent this difficulty.

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.