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

The Semiexperimental Approach at Work: Equilibrium Structure of Radical Species

The so-called semiexperimental (SE) approach is a powerful technique for obtaining highly accurate equilibrium structures for isolated systems. This Featured Article describes its extension to open-shell species, thus providing the first systematic investigation on radical equilibrium geometries to be used for benchmarking purposes. The small yet significant database obtained demonstrates that there is no reduction in accuracy when moving from closed-shell species to radicals. We also provide an extension of the applicability of the SE approach to medium-/large-sized radicals by exploiting the so-called "Lego-brick" approach, which is based on the assumption that a molecular system can be seen as formed by smaller fragments for which the SE equilibrium structure is available. In this Featured Article we show that this model can be successfully applied also to open-shell species.

Toward Accurate Quantum Chemical Methods for Molecules of Increasing Dimension: The New Family of Pisa Composite Schemes

The new versions of the Pisa composite scheme introduced in the present paper are based on the careful selection of different quantum chemical models for energies, geometries, and vibrational frequencies, with the aim of maximizing the accuracy of the overall description while retaining a reasonable cost for all the steps. In particular, the computation of accurate electronic energies has been further improved introducing more reliable complete basis set extrapolations and estimation of core-valence correlation, together with improved basis sets for third-row atoms. Furthermore, the reduced-cost frozen natural orbital (FNO) model has been introduced and validated for large molecules. Accurate molecular structures can be obtained avoiding complete basis set extrapolation and evaluating core-valence correlation at the MP2 level. Unfortunately, analytical gradients are not available for the FNO version of the model. Therefore, for large molecules, an accurate reduced-cost alternative is offered by evaluation of valence contributions with a double-hybrid functional in conjunction with the same MP2 contribution for core-valence correlation or by means of a one-parameter approximation. The same double-hybrid functional and basis set are employed to evaluate zero-point energies and partition functions. After the validation of the new models for small systems, a panel of molecular bricks of life has been used to analyze their performances for problems of current fundamental or technological interest. The fully black-box implementation of the computational workflow paves the way toward the accurate yet not prohibitively expensive study of medium- to large-sized molecules also by experimentally oriented researchers.

Accurate Biomolecular Structures by the Nano-LEGO Approach: Pick the Bricks and Build Your Geometry

The determination of accurate equilibrium molecular structures plays a fundamental role for understanding many physical–chemical properties of molecules, ranging from the precise evaluation of the electronic structure to the analysis of the role played by dynamical and environmental effects in tuning their overall behavior. For small semi-rigid systems in the gas phase, state-of-the-art quantum chemical computations rival the most sophisticated experimental (from, for example, high-resolution spectroscopy) results. For larger molecules, more effective computational approaches must be devised. To this end, we have further enlarged the compilation of available semi-experimental (SE) equilibrium structures, now covering the most important fragments containing H, B, C, N, O, F, P, S, and Cl atoms collected in the new SE100 database. Next, comparison with geometries optimized by methods rooted in the density functional theory showed that the already remarkable results delivered by PW6B95 and, especially, rev-DSDPBEP86 functionals can be further improved by a linear regression (LR) approach. Use of template fragments (taken from the SE100 library) together with LR estimates for the missing interfragment parameters paves the route toward accurate structures of large molecules, as witnessed by the very small deviations between computed and experimental rotational constants. The whole approach has been implemented in a user-friendly tool, termed nano-LEGO, and applied to a number of demanding case studies.

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 N₂O⋯CO. In order to assure a large initial dataset of rotational parameters, we first measured the microwave spectra of the ¹⁵N₂O⋯¹²CO and ¹⁵N₂O⋯¹³CO 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.

Laboratory Observation of, Astrochemical Search for, and Structure of Elusive Erythrulose in the Interstellar Medium

Rotational spectroscopy provides the most powerful means of identifying molecules of biological interest in the interstellar medium (ISM), but despite their importance, the detection of carbohydrates has remained rather elusive. Here, we present a comprehensive Fourier transform rotational spectroscopic study of elusive erythrulose, a sugar building block likely to be present in the ISM, employing a novel method of transferring the hygroscopic oily carbohydrate into the gas phase. The high sensitivity of the experiment allowed the rotational spectra of all monosubstituted isotopologue species of ¹³C-¹²C₃H₈O₄ to be recorded, which, together with quantum chemical calculations, enabled us to determine their equilibrium geometries (r_e^SE) with great precision. Searches employing the new experimental data for erythrulose have been undertaken in different ISM regions, so far including the cold areas Barnard 1, the pre-stellar core TMC-1, Sagittarius B2. Although no lines of erythrulose were found, this data will serve to enable future searches and possible detections in other ISM regions.

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.

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

The symmetrically substituted diallyl disulfide adopts an asymmetric conformation in the gas phase, as observed in the rotational spectrum.

Competing Dispersive Interactions: From Small Energy Differences to Large Structural Effects in Methyl Jasmonate and Zingerone

Modern structural studies of biologically relevant molecules require an exhaustive interplay between experiment and theory. In this work, we present two examples where a poor choice of the theoretical method led to a misinterpretation of experimental results. We do that by performing a rotational spectroscopy study on two large and flexible biomolecules: methyl jasmonate and zingerone. The results show the enormous potential of rotational spectroscopy as a benchmark to evaluate the performance of theoretical methods.

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.

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 N_axial /N_equatorial ≈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.

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.

Axial–equatorial isomerism and semiexperimental equilibrium structures of fluorocyclohexane

Accuracy vs. computational effort: the mixed estimation method provides very accurate semiexperimental equilibrium structures, as illustrated by the subtle differences between axial/equatorial fluorocyclohexane.