1
|
Khanifaev J, Schrader T, Perlt E. The effect of machine learning predicted anharmonic frequencies on thermodynamic properties of fluid hydrogen fluoride. J Chem Phys 2024; 160:124302. [PMID: 38516969 DOI: 10.1063/5.0195386] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/02/2024] [Accepted: 03/02/2024] [Indexed: 03/23/2024] Open
Abstract
Anharmonic effects play a crucial role in determining thermochemical properties of liquids and gases. For such extended phases, the inclusion of anharmonicity in reliable electronic structure methods is computationally extremely demanding, and hence, anharmonic effects are often lacking in thermochemical calculations. In this study, we apply the quantum cluster equilibrium method to transfer density functional theory calculations at the cluster level to the macroscopic, liquid, and gaseous phase of hydrogen fluoride. This allows us to include anharmonicity, either via vibrational self-consistent field calculations for smaller clusters or using a regression model for larger clusters. We obtain the structural composition of the fluid phases in terms of the population of different clusters as well as isobaric heat capacities as an example for thermodynamic properties. We study the role of anharmonicities for these analyses and observe that, in particular, the dominating structural motifs are rather sensitive to the anharmonicity in vibrational frequencies. The regression model proves to be a promising way to get access to anharmonic features, and the extension to more sophisticated machine-learning models is promising.
Collapse
Affiliation(s)
- Jamoliddin Khanifaev
- Otto Schott Institute of Materials Research, Friedrich Schiller University Jena, 07743 Jena, Germany
| | - Tim Schrader
- Otto Schott Institute of Materials Research, Friedrich Schiller University Jena, 07743 Jena, Germany
| | - Eva Perlt
- Otto Schott Institute of Materials Research, Friedrich Schiller University Jena, 07743 Jena, Germany
| |
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
|
de Andrade KN, Peixoto BP, Carneiro JWDM, Fiorot RG. Exploring borderline S N1-S N2 mechanisms: the role of explicit solvation protocols in the DFT investigation of isopropyl chloride. RSC Adv 2024; 14:4692-4701. [PMID: 38318615 PMCID: PMC10841197 DOI: 10.1039/d4ra00066h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/03/2024] [Accepted: 01/30/2024] [Indexed: 02/07/2024] Open
Abstract
Nucleophilic substitution at saturated carbon is a crucial class of organic reactions, playing a pivotal role in various chemical transformations that yield valuable compounds for society. Despite the well-established SN1 and SN2 mechanisms, secondary substrates, particularly in solvolysis reactions, often exhibit a borderline pathway. A molecular-level understanding of these processes is fundamental for developing more efficient chemical transformations. Typically, quantum-chemical simulations of the solvent medium combine explicit and implicit solvation methods. The configuration of explicit molecules can be defined through top-down approaches, such as Monte Carlo (MC) calculations for generating initial configurations, and bottom-up methods that involve user-dependent protocols to add solvent molecules around the substrate. Herein, we investigated the borderline mechanism of the hydrolysis of a secondary substrate, isopropyl chloride (iPrCl), at DFT-M06-2X/aug-cc-pVDZ level, employing explicit and explicit + implicit protocols. Top-down and bottom-up approaches were employed to generate substrate-solvent complexes of varying number (n = 1, 3, 5, 7, 9, and 12) and configurations of H2O molecules. Our findings consistently reveal that regardless of the solvation approach, the hydrolysis of iPrCl follows a loose-SN2-like mechanism with nucleophilic solvent assistance. Increasing the water cluster around the substrate in most cases led to reaction barriers of ΔH‡ ≈ 21 kcal mol-1, with nine water molecules from MC configurations sufficient to describe the reaction. The More O'Ferrall-Jencks plot demonstrates an SN1-like character for all transition state structures, showing a clear merged profile. The fragmentation activation strain analyses indicate that energy barriers are predominantly controlled by solvent-substrate interactions, supported by the leaving group stabilization assessed through CHELPG atomic charges.
Collapse
Affiliation(s)
- Karine Nascimento de Andrade
- Department of Organic Chemistry, Chemistry Institute, Universidade Federal Fluminense (UFF) Outeiro de São João Batista 24020-141 Niterói RJ Brazil
| | - Bárbara Pereira Peixoto
- Department of Organic Chemistry, Chemistry Institute, Universidade Federal Fluminense (UFF) Outeiro de São João Batista 24020-141 Niterói RJ Brazil
| | - José Walkimar de Mesquita Carneiro
- Department of Inorganic Chemistry, Chemistry Institute, Universidade Federal Fluminense (UFF) Outeiro de São João Batista 24020-141 Niterói RJ Brazil
| | - Rodolfo Goetze Fiorot
- Department of Organic Chemistry, Chemistry Institute, Universidade Federal Fluminense (UFF) Outeiro de São João Batista 24020-141 Niterói RJ Brazil
| |
Collapse
|
5
|
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
|
6
|
Hashim FH, Yu F, Izgorodina EI. Appropriate clusterset selection for the prediction of thermodynamic properties of liquid water with QCE theory. Phys Chem Chem Phys 2023; 25:9846-9858. [PMID: 36945858 DOI: 10.1039/d2cp03712b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/05/2023]
Abstract
Evident in many physical and chemical phenomena, thermodynamics is the study of how energy is stored, transformed and transferred in a molecule or material. However, prediction of these properties with simulation techniques is a non-trivial task as several factors such as composition and intermolecular interactions come into play. While molecular dynamics and ab initio molecular dynamics are the most common techniques for the prediction of thermodynamic properties, there exists many shortcomings associated with their use. Therefore, in this work we instead apply QCE theory to predict the thermodynamic properties of liquid water. This theory assumes that a condensed phase system can be represented as a 'mixture' of varying sized clusters rather than as a continuum. As QCE theory relies on first-principle simulations and statistical thermodynamics to determine the thermodynamic behavior of a system, appropriate selection of clusters is a crucial step towards achieving accurate predictions. In this study, we use molecular dynamics and ab initio calculations to obtain initial configurations of 400 water clusters, Wn where n = 3 to 10 and contrast their stability using two different criteria. The role of entropy towards cluster stabilization is investigated by comparing the binding (ΔEBIND/mol) and Gibbs free binding energy per molecule (ΔGBIND/mol) of various Wn at 298.15 K. Initial clustersets are constructed by exploring two-, three-, four and five-combinations of clustersets using the minimum ΔGBIND/mol structures of Wn. We also expand the ΔGBIND/mol criteria for Wn of sizes 3 to 7 to include values larger than 0.0 kJ mol-1 and smaller than 3.0 kJ mol-1 as a means of improving thermodynamic predictions. 37 of the 459 resulting clustersets predicted the correct boiling point of water and its thermodynamic properties with an accuracy of 95%. A scaled population-weighted infrared spectrum was compared to experimental results to validate the composition of the top 5 clustersets. The predicted spectra showed an exact match within the low frequency range (<1000 cm-1) with some discrepancy at the high frequency range (>3400 cm-1). This work highlights that ΔGBIND/mol is so far the best criteria to apply when determining an appropriate clusterset for QCE theory.
Collapse
Affiliation(s)
- Fairuz H Hashim
- 17 Rainforest Walk, School of Chemistry, Monash University, Australia.
| | - Fiona Yu
- 17 Rainforest Walk, School of Chemistry, Monash University, Australia.
| | | |
Collapse
|
7
|
Tang H, Cai J, Zhu CY, Chen GJ, Wang XH, Sun CY. Review on the clustering behavior in aqueous solutions. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.120382] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
8
|
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
|
9
|
Hunger L, Al‐Sheakh L, Zaitsau DH, Verevkin SP, Appelhagen A, Villinger A, Ludwig R. Dissecting Noncovalent Interactions in Carboxyl‐Functionalized Ionic Liquids Exhibiting Double and Single Hydrogens Bonds Between Ions of Like Charge. Chemistry 2022; 28:e202200949. [PMID: 35785500 PMCID: PMC9543318 DOI: 10.1002/chem.202200949] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2022] [Indexed: 11/11/2022]
Affiliation(s)
- Lasse Hunger
- Institut für Chemie Abteilung für Physikalische Chemie Universität Rostock Albert-Einstein-Str. 27 18059 Rostock Germany
| | - Loai Al‐Sheakh
- Institut für Chemie Abteilung für Physikalische Chemie Universität Rostock Albert-Einstein-Str. 27 18059 Rostock Germany
| | - Dzmitry H. Zaitsau
- Institut für Chemie Abteilung für Physikalische Chemie Universität Rostock Albert-Einstein-Str. 27 18059 Rostock Germany
- Department LL&M University of Rostock Albert-Einstein-Str. 25 18059 Rostock Germany
| | - Sergey P. Verevkin
- Institut für Chemie Abteilung für Physikalische Chemie Universität Rostock Albert-Einstein-Str. 27 18059 Rostock Germany
- Department LL&M University of Rostock Albert-Einstein-Str. 25 18059 Rostock Germany
| | - Andreas Appelhagen
- Institut für Chemie Abteilung für Physikalische Chemie Universität Rostock Albert-Einstein-Str. 27 18059 Rostock Germany
| | - 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 Albert-Einstein-Str. 27 18059 Rostock Germany
- 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
| |
Collapse
|
10
|
Weinhold F. High-Density “Windowpane” Coordination Patterns of Water Clusters and Their NBO/NRT Characterization. Molecules 2022; 27:molecules27134218. [PMID: 35807463 PMCID: PMC9268199 DOI: 10.3390/molecules27134218] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2022] [Revised: 06/25/2022] [Accepted: 06/27/2022] [Indexed: 02/05/2023] Open
Abstract
Cluster mixture models for liquid water at higher pressures suggest the need for water clusters of higher coordination and density than those commonly based on tetrahedral H-bonding motifs. We show here how proton-ordered water clusters of increased coordination and density can assemble from a starting cyclic tetramer or twisted bicyclic (Möbius-like) heptamer to form extended Aufbau sequences of stable two-, three-, and four-coordinate “windowpane” motifs. Such windowpane clusters exhibit sharply reduced (~90°) bond angles that differ appreciably from the tetrahedral angles of idealized crystalline ice Ih. Computed free energy and natural resonance theory (NRT) bond orders provide quantitative descriptors for the relative stabilities of clusters and strengths of individual coordinative linkages. The unity and consistency of NRT description is demonstrated to extend from familiar supra-integer bonds of the molecular regime to the near-zero bond orders of the weakest linkages in the present H-bond clusters. Our results serve to confirm that H-bonding exemplifies resonance–covalent (fractional) bonding in the sub-integer range and to further discount the dichotomous conceptions of “electrostatics” for intermolecular bonding vs. “covalency” for intramolecular bonding that still pervade much of freshman-level pedagogy and force-field methodology.
Collapse
Affiliation(s)
- Frank Weinhold
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, Madison, WI 53706, USA
| |
Collapse
|
11
|
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
|
12
|
Kirchner B, Ingenmey J, von Domaros M, Perlt E. The Ionic Product of Water in the Eye of the Quantum Cluster Equilibrium. Molecules 2022; 27:molecules27041286. [PMID: 35209075 PMCID: PMC8877775 DOI: 10.3390/molecules27041286] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2021] [Revised: 02/06/2022] [Accepted: 02/09/2022] [Indexed: 11/16/2022] Open
Abstract
The theoretical description of water properties continues to be a challenge. Using quantum cluster equilibrium (QCE) theory, we combine state-of-the-art quantum chemistry and statistical thermodynamic methods with the almost historical Clausius-Clapeyron relation to study water self-dissociation and the thermodynamics of vaporization. We pay particular attention to the treatment of internal rotations and their impact on the investigated properties by employing the modified rigid-rotor-harmonic-oscillator (mRRHO) approach. We also study a novel QCE parameter-optimization procedure. Both the ionic product and the vaporization enthalpy yield an astonishing agreement with experimental reference data. A significant influence of the mRRHO approach is observed for cluster populations and, consequently, for the ionic product. Thermodynamic properties are less affected by the treatment of these low-frequency modes.
Collapse
Affiliation(s)
- Barbara Kirchner
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, Beringstr. 4, 53115 Bonn, Germany
- Correspondence:
| | - Johannes Ingenmey
- CNRS, Physico-Chimie des Électrolytes et Nanosystèmes Interfaciaux, Sorbonne Université, F-75005 Paris, France;
| | - Michael von Domaros
- Department of Chemistry, Philipps-Universität Marburg, Hans-Meerwein-Straße 4, 35032 Marburg, Germany;
| | - Eva Perlt
- Otto Schott Institute of Materials Research, Faculty of Physics and Astronomy, Friedrich-Schiller-Universität Jena, Löbdergraben 32, 07743 Jena, Germany;
| |
Collapse
|
13
|
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
|
14
|
Kiefer PM, Daschakraborty S, Pines D, Pines E, Hynes JT. Electron Flow Characterization of Charge Transfer for Carbonic Acid to Strong Base Proton Transfer in Aqueous Solution. J Phys Chem B 2021; 125:11473-11490. [PMID: 34623157 DOI: 10.1021/acs.jpcb.1c05824] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Protonation of the strong base methylamine CH3NH2 by carbonic acid H2CO3 in aqueous solution, HOCOOH···NH2CH3 → HOCOO-···+HNH2CH3, has been previously studied ( J. Phys. Chem. B 2016, 109, 2271-2280; J. Phys. Chem. B 2016, 109, 2281-2290) via Car-Parinnello molecular dynamics. This proton transfer (PT) reaction within a hydrogen (H)-bonded complex was found to be barrierless and very rapid, with key reaction coordinates comprising the proton coordinate, the H-bond separation RON, and a solvent coordinate, reflecting the water solvent rearrangement involved in the neutral to ion pair conversion. In the present work, the reaction's charge flow aspects are analyzed in detail, especially a description via Mulliken charge transfer for PT (MCTPT). A natural bond orbital analysis and some extensions of them are employed for the complex's electronic structure during the reaction trajectories. Results demonstrate that consistent with the MCTPT picture, the charge transfer (CT) occurs from a methylamine base nonbonding orbital to a carbonic acid antibonding orbital. A complementary MCTPT reaction product perspective of CT from the antibonding orbital of the HN+ moiety to the nonbonding orbital of the oxygen in the H-bond complex is also presented. σOH and σHN+ bond order expressions show this CT to occur within the H-bond OHN triad, an aspect key for simultaneous bond-breaking and -forming in the PT reaction.
Collapse
Affiliation(s)
- Philip M Kiefer
- Department of Chemistry, University of Colorado, Boulder, Colorado 80309-0215, United States
| | - Snehasis Daschakraborty
- Department of Chemistry, University of Colorado, Boulder, Colorado 80309-0215, United States
| | - Dina Pines
- Department of Chemistry, Ben-Gurion University of the Negev, P.O. Box 653, Beer-Sheva 84105, Israel
| | - Ehud Pines
- Department of Chemistry, Ben-Gurion University of the Negev, P.O. Box 653, Beer-Sheva 84105, Israel
| | - James T Hynes
- Department of Chemistry, University of Colorado, Boulder, Colorado 80309-0215, United States.,PASTEUR, Department of Chemistry, École Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| |
Collapse
|
15
|
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: 6] [Impact Index Per Article: 2.0] [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
|
16
|
Ludwig R. Towards thermodynamically stable anionic dimers with “anti-electrostatic” hydrogen bonds overcoming like-charge repulsion. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.116882] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
17
|
Lemke KH. Structure and solvation dynamics of the hydroxide ion in ice-like water clusters: a CCSD(T) and car-parrinello molecular dynamics study. Phys Chem Chem Phys 2021; 23:18990-18998. [PMID: 34612437 DOI: 10.1039/d1cp02524d] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Using MP2, CCSD(T) electronic structure theory and ab initio molecular dynamics simulations, we explore the structure, solvation dynamics and vibrational spectra of OH-(H2O)n clusters. Our study reports new cubic and fused cubic global minima structures of OH-(H2O)n for n = 8-26 with surface and interior solvation arrangements. In the case of OH-(H2O)26, we show that MP2 and CCSD(T) calculations predict global minima structures with the hydroxide ion occupying the interior region of a densely packed cubic cluster that is secured by ionic hydrogen bonds. More importantly, results from ab initio molecular dynamics simulations of OH-(H2O)26 demonstrate that the hydroxide ion remains in the cluster interior and hexa-coordinated, irrespective of the temperature, up to around 175 K, then incrementally transitions from a surface-exposed penta- (170-200 K), to a tetra- (225 K) to a tri-coordinated OH-(H2O)3 structure at 300 K. Building on our temperature-dependent vibrational power spectra, we are also able to disentangle structure and temperature effects on individual spectral contributions arising from water molecules located in the inner and outer shell of OH-(H2O)26. Some of these theoretical results provide valuable guidance for the interpretation of IRMPD spectra of small hydroxide-water clusters, but there are also several intriguing implications of these results, in particular, for the solvation of the OH- ion at the surface of water nanodroplets and aqueous interfaces.
Collapse
Affiliation(s)
- Kono H Lemke
- Department of Earth Sciences, University of Hong Kong, Pokfulam Road, Hong Kong, SAR.
| |
Collapse
|
18
|
Gehrke S, Ray P, Stettner T, Balducci A, Kirchner B. Water in Protic Ionic Liquid Electrolytes: From Solvent Separated Ion Pairs to Water Clusters. CHEMSUSCHEM 2021; 14:3315-3324. [PMID: 34169663 PMCID: PMC8456901 DOI: 10.1002/cssc.202100660] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 03/30/2021] [Revised: 06/21/2021] [Indexed: 06/13/2023]
Abstract
The large electrochemical and cycling stability of "water-in-salt" systems have rendered promising prospective electrolytes for batteries. The impact of addition of water on the properties of ionic liquids has already been addressed in several publications. In this contribution, we focus on the changes in the state of water. Therefore, we investigated the protic ionic liquid N-butyl-pyrrolidinium bis(trifluoromethanesulfonyl)imide with varying water content at different temperatures with the aid of molecular dynamics simulations. It is revealed that at very low concentrations, the water is well dispersed and best characterized as shared solvent molecules. At higher concentrations, the water forms larger aggregates and is increasingly approaching a bulk-like state. While the librational and rotational dynamics of the water molecules become faster with increasing concentration, the translational dynamics are found to become slower. Further, all dynamics are found to be faster if the temperature increases. The trends of these findings are well in line with the experimental measured conductivities.
Collapse
Affiliation(s)
- Sascha Gehrke
- Mulliken Center for Theoretical ChemistryUniversity of BonnBeringstr. 4+6D-53115BonnGermany
- Department of Physics and AstronomyUniversity College LondonLondonWC1E 6BTUnited Kingdom
| | - Promit Ray
- Mulliken Center for Theoretical ChemistryUniversity of BonnBeringstr. 4+6D-53115BonnGermany
| | - Timo Stettner
- Institute for Technical Chemistry and Environmental ChemistryFriedrich-Schiller-University JenaPhilosophenweg 7aD-07743JenaGermany
- Center for Energy and Environmental Chemistry Jena (CEEC Jena)Friedrich-Schiller-University JenaPhilosophenweg 7aD-07743JenaGermany
| | - Andrea Balducci
- Institute for Technical Chemistry and Environmental ChemistryFriedrich-Schiller-University JenaPhilosophenweg 7aD-07743JenaGermany
- Center for Energy and Environmental Chemistry Jena (CEEC Jena)Friedrich-Schiller-University JenaPhilosophenweg 7aD-07743JenaGermany
| | - Barbara Kirchner
- Mulliken Center for Theoretical ChemistryUniversity of BonnBeringstr. 4+6D-53115BonnGermany
| |
Collapse
|
19
|
Affiliation(s)
- Frank Weinhold
- Theoretical Chemistry and Department of Chemistry University of Wisconsin Madison Wisconsin 53706 USA
| |
Collapse
|
20
|
Thomsen B, Shiga M. Nuclear quantum effects on autoionization of water isotopologs studied by ab initio path integral molecular dynamics. J Chem Phys 2021; 154:084117. [PMID: 33639728 DOI: 10.1063/5.0040791] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
In this study, we investigate the nuclear quantum effects (NQEs) on the acidity constant (pKA) of liquid water isotopologs under the ambient condition by path integral molecular dynamics (PIMD) simulations. We compared simulations using a fully explicit solvent model with a classical polarizable force field, density functional tight binding, and ab initio density functional theory, which correspond to empirical, semiempirical, and ab initio PIMD simulations, respectively. The centroid variable with respect to the proton coordination number of a water molecule was restrained to compute the gradient of the free energy, which measures the reversible work of the proton abstraction for the quantum mechanical system. The free energy curve obtained by thermodynamic integration was used to compute the pKA value based on probabilistic determination. This technique not only reproduces the pKA value of liquid D2O experimentally measured (14.86) but also allows for a theoretical prediction of the pKA values of liquid T2O and aqueous HDO and HTO, which are unknown due to their scarcity. It is also shown that the NQEs on the free energy curve can result in a downshift of 4.5 ± 0.9 pKA units in the case of liquid water, which indicates that the NQEs plays an indispensable role in the absolute determination of pKA. The results of this study can help inform further extensions into the calculation of the acidity constants of isotope substituted species with high accuracy.
Collapse
Affiliation(s)
- Bo Thomsen
- CCSE, Japan Atomic Energy Agency, 178-4-4, Wakashiba, Kashiwa, Chiba 277-0871, Japan
| | - Motoyuki Shiga
- CCSE, Japan Atomic Energy Agency, 178-4-4, Wakashiba, Kashiwa, Chiba 277-0871, Japan
| |
Collapse
|
21
|
Yap ACW, Lee HS, Loo JL, Mohd NS. Electron generation in water induced by magnetic effect and its impact on dissolved oxygen concentration. SUSTAINABLE ENVIRONMENT RESEARCH 2021; 31:7. [PMCID: PMC7875160 DOI: 10.1186/s42834-021-00080-0] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/19/2020] [Accepted: 01/28/2021] [Indexed: 06/15/2023]
Abstract
pH, oxidation-reduction potential (ORP) and dissolved oxygen (DO) concentration are important parameters in water quality surveillance and treatment. The changes of these parameters are associated with electron density in water. Several techniques including electrolysis and catalysis which require redox reactions and electron exchange are employed to improve these parameters. In recent years, studies reported that magnetic effects can impart considerable changes on the pH, ORP and DO concentration of water. However, the correlation between electron density and magnetic effects on these parameters has yet to be disclosed despite the fact that increased electron density in water could improve water’s reductive properties, heat capacity and hydrogen bonding characteristics. In this study, the magnetic effects on pH, ORP and DO concentration were investigated using different magnets arrangements and water flow rates based on reversed electric motor principle. Results showed that the improvement of pH, ORP and DO concentration from 5.40–5.42 to 5.58–5.62 (+ 3.5%), 392 to 365 mV (− 6.9%), and 7.30 to 7.71 mg L− 1 (+ 5.6%), respectively were achieved using combined variables of non-reversed polarity magnet arrangement (1000–1500 G magnetic strength) and water flow rate of 0.1–0.5 mL s− 1. Such decrement in ORP value also corresponded to 8.0 × 1013 number of electron generation in water. Furthermore, Raman analysis revealed that magnetic effect could strengthen the intermolecular hydrogen bonding of water molecules and favor formation of smaller water clusters. The findings of this study could contribute to potential applications in aquaculture, water quality control and treatment of cancer attributed to free radical induced-oxidative stress.
Collapse
Affiliation(s)
- Augustine Chung Wei Yap
- Department of Mechanical and Material Engineering, Universiti Tunku Abdul Rahman, 43000 Kajang, Malaysia
| | - Hwang Sheng Lee
- Department of Mechanical and Material Engineering, Universiti Tunku Abdul Rahman, 43000 Kajang, Malaysia
- Centre for Photonics and Advanced Materials Research, Universiti Tunku Abdul Rahman, 43000 Kajang, Malaysia
| | - Joo Ling Loo
- Department of Mechatronics and Biomedical Engineering, Universiti Tunku Abdul Rahman, 43000 Kajang, Malaysia
| | - Nuruol Syuhadaa Mohd
- Department of Civil Engineering, University of Malaya, 50603 Kuala Lumpur, Malaysia
| |
Collapse
|
22
|
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
|
23
|
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
|
24
|
Marchelli G, Ingenmey J, Kirchner B. Activity coefficients of binary methanol alcohol mixtures from cluster weighting. ChemistryOpen 2020; 9:774-785. [PMID: 32714740 PMCID: PMC7377217 DOI: 10.1002/open.202000171] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2020] [Revised: 06/16/2020] [Indexed: 11/10/2022] Open
Abstract
The hydrogen bond network of different small alcohols is investigated via cluster analysis. Methanol/alcohol mixtures are studied with increasing chain length and branching of the molecule. Those changes can play an important role in different fields, including solvent and metal extraction. The extended tight binding method GFN2-xTB allows the evaluation and geometry optimization of thousands of clusters built via a genetic algorithm. Interaction energies and geometries are evaluated and discussed for the neat systems. Thermodynamic properties, such as vaporization enthalpies and activity coefficients, are calculated with the binary quantum cluster equilibrium (bQCE) approach using our in-house code peacemaker 2.8. Combined distribution functions of the distances against the angles of the hydrogen bonds are evaluated for neat and mixed clusters and weighted by the equilibrium populations achieved from bQCE calculations.
Collapse
Affiliation(s)
- Gwydyon Marchelli
- Mulliken Center for Theoretical ChemistryRheinische Friedrich-Wilhelms-Universität BonnBeringstr. 4+6D-53115BonnGermany
| | - J. Ingenmey
- Mulliken Center for Theoretical ChemistryRheinische Friedrich-Wilhelms-Universität BonnBeringstr. 4+6D-53115BonnGermany
| | - B. Kirchner
- Mulliken Center for Theoretical ChemistryRheinische Friedrich-Wilhelms-Universität BonnBeringstr. 4+6D-53115BonnGermany
| |
Collapse
|
25
|
Turi L, Rodriguez J, Laria D. Combined Effects from Solvation and Nuclear Quantum Fluctuations on Autoionization Mechanisms in Aqueous Clusters. J Phys Chem B 2020; 124:2198-2208. [PMID: 32075372 DOI: 10.1021/acs.jpcb.9b11087] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Using path-integral molecular dynamics simulations, we examine isomerization paths involving collective proton transfers in [H2O]5 and [H2O]8 clusters under cryogenic conditions. We focused attention on combined effects derived from solvation and nuclear quantum fluctuations on the characteristics of free energy barriers and relative stabilities of reactants and products. In particular, we analyzed two different processes: the first one involves the exchange of donor-acceptor hydrogen bond roles along cyclic moieties, whereas the second one corresponds to charge separation leading to stable [H3O]+[OH]- ion pairs. In the first case, the explicit incorporation of quantum tunneling introduces important modifications in the classical free energy profile. The resulting quantum profile presents two main contributions, one corresponding to compressions of O-O distances and a second one ascribed to nuclear tunneling of the light protons. Solvation effects promote a moderate polarization of the cyclic structures and a partial loss of concertedness in the collective modes, most notably, at the onset of tunneling. Still, the latter effects are also sufficiently strong to promote the stabilization of ion pairs along the classical trajectories. In contrast, the explicit incorporation of nuclear quantum fluctuations leads to charge separated configurations that are marginally stable. As such, the latter states could also be regarded as short-lived intermediate states along the reactive exchange path.
Collapse
Affiliation(s)
- László Turi
- Institute of Chemistry, Eötvös Loránd University, P.O. Box 32, Budapest, 112 H-1518, Hungary
| | - Javier Rodriguez
- Departamento de Fı́sica de la Materia Condensada, Comisión Nacional de Energı́a Atómica, Avenida Libertador 8250, 1429 Buenos Aires, Argentina.,ECyT, UNSAM, Martı́n de Irigoyen 3100, 1650 San Martı́n, Provincia de Buenos Aires, Argentina
| | - Daniel Laria
- Departamento de Fı́sica de la Materia Condensada, Comisión Nacional de Energı́a Atómica, Avenida Libertador 8250, 1429 Buenos Aires, Argentina.,Departamento de Quı́mica Inorgánica, Analı́tica y Quı́mica-Fı́sica and INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires Ciudad Universitaria, Pabellón II, 1428 Buenos Aires, Argentina
| |
Collapse
|
26
|
Pongratz T, Kibies P, Eberlein L, Tielker N, Hölzl C, Imoto S, Beck Erlach M, Kurrmann S, Schummel PH, Hofmann M, Reiser O, Winter R, Kremer W, Kalbitzer HR, Marx D, Horinek D, Kast SM. Pressure-dependent electronic structure calculations using integral equation-based solvation models. Biophys Chem 2020; 257:106258. [DOI: 10.1016/j.bpc.2019.106258] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2019] [Accepted: 08/25/2019] [Indexed: 12/18/2022]
|
27
|
Glendening E, Weinhold F. Natural resonance theory of chemical reactivity, with illustrative application to intramolecular Claisen rearrangement. Tetrahedron 2018. [DOI: 10.1016/j.tet.2018.07.054] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
28
|
Ingenmey J, von Domaros M, Perlt E, Verevkin SP, Kirchner B. Thermodynamics and proton activities of protic ionic liquids with quantum cluster equilibrium theory. J Chem Phys 2018; 148:193822. [PMID: 30307228 DOI: 10.1063/1.5010791] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Johannes Ingenmey
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstr. 4, D-53115 Bonn, Germany
| | - Michael von Domaros
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstr. 4, D-53115 Bonn, Germany
| | - Eva Perlt
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstr. 4, D-53115 Bonn, Germany
| | - Sergey P. Verevkin
- Abteilung Physikalische Chemie, Universität Rostock, Dr.-Lorenz-Weg 2, D-18059 Rostock, Germany
| | - Barbara Kirchner
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstr. 4, D-53115 Bonn, Germany
| |
Collapse
|
29
|
Weinhold F. Theoretical Prediction of Robust Second-Row Oxyanion Clusters in the Metastable Domain of Antielectrostatic Hydrogen Bonding. Inorg Chem 2018; 57:2035-2044. [PMID: 29381336 DOI: 10.1021/acs.inorgchem.7b02943] [Citation(s) in RCA: 33] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We provide ab initio and density functional theory evidence for a family of surprisingly robust like-charged clusters of common HSO4- and H2PO4- oxyanions, ranging up to tetramers of net charge 4-. Our results support other recent theoretical and experimental evidence for "antielectrostatic" hydrogen-bonded (AEHB) species that challenge conventional electrostatic conceptions and force-field modeling of closed-shell ion interactions. We provide structural and energetic descriptors of the predicted kinetic well-depths (in the range 3-10 kcal/mol) and barrier widths (in the range 2-4 Å) for simple AEHB dimers, including evidence of extremely strong hydrogen bonding in the fluoride-bisulfate dianion. For more complex polyanionic species, we employ natural-bond-orbital-based descriptors to characterize the electronic features of the cooperative hydrogen-bonding network that are able to successfully defy Coulomb explosion. The computational results suggest a variety of kinetically stable AEHB species that may be suitable for experimental detection as long-lived gas-phase species or structural units of condensed phases, despite the imposing electrostatic barriers that oppose their formation under ambient conditions.
Collapse
Affiliation(s)
- Frank Weinhold
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison , Madison, Wisconsin 53706, United States
| |
Collapse
|
30
|
Weinhold F. Polyion Covalency: Exotic Species from the Unexplored World of Electrostatically Shielded Molecular Ion Chemistry. Angew Chem Int Ed Engl 2017. [DOI: 10.1002/ange.201708691] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Affiliation(s)
- Frank Weinhold
- Theoretical Chemistry Institute and Department of Chemistry; University of Wisconsin; Madison Wisconsin 53706 USA
| |
Collapse
|
31
|
Weinhold F. Polyion Covalency: Exotic Species from the Unexplored World of Electrostatically Shielded Molecular Ion Chemistry. Angew Chem Int Ed Engl 2017; 56:14577-14581. [PMID: 28977727 DOI: 10.1002/anie.201708691] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2017] [Indexed: 12/21/2022]
Abstract
Standard quantum chemical methods have been employed to describe a variety of kinetically stable polyionic molecular species that are trapped in appreciable potential wells by chemical bonding forces, despite powerful electrostatic opposition that challenges conventional chemical detection and characterization. The studied species are covalent or dative analogs of "anti-electrostatic" hydrogen-bonded (AEHB) species, all illustrating how short-range quantum covalency can overcome the powerful "shielding" opposition of long-range electrostatic forces to form highly charged molecular species, analogous to known neutral or singly ionic counterparts. Computational predictions of representative structural, spectroscopic, and NBO-based electronic signatures of multiply charged analogs of common neutral species (CH3 CH3 , CO2 , FeCO) are provided to suggest the unique material properties characteristic of this shielded domain of polyionic chemical phenomena.
Collapse
Affiliation(s)
- Frank Weinhold
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin, Madison, Wisconsin, 53706, USA
| |
Collapse
|