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

Investigating the stability of aromatic carboxylic acids in hydrated magnesium sulfate under UV irradiation to assist detection of organics on Mars

The Scanning Habitable Environments with Raman and Luminescence for Organics and Chemicals (SHERLOC) instrument onboard the Mars 2020 Perseverance rover detected so far some of the most intense fluorescence signals in association with sulfates analyzing abraded patches of rocks at Jezero crater, Mars. To assess the plausibility of an organic origin of these signals, it is key to understand if organics can survive exposure to ambient Martian UV after exposure by the Perseverance abrasion tool and prior to analysis by SHERLOC. In this work, we investigated the stability of organo-sulfate assemblages under Martian-like UV irradiation and we observed that the spectroscopic features of phthalic and mellitic acid embedded into hydrated magnesium sulfate do not change for UV exposures corresponding to at least 48 Martian sols and, thus, should still be detectable in fluorescence when the SHERLOC analysis takes place, thanks to the photoprotective properties of magnesium sulfate. In addition, different photoproduct bands diagnostic of the parent carboxylic acid molecules could be observed. The photoprotective behavior of hydrated magnesium sulfate corroborates the hypothesis that sulfates might have played a key role in the preservation of organics on Mars, and that the fluorescence signals detected by SHERLOC in association with sulfates could potentially arise from organic compounds.

Accurate structures and rotational constants of bicyclic monoterpenes at DFT cost by means of the bond-corrected Pisa composite scheme (BPCS)

The structural, conformational, and spectroscopic properties in the gas phase of 20 bicyclic monoterpenes and monoterpenoids have been analyzed by a new accurate, reduced-cost computational strategy. In detail, the revDSD-PBEP86 double-hybrid functional in conjunction with the D3BJ empirical dispersion corrections and a suitable triple-zeta basis set provides accurate geometrical parameters, whence equilibrium rotational constants, which are further improved by proper account of core-valence correlation. Average deviations within 0.1% between computed and experimental rotational constants are reached when taking into account the vibrational corrections obtained by the B3LYP functional in conjunction with a double-zeta basis set in the framework of second-order vibrational perturbation theory. In addition to their intrinsic interest, the studied terpenes further extend the panel of systems for which the proposed strategy has provided accurate results at density functional theory cost. Therefore, a very accurate yet robust and user-friendly tool is now available for systematic investigations of the role of stereo-electronic effects on the properties of large systems of current technological and/or biological interest by experimentally oriented researchers.

Investigation of the Proton-Bound Dimer of Dihydrogen Phosphate and Formate Using Infrared Spectroscopy in Helium Droplets

Understanding the structural and dynamic properties of proton-bound complexes is crucial for elucidating fundamental aspects of chemical reactivity and molecular interactions. In this work, the proton-bound complex between dihydrogen phosphate and formate, and its deuterated counterparts, is investigated using IR action spectroscopy in helium droplets. Contrary to the initial expectation that the stronger phosphoric acid would donate a proton to formate, both experiment and theory show that all exchangeable protons are located in the phosphate moiety. The experimental spectra show good agreement with both scaled harmonic and VPT2 anharmonic calculations, indicating that anharmonic effects are small. Some H-bending modes of the nondeuterated complex are found to be sensitive to the helium environment. In the case of the partially deuterated complexes, the experiments indicate that internal dynamics leads to isomeric interconversion upon IR excitation.

Undecanoic Acid and L-Phenylalanine in Vermiculite: Detection, Characterization, and UV Degradation Studies for Biosignature Identification on Mars

Solar radiation that arrives on the surface of Mars interacts with organic molecules present in the soil. The radiation can degrade or transform the organic matter and make the search for biosignatures on the planet's surface difficult. Therefore, samples to be analyzed by instruments on board Mars probes for molecular content should be selectively chosen to have the highest organic preservation content. To support the identification of organic molecules on Mars, the behavior under UV irradiation of two organic compounds, undecanoic acid and L-phenylalanine, in the presence of vermiculite and two chloride salts, NaCl and MgCl, was studied. The degradation of the molecule's bands was monitored through IR spectroscopy. Our results show that, while vermiculite acts as a photoprotective mineral with L-phenylalanine, it catalyzes the photodegradation of undecanoic acid molecules. On the other hand, both chloride salts studied decreased the degradation of both organic species acting as photoprotectors. While these results do not allow us to conclude on the preservation capabilities of vermiculite, they show that places where chloride salts are present could be good candidates for in situ analytic experiments on Mars due to their organic preservation capacity under UV radiation.

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.

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.

Mid-IR and CH stretching vibrational circular dichroism spectroscopy to distinguish various sources of chirality: The case of quinophaneoxazoline based ruthenium(II) complexes

Five diastereomers of ruthenium(II) complexes based on quinolinophaneoxazoline ligands were investigated by vibrational circular dichroism (VCD) in the mid-IR and CH stretching regions. Diastereomers differ in three sources of chirality: the planar chirality of the quinolinophane moiety, the central chirality of an asymmetric carbon atom of the oxazoline ring, and the chirality of the ruthenium atom. VCD, allied to DFT calculations, has been found to be effective in disentangling the various forms of chirality. In particular, a VCD band is identified in the CH stretching region directly connected to the chirality of the metal. The analysis of the calculated VCD spectra is carried out by partitioning the complexes into fragments. The anharmonic analysis is also performed with a recently proposed reduced-dimensionality approach: such treatment is particularly important when examining spectroscopic regions highly perturbed by resonances, like the CH stretching region.

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.

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.

How Crystal Symmetry Dictates Non-Local Vibrational Circular Dichroism in the Solid State

Solid-State Vibrational Circular Dichroism (VCD) can be used to determine the absolute structure of chiral crystals, but its interpretation remains a challenge in modern spectroscopy. In this work, we investigate the effect of a twofold screw axis on the solid-state VCD spectrum in a combined experimental and theoretical analysis of P2₁ crystals of (S)-(+)-1-indanol. Even though the space group is achiral, a single proper symmetry operation has an important impact on the VCD spectrum, which reflects the supramolecular chirality of the crystal. Distinguishing between contributions originating from molecular chirality and from chiral crystal packing, we find that while IR absorption hardly depends on the symmetry of the space group, the situation is different for VCD, where completely new non-local patterns emerge. Understanding the two underlying mechanisms, namely gauge transport and direct coupling, will help to use VCD to distinguish polymorphic forms.

Time Resolved Raman Scattering of Molecules: A Quantum Mechanics Approach with Stochastic Schroedinger Equation

Raman scattering is a very powerful tool employed to characterize molecular systems. Here we propose a novel theoretical strategy to calculate the Raman cross-section in time domain, by computing the cumulative Raman signal emitted during the molecular evolution in time. Our model is based on a numerical propagation of the vibronic wave function under the effect of a light pulse of arbitrary shape. This approach can therefore tackle a variety of experimental setups. Both resonance and nonresonance Raman scattering can be retrieved, and also the time-dependent fluorescence emission is computed. The model has been applied to porphyrin considering both resonance and nonresonance conditions and varying the incident field duration. Moreover the effect of the vibrational relaxation, which should be taken into account when its time scale is similar to that of the Raman emission, has been included through the stochastic Schroedinger equation approach.

Anharmonic Vibrational Raman Optical Activity of Methyloxirane: Theory and Experiment Pushed to the Limits

Combining Raman scattering and Raman optical activity (ROA) with computer simulations reveals fine structural and physicochemical properties of chiral molecules. Traditionally, the region of interest comprised fundamental transitions within 200-1800 cm^-1. Only recently, nonfundamental bands could be observed as well. However, theoretical tools able to match the observed spectral features and thus assist their assignment are rather scarce. In this work, we present an accurate and simple protocol based on a three-quanta anharmonic perturbative approach that is fully fit to interpret the observed signals of methyloxirane within 150-4500 cm^-1. An unprecedented agreement even for the low-intensity combination and overtone transitions has been achieved, showing that anharmonic Raman and ROA spectroscopies can be valuable tools to understand vibrations of chiral molecules or to calibrate computational models.

Vibrational Circular Dichroism from DFT Molecular Dynamics: The AWV Method

The paper illustrates the Activity Weighted Velocities (AWV) methodology to compute Vibrational Circular Dichroism (VCD) anharmonic spectra from Density Functional Theory (DFT) molecular dynamics. AWV calculates the spectra by the Fourier Transform of the time correlation functions of velocities, weighted by specific observables: the Atomic Polar Tensors (APTs) and the Atomic Axial Tensors (AATs). Indeed, AWV shows to correctly reproduce the experimental spectra for systems in the gas and liquid phases, both in the case of weakly and strongly interacting systems. The comparison with the experimental spectra is striking especially in the fingerprint region, as demonstrated by the three benchmark systems discussed: (1S)-Fenchone in the gas phase, (S)-(−)-Propylene oxide in the liquid phase, and (R)-(−)-2-butanol in the liquid phase. The time evolution of APTs and AATs can be adequately described by a linear combination of the tensors of a small set of appropriate reference structures, strongly reducing the computational cost without compromising accuracy. Additionally, AWV allows the partition of the spectral signal in its molecular components without any expensive postprocessing and any localization of the charge density or the wave function.

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.

Calculation of absolute Raman scattering cross-sections using vibrational self-consistent field/vibrational configuration interaction wave functions

In the present study, the differential scattering cross-sections, depolarization ratios and Raman shifts of small molecular systems are obtained from configuration iteration wave functions of vibrational self-consistent field (VSCF) states. The transition polarizabilities were modeled using the Placzek approximation, neglecting those contributions not arising from the electric dipole mechanism. This theoretical approach is considered a good approximation for samples that absorb in the UV range if the excitation radiation falls in the visible region, as is the case of the molecules selected for the present study, namely: water, methane, and acetylene. Potential energy and electronic polarizability surfaces are calculated by the CCSD(T) and CC3 methods with aug-cc-p(C)V(T,Q,5)Z basis sets. The vibrational Hamiltonian includes the vibrational angular momentum contribution of the Watson kinetic energy operator. As expected, due to the variational nature of the VSCF and vibrational configuration interaction (VCI) methods, the Raman transition wavenumbers are substantially improved over the harmonic predictions. Surprisingly, the scattering cross-sections obtained using the harmonic approximation or the VSCF method better agrees with the experimental values than those cross-sections predicted using VCI wave functions. The more significant deviations of the VCI results from the experimental reference may be related to the significant uncertainties of the measured cross-sections. Still, it may also indicate that the VCI Raman transition moments may require a more accurate description of the electronic polarizability surface. Finally, the depolarization ratios calculated for H₂ O and C₂ D₂ using harmonic and VCI wave functions have similar accuracy, whereas, for C₂ H₂ and C₂ HD, the VCI results are more accurate.

Insights into the Thermally Activated Cyclization Mechanism in a Linear Phenylalanine-Alanine Dipeptide

Dipeptides, the prototype peptides, exist in both linear (l-) and cyclo (c-) structures. Since the first mass spectrometry experiments, it has been observed that some l-structures may turn into the cyclo ones, likely via a temperature-induced process. In this work, combining several different experimental techniques (mass spectrometry, infrared and Raman spectroscopy, and thermogravimetric analysis) with tight-binding and ab initio simulations, we provide evidence that, in the case of l-phenylalanyl-l-alanine, an irreversible cyclization mechanism, catalyzed by water and driven by temperature, occurs in the condensed phase. This process can be considered as a very efficient strategy to improve dipeptide stability by turning the comparatively fragile linear structure into the robust and more stable cyclic one. This mechanism may have played a role in prebiotic chemistry and can be further exploited in the preparation of nanomaterials and drugs.

Sum-frequency vibrational spectroscopy of methanol at interfaces due to Fermi resonance

We present a theoretical method of studying sum-frequency vibrational spectroscopy for the CH3 group of methanol at interfaces due to Fermi resonance, which provides a novel and untraditional point of view with respect to traditional approaches.

A Holistic Approach to Determining Stereochemistry of Potential Pharmaceuticals by Circular Dichroism with β-Lactams as Test Cases

This paper's main objective is to show that many different factors must be considered when solving stereochemical problems to avoid misleading conclusions and obtain conclusive results from the analysis of spectroscopic properties. Particularly in determining the absolute configuration, the use of chiroptical methods is crucial, especially when other techniques, including X-ray crystallography, fail, are not applicable, or give inconclusive results. Based on various β-lactam derivatives as models, we show how to reliably determine their absolute configuration (AC) and preferred conformation from circular dichroism (CD) spectra. Comprehensive CD analysis, employing both approaches, i.e., traditional with their sector and helicity rules, and state-of-the-art supported by quantum chemistry (QC) calculations along with solvation models for both electronic (ECD) and vibrational (VCD) circular dichroism ranges, allows confident defining stereochemistry of the β-lactams studied. Based on an in-depth analysis of the results, we have shown that choosing a proper chiroptical method/s strictly depends on the specific case and certain structural features.

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

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

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

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

Impact of fluorination and chlorination on the electronic structure, topology and in-plane ring normal modes of pyridines

Seven partially and fully fluorinated/chlorinated pyridines were investigated by means of FT-IR and Raman spectroscopy combined with quantum chemical calculations, mainly aiming to detect how the nature and position of F and Cl substituents affect the in-plane ring normal modes (RNMs) of pyridines in terms of vibrational wavenumbers, force constants, IR intensities and Raman activities. Taking pyridine as the reference, the RNMs and some derived RNMs through coupling with related C-X (X = F, Cl) stretching vibrations were identified on the basis of their composition in terms of internal coordinates. The impact of fluorination and chlorination on these RNMs was also discussed from the perspective of frontier molecular orbitals (MOs), maps of the molecular electrostatic potential (MEP) and the molecular topology. Natural bond orbital (NBO) analysis revealed the consequences of substitutions on the intramolecular charge delocalisation and consequently the ring bond strength. Moreover, the effects of anharmonicity of the potential on vibrational frequencies were presented and discussed.

Quantifying Conformational Isomerism in Chain Molecules by Linear Raman Spectroscopy: The Case of Methyl Esters

The conformational preferences of the ester group have the potential to facilitate the large amplitude folding of long alkyl chains in the gas phase. They are monitored by Raman spectroscopy in supersonic jet expansions for the model system methyl butanoate, after establishing a quantitative relationship with quantum-chemical predictions for methyl methanoate. This requires a careful analysis of experimental details, and a simulation of the rovibrational contours for near-symmetric top molecules. The technique is shown to be complementary to microwave spectroscopy in quantifying coexisting conformations. It confirms that a C-O-C(=O)-C-C chain segment can be collapsed into a single all-trans conformation by collisional cooling, whereas alkyl chain isomerism beyond this five-membered chain largely survives the jet expansion. This sets the stage for the investigation of linear alkyl alkanoates in terms of dispersion-induced stretched-chain to hairpin transitions by Raman spectroscopy.

AC/DC Analysis: Broad and Comprehensive Approach to Analyze Infrared Intensities at the Atomic Level

We present a complete theoretical protocol to partition infrared intensities into terms owing to individual atoms by two different but related approaches: the atomic contributions (ACs) show how the entire molecular vibrational motion affects the electronic structure of a single atom and the total infrared intensity. On the other hand, the dynamic contributions (DCs) show how the displacement of a single atom alters the electronic structure of the entire molecule and the total intensity. The two analyses are complementary ways of partitioning the same total intensity and conserve most of the features of the total intensity itself. Combined, they are called the AC/DC analysis. These can be further partitioned following the CCTDP (or CCT) models according to the population analysis chosen by the researcher. The main conceptual features of the equations are highlighted, and representative numerical results are shown to support the interpretation of the equations. The results are invariant to rotation and translation and can readily be extended to molecules of any size, shape, or symmetry. Although the AC/DC analysis requires the choice of a charge model, all charge models that correctly reproduce the total molecular dipole moment can be used. A fully automated protocol managed by the Placzek program is made available, free of charge and with input examples.

How much anharmonicity is in vibrational spectra of CH3I and CD3I?

This work presents new experimental and theoretical insights on vibrational spectra of CH₃I and CD₃I in the liquid phase. For the first time, we provided the contributions from different vibrational modes to mid-infrared (MIR) and near-infrared (NIR) spectra and estimated the extent of anharmonicity in the MIR region. Direct comparison of the intensities from ATR-IR and NIR transmission spectra was possible due to normalization of ATR-IR spectra. As a reference for normalization, we applied the area of the ν_s(CH₃)/ν_s(CD₃) band recorded in transmission mode. Our results show that the corresponding vibrational modes of CH₃I and CD₃I have similar contributions to the total intensity (MIR + NIR), however, these contributions are distributed in a different way between MIR and NIR regions. As expected, most of intensity in MIR spectra originates from the fundamental transitions (>90%). The fundamental bands together with the first overtones and the binary combinations contribute to more than 99% of MIR intensity for both compounds. Therefore, reliable reconstruction of MIR spectra can be achieved by considering only these vibrational modes. On the other hand, accurate simulation of NIR spectra requires including the higher-order transitions. In the case of CD₃I, the fourth-order transitions contribute to 12.7% of NIR intensity. The contributions from NIR region are significantly smaller than those from MIR range and were estimated to be 6.7% for CH₃I and 2.3% for CD₃I. The theoretical calculations provide a reasonable estimation of the total contribution from the fundamental bands. Yet, the calculated contributions from the anharmonic transitions are different from those obtained from the experimental data. MIR spectra of CH₃I and CD₃I reveal an unexpected increase in the intensity of some overtones and combination bands indicating the presence of Fermi resonances. These resonances are responsible for differences in contributions from the first overtones and binary combinations between CH₃I and CD₃I.

How to VPT2: Accurate and Intuitive Simulations of CH Stretching Infrared Spectra Using VPT2+K with Large Effective Hamiltonian Resonance Treatments

This article primarily discusses the utility of vibrational perturbation theory for the prediction of X-H stretching vibrations with particular focus on the specific variant, second-order vibrational perturbation theory with resonances (VPT2+K). It is written as a tutorial, reprinting most important formulas and providing numerous simple examples. It discusses the philosophy and practical considerations behind vibrational simulations with VPT2+K, including but not limited to computational method selection, cost-saving approximations, approaches to evaluating intensity, resonance identification, and effective Hamiltonian structure. Particular attention is given to resonance treatments, beginning with simple Fermi dyads and gradually progressing to arbitrarily large polyads that describe both Fermi and Darling-Dennison resonances. VPT2+K combined with large effective Hamiltonians is shown to be a reliable framework for modeling the complicated CH stretching spectra of alkenes. An error is also corrected in the published analytic formula for the VPT2 transition moment between the vibrational ground state and triply excited states.

Lactic Acid Spectroscopy: Intra- and Intermolecular Interactions

Lactic acid, a relevant molecule in biology and the environment, is an α-hydroxy acid with a high propensity to form hydrogen bonds, both internally and to other hydrogen-bond-accepting molecules. This work includes the novel recording of infrared spectra of gas-phase lactic acid using Fourier transform infrared spectroscopy, and the vibrational absorption features of lactic acid are assigned with the aid of computationally simulated vibrational spectra with anharmonic corrections. Theoretical chemistry methods are used to relate intramolecular hydrogen-bond strengths to the relative stability of lactic acid conformers. The formation of hydrogen-bonded lactic acid dimers and 1:1 water complexes is investigated by simulated vibrational spectra and calculated thermodynamic parameters for the lactic acid monomer and dimer and its water complex in the gas phase. The results of this study are discussed in the context of environmental chemistry with an emphasis on indoor environments.

Basis set dependence of S[double bond, length as m-dash]O stretching frequencies and its consequences for IR and VCD spectra predictions

Benchmarking functionals and basis sets for the computational prediction of molecular properties is usually done on very small model systems. Larger organic molecules containing heavier second row atoms are not the typical model structures. We herein present the first survey of basis sets and functionals for the prediction of the IR and VCD spectra of chiral tosylates and sulfinates as we noted drastic deviations between computed harmonic frequencies obtained at B3LYP/6-311++G(2d,p) level of theory and those observed in experimental solution phase IR and VCD spectra. We show that the harmonic frequencies of the asymmetric and symmetric S[double bond, length as m-dash]O stretching modes of tosylates are predicted at significantly too low vibrational frequencies if the employed basis set does not provide higher order polarization functions. The results of our benchmarks show that at least the 6-311G(3df,2dp) basis (or equivalent Dunning and Ahlrichs variants) should be used.

Looking for the bricks of the life in the interstellar medium: The fascinating world of astrochemistry

The discovery in the interstellar medium of molecules showing a certain degree of complexity, and in particular those with a prebiotic character, has attracted great interest. A complex chemistry takes place in space, but the processes that lead to the production of molecular species are a matter of intense discussion, the knowledge still being at a rather primitive stage. Debate on the origins of interstellar molecules has been further stimulated by the identification of biomolecular building blocks, such as nucleobases and amino acids, in meteorites and comets. Since many of the molecules found in space play a role in the chemistry of life, the issue of their molecular genesis and evolution might be related to the profound question of the origin of life itself. Understanding the underlying chemical processes, including the production, reactions and destruction of compounds, requires the concomitant study of spectroscopy, gas-phase reactivity, and heterogeneous processes on dust-grains. The aim of this contribution is to provide a general view of a complex and multifaceted challenge, while focusing on the role played by molecular spectroscopy and quantum-chemical computations. In particular, the derivation of the molecular spectroscopic features and the investigation of gas-phase formation routes of prebiotic species in the interstellar medium are addressed from a computational point of view.

Accuracy Meets Interpretability for Computational Spectroscopy by Means of Hybrid and Double-Hybrid Functionals

Accuracy and interpretability are often seen as the devil and holy grail in computational spectroscopy and their reconciliation remains a primary research goal. In the last few decades, density functional theory has revolutionized the situation, paving the way to reliable yet effective models for medium size molecules, which could also be profitably used by non-specialists. In this contribution we will compare the results of some widely used hybrid and double hybrid functionals with the aim of defining the most suitable recipe for all the spectroscopic parameters of interest in rotational and vibrational spectroscopy, going beyond the rigid rotor/harmonic oscillator model. We will show that last-generation hybrid and double hybrid functionals in conjunction with partially augmented double- and triple-zeta basis sets can offer, in the framework of second order vibrational perturbation theory, a general, robust, and user-friendly tool with unprecedented accuracy for medium-size semi-rigid molecules.

Graphene Domain Signature of Raman Spectra of sp2 Amorphous Carbons

The standard D-G-2D pattern of Raman spectra of sp² amorphous carbons is considered from the viewpoint of graphene domains presenting their basic structure units (BSUs) in terms of molecular spectroscopy. The molecular approximation allows connecting the characteristic D-G doublet spectra image of one-phonon spectra with a considerable dispersion of the C=C bond lengths within graphene domains, governed by size, heteroatom necklace of BSUs as well as BSUs packing. The interpretation of 2D two-phonon spectra reveals a particular role of electrical anharmonicity in the spectra formation and attributes this effect to a high degree of the electron density delocalization in graphene domains. A size-stimulated transition from molecular to quasi-particle phonon consideration of Raman spectra was experimentally traced, which allowed evaluation of a free path of optical phonons in graphene crystal.

Absolute Configuration and Conformation of (-)-R-t-Butylphenylphosphinoamidate: Chiroptical Spectroscopy and X-ray Analysis

The absolute configuration and conformations of (−)-tert-butylphenylphosphinoamidate were determined using three different chiroptical spectroscopic methods, namely vibrational circular dichroism (VCD), electronic circular dichroism (ECD), and optical rotatory dispersion (ORD). In each of the spectroscopic methods used, experimental data for the (−)-enantiomer of tert-butylphenylphosphinoamidate were measured in the solution phase. Using the concentration-dependent experimental infrared spectra, the existence of dimers in the solution was investigated, and the monomer–dimer equilibrium constant was determined. Concomitant quantum mechanical predictions of the VCD, ECD, and ORD for monomeric tert-butylphenylphosphinoamidate were carried out using density functional theory (DFT) calculations using the B3LYP functional and the 6-31G(d), 6-311G(2d,2p) and aug-cc-pVDZ basis sets. Similar predictions for dimeric tert-butylphenylphosphinoamidate were also obtained using the B3LYP/6-31G(d) method. A comparison of theoretically predicted data with the corresponding experimental data led to the elucidation of the absolute configuration as (−)-(R)-tert-butylphenylphosphinoamidate with one predominant conformation in the solution. This conclusion was independently supported by X-ray analysis of the complex with (+)-R-2,2′-dihydroxy-1,1′-binaphthol ((+)-R- BINOL).

Vibrations of the guanine-cytosine pair in chloroform: an anharmonic computational study

We compute at the anharmonic level the vibrational spectra of the Watson-Crick dimer formed by guanosine (G) and cytidine (C) in chloroform, together with those of G, C and the most populated GG dimer. The spectra for deuterated and partially deuterated GC are also computed. We use DFT calculations, with B3LYP and CAM-B3LYP as reference functionals. Solvent effects from chloroform are included via the Polarizable Continuum Model (PCM), and by performing tests on models including up two chloroform molecules. Both B3LYP and CAM-B3LYP calculations reproduce the shape of the experimental spectra well in the fingerprint region (1500-1700 cm-1) and in the N-H stretching region (2800-3600 cm-1), with B3LYP providing better quantitative agreement with experiments. According to our calculations, the N-H amido streching mode of G falls at ∼2900 cm-1, while the N-H amino of G and C falls at ∼3100 cm-1 when hydrogen-bonded, or ∼3500 cm-1 when free. Overtone and combination bands strongly contribute to the absorption band at ∼3300 cm-1. Inclusion of bulk solvent effects significantly increases the accuracy of the computed spectra, while solute-solvent interactions have a smaller, though still noticeable, effect. Some key aspects of the anharmonic treatment of strongly vibrationally coupled supermolecular systems and the related methodological issues are also discussed.

Time-Resolved Femtosecond Stimulated Raman Spectra and DFT Anharmonic Vibrational Analysis of an Electronically Excited Rhenium Photosensitizer

Time-resolved femtosecond stimulated Raman spectra (FSRS) of a prototypical organometallic photosensitizer/photocatalyst ReCl(CO)₃(2,2'-bipyridine) were measured in a broad spectral range ∼40-2000 (4000) cm^-1 at time delays from 40 fs to 4 ns after 400 nm excitation of the lowest allowed electronic transition. Theoretical ground- and excited-state Raman spectra were obtained by anharmonic vibrational analysis using second-order vibrational perturbation theory on vibrations calculated by harmonic approximation at density functional theory-optimized structures. A good match with anharmonically calculated vibrational frequencies allowed for assigning experimental Raman features to particular vibrations. Observed frequency shifts upon excitation (ν(ReCl) and ν(CC inter-ring) vibrations upward; ν(CC, CN) and ν(Re-C) downward) are consistent with the bonding/antibonding characters of the highest occupied molecular orbital and the lowest unoccupied molecular orbital involved in excitation and support the delocalized formulation of the lowest triplet state as ReCl(CO)₃ → bpy charge transfer. FSRS spectra show a mode-specific temporal evolution, providing insights into the intersystem crossing (ISC) mechanism and subsequent relaxation. Most of the Raman features are present at ∼40 fs and exhibit small shifts and intensity changes with time. The 1450-1600 cm^-1 group of bands due to CC, CN, and CC(inter-ring) stretching vibrations undergoes extensive restructuring between 40 and ∼150 fs, followed by frequency upshifts and a biexponential (0.38, 21 ps) area growth, indicating progressing charge separation in the course of the formation and relaxation of the lowest triplet state. Early (40-150 fs) restructuring was also observed in the low-frequency range for ν(Re-Cl) and δ(Re-C-O) vibrations that are presumably activated by ISC. FSRS experimental innovations employed to measure low- and high-energy Raman features simultaneously are described and discussed in detail.

Toward Fully Unsupervised Anharmonic Computations Complementing Experiment for Robust and Reliable Assignment and Interpretation of IR and VCD Spectra from Mid-IR to NIR: The Case of 2,3-Butanediol and trans-1,2-Cyclohexanediol

The infrared (IR) and vibrational circular dichroism (VCD) spectra of 2,3-butanediol and trans-1,2-cyclohexanediol from 900 to 7500 cm^–1 (including mid-IR, fundamental CH and OH stretchings, and near-infrared regions) have been investigated by a combined experimental and computational strategy. The computational approach is rooted in density functional theory (DFT) computations of harmonic and leading anharmonic mechanical, electrical, and magnetic contributions, followed by a generalized second-order perturbative (GVPT2) evaluation of frequencies and intensities for all the above regions without introducing any ad hoc scaling factor. After proper characterization of large-amplitude motions, all resonances plaguing frequencies and intensities are taken into proper account. Comparison of experimental and simulated spectra allows unbiased assignment and interpretation of the most interesting features. The reliability of the GVPT2 approach for OH stretching fundamentals and overtones is confirmed by the remarkable agreement with a local mode model purposely tailored for the latter two regions. Together with the specific interest of the studied molecules, our results confirm that an unbiased assignment and interpretation of vibrational spectra for flexible medium-size molecules can be achieved by means of a nearly unsupervised reliable, robust, and user-friendly DFT/GVPT2 model.

Vibrational optical activity for structural characterization of natural products

This review presents the recent progress towards elucidating the structures of chiral natural products and applications using vibrational optical activity (VOA) spectroscopy.

The challenging playground of astrochemistry: an integrated rotational spectroscopy - quantum chemistry strategy

While it is now well demonstrated that the interstellar medium (ISM) is characterized by a diverse and complex chemistry, a significant number of features in radioastronomical spectra are still unassigned and call for new laboratory efforts, which are increasingly based on integrated experimental and computational strategies. In parallel, the identification of an increasing number of molecules containing more than five atoms and at least one carbon atom (the so-called "interstellar" complex organic molecules), which can play a relevant role in the chemistry of life, raises the additional issue of how these species can be produced in the typical harsh conditions of the ISM. On these grounds, this perspective aims to present an integrated rotational spectroscopy - quantum chemistry approach for supporting radioastronomical observations and a computational strategy for contributing to the elucidation of chemical reactivity in the interstellar space.

Spectroscopic evidence of a particular intermolecular interaction in iodomethane–ethanol mixtures: the cooperative effect of halogen bonding, hydrogen bonding, and the solvent effect

The particular intermolecular interaction of an iodomethane–ethanol mixture is revealed by NIR, Raman, DFT calculation, and 2D correlation analysis.

Evaluation of Molecular Polarizability and of Intensity Carrying Modes Contributions in Circular Dichroism Spectroscopies

We re-examine the theory of electronic and vibrational circular dichroism spectroscopy in terms of the formalism of frequency-dependent molecular polarizabilities. We show the link between Fermi’s gold rule in circular dichroism and the trace of the complex electric dipole–magnetic dipole polarizability. We introduce the C++ code polar to compute the molecular polarizability complex tensors from quantum chemistry outputs, thus simulating straightforwardly UV-visible absorption (UV-Vis)/electronic circular dichroism (ECD) spectra, and infrared (IR)/vibrational circular dichroism (VCD) spectra. We validate the theory and the code by referring to literature data of a large group of chiral molecules, showing the remarkable accuracy of density functional theory (DFT) methods. We anticipate the application of this methodology to the interpretation of vibrational spectra in various measurement conditions, even in presence of metal surfaces with plasmonic properties. Our theoretical developments aim, in the long run, at embedding the quantum-mechanical details of the chiroptical spectroscopic response of a molecule into the simulation of the electromagnetic field distribution at the surface of plasmonic devices. Such simulations are also instrumental to the interpretation of the experimental spectra measured from devices designed to enhance chiroptical interactions by the surface plasmon resonance of metal nanostructures.