1
|
Sklenář A, Růžičková L, Schrenková V, Bednárová L, Pazderková M, Chatziadi A, Zmeškalová Skořepová E, Šoóš M, Kaminský J. Solid-state vibrational circular dichroism for pharmaceutical applications: Polymorphs and cocrystal of sofosbuvir. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2024; 318:124478. [PMID: 38788502 DOI: 10.1016/j.saa.2024.124478] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/25/2023] [Revised: 05/06/2024] [Accepted: 05/16/2024] [Indexed: 05/26/2024]
Abstract
X-ray diffraction is a commonly used technique in the pharmaceutical industry for the determination of the atomic and molecular structure of crystals. However, it is costly, sometimes time-consuming, and it requires a considerable degree of expertise. Vibrational circular dichroism (VCD) spectroscopy resolves these limitations, while also exhibiting substantial sensitivity to subtle modifications in the conformation and molecular packaging in the solid state. This study showcases VCD's ability to differentiate between various crystal structures of the same molecule (polymorphs, cocrystals). We examined the most effective approach for producing high-quality spectra and unveiled the intricate link between structure and spectrum via quantum-chemical computations. We rigorously assessed, using alanine as a model compound, multiple experimental conditions on the resulting VCD spectra, with the aim of proposing an optimal and efficient procedure. The proposed approach, which yields reliable, reproducible, and artifact-free results with maximal signal-to-noise ratio, was then validated using a set comprising of three amino acids (serine, alanine, tyrosine), one hydroxy acid (tartaric acid), and a monosaccharide (ribose) to mimic active pharmaceutical components. Finally, the optimized approach was applied to distinguish three polymorphs of the antiviral drug sofosbuvir and its cocrystal with piperazine. Our results indicate that solid-state VCD is a prompt, cost-effective, and easy-to-use technique to identify crystal structures, demonstrating potential for application in pharmaceuticals. We also adapted the cluster and transfer approach to calculate the spectral properties of molecules in a periodic crystal environment. Our findings demonstrate that this approach reliably produces solid-state VCD spectra of model compounds. Although for large molecules with many atoms per unit cell, such as sofosbuvir, this approach has to be simplified and provides only a qualitative match, spectral calculations, and energy analysis helped us to decipher the observed differences in the experimental spectra of sofosbuvir.
Collapse
Affiliation(s)
- Adam Sklenář
- Institute of Organic Chemistry and Biochemistry of the CAS, Flemingovo nám.2, Prague 166 10, Czech Republic; University of Chemistry and Technology, Prague, Technická 5, Prague 166 28, Czech Republic
| | - Lucie Růžičková
- Institute of Organic Chemistry and Biochemistry of the CAS, Flemingovo nám.2, Prague 166 10, Czech Republic; Imperial College London, Department of Life Sciences, South Kensington Campus, London SW7 2AZ, UK
| | - Věra Schrenková
- Institute of Organic Chemistry and Biochemistry of the CAS, Flemingovo nám.2, Prague 166 10, Czech Republic
| | - Lucie Bednárová
- Institute of Organic Chemistry and Biochemistry of the CAS, Flemingovo nám.2, Prague 166 10, Czech Republic
| | - Markéta Pazderková
- Institute of Organic Chemistry and Biochemistry of the CAS, Flemingovo nám.2, Prague 166 10, Czech Republic
| | - Argyro Chatziadi
- University of Chemistry and Technology, Prague, Technická 5, Prague 166 28, Czech Republic
| | - Eliška Zmeškalová Skořepová
- University of Chemistry and Technology, Prague, Technická 5, Prague 166 28, Czech Republic; Institute of Physics of the CAS, Na Slovance 1999/2, Prague 182 21, Czech Republic
| | - Miroslav Šoóš
- University of Chemistry and Technology, Prague, Technická 5, Prague 166 28, Czech Republic
| | - Jakub Kaminský
- Institute of Organic Chemistry and Biochemistry of the CAS, Flemingovo nám.2, Prague 166 10, Czech Republic.
| |
Collapse
|
2
|
Duston T, Tao Z, Bian X, Bhati M, Rawlinson J, Littlejohn RG, Pei Z, Shao Y, Subotnik JE. A Phase-Space Electronic Hamiltonian For Vibrational Circular Dichroism. J Chem Theory Comput 2024. [PMID: 39226223 DOI: 10.1021/acs.jctc.4c00662] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/05/2024]
Abstract
We show empirically that a phase-space non-Born-Oppenheimer electronic Hamiltonian approach to quantum chemistry (where the electronic Hamiltonian is parametrized by both nuclear position and momentum, ĤPS(R,P)) is both a practical and accurate means to recover vibrational circular dichroism spectra. We further hypothesize that such a phase-space approach may lead to very new dynamical physics beyond spectroscopic circular dichroism, with potential implications for understanding chiral induced spin selectivity (CISS), noting that classical phase-space approaches conserve the total nuclear plus electronic momentum, whereas classical Born-Oppenheimer approaches do not (they conserve only the nuclear momentum).
Collapse
Affiliation(s)
- Titouan Duston
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Zhen Tao
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Xuezhi Bian
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Mansi Bhati
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Jonathan Rawlinson
- Department of Mathematics, Nottingham Trent University, Nottingham NG1 4FQ, United Kingdom
| | - Robert G Littlejohn
- Department of Physics, University of California, Berkeley, California 94720, United States
| | - Zheng Pei
- Department of Chemistry, The University of Oklahoma, Norman, Oklahoma 73104, United States
| | - Yihan Shao
- Department of Chemistry, The University of Oklahoma, Norman, Oklahoma 73104, United States
| | - Joseph E Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
3
|
Bloino J, Jähnigen S, Merten C. After 50 Years of Vibrational Circular Dichroism Spectroscopy: Challenges and Opportunities of Increasingly Accurate and Complex Experiments and Computations. J Phys Chem Lett 2024; 15:8813-8828. [PMID: 39167088 DOI: 10.1021/acs.jpclett.4c01700] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/23/2024]
Abstract
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.
Collapse
Affiliation(s)
- Julien Bloino
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Sascha Jähnigen
- Freie Universität Berlin, Institut für Chemie und Biochemie, Arnimallee 22, 14195 Berlin, Germany
| | - Christian Merten
- Ruhr-Universität Bochum, Fakultät für Chemie und Biochemie, Universitätsstraße 150, 44801 Bochum, Germany
| |
Collapse
|
4
|
Drost DA, Merten C. Vibrational circular dichroism spectra of proline in water at different pH values. Phys Chem Chem Phys 2024; 26:17753-17759. [PMID: 38873734 DOI: 10.1039/d4cp01768d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2024]
Abstract
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.
Collapse
Affiliation(s)
- Deborah A Drost
- Ruhr University Bochum, Faculty of Chemistry and Biochemistry, Organic Chemistry II, Universitätsstraße 150, 44801 Bochum, Germany. www.mertenlab.de
| | - Christian Merten
- Ruhr University Bochum, Faculty of Chemistry and Biochemistry, Organic Chemistry II, Universitätsstraße 150, 44801 Bochum, Germany. www.mertenlab.de
| |
Collapse
|
5
|
Schreder L, Luber S. Chiral Spectroscopy of Bulk Systems with Propagated Localized Orbitals. J Chem Theory Comput 2024; 20:3894-3910. [PMID: 38661175 DOI: 10.1021/acs.jctc.4c00235] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/26/2024]
Abstract
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.
Collapse
Affiliation(s)
- Lukas Schreder
- University of Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland
| | - Sandra Luber
- University of Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland
| |
Collapse
|
6
|
Bowles J, Jähnigen S, Agostini F, Vuilleumier R, Zehnacker A, Calvo F, Clavaguéra C. Vibrational Circular Dichroism Spectroscopy with a Classical Polarizable Force Field: Alanine in the Gas and Condensed Phases. Chemphyschem 2024; 25:e202300982. [PMID: 38318765 DOI: 10.1002/cphc.202300982] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2023] [Revised: 02/05/2024] [Accepted: 02/06/2024] [Indexed: 02/07/2024]
Abstract
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.
Collapse
Affiliation(s)
- Jessica Bowles
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405, Orsay, France
| | - Sascha Jähnigen
- PASTEUR Laboratory, Département de Chimie, Ecole Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005, Paris, France
| | - Federica Agostini
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405, Orsay, France
| | - Rodolphe Vuilleumier
- PASTEUR Laboratory, Département de Chimie, Ecole Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005, Paris, France
| | - Anne Zehnacker
- Université Paris-Saclay, CNRS, Institut des Sciences Moléculaires d'Orsay UMR8214, 91405, Orsay, France
| | - Florent Calvo
- Université Grenoble Alpes, CNRS, LIPhy, 38000, Grenoble, France
| | - Carine Clavaguéra
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, 91405, Orsay, France
| |
Collapse
|
7
|
Tao Z, Qiu T, Bhati M, Bian X, Duston T, Rawlinson J, Littlejohn RG, Subotnik JE. Practical phase-space electronic Hamiltonians for ab initio dynamics. J Chem Phys 2024; 160:124101. [PMID: 38526114 DOI: 10.1063/5.0192084] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2023] [Accepted: 02/14/2024] [Indexed: 03/26/2024] Open
Abstract
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.
Collapse
Affiliation(s)
- Zhen Tao
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Tian Qiu
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Mansi Bhati
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Xuezhi Bian
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Titouan Duston
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Jonathan Rawlinson
- Department of Mathematics, University of Manchester, Manchester M13 9PL, United Kingdom
| | - Robert G Littlejohn
- Department of Physics, University of California, Berkeley, California 94720, USA
| | - Joseph E Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| |
Collapse
|
8
|
Littlejohn R, Rawlinson J, Subotnik J. Diagonalizing the Born-Oppenheimer Hamiltonian via Moyal perturbation theory, nonadiabatic corrections, and translational degrees of freedom. J Chem Phys 2024; 160:114103. [PMID: 38501907 DOI: 10.1063/5.0192465] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2023] [Accepted: 02/19/2024] [Indexed: 03/20/2024] Open
Abstract
This article describes a method for calculating higher order or nonadiabatic corrections in Born-Oppenheimer theory and its interaction with the translational degrees of freedom. The method uses the Wigner-Weyl correspondence to map nuclear operators into functions on the classical phase space and the Moyal star product to represent operator multiplication on those functions. These are explained in the body of the paper. The result is a power series in κ2, where κ = (m/M)1/4 is the usual Born-Oppenheimer parameter. The lowest order term is the usual Born-Oppenheimer approximation, while higher order terms are nonadiabatic corrections. These are needed in calculations of electronic currents, momenta, and densities. The separation of nuclear and electronic degrees of freedom takes place in the context of the exact symmetries (for an isolated molecule) of translations and rotations, and these, especially translations, are explicitly incorporated into our discussion. This article presents an independent derivation of the Moyal expansion in molecular Born-Oppenheimer theory. We show how electronic currents and momenta can be calculated within the framework of Moyal perturbation theory; we derive the transformation laws of the electronic Hamiltonian, the electronic eigenstates, and the derivative couplings under translations; we discuss in detail the rectilinear motion of the molecular center of mass in the Born-Oppenheimer representation; and we show how the elimination of the translational components of the derivative couplings leads to a unitary transformation that has the effect of exactly separating the translational degrees of freedom.
Collapse
Affiliation(s)
- Robert Littlejohn
- Department of Physics, University of California, Berkeley, California 94720, USA
| | - Jonathan Rawlinson
- School of Mathematics, University of Manchester, Manchester, United Kingdom
| | - Joseph Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| |
Collapse
|
9
|
Puente AR, Polavarapu PL. Influence of microsolvation on vibrational circular dichroism spectra in dimethyl sulfoxide solvent: A Bottom-Up approach using Quantum cluster growth. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2023; 303:123231. [PMID: 37562213 DOI: 10.1016/j.saa.2023.123231] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/22/2023] [Revised: 07/09/2023] [Accepted: 08/01/2023] [Indexed: 08/12/2023]
Abstract
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.
Collapse
Affiliation(s)
- Andrew R Puente
- Department of Chemistry, Vanderbilt University, Nashville, TN 37235, USA
| | | |
Collapse
|
10
|
Merten C. Modelling solute-solvent interactions in VCD spectra analysis with the micro-solvation approach. Phys Chem Chem Phys 2023; 25:29404-29414. [PMID: 37881890 DOI: 10.1039/d3cp03408a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2023]
Abstract
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.
Collapse
Affiliation(s)
- Christian Merten
- Ruhr Universität Bochum, Fakultät für Chemie und Biochemie, Organische Chemie II, Universitätsstraße 150, 44801 Bochum, Germany.
| |
Collapse
|
11
|
Jähnigen S. Vibrational Circular Dichroism Spectroscopy of Chiral Molecular Crystals: Insights from Theory. Angew Chem Int Ed Engl 2023; 62:e202303595. [PMID: 37071543 DOI: 10.1002/anie.202303595] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2023] [Revised: 04/18/2023] [Accepted: 04/18/2023] [Indexed: 04/19/2023]
Abstract
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.
Collapse
Affiliation(s)
- Sascha Jähnigen
- PASTEUR, Département de Chimie, Ecole Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005, Paris, France
| |
Collapse
|
12
|
Ditler E, Kumar C, Luber S. Vibrational circular dichroism spectra of natural products by means of the nuclear velocity perturbation theory. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2023; 298:122769. [PMID: 37163900 DOI: 10.1016/j.saa.2023.122769] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/22/2022] [Revised: 03/24/2023] [Accepted: 04/18/2023] [Indexed: 05/12/2023]
Abstract
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.
Collapse
Affiliation(s)
- Edward Ditler
- University of Zurich, Winterthurerstrasse 190, 8057, Zurich, Switzerland
| | - Chandan Kumar
- University of Zurich, Winterthurerstrasse 190, 8057, Zurich, Switzerland
| | - Sandra Luber
- University of Zurich, Winterthurerstrasse 190, 8057, Zurich, Switzerland.
| |
Collapse
|
13
|
Jähnigen S, Le Barbu-Debus K, Guillot R, Vuilleumier R, Zehnacker A. How Crystal Symmetry Dictates Non-Local Vibrational Circular Dichroism in the Solid State. Angew Chem Int Ed Engl 2023; 62:e202215599. [PMID: 36441537 PMCID: PMC10107176 DOI: 10.1002/anie.202215599] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2022] [Revised: 11/28/2022] [Accepted: 11/28/2022] [Indexed: 11/29/2022]
Abstract
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 P21 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.
Collapse
Affiliation(s)
- Sascha Jähnigen
- PASTEUR, Département de Chimie, Ecole Normale Supérieure, CNRS, PSL University, Sorbonne Université, 75005, Paris, France
| | - Katia Le Barbu-Debus
- Institut des Sciences Moléculaires d'Orsay (ISMO), CNRS, Université Paris-Saclay, 91405, Orsay, France
| | - Régis Guillot
- Institut de Chimie Moléculaire et des Matériaux d'Orsay (ICMMO), CNRS, Université Paris-Saclay, 91405, Orsay, France
| | - Rodolphe Vuilleumier
- PASTEUR, Département de Chimie, Ecole Normale Supérieure, CNRS, PSL University, Sorbonne Université, 75005, Paris, France
| | - Anne Zehnacker
- Institut des Sciences Moléculaires d'Orsay (ISMO), CNRS, Université Paris-Saclay, 91405, Orsay, France
| |
Collapse
|
14
|
Villaseco Arribas E, Agostini F, Maitra NT. Exact Factorization Adventures: A Promising Approach for Non-Bound States. Molecules 2022; 27:molecules27134002. [PMID: 35807246 PMCID: PMC9267945 DOI: 10.3390/molecules27134002] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2022] [Revised: 06/07/2022] [Accepted: 06/09/2022] [Indexed: 11/29/2022] Open
Abstract
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.
Collapse
Affiliation(s)
| | - Federica Agostini
- Institut de Chimie Physique UMR8000, Université Paris-Saclay, CNRS, 91405 Orsay, France;
| | - Neepa T. Maitra
- Department of Physics, Rutgers University, Newark, NJ 07102, USA;
- Correspondence:
| |
Collapse
|
15
|
Ditler E, Zimmermann T, Kumar C, Luber S. Implementation of Nuclear Velocity Perturbation and Magnetic Field Perturbation Theory in CP2K and Their Application to Vibrational Circular Dichroism. J Chem Theory Comput 2022; 18:2448-2461. [PMID: 35363490 DOI: 10.1021/acs.jctc.2c00006] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
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.
Collapse
Affiliation(s)
- Edward Ditler
- Department of Chemistry, University of Zurich, CH-8057 Zurich, Switzerland
| | - Tomáš Zimmermann
- Department of Chemistry, University of Zurich, CH-8057 Zurich, Switzerland
| | - Chandan Kumar
- Department of Chemistry, University of Zurich, CH-8057 Zurich, Switzerland
| | - Sandra Luber
- Department of Chemistry, University of Zurich, CH-8057 Zurich, Switzerland
| |
Collapse
|
16
|
Le Barbu-Debus K, Zehnacker A. Competition between inter and intramolecular hydrogen bond evidenced by vibrational circular dichroism spectroscopy: The case of (1S,2R)-(-)-cis-1-amino-2-indanol. Chirality 2021; 33:858-874. [PMID: 34570370 DOI: 10.1002/chir.23362] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2021] [Revised: 08/29/2021] [Accepted: 08/31/2021] [Indexed: 02/06/2023]
Abstract
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-d6 . 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.
Collapse
Affiliation(s)
- Katia Le Barbu-Debus
- Institut des Sciences Moléculaires d'Orsay, Université Paris-Saclay, Orsay, France
| | - Anne Zehnacker
- Institut des Sciences Moléculaires d'Orsay, Université Paris-Saclay, Orsay, France
| |
Collapse
|
17
|
Jähnigen S, Sebastiani D, Vuilleumier R. The important role of non-covalent interactions for the vibrational circular dichroism of lactic acid in aqueous solution. Phys Chem Chem Phys 2021; 23:17232-17241. [PMID: 34369531 DOI: 10.1039/d1cp03106f] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
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.
Collapse
Affiliation(s)
- Sascha Jähnigen
- PASTEUR, Département de chimie, École normale supérieure, PSL University, Sorbonne Université, CNRS, 75005, Paris, France.
| | | | | |
Collapse
|
18
|
Jähnigen S, Zehnacker A, Vuilleumier R. Computation of Solid-State Vibrational Circular Dichroism in the Periodic Gauge. J Phys Chem Lett 2021; 12:7213-7220. [PMID: 34310135 DOI: 10.1021/acs.jpclett.1c01682] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
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.
Collapse
Affiliation(s)
- Sascha Jähnigen
- Institut des Sciences Moléculaires d'Orsay (ISMO), CNRS, Université Paris-Saclay, 91405 Orsay, France
- PASTEUR Laboratory, Département de Chimie, Ecole Normale Supérieure, PSL University, Sorbonne Université, CNRS,, 75005 Paris, France
| | - Anne Zehnacker
- Institut des Sciences Moléculaires d'Orsay (ISMO), CNRS, Université Paris-Saclay, 91405 Orsay, France
| | - Rodolphe Vuilleumier
- PASTEUR Laboratory, Département de Chimie, Ecole Normale Supérieure, PSL University, Sorbonne Université, CNRS,, 75005 Paris, France
| |
Collapse
|
19
|
Kirchner B, Blasius J, Esser L, Reckien W. Predicting Vibrational Spectroscopy for Flexible Molecules and Molecules with Non‐Idle Environments. ADVANCED THEORY AND SIMULATIONS 2020. [DOI: 10.1002/adts.202000223] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Barbara Kirchner
- Mulliken Center for Theoretical Chemistry Rheinische Friedrich‐Wilhelms‐Universität Bonn Beringstr. 4+6 D‐53115 Bonn Germany
| | - Jan Blasius
- Mulliken Center for Theoretical Chemistry Rheinische Friedrich‐Wilhelms‐Universität Bonn Beringstr. 4+6 D‐53115 Bonn Germany
| | - Lars Esser
- Mulliken Center for Theoretical Chemistry Rheinische Friedrich‐Wilhelms‐Universität Bonn Beringstr. 4+6 D‐53115 Bonn Germany
| | - Werner Reckien
- Mulliken Center for Theoretical Chemistry Rheinische Friedrich‐Wilhelms‐Universität Bonn Beringstr. 4+6 D‐53115 Bonn Germany
| |
Collapse
|
20
|
Giovannini T, Egidi F, Cappelli C. Theory and algorithms for chiroptical properties and spectroscopies of aqueous systems. Phys Chem Chem Phys 2020; 22:22864-22879. [PMID: 33043930 DOI: 10.1039/d0cp04027d] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
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.
Collapse
Affiliation(s)
- Tommaso Giovannini
- Department of Chemistry, Norwegian University of Science and Technology (NTNU), 7491 Trondheim, Norway
| | | | | |
Collapse
|
21
|
Blasius J, Kirchner B. Cluster-Weighting in Bulk Phase Vibrational Circular Dichroism. J Phys Chem B 2020; 124:7272-7283. [DOI: 10.1021/acs.jpcb.0c06313] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Jan Blasius
- Mulliken Center for Theoretical Chemistry, Rheinische Friedrich-Wilhelms-Universität Bonn, Beringstr. 4+6, D-53115 Bonn, Germany
| | - Barbara Kirchner
- Mulliken Center for Theoretical Chemistry, Rheinische Friedrich-Wilhelms-Universität Bonn, Beringstr. 4+6, D-53115 Bonn, Germany
| |
Collapse
|
22
|
Brehm M, Thomas M, Gehrke S, Kirchner B. TRAVIS—A free analyzer for trajectories from molecular simulation. J Chem Phys 2020; 152:164105. [DOI: 10.1063/5.0005078] [Citation(s) in RCA: 164] [Impact Index Per Article: 41.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023] Open
Affiliation(s)
- M. Brehm
- Institut für Chemie, Martin-Luther-Universität Halle–Wittenberg, von-Danckelmann-Platz 4, D-06120 Halle (Saale), Germany
| | - M. Thomas
- Institut für Chemie, Martin-Luther-Universität Halle–Wittenberg, von-Danckelmann-Platz 4, D-06120 Halle (Saale), Germany
| | - S. Gehrke
- Mulliken Center for Theoretical Chemistry, Rheinische Friedrich-Wilhelms-Universität Bonn, Beringstr. 4+6, D-53115 Bonn, Germany
| | - B. Kirchner
- Mulliken Center for Theoretical Chemistry, Rheinische Friedrich-Wilhelms-Universität Bonn, Beringstr. 4+6, D-53115 Bonn, Germany
| |
Collapse
|
23
|
Nafie LA. Vibrational optical activity: From discovery and development to future challenges. Chirality 2020; 32:667-692. [DOI: 10.1002/chir.23191] [Citation(s) in RCA: 33] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2019] [Revised: 01/21/2020] [Accepted: 01/23/2020] [Indexed: 12/13/2022]
|
24
|
Polavarapu PL, Santoro E. Vibrational optical activity for structural characterization of natural products. Nat Prod Rep 2020; 37:1661-1699. [DOI: 10.1039/d0np00025f] [Citation(s) in RCA: 44] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
This review presents the recent progress towards elucidating the structures of chiral natural products and applications using vibrational optical activity (VOA) spectroscopy.
Collapse
|
25
|
Le Barbu-Debus K, Bowles J, Jähnigen S, Clavaguéra C, Calvo F, Vuilleumier R, Zehnacker A. Assessing cluster models of solvation for the description of vibrational circular dichroism spectra: synergy between static and dynamic approaches. Phys Chem Chem Phys 2020; 22:26047-26068. [DOI: 10.1039/d0cp03869e] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Solvation effects are essential for defining the shape of vibrational circular dichroism (VCD) spectra.
Collapse
Affiliation(s)
- Katia Le Barbu-Debus
- Institut des Sciences Moléculaires d’Orsay (ISMO)
- CNRS
- Université Paris-Saclay
- F-91405 Orsay
- France
| | - Jessica Bowles
- Université Paris-Saclay
- CNRS
- Institut de Chimie Physique
- UMR8000
- 91405 Orsay
| | - Sascha Jähnigen
- PASTEUR
- Département de Chimie
- Ecole Normale Supérieure
- PSL University
- Sorbonne Université
| | - Carine Clavaguéra
- Université Paris-Saclay
- CNRS
- Institut de Chimie Physique
- UMR8000
- 91405 Orsay
| | - Florent Calvo
- Université Grenoble Alpes
- CNRS
- LiPhy
- F-38000 Grenoble
- France
| | - Rodolphe Vuilleumier
- PASTEUR
- Département de Chimie
- Ecole Normale Supérieure
- PSL University
- Sorbonne Université
| | - Anne Zehnacker
- Institut des Sciences Moléculaires d’Orsay (ISMO)
- CNRS
- Université Paris-Saclay
- F-91405 Orsay
- France
| |
Collapse
|
26
|
Johnson J, Polavarapu PL. Chiral Molecular Structures of Substituted Indans: Ring Puckering, Rotatable Substituents, and Vibrational Circular Dichroism. ACS OMEGA 2019; 4:4963-4976. [PMID: 31459680 PMCID: PMC6648482 DOI: 10.1021/acsomega.8b03628] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/26/2018] [Accepted: 02/21/2019] [Indexed: 05/06/2023]
Abstract
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 NH2 group, for each ring puckering conformation, resulting in a total of six conformations. Whereas 1-indanol in the nonhydrogen-bonding solvent CCl4 also has six conformations similar to those of 1-aminoindan, 1-indanol in the hydrogen-bonding solvent DMSO-d 6 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.
Collapse
|
27
|
Sidler D, Bleiziffer P, Riniker S. Beyond the Rosenfeld Equation: Computation of Vibrational Circular Dichroism Spectra for Anisotropic Solutions. J Chem Theory Comput 2019; 15:2492-2503. [PMID: 30802403 DOI: 10.1021/acs.jctc.8b01156] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
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 C2-symmetric molecules can be identified, which allows the importance of additional quadrupole terms in anisotropic VCD spectra calculation to be assessed.
Collapse
Affiliation(s)
- Dominik Sidler
- Laboratory of Physical Chemistry , ETH Zürich , Vladimir-Prelog-Weg 2 , 8093 Zürich , Switzerland
| | - Patrick Bleiziffer
- Laboratory of Physical Chemistry , ETH Zürich , Vladimir-Prelog-Weg 2 , 8093 Zürich , Switzerland
| | - Sereina Riniker
- Laboratory of Physical Chemistry , ETH Zürich , Vladimir-Prelog-Weg 2 , 8093 Zürich , Switzerland
| |
Collapse
|
28
|
Jähnigen S, Scherrer A, Vuilleumier R, Sebastiani D. VCD‐Verstärkung durch chirale Packungseffekte in molekularen Kristallen. Angew Chem Int Ed Engl 2018. [DOI: 10.1002/ange.201805671] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Sascha Jähnigen
- Martin-Luther-Universität Halle-Wittenberg Institut für Chemie von-Danckelmann-Platz 4 Halle 06120 Deutschland
| | - Arne Scherrer
- Martin-Luther-Universität Halle-Wittenberg Institut für Chemie von-Danckelmann-Platz 4 Halle 06120 Deutschland
| | - Rodolphe Vuilleumier
- PASTEUR, Département de chimie, École normale supérieure, PSL University, Sorbonne Université, CNRS 75005 Paris Frankreich
| | - Daniel Sebastiani
- Martin-Luther-Universität Halle-Wittenberg Institut für Chemie von-Danckelmann-Platz 4 Halle 06120 Deutschland
| |
Collapse
|
29
|
Jähnigen S, Scherrer A, Vuilleumier R, Sebastiani D. Chiral Crystal Packing Induces Enhancement of Vibrational Circular Dichroism. Angew Chem Int Ed Engl 2018; 57:13344-13348. [DOI: 10.1002/anie.201805671] [Citation(s) in RCA: 32] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2018] [Revised: 06/25/2018] [Indexed: 01/15/2023]
Affiliation(s)
- Sascha Jähnigen
- Martin-Luther-Universität Halle-WittenbergInstitut für Chemie von-Danckelmann-Platz 4 Halle 06120 Germany
| | - Arne Scherrer
- Martin-Luther-Universität Halle-WittenbergInstitut für Chemie von-Danckelmann-Platz 4 Halle 06120 Germany
| | - Rodolphe Vuilleumier
- PASTEUR, Département de chimie, École normale supérieure, PSL University, Sorbonne Université, CNRS 75005 Paris France
| | - Daniel Sebastiani
- Martin-Luther-Universität Halle-WittenbergInstitut für Chemie von-Danckelmann-Platz 4 Halle 06120 Germany
| |
Collapse
|
30
|
Le Barbu-Debus K, Scherrer A, Bouchet A, Sebastiani D, Vuilleumier R, Zehnacker A. Effect of puckering motion and hydrogen bond formation on the vibrational circular dichroism spectrum of a flexible molecule: the case of (S)-1-indanol. Phys Chem Chem Phys 2018; 20:14635-14646. [DOI: 10.1039/c8cp01695j] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Vibrational circular dichroism spectra of (S)-1-indanol in DMSO and CCl4 are described by cluster-in-the-bulk static calculations and first principles molecular dynamics.
Collapse
Affiliation(s)
- Katia Le Barbu-Debus
- Institut des Sciences Moléculaires d’Orsay (ISMO)
- CNRS
- Univ. Paris-Sud
- Université Paris-Saclay
- F-91405 Orsay
| | - Arne Scherrer
- PASTEUR
- Département de chimie
- École normale supérieure
- PSL University
- Sorbonne Université
| | - Aude Bouchet
- Institut des Sciences Moléculaires d’Orsay (ISMO)
- CNRS
- Univ. Paris-Sud
- Université Paris-Saclay
- F-91405 Orsay
| | - Daniel Sebastiani
- Martin-Luther-Universität Halle-Wittenberg
- Institut für Chemie
- 06120 Halle
- Germany
| | - Rodolphe Vuilleumier
- PASTEUR
- Département de chimie
- École normale supérieure
- PSL University
- Sorbonne Université
| | - Anne Zehnacker
- Institut des Sciences Moléculaires d’Orsay (ISMO)
- CNRS
- Univ. Paris-Sud
- Université Paris-Saclay
- F-91405 Orsay
| |
Collapse
|
31
|
Giovannini T, Del Frate G, Lafiosca P, Cappelli C. Effective computational route towards vibrational optical activity spectra of chiral molecules in aqueous solution. Phys Chem Chem Phys 2018; 20:9181-9197. [DOI: 10.1039/c8cp00487k] [Citation(s) in RCA: 37] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Abstract
A polarizable QM/MM approach to accurately compute the Vibrational Optical Activity (VOA) spectra of chiral systems is proposed and applied to aqueous solutions of (l)-methyl lactate and (S)-glycidol.
Collapse
|
32
|
Brehm M, Thomas M. Computing Bulk Phase Raman Optical Activity Spectra from ab initio Molecular Dynamics Simulations. J Phys Chem Lett 2017; 8:3409-3414. [PMID: 28685571 DOI: 10.1021/acs.jpclett.7b01616] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
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.
Collapse
Affiliation(s)
- Martin Brehm
- Institut für Chemie - Theoretische Chemie, Martin-Luther-Universität Halle-Wittenberg , Von-Danckelmann-Platz 4, 06120 Halle (Saale), Germany
| | - Martin Thomas
- Institut für Chemie - Theoretische Chemie, Martin-Luther-Universität Halle-Wittenberg , Von-Danckelmann-Platz 4, 06120 Halle (Saale), Germany
| |
Collapse
|
33
|
Srebro-Hooper M, Autschbach J. Calculating Natural Optical Activity of Molecules from First Principles. Annu Rev Phys Chem 2017; 68:399-420. [DOI: 10.1146/annurev-physchem-052516-044827] [Citation(s) in RCA: 94] [Impact Index Per Article: 13.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
| | - Jochen Autschbach
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260
| |
Collapse
|
34
|
Reiter K, Kühn M, Weigend F. Vibrational circular dichroism spectra for large molecules and molecules with heavy elements. J Chem Phys 2017; 146:054102. [DOI: 10.1063/1.4974897] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023] Open
|
35
|
Scherrer A, Vuilleumier R, Sebastiani D. Vibrational circular dichroism from ab initio molecular dynamics and nuclear velocity perturbation theory in the liquid phase. J Chem Phys 2016; 145:084101. [DOI: 10.1063/1.4960653] [Citation(s) in RCA: 38] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
36
|
Eich FG, Agostini F. The adiabatic limit of the exact factorization of the electron-nuclear wave function. J Chem Phys 2016; 145:054110. [DOI: 10.1063/1.4959962] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- F. G. Eich
- Max Planck Institute for the Structure and Dynamics of Matter, Luruper Chaussee 149, D-22761 Hamburg, Germany
| | - Federica Agostini
- Max Planck Institute of Microstructure Physics, Weinberg 2, D-06120 Halle, Germany
| |
Collapse
|
37
|
Schild A, Agostini F, Gross EKU. Electronic Flux Density beyond the Born–Oppenheimer Approximation. J Phys Chem A 2016; 120:3316-25. [DOI: 10.1021/acs.jpca.5b12657] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Axel Schild
- Max-Planck Institut für Mikrostrukturphysik, Weinberg 2, D-06120 Halle, Germany
| | - Federica Agostini
- Max-Planck Institut für Mikrostrukturphysik, Weinberg 2, D-06120 Halle, Germany
| | - E. K. U. Gross
- Max-Planck Institut für Mikrostrukturphysik, Weinberg 2, D-06120 Halle, Germany
| |
Collapse
|
38
|
Thomas M, Kirchner B. Classical Magnetic Dipole Moments for the Simulation of Vibrational Circular Dichroism by ab Initio Molecular Dynamics. J Phys Chem Lett 2016; 7:509-513. [PMID: 26771403 DOI: 10.1021/acs.jpclett.5b02752] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
We present a new approach for calculating vibrational circular dichroism spectra by ab initio molecular dynamics. In the context of molecular dynamics, these spectra are given by the Fourier transform of the cross-correlation function of magnetic dipole moment and electric dipole moment. We obtain the magnetic dipole moment from the electric current density according to the classical definition. The electric current density is computed by solving a partial differential equation derived from the continuity equation and the condition that eddy currents should be absent. In combination with a radical Voronoi tessellation, this yields an individual magnetic dipole moment for each molecule in a bulk phase simulation. Using the chiral alcohol 2-butanol as an example, we show that experimental spectra are reproduced very well. Our approach requires knowing only the electron density in each simulation step, and it is not restricted to any particular electronic structure method.
Collapse
Affiliation(s)
- Martin Thomas
- Mulliken Center for Theoretical Chemistry, Rheinische Friedrich-Wilhelms-Universität Bonn , Beringstraße 4, 53115 Bonn, Germany
| | - Barbara Kirchner
- Mulliken Center for Theoretical Chemistry, Rheinische Friedrich-Wilhelms-Universität Bonn , Beringstraße 4, 53115 Bonn, Germany
| |
Collapse
|
39
|
Scherrer A, Agostini F, Sebastiani D, Gross EKU, Vuilleumier R. Nuclear velocity perturbation theory for vibrational circular dichroism: An approach based on the exact factorization of the electron-nuclear wave function. J Chem Phys 2015; 143:074106. [DOI: 10.1063/1.4928578] [Citation(s) in RCA: 55] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Arne Scherrer
- Martin-Luther-University Halle-Wittenberg, von-Danckelmann-Platz 4, D-06120 Halle, Germany
- UMR 8640 ENS-CNRS-UPMC, Département de Chimie, 24 rue Lhomond, École Normale Supérieure, 75005 Paris, France
- UPMC Université Paris 06, 4, Place Jussieu, 75005 Paris, France
| | - Federica Agostini
- Max-Planck-Institute of Microstructure Physics, Weinberg 2, D-06120 Halle, Germany
| | - Daniel Sebastiani
- Martin-Luther-University Halle-Wittenberg, von-Danckelmann-Platz 4, D-06120 Halle, Germany
| | - E. K. U. Gross
- Max-Planck-Institute of Microstructure Physics, Weinberg 2, D-06120 Halle, Germany
| | - Rodolphe Vuilleumier
- UMR 8640 ENS-CNRS-UPMC, Département de Chimie, 24 rue Lhomond, École Normale Supérieure, 75005 Paris, France
- UPMC Université Paris 06, 4, Place Jussieu, 75005 Paris, France
| |
Collapse
|