Lowest-Lying Conformers of Alanine: Pushing Theory to Ascertain Precise Energetics and Semiexperimental Re Structures

Factorized Quadruples and a Predictor of Higher-Level Correlation in Thermochemistry

Coupled cluster theory has had a momentous impact on the ab initio prediction of molecular properties, and remains a staple ingratiate in high-accuracy thermochemical model chemistries. However, these methods require inclusion of at least some connected quadruple excitations, which generally scale at best as O ( N 9 ) with the number of basis functions. It is very difficult to predict, a priori, the effect correlation of past CCSD(T) on a given reaction energy. The purpose of this work is to examine cost-effective quadruple corrections based on the factorization theorem of the many-body perturbation theory that may address these challenges. We show that the O ( N 7 ) factorized CCSD(TQf) method introduces minimal error to predicted correlation and reaction energies as compared to the O ( N 9 ) CCSD(TQ). Further, we examine the performance of Goodson's continued fraction method in the estimation of CCSDT(Q)Λ contributions to reaction energies as well as a "new" method related to %TAE[(T)] that we refer to as a scaled perturbation estimator. We find that the scaled perturbation estimator based upon CCSD(TQf)/cc-pVDZ is capable of predicting CCSDT(Q)Λ/cc-pVDZ contributions to reaction energies with an average error of 0.07 kcal mol-1 and an L2D of 0.52 kcal mol-1 when applied to a test-suite of nearly 3000 reactions. This offers a means by which to reliably "ballpark" how important post-CCSD(T) contributions are to reaction energies while incurring no more than CCSD(T) formal cost and a little mental math.

Sub 20 cm-1 computational prediction of the CH bond energy - a case of systematic error in computational thermochemistry

The bond dissociation energy of methylidyne, D0(CH), is studied using an improved version of the High-Accuracy Extrapolated ab initio Thermochemistry (HEAT) approach as well as the Feller-Peterson-Dixon (FPD) model chemistry. These calculations, which include basis sets up to nonuple (aug-cc-pCV9Z) quality, are expected to be capable of providing results substantially more accurate than the ca. 1 kJ mol-1 level that is characteristic of standard high-accuracy protocols for computational thermochemistry. The calculated 0 K CH bond energy (27 954 ± 15 cm-1 for HEAT and 27 956 ± 15 cm-1 for FPD), along with equivalent treatments of the CH ionization energy and the CH+ dissociation energy (85 829 ± 15 cm-1 and 32 946 ± 15 cm-1, respectively), were compared to the existing benchmarks from Active Thermochemical Tables (ATcT), uncovering an unexpected difference for D0(CH). This has prompted a detailed reexamination of the provenance of the corresponding ATcT benchmark, allowing the discovery and subsequent correction of a systematic error present in several published high-level calculations, ultimately yielding an amended ATcT benchmark for D0(CH). Finally, the current theoretical results were added to the ATcT Thermochemical Network, producing refined ATcT estimates of 27 957.3 ± 6.0 cm-1 for D0(CH), 32 946.7 ± 0.6 cm-1 for D0(CH+), and 85 831.0 ± 6.0 cm-1 for IE(CH).

Linking Experimental and Ab Initio Thermochemistry of Adsorbates with a Generalized Thermochemical Hierarchy

Enthalpies of formation of adsorbates are crucial parameters in the microkinetic modeling of heterogeneously catalyzed reactions since they quantify the stability of intermediates on the catalyst surface. This quantity is often computed using density functional theory (DFT), as more accurate methods are computationally still too expensive, which means that the derived enthalpies have a large uncertainty. In this study, we propose a new error cancellation method to compute the enthalpies of formation of adsorbates from DFT more accurately through a generalized connectivity-based hierarchy. The enthalpy of formation is determined through a hypothetical reaction that preserves atomistic and bonding environments. The method is applied to a data set of 60 adsorbates on Pt(111) with up to 4 heavy (non-hydrogen) atoms. Enthalpies of formation of the fragments required for the bond balancing reactions are based on experimental heats of adsorption for Pt(111). The comparison of enthalpies of formation derived from different DFT functionals using the isodesmic reactions shows that the effect of the functional is significantly reduced due to the error cancellation. Thus, the proposed methodology creates an interconnected thermochemical network of adsorbates that combines experimental with ab initio thermochemistry in a single and more accurate thermophysical database.

High-Accuracy Heats of Formation for Alkane Oxidation: From Small to Large via the Automated CBH-ANL Method

It is generally challenging to obtain high-accuracy predictions for the heat of formation for species with more than a handful of heavy atoms, such as those of importance in standard combustion mechanisms. To this end, we construct the CBH-ANL approach and illustrate that, for a set of 194 alkane oxidation species, it can be used to produce ΔH_f(0 K) values with 2σ uncertainties of 0.2-0.5 kcal mol^-1. This set includes the alkanes, hydroperoxides, and alkyl, peroxy, and hydroperoxyalkyl radicals for 17 representative hydrocarbon fuels containing up to 10 heavy atoms with various degrees of branching in the alkane backbone. The CBH-ANL approach, automated in the QTC and AutoMech software suites, builds balanced chemical equations for the calculation of ΔH_f(0 K), in which the reference species may be up to five heavy atoms. The high-level ANL0 and ANL1 reference ΔH_f(0 K) values are further refined for even the largest of these reference species with a novel laddering approach. We perform a comprehensive quantification of the uncertainties for both the individual reference species (the largest of which is 0.15 kcal mol^-1) and the propagation of those uncertainties when used in the calculation of ΔH_f(0 K) for the 194 target species. We examine the sensitivity of the predicted ΔH_f(0 K) values to (i) electronic energies from various methods, including ωB97X-D/cc-pVTZ, B2PLYP-D3/cc-pVTZ, CCSD(T)-F12b/cc-pVDZ-F12//B2PLYP-D3/cc-pVTZ, and CCSD(T)-F12b/cc-pVTZ-F12//B2PLYP-D3/cc-pVTZ; (ii) the zero-point vibrational energies (ZPVEs), where we consider harmonic ZPVEs as well as two scaling-based estimates of the anharmonic ZPVEs, all implemented for both ωB97X-D/cc-pVTZ and B2PLYP-D3/cc-pVTZ calculations; (iii) the particular CBH-ANL scheme employed; and (iv) the procedure for choosing the reference conformer for the analyses. The discussion concludes with a summary of the estimated overall uncertainty in the predictions and a validation of the predictions for the alkane subset.

Benchmark Structures and Conformational Landscapes of Amino Acids in the Gas Phase: A Joint Venture of Machine Learning, Quantum Chemistry, and Rotational Spectroscopy

The accurate characterization of prototypical bricks of life can strongly benefit from the integration of high resolution spectroscopy and quantum mechanical computations. We have selected a number of representative amino acids (glycine, alanine, serine, cysteine, threonine, aspartic acid and asparagine) to validate a new computational setup rooted in quantum-chemical computations of increasing accuracy guided by machine learning tools. Together with low-lying energy minima, the barriers ruling their interconversion are evaluated in order to unravel possible fast relaxation paths. Vibrational and thermal effects are also included in order to estimate relative free energies at the temperature of interest in the experiment. The spectroscopic parameters of all the most stable conformers predicted by this computational strategy, which do not have low-energy relaxation paths available, closely match those of the species detected in microwave experiments. Together with their intrinsic interest, these accurate results represent ideal benchmarks for more approximate methods.

Resolving Discrepancies between State-of-the-Art Theory and Experiment for HO2 + HO2 via Multiscale Informatics

Recent high-level theoretical calculations predict a mild temperature dependence for HO₂ + HO₂ inconsistent with state-of-the-art experimental determinations that upheld the stronger temperature dependence observed in early experiments. Via MultiScale Informatics analysis of the theoretical and experimental data, we identified an alternative interpretation of the raw experimental data that uses HO₂ + HO₂ rate constants nearly identical to theoretical predictions─implying that the theoretical and experimental data are actually consistent, at least when considering the raw data from experimental studies. Similar analyses of typical signals from low-temperature experiments indicate that an HOOOOH intermediate─identified by recent theory but absent from earlier interpretations─yields modest effects that are smaller than, but may have contributed to, the scatter in data among different experiments. More generally, the findings demonstrate that modern chemical theories and experiments have progressed to a point where meaningful comparison requires joint consideration of their data simultaneously.

Structure, energetics, and spectroscopy of the chromophores of HHe+n, H2He+n, and He+n clusters and their deuterated isotopologues

The linear molecular ions H₂He⁺, HHe+2, and He+3 are the central units (chromophores) of certain He-solvated complexes of the H₂He+n, HHe+n, and He+n families, respectively. These are complexes which do exist, according to mass-spectrometry studies, up to very high n values. Apparently, for some of the H₂He+n and He+n complexes, the linear symmetric tetratomic H₂He+2 and the diatomic He+2 cations, respectively, may also be the central units. In this study, definitive structures, relative energies, zero-point vibrational energies, and (an)harmonic vibrational fundamentals, and, in some cases, overtones and combination bands, are established mostly for the triatomic chromophores. The study is also extended to the deuterated isotopologues D₂He⁺, DHe+2, and D₂He+2. To facilitate and improve the electronic-structure computations performed, new atom-centered, fixed-exponent, Gaussian-type basis sets called MAX, with X = T(3), Q(4), P(5), and H(6), are designed for the H and He atoms. The focal-point-analysis (FPA) technique is employed to determine definitive relative energies with tight uncertainties for reactions involving the molecular ions. The FPA results determined include the 0 K proton and deuteron affinities of the ⁴He atom, 14 875(9) cm^-1 [177.95(11) kJ mol^-1] and 15 229(8) cm^-1 [182.18(10) kJ mol^-1], respectively, the dissociation energies of the He+2 → He⁺ + He, HHe+2 → HHe⁺ + He, and He+3 → He+2 + He reactions, 19 099(13) cm^-1 [228.48(16) kJ mol^-1], 3948(7) cm^-1 [47.23(8) kJ mol^-1], and 1401(12) cm^-1 [16.76(14) kJ mol^-1], respectively, the dissociation energy of the DHe+2 → DHe⁺ + He reaction, 4033(6) cm^-1 [48.25(7) kJ mol^-1], the isomerization energy between the two linear isomers of the [H, He, He]⁺ system, 3828(40) cm^-1 [45.79(48) kJ mol^-1], and the dissociation energies of the H₂He⁺ → H+2 + He and the H₂He+2 → H₂He⁺ + He reactions, 1789(4) cm^-1 [21.40(5) kJ mol^-1] and 435(6) cm^-1 [5.20(7) kJ mol^-1], respectively. The FPA estimates of the first dissociation energy of D₂He⁺ and D₂He+2 are 1986(4) cm^-1 [23.76(5) kJ mol^-1] and 474(5) cm^-1 [5.67(6) kJ mol^-1], respectively. Determining the vibrational fundamentals of the triatomic chromophores with second-order vibrational perturbation theory (VPT2) and vibrational configuration interaction (VCI) techniques, both built around the Eckart-Watson Hamiltonian, proved unusually challenging. For the species studied, VPT2 has difficulties yielding dependable results, in some cases even for the fundamentals of the H-containing molecular cations, while carefully executed VCI computations yield considerably improved spectroscopic results. In a few cases unusually large anharmonic corrections to the fundamentals, on the order of 15% of the harmonic value, have been observed.

Accurate evaluation of combustion enthalpy by ab-intio computations

Accurate evaluation of combustion enthalpy is of high scientific and industrial importance. Although ab-initio computation of the heat of reactions is one of the promising and well-established approaches in computational chemistry, reliable and precise computation of heat of combustion reactions by ab-initio methods is surprisingly scarce in the literature. A handful of works carried out for this purpose report significant inconsistencies between the computed and experimentally determined combustion enthalpies and suggest empirical corrections to improve the accuracy of the ab-initio predicted data. The main aim of the present study is to investigate the reasons behind those reported inconsistencies and propose guidelines for a high-accuracy estimation of heat of reactions via ab-initio computations. We show comparably accurate prediction of combustion enthalpy of 40 organic molecules based on a DSD-PBEP86 double-hybrid density functional theory approach and CCSD(T)-F12 coupled-cluster computations, with mean unsigned errors with respect to experimental data being below 0.5% for both methods.

Probing Molecular Chirality of Ground and Electronically Excited States in the UV-vis and X-ray Regimes: An EOM-CCSD Study

We present several strategies for computing electronic circular dichroism (CD) spectra across different frequency ranges at the equation-of-motion coupled-cluster singles and doubles level of theory. CD spectra of both ground and electronically excited states are discussed. For selected cases, the approach is compared with coupled-cluster linear response results as well as time-dependent density functional theory. The extension of the theory to include the effect of spin-orbit coupling is presented and illustrated by calculations of X-ray CD spectra at the L-edge.

Spiers Memorial Lecture: theory of unimolecular reactions

One hundred years ago, at an earlier Faraday Discussion meeting, Lindemann presented a mechanism that provides the foundation for contemplating the pressure dependence of unimolecular reactions. Since that time, our...

A Comprehensive Picture of the Structures, Energies, and Bonding in the Alanine Dimers

High level quantum mechanical computations and extensive stochastic searches of the potential energy surfaces of the Alanine dimers uncover rich and complex structural and interaction landscapes. A total of 416 strongly bound (up 13.4 kcal mol^-1 binding energies at the DLPNO-CCSD(T)/6-311++G(d,p) level corrected by the basis set superposition error and by the zero point vibrational energies over B3LYP-D3 geometries), close energy equilibrium structures were located, bonded via 32 specific types of intermolecular contacts including Y⋅⋅⋅H-X primary and Y⋅⋅⋅H-C secondary hydrogen bonds, H⋅⋅⋅H dihydrogen contacts, and non conventional anti-electrostatic Y δ - ⋯ X δ - interactions. The putative global minimum is triply degenerate, corresponding to the structure of the common dimer of a carboxylic acid. All quantum descriptors of chemical bonding point to a multitude of weak individual interactions within each dimer, whose cumulative effect results in large binding energies and in an attractive fluxional wall of non-covalent interactions in the interstitial region between the monomers.

Structural and Vibrational Properties of Amino Acids from Composite Schemes and Double-Hybrid DFT: Hydrogen Bonding in Serine as a Test Case

The structures, relative stabilities, and vibrational wavenumbers of the two most stable conformers of serine, stabilized by the O-H···N, O-H···O═C and N-H···O-H intramolecular hydrogen bonds, have been evaluated by means of state-of-the-art composite schemes based on coupled-cluster (CC) theory. The so-called "cheap" composite approach (CCSD(T)/(CBS+CV)^MP2) allowed determination of accurate equilibrium structures and harmonic vibrational wavenumbers, also pointing out significant corrections beyond the CCSD(T)/cc-pVTZ level. These accurate results stand as a reference for benchmarking selected hybrid and double-hybrid, dispersion-corrected DFT functionals. B2PLYP-D3 and DSDPBEP86 in conjunction with a triple-ζ basis set have been confirmed as effective methodologies for structural and spectroscopic studies of medium-sized flexible biomolecules, also showing intramolecular hydrogen bonding. These best performing double-hybrid functionals have been employed to simulate IR spectra by means of vibrational perturbation theory, also considering hybrid CC/DFT schemes. The best overall agreement with experiment, with mean absolute error of 8 cm^-1, has been obtained by combining CCSD(T)/(CBS+CV)^MP2 harmonic wavenumbers with B2PLYP-D3/maug-cc-pVTZ anharmonic corrections. Finally, a composite scheme entirely based on CCSD(T) calculations (CCSD(T)/CBS+CV) has been employed for energetics, further confirming that serine II is the most stable conformer, also when zero-point vibrational energy corrections are included.

A molecular modeling of iodinated organic compounds

Thermodynamic quantities of some iodinated compounds have been studied by means of ab initio and density functional theory calculations. Since there are some limitations for heavy elements in the fifth row of periodic table such as the iodine atom, a high level of ab initio theory is developed to calculate the energy of iodinated compounds. The standard enthalpies of formation for alkyl iodides and other iodinated organic compounds have been computed and compared with available experimental data. The impact of iodine atoms on the change of enthalpies for deprotonation of iodobenzenes has been studied and a mathematical model is proposed for the effect of iodine on the change of deprotonation. The acidity of diiodophenols as a model for thyroxin hormone has been successfully calculated. The barriers for nucleophilic reactions of substituted methyl iodides are investigated and the energy profiles are characterized by double-well potentials which correspond to pre-reaction and post-reaction complexes, separated by a transition state.

Shape and core-excited resonances in electron scattering from alanine

We present detailed ab initio scattering calculations using the R-matrix method for electron collisions with the most stable conformer of α-alanine. The shape resonances that we identify are in good agreement with earlier calculations and experiments. Core-excited and mixed-character resonances are identified and characterized computationally for the first time. Dissociative electron attachment results are discussed in relation to the resonances identified.

Thiols as Hydrogen Bond Acceptors and Donors: Spectroscopy of 2-Phenylethanethiol Complexes

Evidence and understanding of sulfur-centered hydrogen bonding, especially where the donor is a thiol, lags far behind that for conventional OH interactions. To help address this deficiency, conformer specific IR spectra of 2-phenylethanethiol (PET) and associated 1:1 solvent complexes have been measured in SH, OH, and CH stretch regions using resonant-two-photon-ionization (R2PI) and IR-UV ion dip spectroscopic techniques. The aromatic and aliphatic CH stretch regions show signature differences between anti and gauche conformers. Supported by ab initio calculations, a PET-water cluster with an OH···S arrangement and a PET-diethyl ether cluster expressing an SH···O interaction were identified. The SH stretch band of the SH···O complex is red-shifted and undergoes significant intensity enhancement compared to the bare molecule, which is characteristic of hydrogen bonding. These findings offer insight into the nature of the thiol functional group as a potential hydrogen bond donor and acceptor.

Definitive thermochemistry and kinetics of the interconversions among conformers of n-butane and n-pentane

The focal-point analysis (FPA) technique is used for the definitive characterization of conformational interconversion parameters, including activation energy barriers, activation free energies, and kinetic rate coefficients at 298 K, of two n-alkanes, n-butane, and n-pentane, yielding the first complete analysis of their interconversion kinetics. The FPA implementation developed in this study is based on geometry optimizations and harmonic frequency computations carried out with density functional theory methods and single-point energy computations up to the CCSD(T) level of electronic structure theory using atom-centered Gaussian basis sets as large as cc-pV5Z. The anharmonic vibrational computations are carried out, at the MP2/6-31G* level of theory. Reflecting the convergence behavior of the Gibbs free-energy terms and the interconversion parameters, well-defined uncertainties, mostly neglected in previous theoretical studies, are provided. Finally, the effect of these uncertainties on the concentrations of the conformers of n-butane and n-pentane is examined via a global Monte-Carlo uncertainty analysis. © 2017 Wiley Periodicals, Inc.

l-Alanine capping of ZnO nanorods: increased carrier concentration in ZnO/CuI heterojunction diode

ZnO nanorods were capped with a simple amino acid, viz., l-alanine to increase the carrier concentration and improve the performance of ZnO/CuI heterojunction diodes. The effect of l-alanine capping on the morphology, structural, optical, electrochemical and electrical properties of ZnO nanorods had been studied in detail. The stable structure with two equally strong Zn–O coordinate bonds predicted by density functional theory was in agreement with the experimental results of FTIR spectroscopy. Due to the presence of electron-releasing (+I effect) moieties in l-Alanine, the carrier concentration of capped ZnO nanorods was two orders of magnitude higher and the ZnO/CuI heterojunction device showed more than a two-fold increase in the photovoltaic power conversion efficiency.

ZnO nanorods were capped with a simple amino acid, viz., l-Alanine to increase the carrier concentration and improve the performance of ZnO/CuI heterojunction diodes.

Benchmarking Quantum Chemical Methods: Are We Heading in the Right Direction?

Theoreticians and experimentalists should work together more closely to establish reliable rankings and benchmarks for quantum chemical methods. Comparison to carefully designed experimental benchmark data should be a priority. Guidelines to improve the situation for experiments and calculations are proposed.

Impact of Hydrogen Bonding on the Susceptibility of Peptides to Oxidation

The tendency of peptides to be oxidized is intimately connected with their function and even their ability to exist in an oxidative environment. Here we report high-level theoretical studies that show that hydrogen bonding can alter the susceptibility of peptides to oxidation, with complexation to a hydrogen-bond acceptor facilitating oxidation, and vice versa, impacting the feasibility of a diverse range of biological processes. It can even provide an energetically viable mechanistic alternative to direct hydrogen-atom abstraction. We find that hydrogen bonding to representative reactive groups leads to a broad (≈400 kJ mol^-1 ) spectrum of ionization energies in the case of model amide, thiol and phenol systems. While some of the oxidative processes at the extreme ends of the spectrum are energetically prohibitive, subtle environmental and solvent effects could potentially mitigate the situation, leading to a balance between hydrogen bonding and oxidative susceptibility.

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.

Photochemical Chiral Symmetry Breaking in Alanine

We introduce a general theoretical approach for the simulation of photochemical dynamics under the influence of circularly polarized light to explore the possibility of generating enantiomeric enrichment through polarized-light-selective photochemistry. The method is applied to the simulation of the photolysis of alanine, a prototype chiral amino acid. We show that a systematic enantiomeric enrichment can be obtained depending on the helicity of the circularly polarized light that induces the excited-state photochemistry of alanine. By analyzing the patterns of the photoinduced fragmentation of alanine we find an inducible enantiomeric enrichment up to 1.7%, which is also in good correspondence to the experimental findings. Our method is generally applicable to complex systems and might serve to systematically explore the photochemical origin of homochirality.

First-principles data set of 45,892 isolated and cation-coordinated conformers of 20 proteinogenic amino acids

We present a structural data set of the 20 proteinogenic amino acids and their amino-methylated and acetylated (capped) dipeptides. Different protonation states of the backbone (uncharged and zwitterionic) were considered for the amino acids as well as varied side chain protonation states. Furthermore, we studied amino acids and dipeptides in complex with divalent cations (Ca²⁺, Ba²⁺, Sr²⁺, Cd²⁺, Pb²⁺, and Hg²⁺). The database covers the conformational hierarchies of 280 systems in a wide relative energy range of up to 4 eV (390 kJ/mol), summing up to a total of 45,892 stationary points on the respective potential-energy surfaces. All systems were calculated on equal first-principles footing, applying density-functional theory in the generalized gradient approximation corrected for long-range van der Waals interactions. We show good agreement to available experimental data for gas-phase ion affinities. Our curated data can be utilized, for example, for a wide comparison across chemical space of the building blocks of life, for the parametrization of protein force fields, and for the calculation of reference spectra for biophysical applications.

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.

A complex-polarization-propagator protocol for magneto-chiral axial dichroism and birefringence dispersion

A schematic representation of magneto-chiral effects.

Generalized Vibrational Perturbation Theory for Rotovibrational Energies of Linear, Symmetric and Asymmetric Tops: Theory, Approximations, and Automated Approaches to Deal with Medium-to-Large Molecular Systems

Models going beyond the rigid-rotor and the harmonic oscillator levels are mandatory for providing accurate theoretical predictions for several spectroscopic properties. Different strategies have been devised for this purpose. Among them, the treatment by perturbation theory of the molecular Hamiltonian after its expansion in power series of products of vibrational and rotational operators, also referred to as vibrational perturbation theory (VPT), is particularly appealing for its computational efficiency to treat medium-to-large systems. Moreover, generalized (GVPT) strategies combining the use of perturbative and variational formalisms can be adopted to further improve the accuracy of the results, with the first approach used for weakly coupled terms, and the second one to handle tightly coupled ones. In this context, the GVPT formulation for asymmetric, symmetric, and linear tops is revisited and fully generalized to both minima and first-order saddle points of the molecular potential energy surface. The computational strategies and approximations that can be adopted in dealing with GVPT computations are pointed out, with a particular attention devoted to the treatment of symmetry and degeneracies. A number of tests and applications are discussed, to show the possibilities of the developments, as regards both the variety of treatable systems and eligible methods.

Accurate molecular structures and infrared spectra of trans-2,3-dideuterooxirane, methyloxirane, and trans-2,3-dimethyloxirane

Oxirane derivatives are the most used benchmarks for chiroptical spectroscopies in view of their small size and relative rigidity. The molecular structure, vibrational harmonic and anharmonic frequencies, and infrared intensities of the ground electronic states are analyzed in this paper. Equilibrium structure and harmonic force fields have been evaluated by means of high-level quantum-chemical calculations at the coupled-cluster level including single and double excitations together with a perturbative treatment of triples (CCSD(T)). Extrapolation to the complete basis-set limit as well as core-correlation effects have also been taken into account. Anharmonic contributions have been computed at the CCSD(T)/cc-pVTZ level for trans-2,3-dideuterooxirane. These data can serve as references to evaluate the accuracy of less expensive computational approaches rooted in the density functional theory (DFT). The latter have been used within hybrid CC/DFT approaches, which have been applied to simulate fully anharmonic infrared (IR) spectra. Finally, the best theoretical estimates of the equilibrium structures and vibrational wavenumbers are compared to the most accurate experimental data and show in all cases very good agreement, i.e., within 0.001 Å, 0.1 deg, 10 cm(-1), and 0.5 km mol(-1), for bond lengths, angles, wavenumbers, and IR intensities, respectively.

Molecular structure and spectroscopic signatures of acrolein: theory meets experiment

A comprehensive study of the molecular structure and IR spectrum of cis and trans acrolein has been performed by an integrated computational approach coupling methods rooted in the coupled-cluster ansatz and the density functional theory. From the one side, DFT anharmonic force fields allow us to determine very reliable semiexperimental structures for both isomers, which are in remarkable agreement with the geometries issuing from CCSD(T) computations accounting for the extrapolation to the complete basis set and core correlation. The same kind of coupled-cluster computations provide dipole moment, relative energies, and interconversion barrier in remarkable agreement with experiments. Finally, harmonic CCSD(T) results coupled to DFT evaluation of mechanical and electrical anharmonicity allow us, in the framework of second-order perturbative vibrational theory, to confirm most of the experimental assignments of IR spectra, and to suggest some additional interpretations for congested regions including fundamental bands together with overtones and combination bands.

VUV photodynamics and chiral asymmetry in the photoionization of gas phase alanine enantiomers

The valence shell photoionization of the simplest proteinaceous chiral amino acid, alanine, is investigated over the vacuum ultraviolet region from its ionization threshold up to 18 eV. Tunable and variable polarization synchrotron radiation was coupled to a double imaging photoelectron/photoion coincidence (i(2)PEPICO) spectrometer to produce mass-selected threshold photoelectron spectra and derive the state-selected fragmentation channels. The photoelectron circular dichroism (PECD), an orbital-sensitive, conformer-dependent chiroptical effect, was also recorded at various photon energies and compared to continuum multiple scattering calculations. Two complementary vaporization methods-aerosol thermodesorption and a resistively heated sample oven coupled to an adiabatic expansion-were applied to promote pure enantiomers of alanine into the gas phase, yielding neutral alanine with different internal energy distributions. A comparison of the photoelectron spectroscopy, fragmentation, and dichroism measured for each of the vaporization methods was rationalized in terms of internal energy and conformer populations and supported by theoretical calculations. The analytical potential of the so-called PECD-PICO detection technique-where the electron spectroscopy and circular dichroism can be obtained as a function of mass and ion translational energy-is underlined and applied to characterize the origin of the various species found in the experimental mass spectra. Finally, the PECD findings are discussed within an astrochemical context, and possible implications regarding the origin of biomolecular asymmetry are identified.

Conformational analysis and intramolecular interactions of L-proline methyl ester and its N-acetylated derivative through spectroscopic and theoretical studies

This work reports a detailed study regarding the conformational preferences of L-proline methyl ester (ProOMe) and its N-acetylated derivative (AcProOMe) to elucidate the effects that rule their behaviors, through nuclear magnetic resonance (NMR) and infrared (IR) spectroscopies combined with theoretical calculations. These compounds do not present a zwitterionic form in solution, simulating properly amino acid residues in biological media, in a way closer than amino acids in the gas phase. Experimental (3)JHH coupling constants and infrared data showed excellent agreement with theoretical calculations, indicating no variations in conformer populations on changing solvents. Natural bond orbital (NBO) results showed that hyperconjugative interactions are responsible for the higher stability of the most populated conformer of ProOMe, whereas for AcProOMe both hyperconjugative and steric effects rule its conformational equilibrium.

Accurate Characterization of the Peptide Linkage in the Gas Phase: A Joint Quantum-Chemical and Rotational Spectroscopy Study of the Glycine Dipeptide Analogue

Accurate structures of aminoacids in the gas phase have been obtained by joint microwave and quantum-chemical investigations. However, the structure and conformational behavior of α-aminoacids once incorporated into peptide chains are completely different and have not yet been characterized with the same accuracy. To fill this gap, we present here an accurate characterization of the simplest dipeptide analogue (N-acetyl-glycinamide) involving peptidic bonds. State-of-the-art quantum-chemical computations are complemented by a comprehensive study of the rotational spectrum using a combination of Fourier transform microwave spectroscopy with laser ablation. The coexistence of the C7 and C5 conformers has been proved and energetically as well as spectroscopically characterized. This joint theoretical-experimental investigation demonstrated the feasibility of obtaining accurate structures for flexible small biomolecules, thus paving the route to the elucidation of the inherent behavior of peptides.

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.

Accurate molecular structure and spectroscopic properties of nucleobases: a combined computational-microwave investigation of 2-thiouracil as a case study

The computational composite scheme purposely set up for accurately describing the electronic structure and spectroscopic properties of small biomolecules has been applied to the first study of the rotational spectrum of 2-thiouracil. The experimental investigation was made possible thanks to the combination of the laser ablation technique with Fourier transform microwave spectrometers. The joint experimental-computational study allowed us to determine the accurate molecular structure and spectroscopic properties of the title molecule, but more importantly, it demonstrates a reliable approach for the accurate investigation of isolated small biomolecules.

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.