Discovering the Elusive Global Minimum in a Ternary Chiral Cluster: Rotational Spectra of Propylene Oxide Trimer

Prebiotic Astrochemistry from Astronomical Observations and Laboratory Spectroscopy

The discovery of more than 200 gas-phase chemical compounds in interstellar space has led to the speculation that this nonterrestrial synthesis may play a role in the origin of life. These identifications were possible because of laboratory spectroscopy, which provides the molecular fingerprints for astronomical observations. Interstellar chemistry produces a wide range of small, organic molecules in dense clouds, such as NH2COCH3, CH3OCH3, CH3COOCH3, and CH2(OH)CHO. Carbon (C) is also carried in the fullerenes C60 and C70, which can preserve C-C bonds from circumstellar environments for future synthesis. Elusive phosphorus has now been found in molecular clouds, the sites of star formation, in the molecules PO and PN. Such clouds can collapse into solar systems, although the chemical/physical processing of the emerging planetary disk is uncertain. The presence of molecule-rich interstellar starting material, as well as the link to planetary bodies such as meteorites and comets, suggests that astrochemical processes set a prebiotic foundation.

The Heavy Atom Structure, "cis effect" on Methyl Internal Rotation, and 14N Nuclear Quadrupole Coupling of 1-Cyanopropene from Quantum Chemical and Microwave Spectroscopic Analysis

The microwave spectrum of 1-cyanopropene (crotonitrile) was remeasured using two pulsed molecular jet Fourier transform microwave spectrometers operating from 2.0 to 40.0 GHz. The molecule exists in two isomer forms, E and Z, with respect to the orientation between the methyl and the cyano groups. The spectrum of the Z isomer is more intense. Due to internal rotation of the methyl group, doublets containing A and E torsional species were found for all rotational transitions. Hyperfine splittings arising from the 14N nuclear quadrupole coupling were resolved. The heavy atom structure of the Z isomer was determined by observation of 13C and 15N isotopologue spectra in natural abundances. The experimental results were supported by quantum chemistry. The complex spectral patterns were analyzed and fitted globally, and the barriers to methyl internal rotation are determined to be 478.325(28) cm-1 and 674.632(76) cm-1 for the Z and E isomers, respectively. The non-bonded intramolecular electrostatic attraction between the methyl group and the 1-cyano substituent overcomes steric hindrance, leading to higher stability of the Z isomer. The consequence is a slight opening of 3.2° of the C(1)-C(2)-C(3) angle and a radical decrease of the methyl torsional barrier in the Z isomer due to steric repulsion.

Binary conformers of a flexible, long-chain fluoroalcohol: dispersion controlled selectivity and relative abundances in a jet

The complex conformational panorama of binary 4,4,4-trifluoro-1-butanol (TFB) aggregates was investigated using chirped-pulse Fourier transform microwave spectroscopy, aided by conformational searches using CREST (Conformer-Rotamer Ensemble Sampling Tool) and quantum chemistry calculations. From nearly 1500 initial dimer geometries, 16 most stable binary candidates were obtained within a relative energy window of ∼4 kJ mol-1. Rotational spectra of five binary conformers were experimentally observed in supersonic expansion and assigned. Interestingly, three out of the five observed binary conformers are composed solely of monomer conformers, which were not observed in their isolated gas phase forms in jet expansion. In addition, an observed dimer that is made exclusively of the most stable TFB monomer subunits does not correspond to the global minimum. The intricate kinetically and thermodynamically controlled dimer formation mechanisms are discussed, and a modified kinetic-thermodynamic model was developed, providing conformational abundances that are in good agreement with the experiment. Subsequent non-covalent interaction analyses reveal that the observed conformers are held together by one primary O-H⋯O hydrogen bond and secondary intermolecular C-H⋯O, C-H⋯F, and/or O-H⋯F interactions, as well as C-H⋯H-C London dispersion interactions between the methylene groups. Further symmetry-adapted perturbation theory analyses of the TFB dimer conformers and related alcohol dimers reveal a considerable rise in dispersion contributions with increasing n-alkyl carbon chain length and highlight the role of dispersion interactions in preferentially stabilizing the global minimum of the TFB dimer.

Conformational distributions of tetrahydro-2-turoic acid in water at different pH values by their IR and vibrational circular dichroism spectra

Infrared (IR) and vibrational circular dichroism (VCD) spectra of tetrahydro-2-furoic acid (THFA) in aqueous solutions under several different pH conditions were recorded. To interpret the IR and VCD spectra of THFA obtained in highly acidic and basic aqueous solutions, extensive and systematic conformational searches were conducted to acquire the low-energy minima for both the neutral and deprotonated forms of THFA species, as well as their hydrated clusters. This was accomplished by using the conformer-rotamer ensemble sampling tool (CREST) with an implicit solvation model for water. The CREST candidates were further optimized at the B3LYP-D3BJ/def2-TZVP level of theory. The simulated VCD spectra of the neutral THFA conformers in the polarizable continuum model (PCM) of water alone exhibit little agreement with the experimental data under highly acidic conditions. Applying the clusters-in-a-liquid solvation model by considering the monohydrate THFA conformers in the PCM of water, significantly improved agreement with the experimental data. Similarly, the deprotonated THFA species solvated with one to four explicit water molecules in the PCM of water were considered. While the IR and VCD spectra of the deprotonated THFA monohydrate conformers offer the best agreement with the experimental data, other larger hydrated clusters, particularly the dihydrates, also contribute to the spectra. Through the synergistic combined experimental and theoretical approach used in the study, comprehensive conformational distributions of the predominant THFA species across various pH conditions were extracted.

Quantum chemistry meets high-resolution spectroscopy for characterizing the molecular bricks of life in the gas-phase

Computation of accurate geometrical structures and spectroscopic properties of large flexible molecules in the gas-phase is tackled at an affordable cost using a general exploration/exploitation strategy. The most distinctive feature of the approach is the careful selection of different quantum chemical models for energies, geometries and vibrational frequencies with the aim of maximizing the accuracy of the overall description while retaining a reasonable cost for all the steps. In particular, a composite wave-function method is used for energies, whereas a double-hybrid functional (with the addition of core-valence correlation) is employed for geometries and harmonic frequencies and a cheaper hybrid functional for anharmonic contributions. A thorough benchmark based on a wide range of prototypical molecular bricks of life shows that the proposed strategy is close to the accuracy of state-of-the-art composite wave-function methods, and is applicable to much larger systems. A freely available web-utility post-processes the geometries optimized by standard electronic structure codes paving the way toward the accurate yet not prohibitively expensive study of medium- to large-sized molecules by experimentally-oriented researchers.

Synthesis and absolute configuration of cyclic synthetic cathinones derived from α-tetralone

The emergence of new synthetic cathinones continues to be a matter of public health concern. In fact, already known products (drugs) are being rapidly replaced by new structurally related alternatives, whereby modifications in the basic cathinone structure are used by manufacturers to circumvent the legislation. On the other hand, some derivatives of synthetic cathinones represent important pharmaceuticals with antidepressant properties. In the search for pharmaceutically relevant analogs, the main goal of the present study was to design and characterize novel cyclic α-tetralone-based derivatives of synthetic cathinones. We synthesized a series of derivatives and verified their chemical structure. Subsequently, chiral separation has been accomplished by high-performance liquid chromatography (HPLC) equipped with a circular dichroism (CD) detector, which directly provided CD spectra of the enantiomers of the analyzed substances at 252 nm. Using density functional theory calculations, we have obtained stable conformers of selected enantiomers in solution and their relative abundances, which we used to simulate their spectra. The experimental and calculated data have been used to assign the absolute configuration of six as-yet unknown synthetic cathinones.

Wetting vs Droplet Aggregation: A Broadband Rotational Spectroscopic Study of 3-Methylcatechol⋅⋅⋅Water Clusters

Two competing solvation pathways of 3-methylcatechol (MC), an atmospherically relevant aromatic molecule, with up to five water molecules were explored in detail by using a combination of broadband rotational spectroscopy and computational chemistry. Theoretically, two different pathways of solvation emerge: the commonly observed droplet pathway which involves preferential binding among the water molecules while the solute serves as an anchor point for the formation of a water cluster, and an unexpected wetting pathway which involves interactions between the water molecules and the aromatic face of MC, i.e., a wetting of the π-surface. Conclusive identification of the MC hydrate structures, and therefore the wetting pathway, was facilitated by rotational spectra of the parent MC hydrates and several H2 18 O and 13 C isotopologues which exhibit splittings associated with methyl internal rotation and/or water tunneling motions. Theoretical modelling and analyses offer insights into the tunneling and conversion barriers associated with the observed hydrate conformers and the nature of the non-covalent interactions involved in choosing the unusual wetting pathway.

Deciphering the non-covalent interactions in the furan⋯hexane complex using rotational spectroscopy and theoretical analyses

The rotational spectrum of a binary complex formed between furan and n-hexane was investigated using a chirped pulse Fourier transform microwave spectrometer in the range of 2-6 GHz. While furan has only one conformer, n-hexane exists in multiple conformations. The conformational landscape of the binary complex was systematically explored by using a semiempirical conformational search tool, namely CREST. The CREST conformational candidates were subjected to further geometry optimization and harmonic frequency calculations at the B3LYP-D3BJ/def2-TZVP level of theory, resulting in 34 minima within an energy window of 5 kJ mol-1. The three most stable furan⋯hexane minima all contain the most stable n-hexane conformer subunit and are separated by relatively low conformational conversion barriers. Additional calculations were carried out to support the conclusive identification of the global minimum structure responsible for the set of assigned rotational transitions. These include calculations at the B3LYP-D3BJ level with the aug-cc-pVTZ and 6-311++G(d,p) basis sets and the MP2/def2-TZVP level, as well as the single point energy calculations at the CCSD(T)-F12/cc-pVDZ level. Further non-covalent interaction and principal interacting orbital analyses show that the synergy of the πfuran → σ*hexane and σhexane → π*furan interactions plays an important role in stabilizing the observed furan-hexane conformer.

Infrared and vibrational circular dichroism spectra of methyl β-D-glucopyranose in water: The application of the quantum cluster growth and clusters-in-a-liquid solvation models

The infrared (IR) and vibrational circular dichroism (VCD) spectra of methyl β-D-glucopyranose in water were measured. Both implicit and explicit solvation models were utilized to explain the observed spectra. The vast body of existing experimental and theoretical data suggested that about eight explicit water molecules are needed to account for the solvent effects, supported by the current Quantum Cluster Growth (QCG) analysis. Extensive manual and systematic conformational searches of the molecular target and its water clusters were carried out by using a recently developed conformational searching tool, conformer-rotamer ensemble sampling tool (CREST), and the microsolvation model in the associated QCG code. The Boltzmann averaged IR and VCD spectra of the methyl β-D-glucopyranose-(water)n (n = 8) conformers in the PCM of water provide better agreement with the experimental ones than those with n = 0, 1, and 2. The explicit solvation with eight water molecules was shown to greatly modify the conformational preference of methyl β-D-glucopyranose from its monomeric form. Further analyses show that the result is consistent with the existence of long-lived methyl β-D-glucopyranose monohydrates with the additional explicit water effects being accounted for with the quantum mechanical treatment of the other seven close-by water molecules in the PCM of water.

Conformational Distributions of Phenyl β-D-Glucopyranoside and Gastrodin in Solution by Vibrational Optical Activity and Theoretical Calculations

The conformational landscapes of two highly flexible monosaccharide derivatives, namely phenyl β-D-glucopyranoside (ph-β-glu) and 4-(hydroxymethyl)phenyl β-D-glucopyranoside, also commonly known as gastrodin, were explored using a combined experimental and theoretical approach. For the infrared, Raman, and the associated vibrational optical activity (VOA), i.e., vibrational circular dichroism and Raman optical activity, experiments of these two compounds in DMSO and in water were carried out. Extensive and systematic conformational searches were performed using a recently developed conformational searching tool called CREST (conformer-rotamer ensemble sampling tool) in the two solvents. Fourteen and twenty-four low-energy conformers were identified at the DFT level for ph-β-glu and gastrodin, respectively. The spectral simulations of individual conformers were done at the B3LYP-D3BJ/def2-TZVPD level with the polarizable continuum model of the solvents. The VOA spectral features exhibit much higher specificity to conformational differences than their parent infrared and Raman. The excellent agreements achieved between the experimental and simulated VOA spectra allow for the extraction of experimental conformational distributions of these two carbohydrates in solution directly. The experimental percentage abundances based on the hydroxymethyl (at the pyranose ring) conformations G+, G-, and T for ph-β-glu were obtained to be 15%, 75%, and 10% in DMSO and 53%, 40%, and 7% in water, respectively, in comparison to the previously reported gas phase values of 68%, 25%, and 7%, highlighting the important role of solvents in conformational preferences. The corresponding experimental distributions for gastrodin are 56%, 22%, and 22% in DMSO and 70%, 21%, and 9% in water.

Indazole-derived synthetic cannabinoids: Absolute configuration determination and structure characterization by circular dichroism and DFT calculations

An increasing number of products containing synthetic cannabinoids pose a growing crisis to public health worldwide. Recently, a rising number of cases of serious adverse health effects, intoxications, and death cases associated with synthetic cannabinoids were reported. The current study represents the comprehensive structural analysis of three new synthetic cannabinoids (AB-, ADB- and AMB-FUBINACA) in solution investigated by electronic and vibrational circular dichroism together with the conventional methods of infrared and ultraviolet absorption spectroscopy, all supported by the density functional theory (DFT) calculations. The best level of theory to reproduce the experimental wavenumbers and wavelengths was found to be the B3PW91 method with a 6-311++G(d,p) basis set including the implicit solvent effect simulation. Very good agreement between the experimental and simulated spectra allowed us to determine the absolute configuration and a detailed interpretation of the IR absorption, VCD, ECD and UV spectra of AB-, ADB- and AMB-FUBINACA. In addition, the HOMO and LUMO electronic transitions were calculated.

CO2 Aggregation on Monoethanolamine: Observations from Rotational Spectroscopy

The initial stages of the gas-phase nucleation between CO₂ and monoethanolamine were investigated via broadband rotational spectroscopy with the aid of extensive theoretical structure sampling. Sub-nanometer-scale aggregation patterns of monoethanolamine-(CO₂ )_n , n=1-4, were identified. An interesting competition between the monoethanolamine intramolecular hydrogen bond and the intermolecular interactions between monoethanolamine and CO₂ upon cluster growth was discovered, revealing an intriguing CO₂ binding priority to the hydroxyl group over the amine group. These findings are in sharp contrast to the general results for aqueous solutions. In the quinary complex, a cap-like CO₂ tetramer was observed cooperatively surrounding the monoethanolamine. As the cluster approaches the critical size of new particle formation, the contribution of CO₂ self-assembly to the overall stability increases.

Chirality Discrimination at Binary Organic|Water Interfaces Monitored by Interfacial Tension Measurements with Preliminary Comparison with Molecular Dynamics Simulations

Chirality discrimination at a binary toluene (organic)/water(aqueous) interface between R- or S-Tol-BINAP (2,2'-Bis(di-p-tolylphosphino)-1,1'-binaphthyl) molecules and the water-soluble serine chiral specie is examined for the first time, using a combination of interfacial tension measurements and molecular dynamic simulations. Experimental interfacial measurements exhibit a clear chirality-controlled difference when a homochiral versus a heterochiral enantiomeric pairs are introduced at the interfaces. The related molecular dynamics simulations support the experimental results and provide further molecular insight of intermolecular interactions at the interfaces. The results indicate that interfacial tension measurements can capture the preferential interactions which exist between different pairs of enantiomers at the binary interfaces, opening up a new way for probing chirality discrimination at liquid-liquid interfaces.

Connections between the accuracy of rotational constants and equilibrium molecular structures

Rotational spectroscopy is the technique of choice for investigating molecular structures in the gas phase. Indeed, rotational constants are strongly connected to the geometry of the molecular system under consideration. Therefore, they are powerful tools for assessing the accuracy that quantum chemical approaches can reach in structural determinations. In this review article, it is shown how it is possible to measure the accuracy of a computed equilibrium geometry based on the comparison of rotational constants. But, it is also addressed what accuracy is required by computations for providing molecular structures and thus rotational constants that are useful to experiment. Quantum chemical methodologies for obtaining the "0.1% accuracy" for rotational constants are reviewed for systems ranging in size from small molecules to small polycyclic aromatic hydrocarbons. This accuracy for systems containing two dozen or so atoms opens the way towards future applications such as the accurate characterization of non-covalent interactions, which play a key role in several biological and technological processes.

Conformations of Steroid Hormones: Infrared and Vibrational Circular Dichroism Spectroscopy

Steroid hormone molecules may exhibit very different functionalities based on the associated functional groups and their 3D arrangements in space, i.e., absolute configurations and conformations. Infrared (IR) and vibrational circular dichroism (VCD) spectra of four different steroid hormones, namely dehydroepiandrosterone (DHEA), 17α-methyltestosterone (MTTT), (16α,17)-epoxyprogesterone (Epoxy-P4), and dehydroepiandrosterone acetate (AcO-DHEA), were measured in deuterated dimethyl sulfoxide and some also in carbon tetrachloride. Extensive conformational searches were carried out using the recent developed conformer-rotamer ensemble sampling tool (CREST) which also accounts for solvent effects using an implicit solvation model. All the CREST conformational candidates were then reoptimized at the B3LYP-D3BJ/def2-TZVPD with the PCM of solvent. The good agreements between the experimental IR and VCD spectra and the theoretical simulations provide a conclusive information about their conformational distribution and absolute configurations. The experimental and theoretical IR and VCD spectra of AcO-DHEA in the carbonyl and alkene stretching region showed some discrepancies, and the possible causes related to solvent effects, large amplitude motions and levels of theory used in the modelling were explored in detail. As part of the investigation, additional calculations at the B3LYP-D3BJ/6-31++G (2d,p) and B3LYP-D3BJ/cc-pVTZ levels, as well as some 'mixed' calculations with the double-hybrid functional B2PLYP-D3 were also carried out. The results indicate that the double-hybrid functional is important for predicting the correct IR band pattern in the carbonyl and alkene stretching region.

DFT exchange: sharing perspectives on the workhorse of quantum chemistry and materials science

In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners. The format of the paper is that of a roundtable discussion, in which the participants express and exchange views on DFT in the form of 302 individual contributions, formulated as responses to a preset list of 26 questions. Supported by a bibliography of 777 entries, the paper represents a broad snapshot of DFT, anno 2022.

Rotational Spectroscopy of the 2,2,3,3,3-Pentafluoropropanol⋅⋅⋅Water Complex: Conformations and Large Amplitude Motions

The 2,2,3,3,3-pentafluoropropanol (PFP) monomer can exist in five conformations defined by the CCCO and CCOH dihedral angles: four mirror-imaged pairs (G+g+/G-g-, G+g-/G-g+, G+t/G-t, Tg+/Tg-) and an achiral Tt form. We examined the conformational landscape of the PFP⋅⋅⋅water complex using chirped pulsed Fourier transform microwave spectroscopy and theoretical calculations. Rotational spectra of two PFP⋅⋅⋅water conformers, PFPG+g+⋅⋅⋅W_H and PFPTg+⋅⋅⋅W_H , and seven deuterated isotopologues of each, were assigned. Tunneling splittings were observed for both conformers and are attributed to the exchange of the bonded and non-bonded hydrogen atoms of water. On the other hand, the tunneling splitting associated with the OH flipping motion in PFPTg+/Tg- appears to be quenched upon hydrogen bonding with water. The large amplitude motions associated with the water subunits were examined in detail to explain the very different magnitudes of the experimental and theoretical permanent electric dipole moment components. The study highlights the challenge in correctly identifying the conformers observed when large amplitude motions are involved. Quantum theory of atoms in molecules (QTAIM) and non-covalent interaction (NCI) analyses, as well as electrostatic potential (ESP) calculations were carried out to explore the nature of the non-covalent interactions and to appreciate the effects of fluorination.

Sniffing out camphor: the fine balance between hydrogen bonding and London dispersion in the chirality recognition with α-fenchol

Binary complexes between the chiral monoterpenoids camphor and α-fenchol were explored with vibrational and rotational jet spectroscopy as well as density functional theory in order to explore how chirality can influence the binding preferences in gas-phase complexes. The global minimum structures of the two diastereomers were assigned. It is found that chirality recognition leads to different compromises in the fine balance between intermolecular interactions. While one isomer features a stronger hydrogen bond, the other one is more tightly arranged and stabilized by larger London dispersion interactions. These new spectroscopic results help understand the influence of chirality in molecular aggregation and unveil the kind of interactions involved between a chiral alcohol and a chiral ketone with large dispersion contributions.

The First Stages of Nanomicelle Formation Captured in the Sevoflurane Trimer

Self-aggregation of sevoflurane, an inhalable, fluorinated anesthetic, provides a challenge for current state-of-the-art high-resolution techniques due to its large mass and the variety of possible hydrogen bonds between monomers. Here we present the observation of sevoflurane trimer by chirped-pulse Fourier transform microwave spectroscopy, identified through the interplay of experimental and computational methods. The trimer (>600 Da), one of the largest molecular aggregates observed through rotational spectroscopy, does not resemble the binding (C-H···O) motif of the already characterized sevoflurane dimer, instead adapting a new binding configuration created predominantly from 17 CH···F hydrogen bonds that resembles a nanomicellar arrangement. The observation of such a heavy aggregate highlights the potential of rotational spectroscopy to study larger biochemical systems in the limit of spectroscopic congestion but also showcases the challenges ahead as the mass of the system increases.

Raman Optical Activity of N-Acetyl-L-Cysteine in water and in methanol: the "clusters-in-a-liquid" model and ab initio molecular dynamics simulations

Raman and Raman Optical Activity (ROA) spectra of N-acetyl-L-cysteine (NALC), a flexible chiral molecule, were measured in water and in methanol to evaluate the solvent effects. Two different solvation approaches, i.e. the DFT based clusters-in-a-liquid solvent model and the ab initio molecular dynamics (AIMD) simulations, were applied to simulate the Raman and ROA spectra. Systematic conformational searches were carried out using a recently developed conformational searching tool, CREST, with the inclusion of polarizable continuum model of water and of methanol. The CREST candidates of NALC and the NALC-solvent complexes were re-optimized and their Raman and ROA simulations were done at the B3LYP-D3BJ/def2-TZVP and the B3LYP-aug-cc-pVDZ//cc-pVTZ levels. Also, AIMD simulations , which includes some anharmonic effects and all intermolecular interactions in solution, were performed. By empirically weighting the computed Raman and ROA spectra of each conformer, good agreements with the experimental data were achieved with both approaches, while AIMD offered some improvements in the carbonyl and in the low wavenumber regions over the static DFT approach. The pros and cons of these two different approaches for accounting the solvent effects on Raman and ROA of this flexible chiral system will also be discussed.

Modulation of π character upon complexation captured by molecular rotation spectra

Two configurations of the furan–CF3Cl complex have been observed by high-resolution rotational spectroscopy. One is characterized by a dominant Cl lone pairs∙∙∙π*aromatic interaction and the other is stabilized by a...

2,2,3,3,3-pentafluoro-1-propanol and its dimer: structural diversity, conformational conversion, and tunnelling motion

Rotational spectra of the 2,2,3,3,3-pentafluoro-1-propanol (PFP) were measured using cavity and chirped pulse Fourier transform microwave spectrometers. Of the nine possible PFP configurations which include four mirror-imaged pairs and an...

Vibrational Spectroscopy of Homo- and Heterochiral Amino Acid Dimers: Conformational Landscapes

The homo- and heterochiral protonated dimers of asparagine with serine and with valine were investigated using infrared multiple-photon dissociation (IRMPD) spectroscopy. Extensive quantum-chemical calculations were used in a three-tiered strategy to screen the conformational spaces of all four dimer species. The resulting binary structures were further grouped into five different types based on their intermolecular binding topologies and subunit configurations. For each dimer species, there are eight to fourteen final conformational geometries within a 10 kJ mol-1 window of the global minimum structure for each species. The comparison between the experimental IRMPD spectra and the simulated harmonic IR features allowed us to clearly identify the types of structures responsible for the observation. The monomeric subunits of the observed homo- and heterochiral dimers are compared to the corresponding protonated/neutral amino acid monomers observed experimentally in previous IRMDP/rotational spectroscopic studies. Possible chirality and kinetic influences on the experimental IRMPD spectra are discussed.

Rovibronic signatures of molecular aggregation in the gas phase: subtle homochirality trends in the dimer, trimer and tetramer of benzyl alcohol

Molecular aggregation is of paramount importance in many chemical processes, including those in living beings. Thus, characterization of the intermolecular interactions is an important step in its understanding. We describe here the aggregation of benzyl alcohol at the molecular level, a process governed by a delicate equilibrium between OH⋯O and OH⋯π hydrogen bonds and dispersive interactions. Using microwave, FTIR, Raman and mass-resolved double-resonance IR/UV spectroscopic techniques, we explored the cluster growth up to the tetramer and found a complex landscape, partly due to the appearance of multiple stereoisomers of very similar stability. Interestingly, a consistently homochiral synchronization of transiently chiral monomer conformers was observed during cluster growth to converge in the tetramer, where the fully homochiral species dominates the potential energy surface. The data on the aggregation of benzyl alcohol also constitute an excellent playground to fine-tune the parameters of the most advanced functionals.

Stacked but not Stuck: Unveiling the Role of π→π* Interactions with the Help of the Benzofuran-Formaldehyde Complex

The 1:1 benzofuran–formaldehyde complex has been chosen as model system for analyzing π→π* interactions in supramolecular organizations involving heteroaromatic rings and carbonyl groups. A joint “rotational spectroscopy–quantum chemistry” strategy unveiled the dominant role of π→π* interactions in tuning the intermolecular interactions of such adduct. The exploration of the intermolecular potential energy surface led to the identification of 14 low‐energy minima, with 4 stacked isomers being more stable than those linked by hydrogen bond or lone‐pair→π interactions. All energy minima are separated by loose transition states, thus suggesting an effective relaxation to the global minimum under the experimental conditions. This expectation has been confirmed by the experimental detection of only one species, which was unambiguously assigned owing to the computation of accurate spectroscopic parameters and the characterization of 11 isotopologues. The large number of isotopic species opened the way to the determination of the first semi‐experimental equilibrium structure for a molecular complex of such a dimension.

Microwave Spectrum of Two-Top Molecule: 2-Acetyl-3-Methylthiophene

The microwave spectrum of 2-acetyl-3-methylthiophene (2A3MT) was recorded in the frequency range from 2 to 26.5 GHz using a molecular jet Fourier transform microwave spectrometer and could be fully assigned to the anti-conformer of the molecule, while the syn-conformer was not observable. Torsional splittings of all rotational transitions in quintets due to internal rotations of the acetyl methyl and the ring methyl groups were resolved and analyzed, yielding barriers to internal rotation of 306.184(46) cm^-1 and 321.813(64) cm^-1 , respectively. The rotational and centrifugal distortion constants were determined with high accuracy, and the experimental values are compared to those derived from quantum chemical calculations. The experimentally determined inertial defect supports the conclusion that anti-2A3MT is planar, even though a number of MP2 calculations predicted the contrary.

Microwave Study of Triflic Acid Hydrates: Evidence for the Transition from Hydrogen-Bonded Clusters to a Microsolvated Ion Pair

Rotational spectra of the mono-, di-, and trihydrates of triflic acid, CF₃SO₃H···(H₂O)_n=1-3, have been recorded by pulsed nozzle Fourier transform microwave spectroscopy and spectroscopic constants obtained have been compared with values calculated at several levels of theory. The experimental results are consistent with the theoretical predictions presented here and elsewhere, indicating that with only one or two water molecules, triflic acid remains un-ionized in a cold molecular complex. The experiments further concur with theoretical predictions that the addition of a third water molecule transforms the system into what is best regarded as a hydrated hydronium triflate ion pair. Thus, only three water molecules are needed to induce ionization of triflic acid in a cold molecular cluster. This number is somewhat low compared with that for other simple protic acids and likely reflects the superacidity of triflic acid. Simple energetic arguments can be used to rationalize this result.

Structures and hydrogen bonding of 1,7-dioxaspiro[5.5]undecane and its hydrates

The conformations of 1,7DSU and its stepwise solvation by up to 5 water molecules were explored using supersonic-jet Fourier transform microwave spectroscopy with the supplement of theoretical calculations. Experimentally, the rotational spectra of the most stable structures of the monomer, monohydrate and dihydrate were observed and assigned. The characteristics of the stability and intermolecular interaction topologies of the 1,7DSU monomer and its hydrated clusters were obtained by CREST conformational sampling followed by B3LYP-D3(BJ)/def2-TZVP geometrical optimizations and MP2/aug-cc-pVTZ single-point energy calculations. The first water molecule links to the 1,7DSU monomer through an O_wH⋯O hydrogen bond. The water molecules tend to aggregate with each other and form cyclic structures for the n = 2-5 clusters. The interactions between water and the 1,7DSU monomer as well as those between water and water were revealed. The analyses of non-covalent interactions and the natural bond orbital suggest that the O_wH⋯O_1,7DSU, O_wH⋯O_w, and CH⋯O_w hydrogen bonds play a prominent role in structural stability.

Development and Validation of a Parameter-Free Model Chemistry for the Computation of Reliable Reaction Rates

A recently developed model chemistry (jun-Cheap) has been slightly modified and proposed as an effective, reliable, and parameter-free scheme for the computation of accurate reaction rates with special reference to astrochemical and atmospheric processes. Benchmarks with different sets of state-of-the-art energy barriers spanning a wide range of values show that, in the absence of strong multireference contributions, the proposed model outperforms the most well-known model chemistries, reaching a subchemical accuracy without any empirical parameter and with affordable computer times. Some test cases show that geometries, energy barriers, zero point energies, and thermal contributions computed at this level can be used in the framework of the master equation approach based on the ab initio transition-state theory for obtaining accurate reaction rates.

Integration of theory, simulation, artificial intelligence and virtual reality: a four-pillar approach for reconciling accuracy and interpretability in computational spectroscopy

The established pillars of computational spectroscopy are theory and computer based simulations. Recently, artificial intelligence and virtual reality are becoming the third and fourth pillars of an integrated strategy for the investigation of complex phenomena. The main goal of the present contribution is the description of some new perspectives for computational spectroscopy, in the framework of a strategy in which computational methodologies at the state of the art, high-performance computing, artificial intelligence and virtual reality tools are integrated with the aim of improving research throughput and achieving goals otherwise not possible. Some of the key tools (e.g., continuous molecular perception model and virtual multifrequency spectrometer) and theoretical developments (e.g., non-periodic boundaries, joint variational-perturbative models) are shortly sketched and their application illustrated by means of representative case studies taken from recent work by the authors. Some of the results presented are already well beyond the state of the art in the field of computational spectroscopy, thereby also providing a proof of concept for other research fields.

Rotational Spectrum and Molecular Structures of the Binary Aggregates of 1,1,1,3,3,3-Hexafluoro-2-propanol with Ne and Ar

The structures and binding topologies of two binary van der Waals complexes 1,1,1,3,3,3-hexafluoro-2-propanol (HFIP)···Ne and ···Ar were investigated. The rotational spectra of these two complexes including several isotopic species containing ²⁰Ne, ²²Ne, ⁴⁰Ar, ¹³C, and hydroxyl D were measured using a chirped pulse Fourier transform microwave spectrometer and a cavity-based Fourier transform microwave spectrometer. While HFIP was shown to exist in both the gauche and trans configurations based on previous reports, the rare gas atom is predicted to attach to HFIP in several different binding topologies, leading to a total of nine possible structural isomers for each complex. Only one isomer was detected for each species, and it corresponds to the most stable one predicted, based on the comparison of the experimental rotational constants and electric dipole moment components with the theoretical predictions and on the isotopic data. We applied quantum theory of atoms in molecules (QTAIM) and electrostatic potential calculations to examine the different rare gas binding sites and to explore the nature of the interactions in these two complexes and several previously reported alcohol···Ar complexes. The effects of fluorination are also discussed by comparison with the binary complexes of isopropanol···Ne and ···Ar.

Dispersion forces in chirality recognition - a density functional and wave function theory study of diols

In the discussion of chirality recognition, steric considerations and strongly directed interactions such as hydrogen bonds are primarily discussed. However, given the sheer size of biomolecules, it is expected that dispersion forces could also play a determining role for aggregate formation and associated chirality recognition. With the example of diol molecules, we explore different factors in the formation of homo- and hetero-dimers as well as their relative stability. By comparing density functional results with the analysis of local correlation methods, we infer the impact of dispersion not only on the energies but also on the structures of such chiral aggregates. A local orbital based scheme is used to calculate wave function dispersion-free gradients and compare to uncorrected density functional structures.

Matrix Isolation-Vibrational Circular Dichroism Spectroscopic Study of Conformations and Non-Covalent Interactions of Tetrahydro-2-Furoic Acid

The conformational landscape and aggregation behaviour of tetrahydro-2-furoic acid (THFA) were investigated by using matrix isolation-vibrational circular dichroism (MI-VCD). The well-resolved experimental MI-IR and MI-VCD features in an argon matrix at 10 K allow one to identify two dominant monomeric conformations as trans-THFA where the hydroxyl and carbonyl groups of COOH are at opposite sides, as well as one cis-conformer. At 24 K and 30 K deposition temperatures, the experimental IR and VCD spectral features reveal further growth of the binary THFA aggregates. Systematic conformational searches identified three vastly different binary binding topologies, resulting in a few hundred stable (THFA)₂ conformers. Interestingly, the main binary structures observed correspond to an unusual type of structure which is made of two trans-THFA subunits, in contrast to the usual double H-bonded ring binary structures, identified in a previous solution study. The present work showcases the power of MI-VCD spectroscopy in revealing unusual structures formed in a cold rare gas matrix.

Unveiling Bifunctional Hydrogen Bonding with the Help of Quantum Chemistry: The Imidazole-Water Adduct as Test Case

The ubiquitous role of water and its amphiprotic nature call for a deeper insight into the physical-chemical properties of hydrogen-bonded complexes formed with building blocks of biomolecules. In this work, the semiexperimental (SE) approach combined with the template model (TM) protocol allowed the accurate determination of the equilibrium structure of two isomeric forms of the imidazole-water complex. In this procedure, the integration of experiment (thanks to a recent rotational spectroscopy investigation) and theory is exploited, also providing the means of assessing the reliability and accuracy of different quantum-chemical approaches. Overall, this study demonstrated the robustness of the combined SE-TM approach, which can provide accurate results using affordable quantum-chemical methods. Finally, the structural and energetic characteristics of these complexes have been examined in detail and compared with those of analogous heterocycle-water adducts, also exploiting energy decomposition analyses.

Dynamic chiral self-recognition in aromatic dimers of styrene oxide revealed by rotational spectroscopy

Chiral molecular recognition is a pivotal phenomenon in biomolecular science, governed by subtle balances of intermolecular forces that are difficult to quantify. Non-covalent interactions involving aromatic moieties are particularly important in this realm, as recurring motifs in biomolecular aggregation. In this work, we use high-resolution broadband rotational spectroscopy to probe the dynamic conformational landscape enclosing the self-pairing topologies of styrene oxide, a chiral aromatic system. We reach a definite assignment of four homochiral and two heterochiral dimers using auxiliary quantum chemistry calculations as well as structure-solving methods based on experimental isotopic information. A complete picture of the dimer conformational space is obtained, and plausible routes for conformational relaxation are derived. Molecular structures are discussed in terms of conformational flexibility, the concerted effort of weak intermolecular interactions, and their role in the expression of the molecular fit.

Conformational Landscape, Chirality Recognition and Chiral Analyses: Rotational Spectroscopy of Tetrahydro-2-Furoic Acid⋅⋅⋅Propylene Oxide Conformers

A chiral adduct formed between a chiral carboxylic acid, tetrahydro-2-furoic acid (THFA), and a chiral ester, propylene oxide (PO), was investigated using rotational spectroscopy and DFT calculations. Isolated THFA exists dominantly as three different conformers: I, II, and III in a jet, with I and II taking on the trans-COOH configuration and III having the cis-COOH configuration. We utilized CREST, a conformational ensemble space exploration tool, to identify the possible conformations of the binary adduct, THFA⋅⋅⋅PO. Subsequent DFT geometry optimizations predicted about two hundred homochiral and heterochiral binary structures with 28 low energy structures within an energy window of 15 kJ mol^-1 . A rich broadband rotational spectrum was obtained with a mixture of trace amounts of THFA+PO in neon in a supersonic jet expansion. Six THFA⋅⋅⋅PO conformers were identified experimentally. Kinetically favored binary products which contain trans-COOH I dominate among the observed conformers, while thermodynamically more stable adducts were also detected. Detailed analyses of the structures of the observed conformers show interesting chirality-controlled structural preferences. Such non-covalently bound chiral contact pairs are the foundation of chiral-tag rotational spectroscopy, an exciting new analytical application of rotational spectroscopy for determination of enantiomeric excess. Enantiomeric excess analyses were performed and the results are discussed.

Modifying conformational distribution of chiral tetrahydro-2-furoic acid through its interaction with water: a rotational spectroscopic and theoretical investigation

Rotational spectrum of a binary complex formed between tetrahydro-2-furoic acid (THFA) and water was measured using a chirped pulse Fourier transform microwave spectrometer. A comprehensive theoretical conformational search procedure was carried out using CREST, a conformational searching tool, and DFT calculations to aid the spectral assignment and interpretation. The final conformer ensemble is classified into two structural groups: Type 1 conformers showing a classic carboxylic acid monohydrate structure with two strong hydrogen-bonds formed between the COOH group of cis-THFA and water, and the much less stable Type 2 conformers with trans-THFA and weaker intermolecular interactions with water. The 'cis-' and 'trans-' labels refer to the configurations where the carboxylic C[double bond, length as m-dash]O and OH functional groups are on the same or opposite side, respectively. Only the two most stable Type 2 conformers containing trans-THFA I and II were observed experimentally in a neon jet expansion with an abundance ratio of 1 : 1. This relative abundance observation differs greatly from that of the THFA monomer, i.e. with trans-THFA I : trans-THFA II : cis-THFA III of 10 : 1 : 1 in a neon jet expansion, reported previously. The observation indicates a kinetically controlled formation process of different types of the monohydrates in a jet expansion, whereas a thermodynamically controlled process dominates within each type of structures. The relative stability of the THFA ring conformations is altered by interaction with water, showing a noticeable water induced conformational preference.

Discovering the Elusive Global Minimum in a Ternary Chiral Cluster: Rotational Spectra of Propylene Oxide Trimer

The chirality controlled conformational landscape of the trimer of propylene oxide (PO), a prototypical chiral molecule, was investigated using rotational spectroscopy and a range of theoretical tools for conformational searches and for evaluating vibrational contributions to effective structures. Two sets of homochiral (PO)3 rotational transitions were assigned and the associated conformers identified with theoretical support. One set of heterochiral (PO)3 transitions was assigned, but no structures generated by one of the latest, advanced conformational search codes could account for them. With the aid of a Python program, the carbon atom backbone and then the heterochiral (PO)3 structure were generated using 13 C isotopic data. Excellent agreement between theoretical and experimental rotational constants and relative dipole moment components of all three conformers was achieved, especially after applying vibrational corrections to the rotational constants.