A second-order perturbation theory route to vibrational averages and transition properties of molecules: general formulation and application to infrared and vibrational circular dichroism spectroscopies

Hydrogen-atom-assisted thione-thiol tautomerization of thiourea derivatives in para-H2 matrix

Thiourea (TU) and its N-methylated derivative, N-methyl thiourea (NMTU), were exposed to H atoms generated in cryogenic para-H2 matrices. The reactions were followed online by FT-IR spectroscopy. The freshly deposited matrices exclusively contained the more stable thione tautomeric forms. However, upon exposure to H atoms, the peaks belonging to the precursor molecules clearly decreased along with the simultaneous appearance of new signals. These new bands could be attributed to the corresponding higher-energy thiol forms (in the case of TU) and, tentatively, to an intermediate radical (in both the TU and NMTU experiments). The radicals are suggested to be the H-atom-addition products of the TU and NMTU thione precursors, with the addition occurring on the S atom. These intermediates may then react with another free H atom, leading to the formation of the more energetic thiol tautomers, following an H-atom-abstraction process. As such, these radicals act as the centerpiece of the reaction scheme, enabling the thione-thiol tautomerization. This H-atom-assisted process is similar to that observed for the related molecule, thioacetamide. The interpretation of the experimental results was supplemented by quantum-chemical computations, which predicted all the above-mentioned reactions to be barrierless. The presence of H atoms opens a barrierless pathway; thus, the process does not necessarily require activation through irradiation (e.g., broadband UV). These findings point to the ubiquitous nature of the facile hydrogenation/dehydrogenation of the S atom, implying that thione-thiol tautomerization may occur easily.

Gas-phase, conformer-specific infrared spectra of 3-chlorophenol and 3-fluorophenol

Conformational isomerism of phenol derivatives has been a subject of extensive spectroscopic study. Combining the capabilities of the widely tuneable infrared free-electron laser FELIX with molecular beam technologies allows for revisiting existing data and gaining additional insights into far-IR spectroscopy of halogenated phenols. Here we present conformer-resolved infrared spectra of the syn and anti conformers of 3-chlorophenol and 3-fluorophenol recorded via IR-UV ion-dip spectroscopy. The experimental work is complemented by density functional theory calculations to aid in assignment of the observed bands. The experimental spectra for the two conformers of each molecule show overall a great similarity, but also include some distinct conformer-specific bands in the spectral range investigated. Our spectra confirm previously reported OH torsional mode frequencies for the syn and anti conformers of 3-chlorophenol (3CP) at 315 cm-1, (Manocha et al., J. Phys. Chem., 1973, 77, 2094) but reverse their assignment of the 311 and 319 cm-1 bands for 3-fluorophenol. 1D torsional mode calculations were performed for 3CP to help assign possible OH torsion overtones.

Accelerating Quantum Anharmonic Vibrational Calculations by Atom-Specific Hybrid Basis Set-Based Potential Energy Surface Approach

The development of accurate yet fast quantum mechanical methods to calculate the anharmonic vibrational spectra of large molecules is one of the major goals of ongoing developments in this field. This study extensively explores and validates a hybrid electronic basis set approach for anharmonic vibrational calculations, where the molecule is segregated into different computational layers, and such layers are then treated with different levels of electronic basis sets. Following the system-bath model, the atoms corresponding to the active sites are treated in more accurate but computationally slower, large basis set and the rest of the atoms in less accurate but computationally faster, small basis set to construct the anharmonic hybrid potential energy surface (PES). Such a hybrid protocol for constructing an anharmonic PES is named as the atom-specific hybrid basis set (ASHBS) approach. The accuracy of the ASHBS approach is tested and established by evaluating the harmonic and anharmonic frequencies of a set of four prototype molecules. Following the ASHBS approach, for a chosen active site, the transitions of the corresponding modes are found to be closer to that of a high basis set, with a majority of target modes displaying a mean absolute error of around 3.3 cm-1, while achieving the computational acceleration of 2-3 times. This study also provides insights into determining the optimal layer size to balance computational efficiency and accuracy, offering a suitable alternative, particularly for large molecules.

Constrained Nuclear-Electronic Orbital Transition State Theory Using Energy Surfaces with Nuclear Quantum Effects

Hydrogen-atom transfer is crucial in a myriad of chemical and biological processes, yet the accurate and efficient description of hydrogen-atom transfer reactions and kinetic isotope effects remains challenging due to significant quantum effects on hydrogenic motion, especially tunneling and zero-point energy. In this paper, we combine transition state theory (TST) with the recently developed constrained nuclear-electronic orbital (CNEO) theory to propose a new transition state theory denoted CNEO-TST. We use CNEO-TST with CNEO density functional theory (CNEO-DFT) to predict reaction rate constants for two prototypical gas-phase hydrogen-atom transfer reactions and their deuterated isotopologic reactions. CNEO-TST is similar to conventional TST except that it employs constrained minimized energy surfaces to include zero-point energy and shallow tunneling effects in the effective potential. We find that the new theory predicts reaction rates quite accurately at room temperature. The effective potential surface must be generated by CNEO theory rather than by ordinary electronic structure theory, but because of the favorable computational scaling of CNEO-DFT, the cost is economical even for large systems. Our results show that dynamics calculations with this approach achieve accuracy comparable to variational TST with a semiclassical multidimensional tunneling transmission coefficient at and above room temperature. Therefore, CNEO-TST can be a useful tool for rate prediction, even for reactions involving highly quantal motion, such as many chemical and biochemical reactions involving transfers of hydrogen atoms, protons, or hydride ions.

Phenyl Radical Activates Molecular Hydrogen Through Protium and Deuterium Tunneling

Activating dihydrogen, H2, is a challenging endeavor typically achieved using transition metal centers. Pure main-group compounds capable of this are rare and have emerged in recent decades. These systems rely on synergistic donor-acceptor interactions with H2's antibonding σ* and bonding σ orbital. An alternative (hydrocarbon) radical-mediated activation is problematic because the H-H bond is stronger (104.2 kcal mol-1) than most C-H bonds. Here, we explore using the phenyl radical to tackle this, forming benzene with a C-H bond energy (112.9 kcal mol-1) that provides a thermodynamic driving force. We mainly observe a benzene-HI complex upon photolysis of iodobenzene in an H2-doped neon matrix at 4.4 K despite a barrier of 7.6 kcal mol-1, while phenyl radical forms in case of the heavier D2 isotopologue. When D2 molecules are allowed to diffuse, mono-deuterated benzene accumulates within hours. Computations using path integral-based instanton theory highlight that primarily the transferred hydrogen atom is moving during the reaction which greatly increases the tunneling probability. In excellent agreement with the experimental results, we predict significant tunneling rate constants for both isotopologues, H2 and D2, featuring a strong kinetic isotope effect of up to four orders of magnitude at the lowest temperatures.

Accurate Structure and Spectroscopic Properties of Azulene and Its Derivatives by Means of Pisa Composite Schemes and Vibrational Perturbation Theory to Second Order

The structural and spectroscopic properties in the gas phase of azulene and some of its N-bearing derivatives have been analyzed by a general computational strategy based on the recent Pisa composite schemes (PCSs). First of all, an accurate semiexperimental equilibrium structure has been derived for azulene and employed to validate the geometrical parameters delivered by different quantum chemical methods. Next, different isomerization energies (azulene to naphthalene, 1-aza-azulene to quinoline and to other isomers) have been computed by an explicitly correlated PCS version employing frozen natural orbitals. Accurate geometries have been obtained by a cheaper PCS variant based on a double-hybrid functional improved by one-parameter bond corrections, with the same functional providing also remarkable harmonic frequencies. The corresponding equilibrium rotational constants show average deviations within 0.1% from experimental results when taking into account anharmonic vibrational corrections obtained by a global hybrid functional. Therefore, reliable computational estimates have been produced for the rotational constants of several nitrogen derivatives (isomeric aza-azulenes and guaiazulene), whose non-negligible dipole moments could allow experimental microwave characterizations. An analogous approach delivers infrared spectra in remarkable agreement with their experimental counterparts for naphthalene, quinoline, and azulene, together with reliable predictions for the still-unknown spectrum of 1-aza-azulene. In addition to their intrinsic interest, the results of this paper further confirm that a very accurate yet robust and user-friendly tool is now available for aiding high-resolution spectroscopic studies of quite large systems of current technological and/or biological interest.

Assessing the Partial Hessian Approximation in QM/MM-Based Vibrational Analysis

The partial Hessian approximation is often used in vibrational analysis of quantum mechanics/molecular mechanics (QM/MM) systems because calculating the full Hessian matrix is computationally impractical. This approach aligns with the core concept of QM/MM, which focuses on the QM subsystem. Thus, using the partial Hessian approximation implies that the main interest is in the local vibrational modes of the QM subsystem. Here, we investigate the accuracy and applicability of the partial Hessian vibrational analysis (PHVA) approach as it is typically used within QM/MM, i.e., only the Hessian belonging to the QM subsystem is computed. We focus on solute-solvent systems with small, rigid solutes. To separate two of the major sources of errors, we perform two separate analyses. First, we study the effects of the partial Hessian approximation on local normal modes, harmonic frequencies, and harmonic IR and Raman intensities by comparing them to those obtained using full Hessians, where both partial and full Hessians are calculated at the QM level. Then, we quantify the errors introduced by QM/MM used with the PHVA by comparing normal modes, frequencies, and intensities obtained using partial Hessians calculated using a QM/MM-type embedding approach to those obtained using partial Hessians calculated at the QM level. Another aspect of the PHVA is the appearance of normal modes resembling the translation and rotation of the QM subsystem. These pseudotranslational and pseudorotational modes should be removed as they are collective vibrations of the atoms in the QM subsystem relative to a frozen MM subsystem and, thus, not well-described. We show that projecting out translation and rotation, usually done for systems in isolation, can adversely affect other normal modes. Instead, the pseudotranslational and pseudorotational modes can be identified and removed.

Observation of metastable structures of the ethylene glycol-water dimer in helium nanodroplets

Ethylene glycol (EG) is the simplest organic diol. Here we measure infrared spectra of the EG monomer and its dimer with water, the complex, EG(H2O), embedded in superfluid helium nanodroplets. For the monomer, only a single, gauche, conformation is observed. For EG(H2O), no trace of the global energy minimum is seen, a structure that would maximize the hydrogen bonding contacts. Instead, only metastable structures are formed, suggesting that dimerization in a superfluid environment leads to kinetic trapping in local energy minima. In addition, we obtain evidence for a dimer where the conformation of EG switches from gauche to trans on account of dimerization with a water molecule. This observation is assumed to be driven over an energy barrier by utilizing the energy released as hydrogen bonding occurs.

Vibrational Features of Oxyamines: A Comparative Study of N,N-Diethylhydroxylamine and N,N-Diethylacetyloxyamine

The ability of oxyamines to undergo homolytic cleavage of the O-C bond, leading to the formation of stable radicals, is widely used in polymerization processes and to prevent oxidative stress in materials. We present a mid and near-infrared spectroscopy study on two model compounds, N,N-diethylhydroxyloxyamine (C4H11NO) and its acetyl derivative N,N-diethylacetyloxyamine (C6H13NO2) in the liquid phase. The analysis of the spectra is based on a complete exploration of the conformational space, coupled to harmonic and anharmonic calculations performed using the generalized second-order vibrational perturbation theory (GVPT2) formalism at the B3LYP-D3(BJ)/Def2-TZVP level of calculation and potential energy distribution analysis. In the most stable species out of 25, the three amine chains present an all-anti arrangement, with the carbonyl oxygen atom pointing towards the nitrogen lone pair. The simulated spectra are in overall good agreement with the experimental ones, and suitable for the assignment of the main observed bands. Furthermore, similarities and divergences between the two molecules are discussed.

Toward the identification of cyano-astroCOMs via vibrational features: benzonitrile as a test case

The James Webb Space Telescope (JWST) opened a new era for the identification of molecular systems in the interstellar medium (ISM) by vibrational features. One group of molecules of increasing interest is cyano-derivatives of aromatic organic molecules, which have already been identified in the ISM on the basis of the analysis of rotational signatures, and so, are plausible candidates for the detection by the JWST. Benzonitrile considered in this work represents a suitable example for the validation of a computational strategy, which can be further applied for different, larger, and not-yet observed molecules. For this purpose, anharmonic simulations of infrared (IR) spectra have been compared with recent FTIR experimental studies. The anharmonic computations using the generalized second-order vibrational perturbation theory (GVPT2) in conjunction with a hybrid force field combining the harmonic part of revDSD-PBEP86-D3/jun-cc-pVTZ with anharmonic corrections from B3LYP-D3/SNSD show very good agreement with those in the experiment, with a mean error of 11 c m - 1 for all fundamental transitions overall and only 2 c m - 1 for the C ≡ N stretching fundamental at 4.49 μ m . The inclusion of overtones up to three-quanta transitions also allowed the prediction of spectra in the near-infrared region, which shows distinct features due to C ≡ N overtones at the 2.26 μ m and 1.52 μ m . The remarkable accuracy of the GVPT2 results opens a pathway for the reliable prediction of spectra for a broader range of cyano-astroCOMs.

Accurate Vibrational and Ro-Vibrational Contributions to the Properties of Large Molecules by a New Engine Employing Curvilinear Internal Coordinates and Vibrational Perturbation Theory to Second Order

The unbiased comparison between theory and experiment requires approaches more sophisticated than the basic harmonic-oscillator rigid-rotor model, for taking into account vibrational averaging effects and ro-vibrational couplings in molecules of increasing size. Second-order vibrational perturbation theory based on curvilinear internal coordinates (ICs) offers a remarkable compromise between accuracy and computational cost, thanks to the reduction of mode-mode couplings with respect to their counterparts based on Cartesian coordinates. Therefore, we have developed, implemented, and validated a general engine employing ICs, which allows the accurate evaluation of vibrational averages and ro-vibrational couplings for molecules containing up to about 50 atoms beyond the harmonic approximation. After validation of the new tool for relatively small molecules, the effectiveness of ICs has been demonstrated for some flexible and/or quite large molecular bricks of life.

After 50 Years of Vibrational Circular Dichroism Spectroscopy: Challenges and Opportunities of Increasingly Accurate and Complex Experiments and Computations

VCD research continues to thrive, driven by ongoing experimental and theoretical advances. Modern studies deal with increasingly complex samples featuring weak intermolecular interactions and shallow potential energy surfaces. Likewise, the combination of VCD measurements with, for instance, cryo-spectroscopic techniques has significantly increased their sensitivity. The extent to which such modern measurements enhance the informative value of VCD depends significantly on the quality of the theoretical models, which must adequately account for anharmonicity, solvation and molecular dynamics. We herein discuss how experimental advancements engage in a stimulating interplay with recent theoretical developments, pursuing either the static or the dynamic computational route. Both paths have their own strengths and limitations, each addressing fundamentally different problems. We give an outlook on future challenges of VCD research, including the possibility to combine static and dynamic approaches to obtain a full picture of the sample.

Hydrogen-atom-assisted processes on thioacetamide in para-H2 matrix - formation of thiol tautomers

TA has been isolated in low-temperature para-H2 matrices and it has been exposed to H atoms. In accordance with previous experimental results, TA exclusively exists in its more stable thione tautomeric form in the freshly deposited matrix. However, upon H atom generation, the bands belonging to the precursor start decreasing with the simultaneous appearance of new bands. By comparing the position of these new peaks with earlier results obtained in para-H2, they can be undoubtedly attributed to the presence of the higher-energy thiol tautomeric forms of TA. No other products could be observed, except for NH3. Quantum-chemical computations have been invoked to understand the mechanism behind the observed thione-thiol tautomerization assisted by H atoms. Accordingly, tautomerization starts with barrierless H-atom addition on the S atom of the thione resulting in the formation of the intermediate 1-amino-1-mercapto-ethyl radical (H3C-Ċ(-SH)-NH2), which has been detected tentatively for the first time. The second step is a barrierless H-atom induced H-abstraction from the -NH2 moiety of the H3C-Ċ(-SH)-NH2 radical. A comprehensive mechanism for tautomerization is proposed based on the experimental and theoretical results. Although earlier studies showed the possibility of TA thione-thiol tautomerization in cryogenic matrices achieved by energetic UV irradiation, the present study points out that this can also take place through a barrierless pathway by simply exposing the TA molecules to H atoms. As such, this is the first evidence for the occurrence of such a reaction in a matrix-isolated environment. The current results also help us better understand the mystery behind thione-thiol tautomerization in low-temperature environments.

Harmonic and anharmonic vibrational computations for biomolecular building blocks: Benchmarking DFT and basis sets by theoretical and experimental IR spectrum of glycine conformers

Advanced vibrational spectroscopic experiments have reached a level of sophistication that can only be matched by numerical simulations in order to provide an unequivocal analysis, a crucial step to understand the structure-function relationship of biomolecules. While density functional theory (DFT) has become the standard method when targeting medium-size or larger systems, the problem of its reliability and accuracy are well-known and have been abundantly documented. To establish a reliable computational protocol, especially when accuracy is critical, a tailored benchmark is usually required. This is generally done over a short list of known candidates, with the basis set often fixed a priori. In this work, we present a systematic study of the performance of DFT-based hybrid and double-hybrid functionals in the prediction of vibrational energies and infrared intensities at the harmonic level and beyond, considering anharmonic effects through vibrational perturbation theory at the second order. The study is performed for the six-lowest energy glycine conformers, utilizing available "state-of-the-art" accurate theoretical and experimental data as reference. Focusing on the most intense fundamental vibrations in the mid-infrared range of glycine conformers, the role of the basis sets is also investigated considering the balance between computational cost and accuracy. Targeting larger systems, a broad range of hybrid schemes with different computational costs is also tested.

Importance of Polarizable Embedding for Computing Optical Rotation: The Case of Camphor in Ethanol

Solvation effects on optical rotation are notoriously challenging to model for computational chemistry, as the specific rotatory power of a molecule can vary wildly going from apolar to polar or even protic solvents. To address such a problem, we present a polarizable embedding implementation of an electric and magnetic response property based on density functional theory and the AMOEBA polarizable force field, and apply such an implementation to the study of the optical rotation of camphor in ethanol. By comparing a continuum model, and electrostatic and polarizable embedding QM/MM models, we observe that accounting for the environment's polarization gives rise to not only a different quantitative prediction, in very good agreement with experiments for the QM/AMOEBA model, but also to a very different qualitative picture, with the values of the optical rotation computed along a classical molecular dynamics trajectory with electrostatic embedding being statistically uncorrelated to the ones obtained with the polarizable description.

Photochemical Generation and Characterization of C-Aminophenyl-Nitrilimines: Insights on Their Bond-Shift Isomers by Matrix-Isolation IR Spectroscopy and Density Functional Theory Calculations

The intriguing ability of C-phenyl-nitrilimine to co-exist as allenic and propargylic bond-shift isomers motivated us to investigate how substituents in the phenyl ring influence this behavior. Building on our previous work on the meta- and para-OH substitution, here we extended this investigation to explore the effect of the NH2 substitution. For this purpose, C-(4-aminophenyl)- and C-(3-aminophenyl)-nitrilimines were photogenerated in an argon matrix at 15 K by narrowband UV-light irradiation (λ = 230 nm) of 5-(4-aminophenyl)- and 5-(3-aminophenyl)-tetrazole, respectively. The produced nitrilimines were further photoisomerized to carbodiimides via 1H-diazirines by irradiations at longer wavelengths (λ = 380 or 330 nm). Combining IR spectroscopic measurements and DFT calculations, it was found that the para-NH2-substituted nitrilimine exists as a single isomeric structure with a predominant allenic character. In contrast, the meta-NH2-substituted nitrilimine co-exists as two bond-shift isomers characterized by propargylic and allenic structures. To gain further understanding of the effects of phenyl substitution on the bond-shift isomerism of the nitrilimine fragment, we compared geometric and charge distribution data derived from theoretical calculations performed for C-phenyl-nitrilimine with those performed for the derivatives resulting from NH2 (electron-donating group) and NO2 (electron-withdrawing group) phenyl substitutions.

Performance of Effective Harmonic Oscillator Approach for the Calculations of Vibrational Transition Energies of Large Molecules

The accuracy and performance of the effective harmonic oscillator approximation for the description of anharmonic vibrational structure calculations are tested for large molecular systems and compared with experimental values along with vibrational self-consistent field and second-order perturbation theories. The effective harmonic oscillator approach is an effective single-particle approximation where the variational parameters are the centroids and widths of the multidimensional Gaussian product functions posited as the vibrational wave functions. A comprehensive calculation for 849 transitions that include the fundamentals, two and three quanta overtone transitions, and several combination bands of three polyaromatic hydrocarbons and one DNA nucleobase with a total of 231 normal modes are assessed. A comparison of EHO results with the experimental values is done for the polyaromatic hydrocarbons, and a close agreement is found between the two results. It also offers anharmonic eigenstates and eigenfunctions that are nearly identical with vibrational self-consistent field theory. An extensive analysis on the resultant wave functions of the excited states is performed. The overall root-mean-square deviation (RMSD) between these two methods for 849 transitions understudy is only about 8.3 cm-1, suggesting the effective harmonic oscillator as a viable alternative for the reliable calculations of transition energies of large molecular systems.

Vibrational Circular Dichroism Spectroscopy with a Classical Polarizable Force Field: Alanine in the Gas and Condensed Phases

Polarizable force fields are an essential component for the chemically accurate modeling of complex molecular systems with a significant degree of fluxionality, beyond harmonic or perturbative approximations. In this contribution we examine the performance of such an approach for the vibrational spectroscopy of the alanine amino acid, in the gas and condensed phases, from the Fourier transform of appropriate time correlation functions generated along molecular dynamics (MD) trajectories. While the infrared (IR) spectrum only requires the electric dipole moment, the vibrational circular dichroism (VCD) spectrum further requires knowledge of the magnetic dipole moment, for which we provide relevant expressions to be used with polarizable force fields. The AMOEBA force field was employed here to model alanine in the neutral and zwitterionic isolated forms, solvated by water or nitrogen, and as a crystal. Within this framework, comparison of the electric and magnetic dipole moments to those obtained with nuclear velocity perturbation theory based on density-functional theory for the same MD trajectories are found to agree well with one another. The statistical convergence of the IR and VCD spectra is examined and found to be more demanding in the latter case. Comparisons with experimental frequencies are also provided for the condensed phases.

Scaling-up VPT2: A feasible route to include anharmonic correction on large molecules

Vibrational analysis plays a crucial role in the investigation of molecular systems. Though methodologies like second-order vibrational perturbation theory (VPT2) have paved the way to more accurate simulations, the computational cost remains a difficult barrier to overcome when the molecular size increases. Building upon recent advances in the identification of resonances, we propose an approach making anharmonic simulations possible for large-size systems, typically unreachable by standard means. This relies on the fact that, often, only portions of the whole spectra are of actual interest. Therefore, the anharmonic corrections can be included selectively on subsets of normal modes directly related to the regions of interest. Starting from the VPT2 equations, we evaluate rigorously and systematically the impact of the truncated anharmonic treatment onto simulations. The limit and feasibility of the reduced-dimensionality approach are detailed, starting on a smaller model system. The methodology is then challenged on the IR absorption and vibrational circular dichroism spectra of an organometallic complex in three different spectral ranges.

Multiple-time scale integration method based on an interpolated potential energy surface for ab initio path integral molecular dynamics

A new multiple-time scale integration method is presented that propagates ab initio path integral molecular dynamics (PIMD). This method uses a large time step to generate an approximate geometrical configuration whose energy and gradient are evaluated at the level of an ab initio method, and then, a more precise integration scheme, e.g., the Bulirsch-Stoer method or velocity Verlet integration with a smaller time step, is used to integrate from the previous step using the computationally efficient interpolated potential energy surface constructed from two consecutive points. This method makes the integration of PIMD more efficient and accurate compared with the velocity Verlet integration. A Nosé-Hoover chain thermostat combined with this new multiple-time scale method has good energy conservation even with a large time step, which is usually challenging in velocity Verlet integration for PIMD due to the very small chain mass when a large number of beads are used. The new method is used to calculate infrared spectra and free energy profiles to demonstrate its accuracy and capabilities.

Performance of EOM-CCSD(T)(a)*-Based Quartic Force Fields in Computing Fundamental, Anharmonic Vibrational Frequencies of Molecular Electronically Excited States with Application to the Ã1A″ State of :CCH2 (Vinylidene)

Highly accurate anharmonic vibrational frequencies of electronically excited states are not as easily computed as their ground electronic state counterparts, but recently developed approximate triple excited state methods may be changing that. One emerging excited state method is equation of motion coupled cluster theory at the singles and doubles level with perturbative triples computed via the (a)* formalism, EOMEE-CCSD(T)(a)*. One of the most employed means for the ready computation of vibrational anharmonic frequencies for ground electronic states is second-order vibrational perturbation theory (VPT2), a theory based on quartic force fields (QFFs),fourth-order Taylor series expansions of the potential portion of the internuclear Watson Hamiltonian. The combination of these two is herein benchmarked for its performance for use as a means of computing rovibrational spectra of electronically excited states. Specifically, the EOMEE-CCSD(T)(a)* approach employing a complete basis set extrapolation along with core electron inclusion and relativity (the so-called "CcCR" approach) defining the QFF produces anharmonic fundamental vibrational frequencies within 2.83%, on the average, of reported gas-phase experimentally assigned values for the test set including the A ~ 1 A ″ states of HCF, HCCl, HSiF, HNO, and HPO. However, some states have exceptional accuracy in the fundamentals, most notably for ν2 of A ~ 1 A ″ HCCl in which the CcCR QFF value is within 1.8 cm-1 at 927.9 cm-1 (or 0.2%) of the experiment. Additionally, this approach produces rotational constants to, on the absolute average, within 0.41% of available experimental data, showcasing notable accuracy in the computation of rovibronic spectral data. Furthermore, utilizing a hybrid approach composed of harmonic CcCR force constants along with a set of simple EOMEE-CCSD(T)(a)*/aug-cc-pVQZ QFF cubic and quartic force constants is faster than using pure CcCR and better represents those modes that suffer from numerical instability in the anharmonic portion of the QFF, implying that this so-called "CcCR + QZ" QFF approach may be the best for future applications. Finally, complete, rovibrational spectral data are provided for A ~ 1 A 2 :CCH2, a molecule of potential astrochemical interest, in order to aid in its potential future experimental rovibronic characterization.

Methoxyacetone revisited

Conformational space of methoxyacetone (MA) was studied at the MP2/6-311++G(d,p) and DFT(B3LYP)/6-311++G(d,p) levels of theory. Computations predict MA to adopt four conformations, resulting from internal rotations around the O=C-C-O (Trans, Cis) and C-C-O-C (trans, gauche) dihedral angles. The Tt (Trans-trans) conformer is the most stable. The computed energies of two gauche (Tg and Cg) conformers fall in the 3-8 kJ mol-1 range above Tt and should account for 1/3 of the room-temperature gas-phase equilibrium. The energy of Ct form is 11 kJ mol-1 above Tt, and its expected population is negligible (below 1 %). In our earlier work, MA monomers were isolated in cryogenic argon matrices and characterized by infrared spectroscopy. In the experiment, only the most stable Tt conformer was detected in the sample. Signatures of the other conformers were not detected, either in freshly deposited samples, or in samples subjected to different UV irradiations. We rationalize those observations in terms of computed barriers for intramolecular torsions, indicating occurrence of conformational cooling during deposition. The experimental infrared spectrum of the Tt form is now assigned with the aid of anharmonic DFT computations. Exposure of MA to UV irradiation in the 300-260 nm range led to photolysis, according to the Norrish type II mechanism, resulting in dimer between enol acetone and formaldehyde observed as a cage-confined intermediate photoproduct. The subsequent photolysis resulted in the formation of carbon monoxide as the dominating photoproduct, formed in the Norrish type I photoreaction. Mechanistic interpretation of this photo decarbonylation reaction is presented.

Infrared Spectroscopy of Fluorenyl Cations at Cryogenic Temperatures

The notion of (anti)aromaticity is a successful concept in chemistry to explain the structure and stability of polycyclic hydrocarbons. Cyclopentadienyl and fluorenyl cations are among the most studied classical antiaromatic systems. In this work, fluorenyl cations are investigated by high-resolution gas-phase infrared spectroscopy in helium droplets. Bare fluorenyl cations are generated in the gas phase by electrospray ionization. After mass-to-charge selection, ions are captured in ultracold helium nanodroplets and probed by infrared spectroscopy using a widely tunable free-electron laser in the 600-1700 cm-1 range. The highly resolved cryogenic infrared spectra confirm, in combination with DFT computations, that all cations are present in their singlet states.

Reassignment of the Photoelectron Spectrum of Methylketene Using a Hybrid Model of Harmonic and Anharmonic Oscillators to Compute Franck-Condon Factors

We constructed a hybrid model of harmonic and anharmonic oscillators to compute Franck-Condon factors and interpret the photoelectron spectrum of methylketene. The equilibrium structures of methylketene and its cation were optimized, and then, the harmonic and anharmonic vibrational frequencies were computed using the B3LYP, PBE0, APFD, and ωB97XD approaches of the density functional theory. The photoelectron spectrum of methylketene was simulated by computing the Franck-Condon factors with both the harmonic and hybrid models. The adiabatic ionization energy of methylketene was computed by using the CCSD(T) approach extrapolating to the complete basis set limit. The simulated photoelectron spectra are consistent with those from the experiment for both the harmonic and hybrid models. However, the error in band positions is reduced by using the hybrid model. The computed adiabatic ionization energies of methylketene are in agreement with the experiment, with the smallest error being 0.017 eV. Our interpretation based on the theoretical spectrum led to the reassignment of the experimental photoelectron spectrum of methylketene.

Comparison of different density functional theory methods for the calculation of vibrational circular dichroism spectra

The determination of the absolute configuration (AC) of an organic molecule is still a challenging task for which the combination of spectroscopic with quantum-mechanical methods has become a promising approach. In this study, we investigated the accuracy of DFT methods (480 overall combinations of 15 functionals, 16 basis sets, and 2 solvation models) to calculate the VCD spectra of six chiral organic molecules in order to benchmark their capability to facilitate the determination of the AC.

Free Ethynylarsinidene and Ethynylstibinidene: Heavier Analogues of Nitrenes and Phosphinidenes

Until now, there has been very little experimental evidence for the existence of free arsinidenes and stibinidenes, apart from the hydrides, AsH and SbH. Here, we report on photogeneration of triplet ethynylarsinidene, HCCAs, and triplet ethynylstibinidene, HCCSb, from ethynylarsine and ethynylstibine, respectively, in solid argon matrices. The products were identified using infrared spectroscopy and the associated UV absorption spectra are interpreted with the aid of theoretical predictions.

Intramolecular hole-transfer in protonated anthracene

Excitation spectra of protonated and deuteronated anthracene are obtained by triple-resonance dissociation spectroscopy. Very cold cations, protonated/deuteronated exclusively at the 9-position, are generated from two-colour two-photon threshold ionisation of 9-dihydroanthracenyl radicals (C14H11). The excitation spectra reveal rich structure, not resolved in previous studies, that is assigned based on anharmonic and Herzberg-Teller coupling calculations. This work reveals that the excitation of protonated anthracene induces a symmetry-breaking intramolecular charge-transfer process along a Marcus-Hush coordinate, where the positively charged hole hops from the central bridging sp2 carbon, onto one of the aromatic rings. Signatures of this charge-transfer event are observed in the excitation spectrum, through active Herzberg-Teller progressions.

Computational infrared and Raman spectra by hybrid QM/MM techniques: a study on molecular and catalytic material systems

Vibrational spectroscopy is one of the most well-established and important techniques for characterizing chemical systems. To aid the interpretation of experimental infrared and Raman spectra, we report on recent theoretical developments in the ChemShell computational chemistry environment for modelling vibrational signatures. The hybrid quantum mechanical and molecular mechanical approach is employed, using density functional theory for the electronic structure calculations and classical forcefields for the environment. Computational vibrational intensities at chemical active sites are reported using electrostatic and fully polarizable embedding environments to achieve more realistic vibrational signatures for materials and molecular systems, including solvated molecules, proteins, zeolites and metal oxide surfaces, providing useful insight into the effect of the chemical environment on the signatures obtained from experiment. This work has been enabled by the efficient task-farming parallelism implemented in ChemShell for high-performance computing platforms.  This article is part of a discussion meeting issue 'Supercomputing simulations of advanced materials'.

pbqff: Push-Button Quartic Force Fields

pbqff is an open-source program for fully automating the production of quartic force fields (QFFs) and their corresponding anharmonic spectroscopic data. Rather than being a monolithic piece of code, it consists of several key modules including a generic interface to quantum chemistry codes and, notably, queuing systems; a molecular point group symmetry library; an internal-to-Cartesian coordinate conversion module; a module for the ordinary least-squares fitting of potential energy surfaces; and an improved second-order rotational and vibrational perturbation theory package for asymmetric and symmetric tops that handles type-1 and -2 Fermi resonances, Fermi resonance polyads, and Coriolis resonances. All of these pieces are written in Rust, a modern, safe, and performant programming language that has much to offer for scientific programming. This work introduces pbqff and its surrounding ecosystem, in addition to reporting new anharmonic vibrational data for c-(C)C₃H₂ and describing how the components of pbqff can be leveraged in other projects.

Infrared spectra of isoquinolinium (iso-C9H7NH+) and isoquinolinyl radicals (iso-C9H7NH and 1-, 3-, 4-, 5-, 6-, 7- and 8-iso-HC9H7N) isolated in solid para-hydrogen

Protonated polycyclic aromatic nitrogen heterocycles (H⁺PANH) are prospective candidates that may contribute to interstellar unidentified infrared (UIR) emission bands because protonation enhances the relative intensities of the bands near 6.2, 7.7 and 8.6 μm, and the presence of the N atom induces a blue shift of the ring-stretching modes so that the spectra of H⁺PANH match better with the 6.2 μm feature in class-A UIR spectra. We report the infrared (IR) spectra of protonated isoquinoline (the 2-isoquinolinium cation, iso-C₉H₇NH⁺), its neutral counterpart (the 2-isoquinolinyl radical, iso-C₉H₇NH), and another mono-hydrogenated product (the 6-isoquinolinyl radical, 6-iso-HC₉H₇N), produced on the electron-bombardment of a mixture of isoquinoline (iso-C₉H₇N) with excess para-hydrogen (p-H₂) during matrix deposition at 3.2 K. To generate additional isomers of hydrogenated isoquinoline, we irradiated iso-C₉H₇N/Cl₂/p-H₂ matrices at 365 nm to generate Cl atoms, followed by IR irradiation to generate H atoms via Cl + H₂ (v = 1) → HCl + H; the H atoms thus generated reacted with iso-C₉H₇N. In addition to iso-C₉H₇NH and 6-iso-HC₉H₇N observed in the electron-bombardment experiments, we identified six additional hydrogenated isoquinoline species, 1-, 3-, 4-, 5-, 7- and 8-iso-HC₉H₇N, via their IR spectra; hydrogenation on the N atom and all available carbon atoms except for the two sharing carbon atoms on the fused ring was observed. Spectral groupings were achieved according to their behaviors after maintenance of the matrix in darkness and on secondary photolysis at various wavelengths. The assignments were supported via comparison of the experimental results with the vibrational wavenumbers and IR intensities of possible isomers predicted using the B3LYP/6-311++G(d,p) method. The implications in the identification of the UIR band are discussed.

Dinitrogen Activation in the Gas Phase: Spectroscopic Characterization of C-N Coupling in the V3 C+ +N2 Reaction

We report on cluster-mediated C-N bond formation in the gas phase using N₂ as a nitrogen source. The V₃ C⁺ +N₂ reaction is studied by a combination of ion-trap mass spectrometry with infrared photodissociation (IRPD) spectroscopy and complemented by electronic structure calculations. The proposed reaction mechanism is spectroscopically validated by identifying the structures of the reactant and product ions. V₃ C⁺ exhibits a pyramidal structure of C₁ -symmetry. N₂ activation is initiated by adsorption in an end-on fashion at a vanadium site, followed by spontaneous cleavage of the N≡N triple bond and subsequent C-N coupling. The IRPD spectrum of the metal nitride product [NV₃ (C=N)]⁺ exhibits characteristic C=N double bond (1530 cm^-1 ) and V-N single bond (770, 541 and 522 cm^-1 ) stretching bands.

Vibrational resonance analysis of linear molecules using resummation of divergent Rayleigh-Schrödinger perturbation theory series

The large order Rayleigh-Schrödinger perturbation theory (RSPT) was applied for calculating vibrational states of linear molecules. Two molecules (CO₂ and C₂H₂) were used as test cases with using of isomorphic Watson Hamiltonian and quartic force fields. For CO₂ the Sayvetz condition can remove all degeneracies for purely vibrational states and the non-degenerate perturbation theory can be applied. However, an existence of two degenerate modes in C₂H₂ requires using the upgraded degenerate version of RSPT that was employed in this context for the first time. The dominating divergent behavior of such series requires the resummation technique that mimics the multivalued nature of the underlying solutions, and the applied quartic Padé-Hermite approximants (QPHA) provided full solution of the problem. Moreover, some mathematical properties of QPHA proved to be an efficient tool for studying resonance effects through the Katz theorem that controls the singular points of the eigenvalues on the complex plane. In the case of C₂H₂, not only all earlier observed classical resonances were confirmed and quantified, but also subtle interpolyad resonances (K_2/55,K_3/4555), proposed recently by Herman (2011) were described as well. Following the analysis, we found several novel resonances, of which we proposed one independent interpolyad resonance K_2/4444. The complete analysis of such critical points provided the full resonance picture of all studied molecules.

Dual Photochemistry of Benzimidazole

Monomers of benzimidazole trapped in an argon matrix at 15 K were characterized by vibrational spectroscopy and identified as 1H-tautomers exclusively. The photochemistry of matrix-isolated 1H-benzimidazole was induced by excitations with a frequency-tunable narrowband UV light and followed spectroscopically. Hitherto unobserved photoproducts were identified as 4H- and 6H-tautomers. Simultaneously, a family of photoproducts bearing the isocyano moiety was identified. Thereby, the photochemistry of benzimidazole was hypothesized to follow two reaction pathways: the fixed-ring and the ring-opening isomerizations. The former reaction channel results in the cleavage of the NH bond and formation of a benzimidazolyl radical and an H-atom. The latter reaction channel involves the cleavage of the five-membered ring and concomitant shift of the H-atom from the CH bond of the imidazole moiety to the neighboring NH group, leading to 2-isocyanoaniline and subsequently to the isocyanoanilinyl radical. The mechanistic analysis of the observed photochemistry suggests that detached H-atoms, in both cases, recombine with the benzimidazolyl or isocyanoanilinyl radicals, predominantly at the positions with the largest spin density (revealed using the natural bond analysis computations). The photochemistry of benzimidazole therefore occupies an intermediate position between the earlier studied prototype cases of indole and benzoxazole, which exhibit exclusively the fixed-ring and the ring-opening photochemistries, respectively.

Matrix-isolation and cryosolution-VCD spectra of α-pinene as benchmark for anharmonic vibrational spectra calculations

The inclusion of anharmonicity in vibrational spectral analyis remains associated to small molecular systems with up to a dozen of atoms, with half a dozen of non-hydrogen atoms, typically thesize of propylene oxide. One may see two reasons for this: first of all, larger systems are often thought to be computationally too demanding (high computational costs) for a full anharmonic vibrational analysis. Second, the identification of resonances and their correction is often considered something only expert theoreticians could address because of the lack of unequivocal criteria. In this contribution, we illustrate that resonances can indeed become a complex problem, which can be handled almost transparently thanks to recent advances in vibrational perturbation theory (VPT2). The study also emphasizes the importance and the central role played by experiment in benchmarking novel theoretical approaches. In fact, we herein provide the currently highest resolution VCD spectra available for α- and β-pinene obtained under matrix-isolation conditions and in liquid Xenon as solvent. They are interpreted by VPT2 with novel tests for the identification of resonances. Hence, the study demonstrates the mutual stimulation of advances in both experimental techniques and computational models.

Impact of anion polarizability on ion pairing in microhydrated salt clusters

Despite longstanding interest in the mechanism of salt dissolution in aqueous media, a molecular level understanding remains incomplete. Here, cryogenic ion trap vibrational action spectroscopy is combined with electronic structure calculations to track salt hydration in a gas phase model system one water molecule at a time. The infrared photodissociation spectra of microhydrated lithium dihalide anions [LiXX'(H2O) n ]- (XX' = I2, ClI and Cl2; n = 1-3) in the OH stretching region (3800-2800 cm-1) provide a detailed picture of how anion polarizability influences the competition among ion-ion, ion-water and water-water interactions. While exclusively contact ion pairs are observed for n = 1, the formation of solvent-shared ion pairs, identified by markedly red-shifted OH stretching bands (<3200 cm-1), originating from the bridging water molecules, is favored already for n = 2. For n = 3, Li+ reaches its maximum coordination number of four only in [LiI2(H2O)3]-, in accordance with the hard and soft Lewis acid and base principle. Water-water hydrogen bond formation leads to a different solvent-shared ion pair motif in [LiI2(H2O)3]- and network formation even restabilizes the contact ion pair motif in [LiCl2(H2O)3]-. Structural assignments are exclusively possible after the consideration of anharmonic effects. Molecular dynamics simulations confirm that the significance of large amplitude motion (of the water molecules) increases with increasing anion polarizability and that needs to be considered already at cryogenic temperatures.

Spectroscopic and quantum mechanical study of a scavenger molecule: N,N-diethylhydroxylamine

We report a combination of quantum mechanical calculations and a range of spectroscopic measurements in the gas phase of N,N-diethylhydroxylamine, an important scavenger compound. Three conformers were observed by pulsed jet Fourier transform microwave spectroscopy in the 6.5-18.5 GHz frequency range. They are characterized by the hydroxyl hydrogen atom being in trans orientation with respect to the bisector of the CNC angle while the side alkyl chains can be both trans (global minimum, C_s symmetry, A = 7608.1078(4), B = 2020.2988(2) and C = 1760.5423(2) MHz) or one trans and the other gauche (second energy minimum, A = 5302.896(1), B = 2395.9822(4) and C = 1804.8567(3) MHz) or gauche' (third energy minimum, A = 5960.8025(6), B = 2273.6627(4) and C = 1975.8074(4) MHz). For the global minimum, the ¹³C_α,¹³C_β and ¹⁵N isotopologues were observed in natural abundance, allowing for an accurate partial structure determination. Moreover, several lines were detected by free jet absorption millimeter wave spectroscopy in the 59.6-74.4 GHz spectral range. The electron binding energies of the highest occupied molecular orbital and the next-to-highest occupied molecular orbital, determined by photoelectron spectroscopy, are 8.95 and 10.76 eV, respectively. Supporting calculations evidence that, (i) upon ionization of the HOMO, the molecular structure changes from an amine to an N-oxoammonium arrangement and (ii) the 0-0 of the HOMO-1 photoionization is 10.46 eV. The K-shell binding energies, determined by X-ray photoelectron spectroscopy, are 290.42 eV (C_β), 291.45 eV (C_α), 405.98 eV (N) and 538.75 eV (O). The Fourier transform near infrared spectrum is reported and a tentative assignment is proposed. The equilibrium wavenumber (ω̃ = 3811 cm^-1) and the anharmonicity constant (ω̃χ = -87.5 cm^-1) of the hydroxyl stretching mode were estimated using a quadratic model.

Anharmonic Aspects in Vibrational Circular Dichroism Spectra from 900 to 9000 cm-1 for Methyloxirane and Methylthiirane

Vibrational circular dichroism (VCD) spectra and the corresponding IR spectra of the chiral isomers of methyloxirane and of methylthiirane have been reinvestigated, both experimentally and theoretically, with particular attention to accounting for anharmonic corrections, as calculated by the GVPT2 approach. De novo recorded VCD spectra in the near IR (NIR) range regarding CH-stretching overtone transitions, together with the corresponding NIR absorption spectra, were also considered and accounted for, both with the GVPT2 and with the local mode approaches. Comparison of the two methods has permitted us to better describe the nature of active “anharmonic” modes in the two molecules and the role of mechanical and electrical anharmonicity in determining the intensities of VCD and IR/NIR data. Finally, two nonstandard IR/NIR regions have been investigated: the first one about ≈2000 cm^–1, involving mostly two-quanta bending mode transitions, the second one between 7000 and 7500 cm^–1 involving three-quanta transitions containing CH-stretching overtones and HCC/HCH bending modes.

Molecules of life: studying the interaction between water and phosphine in argon matrices

The interaction between water and phosphine isolated in solid argon matrices has been investigated. Infrared spectra of matrices containing water and phosphine, as well as their deuterium and ¹⁸O isotopologues, revealed that a weak hydrogen bonded complex is formed, with an O-H⋯P bonding motif. This was supported by high level ab initio calculations, which predict a binding energy of only 5.1 kJ mol^-1 (430 cm^-1).

A Computational Protocol for Vibrational Circular Dichroism Spectra of Cyclic Oligopeptides

Cyclic peptides are a promising class of compounds for next-generation antibiotics as they may provide new ways of limiting antibiotic resistance development. Although their cyclic structure will introduce some rigidity, their conformational space is large and they usually have multiple chiral centers that give rise to a wide range of possible stereoisomers. Chiroptical spectroscopies such as vibrational circular dichroism (VCD) are used to assign stereochemistry and discriminate enantiomers of chiral molecules, often in combination with electronic structure methods. The reliable determination of the absolute configuration of cyclic peptides will require robust computational methods than can identify all significant conformers and their relative population and reliably assign their stereochemistry from their chiroptical spectra by comparison with ab initio calculated spectra. We here present a computational protocol for the accurate calculation of the VCD spectra of a series of flexible cyclic oligopeptides. The protocol builds on the Conformer-Rotamer Ensemble Sampling Tool (CREST) developed by Grimme and co-workers ( Phys. Chem. Chem. Phys. 2020, 22, 7169-7192 and J. Chem. Theory. Comput. 2019, 15, 2847-2862) in combination with postoptimizations using B3LYP and moderately sized basis sets. Our recommended computational protocol for the computation of VCD spectra of cyclic oligopeptides consists of three steps: (1) conformational sampling with CREST and tight-binding density functional theory (xTB); (2) energy ranking based on single-point energy calculations as well as geometry optimization and VCD calculations of conformers that are within 2.5 kcal/mol of the most stable conformer using B3LYP/6-31+G*/CPCM; and (3) VCD spectra generation based on Boltzmann weighting with Gibbs free energies. Our protocol provides a feasible basis for generating VCD spectra also for larger cyclic peptides of biological/pharmaceutical interest and can thus be used to investigate promising compounds for next-generation antibiotics.

Matrix Isolation Study of Fumaric and Maleic Acids in Solid Nitrogen

Fumaric and maleic acids ((E)- and (Z)-HOOC-CH═CH-COOH, FA and MA) were studied experimentally by infrared spectroscopy in nitrogen matrixes and theoretically by quantum chemical calculations. The calculations, carried out at the DFT(B3LYP) and MP2 levels of theory, predicted the existence of at least 5 conformers of maleic acid and 10 conformers of fumaric acid. After the deposition of the matrixes, two conformers of maleic acid (I and II) and three conformers (I-III) of fumaric acid were observed and characterized vibrationally. Selective narrowband near-infrared (NIR) excitation of the first OH stretching overtones of the different conformers of maleic and fumaric acids initially present in the matrixes allowed the generation of higher-energy forms, never before observed experimentally. In the case of maleic acid, conformers I (a cis-trans form, where cis and trans designate the conformation of the carboxylic groups of the molecule) and II (cis-cis) were found to generate the novel conformers VI (trans-trans) and VII (cis-trans), respectively. The conversion of conformer II into the most stable conformer I was also observed. For fumaric acid, the cis-cis conformers I-III were found to give rise to the new cis-trans conformers IV-VII, respectively. The tunneling decay of the new conformers produced upon NIR excitation of the lowest-energy conformers initially trapped in the matrixes was observed, and their lifetimes in solid N₂ were determined. The increased stability of all of the observed high-energy conformers of the studied acids in the N₂ matrix, compared to the argon matrix, where they could not be observed experimentally, demonstrates the stabilizing effect of the interaction between the OH groups of the acids with the matrix N₂ molecules, in line with previous observations for other carboxylic acids. In addition, the photochemistry of matrix-isolated maleic and fumaric acids upon broad-band UV irradiation (λ > 235 nm) was also investigated. UV-induced isomerization of both acids around the C═C double bond was observed, together with their decarboxylation to acrylic acid.

Studying the Key Intermediate of RNA Autohydrolysis by Cryogenic Gas-Phase Infrared Spectroscopy

Over the course of the COVID-19 pandemic, mRNA-based vaccines have gained tremendous importance. The development and analysis of modified RNA molecules benefit from advanced mass spectrometry and require sufficient understanding of fragmentation processes. Analogous to the degradation of RNA in solution by autohydrolysis, backbone cleavage of RNA strands was equally observed in the gas phase; however, the fragmentation mechanism remained elusive. In this work, autohydrolysis-like intermediates were generated from isolated RNA dinucleotides in the gas phase and investigated using cryogenic infrared spectroscopy in helium nanodroplets. Data from both experiment and density functional theory provide evidence for the formation of a five-membered cyclic phosphate intermediate and rule out linear or six-membered structures. Furthermore, the experiments show that another prominent condensed-phase reaction of RNA nucleotides can be induced in the gas phase: the tautomerization of cytosine. Both observed reactions are therefore highly universal and intrinsic properties of the investigated molecules.

F12-TZ-cCR: A Methodology for Faster and Still Highly Accurate Quartic Force Fields

The F12-TZ-cCR quartic force field (QFF) methodology, defined here as CCSD(T)-F12b/cc-pCVTZ-F12 with further corrections for relativity, is introduced as a cheaper and even more accurate alternative to more costly composite QFF methods like those containing complete basis set extrapolations within canonical coupled cluster theory. F12-TZ-cCR QFFs produce B₀ and C₀ vibrationally averaged principal rotational constants within 7.5 MHz of gas-phase experimental values for tetraatomic and larger molecules, offering higher accuracy in these constants than the previous composite methods. In addition, F12-TZ-cCR offers an order of magnitude decrease in the computational cost of highly accurate QFF methodologies accompanying this increase in accuracy. An additional order of magnitude in cost reduction is achieved in the F12-DZ-cCR method, while also matching the accuracy of the traditional composite method's B₀ and C₀ constants. Finally, F12-DZ and F12-TZ are benchmarked on the same test set, revealing that both methods can provide anharmonic vibrational frequencies that are comparable in accuracy to all three of the more expensive methodologies, although their rotational constants lag behind. Hence, the present work demonstrates that highly accurate theoretical rovibrational spectral data can be obtained for a fraction of the cost of conventional QFF methodologies, extending the applicability of QFFs to larger molecules.

Infrared spectra of carbocations and CH4+ in helium

Infrared (IR) spectra of several hydrocarbon cations are reported, namely CH₃⁺, CH₄⁺, CH₅⁺, CH₅⁺(CH₄) and C₂H₅⁺. The spectra were generated from weakly-bound helium-cation complexes formed by electron ionization of helium nanodroplets doped with a neutral hydrocarbon precursor. Spectroscopic transitions were registered by photoexcitation of the complexes coupled with mass spectrometric detection of the bare ions. For CH₃⁺, we provide evidence showing that the helium-bound complexes contain 10-20 helium atoms (on average) and have a rotational temperature of ∼5 K. We show that this technique is well-suited to the study of highly symmetric or fluxional ionic species, as these intrinsic properties are preserved in the helium environment. This is in contrast to conventional tagging methods that use a single atom or molecule, which can change the point group or rigidity of the core ion and therefore the spectral profile. We demonstrate this for the highly fluxional molecular ion CH₅⁺, whose spectrum in the current study matches that of the gas phase ion, whereas the fluxionality is lost when a methane tag is added. Finally, we present the first IR spectrum of methane cation, CH₄⁺. The spectrum of this fundamental organic ion shows CH stretching bands consistent with a non-tetrahedral structure, a consequence of Jahn-Teller distortion.

Structural and Energetic Properties of Amino Acids and Peptides Benchmarked by Accurate Theoretical and Experimental Data

Structural, energetic, and spectroscopic data derived in this work aim at the setup of an "experimentally validated" database for amino acids and polypeptides conformers. First, the "cheap" composite scheme (ChS, CCSD(T)/(CBS+CV)^MP2) is tested for evaluation of conformational energies of all eight stable conformers of glycine, by comparing to the more accurate CCSD(T)/CBS+CV computations (Phys. Chem. Chem. Phys. 2013, 15, 10094-10111 and J Mol. Model. 2020, 26, 129). The recently proposed jun-ChS (J. Chem. Theory and Comput. 2020, 16, 988-1006), employing the jun-cc-pVnZ basis set family for CCSD(T) computations and CBS extrapolation, yields conformational energies accurate to 0.2 kJ·mol^-1, at reduced computational cost with respect to aug-ChS employing aug-cc-pVnZ basis sets. The jun-ChS composite scheme is further applied to derive conformational energies for three dipeptide analogues Ac-Gly-NH₂, Ac-Ala-NH₂, and Gly-Gly. Finally, dipeptide conformational energies and semiexperimental equilibrium rotational constants along with the CCSD(T)/(CBS+CV)^MP2 structural parameters (J. Phys. Chem. Lett. 2014, 5, 534-540) stand as the reference for benchmarking of selected density functional methodologies. The double-hybrid functionals B2-PLYP-D3(BJ) and DSD-PBEP86, perform best for structural and energetic characterization of all dipeptide analogues. From hybrid functionals CAM-B3LYP-D3(BJ) and ωB97X-D3(BJ) represent promising methods applicable for larger peptide-based systems for which computations with double-hybrid functionals are not feasible.