1
|
Fusè M, Mazzeo G, Longhi G, Abbate S, Yang Q, Bloino J. Scaling-up VPT2: A feasible route to include anharmonic correction on large molecules. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2024; 311:123969. [PMID: 38330757 DOI: 10.1016/j.saa.2024.123969] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/02/2023] [Revised: 01/18/2024] [Accepted: 01/24/2024] [Indexed: 02/10/2024]
Abstract
Vibrational analysis plays a crucial role in the investigation of molecular systems. Though methodologies like second-order vibrational perturbation theory (VPT2) have paved the way to more accurate simulations, the computational cost remains a difficult barrier to overcome when the molecular size increases. Building upon recent advances in the identification of resonances, we propose an approach making anharmonic simulations possible for large-size systems, typically unreachable by standard means. This relies on the fact that, often, only portions of the whole spectra are of actual interest. Therefore, the anharmonic corrections can be included selectively on subsets of normal modes directly related to the regions of interest. Starting from the VPT2 equations, we evaluate rigorously and systematically the impact of the truncated anharmonic treatment onto simulations. The limit and feasibility of the reduced-dimensionality approach are detailed, starting on a smaller model system. The methodology is then challenged on the IR absorption and vibrational circular dichroism spectra of an organometallic complex in three different spectral ranges.
Collapse
Affiliation(s)
- Marco Fusè
- Dipartimento di Medicina Molecolare e Traslazionale, Università di Brescia, Viale Europa 11, 25123, Brescia, Italy.
| | - Giuseppe Mazzeo
- Dipartimento di Medicina Molecolare e Traslazionale, Università di Brescia, Viale Europa 11, 25123, Brescia, Italy
| | - Giovanna Longhi
- Dipartimento di Medicina Molecolare e Traslazionale, Università di Brescia, Viale Europa 11, 25123, Brescia, Italy; Istituto Nazionale di Ottica (INO), CNR, Research Unit of Brescia, c/o CSMT, VIA Branze 45, 25123, Brescia, Italy
| | - Sergio Abbate
- Dipartimento di Medicina Molecolare e Traslazionale, Università di Brescia, Viale Europa 11, 25123, Brescia, Italy; Istituto Nazionale di Ottica (INO), CNR, Research Unit of Brescia, c/o CSMT, VIA Branze 45, 25123, Brescia, Italy
| | - Qin Yang
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo náměstí 542/2, 160 00, Prague, Czech Republic
| | - Julien Bloino
- Scuola Normale Superiore, Piazza dei Cavalieri, 56125, Pisa, Italy.
| |
Collapse
|
2
|
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
|
3
|
Galimberti DR. Vibrational Circular Dichroism from DFT Molecular Dynamics: The AWV Method. J Chem Theory Comput 2022; 18:6217-6230. [PMID: 36112978 PMCID: PMC9558311 DOI: 10.1021/acs.jctc.2c00736] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2022] [Indexed: 11/29/2022]
Abstract
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.
Collapse
Affiliation(s)
- Daria Ruth Galimberti
- Institute
for Molecules and Materials, Radboud University, Heyendaalseweg 135, 6525 AJ Nijmegen, The Netherlands
| |
Collapse
|
4
|
Kirchner B, Blasius J, Alizadeh V, Gansäuer A, Hollóczki O. Chemistry Dissolved in Ionic Liquids. A Theoretical Perspective. J Phys Chem B 2022; 126:766-777. [PMID: 35034453 DOI: 10.1021/acs.jpcb.1c09092] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
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.
Collapse
Affiliation(s)
- Barbara Kirchner
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstraße 4+6, D-53115 Bonn, Germany
| | - Jan Blasius
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstraße 4+6, D-53115 Bonn, Germany
| | - Vahideh Alizadeh
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstraße 4+6, D-53115 Bonn, Germany
| | - Andreas Gansäuer
- Kekulé-Institut für Organische Chemie und Biochemie, University of Bonn, Gerhard-Domagk-Straße 1, D-53121 Bonn, Germany
| | - Oldamur Hollóczki
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstraße 4+6, D-53115 Bonn, Germany.,Department of Physical Chemistry, Faculty of Science and Technology, University of Debrecen, Egyetem tér 1, H-4010 Debrecen, Hungary
| |
Collapse
|
5
|
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
|