Vibrational circular dichroism from ab initio molecular dynamics and nuclear velocity perturbation theory in the liquid phase

Vibrational circular dichroism spectra of proline in water at different pH values

Recording VCD spectra of aqueous solution poses a particular challenge as water is a strong infrared absorber. Likewise, the computational analysis of VCD spectra by means of DFT-based spectral calculations requires the consideration of explicit solvent molecules, thus posing an even greater challenge. Several studies suggested that by modeling the solvent environment with a few water molecules in a micro-solvation approach would be sufficient to describe experimental spectra. For example, using proline at different pH values, we herein show that a change in the relative spatial orientation of a single water molecule in five-fold solvated structures strongly affects the computed VCD spectral signatures and that Boltzmann-weighted spectra do not correctly reproduce the experiment. We thus explored an approach based on molecular dynamics and subsequent DFT-calculations, in which we considered 30 water molecules (about 1.5 solvation shells). Once again, it was found that the Boltzmann-weighted spectra obtained on the basis of several hundred structures did not correctly reproduce experimental signatures, and a simple averaging scheme resulted in well-matching spectra with comparable bandwidths. The rationale behind the procedure was that sampling the configurational space of the solvent molecules is as equally important as the conformational sampling of the solute. For conformationally more flexible molecules, it is assumed that a much larger set of structures will have to be computed in order to properly sample the conformational space of both solute and solvent.

Chiral Spectroscopy of Bulk Systems with Propagated Localized Orbitals

We present approaches for the simulation of electronic circular dichroism, Raman, and Raman optical activity (ROA) spectra for isolated and periodic systems as well as subsystem analysis thereof. The method is based on the use of time-dependent maximally localized Wannier functions in the CP2K package and accounts for origin dependencies inherent to the Gaussian and plane wave with pseudopotentials approach as well as the origin dependence of the magnetic dipole and electric quadrupole operators. Tests on the H-bonded enantiomers of alanine by harmonic normal-mode analysis and on an aqueous solution of l-alanine by ab initio molecular dynamics obeying periodic boundary conditions (PBCs) are presented as total and subsystem-resolved spectra. To our knowledge, this is the first instance of an ROA spectrum derived from real-time propagation obeying PBCs and the first ROA simulation considering off-, pre-, and on-resonance effects within PBCs.

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.

Practical phase-space electronic Hamiltonians for ab initio dynamics

Modern electronic structure theory is built around the Born-Oppenheimer approximation and the construction of an electronic Hamiltonian Ĥel(X) that depends on the nuclear position X (and not the nuclear momentum P). In this article, using the well-known theory of electron translation (Γ') and rotational (Γ″) factors to couple electronic transitions to nuclear motion, we construct a practical phase-space electronic Hamiltonian that depends on both nuclear position and momentum, ĤPS(X,P). While classical Born-Oppenheimer dynamics that run along the eigensurfaces of the operator Ĥel(X) can recover many nuclear properties correctly, we present some evidence that motion along the eigensurfaces of ĤPS(X,P) can better capture both nuclear and electronic properties (including the elusive electronic momentum studied by Nafie). Moreover, only the latter (as opposed to the former) conserves the total linear and angular momentum in general.

Influence of microsolvation on vibrational circular dichroism spectra in dimethyl sulfoxide solvent: A Bottom-Up approach using Quantum cluster growth

Chiroptical spectroscopic measurements serve as routine methods to assign the absolute configuration of chiral compounds and interpret their conformational behavior in solution. One common challenge is the use of strongly hydrogen-bonding solvents, which can significantly bias the conformational ensemble and affect the vibrational circular dichroism (VCD) active bands in solution. One such solvent is dimethyl sulfoxide (DMSO)-an excellent solvent for stubborn compounds-that must be explicitly considered in VCD analysis. Explicit consideration of solvent remains a critical challenge in chiroptical spectroscopy due to the need to explore solute-solvent conformational space and the computational expense in modeling these clusters. Interested in the recent development of the Quantum Cluster Growth (QCG) program by the Grimme lab, we set out to model and interpret previously reported VCD spectra for several molecules using their efficient program. Our purposes are two-fold: (1) to investigate the applicability of the QCG program to the problem of reproducing VCD spectra in DMSO solvent and (2) to identify limitations in using this approach. We find that we can conveniently model and analyze the VCD spectra of investigated molecules in DMSO. However, the final set of conformers used for VCD calculations are functional dependent and different sets of conformers can provide satisfactory quantitative agreement between experimental and predicted VCD spectra. We hope that this study provides guidance for future chiroptical studies in the challenging DMSO solvent.

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.

Selective Chirality Transfer to the Bis(trifluoromethylsulfonyl)imide Anion of an Ionic Liquid

Chirality transfer from the chiral molecule (R)-1,2-propylene oxide to the achiral anion of the 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ionic liquid is observed. The chiral probe selectively affects one part of the binary ionic liquid, i. e., it has previously been shown experimentally and theoretically that this particular imidazolium cation can be affected by chirality transfer, but in the present system chirality is almost exclusively transferred to the anion and not to both parts of the solvent (anion and cation). This observation is of high relevance because of its selectivity and because anion effects are usually much more important in ionic liquid research than cation effects. From ab initio molecular dynamics simulations, a conformational analysis and dissected vibrational circular dichroism spectra are obtained to study the chirality transfer. While in the neat ionic liquid two mirror imaged trans conformers of the anion occur almost equally, we observe an excess of one of these conformers in the presence of the chiral solute, causing optical activity of the anion. Although the cis conformers are not tremendously affected by the chirality transfer, they gain in total population when (R)-1,2-propylene oxide is dissolved in the ionic liquid.

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.

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.

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.

Exact Factorization Adventures: A Promising Approach for Non-Bound States

Modeling the dynamics of non-bound states in molecules requires an accurate description of how electronic motion affects nuclear motion and vice-versa. The exact factorization (XF) approach offers a unique perspective, in that it provides potentials that act on the nuclear subsystem or electronic subsystem, which contain the effects of the coupling to the other subsystem in an exact way. We briefly review the various applications of the XF idea in different realms, and how features of these potentials aid in the interpretation of two different laser-driven dissociation mechanisms. We present a detailed study of the different ways the coupling terms in recently-developed XF-based mixed quantum-classical approximations are evaluated, where either truly coupled trajectories, or auxiliary trajectories that mimic the coupling are used, and discuss their effect in both a surface-hopping framework as well as the rigorously-derived coupled-trajectory mixed quantum-classical approach.

Impact on the Raman spectra of liquids when a polarized light source is used

The polarization state of the excitation light used in two Raman systems was controlled to study its effect in the unpolarized Raman spectra of unstructured samples. Both systems work in different regions of the electromagnetic spectrum (NIR and visible). Four polarization states (linear, linear at 45° and 90°, and circular) were used to excite liquid samples (ethanol, acetone, and their mixture). The results show that the Raman peaks intensities' ratio varies according to the polarization state of the excitation light. Peaks related to functional groups and C-H stretching modes increase their intensity when circular polarization (CP) is applied. The latter may help to study liquid mixtures with low concentrations. Different polarizing light states give a more detailed spectroscopic analysis since it gathers more structural information of the samples tested in this work with an undefined structure.

Implementation of Nuclear Velocity Perturbation and Magnetic Field Perturbation Theory in CP2K and Their Application to Vibrational Circular Dichroism

We present the implementation of nuclear velocity perturbation theory (NVPT), using a pioneering combination of atom-centered (velocity-dependent) Gaussian basis functions and plane waves in the CP2K package. The atomic polar tensors (APTs) and atomic axial tensors (AATs) are evaluated in the velocity representation using efficient density functional perturbation theory. The presence of nonlocal pseudopotentials, the representation of potentials on numerical integration grids, and effects arising from the basis functions being centered on the atoms have been considered in the implementation. The Magnetic Field Perturbation Theory (MFPT) using gauge-including atomic orbitals is implemented in the same code and compared to the NVPT. Our implementation is the first to compare both approaches (MFPT and NVPT) in the same code. The implementation has been verified via sum rules and by investigating the gauge origin dependence of the AATs for a set of small molecules, oxirane, and fluoro-oxirane. We also present vibrational circular dichroism spectra that are related to the APTs and AATs, applying both theories.

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.

Solvation of the Boc-Val-Phe-nPr peptide characterized by VCD spectroscopy and DFT calculations

The conformational preferences of peptides are strongly determined by hydrogen bonding interactions. Intermolecular solute-solvent interactions compete with intramolecular interactions, which typically stabilize the secondary structure of the peptide. The analysis of vibrational circular dichroism (VCD) spectra can give insights into solvation-induced changes in the conformational distribution of small peptides. Here we describe the VCD spectroscopic characterization of the model peptide Boc-Val-Phe-nPr in chloroform as representative for a weakly interacting solvent and dimethyl sulfoxide (DMSO-d₆) as a strongly hydrogen bonding solvent. We show that the conformational preferences of the peptide in chloroform are well-described by the computationally predicted distribution of the isolated molecule assuming only implicit solvation effects through a continuum solvation model. In order to simulate the spectra recorded in DMSO-d₆, solvation was accounted for explicitly by computed microsolvated structures containing one to three solvent molecules. A good match of the computed spectra with the experimental data is obtained by this method. Comparing the conformational distributions in deuterated chloroform-d₁ and DMSO-d₆, structures with intramolecular hydrogen bonds such as the (δ,δ)-conformer family contribute to the conformational distribution only when there is no strong interaction with the solvent. This is in contrast to the results for the related Boc-Pro-Phe-nPr studied before, for which the intramolecular interaction was found to persist in DMSO-d₆. Furthermore, we discuss the influence of hydrogen bonding to different numbers of solvent molecules on the spectral signatures and show that the structure of the peptide in DMSO-d₆ is best described as a mixture of twofold-solvated (δ,β)- and threefold-solvated (β,β)-conformers.

Classical Trajectory of Molecules in Electromagnetic Field: A Handy Method to Simulate Molecular Vibrational Spectra

Within harmonic approximations, molecular vibrational spectra are simulated in a standard way through force field diagonalization and following transformation of Cartesian to normal-mode tensor derivatives. This may become tedious for large systems of many thousands of atoms and also not necessary because of a limited resolution required to interpret an experiment. We developed an algorithm based on the real-time real-field molecular dynamics, effectively at zero temperature, invoked in a molecule by the electromagnetic field of light. The algorithm is simple to implement and suitable for parallel computing, and it can be potentially extended to more complicated molecular-light interaction modes. It circumvents the diagonalization and is suitable to model vibrational optical activity (vibrational circular dichroism and, to a lesser extent, Raman optical activity). For large molecules, it becomes faster than diagonalization, but it also enables the assignment of vibrational spectral bands to local molecular motions.

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.

Chiral Molecular Structure Determination for a Desired Compound Just from Its Molecular Formula and Vibrational Optical Activity Spectra

A novel proof-of-concept model for chiral molecular structure determination using just the molecular formula and vibrational optical activity (VOA) spectra is presented. To verify this concept, the molecular formula of a desired compound is used to generate all possible chiral structural isomers and their VOA spectra are predicted. The similarity analyses of predicted VOA spectra were then carried out in two different ways: (a) similarity between VOA spectrum of one structural isomer with those of the rest, referred to as cross-correlations; (b) similarity between VOA spectra of all chiral structural isomers with the experimental VOA spectra of the desired compound. Three different molecular formulae, C₄H₈O, C₃H₅ClO, and C₆H₁₀O, and their chiral structural isomers (6, 9, and 75 respectively), were investigated. In each case, the correct chiral molecular structure of the desired compound was identified without ambiguity. Cross-correlation analysis revealed the uniqueness of VOA spectra in deducing the chiral molecular structure solely from its molecular formula. Different chiral structural isomers associated with the molecular formula CH₃NO₂ were also found to have no significant cross-correlations between their VOA spectra, opening a pathway to detect and identify the elusive chiral N-hydroxyoxaziridine from its VOA spectra.

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.

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.

Theory and algorithms for chiroptical properties and spectroscopies of aqueous systems

Chiroptical properties and spectroscopies are valuable tools to study chiral molecules and assign absolute configurations. The spectra that result from chiroptical measurements may be very rich and complex, and hide much of their information content. For this reason, the interplay between experiments and calculations is especially useful, provided that all relevant physico-chemical interactions that are present in the experimental sample are accurately modelled. The inherent difficulty associated to the calculation of chiral signals of systems in aqueous solutions requires the development of specific tools, able to account for the peculiarities of water-solute interactions, and especially its ability to form hydrogen bonds. In this perspective we discuss a multiscale approach, which we have developed and challenged to model the most used chiroptical techniques.

Solvent Effects and Aggregation Phenomena Studied by Vibrational Optical Activity and Molecular Dynamics: The Case of Pantolactone

Raman and Raman optical activity (ROA), IR, and vibrational circular dichroism (VCD) spectra of (R)- and (S)-pantolactone have been recorded in three solvents. ROA has been employed on water and DMSO solutions, VCD on DMSO and CCl₄ solutions. In the last solvent, monomer-dimer equilibrium is present. Due to the low conformational flexibility of the isolated molecule and to the possibility of aggregation, this compound has been used here to test different protocols for computation of the spectroscopic responses taking into account solvent effects. Molecular dynamics (MD) simulations have been carried out together with statistical clustering methods based on collective variables to extract the structures needed to calculate the spectra. Quantum mechanical DFT calculations based on PCM are compared with approaches based on different representations of the solvent shell (MM or QM level). Appropriate treatment of the solvent permits obtaining of good band-shapes, with the added advantage that the MD analysis allows one to take into account flexibility of dimeric structures justifying the broadness of observed bands and the absence of intense VCD couplets in the carbonyl and OH stretching regions.

Glucose in dry and moist ionic liquid: vibrational circular dichroism, IR, and possible mechanisms

Ionic liquids and their mixtures with water show remarkable features in cellulose processing. For this reason, understanding the behavior of carbohydrates in ionic liquids is important. In the present study, we investigated three d-glucose isomers (α, β and open-chain) in 1-ethyl-3-methylimidazolium acetate in the presence and absence of water, through ab initio molecular dynamics simulations. In the complex hydrogen bonding network of these mixtures, the most interesting observation is that upon water addition every hydrogen bond elongates, except the glucose-glucose hydrogen bond for the open-chain and the α-form which shortens, clearly showing the beginning of the crystallization process. The ring glucose rearranges from on-top to in-plane and the open form changes from a coiled to a more linear arrangement when adding water which explains the contradiction that the center of mass distances of the glucose molecules with other glucose molecules grow while the hydrogen bonds shorten. The appearance of coiled open forms indicates that the previously suggested isomerization between these forms is possible and might play a role in the solubility of the related carbohydrates. The calculated IR and VCD spectra reveal insight into the intermolecular interactions, with good to excellent agreements with experimental spectra. Investigating the role of the cation, distances between the acidic carbon atom of the cation and the glucose carbon atom where ring closure and opening occurs are found, which are way shorter than dispersion-like interactions between aliphatic hydrocarbons.

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.

More complex, less complicated? Explicit solvation of hydroxyl groups for the analysis of VCD spectra

With increasing size of the molecules, hydrogen bonding induced solvent effects on the IR and VCD spectra become more negligible.

How many solvent molecules are required to solvate chiral 1,2-diols with hydrogen bonding solvents? A VCD spectroscopic study

A VCD spectroscopic analysis of selected model systems for solute–solvent interactions of chiral diols with hydrogen bonding solvents DMSO and ACN.

Assessing cluster models of solvation for the description of vibrational circular dichroism spectra: synergy between static and dynamic approaches

Solvation effects are essential for defining the shape of vibrational circular dichroism (VCD) spectra.

Chiral Molecular Structures of Substituted Indans: Ring Puckering, Rotatable Substituents, and Vibrational Circular Dichroism

The chiral molecular structures of four different substituted indans, namely, (S)-1-methylindan, (R)-1-methylindan-1-d, (R)-1-aminoindan, and (S)-1-indanol, were investigated using experimental vibrational absorption and vibrational circular dichroism spectra and corresponding spectra predicted using quantum chemical (QC) calculations. All of these molecules possess two ring puckering conformations, with ring puckering leading to the pseudoequatorial substituent being approximately four times more abundant over that leading to the pseudoaxial substituent. The amino group in 1-aminoindan has three conformations arising from the rotation of NH₂ group, for each ring puckering conformation, resulting in a total of six conformations. Whereas 1-indanol in the nonhydrogen-bonding solvent CCl₄ also has six conformations similar to those of 1-aminoindan, 1-indanol in the hydrogen-bonding solvent DMSO-d ₆ adopts numerous conformations, of which 30 conformers are considered to have at least ∼1% or more population. In DMSO solution, ring puckering leading to pseudoequatorial substituent accounts for 77% population and 23% for pseudoaxial substituent. The QC spectra predicted for the geometry optimized conformers are found to be in excellent quantitative agreement with corresponding experimental spectra in all of the molecules considered. The procedures suggested in this work are hoped to provide successful pathways for future chiral molecular structural analyses.

Beyond the Rosenfeld Equation: Computation of Vibrational Circular Dichroism Spectra for Anisotropic Solutions

The difference in absorption of left and right circularly polarized light by chiral molecules can be described by the Rosenfeld equation for isotropic samples. It allows the assignment of the absolute stereochemistry by comparing experimental and computationally derived spectra. Despite the simple form of the Rosenfeld equation, its evaluation in the infrared regime remained challenging, as the contribution from the magnetic dipole operator is zero within the Born-Oppenheimer (BO) approximation. In order to resolve this issue, "beyond BO" theories had to be developed, among which Stephen's magnetic field perturbation (MFP) approach offers a computationally easily accessible form. In this work, optical activity is discussed for cylindrically symmetric solutions, which cannot be described anymore by Rosenfeld's equation due to broken spherical symmetry. Mathematical properties of natural and electric-field induced anisotropies are discussed on the basis of the gauge-independent theoretical framework of Buckingham and Dunn. The issue of achiral noise arising from external field perturbations is considered, and potential remedies are introduced. Natural anisotropic vibrational circular dichroism (VCD) equations are solved numerically by applying the MFP approach within the Hartree-Fock (HF) formalism. Properties of anisotropic VCD spectra are discussed for R-(+)-methyloxirane and (1 S,2 S)-cyclopropane-1,2-dicarbonitrile. In particular, by using a group theoretical argument, a gauge-independent lower bound for the quadrupole contribution of C₂-symmetric molecules can be identified, which allows the importance of additional quadrupole terms in anisotropic VCD spectra calculation to be assessed.

Solvation and self-aggregation of chiral alcohols: how hydrogen bonding affects their VCD spectral signatures

IR and VCD spectra of chiral alcohols in different solvents are analyzed with DFT spectra calculations. We show that for ACN or DMSO explicit solvation is needed to reproduce experimental spectra.

Regional Susceptibility in VCD Spectra to Dynamic Molecular Motions: The Case of a Benzyl α-Hydroxysilane

Experimental and theoretical studies of the vibrational circular dichroism (VCD) spectrum of 3-methyl-1-(methyldiphenlsilyl)-1-phenylbutan-1-ol, whose absolute configuration is key to elucidating the Brook rearrangement of tertiary benzylic α-hydroxylsilanes, are presented. It is found that the entire OH-bending region in this spectrum-a region that provides important marker bands-cannot be reproduced at all by standard theoretical approaches even though other regions are well described. Using a novel approach to disentangle contributions to the rotational strength of these bands, internal coordinates are identified that critically influence the appearance of this part of the spectrum. We show that the agreement between experiment and theory is greatly improved when structural dynamics along these coordinates are explicitly taken into account. The general applicability of the approach underlines its usefulness for structurally flexible chiral systems, a situation that is more the rule rather than the exception.

Chiroptical properties of 2,2'-bioxirane

The two enantiomers of 2,2'-bioxirane were synthesized, and their chiroptical properties were thoroughly investigated in various solvents by polarimetry, vibrational circular dichroism (VCD), and Raman optical activity (ROA). Density functional theory (DFT) calculations at the B3LYP/aug-cc-pVTZ level revealed the presence of three conformers (G₊ , G_- , and cis) with Gibbs populations of 51, 44, and 5% for the isolated molecule, respectively. The population ratios of the two main conformers were modified for solvents exhibiting higher dielectric constants (G_- form decreases whereas G₊ form increases). The behavior of the specific optical rotation values with the different solvents was correctly reproduced by time-dependent DFT calculations using the polarizable continuum model (PCM), except for the benzene for which explicit solvent model should be necessary. Finally, VCD and ROA spectra were perfectly reproduced by the DFT/PCM calculations for the Boltzmann-averaged G₊ and G_- conformers.

Effect of puckering motion and hydrogen bond formation on the vibrational circular dichroism spectrum of a flexible molecule: the case of (S)-1-indanol

Vibrational circular dichroism spectra of (S)-1-indanol in DMSO and CCl₄ are described by cluster-in-the-bulk static calculations and first principles molecular dynamics.

Computing Bulk Phase Raman Optical Activity Spectra from ab initio Molecular Dynamics Simulations

We present our novel methodology for computing Raman optical activity (ROA) spectra of liquid systems from ab initio molecular dynamics (AIMD) simulations. The method is built upon the recent developments to obtain magnetic dipole moments from AIMD and to integrate molecular properties by using radical Voronoi tessellation. These techniques are used to calculate optical activity tensors for large and complex periodic bulk phase systems. Only AIMD simulations are required as input, and no time-consuming perturbation theory is involved. The approach relies only on the total electron density in each time step and can readily be combined with a wide range of electronic structure methods. To the best of our knowledge, these are the first computed ROA spectra for a periodic bulk phase system. As an example, the experimental ROA spectrum of liquid (R)-propylene oxide is reproduced very well.

Vibrational optical activity as probe for intermolecular interactions

A detailed VCD spectroscopic analysis of well-selected chiral model systems can give valuable and unprecedented insights into intermolecular interactions such as solvation or reactant–substrate binding in catalysis.

Effective Fully Polarizable QM/MM Approach To Model Vibrational Circular Dichroism Spectra of Systems in Aqueous Solution

We propose a methodology, based on the combination of classical Molecular Dynamics (MD) simulations with a fully polarizable Quantum Mechanical (QM)/Molecular Mechanics (MM)/Polarizable Continuum Model (PCM) Hamiltonian, to calculate Vibrational Circular Dichroism (VCD) spectra of chiral systems in aqueous solution. Polarization effects are included in the MM force field by exploiting an approach based on Fluctuating Charges (FQ). By performing the MD, the description of the solvating environment is enriched by taking into account the dynamical aspects of the solute-solvent interactions. On the other hand, the QM/FQ/PCM calculation of the VCD spectrum ensures an accurate description of the electronic density of the solute and a proper account for the specific interactions in solution. The application of our approach to (R)-methyloxirane and (l)-alanine in aqueous solution gives calculated spectra in remarkable agreement with their experimental counterparts and a substantial improvement with respect to the same spectra calculated with the PCM.