Accurate structure, thermodynamic and spectroscopic parameters from CC and CC/DFT schemes: the challenge of the conformational equilibrium in glycine

Two-level stochastic search of low-energy conformers for molecular spectroscopy: implementation and validation of MM and QM models

The search for stationary points in the molecular potential energy surfaces (PES) is a problem of increasing relevance in molecular sciences especially for large, flexible systems featuring several large-amplitude internal motions.

Molecular synthons for accurate structural determinations: the equilibrium geometry of 1-chloro-1-fluoroethene

The equilibrium structure for 1-chloro-1-fluoroethene is obtained by high-level quantum chemical calculations and the semi-experimental approach.

Diving for Accurate Structures in the Ocean of Molecular Systems with the Help of Spectroscopy and Quantum Chemistry

The prediction and interpretation of structural properties are the starting points for a deep understanding of thermochemistry, kinetics, and spectroscopic signatures of molecular systems. To give an example, detailed knowledge of the conformational behavior of the main building blocks of biomolecules in the gas phase (i.e., without the perturbing effect of the environment) is a mandatory prerequisite toward the understanding of the role played by different interactions in determining the biological activity in terms of structure-activity relationships. The first step to take is an unambiguous definition of molecular structure. We address the so-called Born-Oppenheimer equilibrium structure, which is defined in a rigorous manner and isotopically independent, and the target accuracy. For the latter, we aim at so-called "spectroscopic" accuracy, which implies uncertainties of a few milliangstroms for bond lengths and smaller than a tenth of degree for angles. If on one side the continuous enhancements of the experimental techniques give access to new and unprecedented spectroscopic determinations, on the other side they require increasing efforts for an unbiased interpretation and analysis. Among the pieces of information, accurate molecular structures play a particularly important role. Indeed, there is a strong relationship between the experimental outcome and the electronic structure of the system. Spectroscopic techniques, in particular those exploited in the gas phase, are therefore accurate and reliable sources for structural information. However, it is seldom straightforward to derive molecular structures directly from the experimental information. Indeed, even in the favorable case of investigations in the gas phase, vibrational effects are always present, and disentangling their contributions requires collection of information for all vibrational modes, a nearly impossible task. To overcome these limitations, joint theory-spectroscopy strategies can be identified, which are referred to as "top-down" and "bottom-up". The first approach, denoted as the semiexperimental approach, relies on extracting from experimental outcomes the equilibrium structure by using quantum-chemical computations to recover vibrational effects. The bottom-up approach consists in verifying the computed equilibrium geometry by means of a comparison between calculated and experimental spectroscopic parameters that probe structural characteristics. In this contribution, we try to review the most important challenges in accurate molecular structure determinations, with particular emphasis on the "solution" provided by a joint theoretical-experimental approach and on the current state of the art. Starting from the illustration of different strategies, we proceed by addressing the increasing complexity in the derivation of equilibrium geometries: we start from the construction of a database of accurate structures, we then face the problem of extending the dimension of the systems amenable to accurate structural determinations, and finally we move to the challenge of understanding the nature of intermolecular interactions.

Correct Modeling of Cisplatin: a Paradigmatic Case

Quantum chemistry is a useful tool in modern approaches to drug and material design, but only when the adopted model reflects a correct physical picture. Paradigmatic is the case of cis‐diaminodichloroplatinum(II), cis‐[Pt(NH₃)₂Cl₂], for which the correct simulation of the structural and vibrational properties measured experimentally still remains an open question. By using this molecule as a proof of concept, it is shown that state‐of‐the‐art quantum chemical calculations and a simple model, capturing the basic physical flavors, a cis‐[Pt(NH₃)₂Cl₂] dimer, can provide the accuracy required for interpretative purposes. The present outcomes have fundamental implications for benchmark studies aiming at assessing the accuracy of a given computational protocol.

VMS-ROT: A New Module of the Virtual Multifrequency Spectrometer for Simulation, Interpretation, and Fitting of Rotational Spectra

The Virtual Multifrequency Spectrometer (VMS) is a tool that aims at integrating a wide range of computational and experimental spectroscopic techniques with the final goal of disclosing the static and dynamic physical-chemical properties "hidden" in molecular spectra. VMS is composed of two parts, namely, VMS-Comp, which provides access to the latest developments in the field of computational spectroscopy, and VMS-Draw, which provides a powerful graphical user interface (GUI) for an intuitive interpretation of theoretical outcomes and a direct comparison to experiment. In the present work, we introduce VMS-ROT, a new module of VMS that has been specifically designed to deal with rotational spectroscopy. This module offers an integrated environment for the analysis of rotational spectra: from the assignment of spectral transitions to the refinement of spectroscopic parameters and the simulation of the spectrum. While bridging theoretical and experimental rotational spectroscopy, VMS-ROT is strongly integrated with quantum-chemical calculations, and it is composed of four independent, yet interacting units: (1) the computational engine for the calculation of the spectroscopic parameters that are employed as a starting point for guiding experiments and for the spectral interpretation, (2) the fitting-prediction engine for the refinement of the molecular parameters on the basis of the assigned transitions and the prediction of the rotational spectrum of the target molecule, (3) the GUI module that offers a powerful set of tools for a vis-à-vis comparison between experimental and simulated spectra, and (4) the new assignment tool for the assignment of experimental transitions in terms of quantum numbers upon comparison with the simulated ones. The implementation and the main features of VMS-ROT are presented, and the software is validated by means of selected test cases ranging from isolated molecules of different sizes to molecular complexes. VMS-ROT therefore offers an integrated environment for the analysis of the rotational spectra, with the innovative perspective of an intimate connection to quantum-chemical calculations that can be exploited at different levels of refinement, as an invaluable support and complement for experimental studies.

Can THF hydrogen bond to glycine as strong as water?

Hydrogen bond complexation between glycine and THF and between glycine and water involving four lowest-energy glycine conformers have been studied. The complexes have been investigated in the gas phase at the ab initio molecular orbital theory (MP2) with aug-cc-pVDZ basis set and density functional theory (B3LYP) with aug-cc-pVTZ basis set. Bader’s theory of atoms in molecules (AIM), natural bond orbital (NBO), and symmetry adapted perturbation theory (SAPT) analyses are employed to elucidate the interaction characteristics in the complexes. The premise that the hydrogen bond donor ability of the O–H group of the carboxyl group dominates the interaction between glycine and THF and between glycine and water is confirmed. It is found that in comparison with water, THF binds more strongly to glycine. The quantum studies indicate that contribution of N–H···O and C–H···O hydrogen bonds in the complexes, although lower in magnitude to O–H···O interactions, play an important role in the stability of complexes. The blue and red shifts in the stretching frequencies of the hydrogen bond donors X–H (X = O, C, N) have also been related to stabilization energies. Decomposition of the stabilization energy based on the SAPT method clearly indicates the dominant role of the electrostatic interactions in all the complexes under study; however, induction and dispersion interaction terms are relatively higher in glycine–THF complexes.

Development and Implementation of Advanced Fitting Methods for the Calculation of Accurate Molecular Structures

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 dynamical and environmental effects in tuning their overall behavior. For this purpose the so-called semiexperimental approach, based on a nonlinear least-squares fit of the moments of inertia associated with a set of available isotopologues, allows one to obtain very accurate results, without the unfavorable computational cost characterizing high-level quantum chemical methods. In the present work the MSR (Molecular Structure Refinement) software for the determination of equilibrium structures by means of the semiexperimental approach is presented, and its implementation is discussed in some detail. The software, which is interfaced with a powerful graphical user interface, includes different optimization algorithms, an extended error analysis, and a number of advanced features, the most remarkable ones concerning the choice of internal coordinates and the method of predicate observations. In particular, a new black-box scheme for defining automatically a suitable set of nonredundant internal coordinates of A₁ symmetry in place of the customary Z-matrix has been designed and tested. Finally, the implementation of the method of the predicate observations is discussed and validated for a set of test molecules. As an original application, the method is employed for the determination of the semiexperimental structure for the most stable conformer of glycine.

Accurate Vibrational-Rotational Parameters and Infrared Intensities of 1-Bromo-1-fluoroethene: A Joint Experimental Analysis and Ab Initio Study

The medium-resolution gas-phase infrared (IR) spectra of 1-bromo-1-fluoroethene (BrFC═CH₂, 1,1-C₂H₂BrF) were investigated in the range 300-6500 cm^-1, and the vibrational analysis led to the assignment of all fundamentals as well as many overtone and combination bands up to three quanta, thus giving an accurate description of its vibrational structure. Integrated band intensity data were determined with high precision from the measurements of their corresponding absorption cross sections. The vibrational analysis was supported by high-level ab initio investigations. CCSD(T) computations accounting for extrapolation to the complete basis set and core correlation effects were employed to accurately determine the molecular structure and harmonic force field. The latter was then coupled to B2PLYP and MP2 computations in order to account for mechanical and electrical anharmonicities. Second-order perturbative vibrational theory was then applied to the thus obtained hybrid force fields to support the experimental assignment of the IR spectra.

Structural features of the carbon-sulfur chemical bond: a semi-experimental perspective

In this work semi-experimental and theoretical equilibrium geometries of 10 sulfur-containing organic molecules, as well as 4 oxygenated ones, are determined by means of a computational protocol based on density functional theory. The results collected in the present paper further enhance our online database of accurate semi-experimental equilibrium molecular geometries, adding 13 new molecules containing up to 8 atoms, for 12 of which the first semi-experimental equilibrium structure is reported, to the best of our knowledge. We focus in particular on sulfur-containing compounds, aiming both to provide new accurate data on some rather important chemical moieties, only marginally represented in the literature of the field, and to examine the structural features of carbon-sulfur bonds in the light of the previously presented linear regression approach. The structural changes issuing from substitution of oxygen by sulfur are discussed to get deeper insights on how modifications in electronic structure and nuclear potential can affect equilibrium geometries. With respect to our previous works, we perform non-linear constrained optimizations of equilibrium SE structures with a new general and user-friendly software under development in our group with updated definition of useful statistical indicators.

State-of-the-Art Thermochemical and Kinetic Computations for Astrochemical Complex Organic Molecules: Formamide Formation in Cold Interstellar Clouds as a Case Study

We describe an integrated computational strategy aimed at providing reliable thermochemical and kinetic information on the formation processes of astrochemical complex organic molecules. The approach involves state-of-the-art quantum-mechanical computations, second-order vibrational perturbation theory, and kinetic models based on capture and transition state theory together with the master equation approach. Notably, tunneling, quantum reflection, and leading anharmonic contributions are accounted for in our model. Formamide has been selected as a case study in view of its interest as a precursor in the abiotic amino acid synthesis. After validation of the level of theory chosen for describing the potential energy surface, we have investigated several pathways of the OH + CH2NH and NH2 + H2CO reaction channels. Our results show that both reaction channels are essentially barrierless (in the sense that all relevant transition states lie below or only marginally above the reactants) and once tunneling is taken into the proper account indicate that the reaction can occur under the low temperature conditions of interstellar environments.

Reliable vibrational wavenumbers for C=O and N-H stretchings of isolated and hydrogen-bonded nucleic acid bases

The accurate prediction of vibrational wavenumbers for functional groups involved in hydrogen-bonded bridges remains an important challenge for computational spectroscopy. For the specific case of the C=O and N-H stretching modes of nucleobases and their oligomers, the paucity of experimental reference values needs to be compensated by reliable computational data, which require the use of approaches going beyond the standard harmonic oscillator model. Test computations performed for model systems (formamide, acetamide and their cyclic homodimers) in the framework of the second order vibrational perturbation theory (VPT2) confirmed that anharmonic corrections can be safely computed by global hybrid (GHF) or double hybrid (DHF) functionals, whereas the harmonic part is particularly challenging. As a matter of fact, GHFs perform quite poorly and even DHFs, while fully satisfactory for C=O stretchings, face unexpected difficulties when dealing with N-H stretchings. On these grounds, a linear regression for N-H stretchings has been obtained and validated for the heterodimers formed by 4-aminopyrimidine with 6-methyl-4-pyrimidinone (4APM-M4PMN) and by uracil with water. In view of the good performance of this computational model, we have built a training set of B2PLYP-D3/maug-cc-pVTZ harmonic wavenumbers (including linear regression scaling for N-H) for six-different uracil dimers and a validation set including 4APM-M4PMN, one of the most stable hydrogen-bonded adenine homodimers, as well as the adenine-uracil, adenine-thymine, guanine-cytosine and adenine-4-thiouracil heterodimers. Because of the unfavourable scaling of DHF harmonic wavenumbers with the dimensions of the investigated systems, we have optimized a linear regression of B3LYP-D3/N07D harmonic wavenumbers for the training set, which has been next checked against the validation set. This relatively cheap model, which shows very good agreement with experimental data (average errors of about 10 cm(-1)), paves the route toward a reliable analysis of spectroscopic signatures for larger polynucleotides.

The Virtual Multifrequency Spectrometer: a new paradigm for spectroscopy

On going developments of hardware and software are changing computational spectroscopy from a strongly specialized research area to a general tool in the inventory of most researchers. Increased interactions between experimentally-oriented users and theoretically-oriented developers of new methods and models would result in more robust, flexible and reliable tools and studies for the systems of increasing complexity, which are of current scientific and technological interest. This is the philosophy behind this review, which presents the development of a so-called virtual multi-frequency spectrometer (VMS) including state-of-the-art approaches in a user-friendly frame. The current status of the VMS tool will be illustrated by a number of case studies with special reference to infrared and UV-vis regions of the electro-magnetic spectrum including also chiral spectroscopies. Only the basic theoretical background will be provided avoiding explicit equations as far as possible, and pointing out the most recent advancements beyond the standard rigid-rotor harmonic-oscillator model coupled to vertical electronic excitation energies.

Anharmonic Effects on Vibrational Spectra Intensities: Infrared, Raman, Vibrational Circular Dichroism, and Raman Optical Activity

The aim of this paper is 2-fold. First, we want to report the extension of our virtual multifrequency spectrometer (VMS) to anharmonic intensities for Raman optical activity (ROA) with the full inclusion of first- and second-order resonances for both frequencies and intensities in the framework of the generalized second-order vibrational perturbation theory (GVPT2) for all kinds of vibrational spectroscopies. Then, from a more general point of view, we want to present and validate the performance of VMS for the parallel analysis of different vibrational spectra for medium-sized molecules (IR, Raman, VCD, ROA) including both mechanical and electric/magnetic anharmonicity. For the well-known methyloxirane benchmark, careful selection of density functional, basis set, and resonance thresholds permitted us to reach qualitative and quantitative agreement between experimental and computed band positions and shapes. Next, the whole series of halogenated azetidinones is analyzed, showing that it is now possible to interpret different spectra in terms of mass, electronegativity, polarizability, and hindrance variation between closely related substituents, chiral spectroscopies being particular effective in this connection.

CC/DFT Route toward Accurate Structures and Spectroscopic Features for Observed and Elusive Conformers of Flexible Molecules: Pyruvic Acid as a Case Study

The structures and relative stabilities as well as the rotational and vibrational spectra of the three low-energy conformers of pyruvic acid (PA) have been characterized using a state-of-the-art quantum-mechanical approach designed for flexible molecules. By making use of the available experimental rotational constants for several isotopologues of the most stable PA conformer, Tc-PA, the semiexperimental equilibrium structure has been derived. The latter provides a reference for the pure theoretical determination of the equilibrium geometries for all conformers, thus confirming for these structures an accuracy of 0.001 Å and 0.1 deg for bond lengths and angles, respectively. Highly accurate relative energies of all conformers (Tc-, Tt-, and Ct-PA) and of the transition states connecting them are provided along with the thermodynamic properties at low and high temperatures, thus leading to conformational enthalpies accurate to 1 kJ mol(-1). Concerning microwave spectroscopy, rotational constants accurate to about 20 MHz are provided for the Tt- and Ct-PA conformers, together with the computed centrifugal-distortion constants and dipole moments required to simulate their rotational spectra. For Ct-PA, vibrational frequencies in the mid-infrared region accurate to 10 cm(-1) are reported along with theoretical estimates for the transitions in the near-infrared range, and the corresponding infrared spectrum including fundamental transitions, overtones, and combination bands has been simulated. In addition to the new data described above, theoretical results for the Tc- and Tt-PA conformers are compared with all available experimental data to further confirm the accuracy of the hybrid coupled-cluster/density functional theory (CC/DFT) protocol applied in the present study. Finally, we discuss in detail the accuracy of computational models fully based on double-hybrid DFT functionals (mainly at the B2PLYP/aug-cc-pVTZ level) that avoid the use of very expensive CC calculations.

Semiexperimental equilibrium structures for building blocks of organic and biological molecules: the B2PLYP route

The B2PLYP double hybrid functional, coupled with the correlation-consistent triple-ζ cc-pVTZ (VTZ) basis set, has been validated in the framework of the semiexperimental (SE) approach for deriving accurate equilibrium structures of molecules containing up to 15 atoms. A systematic comparison between new B2PLYP/VTZ results and several equilibrium SE structures previously determined at other levels, in particular B3LYP/SNSD and CCSD(T) with various basis sets, has put in evidence the accuracy and the remarkable stability of such model chemistry for both equilibrium structures and vibrational corrections. New SE equilibrium structures for phenylacetylene, pyruvic acid, peroxyformic acid, and phenyl radical are discussed and compared with literature data. Particular attention has been devoted to the discussion of systems for which lack of sufficient experimental data prevents a complete SE determination. In order to obtain an accurate equilibrium SE structure for these situations, the so-called templating molecule approach is discussed and generalized with respect to our previous work. Important applications are those involving biological building blocks, like uracil and thiouracil. In addition, for more general situations the linear regression approach has been proposed and validated.

Quantum Chemistry Meets Spectroscopy for Astrochemistry: Increasing Complexity toward Prebiotic Molecules

For many years, scientists suspected that the interstellar medium was too hostile for organic species and that only a few simple molecules could be formed under such extreme conditions. However, the detection of approximately 180 molecules in interstellar or circumstellar environments in recent decades has changed this view dramatically. A rich chemistry has emerged, and relatively complex molecules such as C60 and C70 are formed. Recently, researchers have also detected complex organic and potentially prebiotic molecules, such as amino acids, in meteorites and in other space environments. Those discoveries have further stimulated the debate on the origin of the building blocks of life in the universe. Many efforts continue to focus on the physical, chemical, and astrophysical processes by which prebiotic molecules can be formed in the interstellar dust and dispersed to Earth or to other planets.Spectroscopic techniques, which are widely used to infer information about molecular structure and dynamics, play a crucial role in the investigation of planetary atmosphere and the interstellar medium. Increasingly these astrochemical investigations are assisted by quantum-mechanical calculations of structures as well as spectroscopic and thermodynamic properties, such as transition frequencies and reaction enthalpies, to guide and support observations, line assignments, and data analysis in these new and chemically complicated situations. However, it has proved challenging to extend accurate quantum-chemical computational approaches to larger systems because of the unfavorable scaling with the number of degrees of freedom (both electronic and nuclear).In this Account, we show that it is now possible to compute physicochemical properties of building blocks of biomolecules with an accuracy rivaling that of the most sophisticated experimental techniques, and we summarize specific contributions from our groups. As a test case, we present the underlying computational machinery through the investigation of oxirane. We describe how we determine the molecular structure and then how we characterize the rotational and IR spectra, the most important issues for a correct theoretical description and a proper comparison with experiment. Next, we analyze the spectroscopic properties of representative building blocks of DNA bases (uracil and pyrimidine) and of proteins (glycine and glycine dipeptide analogue).Solvation, surface chemistry (dust fraction, adsorption, desorption), and inter- and intramolecular interactions, such as self-organization and self-interaction, are important molecular processes for understanding astrochemistry. Using the specific cases of uracil dimers and glycine adsorbed on silicon grains, we also illustrate approaches in which we treat different regions, interactions, or effects at different levels of sophistication.

Criteria for first- and second-order vibrational resonances and correct evaluation of the Darling-Dennison resonance coefficients using the canonical Van Vleck perturbation theory

The second-order vibrational Hamiltonian of a semi-rigid polyatomic molecule when resonances are present can be reduced to a quasi-diagonal form using second-order vibrational perturbation theory. Obtaining exact vibrational energy levels requires subsequent numerical diagonalization of the Hamiltonian matrix including the first- and second-order resonance coupling coefficients. While the first-order Fermi resonance constants can be easily calculated, the evaluation of the second-order Darling-Dennison constants requires more complicated algebra for seven individual cases with different numbers of creation-annihilation vibrational quanta. The difficulty in precise evaluation of the Darling-Dennison coefficients is associated with the previously unrecognized interference with simultaneously present Fermi resonances that affect the form of the canonically transformed Hamiltonian. For the first time, we have presented the correct form of the general expression for the evaluation of the Darling-Dennison constants that accounts for the underlying effect of Fermi resonances. The physically meaningful criteria for selecting both Fermi and Darling-Dennison resonances are discussed and illustrated using numerical examples.

Quantum chemical study on HKrC₅N, HXeC₅N, and related rare gas compounds

The recent identification of HRgC5N (Rg = Kr, Xe) in a cryogenic matrix calls for an in-depth theoretical study on these compounds. Here we present the results of CCSD(T), MP2, and DFT calculations concerning the molecular structure, stability, and vibrational spectroscopy. The procedure combining CCSD(T) calculations for variable H-Rg distances with the anharmonic description of the corresponding stretching vibration, based on a Morse-type potential energy function, was proposed and has led to good agreement between computational and experimental values for H-Rg stretching frequencies, at relatively low computational costs. High Raman scattering activity of HRgC5N and of its isomers, predicted at the DFT level, gives some prospects for the detection of these molecules with a method alternative to the IR absorption spectroscopy.

Conformational preferences of Ac-Gly-NHMe in solution

The conformational behaviour of Ac-Gly-NHMe and its fluorinated [CF₃-C(O)-Gly-NHMe] and N-methyl[Ac-Gly-N(Me)₂] derivatives is investigated in nonpolar, polar and polar protic solutions by NMR and IR spectroscopies and theoretical calculations.

Microsolvation of 2-thiouracil: molecular structure and spectroscopic parameters of the thiouracil-water complex

State-of-the-art quantum-chemical computations have been employed to accurately determine the equilibrium structure and interaction energy of the 2-thiouracil-water complex, thus extending available reference data for biomolecule solvation patterns. The coupled-cluster level of theory in conjunction with a triple-ζ basis set has been considered together with extrapolation to the basis set limit, performed by employing second-order Møller-Plesset perturbation theory, and inclusion of core-correlation and diffuse-function corrections. On the basis of the comparison of experiment and theory for 2-thiouracil [ Puzzarini et al. Phys. Chem. Chem. Phys. 2013 , 15 , 16965 - 16975 ], structural changes due to water complexation have been pointed out. Molecular and spectroscopic properties of the 2-thiouracil-water complex have then been studied by means of the composite computational approach introduced for the molecular structure evaluation. Among the results achieved, we mention the accurate determination of the molecular dipole moment and of the spectroscopic parameters required for predicting the rotational spectrum.

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.

A new Gaussian MCTDH program: implementation and validation on the levels of the water and glycine molecules

We report the main features of a new general implementation of the Gaussian Multi-Configuration Time-Dependent Hartree model. The code allows effective computations of time-dependent phenomena, including calculation of vibronic spectra (in one or more electronic states), relative state populations, etc. Moreover, by expressing the Dirac-Frenkel variational principle in terms of an effective Hamiltonian, we are able to provide a new reliable estimate of the representation error. After validating the code on simple one-dimensional systems, we analyze the harmonic and anharmonic vibrational spectra of water and glycine showing that reliable and converged energy levels can be obtained with reasonable computing resources. The data obtained on water and glycine are compared with results of previous calculations using the vibrational second-order perturbation theory method. Additional features and perspectives are also shortly discussed.

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.

Dispersion corrected DFT approaches for anharmonic vibrational frequency calculations: nucleobases and their dimers

Computational spectroscopy techniques have become in the last few years an effective means to analyze and assign infrared (IR) spectra of molecular systems of increasing dimensions and in different environments. However, transition from compilation of harmonic data to fully anharmonic simulations of spectra is still underway. The most promising results for large systems have been obtained, in our opinion, by perturbative vibrational approaches based on potential energy surfaces computed by hybrid (especially B3LYP) density functionals and medium size (e.g. SNSD) basis sets. In this framework, we are actively developing a comprehensive and robust computational protocol aimed at quantitative reproduction of the spectra of nucleic acid base complexes and their adsorption on solid supports (organic/inorganic). In this contribution we report the essential results of the first step devoted to isolated monomers and dimers. It is well known that in order to model the vibrational spectra of weakly bound molecular complexes dispersion interactions should be taken into proper account. In this work we have chosen two popular and inexpensive approaches to model dispersion interactions, namely the semi-empirical dispersion correction (D3) and pseudopotential based (DCP) methodologies both in conjunction with the B3LYP functional. These have been used for simulating fully anharmonic IR spectra of nucleobases and their dimers through generalized second order vibrational perturbation theory (GVPT2). We have studied, in particular, isolated adenine, hypoxanthine, uracil, thymine and cytosine, the hydrogen-bonded and stacked adenine and uracil dimers, and the stacked adenine-naphthalene heterodimer. Anharmonic frequencies are compared with standard B3LYP results and experimental findings, while the computed interaction energies and structures of complexes are compared to the best available theoretical estimates.

Fully anharmonic IR and Raman spectra of medium-size molecular systems: accuracy and interpretation

Computation of full infrared (IR) and Raman spectra (including absolute intensities and transition energies) for medium- and large-sized molecular systems beyond the harmonic approximation is one of the most interesting challenges of contemporary computational chemistry. Contrary to common beliefs, low-order perturbation theory is able to deliver results of high accuracy (actually often better than those issuing from current direct dynamics approaches) provided that anharmonic resonances are properly managed. This perspective sketches the recent developments in our research group toward the development of a robust and user-friendly virtual spectrometer rooted in second-order vibrational perturbation theory (VPT2) and usable also by non-specialists essentially as a black-box procedure. Several examples are explicitly worked out in order to illustrate the features of our computational tool together with the most important ongoing developments.

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.

An integrated experimental and quantum-chemical investigation on the vibrational spectra of chlorofluoromethane

The vibrational analysis of the gas-phase infrared spectra of chlorofluoromethane (CH2ClF, HCFC-31) was carried out in the range 200-6200 cm(-1). The assignment of the absorption features in terms of fundamental, overtone, combination, and hot bands was performed on the medium-resolution (up to 0.2 cm(-1)) Fourier transform infrared spectra. From the absorption cross section spectra accurate values of the integrated band intensities were derived and the global warming potential of this compound was estimated, thus obtaining values of 323, 83, and 42 on a 20-, 100-, and 500-year horizon, respectively. The set of spectroscopic parameters here presented provides the basic data to model the atmospheric behavior of this greenhouse gas. In addition, the obtained vibrational properties were used to benchmark the predictions of state-of-the-art quantum-chemical computational strategies. Extrapolated complete basis set limit values for the equilibrium geometry and harmonic force field were obtained at the coupled-cluster singles and doubles level of theory augmented by a perturbative treatment of triple excitations, CCSD(T), in conjunction with a hierarchical series of correlation-consistent basis sets (cc-pVnZ, with n = T, Q, and 5), taking also into account the core-valence correlation effects and the corrections due to diffuse (aug) functions. To obtain the cubic and quartic semi-diagonal force constants, calculations employing second-order Møller-Plesset perturbation (MP2) theory, the double-hybrid density functional B2PLYP as well as CCSD(T) were performed. For all anharmonic force fields the performances of two different perturbative approaches in computing the vibrational energy levels (i.e., the generalized second order vibrational treatment, GVPT2, and the recently proposed hybrid degeneracy corrected model, HDCPT2) were evaluated and the obtained results allowed us to validate the spectroscopic predictions yielded by the HDCPT2 approach. The predictions of the deperturbed second-order perturbation approach, DVPT2, applied to the computation of infrared intensities beyond the double-harmonic approximation were compared to the accurate experimental values here determined. Anharmonic DFT and MP2 corrections to CCSD(T) intensities led to a very good agreement with the absorption cross section measurements over the whole spectral range here analysed.

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.

Anharmonic theoretical simulations of infrared spectra of halogenated organic compounds

The recent implementation of the computation of infrared (IR) intensities beyond the double-harmonic approximation [J. Bloino and V. Barone, J. Chem. Phys. 136, 124108 (2012)] paved the route to routine calculations of infrared spectra for a wide set of molecular systems. Halogenated organic compounds represent an interesting class of molecules, from both an atmospheric and computational point of view, due to the peculiar chemical features related to the halogen atoms. In this work, we simulate the IR spectra of eight halogenated molecules (CH2F2, CHBrF2, CH2DBr, CF3Br, CH2CHF, CF2CFCl, cis-CHFCHBr, cis-CHFCHI), using two common hybrid and double-hybrid density functionals in conjunction with both double- and triple-ζ quality basis sets (SNSD and cc-pVTZ) as well as employing the coupled-cluster theory with basis sets of at least triple-ζ quality. Finally, we compare our results with available experimental spectra, with the aim of checking the accuracy and the performances of the computational approaches.

Accurate structure, thermodynamics, and spectroscopy of medium-sized radicals by hybrid coupled cluster/density functional theory approaches: the case of phenyl radical

The coupled-cluster singles doubles model with perturbative treatment of triples (CCSD(T)) coupled with extrapolation to the complete basis-set limit and additive approaches represent the "golden standard" for the structural and spectroscopic characterization of building blocks of biomolecules and nanosystems. However, when open-shell systems are considered, additional problems related to both specific computational difficulties and the need of obtaining spin-dependent properties appear. In this contribution, we present a comprehensive study of the molecular structure and spectroscopic (IR, Raman, EPR) properties of the phenyl radical with the aim of validating an accurate computational protocol able to deal with conjugated open-shell species. We succeeded in obtaining reliable and accurate results, thus confirming and, partly, extending the available experimental data. The main issue to be pointed out is the need of going beyond the CCSD(T) level by including a full treatment of triple excitations in order to fulfil the accuracy requirements. On the other hand, the reliability of density functional theory in properly treating open-shell systems has been further confirmed.