Conformer-Selective Photoelectron Circular Dichroism

Conformational flexibility and chirality both play a key role in molecular recognition. It is therefore very useful to develop spectroscopic methods that simultaneously probe both properties. It has been theoretically predicted that photoelectron circular dichroism (PECD) should be very sensitive to conformational isomerism. However, experimental proof has been less forthcoming and only exists for a very few favorable cases. Here, we present a new PECD scheme based on resonance-enhanced two-photon ionization (RE2PI) using UV/Vis nanosecond laser excitations. The spectral resolution obtained thereby guarantees conformer-selectivity by inducing resonant conformer-specific ππ* S1←S0 transitions. We apply this experimental scheme to the study of chiral 1-indanol, which exists in two conformers linked by a ring inversion and defined by the position of the hydroxyl group, namely axial and equatorial. We show that the PECD of the equatorial and axial forms considerably differ in sign, magnitude and shape. We also discuss the influence of the total ionization energy, vibronic excitation of intermediate and final states, and relative polarization of the excitation and ionization lasers. Conformer-specificity adds a new dimension to the applications of PECD in analytical chemistry addressing now the general case of floppy systems.

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

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

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.

Modelling solute-solvent interactions in VCD spectra analysis with the micro-solvation approach

Vibrational circular dichroism (VCD) spectroscopy has become an important part of the (stereo-)chemists' toolbox as a reliable method for the determination of absolute configurations. Being the chiroptical version of infrared spectroscopy, it has also been recognized as being very sensitive to conformational changes and intermolecular interactions. This sensitivity originates from the fact that the VCD spectra of individual conformers are often more different than their IR spectra, so that changes in conformational distributions or band positions and intensities become more pronounced. What is an advantage for studies focussing on intermolecular interactions can, however, quickly turn into a major obstacle during AC determinations: solute-solvent interactions can have a strong influence on spectral signatures and they must be accurately treated when simulating VCD and IR spectra. In this perspective, we showcase selected examples which exhibit particularly pronounced solvent effects. It is demonstrated that it is typically sufficient to model solute-solvent interactions by placing single solvent molecules near hydrogen bonding sites of the solute and subsequently use the optimized structures for spectra simulations. This micro-solvation approach works reasonably well for medium-sized, not too conformationally flexible molecules. We thus also discuss its limitations and outline the next steps that method development needs to take in order to further improve the workflows for VCD spectra predictions.

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.

Vibrational Circular Dichroism Spectroscopy of Chiral Molecular Crystals: Insights from Theory

Chirality is a curious phenomenon that appears in various forms. While the concept of molecular (RS-)chirality is ubiquitous in chemistry, there are also more intricate forms of structural chirality. One of them is the enantiomorphism of crystals, especially molecular crystals, that describes the lack of mirror symmetry in the unit cell. Its relation to molecular chirality is not obvious, but still an open question, which can be addressed with chiroptical tools. Vibrational circular dichroism (VCD) denotes chiral infrared (IR) spectroscopy that is susceptible to both, the molecular as well as the intermolecular space by means of vibrational transitions. When carried out in the solid state, VCD delivers a very rich set of non-local contributions that are determined by crystal packing and collective motion. Since its discovery in the 1970s, VCD has become the method of choice for the determination of absolute configurations, but its applicability reaches beyond towards the study of different crystal forms and polymorphism. This brief review summarises the theoretical concepts of crystal chirality and how computations of solid-state VCD can shed light into the intimate connection of chiral structure and vibrational optical activity.

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.

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.

Isolated and solvated chiroptical behavior in conformationally flexible butanamines

The nonresonant optical activity of two highly flexible aliphatic amines, (2R)-3-methyl-2-butanamine (R-MBA) and (2R)-(3,3)-dimethyl-2-butanamine (R-DMBA), has been probed under isolated and solvated conditions to examine the roles of conformational isomerism and to explore the influence of extrinsic perturbations. The optical rotatory dispersion (ORD) measured in six solvents presented uniformly negative rotatory powers over the 320-590 nm region, with the long-wavelength magnitude of chiroptical response growing nearly monotonically as the dielectric constant of the surroundings diminished. The intrinsic specific optical rotation,α λ T (in deg dm-1 [g/mL]-1 ), extracted for ambient vapor-phase samples of R-MBA [-11.031(98) and -2.29 (11)] and R-DMBA [-9.434 (72) and -1.350 (48)] at 355 and 633 nm were best reproduced by counterintuitive solvents of high polarity (yet low polarizability) like acetonitrile and methanol. Attempts to interpret observed spectral signatures quantitatively relied on the linear-response frameworks of density-functional theory (B3LYP, cam-B3LYP, and dispersion-corrected analogs) and coupled-cluster theory (CCSD), with variants of the polarizable continuum model (PCM) deployed to account for the effects of implicit solvation. Building on the identification of several low-lying equilibrium geometries (nine for R-MBA and three for R-DMBA), ensemble-averaged ORD profiles were calculated at T = 300 K by means of the independent-conformer ansatz, which enabled response properties predicted for the optimized structure of each isomer to be combined through Boltzmann-weighted population fractions derived from corresponding relative internal-energy or free-energy values, the latter of which stemmed from composite CBS-APNO and G4 analyses. Although reasonable accord between theory and experiment was realized for the isolated (vapor-phase) species, the solution-phase results were less satisfactory and tended to degrade progressively as the solvent polarity increased. These trends were attributed to solvent-mediated changes in structural parameters and energy metrics for the transition states that separate and putatively isolate the equilibrium conformations supported by the ground electronic potential-energy surface, with the resulting displacement of barrier locations and/or decrease of barrier heights compromising the underlying premise of the independent-conformer ansatz.

Influence of the environment on the infrared spectrum of alanine: An effective mode analysis

The vibrational spectrum of the alanine amino acid was computationally determined in the infrared range 1000-2000 cm^-1, under various environments encompassing the gas, hydrated, and crystalline phases, by means of classical molecular dynamics trajectories, carried out with the Atomic Multipole Optimized Energetics for Biomolecular Simulation polarizable force field. An effective mode analysis was performed, in which the spectra are optimally decomposed into different absorption bands arising from well-defined internal modes. In the gas phase, this analysis allows us to unravel the significant differences between the spectra obtained for the neutral and zwitterionic forms of alanine. In condensed phases, the method provides invaluable insight into the molecular origins of the vibrational bands and further shows that peaks with similar positions can be traced to rather different molecular motions.

Quantum Mechanical-Cluster Approach to Solve the Bioisosteric Replacement Problem in Drug Design

Bioisosteres are molecules that differ in substituents but still have very similar shapes. Bioisosteric replacements are ubiquitous in modern drug design, where they are used to alter metabolism, change bioavailability, or modify activity of the lead compound. Prediction of relative affinities of bioisosteres with computational methods is a long-standing task; however, the very shape closeness makes bioisosteric substitutions almost intractable for computational methods, which use standard force fields. Here, we design a quantum mechanical (QM)-cluster approach based on the GFN2-xTB semi-empirical quantum-chemical method and apply it to a set of H → F bioisosteric replacements. The proposed methodology enables advanced prediction of biological activity change upon bioisosteric substitution of -H with -F, with the standard deviation of 0.60 kcal/mol, surpassing the ChemPLP scoring function (0.83 kcal/mol), and making QM-based ΔΔG estimation comparable to ∼0.42 kcal/mol standard deviation of in vitro experiment. The speed of the method and lack of tunable parameters makes it affordable in current drug research.

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.

Alternating 1-Phenyl-2,2,2-Trifluoroethanol Conformational Landscape With the Addition of One Water: Conformations and Large Amplitude Motions

The 1:1 adduct of 1-phenyl-2,2,2-trifluoroethanol (PhTFE), a chiral fluoroalcohol, with water was investigated using chirped pulse Fourier transform microwave spectroscopy and computational methods. While PhTFE itself was predicted to have three minima, I (gauche+), II (trans), and III (gauche-), only I and II were stable and only I was observed experimentally. A systematic search of the PhTFE···H₂O conformational landscape identified 110 stable minima, 14 of which are within a 15 kJ mol^-1 energy window. Rotational spectra of the two PhTFE···H₂O conformers along with several deuterium and ¹⁸O isotopologues were assigned, and the isotopic data were used to verify the corresponding structures. In the two observed monohydrate conformers, one contains PhTFE I where the water subunit is inserted into the existing intramolecular OH···F contact of I, and the binary adduct is stabilized by two intermolecular contacts: OH···O_W and H_W···F, whereas the other contains PhTFE II where the water subunit interacts with both the alcohol hydrogen and phenyl ring of II, demonstrating that interaction with water sufficiently stabilizes II for its observation in a jet expansion. Interestingly, the predicted electric dipole moment components at the identified minima deviate considerably from the experimental ones. Such deviations were analyzed in terms of dynamic effects associated with the large amplitude motions of the unbound H_W. In addition, tunnelling effects associated with the exchange of the bonded and nonbonded H_W were also discussed.

Dynamic Stereochemistry of a Biphenyl‐Bisprolineamide Model Catalyst and its Imidazolidinone Intermediates

In this study, we characterize the dynamic stereochemistry of a biphenyl‐2,2’‐bis(proline amide) catalyst in chloroform and DMSO as representative weakly and strongly hydrogen bonding solvents. Using vibrational circular dichroism (VCD) spectroscopy and density functional theory (DFT) based spectra calculations, we show that the preferred axial stereochemistry of the catalyst is determined by solute‐solvent interactions. Explicitly considering solvation with DMSO molecules is found to be essential to correctly predict the conformational preferences of the catalyst. Furthermore, we investigate the stereochemistry of the corresponding enamines and imidazolidinones that are formed upon reaction with isovaleraldehyde. The enamines are found to rapidly convert to endo‐imidazolidinones and the thermodynamically favored exo‐imidazolidinones are formed only slowly. The present study demonstrates that the stereochemistry of these imidazolidinones can be deduced directly from the VCD spectra analysis without any further detailed analysis of NMR spectra. Hence, we herein exemplify the use of VCD spectroscopy for an in situ characterization of intermediates relevant in asymmetric catalysts.

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.

VCD spectroscopy reveals conformational changes of chiral crown ethers upon complexation of potassium and ammonium cations

Two chiral derivatives of 18-crown-6, namely the host molecules 2,3-diphenyl- and 2-phenyl-18c6, serve as model systems to investigate whether VCD spectroscopy can be used to monitor conformational changes occurring upon complexation of guests. Host-guest complexes of both crown ethers were prepared by addition of KNO₃. The more bulky 2,3-diphenyl-18c6 is found to undergo major conformational changes upon encapsulation of K⁺, which are revealed as characteristic changes of the VCD spectral signatures. In contrast, while 2-phenyl-18c6 also incorporates K⁺ into the macrocycle, strong conformational changes are not occurring and thus spectral changes are negligible. With an octyl ammonium cation as guest molecule, 2,3-diphenyl-18c6 shows the same conformational and spectral changes that were observed for K⁺-complexes. In addition, the asymmetric NH₃-deformation modes are found to gain VCD intensity through an induced VCD process. An analysis of the vibrational spectra enables a differentiation of VCD active and inactive guest modes: There appears to be a correlation between the symmetry of the vibrational mode and the induced VCD intensity. While this finding makes the host-guest complexes interesting systems for future theoretical studies on the origin of induced VCD signatures, the observations described in this study demonstrate that VCD spectroscopy is indeed a suitable technique for the characterization of supramolecular host-guest complexes.

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.

Photoelectron Circular Dichroism as a Signature of Subtle Conformational Changes: The Case of Ring Inversion in 1-Indanol

Chirality plays a fundamental role in the molecular recognition processes. Molecular flexibility is also crucial in molecular recognition, allowing the interacting molecules to adjust their structures and hence optimize the interaction. Methods probing simultaneously chirality and molecular conformation are therefore crucially needed. Taking advantage of a possible control in the gas phase of the conformational distribution between the equatorial and axial conformers resulting from a ring inversion in jet-cooled 1-indanol, we demonstrate here the sensitivity of valence-shell photoelectron circular dichroism (PECD) to both chirality and subtle conformational changes, in a case where the photoelectron spectra of the two conformers are identical. For the highest occupied orbital, we observe a dramatic inversion of the PECD-induced photoelectron asymmetries, while the photoionization cross-section and usual anisotropy (β) parameter are completely insensitive to conformational isomerism. Such a sensitivity is a major asset for the ongoing developments of PECD-based techniques as a sensitive chiral (bio)chemical analytical tool in the gas phase.

The solid state VCD of a novel N-acylhydrazone trifluoroacetate

A novel N-acylhydrazone with pharmaceutical importance was subject of structural and IR/VCD investigations in the solid state. In the crystal structure, dimers of anion-cation pairs are stabilized by H-bonding and ionic interactions. Some less common interaction types, like C=N···C-NH₃⁺ (σ-hole) interactions, hydrazone-aromatic interactions and dispersive contacts of the CF₃ groups are also present in the crystal. Satisfactory reproduction of the solid state IR and VCD spectra required that quantum-chemical calculations be done on a tetramer (four cation-anion pairs) cut out from the crystal structure, exhibiting key intermolecular interactions. Ten DFT functionals were assessed as to the agreement between the calculated and experimental spectra. Various approaches to scaling of the calculated frequencies were applied. The best results were yielded with individual (optimized) frequency scaling factors (FSFs) and band half-widths at half maximum-(HWHM) for four separate spectral subregions. The best matching between the experimental and theoretical spectra (according to SimIR, SimVCD and SimVDF indices) was found for the B3PW91 functional, however, a few other functionals follow closely in the ranking. Based on the quantum chemical calculations, spectral assignments have been made.

Chemistry Dissolved in Ionic Liquids. A Theoretical Perspective

The theoretical treatment of ionic liquids must focus now on more realistic models while at the same time keeping an accurate methodology when following recent ionic liquids research trends or allowing predictability to come to the foreground. In this Perspective, we summarize in three cases of advanced ionic liquid research what methodological progress has been made and point out difficulties that need to be overcome. As particular examples to discuss we choose reactions, chirality, and radicals in ionic liquids. All these topics have in common that an explicit or accurate treatment of the electronic structure and/or intermolecular interactions is required (accurate methodology), while at the same time system size and complexity as well as simulation time (realistic model) play an important role and must be covered as well.

VCD spectroscopy distinguishes the enamine and iminium ion of a 1,1’-binaphthyl azepine

We present a VCD spectroscopic characterization of a chiral 1,1’-binaphthyl azepine catalyst and show that the VCD spectra of an in-situ generated enamine and an ex-situ prepared iminium ion are...

Competition between inter and intramolecular hydrogen bond evidenced by vibrational circular dichroism spectroscopy: The case of (1S,2R)-(-)-cis-1-amino-2-indanol

The infrared (IR) absorption and vibrational circular dichroism (VCD) spectra of an intramolecularly hydrogen-bonded chiral amino-alcohol, (1S,2R)-(-)-cis-1-amino-2-indanol, are studied in DMSO-d₆ . The spectra are simulated at the density functional theory (DFT) level within the frame of the cluster-in-the-liquid model. Both IR and VCD spectra show a clear signature of the formation of intermolecular hydrogen bonds at the detriment of the intramolecular OH … N interaction present in the isolated molecule. Two solvent molecules are necessary to reproduce the experimental spectra. Whereas the first DMSO molecule captures the main spectral modifications due to hydrogen bond formation between the solute and the solvent, the second DMSO molecule is necessary for a good description of the Boltzmann contribution of the different complexes, based on their Gibbs free energy.

The important role of non-covalent interactions for the vibrational circular dichroism of lactic acid in aqueous solution

We present a computational study of vibrational circular dichroism (VCD) in solutions of (S)-lactic acid, relying on ab initio molecular dynamics (AIMD) and full solvation with bulk water. We discuss the effect of the hydrogen bond network on the aggregation behaviour of the acid: while aggregates of the solute represent conditions encountered in a weakly interacting solvent, the presence of water drastically interferes with the clusters - more strongly than originally anticipated. For both scenarios we computed the VCD spectra by means of nuclear velocity perturbation theory (NVPT). The comparison with experimental data allows us to establish a VCD-structure relationship that includes the solvent network around the chiral solute. We suggest that fundamental modes with strong polarisation such as the carbonyl stretching vibration can borrow VCD from the chirally restructured solvent cage, which extends the common explanatory models of VCD generation in aqueous solution.

Computation of Solid-State Vibrational Circular Dichroism in the Periodic Gauge

We introduce a new theoretical formalism to compute solid-state vibrational circular dichroism (VCD) spectra from molecular dynamics simulations. Having solved the origin-dependence problem of the periodic magnetic gauge, we present IR and VCD spectra of (1S,2S)-trans-1,2-cyclohexanediol obtained from first-principles molecular dynamics calculations and nuclear velocity perturbation theory, along with the experimental results. Because the structure model imposes periodic boundary conditions, the common origin of the rotational strength has hitherto been ill-defined and was approximated by means of averaging multiple origins. The new formalism reconnects the periodic model with the finite physical system and restores gauge freedom. It nevertheless fully accounts for nonlocal spatial couplings from the gauge transport term. We show that even for small simulation cells the rich nature of solid-state VCD spectra found in experiments can be reproduced to a very satisfactory level.

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.

VCD spectroscopy reveals that a water molecule determines the conformation of azithromycin in solution

We report the IR and VCD spectra of azithromycin, a macrolide antibiotic with a total of 18 stereogenic centers. The computational analysis of the spectra reveals that a single water molecule has to be considered in the conformational search. Its key role is the stabilization of an extended hydrogen bonding network and an otherwise unstable conformation that determines the VCD spectral signatures.