1
|
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
|
2
|
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
|
3
|
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
|
4
|
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
|
5
|
Perlt E, Berger SA, Kelterer AM, Kirchner B. Anharmonicity of Vibrational Modes in Hydrogen Chloride-Water Mixtures. J Chem Theory Comput 2019; 15:2535-2547. [PMID: 30811198 DOI: 10.1021/acs.jctc.8b01070] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A thorough analysis of molecular vibrations in the binary system hydrogen chloride/water is presented considering a set of small mixed and pure clusters. In addition to the conventional normal-mode analysis based on the diagonalization of the Hessian, anharmonic frequencies were obtained from the perturbative VPT2 and PT2-VSCF method using hybrid density functional theory. For all normal modes, potential energy curves were modeled by displacing the atoms from the minimum geometry along the normal mode vectors. Three model potentials, a harmonic potential, a Morse potential, and a fourth order polynomial, were applied to fit these curves. From these data, it was possible not only to characterize distinct vibrations as mainly harmonic, anharmonic, or involving higher order terms but also to extract force constants, k, and anharmonicity constants, xe. By investigating all different types of intramolecular vibrations including covalent stretching or bending vibrations and intermolecular vibrations such as librations, we could demonstrate that while vibrational frequencies can be obtained applying scaling factors to harmonic results, useful anharmonicity constants cannot be predicted in such a way and the usage of more elaborate vibrational methods is necessary. For each particular type of molecular vibration, we could however determine a relationship between the wavenumber or wavenumber shift and the anharmonicity constant, which allows us to estimate mode dependent anharmonicity constants for larger clusters in the future.
Collapse
Affiliation(s)
- Eva Perlt
- Department of Chemistry , University of California, Irvine , 1102 Natural Sciences II , Irvine , California 92697-2025 , United States
| | - Sarah A Berger
- Institute of Physical and Theoretical Chemistry, NAWI Graz , Graz University of Technology , Stremayrgasse 9 , 8010 Graz , Austria
| | - Anne-Marie Kelterer
- Institute of Physical and Theoretical Chemistry, NAWI Graz , Graz University of Technology , Stremayrgasse 9 , 8010 Graz , Austria
| | - Barbara Kirchner
- Mulliken Center for Theoretical Chemistry , University of Bonn , Beringstrasse 4 , D-53115 Bonn , Germany
| |
Collapse
|
6
|
Abstract
One of the most important electrostatic interactions between molecules is most definitely the hydrogen bond. Understanding the basis of this interaction may offer us the insight needed to understand its effect on the macroscopic scale. Hydrogen bonding is for example the reason for anomalous properties in compounds like water and naturally life as we know it. The strength of the bond depends on numerous factors, among them the electronegativity of participating atoms. In this work we calculated the strength of hydrogen bonds between hydrides of the upper-right part of the periodic table (C, N, O, F, P, S, Cl, As, Se, Br) using quantum-chemical methods. The aim was to determine what influences the strength of strong and weak hydrogen bonds in simple hydrides. Various relationships were checked. A relation between the strength of the bond and the electronegativity of the participating atoms was found. We also observed a correlation between the strength of hydrogen bonds and the inter-atomic distances, along with the dependence on the charge transfer on the atom of the donor. We also report characteristic geometries of different dimers.
Collapse
|
7
|
von Domaros M, Perlt E. Anharmonic effects in the quantum cluster equilibrium method. J Chem Phys 2017; 146:124114. [PMID: 28388115 DOI: 10.1063/1.4978958] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The well-established quantum cluster equilibrium (QCE) model provides a statistical thermodynamic framework to apply high-level ab initio calculations of finite cluster structures to macroscopic liquid phases using the partition function. So far, the harmonic approximation has been applied throughout the calculations. In this article, we apply an important correction in the evaluation of the one-particle partition function and account for anharmonicity. Therefore, we implemented an analytical approximation to the Morse partition function and the derivatives of its logarithm with respect to temperature, which are required for the evaluation of thermodynamic quantities. This anharmonic QCE approach has been applied to liquid hydrogen chloride and cluster distributions, and the molar volume, the volumetric thermal expansion coefficient, and the isobaric heat capacity have been calculated. An improved description for all properties is observed if anharmonic effects are considered.
Collapse
Affiliation(s)
- Michael von Domaros
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstr. 6, D-53115 Bonn, Germany
| | - Eva Perlt
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstr. 6, D-53115 Bonn, Germany
| |
Collapse
|
8
|
von Domaros M, Jähnigen S, Friedrich J, Kirchner B. Quantum cluster equilibrium model of N-methylformamide–water binary mixtures. J Chem Phys 2016; 144:064305. [DOI: 10.1063/1.4941278] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Michael von Domaros
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstr. 4, D-53115 Bonn, Germany
| | - Sascha Jähnigen
- Martin-Luther-Universität Halle-Wittenberg, von-Danckelmann-Platz 4, D-06120 Halle, Germany
| | - Joachim Friedrich
- Technische Universität Chemnitz, Straße der Nationen 62, D-09111 Chemnitz, Germany
| | - Barbara Kirchner
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstr. 4, D-53115 Bonn, Germany
| |
Collapse
|
9
|
Matisz G, Kelterer AM, Fabian WMF, Kunsági-Máté S. Structural properties of methanol–water binary mixtures within the quantum cluster equilibrium model. Phys Chem Chem Phys 2015; 17:8467-79. [DOI: 10.1039/c4cp05836d] [Citation(s) in RCA: 48] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The Quantum Cluster Equilibrium (QCE) method computes cluster distributions and thermodynamic properties of binary methanol–water mixtures in agreement with experiments.
Collapse
Affiliation(s)
- G. Matisz
- Department of General and Physical Chemistry
- University of Pécs
- 7624 Pécs
- Hungary
- János Szentágothai Research Center
| | - A.-M. Kelterer
- Institute of Physical and Theoretical Chemistry
- Graz University of Technology
- NAWI Graz
- 8010 Graz
- Austria
| | - W. M. F. Fabian
- Institute of Chemistry
- University of Graz
- NAWI Graz
- 8010 Graz
- Austria
| | - S. Kunsági-Máté
- Department of General and Physical Chemistry
- University of Pécs
- 7624 Pécs
- Hungary
- János Szentágothai Research Center
| |
Collapse
|