Classical Magnetic Dipole Moments for the Simulation of Vibrational Circular Dichroism by ab Initio Molecular Dynamics

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

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

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.

Exploring the Influence of the Phosphorus-Heteroatom Substitution in Nicotine on Its Electronic and Vibrational Spectroscopic Properties

The influence of phosphorus substitution of nitrogen in heterocyclic compounds on the vibrational spectroscopy as well as frontier molecular orbitals are analyzed. Nicotine with two nitrogen atoms in its structure is taken as the sample system to be studied computationally. By replacing the nitrogen atom in one or both rings of this molecule with phosphorus, three nicotine derivatives are created. The vibrational circular dichroism and infrared spectra of these four molecules in their monomer state, as well as the assemblies up to trimers are determined. The aforementioned spectra are calculated using static quantum chemical calculations employing a cluster-weighted approach. The calculated gas phase spectra of nicotine are compared to their respective experimental spectra. It is observed that the nicotine derivatives with phosphorus in the methylpyrrolidine ring have considerably different gas phase and bulk phase vibrational circular dichroism spectra when compared to nicotine. The phosphorus substitution reduces the energy gap between the highest occupied molecular orbital and the lowest unoccupied molecular orbital as well as altering the polarizability and reactivity of the investigated molecules.

Efficient Prediction of Mole Fraction Related Vibrational Frequency Shifts

While so far it has been possible to calculate vibrational spectra of mixtures at a particular composition, we present here a novel cluster approach for a fast and robust calculation of mole fraction dependent infrared and vibrational circular dichroism spectra at the example of acetonitrile/(R)-butan-2-ol mixtures. By assigning weights to a limited number of quantum chemically calculated clusters, vibrational spectra can be obtained at any desired composition by a weighted average of the single cluster spectra. In this way, peak positions carrying information about intermolecular interactions can be predicted. We show that mole fraction dependent peak shifts can be accurately modeled and, that experimentally recorded infrared spectra can be reproduced with high accuracy over the entire mixing range. Because only a very limited number of clusters is required, the presented approach is a valuable and computationally efficient tool to access mole fraction dependent spectra of mixtures on a routine basis.

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.

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.

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.

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.

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.

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.

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.

Optimized Atomic Partial Charges and Radii Defined by Radical Voronoi Tessellation of Bulk Phase Simulations

We present a novel method for the computation of well-defined optimized atomic partial charges and radii from the total electron density. Our method is based on a two-step radical Voronoi tessellation of the (possibly periodic) system and subsequent integration of the total electron density within each Voronoi cell. First, the total electron density is partitioned into the contributions of each molecule, and subsequently the electron density within each molecule is assigned to the individual atoms using a second set of atomic radii for the radical Voronoi tessellation. The radii are optimized on-the-fly to minimize the fluctuation (variance) of molecular and atomic charges. Therefore, our method is completely free of empirical parameters. As a by-product, two sets of optimized atomic radii are produced in each run, which take into account many specific properties of the system investigated. The application of an on-the-fly interpolation scheme reduces discretization noise in the Voronoi integration. The approach is particularly well suited for the calculation of partial charges in periodic bulk phase systems. We apply the method to five exemplary liquid phase simulations and show how the optimized charges can help to understand the interactions in the systems. Well-known effects such as reduced ion charges below unity in ionic liquid systems are correctly predicted without any tuning, empiricism, or rescaling. We show that the basis set dependence of our method is very small. Only the total electron density is evaluated, and thus, the approach can be combined with any electronic structure method that provides volumetric total electron densities-it is not limited to Hartree-Fock or density functional theory (DFT). We have implemented the method into our open-source software tool TRAVIS.

Quantitative Comparison of Experimental and Computed IR-Spectra Extracted from Ab Initio Molecular Dynamics

Experimentally measured infrared spectra are often compared to their computed equivalents. However, the accordance is typically characterized by visual inspection, which is prone to subjective judgment. The primary challenge for a similarity-based analysis is that the artifacts introduced by each approach are very different and, therefore, may require preprocessing steps to determine and correct impeding irregularities. To allow for automated objective assessment, we propose a practical and comprehensive workflow involving scaling factors, a novel baseline correction scheme, and peak smoothing. The resulting spectra can then easily be compared quantitatively using similarity measures, for which we found the Pearson correlation coefficient to be the most suitable. The proposed procedure is then applied to compare the agreement of the experimental infrared spectra from the NIST Chemistry Web book with the calculated spectra using standard harmonic frequency analysis and spectra extracted from ab initio molecular dynamics simulations at different levels of theory. We conclude that the direct, quantitative comparison of calculated and measured IR spectra might become a novel, sophisticated approach to benchmark quantum-chemical methods. In the present benchmark, simulated spectra based on ab initio molecular dynamics show in general better agreement with the experiment than static calculations.

A Systematic Study of DFT Performance for Geometry Optimizations of Ionic Liquid Clusters

Clusters of two ion pairs of imidazolium-based ionic liquids were optimized with 43 different levels of theory, including DFT functionals and MP2-based methods combined with varying Dunning's basis sets, and added dispersion corrections. Better preforming DFT functionals were then applied to clusters consisting of four ion pairs. Excellent performance of some DFT functionals for the two ion pair clusters did not always match that of the four ion-paired clusters despite interionic distances remaining constant between the optimized two and four ion-paired clusters of the same ionic liquid. Combinations of DFT functional and basis set such as ωB97X-D/cc-pVDZ, M06-2X/aug-cc-pVDZ, B3LYP-D3/cc-pVTZ, and TPSS-D3/cc-pVTZ gave excellent results for geometry optimization of two ion-paired clusters of imidazolium ionic liquids but gave larger deviations when applied to the four ion-paired clusters of varying ionic liquids. Empirical dispersion corrections were seen to be crucial in correctly capturing correlation effects in the studied ionic liquid clusters, becoming more important in larger clusters. Dunning's double-ζ basis set, cc-pVDZ, is associated with the smallest root mean squared deviations for geometries; however, it also produces the largest deviations in total electronic energies. ωB97X-D and M06-2X produced the best performance with the augmented version of this basis set. The triple-ζ basis set, cc-pVTZ, leads to the best performance of most of the DFT functionals (especially the dispersion-corrected ones) used, whereas its augmented version, aug-cc-pVTZ, was not seen to improve results. The combinations of functional and basis set that gave the best geometry and energetics in both two and four ion-paired clusters were PBE-D3/cc-pVTZ, ωB97X-D/aug-cc-pVDZ, and BLYP-D3/cc-pVTZ. All three combinations are recommended for geometry optimizations of larger clusters of ionic liquids. PBE-D3/cc-pVTZ performed the best with an average deviation of 2.3 kJ mol^-1 and a standard deviation of 3.4 kJ mol^-1 for total electronic energy when applied to four ion-paired clusters. Geometries optimized with FMO2-SRS-MP2/cc-pVTZ produced total energy within 2.0 kJ mol^-1 off the benchmark in two ion-paired clusters, with the cc-pVDZ basis set performing unsurprisingly poorly with the same method. The error increased to 4.8 kJ mol^-1 on average in four ion-paired clusters, with the smallest RMSD deviations in geometries when compared to the benchmark ones. This study is the first report that investigated the performance of DFT functionals for two and four ion-paired clusters of a wide range of ionic liquids consisting of commonly used cations such as pyrrolidinium, imidazolium, pyridinium, and ammonium. It also identified the importance of assessing the performance of quantum chemical methods for ionic liquids on a variety of cation-anion combinations.

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.

Anharmonicity of Vibrational Modes in Hydrogen Chloride-Water Mixtures

A thorough analysis of molecular vibrations in the binary system hydrogen chloride/water is presented considering a set of small mixed and pure clusters. In addition to the conventional normal-mode analysis based on the diagonalization of the Hessian, anharmonic frequencies were obtained from the perturbative VPT2 and PT2-VSCF method using hybrid density functional theory. For all normal modes, potential energy curves were modeled by displacing the atoms from the minimum geometry along the normal mode vectors. Three model potentials, a harmonic potential, a Morse potential, and a fourth order polynomial, were applied to fit these curves. From these data, it was possible not only to characterize distinct vibrations as mainly harmonic, anharmonic, or involving higher order terms but also to extract force constants, k, and anharmonicity constants, x_e. By investigating all different types of intramolecular vibrations including covalent stretching or bending vibrations and intermolecular vibrations such as librations, we could demonstrate that while vibrational frequencies can be obtained applying scaling factors to harmonic results, useful anharmonicity constants cannot be predicted in such a way and the usage of more elaborate vibrational methods is necessary. For each particular type of molecular vibration, we could however determine a relationship between the wavenumber or wavenumber shift and the anharmonicity constant, which allows us to estimate mode dependent anharmonicity constants for larger clusters in the future.

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.

Simulating structure and dynamics in small droplets of 1-ethyl-3-methylimidazolium acetate

To investigate the structure and dynamics of small ionic liquid droplets in gas phase, we performed a DFT-based ab initio molecular dynamics study of several 1-ethyl-3-methylimidazolium acetate clusters in vacuum as well as a bulk phase simulation. We introduce an unbiased criterion for average droplet diameter and density. By extrapolation of the droplet densities, we predict the experimental bulk phase density with a deviation of only a few percent. The hydrogen bond geometry between cations and anions is very similar in droplets and bulk, but the hydrogen bond dynamics is significantly slower in the droplets, becoming slower with increasing system size, with hydrogen bond lifetimes up to 2000 ps. From a normal mode analysis of the trajectories, we identify the modes of the ring proton C-H stretching, which are strongly affected by hydrogen bonding. From analyzing these, we find that the hydrogen bond becomes weaker with increasing system size. The cations possess an increased concentration inside the clusters, whereas the anions show an excess concentration on the outside. Almost all anions point towards the droplet center with their carboxylic groups. Ring stacking is found to be a very important structural motif in the droplets (as in the bulk), but side chain interactions are only of minor importance. By using Voronoi tessellation, we define the exposed droplet surface and find that it consists mainly of hydrogen atoms from the cation's and anion's methyl and ethyl groups. Polar atoms are rarely found on the surface, such that the droplets appear completely hydrophobic on the outside.

An Efficient Lossless Compression Algorithm for Trajectories of Atom Positions and Volumetric Data

We present our newly developed and highly efficient lossless compression algorithm for trajectories of atom positions and volumetric data. The algorithm is designed as a two-step approach. In the first step, efficient polynomial extrapolation schemes reduce the information entropy of the data by exploiting both spatial and temporal continuity. The second step processes the data by a series of transformations (Burrows-Wheeler, move-to-front, run length encoding) and finally compresses the stream with multitable canonical Huffman coding. Our approach reaches a compression ratio of around 15:1 for typical position trajectories in the XYZ format. For volumetric data trajectories in Gaussian Cube format (such as electron density), even a compression ratio of around 35:1 is yielded, which is by far the smallest size of all formats compared here. At the same time, compression and decompression are still reasonably fast for everyday use. The precision of the data can be selected by the user. For storage of the compressed data, we introduce the BQB file format, which is very robust, flexible, and efficient. In contrast to most archiving formats, it allows fast random access to individual trajectory frames. Our method is implemented in C++ and provided as free software under the GNU LGPL license. It has been included in the TRAVIS program package but is also available as stand-alone tool and as a library ("libbqb") for use in other projects.

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.

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.