1
|
Rojas-González FE, Castillo-Quevedo C, Rodríguez-Kessler PL, Jimenez-Halla JOC, Vásquez-Espinal A, Eithiraj RD, Cortez-Valadez M, Cabellos JL. Exploration of Free Energy Surface of the Au 10 Nanocluster at Finite Temperature. Molecules 2024; 29:3374. [PMID: 39064952 PMCID: PMC11279810 DOI: 10.3390/molecules29143374] [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: 05/26/2024] [Revised: 07/08/2024] [Accepted: 07/15/2024] [Indexed: 07/28/2024] Open
Abstract
The first step in comprehending the properties of Au10 clusters is understanding the lowest energy structure at low and high temperatures. Functional materials operate at finite temperatures; however, energy computations employing density functional theory (DFT) methodology are typically carried out at zero temperature, leaving many properties unexplored. This study explored the potential and free energy surface of the neutral Au10 nanocluster at a finite temperature, employing a genetic algorithm coupled with DFT and nanothermodynamics. Furthermore, we computed the thermal population and infrared Boltzmann spectrum at a finite temperature and compared it with the validated experimental data. Moreover, we performed the chemical bonding analysis using the quantum theory of atoms in molecules (QTAIM) approach and the adaptive natural density partitioning method (AdNDP) to shed light on the bonding of Au atoms in the low-energy structures. In the calculations, we take into consideration the relativistic effects through the zero-order regular approximation (ZORA), the dispersion through Grimme's dispersion with Becke-Johnson damping (D3BJ), and we employed nanothermodynamics to consider temperature contributions. Small Au clusters prefer the planar shape, and the transition from 2D to 3D could take place at atomic clusters consisting of ten atoms, which could be affected by temperature, relativistic effects, and dispersion. We analyzed the energetic ordering of structures calculated using DFT with ZORA and single-point energy calculation employing the DLPNO-CCSD(T) methodology. Our findings indicate that the planar lowest energy structure computed with DFT is not the lowest energy structure computed at the DLPN0-CCSD(T) level of theory. The computed thermal population indicates that the 2D elongated hexagon configuration strongly dominates at a temperature range of 50-800 K. Based on the thermal population, at a temperature of 100 K, the computed IR Boltzmann spectrum agrees with the experimental IR spectrum. The chemical bonding analysis on the lowest energy structure indicates that the cluster bond is due only to the electrons of the 6 s orbital, and the Au d orbitals do not participate in the bonding of this system.
Collapse
Affiliation(s)
| | - César Castillo-Quevedo
- Departamento de Fundamentos del Conocimiento, Centro Universitario del Norte, Universidad de Guadalajara, Carretera Federal No. 23, km. 191, Colotlán 46200, Jalisco, Mexico;
| | | | - José Oscar Carlos Jimenez-Halla
- Departamento de Química, División de Ciencias Exactas y Naturales, Universidad de Guanajuato, Noria Alta s/n, Guanajuato 36050, Guanajuato, Mexico;
| | - Alejandro Vásquez-Espinal
- Química y Farmacia, Facultad de Ciencias de la Salud, Universidad Arturo Prat. Casilla 121, Iquique 1100000, Chile;
| | | | - Manuel Cortez-Valadez
- CONAHCYT-Departamento de Investigación en Física, Universidad de Sonora, Apdo. Postal 5-88, Hermosillo 83190, Sonora, Mexico;
| | - José Luis Cabellos
- Coordinación de Investigación y Desarrollo Tecnológico, Universidad Politécnica de Tapachula, Carretera Tapachula a Puerto Madero km. 24, Tapachula 30830, Chiapas, Mexico
| |
Collapse
|
2
|
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
|
3
|
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
|
4
|
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
|
5
|
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
|
6
|
Van Benschoten W, Petras HR, Shepherd JJ. Electronic Free Energy Surface of the Nitrogen Dimer Using First-Principles Finite Temperature Electronic Structure Methods. J Phys Chem A 2023; 127:6842-6856. [PMID: 37535315 PMCID: PMC10440793 DOI: 10.1021/acs.jpca.3c01741] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2023] [Revised: 07/20/2023] [Indexed: 08/04/2023]
Abstract
We use full configuration interaction and density matrix quantum Monte Carlo methods to calculate the electronic free energy surface of the nitrogen dimer within the free-energy Born-Oppenheimer approximation. As the temperature is raised from T = 0, we find a temperature regime in which the internal energy causes bond strengthening. At these temperatures, adding in the entropy contributions is required to cause the bond to gradually weaken with increasing temperature. We predict a thermally driven dissociation for the nitrogen dimer between 22,000 to 63,200 K depending on symmetries and basis set. Inclusion of more spatial and spin symmetries reduces the temperature required. The origin of these observations is explored using the structure of the density matrix at various temperatures and bond lengths.
Collapse
Affiliation(s)
| | - Hayley R. Petras
- Department
of Chemistry, University of Iowa, Iowa City, Iowa 52242, United States
| | - James J. Shepherd
- Department
of Chemistry, University of Iowa, Iowa City, Iowa 52242, United States
| |
Collapse
|
7
|
Hu Y, Wang K, Müller M, Wessel E, Spatschek R. Theoretical Prediction of the Sublimation Behavior by Combining Ab Initio Calculations with Statistical Mechanics. MATERIALS (BASEL, SWITZERLAND) 2023; 16:2826. [PMID: 37049120 PMCID: PMC10095645 DOI: 10.3390/ma16072826] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/28/2023] [Revised: 03/29/2023] [Accepted: 03/30/2023] [Indexed: 06/19/2023]
Abstract
We develop a theoretical model to predict the sublimation vapor pressure of pure substances. Moreover, we present a simple monoatomic molecule approximation, which reduces the complexity of the vapor pressure expression for polyatomic gaseous molecules at a convincing level of accuracy, with deviations of the Arrhenius prefactor for NaCl and NaF being 5.02% and 7.08%, respectively. The physical model is based on ab initio calculations, statistical mechanics, and thermodynamics. We illustrate the approach for Ni, Cr, Cu (metallic bond), NaCl, NaF, ZrO2 (ionic bond) and SiO2 (covalent bond). The results are compared against thermodynamic databases, which show high accuracy of our theoretical predictions, and the deviations of the predicted sublimation enthalpy are typically below 10%, for Cu even only 0.1%. Furthermore, the partial pressures caused by gas phase reactions are also explored, showing good agreement with experimental results.
Collapse
Affiliation(s)
- Yang Hu
- Institute of Energy and Climate Research IEK-2, Forschungszentrum Jülich GmbH, 52428 Jülich, Germany
| | - Kai Wang
- Institute of Energy and Climate Research IEK-2, Forschungszentrum Jülich GmbH, 52428 Jülich, Germany
| | - Michael Müller
- Institute of Energy and Climate Research IEK-2, Forschungszentrum Jülich GmbH, 52428 Jülich, Germany
| | - Egbert Wessel
- Institute of Energy and Climate Research IEK-2, Forschungszentrum Jülich GmbH, 52428 Jülich, Germany
| | - Robert Spatschek
- Institute of Energy and Climate Research IEK-2, Forschungszentrum Jülich GmbH, 52428 Jülich, Germany
- JARA Energy, 52428 Jülich, Germany
| |
Collapse
|
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
|
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
|
10
|
Spicher S, Plett C, Pracht P, Hansen A, Grimme S. Automated Molecular Cluster Growing for Explicit Solvation by Efficient Force Field and Tight Binding Methods. J Chem Theory Comput 2022; 18:3174-3189. [PMID: 35482317 DOI: 10.1021/acs.jctc.2c00239] [Citation(s) in RCA: 35] [Impact Index Per Article: 17.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
An automated and broadly applicable workflow for the description of solvation effects in an explicit manner is introduced. This method, termed quantum cluster growth (QCG), is based on the semiempirical GFN2-xTB/GFN-FF methods, enabling efficient geometry optimizations and MD simulations. Fast structure generation is provided using the intermolecular force field xTB-IFF. Additionally, the approach uses an efficient implicit solvation model for the electrostatic embedding of the growing clusters. The novel QCG procedure presents a robust cluster generation tool for subsequent application of higher-level (e.g., DFT) methods to study solvation effects on molecular geometries explicitly or to average spectroscopic properties over cluster ensembles. Furthermore, the computation of the solvation free energy with a supermolecular approach can be carried out with QCG. The underlying growing process is physically motivated by computing the leading-order solute-solvent interactions first and can account for conformational and chemical changes due to solvation for low-energy barrier processes. The conformational space is explored with the NCI-MTD algorithm as implemented in the CREST program, using a combination of metadynamics and MD simulations. QCG with GFN2-xTB yields realistic solution geometries and reasonable solvation free energies for various systems without introducing many empirical parameters. Computed IR spectra of some solutes with QCG show a better match to the experimental data compared to well-established implicit solvation models.
Collapse
Affiliation(s)
- Sebastian Spicher
- Mulliken Center for Theoretical Chemistry, Institute of Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| | - Christoph Plett
- Mulliken Center for Theoretical Chemistry, Institute of Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| | - Philipp Pracht
- Mulliken Center for Theoretical Chemistry, Institute of Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Institute of Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Institute of Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| |
Collapse
|
11
|
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
|
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
|