1
|
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
|
2
|
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
|
3
|
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
|
4
|
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
|
5
|
Complex formation in methanol-chloroform solutions: Vibrational spectroscopy and quantum cluster equilibrium study. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.120499] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
6
|
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]
|
7
|
Taherivardanjani S, Blasius J, Brehm M, Dötzer R, Kirchner B. Conformer Weighting and Differently Sized Cluster Weighting for Nicotine and Its Phosphorus Derivatives. J Phys Chem A 2022; 126:7070-7083. [PMID: 36170053 DOI: 10.1021/acs.jpca.2c03133] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Weighting methods applied to systems with many conformers have been broadly employed to calculate thermodynamic properties, structural characteristics, and populations. To better understand and test the sensitivity of conventional weighting methods, the conformational distributions of nicotine and its phosphorus-substituted derivatives are investigated. The weighting schemes used for this are all based on Boltzmann statistics. Classical Boltzmann factors based on the electronic energy and the Gibbs free energy are calculated at different quantum chemical levels of theory and compared to cluster weights obtained by the quantum cluster equilibrium method. Furthermore, the influence of the modified rigid-rotor-harmonic-oscillator (mRRHO) approximation on the cluster weights is investigated. The substitution of the nitrogen atom in the methylpyrrolidine ring by a phosphorus atom results in more monomer conformers and clusters being populated. The conformational distribution of the monomers remained stable at different levels of theory and weighting methods. However, going to dimers and trimers, we observe a significant influence of the level of theory, weighting method, and mRRHO cutoff on the populations of these clusters. We show that mRRHO cutoff values of 50 and 100 cm-1 yield similar results, which is why 50 cm-1 is recommended as a robust choice. Furthermore, we observe that the global minimum for ΔE0 and ΔG varies in a few cases and that the global minimum is not always the dominantly occupied structure.
Collapse
Affiliation(s)
- Shima Taherivardanjani
- Mulliken Center for Theoretical Chemistry, Clausius Institute of Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4 + 6, D-53115 Bonn, Germany
| | - Jan Blasius
- Mulliken Center for Theoretical Chemistry, Clausius Institute of Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4 + 6, D-53115 Bonn, Germany
| | - Martin Brehm
- Institut für Chemie, Martin-Luther-Universität Halle-Wittenberg, D-06120 Halle (Saale), Germany
| | - Reinhard Dötzer
- Competence Center Analytics, BASF SE, D-67056 Ludwigshafen, Germany
| | - Barbara Kirchner
- Mulliken Center for Theoretical Chemistry, Clausius Institute of Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4 + 6, D-53115 Bonn, Germany
| |
Collapse
|
8
|
Weinhold F. Chlorine Dioxide: An Exception that Proves the Rules of Localized Chemical Bonding. J Chem Phys 2022; 156:124303. [DOI: 10.1063/5.0084739] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We employ natural bond orbital (NBO) and natural resonance theory (NRT) tools to analyze the enigmatic properties of the C2v-symmetric isomer of chlorine dioxide radical (ClO2), whose many challenges to Pauling-type localized bonding concepts were recognized by Linus Pauling himself. Although spin-contamination is minimal in this species, ClO2 exhibits an unusually strong form of "different Lewis structures for different spins" bonding pattern, intrinsically outside the framework of "maximal pairing" concepts. We show how the novel spin-unpaired donor-acceptor interactions lead to weakened bonding in the supramolecular domain of polyradical (ClO2)n homoclusters and aqueous ClO2(H2O)n heteroclusters. Despite feeble binding energies and large inter-radical separations, the polyradical clusters are found to maintain coherent spin patterns in each cluster component, attesting to the quantal donor-acceptor nature of their interactions and the cooperative and anticooperative couplings that govern intra- and intermolecular spin distributions in such spin-clusters.
Collapse
Affiliation(s)
- Frank Weinhold
- Department of Chemistry, University of Wisconsin-Madison, United States of America
| |
Collapse
|
9
|
|
10
|
Marchelli G, Ingenmey J, Hollóczki O, Chaumont A, Kirchner B. Hydrogen Bonding and Vaporization Thermodynamics in Hexafluoroisopropanol-Acetone and -Methanol Mixtures. A Joined Cluster Analysis and Molecular Dynamic Study. Chemphyschem 2021; 23:e202100620. [PMID: 34632686 PMCID: PMC9298724 DOI: 10.1002/cphc.202100620] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2021] [Revised: 09/28/2021] [Indexed: 12/23/2022]
Abstract
Binary mixtures of hexafluoroisopropanol with either methanol or acetone are analyzed via classical molecular dynamics simulations and quantum cluster equilibrium calculations. In particular, their populations and thermodynamic properties are investigated with the binary quantum cluster equilibrium method, using our in‐house code peacemaker 2.8, upgraded with temperature‐dependent parameters. A novel approach, where the final density from classical molecular dynamics, has been used to generate the necessary reference isobars. The hydrogen bond network in both type of mixtures at molar fraction of hexafluoroisopropanol of 0.2, 0.5, and 0.8 respectively is investigated via the molecular dynamics trajectories and the cluster results. In particular, the populations show that mixed clusters are preferred in both systems even at 0.2 molar fractions of hexafluoroisopropanol. Enthalpies and entropies of vaporization are calculated for the neat and mixed systems and found to be in good agreement with experimental values.
Collapse
Affiliation(s)
- Gwydyon Marchelli
- Mulliken Center for Theoretical Chemistry, Rheinische Friedrich-Wilhelms-Universität Bonn, Beringstr. 4+6, D-53115, Bonn, Germany
| | - Johannes Ingenmey
- Mulliken Center for Theoretical Chemistry, Rheinische Friedrich-Wilhelms-Universität Bonn, Beringstr. 4+6, D-53115, Bonn, Germany
| | - Oldamur Hollóczki
- Mulliken Center for Theoretical Chemistry, Rheinische Friedrich-Wilhelms-Universität Bonn, Beringstr. 4+6, D-53115, Bonn, Germany
| | - Alain Chaumont
- Université de Strasbourg, CNRS, CMC UMR 7140, Laboratoire MSM, F-67000, Strasbourg, France
| | - Barbara Kirchner
- Mulliken Center for Theoretical Chemistry, Rheinische Friedrich-Wilhelms-Universität Bonn, Beringstr. 4+6, D-53115, Bonn, Germany
| |
Collapse
|
11
|
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
|
12
|
Malloum A, Conradie J. Hydrogen bond networks of ammonia clusters: What we know and what we don’t know. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.116199] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
|
13
|
Teh S, Hsu PJ, Kuo JL. Size of the hydrogen bond network in liquid methanol: a quantum cluster equilibrium model with extensive structure search. Phys Chem Chem Phys 2021; 23:9166-9175. [PMID: 33885093 DOI: 10.1039/d1cp00427a] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Studies have debated what is a favorable cluster size in liquid methanol. Applications of the quantum cluster equilibrium (QCE) model on a limited set of cluster structures have demonstrated the dominance of cyclic hexamers in liquid methanol. In this study, we examined the aforementioned question by integrating our implementation of QCE with a molecular-dynamics-based structural searching scheme. QCE simulations were performed using a database comprising extensively searched stable conformers of (MeOH)n for n = 2-14, which were optimized by B3LYP/6-31+G(d,p) with and without the dispersion correction. Our analysis indicated that an octamer structure can contribute significantly to cluster probability. By reoptimizing selected conformers with high probability at the MP2 level, we found that the aforementioned octamer became the dominant species due to favorable vibrational free energy, which was attributed to modes of intermolecular vibration.
Collapse
Affiliation(s)
- Soon Teh
- Institute of Atomic and Molecular Sciences, Academia Sinica, Taipei 10617, Taiwan.
| | | | | |
Collapse
|
14
|
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
|
15
|
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
|