1
|
Maya J, Malloum A, Fifen JJ, Dhaouadi Z, Fouda HPE, Conradie J. Quantum cluster equilibrium theory applied to liquid ammonia. J Comput Chem 2024; 45:1279-1288. [PMID: 38353541 DOI: 10.1002/jcc.27327] [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: 11/06/2023] [Revised: 01/21/2024] [Accepted: 01/29/2024] [Indexed: 04/19/2024]
Abstract
Through this paper, the authors propose using the quantum cluster equilibrium (QCE) theory to reinvestigate ammonia clusters in the liquid phase. The ammonia clusters from size monomer to hexadecamer were considered to simulate the liquid ammonia in this approach. The clusterset used to model the liquid ammonia is an ensemble of different structures of ammonia clusters. After studious research of the representative configurations of ammonia clusters through the cluster research program ABCluster, the configurations have been optimized at the MN15/6-31++G(d,p) level of theory. These optimizations lead to geometries and frequencies as inputs for the Peacemaker code. The QCE study of this molecular system permits us to get the liquid phase populations in a temperature range of 190-260 K, covering the temperatures from the melting point to the boiling point. The results show that the population of liquid ammonia comprises mainly the ammonia hexadecamer followed by pentadecamer, tetradecamer, and tridecamer. We noted that the small-sized ammonia clusters do not contribute to the population of liquid ammonia. In addition, the thermodynamic properties, such as heat of vaporization, heat capacity, entropy, enthalpy, and free energies, obtained by the QCE theory have been compared to the experiment given some relatively good agreements in the gas phase and show considerable discrepancies in liquid phase except the density. Finally, based on the predicted population, we calculated the infrared spectrum of liquid ammonia at 215 K temperature. It comes out that the calculated infrared spectrum qualitatively agrees with the experiment.
Collapse
Affiliation(s)
- Josué Maya
- Department of Physics, Faculty of Science, University of Ngaoundere, Ngaoundere, Cameroon
- National Radiation Protection Agency, Yaounde, Cameroon
| | - Alhadji Malloum
- Department of Physics, Faculty of Science, University of Maroua, Maroua, Cameroon
- Department of Chemistry, University of the Free State, Bloemfontein, South Africa
| | - Jean Jules Fifen
- Department of Physics, Faculty of Science, University of Ngaoundere, Ngaoundere, Cameroon
| | - Zoubeida Dhaouadi
- Laboratoire de Spectroscopie Atomique Moléculaire et Application, Université de Tunis El Manar, Tunis, Tunisie
| | | | - Jeanet Conradie
- Department of Chemistry, University of the Free State, Bloemfontein, South Africa
| |
Collapse
|
2
|
Frömbgen T, Drysch K, Zaby P, Dölz J, Ingenmey J, Kirchner B. Quantum Cluster Equilibrium Theory for Multicomponent Liquids. J Chem Theory Comput 2024; 20:1838-1846. [PMID: 38372002 DOI: 10.1021/acs.jctc.3c00799] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/20/2024]
Abstract
In this work, we present a new theory to treat multicomponent liquids based on quantum-chemically calculated clusters. The starting point is the binary quantum cluster equilibrium theory, which is able to treat binary systems. The theory provides one equation with two unknowns. In order to obtain another linearly independent equation, the conservation of mass is used. However, increasing the number of components leads to more unknowns, and this requires linearly independent equations. We address this challenge by introducing a generalization of the conservation of arbitrary quantities accompanied by a comprehensive mathematical proof. Furthermore, a case study for the application of the new theory to ternary mixtures of chloroform, methanol, and water is presented. Calculated enthalpies of vaporization for the whole composition range are given, and the populations or weights of the different clusters are visualized.
Collapse
Affiliation(s)
- Tom Frömbgen
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstraße 4 + 6, Bonn D-53115, Germany
- Max-Planck-Institut Für Chemische Energiekonversion, Stiftstrasse 34-36, Mülheim an der Ruhr D-45470, Germany
| | - Katrin Drysch
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstraße 4 + 6, Bonn D-53115, Germany
| | - Paul Zaby
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstraße 4 + 6, Bonn D-53115, Germany
| | - Jürgen Dölz
- Institute for Numerical Simulation, University of Bonn, Friedrich-Hirzebruch-Allee 7, Bonn D-53115, Germany
| | - Johannes Ingenmey
- CNRS, Physico-Chimie des Électrolytes et Nanosystèmes Interfaciaux, Sorbonne Université, Paris F-75005, France
| | - Barbara Kirchner
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstraße 4 + 6, Bonn D-53115, Germany
| |
Collapse
|
3
|
Taherivardanjani S, Wylie L, Dötzer R, Kirchner B. Exploring the Influence of the Phosphorus-Heteroatom Substitution in Nicotine on Its Electronic and Vibrational Spectroscopic Properties. Chemistry 2024; 30:e202302534. [PMID: 37984418 DOI: 10.1002/chem.202302534] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2023] [Revised: 11/10/2023] [Accepted: 11/13/2023] [Indexed: 11/22/2023]
Abstract
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.
Collapse
Affiliation(s)
- Shima Taherivardanjani
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, Beringstr. 4, D-53115, Bonn, Germany
| | - Luke Wylie
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, Beringstr. 4, D-53115, Bonn, Germany
| | | | - Barbara Kirchner
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, Beringstr. 4, D-53115, Bonn, Germany
| |
Collapse
|
4
|
Malloum A, Conradie J. Structures of DMSO clusters and quantum cluster equilibrium (QCE). J Mol Graph Model 2024; 126:108661. [PMID: 37913567 DOI: 10.1016/j.jmgm.2023.108661] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2023] [Revised: 10/25/2023] [Accepted: 10/26/2023] [Indexed: 11/03/2023]
Abstract
Dimethylsulfoxide (DMSO) clusters are crucial for understanding processes in liquid DMSO. Despite its importance, DMSO clusters have received negligible attention due to the complexity of their potential energy surfaces (PESs). In this work, we explored the PESs of the DMSO clusters from dimer to decamer, starting with classical molecular dynamics, followed by full optimizations at the PW6B95-D3/def2-TZVP level of theory. In addition, the binding energies, the binding enthalpy per DMSO, and the quantum theory of atoms in molecules (QTAIM) analysis of the most stable isomers are reported. Temperature effects on the stability of the isomers have also been assessed. After thoroughly exploring the PESs of the DMSO clusters, 159 configurations have been used to apply the quantum cluster equilibrium (QCE) theory to liquid DMSO. The quantum cluster equilibrium theory has been applied to determine the liquid properties of DMSO from DMSO clusters. Thus, using the QCE, the population of the liquid DMSO, its infrared spectrum, and some thermodynamic properties of the liquid DMSO are predicted. The QCE results show that the population of the liquid DMSO is mainly dominated by the DMSO dimer and decamer, with the contribution in trace of the DMSO monomer, trimer, tetramer, pentamer, and octamer. More interestingly, the predicted infrared spectrum of liquid DMSO is in qualitative agreement with the experiment.
Collapse
Affiliation(s)
- Alhadji Malloum
- Department of Chemistry, University of the Free State, PO BOX 339, Bloemfontein 9300, South Africa; Department of Physics, Faculty of Science, University of Maroua, PO BOX 46, Maroua, Cameroon.
| | - Jeanet Conradie
- Department of Chemistry, University of the Free State, PO BOX 339, Bloemfontein 9300, South Africa; Department of Chemistry, UiT - The Arctic University of Norway, N-9037 Tromsø, Norway
| |
Collapse
|
5
|
Blasius J, Drysch K, Pilz FH, Frömbgen T, Kielb P, Kirchner B. Efficient Prediction of Mole Fraction Related Vibrational Frequency Shifts. J Phys Chem Lett 2023; 14:10531-10536. [PMID: 37972218 DOI: 10.1021/acs.jpclett.3c02761] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2023]
Abstract
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.
Collapse
Affiliation(s)
- Jan Blasius
- Mulliken Center for Theoretical Chemistry, Clausius Institute of Physical and Theoretical Chemistry, University of Bonn, Beringstraße 4-6, D-53115 Bonn, Germany
| | - Katrin Drysch
- Mulliken Center for Theoretical Chemistry, Clausius Institute of Physical and Theoretical Chemistry, University of Bonn, Beringstraße 4-6, D-53115 Bonn, Germany
| | - Frank Hendrik Pilz
- Clausius Institute of Physical and Theoretical Chemistry, University of Bonn, Wegelerstraße 12, D-53115 Bonn, Germany
| | - Tom Frömbgen
- Mulliken Center for Theoretical Chemistry, Clausius Institute of Physical and Theoretical Chemistry, University of Bonn, Beringstraße 4-6, D-53115 Bonn, Germany
- Max-Planck-Institut für Chemische Energiekonversion, Stiftstrasse 34-36, D-45470 Mülheim an der Ruhr, Germany
| | - Patrycja Kielb
- Clausius Institute of Physical and Theoretical Chemistry, University of Bonn, Wegelerstraße 12, D-53115 Bonn, Germany
- Transdisciplinary Research Area "Building Blocks of Matter and Fundamental Interactions" (TRA Matter), University of Bonn, D-53115 Bonn, Germany
| | - Barbara Kirchner
- Mulliken Center for Theoretical Chemistry, Clausius Institute of Physical and Theoretical Chemistry, University of Bonn, Beringstraße 4-6, D-53115 Bonn, Germany
| |
Collapse
|
6
|
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
|
7
|
Taherivardanjani S, Blasius J, Brehm M, Dötzer R, Kirchner B. Conformer Weighting and Differently Sized Cluster Weighting for Nicotine and Its Phosphorus Derivatives. J Phys Chem A 2022; 126:7070-7083. [PMID: 36170053 DOI: 10.1021/acs.jpca.2c03133] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Weighting methods applied to systems with many conformers have been broadly employed to calculate thermodynamic properties, structural characteristics, and populations. To better understand and test the sensitivity of conventional weighting methods, the conformational distributions of nicotine and its phosphorus-substituted derivatives are investigated. The weighting schemes used for this are all based on Boltzmann statistics. Classical Boltzmann factors based on the electronic energy and the Gibbs free energy are calculated at different quantum chemical levels of theory and compared to cluster weights obtained by the quantum cluster equilibrium method. Furthermore, the influence of the modified rigid-rotor-harmonic-oscillator (mRRHO) approximation on the cluster weights is investigated. The substitution of the nitrogen atom in the methylpyrrolidine ring by a phosphorus atom results in more monomer conformers and clusters being populated. The conformational distribution of the monomers remained stable at different levels of theory and weighting methods. However, going to dimers and trimers, we observe a significant influence of the level of theory, weighting method, and mRRHO cutoff on the populations of these clusters. We show that mRRHO cutoff values of 50 and 100 cm-1 yield similar results, which is why 50 cm-1 is recommended as a robust choice. Furthermore, we observe that the global minimum for ΔE0 and ΔG varies in a few cases and that the global minimum is not always the dominantly occupied structure.
Collapse
Affiliation(s)
- Shima Taherivardanjani
- Mulliken Center for Theoretical Chemistry, Clausius Institute of Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4 + 6, D-53115 Bonn, Germany
| | - Jan Blasius
- Mulliken Center for Theoretical Chemistry, Clausius Institute of Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4 + 6, D-53115 Bonn, Germany
| | - Martin Brehm
- Institut für Chemie, Martin-Luther-Universität Halle-Wittenberg, D-06120 Halle (Saale), Germany
| | - Reinhard Dötzer
- Competence Center Analytics, BASF SE, D-67056 Ludwigshafen, Germany
| | - Barbara Kirchner
- Mulliken Center for Theoretical Chemistry, Clausius Institute of Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4 + 6, D-53115 Bonn, Germany
| |
Collapse
|
8
|
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] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [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
|
9
|
Krupová M, Leszczenko P, Sierka E, Emma Hamplová S, Pelc R, Andrushchenko V. Vibrational Circular Dichroism Unravels Supramolecular Chirality and Hydration Polymorphism of Nucleoside Crystals. Chemistry 2022; 28:e202201922. [DOI: 10.1002/chem.202201922] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2022] [Indexed: 11/07/2022]
Affiliation(s)
- Monika Krupová
- Institute of Organic Chemistry and Biochemistry Czech Academy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
- Hylleraas Centre for Quantum Molecular Sciences Department of Chemistry UiT The Arctic University of Norway N-9037 Tromsø Norway
| | - Patrycja Leszczenko
- Institute of Organic Chemistry and Biochemistry Czech Academy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
- Faculty of Chemistry Jagiellonian University Gronostajowa 2 30-387 Kraków Poland
| | - Ewa Sierka
- Institute of Organic Chemistry and Biochemistry Czech Academy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
- Faculty of Chemistry Jagiellonian University Gronostajowa 2 30-387 Kraków Poland
| | - Sára Emma Hamplová
- Institute of Organic Chemistry and Biochemistry Czech Academy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
- Department of Biology and Biochemistry University of Bath Claverton Down Bath BA2 7AY United Kingdom
| | - Radek Pelc
- Institute of Organic Chemistry and Biochemistry Czech Academy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
- Third Faculty of Medicine Charles University Ruská 87 10000 Prague Czech Republic
| | - Valery Andrushchenko
- Institute of Organic Chemistry and Biochemistry Czech Academy of Sciences Flemingovo náměstí 2 16610 Prague Czech Republic
| |
Collapse
|
10
|
Blasius J, Zaby P, Dölz J, Kirchner B. Uncertainty quantification of phase transition quantities from cluster weighting calculations. J Chem Phys 2022; 157:014505. [DOI: 10.1063/5.0093057] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In this work, we investigate how uncertainties in experimental input data influence the results of quantum cluster equilibrium calculations. In particular, we focus on the calculation of vaporization enthalpies and entropies of seven organic liquids, compare two computational approaches for their calculation and investigate how these properties are affected by changes in the experimental input data. It is observed that the vaporization enthalpies and entropies show a smooth dependence on changes in the reference density and boiling point. The reference density is found to have only a small influence on the vaporization thermodynamics, whereas the boiling point has a large influence on the vaporization enthalpy but only a small influence on the vaporization entropy. Furthermore we employed the Gauss--Hermite estimator in order to quantify the uncertainty in the thermodynamic functions that stems from inaccuracies in the experimental reference data at the example of the vaporization enthalpy of (\textit{R})-butan-2-ol. We quantify the uncertainty as 30.95~$\cdot$10$^{-3}$~kJ~mol$^{-1}$. Additionally we compare the convergence behaviour and computational effort of the Gauss--Hermite estimator with the Monte Carlo approach and show the superiority of the former. By this, we present how uncertainty quantification can be applied to examples from theoretical chemistry.
Collapse
Affiliation(s)
- Jan Blasius
- University of Bonn Institute of Physical and Theoretical Chemistry, Germany
| | - Paul Zaby
- University of Bonn Institute of Physical and Theoretical Chemistry, Germany
| | - Jürgen Dölz
- Institute for Numerical Simulation, University of Bonn, Friedrich-Hirzebruch-Allee 7 53115 Bonn, Germany, Germany
| | - Barbara Kirchner
- Theoretical and Physical Chemistry, University of Bonn Institute of Physical and Theoretical Chemistry, Germany
| |
Collapse
|
11
|
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
|
12
|
Al-Sheakh L, Fritsch S, Appelhagen A, Villinger A, Ludwig R. Thermodynamically Stable Cationic Dimers in Carboxyl-Functionalized Ionic Liquids: The Paradoxical Case of "Anti-Electrostatic" Hydrogen Bonding. Molecules 2022; 27:molecules27020366. [PMID: 35056680 PMCID: PMC8778807 DOI: 10.3390/molecules27020366] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2021] [Revised: 12/28/2021] [Accepted: 01/03/2022] [Indexed: 11/18/2022] Open
Abstract
We show that carboxyl-functionalized ionic liquids (ILs) form doubly hydrogen-bonded cationic dimers (c+=c+) despite the repulsive forces between ions of like charge and competing hydrogen bonds between cation and anion (c+–a−). This structural motif as known for formic acid, the archetype of double hydrogen bridges, is present in the solid state of the IL 1−(carboxymethyl)pyridinium bis(trifluoromethylsulfonyl)imide [HOOC−CH2−py][NTf2]. By means of quantum chemical calculations, we explored different hydrogen-bonded isomers of neutral (HOOC–(CH2)n–py+)2(NTf2−)2, single-charged (HOOC–(CH2)n–py+)2(NTf2−), and double-charged (HOOC– (CH2)n−py+)2 complexes for demonstrating the paradoxical case of “anti-electrostatic” hydrogen bonding (AEHB) between ions of like charge. For the pure doubly hydrogen-bonded cationic dimers (HOOC– (CH2)n−py+)2, we report robust kinetic stability for n = 1–4. At n = 5, hydrogen bonding and dispersion fully compensate for the repulsive Coulomb forces between the cations, allowing for the quantification of the two equivalent hydrogen bonds and dispersion interaction in the order of 58.5 and 11 kJmol−1, respectively. For n = 6–8, we calculated negative free energies for temperatures below 47, 80, and 114 K, respectively. Quantum cluster equilibrium (QCE) theory predicts the equilibria between cationic monomers and dimers by considering the intermolecular interaction between the species, leading to thermodynamic stability at even higher temperatures. We rationalize the H-bond characteristics of the cationic dimers by the natural bond orbital (NBO) approach, emphasizing the strong correlation between NBO-based and spectroscopic descriptors, such as NMR chemical shifts and vibrational frequencies.
Collapse
Affiliation(s)
- Loai Al-Sheakh
- Institut für Chemie, Abteilung für Physikalische Chemie, Universität Rostock, Dr.-Lorenz-Weg 2, 18059 Rostock, Germany; (L.A.-S.); (S.F.); (A.A.)
| | - Sebastian Fritsch
- Institut für Chemie, Abteilung für Physikalische Chemie, Universität Rostock, Dr.-Lorenz-Weg 2, 18059 Rostock, Germany; (L.A.-S.); (S.F.); (A.A.)
| | - Andreas Appelhagen
- Institut für Chemie, Abteilung für Physikalische Chemie, Universität Rostock, Dr.-Lorenz-Weg 2, 18059 Rostock, Germany; (L.A.-S.); (S.F.); (A.A.)
| | - Alexander Villinger
- Institut für Chemie, Abteilung für Anorganische Chemie, Universität Rostock, Albert-Einstein-Str. 3a, 18059 Rostock, Germany;
| | - Ralf Ludwig
- Institut für Chemie, Abteilung für Physikalische Chemie, Universität Rostock, Dr.-Lorenz-Weg 2, 18059 Rostock, Germany; (L.A.-S.); (S.F.); (A.A.)
- Department LL&M, University of Rostock, Albert-Einstein−Str. 25, 18059 Rostock, Germany
- Leibniz−Institut für Katalyse an der Universität Rostock e.V., Albert-Einstein−Str. 29a, 18059 Rostock, Germany
- Correspondence:
| |
Collapse
|
13
|
Marchelli G, Ingenmey J, Hollóczki O, Chaumont A, Kirchner B. Hydrogen Bonding and Vaporization Thermodynamics in Hexafluoroisopropanol-Acetone and -Methanol Mixtures. A Joined Cluster Analysis and Molecular Dynamic Study. Chemphyschem 2021; 23:e202100620. [PMID: 34632686 PMCID: PMC9298724 DOI: 10.1002/cphc.202100620] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2021] [Revised: 09/28/2021] [Indexed: 12/23/2022]
Abstract
Binary mixtures of hexafluoroisopropanol with either methanol or acetone are analyzed via classical molecular dynamics simulations and quantum cluster equilibrium calculations. In particular, their populations and thermodynamic properties are investigated with the binary quantum cluster equilibrium method, using our in‐house code peacemaker 2.8, upgraded with temperature‐dependent parameters. A novel approach, where the final density from classical molecular dynamics, has been used to generate the necessary reference isobars. The hydrogen bond network in both type of mixtures at molar fraction of hexafluoroisopropanol of 0.2, 0.5, and 0.8 respectively is investigated via the molecular dynamics trajectories and the cluster results. In particular, the populations show that mixed clusters are preferred in both systems even at 0.2 molar fractions of hexafluoroisopropanol. Enthalpies and entropies of vaporization are calculated for the neat and mixed systems and found to be in good agreement with experimental values.
Collapse
Affiliation(s)
- Gwydyon Marchelli
- Mulliken Center for Theoretical Chemistry, Rheinische Friedrich-Wilhelms-Universität Bonn, Beringstr. 4+6, D-53115, Bonn, Germany
| | - Johannes Ingenmey
- Mulliken Center for Theoretical Chemistry, Rheinische Friedrich-Wilhelms-Universität Bonn, Beringstr. 4+6, D-53115, Bonn, Germany
| | - Oldamur Hollóczki
- Mulliken Center for Theoretical Chemistry, Rheinische Friedrich-Wilhelms-Universität Bonn, Beringstr. 4+6, D-53115, Bonn, Germany
| | - Alain Chaumont
- Université de Strasbourg, CNRS, CMC UMR 7140, Laboratoire MSM, F-67000, Strasbourg, France
| | - Barbara Kirchner
- Mulliken Center for Theoretical Chemistry, Rheinische Friedrich-Wilhelms-Universität Bonn, Beringstr. 4+6, D-53115, Bonn, Germany
| |
Collapse
|
14
|
Zaby P, Ingenmey J, Kirchner B, Grimme S, Ehlert S. Calculation of improved enthalpy and entropy of vaporization by a modified partition function in quantum cluster equilibrium theory. J Chem Phys 2021; 155:104101. [PMID: 34525835 DOI: 10.1063/5.0061187] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023] Open
Abstract
In this work, we present an altered partition function that leads to an improved calculation of the enthalpy and entropy of vaporization in the framework of quantum cluster equilibrium theory. The changes are based on a previously suggested modification [S. Grimme, Chem. Eur. J. 18, 9955-9964 (2012)] of the molecular entropy calculation in the gas phase. Here, the low energy vibrational frequencies in the vibrational partition function are treated as hindered rotations instead of vibrations. The new scheme is tested on a set of nine organic solvents for the calculation of the enthalpy and entropy of vaporization. The enthalpies and entropies of vaporization show improvements from 6.5 error to 3.3 kJ mol-1 deviation to experiment and from 28.4 error to 13.5 J mol-1 K-1 deviation to experiment, respectively. The effect of the corrected partition function is visible in the different populations of clusters, which become physically more meaningful in that larger clusters are higher populated in the liquid phase and the gas phase is mainly populated by the monomers. Furthermore, the corrected partition function also overcomes technical difficulties and leads to an increased stability of the calculations in regard to the size of the cluster set.
Collapse
Affiliation(s)
- Paul Zaby
- Mulliken Center for Theoretical Chemistry, University of Bonn, Bonn, Germany
| | - Johannes Ingenmey
- Mulliken Center for Theoretical Chemistry, University of Bonn, Bonn, Germany
| | - Barbara Kirchner
- Mulliken Center for Theoretical Chemistry, University of Bonn, Bonn, Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, University of Bonn, Bonn, Germany
| | - Sebastian Ehlert
- Mulliken Center for Theoretical Chemistry, University of Bonn, Bonn, Germany
| |
Collapse
|
15
|
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
|
16
|
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
|
17
|
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
|