Raman Optical Activity Spectra from Density Functional Perturbation Theory and Density-Functional-Theory-Based Molecular Dynamics

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.

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.

Vibrational circular dichroism spectra of natural products by means of the nuclear velocity perturbation theory

We present the application of the recently implemented nuclear velocity perturbation theory, using the combined Gaussian and plane waves approach in CP2K, to the vibrational circular dichroism (VCD) spectra of a set of natural products. Even though the calculations were carried out for isolated molecules in the gas-phase limit, neglecting inter-molecular interactions and anharmonic effects, the match between simulated and experimental spectra is reasonable. We also study the influence of different density functionals on the conformational search and the resulting VCD spectra via group coupling matrices (GCMs). The GCM analysis reveals that the VCD signal can in some cases arise from moieties which are close to each other and in other cases from moieties far from each other. Differences in spectra obtained using different exchange-correlation density functionals can be attributed to interaction terms between different moieties in the molecules changing their sign.

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.

A Concise Review on Recent Developments of Machine Learning for the Prediction of Vibrational Spectra

Machine learning has become more and more popular in computational chemistry, as well as in the important field of spectroscopy. In this concise review, we walk the reader through a short summary of machine learning algorithms and a comprehensive discussion on the connection between machine learning methods and vibrational spectroscopy, particularly for the case of infrared and Raman spectroscopy. We also briefly discuss state-of-the-art molecular representations which serve as meaningful inputs for machine learning to predict vibrational spectra. In addition, this review provides an overview of the transferability and best practices of machine learning in the prediction of vibrational spectra as well as possible future research directions.

Analytic calculation and analysis of atomic polar tensors for molecules and materials using the Gaussian and plane waves approach

The evaluation of atomic polar tensors and Born Effective Charge (BEC) tensors from Density Functional Perturbation Theory (DFPT) has been implemented in the CP2K code package. This implementation is based on a combination of the Gaussian and plane wave approach for the description of basis functions and arising potentials. The presence of non-local pseudo-potentials has been considered, as well as contributions arising from the basis functions being centered on the atoms. Simulations of both periodic and non-periodic systems have been implemented and carried out. Dipole strengths and infrared absorption spectra have been calculated for two isomers of the tripeptide Ser-Pro-Ala using DFPT and are compared to the results of standard vibrational analyses using finite differences. The spectra are then decomposed into five subsets by employing localized molecular orbitals/maximally localized Wannier functions, and the results are discussed. Moreover, group coupling matrices are employed for visualization of results. Furthermore, the BECs and partial charges of the surface atoms of a periodic (101) anatase (TiO₂) slab have been investigated in a periodic framework.

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.

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.

Improving Vibrational Mode Interpretation Using Bayesian Regression

To streamline the interpretation of vibrational spectra, this work introduces the use of Bayesian linear regression with automatic relevance determination as a viable approach to decompose the atomic motions along any vibrational mode as a weighted combination of displacements along chemically meaningful internal coordinates. This novel approach denominated vibrational mode automatic relevance determination (VMARD) is presented and compared with the well-established potential energy decomposition (PED) scheme. Good agreement is generally attained between the two methods. VMARD returns a decomposition of the atomic displacement using only a small number of internal coordinates, thus aiding the interpretation of the vibrational spectra. Moreover, the results show that the VMARD descriptions are resilient toward the addition of additional internal coordinates, achieving a concise description of the vibrational modes despite the use of redundant internal coordinates. Potential applications of VMARD involving the gathering of physical insights on the atomic motions along the reaction coordinate at transition state structures, as well as the improvement of theoretically predicted vibrational frequencies, are also presented under a proof-of-concept perspective.

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.

Localized molecular orbitals for calculation and analysis of vibrational Raman optical activity

First calculations of vibrational Raman optical activity based on localized molecular orbitals are presented, which pave the way for novel insight into spectroscopic signatures of chiral systems.

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.