1
|
Mishra R, Bhawnani R, Sartape R, Chauhan R, Thorat AS, Singh MR, Shah JK. Role of Intermolecular Interactions in Deep Eutectic Solvents for CO 2 Capture: Vibrational Spectroscopy and Quantum Chemical Studies. J Phys Chem B 2024. [PMID: 39381893 DOI: 10.1021/acs.jpcb.4c04509] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/10/2024]
Abstract
Recent research and reviews on CO2 capture methods, along with advancements in industry, have highlighted high costs and energy-intensive nature as the primary limitations of conventional direct air capture and storage (DACS) methods. In response to these challenges, deep eutectic solvents (DESs) have emerged as promising absorbents due to their scalability, selectivity, and lower environmental impact compared to other absorbents. However, the molecular origins of their enhanced thermal stability and selectivity for DAC applications have not been explored before. Therefore, the current study focuses on a comprehensive investigation into the molecular interactions within an alkaline DES composed of potassium hydroxide (KOH) and ethylene glycol (EG). Combining Fourier transform infrared (FT-IR) and quantum chemical calculations, the study reports structural changes and intermolecular interactions induced in EG upon addition of KOH and its implications on CO2 capture. Experimental and computational spectroscopic studies confirm the presence of noncovalent interactions (hydrogen bonds) within both EG and the KOH-EG system and point to the aggregation of ions at higher KOH concentrations. Additionally, molecular electrostatic potential (MESP) surface analysis, natural bond orbital (NBO) analysis, quantum theory of atoms-in-molecules (QTAIM) analysis, and reduced density gradient-noncovalent interaction (RDG-NCI) plot analysis elucidate changes in polarizability, charge distribution, hydrogen bond types, noncovalent interactions, and interaction strengths, respectively. Evaluation of explicit and hybrid models assesses their effectiveness in representing intermolecular interactions. This research enhances our understanding of molecular interactions in the KOH-EG system, which are essential for both the absorption and desorption of CO2. The study also aids in predicting and selecting DES components, optimizing their ratios with salts, and fine-tuning the properties of similar solvents and salts for enhanced CO2 capture efficiency.
Collapse
Affiliation(s)
- Rashmi Mishra
- School of Chemical Engineering, Oklahoma State University, 420 Engineering North, Stillwater, Oklahoma 74078, United States
| | - Rajan Bhawnani
- Department of Chemical Engineering, University of Illinois at Chicago, 929 W. Taylor St., Chicago, Illinois 60607, United States
| | - Rohan Sartape
- Department of Chemical Engineering, University of Illinois at Chicago, 929 W. Taylor St., Chicago, Illinois 60607, United States
| | - Rohit Chauhan
- Department of Chemical Engineering, University of Illinois at Chicago, 929 W. Taylor St., Chicago, Illinois 60607, United States
| | - Amey S Thorat
- School of Chemical Engineering, Oklahoma State University, 420 Engineering North, Stillwater, Oklahoma 74078, United States
| | - Meenesh R Singh
- Department of Chemical Engineering, University of Illinois at Chicago, 929 W. Taylor St., Chicago, Illinois 60607, United States
| | - Jindal K Shah
- School of Chemical Engineering, Oklahoma State University, 420 Engineering North, Stillwater, Oklahoma 74078, United States
| |
Collapse
|
2
|
Ahmed U, Daub CD, Sundholm D, Johansson MP. Attractive acceptor-acceptor interactions in self-complementary quadruple hydrogen bonds for molecular self-assembly. Phys Chem Chem Phys 2024; 26:24470-24476. [PMID: 39264175 DOI: 10.1039/d4cp02361g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/13/2024]
Abstract
Molecular self-assembly provides the means for creating large supramolecular structures, extending beyond the capability of standard chemical synthesis. To harness the power of self-assembly, it is necessary to understand its driving forces. A potent method is to exploit self-complementary hydrogen bonding, where a molecule interacts with its own copy by suitable positions of hydrogen-bond donor (D) and acceptor (A) groups. With four hydrogen bonds, there are two possible self complementary patterns: the DDAA/AADD and the DADA/ADAD motifs. Of these, the DDAA pattern is usually more stable. The traditional explanation assumes that the secondary interactions between equal groups, that is, between donors (D⋯D) or acceptors (A⋯A), are repulsive. DDAA arrays would then have two, and DADA arrays six repulsive interactions. Here, using high-end quantum chemical analysis, we show that contrary to the traditional explanation, the secondary A⋯A interactions are, in fact, attractive. We revise the model of secondary interactions accordingly.
Collapse
Affiliation(s)
- Usman Ahmed
- Department of Chemistry, Faculty of Science, University of Helsinki, P. O. Box 55 (A. I. Virtasen aukio 1), FI-00014, Helsinki, Finland.
| | - Christopher D Daub
- Department of Chemistry, Faculty of Science, University of Helsinki, P. O. Box 55 (A. I. Virtasen aukio 1), FI-00014, Helsinki, Finland.
| | - Dage Sundholm
- Department of Chemistry, Faculty of Science, University of Helsinki, P. O. Box 55 (A. I. Virtasen aukio 1), FI-00014, Helsinki, Finland.
| | | |
Collapse
|
3
|
Wallace AM, Sherrill CD. Optimization of damping function parameters for -D3 and -D4 dispersion models for Hartree-Fock based symmetry-adapted perturbation theory. J Chem Phys 2024; 161:114115. [PMID: 39291687 DOI: 10.1063/5.0219185] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2024] [Accepted: 09/02/2024] [Indexed: 09/19/2024] Open
Abstract
Symmetry-adapted perturbation theory (SAPT) directly computes intermolecular interaction energy in terms of electrostatics, exchange-repulsion, induction/polarization, and London dispersion components. In SAPT based on Hartree-Fock ("SAPT0") or based on density functional theory, the most time-consuming step is the computation of the dispersion terms. Previous work has explored the replacement of these expensive dispersion terms with simple damped asymptotic models. We recently examined [Schriber et al. J. Chem. Phys. 154, 234107 (2021)] the accuracy of SAPT0 when replacing its dispersion term with Grimme's popular -D3 correction, reducing the computational cost scaling from O(N5) to O(N3). That work optimized damping function parameters for SAPT0-D3/jun-cc-pVDZ using estimates of the coupled-cluster complete basis set limit [CCSD(T)/CBS] on a 8299 dimer dataset. Here, we explore the accuracy of SAPT0-D3 with additional basis sets, along with an analogous model using -D4. Damping parameters are rather insensitive to basis sets, and the resulting SAPT0-D models are more accurate on average for total interaction energies than SAPT0. Our results are surprising in several respects: (1) improvement of -D4 over -D3 is negligible for these systems, even charged systems where -D4 should, in principle, be more accurate; (2) addition of Axilrod-Teller-Muto terms for three-body dispersion does not improve error statistics for this test set; and (3) SAPT0-D is even more accurate on average for total interaction energies than the much more computationally costly density functional theory based SAPT [SAPT(DFT)] in an aug-cc-pVDZ basis. However, SAPT0 and SAPT0-D3/D4 interaction energies benefit from significant error cancellation between exchange and dispersion terms.
Collapse
Affiliation(s)
- Austin M Wallace
- Center for Computational Molecular Science and Technology and School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| |
Collapse
|
4
|
Gao YJ, Chen WK, Fang WH, Cui G. A Green Host-Guest Protocol to Improve Water Solubility of Fluorescent Dyes. J Phys Chem Lett 2024; 15:9380-9387. [PMID: 39240654 DOI: 10.1021/acs.jpclett.4c02107] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/07/2024]
Abstract
Improving fluorescence emission efficiency is essential to develop novel luminescent materials. However, the low water solubility of conventional fluorescent dyes is a serious obstacle to broadening the application scope. Herein, a green protocol have been proposed: Two poorly water-soluble naphthalimide derivatives MONI and MANI with high fluorescent quantum yields (larger than 0.95 in toluene solution) were loaded in three different sizes of cyclodextrin (CD; α, β, γ-CD) with high water solubility. To further check the feasibility of the proposal, density functional theory (DFT) and time dependent-DFT (TD-DFT) methods combining the Own N-layer Integrated molecular Orbital molecular Mechanics (ONIOM) model with dispersion correction were employed to investigate the geometric and electronic structures of complexes CD·MXNI (X = N, O) in the excited-state process. TD-DFT calculations predict that the fantastic emission behavior of MXNI can be reserved after binding with CD, even improving fluorescent intensity in aqueous solution. Basis set superposition error (BSSE) correction and symmetry adapted perturbation theory (SAPT) were adopted to estimate the complexation energies and weak noncovalent interactions. The middle-sized β-CD is the perfect candidate to allow fluorescent molecules to settle into its cavity, forming an inclusion complex. Energy decomposition analysis (EDA) indicates that dispersion is superior to electrostatics interaction in embedding-type β-CD·MXNI, while it is contrary in α,γ-CD·MXNI. NMR calculations further prove the existence of a strong intermolecular hydrogen bond interaction between host and guest. Weak interactions that limited molecular vibration and hampered the nonradiative inactivation channel are conducive to the enhanced emission intensity.
Collapse
Affiliation(s)
- Yuan-Jun Gao
- State Key Laboratory of Organometallic Chemistry, Shanghai Institute of Organic Chemistry, Chinese Academy of Sciences, Shanghai 200032, China
| | - Wen-Kai Chen
- College of Chemistry and Materials Science, Hebei Key Laboratory of Inorganic Nano-Materials, Hebei Normal University, Shijiazhuang 050024, China
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, Chemistry College, Beijing Normal University, Beijing 100875, P. R. China
| | - Ganglong Cui
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, Chemistry College, Beijing Normal University, Beijing 100875, P. R. China
| |
Collapse
|
5
|
Barclay AJ, McKellar ARW, Pietropolli Charmet A, Moazzen-Ahmadi N. Water-carbon disulfide dimers: observation of a new isomer and ab initio structure theory. Phys Chem Chem Phys 2024; 26:23053-23061. [PMID: 39172124 DOI: 10.1039/d4cp02281e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/23/2024]
Abstract
The weakly bound dimer water-carbon disulfide is studied by ab initio structure theory and high-resolution infrared spectroscopy. The calculations yield three stable minima in the potential energy surface, all planar. The most stable, isomer 1, was observed previously by microwave spectroscopy. It has C2v symmetry, an S-O bond, and a linear O-S-C-S backbone. Isomer 2, the next most stable, has not been considered previously by theory or experiment. It also has C2v symmetry, but with a side-by-side structure. Isomer 3, which is slightly less stable, is side-by-side with one O-H bond pointing toward an S atom. The first gas phase water-CS2 infrared spectra are reported. For isomer 1, H2O-, HDO-, and D2O-CS2 are observed in the CS2ν1 + ν3 region, H2O-CS2 in the H2O ν2 bend region, and D2O-CS2 in the D2O ν1 and ν3 stretch regions. The latter ν3 spectrum enables an experimental determination of the A rotational parameter, which turns out to be larger than expected. The new isomer 2 is observed by means of a band in the H2O ν2 region, confirming its C2v symmetry and calculated structure.
Collapse
Affiliation(s)
- A J Barclay
- Department of Physics and Astronomy, University of Calgary, 2500 University Drive North West, Calgary, Alberta T2N 1N4, Canada.
| | - A R W McKellar
- National Research Council of Canada, Ottawa, Ontario K1A 0R6, Canada
| | - A Pietropolli Charmet
- Dipartimento di Scienze Molecolari e Nanosistemi, Università Ca' Foscari Venezia, Via Torino 155, I-30172, Mestre, Venezia, Italy
| | - N Moazzen-Ahmadi
- Department of Physics and Astronomy, University of Calgary, 2500 University Drive North West, Calgary, Alberta T2N 1N4, Canada.
| |
Collapse
|
6
|
Fishman V, Semidalas E, Martin JML. Basis Set Extrapolation from the Vanishing Counterpoise Correction Condition. J Phys Chem A 2024; 128:7462-7470. [PMID: 39167776 PMCID: PMC11382269 DOI: 10.1021/acs.jpca.4c03012] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/23/2024]
Abstract
Basis set extrapolations are typically rationalized either from analytical arguments involving the partial-wave or principal expansions of the correlation energy in helium-like systems or from fitting extrapolation parameters to reference energetics for a small(ish) training set. Seeking to avoid both, we explore a third alternative: extracting extrapolation parameters from the requirement that the BSSE (basis set superposition error) should vanish at the complete basis set limit. We find this to be a viable approach provided that the underlying basis sets are not too small and reasonably well balanced. For basis sets not augmented by diffuse functions, BSSE minimization and energy fitting yield quite similar parameters.
Collapse
Affiliation(s)
- Vladimir Fishman
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| | - Emmanouil Semidalas
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| | - Jan M L Martin
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
- On sabbatical at Quantum Theory Project, University of Florida, Gainesville, Florida 32611, United States
| |
Collapse
|
7
|
Xu Y, Biczysko M. Toward the identification of cyano-astroCOMs via vibrational features: benzonitrile as a test case. Front Chem 2024; 12:1439194. [PMID: 39296366 PMCID: PMC11408737 DOI: 10.3389/fchem.2024.1439194] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2024] [Accepted: 08/12/2024] [Indexed: 09/21/2024] Open
Abstract
The James Webb Space Telescope (JWST) opened a new era for the identification of molecular systems in the interstellar medium (ISM) by vibrational features. One group of molecules of increasing interest is cyano-derivatives of aromatic organic molecules, which have already been identified in the ISM on the basis of the analysis of rotational signatures, and so, are plausible candidates for the detection by the JWST. Benzonitrile considered in this work represents a suitable example for the validation of a computational strategy, which can be further applied for different, larger, and not-yet observed molecules. For this purpose, anharmonic simulations of infrared (IR) spectra have been compared with recent FTIR experimental studies. The anharmonic computations using the generalized second-order vibrational perturbation theory (GVPT2) in conjunction with a hybrid force field combining the harmonic part of revDSD-PBEP86-D3/jun-cc-pVTZ with anharmonic corrections from B3LYP-D3/SNSD show very good agreement with those in the experiment, with a mean error of 11 c m - 1 for all fundamental transitions overall and only 2 c m - 1 for the C ≡ N stretching fundamental at 4.49 μ m . The inclusion of overtones up to three-quanta transitions also allowed the prediction of spectra in the near-infrared region, which shows distinct features due to C ≡ N overtones at the 2.26 μ m and 1.52 μ m . The remarkable accuracy of the GVPT2 results opens a pathway for the reliable prediction of spectra for a broader range of cyano-astroCOMs.
Collapse
Affiliation(s)
- Yanting Xu
- International Centre for Quantum and Molecular Structures, Department of Physics, College of Sciences, Shanghai University, Shanghai, China
| | - Malgorzata Biczysko
- International Centre for Quantum and Molecular Structures, Department of Physics, College of Sciences, Shanghai University, Shanghai, China
| |
Collapse
|
8
|
Shahidi N, Laub JA, Vogiatzis KD, Doxastakis M. Multiscale Modeling of Vinyl-Addition Polynorbornenes: The Effect of Stereochemistry. Polymers (Basel) 2024; 16:2243. [PMID: 39204463 PMCID: PMC11358979 DOI: 10.3390/polym16162243] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2024] [Revised: 08/01/2024] [Accepted: 08/03/2024] [Indexed: 09/04/2024] Open
Abstract
Vinyl-addition polynorbornenes are candidates for designing high-performance polymers due to unique characteristics, which include a high glass transition temperature associated with a rigid backbone. Recent studies have established that the processability and properties of these polymers can be fine-tuned by using targeted substitutions. However, synthesis with different catalysts results in materials with distinct properties, potentially due to the presence of various stereoisomers that are difficult to quantify experimentally. Herein, we develop all-atom models of polynorbornene oligomers based on classical force fields and density functional theory. To establish the relationship between chemical architecture, chain conformations, and melt structure, we perform detailed molecular dynamics simulations with the fine-tuned atomistic force field and propose simpler coarse-grained descriptions to address the high molecular weight limit. All-atom simulations of oligomers suggest high glass transition temperatures in the range of 550-600 K. In the melt state (800 K), meso chains form highly rigid extended coils (C∞≈11) with amorphous structural characteristics similar to the X-ray diffraction data observed in the literature. In contrast, simulations with racemo chains predict highly helical tubular chain conformations that could promote assembly into crystalline structures.
Collapse
Affiliation(s)
- Nobahar Shahidi
- Department of Chemical and Biomolecular Engineering, University of Tennessee, Knoxville, TN 37996, USA
| | - Jeffrey A. Laub
- Department of Chemistry, University of Tennessee, Knoxville, TN 37996, USA
| | | | - Manolis Doxastakis
- Department of Chemical and Biomolecular Engineering, University of Tennessee, Knoxville, TN 37996, USA
| |
Collapse
|
9
|
Gray M, Herbert JM. Assessing the domain-based local pair natural orbital (DLPNO) approximation for non-covalent interactions in sizable supramolecular complexes. J Chem Phys 2024; 161:054114. [PMID: 39105555 PMCID: PMC11305816 DOI: 10.1063/5.0206533] [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/02/2024] [Accepted: 07/18/2024] [Indexed: 08/07/2024] Open
Abstract
The titular domain-based local pair natural orbital (DLPNO) approximation is the most widely used method for extending correlated wave function models to large molecular systems, yet its fidelity for intermolecular interaction energies in large supramolecular complexes has not been thoroughly vetted. Non-covalent interactions are sensitive to tails of the electron density and involve nonlocal dispersion that is discarded or approximated if the screening of pair natural orbitals (PNOs) is too aggressive. Meanwhile, the accuracy of the DLPNO approximation is known to deteriorate as molecular size increases. Here, we test the DLPNO approximation at the level of second-order Møller-Plesset perturbation theory (MP2) and coupled-cluster theory with singles, doubles, and perturbative triples [CCSD(T)] for a variety of large supramolecular complexes. DLPNO-MP2 interaction energies are within 3% of canonical values for small dimers with ≲10 heavy atoms, but for larger systems, the DLPNO approximation is often quite poor unless the results are extrapolated to the canonical limit where the threshold for discarding PNOs is taken to zero. Counterpoise correction proves to be essential in reducing errors with respect to canonical results. For a sequence of nanoscale graphene dimers up to (C96H24)2, extrapolated DLPNO-MP2 interaction energies agree with canonical values to within 1%, independent of system size, provided that the basis set does not contain diffuse functions; these cause the DLPNO approximation to behave erratically, such that results cannot be extrapolated in a meaningful way. DLPNO-CCSD(T) calculations are typically performed using looser PNO thresholds as compared to DLPNO-MP2, but this significantly impacts accuracy for large supramolecular complexes. Standard DLPNO-CCSD(T) settings afford errors of 2-6 kcal/mol for dimers involving coronene (C24H12) and circumcoronene (C54H18), even at the DLPNO-CCSD(T1) level.
Collapse
Affiliation(s)
- Montgomery Gray
- Department of Chemistry & Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - John M. Herbert
- Department of Chemistry & Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|
10
|
Xu R, Jiang Z, Yang Q, Bloino J, Biczysko M. Harmonic and anharmonic vibrational computations for biomolecular building blocks: Benchmarking DFT and basis sets by theoretical and experimental IR spectrum of glycine conformers. J Comput Chem 2024; 45:1846-1869. [PMID: 38682874 DOI: 10.1002/jcc.27377] [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: 02/12/2024] [Revised: 04/01/2024] [Accepted: 04/02/2024] [Indexed: 05/01/2024]
Abstract
Advanced vibrational spectroscopic experiments have reached a level of sophistication that can only be matched by numerical simulations in order to provide an unequivocal analysis, a crucial step to understand the structure-function relationship of biomolecules. While density functional theory (DFT) has become the standard method when targeting medium-size or larger systems, the problem of its reliability and accuracy are well-known and have been abundantly documented. To establish a reliable computational protocol, especially when accuracy is critical, a tailored benchmark is usually required. This is generally done over a short list of known candidates, with the basis set often fixed a priori. In this work, we present a systematic study of the performance of DFT-based hybrid and double-hybrid functionals in the prediction of vibrational energies and infrared intensities at the harmonic level and beyond, considering anharmonic effects through vibrational perturbation theory at the second order. The study is performed for the six-lowest energy glycine conformers, utilizing available "state-of-the-art" accurate theoretical and experimental data as reference. Focusing on the most intense fundamental vibrations in the mid-infrared range of glycine conformers, the role of the basis sets is also investigated considering the balance between computational cost and accuracy. Targeting larger systems, a broad range of hybrid schemes with different computational costs is also tested.
Collapse
Affiliation(s)
- Ruiqin Xu
- Department of Physics, College of Sciences, Shanghai University, Shanghai, China
| | | | - Qin Yang
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Science, Prague, Czechia
| | - Julien Bloino
- Classe di Scienze, Scuola Normale Superiore, Pisa, Italy
| | - Malgorzata Biczysko
- Department of Physics, College of Sciences, Shanghai University, Shanghai, China
| |
Collapse
|
11
|
Scaglione N, Avila J, Padua A, Costa Gomes M. Tailored carbon dioxide capacity in carboxylate-based ionic liquids. Faraday Discuss 2024. [PMID: 39099430 DOI: 10.1039/d4fd00052h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/06/2024]
Abstract
We have used a library of thermally stable tetraalkylphosphonium carboxylate ionic liquids that were easily prepared from available carboxylic acids. Depending on the pKa in water of the precursor acids, the resulting ionic liquids either dissolve or reversibly chemically absorb CO2, with some exhibiting notable gas capacities, reaching a CO2 mole fraction of 0.2 at 1 bar and 343 K. While equilibrium constants and ionic liquid capacities generally correlate with the pKa of the acids, certain exceptions underscore the influence of liquid structure and physical properties of the ionic liquids, elucidated through molecular dynamics simulations and density functional theory calculations. Unlike the trends observed in other CO2-absorbing ILs, phosphonium carboxylates do not experience increased viscosity upon gas absorption; instead, enhanced diffusivities are observed, facilitating efficient gas-liquid transfer.
Collapse
Affiliation(s)
- Nicolas Scaglione
- Laboratoire de Chimie de l'ENS Lyon, CNRS, Université de Lyon, 46 allée d'Italie, 69364 Lyon, France.
| | - Jocasta Avila
- Laboratoire de Chimie de l'ENS Lyon, CNRS, Université de Lyon, 46 allée d'Italie, 69364 Lyon, France.
| | - Agilio Padua
- Laboratoire de Chimie de l'ENS Lyon, CNRS, Université de Lyon, 46 allée d'Italie, 69364 Lyon, France.
| | - Margarida Costa Gomes
- Laboratoire de Chimie de l'ENS Lyon, CNRS, Université de Lyon, 46 allée d'Italie, 69364 Lyon, France.
| |
Collapse
|
12
|
Shang Y, Yan G, Cai Y, Lu L, Zhao H, Sun R. Theoretical Investigation on Water-Free, Water- and Self-Assisted H-Abstraction Reactions from Dimethylamine by Hydroxy Radicals. J Phys Chem A 2024; 128:6264-6273. [PMID: 39034617 DOI: 10.1021/acs.jpca.4c02732] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/23/2024]
Abstract
Accurate branching ratios of the H-abstraction reactions from dimethylamine (DMA) by OH radicals are important in understanding the atmospheric fate of DMA. In this work, the reaction kinetics of the water-free, water-assisted, and self-assisted H-abstraction reactions between DMA and OH radicals are accurately determined using the multipath canonical variational theory with the small-curvature tunneling correction, to explore the catalytic effects of the reactant (DMA) and product (water). To choose a suitable method that well describes the current reaction systems, various combinations with seven DFT methods and six basis sets are first evaluated, and the M08-HX/ma-TZVP method is identified as the most appropriate, with a mean unsigned deviation of 0.9 kcal mol-1 against the gold-standard CCSD(T)/CBS(T-Q) method. Based on the determined potential energy surfaces with the considerations of ground-state structures and specific-reaction parameters of zero-point energies, rate constants and branching ratios are calculated in a wide temperature range. The calculations show that the participation of water and DMA can lead to three-body complexes with a lower energy and influence the energy barriers, but neither of them shows the catalytic effect on the H-abstraction reactions in terms of kinetics. Additionally, the branching ratio analysis demonstrates that the product distribution is significantly altered in the presence of DMA and water.
Collapse
Affiliation(s)
- Yanlei Shang
- Energy Research Institute, Qilu University of Technology (Shandong Academy of Sciences), Jinan, Shandong 250014, P. R. China
- Shandong Technology Innovation Center of Carbon Neutrality, Jinan, Shandong 250014, P. R. China
- Key Laboratory of Extreme Material Dynamics Technology, Chengdu, Sichuan 610031, P. R. China
| | - Guihuan Yan
- Shandong Technology Innovation Center of Carbon Neutrality, Jinan, Shandong 250014, P. R. China
- Ecology Institute of Shandong Academy of Sciences, Jinan, Shandong 250014, P. R. China
| | - Yang Cai
- Key Laboratory of Extreme Material Dynamics Technology, Chengdu, Sichuan 610031, P. R. China
- The Peac Institute of Multiscale Sciences, Chengdu, Sichuan 610031, P. R. China
| | - Lei Lu
- Key Laboratory of Extreme Material Dynamics Technology, Chengdu, Sichuan 610031, P. R. China
- School of Materials Science and Engineering, Dynamic Materials Data Science Center, Southwest Jiaotong University, Chengdu, Sichuan 610031, P. R. China
| | - Haiyong Zhao
- Xiling DigitIntel Institute, Chengdu, Sichuan 610000, P. R. China
| | - Rongfeng Sun
- Energy Research Institute, Qilu University of Technology (Shandong Academy of Sciences), Jinan, Shandong 250014, P. R. China
- Shandong Technology Innovation Center of Carbon Neutrality, Jinan, Shandong 250014, P. R. China
| |
Collapse
|
13
|
Rais N, Salta Z, Tasinato N. Theoretical investigation of the OH-initiated atmospheric degradation mechanism of CX 2CHX (X = H, F, Cl) by advanced quantum chemical and transition state theory methods. Phys Chem Chem Phys 2024; 26:19976-19991. [PMID: 38995148 DOI: 10.1039/d4cp01453g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/13/2024]
Abstract
Halogenated olefins are anthropogenic compounds with many industrial applications but at the same time raising many environmental and health concerns. Gas-phase electrophilic addition of the OH radical to the olefinic CC bond represents the primary sink for these chemicals in the atmosphere, with the degree and type of halogenation playing a significant role in their overall reactivity. In this work, we present a theoretical investigation of the reaction mechanisms and kinetics for the reactions between the OH radical and CH2CH2 (ethylene, ETH), CF2CHF (trifluoroethylene, TFE) and CCl2CHCl (trichloroethylene, TCE), simulated by state-of-the-art protocols and methods, with the aim of providing a detailed interpretation of the available experimental results, as well as new data of relevance to tropospheric chemistry. Specifically, potential energy surfaces (PESs) are obtained using the jun-Cheap (jChS) composite scheme, whereas temperature and pressure dependent rate coefficients and product distributions in the 100-600 K temperature range are calculated within the Rice-Ramsperger-Kassel-Marcus/master equation (RRKM/ME) framework. The rates for barrierless channels are obtained from variable reaction coordinate-variational transition state theory (VRC-VTST) combined with the two transition state model. While the reactions with ETH and TFE proceed mainly via the formation of addition adducts at P = 1 atm and T = 298 K, the dominant channel for TCE is the Cl-elimination reaction. Global rate constants for the two halogenated olefins, TFE and TCE, are found to be pressure-independent, contrary to the case of ETH. The computed rate constants, as well as their temperature and pressure dependence, are in remarkable agreement with the available experimental data, and they are used to derive atmospheric lifetimes (τ) for both TFE and TCE as a function of altitude (h) in the atmosphere, by taking into account variations in the rate coefficients (k (T, P)) and [OH] concentration.
Collapse
Affiliation(s)
- Nadjib Rais
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126, Pisa, Italy.
- IUSS Scuola Universitaria Superiore, Piazza della Vittoria 15, I-27100, Pavia, Italy
| | - Zoi Salta
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126, Pisa, Italy.
| | - Nicola Tasinato
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126, Pisa, Italy.
| |
Collapse
|
14
|
Falbo E, Lavecchia A. HessFit: A Toolkit to Derive Automated Force Fields from Quantum Mechanical Information. J Chem Inf Model 2024; 64:5634-5645. [PMID: 38897917 DOI: 10.1021/acs.jcim.4c00540] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/21/2024]
Abstract
In this study, we introduce a novel approach to enhance the accuracy of molecular dynamics simulations by refining the force fields (FFs) through a combination of transferable parameters and molecule-specific characteristics derived from quantum mechanical (QM) calculations. Traditional FFs often prioritize generality over precision, leading to limitations in the accuracy of accurately capturing intra- and intermolecular interactions. To address this, we present an open-source toolkit, called HessFit, designed to integrate QM-derived bonded parameters and atomic charges into existing FFs. In combination with bond, angle, torsional, and nonbonded parameters derivation, HessFit can easily extract multiple barrier terms of dihedrals from QM Hessian and gradient or return all terms through a fitting procedure scheme of QM potential energy surface. We showcase the effectiveness of HessFit through comprehensive evaluations of vibrational properties across a diverse set of small molecules and demonstrate that experimental results support its ability in predicting thermodynamic properties of organic molecules compared to previous state-of-the-art approaches. We further explore its application to Zn2+ metalloprotein models, which are hard systems to treat with automatic approaches. Our results demonstrate that HessFit parameters compete with GAFF2 and OPLS parameters to describing small organic molecules, and its feasibility is also comparable to current FFs used to modeling nonstandard residues in Zn proteins for molecular dynamics simulations. The effectiveness of the HessFit protocol makes it a valuable tool for deriving or extending force field parameters for novel compounds in several molecular modeling applications.
Collapse
Affiliation(s)
- Emanuele Falbo
- Department of Pharmacy, Drug Discovery Lab, University of Naples Federico II, Naples 80131, Italy
| | - Antonio Lavecchia
- Department of Pharmacy, Drug Discovery Lab, University of Naples Federico II, Naples 80131, Italy
| |
Collapse
|
15
|
Shang Y, Luo SN. Insights into the role of the H-abstraction reaction kinetics of amines in understanding their degeneration fates under atmospheric and combustion conditions. Phys Chem Chem Phys 2024. [PMID: 39028293 DOI: 10.1039/d4cp02187h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/20/2024]
Abstract
Amines, a class of prototypical volatile organic compounds, have garnered considerable interest within the context of atmospheric and combustion chemistry due to their substantial contributions to the formation of hazardous pollutants in the atmosphere. In the current energy landscape, the implementation of carbon-neutral energy and strategic initiatives leads to generation of new amine sources that cannot be overlooked in terms of the emission scale. To reduce the emission level of amines from their sources and mitigate their impact on the formation of harmful substances, a comprehensive understanding of the fundamental reaction kinetics during the degeneration process of amines is imperative. This perspective article first presents an overview of both traditional amine sources and emerging amine sources within the context of carbon peaking and carbon neutrality and then highlights the importance of H-abstraction reactions in understanding the atmospheric and combustion chemistry of amines from the perspective of reaction kinetics. Subsequently, the current experimental and theoretical techniques for investigating the H-abstraction reactions of amines are introduced, and a concise summary of research endeavors made in this field over the past few decades is provided. In order to provide accurate kinetic parameters of the H-abstraction reactions of amines, advanced kinetic calculations are performed using the multi-path canonical variational theory combined with the small-curvature tunneling and specific-reaction parameter methods. By comparing with the literature data, current kinetic calculations are comprehensively evaluated, and these validated data are valuable for the development of the reaction mechanism of amines.
Collapse
Affiliation(s)
- Yanlei Shang
- Energy Research Institute, Qilu University of Technology (Shandong Academy of Sciences), Jinan, Shandong, 250014, P. R. China.
- School of Materials Science and Engineering, Southwest Jiaotong University, Chengdu, Sichuan, 610031, P. R. China
- Key Laboratory of Extreme Material Dynamics Technology, Chengdu, Sichuan 610031, P. R. China
| | - S N Luo
- School of Materials Science and Engineering, Southwest Jiaotong University, Chengdu, Sichuan, 610031, P. R. China
- Key Laboratory of Extreme Material Dynamics Technology, Chengdu, Sichuan 610031, P. R. China
| |
Collapse
|
16
|
Bovolenta GM, Silva-Vera G, Bovino S, Molpeceres G, Kästner J, Vogt-Geisse S. In-depth exploration of catalytic sites on amorphous solid water: I. The astrosynthesis of aminomethanol. Phys Chem Chem Phys 2024; 26:18692-18706. [PMID: 38922674 DOI: 10.1039/d4cp01865f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/27/2024]
Abstract
Chemical processes taking place on ice-grain mantles are pivotal to the complex chemistry of interstellar environments. In this study, we conducted a comprehensive analysis of the catalytic effects of an amorphous solid water (ASW) surface on the reaction between ammonia (NH3) and formaldehyde (H2CO) to form aminomethanol (NH2CH2OH) using density functional theory. We identified potential catalytic sites based on the binding energy distribution of NH3 and H2CO reactants, on a set-of-clusters surface model composed of 22 water molecules and found a total of 14 reaction paths. Our results indicate that the catalytic sites can be categorized into four groups, depending on the interactions of the carbonyl oxygen and the amino group with the ice surface in the reactant complex. A detailed analysis of the reaction mechanism using Intrinsic Reaction Coordinate and reaction force analysis, revealed three distinct chemical events for this reaction: formation of the C-N bond, breaking of the N-H bond, and formation of the O-H hydroxyl bond. Depending on the type of catalytic site, these events can occur within a single, concerted, albeit asynchronous, step, or can be isolated in a step-wise mechanism, with the lowest overall transition state energy observed at 1.3 kcal mol-1. A key requirement for the low-energy mechanism is the presence of a pair of dangling OH bonds on the surface, found at 5% of the potential catalytic sites on an ASW porous surface.
Collapse
Affiliation(s)
- Giulia M Bovolenta
- Departamento de Físico-Química, Facultad de Ciencias Químicas, Universidad de Concepción, Concepción, Chile.
- Atomistic Simulations, Italian Institute of Technology, 16152 Genova, Italy
| | - Gabriela Silva-Vera
- Departamento de Físico-Química, Facultad de Ciencias Químicas, Universidad de Concepción, Concepción, Chile.
| | - Stefano Bovino
- Chemistry Department, Sapienza University of Rome, P.le A. Moro, 00185 Rome, Italy
- INAF - Osservatorio Astrofisico di Arcetri, Largo E. Fermi 5, 50125 Firenze, Italy
- Departamento de Astronomía, Facultad de Ciencias Físicas y Matemáticas, Universidad de Concepción, Av. Esteban Iturra s/n Barrio Universitario, Concepción, Chile
| | - German Molpeceres
- Departamento de Astrofísica Molecular Instituto de Física Fundamental (IFF-CSIC), Madrid, Spain
| | - Johannes Kästner
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Stefan Vogt-Geisse
- Departamento de Físico-Química, Facultad de Ciencias Químicas, Universidad de Concepción, Concepción, Chile.
| |
Collapse
|
17
|
Salta Z, Schaefer T, Tasinato N, Kieninger M, Katz A, Herrmann H, Ventura ON. Energetics of the OH radical H-abstraction reactions from simple aldehydes and their geminal diol forms. J Mol Model 2024; 30:253. [PMID: 38970670 DOI: 10.1007/s00894-024-06058-0] [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/02/2024] [Accepted: 07/01/2024] [Indexed: 07/08/2024]
Abstract
CONTEXT Carbonyl compounds, especially aldehydes, emitted to the atmosphere, may suffer hydration in aerosols or water droplets in clouds. At the same time, they can react with hydroxyl radicals which may add or abstract hydrogen atoms from these species. The interplay between hydration and hydrogen abstraction is studied using density functional and quantum composite theoretical methods, both in the gas phase and in simulated bulk water. The H-abstraction from the aldehydic and geminal diol forms of formaldehyde, acetaldehyde, glycolaldehyde, glyoxal, methylglyoxal, and acrolein is studied to determine whether the substituent has any noticeable effect in the preference for the abstraction of one form or another. It is found that abstraction of the H-atom adjacent to the carbonyl group gives a more stable radical than same abstraction from the geminal diol in the case of formaldehyde, acetaldehyde, and glycolaldehyde. The presence of a delocalizing group in the Cα (a carbonyl group in glyoxal and methylglyoxal, and a vinyl group in acrolein), reverts this trend, and now the abstraction of the H-atom from the geminal diol gives more stable radicals. A further study was conducted abstracting hydrogen atoms from the other different positions in the species considered, both in the aldehydic and geminal diol forms. Only in the case of glycolaldehyde, the radical formed by H-abstraction from the -CH2OH group is more stable than any of the other radical species. Abstraction of the hydrogen atom in one of the hydroxyl groups in the geminal diol is equivalent to the addition of the •OH radical to the aldehyde. It leads, in some cases, to decomposition into a smaller radical and a neutral molecule. In these cases, some interesting theoretical differences are observed between the results in gas phase and (simulated) bulk solvent, as well as with respect to the method of calculation chosen. METHODS DFT (M06-2X, B2PLYP, PW6B95), CCSD(T), and composite (CBS-QB3, jun-ChS, SCVECV-f12) methods using Dunning basis sets and extrapolation to the CBS limit were used to study the energetics of closed shell aldehydes in their keto and geminal-diol forms, as well as the radical derived from them by hydrogen abstraction. Both gas phase and simulated bulk solvent calculations were performed, in the last case using the Polarizable Continuum Model.
Collapse
Affiliation(s)
- Zoi Salta
- Scuola Normale Superiore, Piazza Dei Cavalieri 7, 56126, Pisa, Italy
| | - Thomas Schaefer
- Atmospheric Chemistry Department (ACD), Leibniz Institute for Tropospheric Research (TROPOS), Permoserstr. 15, 04318, Leipzig, Germany
| | - Nicola Tasinato
- Scuola Normale Superiore, Piazza Dei Cavalieri 7, 56126, Pisa, Italy
| | - Martina Kieninger
- Computational Chemistry and Biology Group, Facultad de Química, CCBG, Universidad de La República, 11400, Montevideo, DETEMA, Uruguay
| | - Aline Katz
- Computational Chemistry and Biology Group, Facultad de Química, CCBG, Universidad de La República, 11400, Montevideo, DETEMA, Uruguay
| | - Hartmut Herrmann
- Atmospheric Chemistry Department (ACD), Leibniz Institute for Tropospheric Research (TROPOS), Permoserstr. 15, 04318, Leipzig, Germany
| | - Oscar N Ventura
- Computational Chemistry and Biology Group, Facultad de Química, CCBG, Universidad de La República, 11400, Montevideo, DETEMA, Uruguay.
| |
Collapse
|
18
|
Salta Z, Ventura ON, Rais N, Tasinato N, Barone V. A new chapter in the never ending story of cycloadditions: The puzzling case of SO 2 and acetylene. J Comput Chem 2024; 45:1587-1602. [PMID: 38517313 DOI: 10.1002/jcc.27350] [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: 01/30/2024] [Revised: 03/04/2024] [Accepted: 03/08/2024] [Indexed: 03/23/2024]
Abstract
A comprehensive study of the different classes of cycloaddition reactions ([3+2], [2+2], and [2+1]) of SO2 to acetylene and ethylene has been performed using density functional theory (DFT) and composite wavefunction methods. The [3+2] cycloaddition reaction, that was previously explored in the context of the cycloaddition of thioformaldehyde S-methylide (TSM) to ethylene and acetylene, proceeds in a concerted way to the formation of stable heterocycles. In this paper, we extend our study to the [2+2] and [2+1] cycloadditions of SO2 to acetylene, which would produce 1,1-oxathiete-2-oxide and thiirene-1,1-dioxide, respectively. One of the main conclusions is that cyclic 1,1-oxathiete-2-oxide can open through a relatively easy breaking of the SO single bond and rearrange toward sulfinyl acetaldehyde (SA). The SA molecule can easily undergo several internal rearrangements, which eventually lead to sulfenic acid and sulfoxide derivatives of ethenone, 1,2,3-dioxathiole, and CO plus sulfinylmethane. The most probable path, however, produces 2-thioxoacetic acid, whose derivatives (or those of the corresponding acetate) are usually obtained by Willgerodt-Kindler-type sulfuration of acetates. This product can in turn decompose, leading to the final products CO2 and H2CS. Comparison of this decomposition path with that of 2-amino-2-thioxoacetic acid shows that the process occurs through different H-transfer processes.
Collapse
Affiliation(s)
- Zoi Salta
- Scuola Normale Superiore, Pisa, Italy
| | - Oscar N Ventura
- Computational Chemistry and Biology Group, CCBG, DETEMA, Facultad de Química, Universidad de la República, Montevideo, Uruguay
| | - Nadjib Rais
- Scuola Normale Superiore, Pisa, Italy
- IUSS Scuola Universitaria Superiore, Pavia, Italy
| | | | | |
Collapse
|
19
|
Martín-Fernández C, Elguero J, Alkorta I. Beryllium as a Base: Complexes of Be(CO) 3 with HX (X=F, Cl, Br, CN, NC, CCH, OH). Chemphyschem 2024:e202400608. [PMID: 38950128 DOI: 10.1002/cphc.202400608] [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: 05/31/2024] [Revised: 06/28/2024] [Accepted: 07/01/2024] [Indexed: 07/03/2024]
Abstract
Beryllium chemistry is typically governed by its electron deficient character, but in some compounds it can act as a base. In order to understand better the unusual basicity of Be, we have systematically explored the complexes of one such compound, Be(CO)3, towards several hydrogen bond donors HX (X=F, Cl, Br, CN, NC, CCH, OH). For all complexes we find three different minima, two hydrogen bonded minima (to the Be or O atoms), and one weak beryllium bonded minimum. Further characterization of the interactions using a topological analysis of the electron density and Symmetry Adapted Perturbation Theory (SAPT) provide insight into the nature of these interactions. Overall these results highlight the capability of certain beryllium compounds to act as either a weak Lewis acid or, unconventionally, a Lewis base whose basicity towards hydrogen bonding is comparable to that of π systems.
Collapse
Affiliation(s)
| | - José Elguero
- Instituto de Química Médica (CSIC), Juan de la Cierva, 3, 28006, Madrid, Spain
| | - Ibon Alkorta
- Instituto de Química Médica (CSIC), Juan de la Cierva, 3, 28006, Madrid, Spain
| |
Collapse
|
20
|
Huang Y, Shen M, Yan H, He Y, Xu J, Zhu F, Yang X, Ye YX, Ouyang G. Achieving a solar-to-chemical efficiency of 3.6% in ambient conditions by inhibiting interlayer charges transport. Nat Commun 2024; 15:5406. [PMID: 38926358 PMCID: PMC11208529 DOI: 10.1038/s41467-024-49373-z] [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: 01/03/2024] [Accepted: 06/04/2024] [Indexed: 06/28/2024] Open
Abstract
Efficiently converting solar energy into chemical energy remains a formidable challenge in artificial photosynthetic systems. To date, rarely has an artificial photosynthetic system operating in the open air surpassed the highest solar-to-biomass conversion efficiency (1%) observed in plants. In this study, we present a three-dimension polymeric photocatalyst achieving a solar-to-H2O2 conversion efficiency of 3.6% under ambient conditions, including real water, open air, and room temperature. The impressive performance is attributed to the efficient storage of electrons inside materials via expeditious intramolecular charge transfer, and the fast extraction of the stored electrons by O2 that can diffuse into the internal pores of the self-supporting three-dimensional material. This construction strategy suppresses the interlayer transfer of excitons, polarizers and carriers, effectively increases the utilization of internal excitons to 82%. This breakthrough provides a perspective to substantially enhance photocatalytic performance and bear substantial implications for sustainable energy generation and environmental remediation.
Collapse
Affiliation(s)
- Yuyan Huang
- Key Laboratory of Bioinorganic and Synthetic Chemistry of Ministry of Education, LIFM, School of Chemistry, IGCME, Sun Yat-sen University, Guangzhou, 510275, China
| | - Minhui Shen
- Key Laboratory of Bioinorganic and Synthetic Chemistry of Ministry of Education, LIFM, School of Chemistry, IGCME, Sun Yat-sen University, Guangzhou, 510275, China
| | - Huijie Yan
- School of Chemical Engineering and Technology, IGCME, Sun Yat-sen University, Zhuhai, 519082, China
| | - Yingge He
- School of Environmental Science and Engineering, Guangdong Provincial Key Laboratory of Environmental Pollution Control and Remediation Technology, Sun Yat-sen University, Guangzhou, 510275, China
| | - Jianqiao Xu
- Key Laboratory of Bioinorganic and Synthetic Chemistry of Ministry of Education, LIFM, School of Chemistry, IGCME, Sun Yat-sen University, Guangzhou, 510275, China
| | - Fang Zhu
- Key Laboratory of Bioinorganic and Synthetic Chemistry of Ministry of Education, LIFM, School of Chemistry, IGCME, Sun Yat-sen University, Guangzhou, 510275, China
| | - Xin Yang
- School of Environmental Science and Engineering, Guangdong Provincial Key Laboratory of Environmental Pollution Control and Remediation Technology, Sun Yat-sen University, Guangzhou, 510275, China
| | - Yu-Xin Ye
- School of Chemical Engineering and Technology, IGCME, Sun Yat-sen University, Zhuhai, 519082, China.
- Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai), Zhuhai, Guangdong, 519082, China.
| | - Gangfeng Ouyang
- Key Laboratory of Bioinorganic and Synthetic Chemistry of Ministry of Education, LIFM, School of Chemistry, IGCME, Sun Yat-sen University, Guangzhou, 510275, China.
- School of Chemical Engineering and Technology, IGCME, Sun Yat-sen University, Zhuhai, 519082, China.
- Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai), Zhuhai, Guangdong, 519082, China.
| |
Collapse
|
21
|
Zhuang L, Wang J, Wan J, Huang C. Why do dipole moments of HCl-water clusters fail to determine acid dissociation? Phys Chem Chem Phys 2024; 26:17910-17917. [PMID: 38888219 DOI: 10.1039/d4cp01316f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/20/2024]
Abstract
This paper quantitatively examines why dipole moments of HCl(H2O)n=1-8 cannot serve as the dissociation criterion for acid molecules using the Hirshfeld-I approach. Also, we propose the possible experimental parameter 〈P(HCl)〉, whose statistical average enables the assessment of acid dissociation in mixed clusters. Furthermore, our calculations reveal that a minimum of four water molecules are necessary to dissociate an HCl molecule.
Collapse
Affiliation(s)
- Lei Zhuang
- School of Materials Science and Physics, China University of Mining and Technology, Xuzhou, Jiangsu 221116, China.
| | - Jing Wang
- School of Materials Science and Physics, China University of Mining and Technology, Xuzhou, Jiangsu 221116, China.
| | - Jianguo Wan
- School of Physics, Nanjing University, Nanjing, Jiangsu 210093, China
| | - Chuanfu Huang
- School of Materials Science and Physics, China University of Mining and Technology, Xuzhou, Jiangsu 221116, China.
| |
Collapse
|
22
|
Yang M, Wang J. Comprehensive Multipath Variational Kinetics Study on Hydrogen Abstraction Reactions from Three Typical Dimethylcyclohexane Isomers by Hydroxyl Radicals: from the Electronic Structure to Model Applications. J Phys Chem A 2024; 128:4517-4531. [PMID: 38804972 DOI: 10.1021/acs.jpca.4c00480] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/29/2024]
Abstract
Cycloalkanes serve as an important class of chemical components in both fossil and alternative transportation fuels and have attracted considerable attention from the combustion community. Hydrogen abstractions from cycloalkanes by hydroxyl radicals initiate the fuel decomposition process and trigger off the subsequent chain reactions and thus play an important role in both combustion and atmospheric chemistry. The target of this study is to fill the vacancy in kinetics data toward the H-abstraction reactions by hydroxyl radical from three typical dimethylcyclohexane isomers through first-principles and direct dynamics. The rate constants involving 18 elementary reactions in total were accurately determined by the multipath canonical variational transition state theory with the multidimensional small-curvature correction for tunneling (MP-CVT/SCT), over a broad temperature range of 200-2000 K. The significant roles of multistructural torsional anharmonicity and recrossing effects were stressed per abstraction site, while the quantum tunneling effect was found to be slight at temperatures of interest in combustion. The discrepancies observed among different reaction systems at a similar abstraction site highlight the fuel molecular effects on site-specific rate constants. The comparison results of total rate constants given by different dynamics approaches prove the importance of considering the torsional anharmonicity, recrossing, and tunneling effects, and the robust feature of the simplified MS-CVT/SCT. The calculated total constants for dimethylcyclohexane isomers by OH are consistent with those measured for methylcyclohexane and 1,4-dimethylcyclohexane at low temperatures. The branching ratio analysis confirms the predominant role of the tertiary abstraction at low-to-intermediate temperatures and its growing competition with distinct secondary abstractions as temperature increases. The calculated rate constants were eventually fitted into the analytical expressions and incorporated into the kinetic models to learn about the influences on modeling performance.
Collapse
Affiliation(s)
- Mo Yang
- National Key Laboratory of Science and Technology on Aero-Engine Aero-Thermodynamics, School of Energy and Power Engineering, Beihang University, Beijing 100191, PR China
| | - Juan Wang
- National Key Laboratory of Science and Technology on Aero-Engine Aero-Thermodynamics, School of Energy and Power Engineering, Beihang University, Beijing 100191, PR China
| |
Collapse
|
23
|
Torres-Boy A, Taccone MI, Kirschbaum C, Ober K, Stein T, Meijer G, von Helden G. Investigation of the Proton-Bound Dimer of Dihydrogen Phosphate and Formate Using Infrared Spectroscopy in Helium Droplets. J Phys Chem A 2024; 128:4456-4466. [PMID: 38771224 PMCID: PMC11163467 DOI: 10.1021/acs.jpca.4c01632] [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/12/2024] [Revised: 05/03/2024] [Accepted: 05/06/2024] [Indexed: 05/22/2024]
Abstract
Understanding the structural and dynamic properties of proton-bound complexes is crucial for elucidating fundamental aspects of chemical reactivity and molecular interactions. In this work, the proton-bound complex between dihydrogen phosphate and formate, and its deuterated counterparts, is investigated using IR action spectroscopy in helium droplets. Contrary to the initial expectation that the stronger phosphoric acid would donate a proton to formate, both experiment and theory show that all exchangeable protons are located in the phosphate moiety. The experimental spectra show good agreement with both scaled harmonic and VPT2 anharmonic calculations, indicating that anharmonic effects are small. Some H-bending modes of the nondeuterated complex are found to be sensitive to the helium environment. In the case of the partially deuterated complexes, the experiments indicate that internal dynamics leads to isomeric interconversion upon IR excitation.
Collapse
Affiliation(s)
| | - Martín I. Taccone
- Fritz
Haber Institute of the Max Planck Society, 14195 Berlin, Germany
| | - Carla Kirschbaum
- Fritz
Haber Institute of the Max Planck Society, 14195 Berlin, Germany
- Institute
of Chemistry and Biochemistry, Freie Universität
Berlin, 14195 Berlin, Germany
| | - Katja Ober
- Fritz
Haber Institute of the Max Planck Society, 14195 Berlin, Germany
| | - Tamar Stein
- Institute
of Chemistry and Fritz Haber Center for Molecular Dynamics, Hebrew University of Jerusalem, 91904 Jerusalem, Israel
| | - Gerard Meijer
- Fritz
Haber Institute of the Max Planck Society, 14195 Berlin, Germany
| | - Gert von Helden
- Fritz
Haber Institute of the Max Planck Society, 14195 Berlin, Germany
| |
Collapse
|
24
|
Liu Q, Huang Y, Ye YX. Construction of Conjugated Organic Polymers for Efficient Photocatalytic Hydrogen Peroxide Generation with Adequate Utilization of Water Oxidation. MATERIALS (BASEL, SWITZERLAND) 2024; 17:2709. [PMID: 38893973 PMCID: PMC11173575 DOI: 10.3390/ma17112709] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/22/2024] [Revised: 05/15/2024] [Accepted: 05/21/2024] [Indexed: 06/21/2024]
Abstract
The visible-light-driven photocatalytic production of hydrogen peroxide (H2O2) is currently an emerging approach for transforming solar energy into chemical energy. In general, the photocatalytic process for producing H2O2 includes two pathways: the water oxidation reaction (WOR) and the oxygen reduction reaction (ORR). However, the utilization efficiency of ORR surpasses that of WOR, leading to a discrepancy with the low oxygen levels in natural water and thereby impeding their practical application. Herein, we report a novel donor-bridge-acceptor (D-B-A) organic polymer conjugated by the Sonogashira-Hagihara coupling reaction with tetraphenylethene (TPE) units as the electron donors, acetylene (A) as the connectors and pyrene (P) moieties as the electron acceptors. Notably, the resulting TPE-A-P exhibits a remarkable solar-to-chemical conversion of 1.65% and a high BET-specific surface area (1132 m2·g-1). Furthermore, even under anaerobic conditions, it demonstrates an impressive H2O2 photosynthetic efficiency of 1770 μmol g-1 h-1, exceeding the vast majority of previously reported photosynthetic systems of H2O2. The outstanding performance is attributed to the effective separation of electrons and holes, along with the presence of sufficient reaction sites facilitated by the incorporation of alkynyl electronic bridges. This protocol presents a successful method for generating H2O2 via a water oxidation reaction, signifying a significant advancement towards practical applications in the natural environment.
Collapse
Affiliation(s)
- Qinzhe Liu
- Key Laboratory of Bioinorganic and Synthetic Chemistry of Ministry of Education, LIFM—Lehn Institute of Functional Materials, School of Chemistry, IGCME—Institute of Green Chemistry and Molecular Engineering, Sun Yat-sen University, Guangzhou 510275, China (Y.H.)
| | - Yuyan Huang
- Key Laboratory of Bioinorganic and Synthetic Chemistry of Ministry of Education, LIFM—Lehn Institute of Functional Materials, School of Chemistry, IGCME—Institute of Green Chemistry and Molecular Engineering, Sun Yat-sen University, Guangzhou 510275, China (Y.H.)
| | - Yu-xin Ye
- School of Chemical Engineering and Technology, IGCME, Sun Yat-sen University, Zhuhai 519082, China
- Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai), Zhuhai 519082, China
| |
Collapse
|
25
|
Mehta N, Martin JML. On the sensitivity of computed partial charges toward basis set and (exchange-)correlation treatment. J Comput Chem 2024; 45:1017-1032. [PMID: 38216516 DOI: 10.1002/jcc.27294] [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: 09/25/2023] [Revised: 12/05/2023] [Accepted: 12/10/2023] [Indexed: 01/14/2024]
Abstract
Partial charges are a central concept in general chemistry and chemical biology, yet dozens of different computational definitions exist. In prior work [Cho et al., ChemPhysChem 21, 688-696 (2020)], we showed that these can be reduced to at most three 'principal components of ionicity'. The present study addressed the dependence of computed partial charges q on 1-particle basis set and (for WFT methods) n -particle correlation treatment or (for DFT methods) exchange-correlation functional, for several representative partial charge definitions such as QTAIM, Hirshfeld, Hirshfeld-I, HLY (electrostatic), NPA, and GAPT. Our findings show that semi-empirical double hybrids can closely approach the CCSD(T) 'gold standard' for this property. In fact, owing to an error compensation in MP2, CCSD partial charges are further away from CCSD(T) than is MP2. The nonlocal correlation is important, especially when there is a substantial amount of nonlocal exchange. Employing range separation proves to be "mostly" not advantageous, while global hybrids perform optimally for 20%-30% Hartree-Fock exchange across all charge types. Basis set convergence analysis shows that an augmented triple-zeta heavy-aug-cc-pV(T+d)Z basis set or a partially augmented jun-cc-pV(T+d)Z basis set is sufficient for Hirshfeld, Hirshfeld-I, HLY, and GAPT charges. In contrast, QTAIM and NPA display slower basis set convergence. It is noteworthy that for both NPA and QTAIM, HF exhibits markedly slower basis set convergence than the correlation components of MP2 and CCSD. Triples corrections in CCSD(T), denoted as CCSD(T)-CCSD, exhibit even faster basis set convergence.
Collapse
Affiliation(s)
- Nisha Mehta
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Reḥovot, Israel
| | - Jan M L Martin
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Reḥovot, Israel
| |
Collapse
|
26
|
Chen W, Zhang D, Fu H, Li J, Yu X, Zhou J, Lu B. Restructuring Electrolyte Solvation by a Partially and Weakly Solvating Cosolvent toward High-Performance Potassium-Ion Batteries. ACS NANO 2024; 18:12512-12523. [PMID: 38701404 DOI: 10.1021/acsnano.4c02108] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/05/2024]
Abstract
Ether-based electrolytes are among the most important electrolytes for potassium-ion batteries (PIBs) due to their low polarization voltage and notable compatibility with potassium metal. However, their development is hindered by the strong binding between K+ and ether solvents, leading to [K+-solvent] cointercalation on graphite anodes. Herein, we propose a partially and weakly solvating electrolyte (PWSE) wherein the local solvation environment of the conventional 1,2-dimethoxyethane (DME)-based electrolyte is efficiently reconfigured by a partially and weakly solvating diethoxy methane (DEM) cosolvent. For the PWSE in particular, DEM partially participates in the solvation shell and weakens the chelation between K+ and DME, facilitating desolvation and suppressing cointercalation behavior. Notably, the solvation structure of the DME-based electrolyte is transformed into a more cation-anion-cluster-dominated structure, consequently promoting thin and stable solid-electrolyte interphase (SEI) generation. Benefiting from optimized solvation and SEI generation, the PWSE enables a graphite electrode with reversible K+ (de)intercalation (for over 1000 cycles) and K with reversible plating/stripping (the K||Cu cell with an average Coulombic efficiency of 98.72% over 400 cycles) and dendrite-free properties (the K||K cell operates over 1800 h). We demonstrate that rational PWSE design provides an approach to tailoring electrolytes toward stable PIBs.
Collapse
Affiliation(s)
- Weijie Chen
- School of Physics and Electronics, Hunan University, Changsha 410082, P. R. China
| | - Dianwei Zhang
- School of Physics and Electronics, Hunan University, Changsha 410082, P. R. China
| | - Hongwei Fu
- School of Physics and Electronics, Hunan University, Changsha 410082, P. R. China
| | - Jinfan Li
- School of Physics and Electronics, Hunan University, Changsha 410082, P. R. China
| | - Xinzhi Yu
- School of Physics and Electronics, Hunan University, Changsha 410082, P. R. China
- Greater Bay Area Institute for Innovation, Hunan University, Guangzhou, Guangdong Province 511300, China
| | - Jiang Zhou
- School of Materials Science and Engineering, Central South University, Changsha 410082, P. R. China
| | - Bingan Lu
- School of Physics and Electronics, Hunan University, Changsha 410082, P. R. China
| |
Collapse
|
27
|
Molpeceres G, Tsuge M, Furuya K, Watanabe N, San Andrés D, Rivilla VM, Colzi L, Aikawa Y. Carbon Atom Condensation on NH 3-H 2O Ices. An Alternative Pathway to Interstellar Methanimine and Methylamine. J Phys Chem A 2024. [PMID: 38709949 DOI: 10.1021/acs.jpca.3c08286] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/08/2024]
Abstract
The recent discovery of the nature and behavior of carbon atoms interacting with interstellar ices has prompted a number of investigations on the chemistry initiated by carbon accretion on icy interstellar dust. In this work, we expand the range of processes promoted by carbon accretion to the chemistry initiated by the interaction of this atom with ammonia (NH3) using quantum chemical calculations. We found that carbon addition to the ammonia molecule forms a rather stable radical, CNH3, that is easily hydrogenated. The complete hydrogenation network is later studied. Our calculations reveal that while conversion to simpler molecules like HCN and HNC is indeed a possible outcome promoted by H-abstraction reactions, methylamine is also easily formed (CH3NH2). In fact, the stability of methylamine against hydrogen abstraction makes this molecule the preferred product of the reaction network. Our results serve as a stepping stone toward the accurate modeling of C-addition reactions in realistic astrochemical kinetic models.
Collapse
Affiliation(s)
- Germán Molpeceres
- Department of Astronomy, Graduate School of Science, The University of Tokyo, Tokyo 113 0033, Japan
- Departamento de Astrofísica Molecular, Instituto de Física Fundamental (IFF-CSIC), C/Serrano 121, 28006 Madrid, Spain
| | - Masashi Tsuge
- Institute of Low Temperature Science, Hokkaido University, Sapporo, Hokkaido 060-0819, Japan
| | - Kenji Furuya
- Department of Astronomy, Graduate School of Science, The University of Tokyo, Tokyo 113 0033, Japan
- National Astronomical Observatory of Japan, Mitaka, Tokyo 181-8588, Japan
| | - Naoki Watanabe
- Institute of Low Temperature Science, Hokkaido University, Sapporo, Hokkaido 060-0819, Japan
| | - David San Andrés
- Centro de Astrobiología (CAB), INTA-CSIC, Carretera de Ajalvir km 4, Torrejón de Ardoz, 28850 Madrid, Spain
| | - Víctor M Rivilla
- Centro de Astrobiología (CAB), INTA-CSIC, Carretera de Ajalvir km 4, Torrejón de Ardoz, 28850 Madrid, Spain
| | - Laura Colzi
- Centro de Astrobiología (CAB), INTA-CSIC, Carretera de Ajalvir km 4, Torrejón de Ardoz, 28850 Madrid, Spain
| | - Yuri Aikawa
- Department of Astronomy, Graduate School of Science, The University of Tokyo, Tokyo 113 0033, Japan
| |
Collapse
|
28
|
Cao W, Warneke J, Wang XB. Probing the Electronic Structure of [B 10H 10] 2- Dianion Encapsulated by an Octamethylcalix[4]pyrrole Molecule. J Phys Chem A 2024; 128:3361-3369. [PMID: 38651632 DOI: 10.1021/acs.jpca.4c01736] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/25/2024]
Abstract
Despite being an important closo-borate in condensed phase boron chemistry, isolated [B10H10]2- is electronically unstable and has never been detected in the gas phase. Herein, we report a successful capture of this fleeting species through binding with an octamethylcalix[4]pyrrole (omC4P) molecule to form a stable gaseous omC4P·[B10H10]2- complex and its characterizations utilizing negative ion photoelectron spectroscopy (NIPES). The recorded NIPE spectrum, contributed by both omC4P and [B10H10]2-, is deconvoluted by subtracting the omC4P contribution to yield a [B10H10]2- spectrum. The obtained [B10H10]2- spectrum consists of four major bands spanning the electron binding energy (EBE) range from 1 to 5 eV, with the EBE gaps matching excellently with the energy intervals of computed high-lying occupied molecular orbitals of the [B10H10]2- dianion. This study showcases a generic method to utilize omC4P to capture unstable multiply charged anions in the gas phase for experimental determination of their electronic structures.
Collapse
Affiliation(s)
- Wenjin Cao
- Physical Sciences Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, Richland, Washington 99352, United States
| | - Jonas Warneke
- Wilhelm-Ostwald-Institute for Physical and Theoretical Chemistry, Leipzig University, Leipzig 04103, Germany
- Leibniz Institute of Surface Engineering (IOM), Permoserstraße 15, Leipzig 04318, Germany
| | - Xue-Bin Wang
- Physical Sciences Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, Richland, Washington 99352, United States
| |
Collapse
|
29
|
Hou T, Yuan X, Jiang S, Xu Z, Zhang X, Lu M, Xu Y. Experimental detection of the diamino-pentazolium cation and theoretical exploration of derived high energy materials. Sci Rep 2024; 14:10120. [PMID: 38698073 PMCID: PMC11065884 DOI: 10.1038/s41598-024-60741-z] [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: 12/29/2023] [Accepted: 04/26/2024] [Indexed: 05/05/2024] Open
Abstract
In this work, we realized the detection of diamino-pentazolium cation (DAPZ+) in the reaction solution experimentally and proved it to be meta-diamino-pentazole based on the transition state theory. Quantum chemical methods were used to predict its spectral properties, charge distribution, stability and aromaticity. Considering that DAPZ+ has excellent detonation properties, it was further explored by assembling it with N5-, N3- and C(NO2)3- anions, respectively. The results show a strong interaction between DAPZ+ and the three anions, which will have a positive effect on its stability. Thanks to the high enthalpy of formation and density, the calculated detonation properties of the three systems are exciting, especially [DAPZ+][N5-] (D: 10,016 m·s-1; P: 37.94 GPa), whose actual detonation velocity may very likely exceed CL-20 (D: 9773 m·s-1). There is no doubt that this work will become the precursor for the theoretical exploration of new polynitrogen ionic compounds.
Collapse
Affiliation(s)
- Tianyang Hou
- School of Chemistry and Chemical Engineering, Nanjing University of Science and Technology, Nanjing, 210094, China
| | - Xiaofeng Yuan
- School of Chemistry and Chemical Engineering, Nanjing University of Science and Technology, Nanjing, 210094, China
| | - Shuaijie Jiang
- School of Chemistry and Chemical Engineering, Nanjing University of Science and Technology, Nanjing, 210094, China
| | - Ze Xu
- School of Chemistry and Chemical Engineering, Nanjing University of Science and Technology, Nanjing, 210094, China
| | - Xiaopeng Zhang
- School of Chemistry and Chemical Engineering, Nanjing University of Science and Technology, Nanjing, 210094, China
| | - Ming Lu
- School of Chemistry and Chemical Engineering, Nanjing University of Science and Technology, Nanjing, 210094, China.
| | - Yuangang Xu
- School of Chemistry and Chemical Engineering, Nanjing University of Science and Technology, Nanjing, 210094, China.
| |
Collapse
|
30
|
Tang Y, Yang Y, Cao Q, Wang K, Ma D, Liu B, Zhu Q, Ma J. Theoretical Study of the Temperature- and Pressure-Dependent Rate Constants for Nine Reactions between CO n ( n = 0-4), O m ( m = 1-3), C 2O, and C 3O 2 during the Radiolysis of Carbon Dioxide. J Phys Chem A 2024. [PMID: 38691833 DOI: 10.1021/acs.jpca.4c00681] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/03/2024]
Abstract
We investigate the reaction pathways of nine important CO2-related reactions using the revDSD-PBEP86-D3(BJ)/jun-cc-pV(T+d)Z level and simultaneously employ an accurate composite method (jun-Cheap) based on coupled-cluster (CC) theory. Subsequently, the Rice-Ramsperger-Kassel-Marcus/master equation (RRKM/ME) is solved to calculate the temperature- and pressure-dependent rate constants. This work investigates reactions involving transition states that have been overlooked in previous literature, including the dissociation of singlet-state C3O2, the triple channel formation of C2O + CO to form C3O2, and the formation of O3 + CO. The results show that CO3 is highly prone to dissociation at high temperatures. Finally, the kinetic data show that over a wide temperature range, our calculations are consistent with previous experimental measurements. The majority of the reaction rate constants studied exhibit significant pressure dependence, while the O3 + CO reaction is pressure-independent at low temperatures. These results are instrumental in the development of detailed kinetic models for the CO2 radiolysis reaction network.
Collapse
Affiliation(s)
- Yuhong Tang
- Institute of Atomic and Molecular Physics, Sichuan University, Chengdu 610065, China
| | - Yu Yang
- Reactor Operation and Application Research Sub-Institute, Nuclear Power Institute of China, Chengdu 610041, China
| | - Qi Cao
- Reactor Operation and Application Research Sub-Institute, Nuclear Power Institute of China, Chengdu 610041, China
| | - Keli Wang
- Institute of Atomic and Molecular Physics, Sichuan University, Chengdu 610065, China
| | - Dandan Ma
- Institute of Atomic and Molecular Physics, Sichuan University, Chengdu 610065, China
| | - Beibei Liu
- Institute of Atomic and Molecular Physics, Sichuan University, Chengdu 610065, China
| | - Quan Zhu
- School of Chemical Engineering, Sichuan University, Chengdu 610065, PR China
- Engineering Research Center of Combustion and Cooling for Aerospace Power, Ministry of Education, Sichuan University, Chengdu 610065, PR China
| | - Jianyi Ma
- Institute of Atomic and Molecular Physics, Sichuan University, Chengdu 610065, China
- Engineering Research Center of Combustion and Cooling for Aerospace Power, Ministry of Education, Sichuan University, Chengdu 610065, PR China
| |
Collapse
|
31
|
Xiong X, Friedman R, Wu W, Su P. QM/MM-Based Energy Decomposition Analysis Method for Large Systems. J Phys Chem A 2024. [PMID: 38687960 DOI: 10.1021/acs.jpca.4c00183] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/02/2024]
Abstract
In this work, a QM/MM-based EDA method, called GKS-EDA(QM/MM), is proposed. As an extension of GKS-EDA, this scheme divides the total interaction energy into electrostatic, exchange-repulsion, polarization, and correlation/dispersion terms. GKS-EDA(QM/MM) can be applied to describe the interactions of large-scale systems combined with various QM/MM platforms. By using the examples of a hydrated hydronium ion complex in water solution, the barnase-barstar complex, and MMP-13-pyrimidinetrione in a metalloprotein, the capability of GKS-EDA(QM/MM) for various interactions in large systems is validated.
Collapse
Affiliation(s)
- Xuewei Xiong
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
- Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen 361005, China
| | - Ran Friedman
- Department of Chemistry and Biomedical Sciences, Linnaeus University, 39182 Kalmar, Sweden
| | - Wei Wu
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
- Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen 361005, China
| | - Peifeng Su
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
- Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen 361005, China
| |
Collapse
|
32
|
Mendolicchio M, Barone V. Unbiased Comparison between Theoretical and Experimental Molecular Structures and Properties: Toward an Accurate Reduced-Cost Evaluation of Vibrational Contributions. J Chem Theory Comput 2024; 20:2842-2857. [PMID: 38556752 DOI: 10.1021/acs.jctc.4c00023] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/02/2024]
Abstract
The tremendous development of hardware and software is constantly increasing the role of quantum chemical (QC) computations in the assignment and interpretation of experimental results. However, an unbiased comparison between theory and experiment requires the proper account of vibrational averaging effects. In particular, high-resolution spectra in the gas phase are now available for molecules containing up to about 50 atoms, which are too large for a brute-force approach with the available QC methods of sufficient accuracy. In the present paper, we introduce hybrid approaches, which allow the accurate evaluation of vibrational averaging effects for molecules of this size beyond the harmonic approximation, with special attention being devoted to rotational constants. After the validation of new tools for relatively small molecules, the β-estradiol hormone and a prototypical molecular motor have been considered to witness the feasibility of accurate computations for large molecules.
Collapse
|
33
|
Cao W, Wang XB. Organic Molecules Mimic Alkali Metals Enabling Spontaneous Harpoon Reactions with Halogens. Chemistry 2024; 30:e202400038. [PMID: 38287792 DOI: 10.1002/chem.202400038] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2024] [Revised: 01/29/2024] [Accepted: 01/30/2024] [Indexed: 01/31/2024]
Abstract
The harpoon mechanism has been a milestone in molecular reaction dynamics. Until now, the entity from which electron harpooning occurs has been either alkali metal atoms or non-metallic analogs in their excited states. In this work, we demonstrate that a common organic molecule, octamethylcalix[4] pyrrole (omC4P), behaves just like alkali metal atoms, enabling the formation of charge-separated ionic bonding complexes with halogens omC4P+ ⋅ X- (X=F-I, SCN) via the harpoon mechanism. Their electronic structures and chemical bonding were determined by cryogenic photoelectron spectroscopy of the corresponding anions and confirmed by theoretical analyses. The omC4P+ ⋅ X- could be visualized to form from the reactants omC4P+X via electron harpooning from omC4P to X at a distance defined by the energy difference between the ionization potential of omC4P and electron affinity of X.
Collapse
Affiliation(s)
- Wenjin Cao
- Physical Sciences Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P. O. Box 999, MS J7-10, Richland, WA, 99352, USA
| | - Xue-Bin Wang
- Physical Sciences Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P. O. Box 999, MS J7-10, Richland, WA, 99352, USA
| |
Collapse
|
34
|
Zhang J, Gulyás KV, Li J, Ma M, Zhou L, Wu L, Xiong R, Erdelyi M, Zhu W, Xu Z. Unexpected effect of halogenation on the water solubility of small organic compounds. Comput Biol Med 2024; 172:108209. [PMID: 38460313 DOI: 10.1016/j.compbiomed.2024.108209] [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: 01/15/2024] [Revised: 02/06/2024] [Accepted: 02/25/2024] [Indexed: 03/11/2024]
Abstract
Halogenation is an indispensable method in the structural modification of lead compounds. It is known to increase lipophilicity and is hence used to improve membrane permeability and thus bioavailability. In this study, we compare the water solubility (logS) of organohalogen compounds and their non-halogenated parent compounds using the molecular matched pair (MMP) analysis method. Unexpectedly, 19.9% of the compounds increased their water solubility upon halogenation. Iodination was observed to have the greatest effect on solubility, followed by chlorination, bromination, and fluorination. Introducing amino, hydroxyl and carboxyl groups into organohalogens improves their aqueous solubilities, whereas introducing a trifluoromethyl group has the opposite effect. According to our quantum chemical calculations, the increased water solubility upon halogenation is, at least partially, attributed to an increased polarity and polarizability. These results improve our understanding of the influence of halogenation on bioactivity.
Collapse
Affiliation(s)
- Jianfang Zhang
- School of Chinese Materia Medica, Nanjing University of Chinese Medicine, Nanjing, 210046, China; State Key Laboratory of Drug Research; Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, 201203, China
| | - Kinga Virág Gulyás
- Department of Chemistry - BMC, Uppsala University, SE-751 23, Uppsala, Sweden
| | - Jintian Li
- State Key Laboratory of Drug Research; Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, 201203, China
| | - Minfei Ma
- State Key Laboratory of Drug Research; Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, 201203, China
| | - Liping Zhou
- State Key Laboratory of Drug Research; Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, 201203, China
| | - Leyun Wu
- State Key Laboratory of Drug Research; Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, 201203, China
| | - Ruisheng Xiong
- Department of Chemistry - BMC, Uppsala University, SE-751 23, Uppsala, Sweden
| | - Mate Erdelyi
- Department of Chemistry - BMC, Uppsala University, SE-751 23, Uppsala, Sweden
| | - Weiliang Zhu
- School of Chinese Materia Medica, Nanjing University of Chinese Medicine, Nanjing, 210046, China; State Key Laboratory of Drug Research; Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, 201203, China.
| | - Zhijian Xu
- School of Chinese Materia Medica, Nanjing University of Chinese Medicine, Nanjing, 210046, China; State Key Laboratory of Drug Research; Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, 201203, China.
| |
Collapse
|
35
|
Deng CL, Obi AD, Tra BYE, Sarkar SK, Dickie DA, Gilliard RJ. Air- and photo-stable luminescent carbodicarbene-azaboraacenium ions. Nat Chem 2024; 16:437-445. [PMID: 38052948 DOI: 10.1038/s41557-023-01381-0] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2023] [Accepted: 10/23/2023] [Indexed: 12/07/2023]
Abstract
Substitution of a C=C bond by an isoelectronic B-N bond is a well-established strategy to alter the electronic structure and stability of acenes. BN-substituted acenes that possess narrow energy gaps have attractive optoelectronic properties. However, they are susceptible to air and/or light. Here we present the design, synthesis and molecular structures of fully π-conjugated cationic BN-doped acenes stabilized by carbodicarbene ligands. They are luminescent in the solution and solid states and show high air and moisture stability. Compared with their neutral BN-substituted counterparts as well as the parent all-carbon acenes, these species display improved quantum yields and small optical gaps. The electronic structures of the azabora-anthracene and azabora-tetracene cations resemble higher-order acenes while possessing high photo-oxidative resistance. Investigations using density functional theory suggest that the stability and photo-physics of these conjugated systems may be ascribed to their cationic nature and the electronic properties of the carbodicarbene.
Collapse
Affiliation(s)
- Chun-Lin Deng
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Akachukwu D Obi
- Department of Chemistry, University of Virginia, Charlottesville, VA, USA
| | - Bi Youan E Tra
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Samir Kumar Sarkar
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Diane A Dickie
- Department of Chemistry, University of Virginia, Charlottesville, VA, USA
| | - Robert J Gilliard
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA, USA.
| |
Collapse
|
36
|
Lazzari F, Mendolicchio M, Barone V. Accurate Geometries of Large Molecules by Integration of the Pisa Composite Scheme and the Templating Synthon Approach. J Phys Chem A 2024; 128:1385-1395. [PMID: 38347709 DOI: 10.1021/acs.jpca.3c08382] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/23/2024]
Abstract
An effective yet reliable computational workflow is proposed, which permits the computation of accurate geometrical structures for large flexible molecules at an affordable cost thanks to the integration of machine learning tools and DFT models together with reduced scaling computations of vibrational averaging effects. After validation of the different components of the overall strategy, a panel of molecules of biological interest have been analyzed. The results confirm that very accurate geometrical parameters can be obtained at reasonable cost for molecules including up to about 50 atoms, which are the largest ones for which comparison with high-resolution rotational spectra is possible. Since the whole computational workflow can be followed employing standard electronic structure codes, accurate results for large-sized molecules can be obtained at DFT cost also by nonspecialists.
Collapse
Affiliation(s)
- Federico Lazzari
- Scuola Normale Superiore di Pisa, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Marco Mendolicchio
- Scuola Normale Superiore di Pisa, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | | |
Collapse
|
37
|
Yu X, Liu H, Ling S, Wu X, Lian C, Xu J. Microfluidic Printing of Vertically-Oriented Nanosheets/MOFs Hetero-Interface for Intensive Pseudocapacitive Storage. SMALL (WEINHEIM AN DER BERGSTRASSE, GERMANY) 2024; 20:e2305396. [PMID: 37797184 DOI: 10.1002/smll.202305396] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/28/2023] [Revised: 09/16/2023] [Indexed: 10/07/2023]
Abstract
Efficient manufacture of electroactive vertically-oriented nanosheets with enhanced electrolyte mass diffusion and strong interfacial redox dynamics is critical for realizing high energy density of miniature supercapacitor (SC), but still challenging. Herein, microfluidic droplet printing is developed to controllably construct vertically-oriented graphene/ZIF-67 hetero-microsphere (VAGS/ZIF-67), where the ZIF-67 is coordinately grown on vertically-oriented graphene framework via Co─O─C bonds. The VAGS/ZIF-67 shows ordered porous channel, high electroactivity and strong interfacial interaction, providing rapid electrolyte diffusion dynamics and high faradaic capacitance in KOH solution (1674 F g-1 , 1004 C g-1 ), which are verified by computational fluid dynamics (CFD) and density functional theory (DFT). Moreover, the VAGS/ZIF-67 based SC exhibits large energy density (100 Wh kg-1 ), excellent durability (10 000 cycles and high/low temperature), and robust power-supply applications in portable electronics.
Collapse
Affiliation(s)
- Xude Yu
- The State Key Laboratory of Chemical Engineering, Department of Chemical Engineering, Tsinghua University, Beijing, 100084, P. R. China
| | - Hengyuan Liu
- The State Key Laboratory of Chemical Engineering, Department of Chemical Engineering, Tsinghua University, Beijing, 100084, P. R. China
| | - Sida Ling
- The State Key Laboratory of Chemical Engineering, Department of Chemical Engineering, Tsinghua University, Beijing, 100084, P. R. China
| | - Xingjiang Wu
- The State Key Laboratory of Chemical Engineering, Department of Chemical Engineering, Tsinghua University, Beijing, 100084, P. R. China
| | - Cheng Lian
- The State Key Laboratory of Chemical Engineering and Shanghai Engineering Research Center of Hierarchical Nanomaterials, School of Chemical Engineering, School of Chemistry and Molecular Engineering, East China University of Science and Technology, Shanghai, 200237, China
| | - Jianhong Xu
- The State Key Laboratory of Chemical Engineering, Department of Chemical Engineering, Tsinghua University, Beijing, 100084, P. R. China
| |
Collapse
|
38
|
Luo T, Wang Y, Elander B, Goldstein M, Mu Y, Wilkes J, Fahrenbruch M, Lee J, Li T, Bao JL, Mohanty U, Wang D. Polysulfides in Magnesium-Sulfur Batteries. ADVANCED MATERIALS (DEERFIELD BEACH, FLA.) 2024; 36:e2306239. [PMID: 37740905 DOI: 10.1002/adma.202306239] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/28/2023] [Revised: 09/08/2023] [Indexed: 09/25/2023]
Abstract
Mg-S batteries hold great promise as a potential alternative to Li-based technologies. Their further development hinges on solving a few key challenges, including the lower capacity and poorer cycling performance when compared to Li counterparts. At the heart of the issues is the lack of knowledge on polysulfide chemical behaviors in the Mg-S battery environment. In this Review, a comprehensive overview of the current understanding of polysulfide behaviors in Mg-S batteries is provided. First, a systematic summary of experimental and computational techniques for polysulfide characterization is provided. Next, conversion pathways for Mg polysulfide species within the battery environment are discussed, highlighting the important role of polysulfide solubility in determining reaction kinetics and overall battery performance. The focus then shifts to the negative effects of polysulfide shuttling on Mg-S batteries. The authors outline various strategies for achieving an optimal balance between polysulfide solubility and shuttling, including the use of electrolyte additives, polysulfide-trapping materials, and dual-functional catalysts. Based on the current understanding, the directions for further advancing knowledge of Mg polysulfide chemistry are identified, emphasizing the integration of experiment with computation as a powerful approach to accelerate the development of Mg-S battery technology.
Collapse
Affiliation(s)
- Tongtong Luo
- Department of Chemistry, Boston College, Chestnut Hill, MA, 02467, USA
| | - Yang Wang
- Department of Chemistry, Boston College, Chestnut Hill, MA, 02467, USA
| | - Brooke Elander
- Department of Chemistry, Boston College, Chestnut Hill, MA, 02467, USA
| | - Michael Goldstein
- Department of Chemistry, Boston College, Chestnut Hill, MA, 02467, USA
| | - Yu Mu
- Department of Chemistry, Boston College, Chestnut Hill, MA, 02467, USA
| | - James Wilkes
- Department of Chemistry, Boston College, Chestnut Hill, MA, 02467, USA
| | | | - Justin Lee
- Department of Chemistry, Boston College, Chestnut Hill, MA, 02467, USA
| | - Tevin Li
- Department of Chemistry, Boston College, Chestnut Hill, MA, 02467, USA
| | - Junwei Lucas Bao
- Department of Chemistry, Boston College, Chestnut Hill, MA, 02467, USA
| | - Udayan Mohanty
- Department of Chemistry, Boston College, Chestnut Hill, MA, 02467, USA
| | - Dunwei Wang
- Department of Chemistry, Boston College, Chestnut Hill, MA, 02467, USA
| |
Collapse
|
39
|
Monge-Palacios M, Wang Q, Alshaarawi A, Sepulveda ACC, Sarathy SM. Quantum chemistry and kinetics of hydrogen sulphide oxidation. Phys Chem Chem Phys 2024; 26:3219-3228. [PMID: 38193631 DOI: 10.1039/d3cp04535h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2024]
Abstract
A fundamental understanding of the acid gas (H2S and CO2) chemistry is key to efficiently implement the desulphurisation process and even the production of clean fuels such as hydrogen or syngas. In this work, we developed a new kinetic model for the pyrolysis and oxidation of hydrogen sulphide by merging two previously reported models with the goal of covering a wider range of conditions and including the effect of carbon dioxide. The resulting model, which consists of 75 species and 514 reactions, was used to conduct rate of production and sensitivity analysis in plug flow reactor simulations, and the results were used to determine the most prominent reactions in which hydrogen sulphide, molecular hydrogen, and sulphur monoxide are involved. The resulting list of important reactions was screened and the kinetics of three of them, i.e., SO2 + S2 → S2O + SO, S2O + S2 → S3 + SO, and SO + SH → S2 + OH, was found to warrant further investigation. With the goal of improving the accurancy of our new kinetic model, we carried out a robust quantum chemistry and Rice-Ramsperger-Kassel-Marcus master equation study to obtain, for the first time, the forward and reverse rate constants for those three reactions at temperatures and pressures of interest for combustion and atmospheric chemistry. This work is the first step of a kinetic study that is aimed at improving the understanding of the chemistry of the pyrolysis and oxidation of H2S, highlighting the importance of sulphur-sulphur interactions and providing a fundamental basis for future kinetic models of H2S not only in the field of combustion, but also in atmospheric chemistry.
Collapse
Affiliation(s)
- M Monge-Palacios
- King Abdullah University of Science and Technology, Thuwal (Makkah) 23955, Saudi Arabia.
| | - Q Wang
- King Abdullah University of Science and Technology, Thuwal (Makkah) 23955, Saudi Arabia.
| | - A Alshaarawi
- Exploration and Petroleum Engineering Center-Advanced Research Center (EXPEC ARC), Saudi Aramco, Dhahran 34465, Saudi Arabia
| | - A C Cavazos Sepulveda
- Exploration and Petroleum Engineering Center-Advanced Research Center (EXPEC ARC), Saudi Aramco, Dhahran 34465, Saudi Arabia
| | - S M Sarathy
- King Abdullah University of Science and Technology, Thuwal (Makkah) 23955, Saudi Arabia.
| |
Collapse
|
40
|
Shiels OJ, Marlton SJP, Poad BLJ, Blanksby SJ, da Silva G, Trevitt AJ. Gas-Phase Phenyl Radical + O 2 Reacts via a Submerged Transition State. J Phys Chem A 2024; 128:413-419. [PMID: 38174881 DOI: 10.1021/acs.jpca.3c06878] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2024]
Abstract
In the gas-phase chemistry of the atmosphere and automotive fuel combustion, peroxyl radical intermediates are formed following O2 addition to carbon-centered radicals which then initiate a complex network of radical reactions that govern the oxidative processing of hydrocarbons. The rapid association of the phenyl radical-a fundamental radical related to benzene-with O2 has hitherto been modeled as a barrierless process, a common assumption for peroxyl radical formation. Here, we provide an alternate explanation for the kinetics of this reaction by deploying double-hybrid density functional theory (DFT), at the DSD-PBEP86-D3(BJ)/aug-cc-pVTZ level of theory, and locate a submerged adiabatic transition state connected to a prereaction complex along the reaction entrance pathway. Using this potential energy scheme, experimental rate coefficients k(T) for the addition of O2 to the phenyl radical are accurately reproduced within a microcanonical kinetic model. This work highlights that purportedly barrierless radical oxidation reactions may instead be modeled using stationary points, which in turn provides insight into pressure and temperature dependence.
Collapse
Affiliation(s)
- Oisin J Shiels
- Molecular Horizons and School of Chemistry and Molecular Bioscience, University of Wollongong, Wollongong 2522, Australia
| | - Samuel J P Marlton
- Molecular Horizons and School of Chemistry and Molecular Bioscience, University of Wollongong, Wollongong 2522, Australia
| | - Berwyck L J Poad
- School of Chemistry and Physics and the Central Analytical Research Facility, Queensland University of Technology, Brisbane 4001, Australia
| | - Stephen J Blanksby
- School of Chemistry and Physics and the Central Analytical Research Facility, Queensland University of Technology, Brisbane 4001, Australia
| | - Gabriel da Silva
- Department of Chemical Engineering, the University of Melbourne, Melbourne 3010, Victoria, Australia
| | - Adam J Trevitt
- Molecular Horizons and School of Chemistry and Molecular Bioscience, University of Wollongong, Wollongong 2522, Australia
| |
Collapse
|
41
|
Fusè M, Mazzeo G, Bloino J, Longhi G, Abbate S. Pushing measurements and interpretation of VCD spectra in the IR, NIR and visible ranges to the detectability and computational complexity limits. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2024; 305:123496. [PMID: 37832448 DOI: 10.1016/j.saa.2023.123496] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/05/2023] [Revised: 08/29/2023] [Accepted: 10/04/2023] [Indexed: 10/15/2023]
Abstract
(R)-Limonene VCD and IR absorption spectra for neat liquid samples are considered from 900 to 16,000 cm-1, using mostly data by Nafie et al. up to 10,000 cm-1 and from previous investigations of the Brescia group. New VCD data are recorded in the merely overtone and combination region between 1800 and 2400 cm-1 and for the Δn= 6 overtone CH-stretching region above 15,000 cm-1. The GVPT2 anharmonic DFT calculations permit satisfactory interpretation of the fundamental + overtone/combination of deformation modes in the mid-IR up to 3500 cm-1. The GVPT2 approach is also used for the first CH-stretching overtone regions together with their combination with deformation modes up to 9000 cm-1. Then the local-mode approach developed within the DFT protocol is employed in all the CH-stretching regions (fundamental + overtones) and is found to satisfactorily account for the observed spectra, justifying the constant VCD pattern observed for all overtones. On the basis of the local-mode model the components of the bisignate VCD spectrum are attributed to the stretchings of the axial and equatorial CH bonds in α-position with respect to the ring CC double bond.
Collapse
Affiliation(s)
- Marco Fusè
- Dipartimento di Medicina Molecolare e Traslazionale, Università di Brescia, Viale Europa 11, 25123, Brescia, Italy
| | - Giuseppe Mazzeo
- Dipartimento di Medicina Molecolare e Traslazionale, Università di Brescia, Viale Europa 11, 25123, Brescia, Italy
| | - Julien Bloino
- Scuola Normale Superiore, Piazza dei Cavalieri, 56125, Pisa, Italy
| | - Giovanna Longhi
- Dipartimento di Medicina Molecolare e Traslazionale, Università di Brescia, Viale Europa 11, 25123, Brescia, Italy; Istituto Nazionale di Ottica (INO), CNR, Research Unit of Brescia, c/o CSMT, VIA Branze 45, 25123, Brescia, Italy
| | - Sergio Abbate
- Dipartimento di Medicina Molecolare e Traslazionale, Università di Brescia, Viale Europa 11, 25123, Brescia, Italy; Istituto Nazionale di Ottica (INO), CNR, Research Unit of Brescia, c/o CSMT, VIA Branze 45, 25123, Brescia, Italy.
| |
Collapse
|
42
|
Deng DD, Long B. Quantitative kinetics of the atmospheric reaction between isocyanic acid and hydroxyl radicals: post-CCSD(T) contribution, anharmonicity, recrossing effects, torsional anharmonicity, and tunneling. Phys Chem Chem Phys 2023; 26:485-492. [PMID: 38079149 DOI: 10.1039/d3cp04385a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2023]
Abstract
Hydroxyl radicals (OH) are the most important atmospheric oxidant, initiating atmospheric reactions for the chemical transformation of volatile organic compounds. Here, we choose the HNCO + OH reaction as a prototype reaction because it contains the fundamental reaction processes for OH radicals: H-abstraction reaction by OH and OH addition reaction. However, its kinetics are unknown under atmospheric conditions. We investigate the reaction of HNCO with OH by using the GMM(P).L method close to the accuracy of single, double, triple, and quadruple excitations and noniterative quintuple excitations with a complete basis set (CCSDTQ(P)/CBS) as benchmark results and a dual-level strategy for kinetics calculations. The calculated rate constant of HNCO + OH is in good agreement with the experimental data available at the temperatures between 620 and 2500 K. We find that the rate constant cannot be correctly obtained by using experimental data to extrapolate the atmospheric temperature ranges. We find that the post-CCSD(T) contribution is very large for the barrier height with the value of -0.85 kcal mol-1 for the H-abstraction reaction, while the previous investigations were done up to the CCSD(T) level. Moreover, we also find that recrossing effects, tunneling, torsional anharmonicity, and anharmonicity are important for obtaining quantitative kinetics in the OH + HNCO reaction.
Collapse
Affiliation(s)
- Dai-Dan Deng
- College of Physics and Mechatronic Engineering, Guizhou Minzu University, Guiyang 550025, China.
| | - Bo Long
- College of Physics and Mechatronic Engineering, Guizhou Minzu University, Guiyang 550025, China.
- College of Materials Science and Engineering, Guizhou Minzu University, Guiyang 550025, China
| |
Collapse
|
43
|
Barone V, Lazzari F. Hunting for Complex Organic Molecules in the Interstellar Medium: The Role of Accurate Low-Cost Theoretical Geometries and Rotational Constants. J Phys Chem A 2023; 127:10517-10527. [PMID: 38033327 PMCID: PMC10726368 DOI: 10.1021/acs.jpca.3c06649] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2023] [Revised: 11/07/2023] [Accepted: 11/08/2023] [Indexed: 12/02/2023]
Abstract
A new approach to computation at affordable cost of accurate geometrical structures and rotational constants for medium-sized molecules in the gas phase is further improved and applied to a large panel of interstellar complex organic molecules. The most distinctive feature of the new model is the effective inclusion of core-valence correlation and vibrational averaging effects in the framework of density functional theory (DFT). In particular, a double-hybrid functional in conjunction with a quadruple-ζ valence/triple-ζ polarization basis set is employed for geometry optimizations, whereas a cheaper hybrid functional in conjunction with a split-valence basis set is used for the evaluation of vibrational corrections. A thorough benchmark based on a wide range of prototypical systems shows that the new scheme approaches the accuracy of state-of-the-art wave function methods with the computational cost of the standard methods (DFT or MP2) routinely employed in the interpretation of microwave spectra. Since the whole computational workflow involves the postprocessing of the output of standard electronic structure codes by a new freely available web utility, the way is paved for the accurate yet not prohibitively expensive study of medium- to large-sized molecules also by nonspecialists.
Collapse
Affiliation(s)
- Vincenzo Barone
- Scuola Normale Superiore, di Pisa, Piazza dei Cavalieri 7, Pisa 56125, Italy
| | - Federico Lazzari
- Scuola Normale Superiore, di Pisa, Piazza dei Cavalieri 7, Pisa 56125, Italy
| |
Collapse
|
44
|
Hancock AC, Goerigk L. Noncovalently bound excited-state dimers: a perspective on current time-dependent density functional theory approaches applied to aromatic excimer models. RSC Adv 2023; 13:35964-35984. [PMID: 38090083 PMCID: PMC10712016 DOI: 10.1039/d3ra07381e] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2023] [Accepted: 11/21/2023] [Indexed: 05/12/2024] Open
Abstract
Excimers are supramolecular systems whose binding strength is influenced by many factors that are ongoing challenges for computational methods, such as charge transfer, exciton coupling, and London dispersion interactions. Treating the various intricacies of excimer binding at an adequate level is expected to be particularly challenging for time-dependent Density Functional Theory (TD-DFT) methods. In addition to well-known limitations for some TD-DFT methods in the description of charge transfer or exciton coupling, the inherent London dispersion problem from ground-state DFT translates to TD-DFT. While techniques to appropriately treat dispersion in DFT are well-developed for electronic ground states, these dispersion corrections remain largely untested for excited states. Herein, we aim to shed light on current TD-DFT methods, including some of the newest developments. The binding of four model excimers is studied across nine density functionals with and without the application of additive dispersion corrections against a wave function reference of SCS-CC2/CBS(3,4) quality, which approximates select CCSDR(3)/CBS data adequately. To our knowledge, this is the first study that presents single-reference wave function dissociation curves at the complete basis set level for the assessed model systems. It is also the first time range-separated double-hybrid density functionals are applied to excimers. In fact, those functionals turn out to be the most promising for the description of excimer binding followed by global double hybrids. Range-separated and global hybrids-particularly with large fractions of Fock exchange-are outperformed by double hybrids and yield worse dissociation energies and inter-molecular equilibrium distances. The deviation between each assessed functional and reference increases with system size, most likely due to missing dispersion interactions. Additive dispersion corrections of the DFT-D3(BJ) and DFT-D4 types reduce the average errors for TD-DFT methods but do so inconsistently and therefore do not offer a black-box solution in their ground-state parametrised form. The lack of appropriate description of dispersion effects for TD-DFT methods is likely hindering the practical application of the herein identified more efficient methods. Dispersion corrections parametrised for excited states appear to be an important next step to improve the applicability of TD-DFT methods and we hope that our work assists with the future development of such corrections.
Collapse
Affiliation(s)
- Amy C Hancock
- School of Chemistry, The University of Melbourne Parkville Australia +61-(0)3-8344 6784
| | - Lars Goerigk
- School of Chemistry, The University of Melbourne Parkville Australia +61-(0)3-8344 6784
| |
Collapse
|
45
|
Pitman SJ, Evans AK, Ireland RT, Lempriere F, McKemmish LK. Benchmarking Basis Sets for Density Functional Theory Thermochemistry Calculations: Why Unpolarized Basis Sets and the Polarized 6-311G Family Should Be Avoided. J Phys Chem A 2023; 127:10295-10306. [PMID: 37982604 DOI: 10.1021/acs.jpca.3c05573] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2023]
Abstract
Basis sets are a crucial but often largely overlooked choice in setting up quantum chemistry calculations. The choice of the basis set can be critical in determining the accuracy and calculation time of your quantum chemistry calculations. Clear recommendations based on thorough benchmarking are essential but not readily available currently. This study investigates the relative quality of basis sets for general properties by benchmarking basis set performance for a diverse set of 139 reactions (from the diet-150-GMTKN55 data set). In our analysis, we find the distributions of errors are often significantly non-Gaussian, meaning that the joint consideration of median errors, mean absolute errors, and outlier statistics is helpful to provide a holistic understanding of basis set performance. Our direct comparison of performance between most modern basis sets provides quantitative evidence for basis set recommendations that broadly align with the established understanding of basis set experts and is evident in the design of modern basis sets. For example, while zeta is a good measure of quality, it is not the only determining factor for an accurate calculation with unpolarized double- and triple-ζ basis sets (like 6-31G and 6-311G) having very poor performance. Appropriate use of polarization functions (e.g., 6-31G*) is essential to obtain the accuracy offered by double- or triple-ζ basis sets. In our study, the best performances for double- and triple-ζ basis sets are 6-31++G** and pcseg-2, respectively. However, the performances of singly polarized double-ζ and doubly polarized triple-ζ basis sets are quite similar with one key exception: the polarized 6-311G basis set family has poor parametrization, which means its performance is more like a double-ζ than a triple-ζ basis set. All versions of the 6-311G basis set family should be avoided entirely for valence chemistry calculations moving forward.
Collapse
Affiliation(s)
- Samuel J Pitman
- School of Chemistry, University of New South Wales, Sydney, NSW 2052, Australia
| | - Alicia K Evans
- School of Chemistry, University of New South Wales, Sydney, NSW 2052, Australia
| | - Robbie T Ireland
- School of Chemistry, University of New South Wales, Sydney, NSW 2052, Australia
| | - Felix Lempriere
- School of Chemistry, University of New South Wales, Sydney, NSW 2052, Australia
| | - Laura K McKemmish
- School of Chemistry, University of New South Wales, Sydney, NSW 2052, Australia
| |
Collapse
|
46
|
King DS, Truhlar DG, Gagliardi L. Variational Active Space Selection with Multiconfiguration Pair-Density Functional Theory. J Chem Theory Comput 2023; 19:8118-8128. [PMID: 37905518 DOI: 10.1021/acs.jctc.3c00792] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2023]
Abstract
The selection of an adequate set of active orbitals for modeling strongly correlated electronic states is difficult to automate because it is highly dependent on the states and molecule of interest. Although many approaches have shown some success, no single approach has worked well in all cases. In light of this, we present the "discrete variational selection" (DVS) approach to active space selection, in which one generates multiple trial wave functions from a diverse set of systematically constructed active spaces and then selects between these wave functions variationally. We apply this DVS approach to 207 vertical excitations of small-to-medium-sized organic and inorganic molecules (with 3 to 18 atoms) in the QUESTDB database by (i) constructing various sets of active space orbitals through diagonalization of parametrized operators and (ii) choosing the result with the lowest average energy among the states of interest. This approach proves ineffective when variationally selecting between wave functions using the density matrix renormalization group (DMRG) or complete active space self-consistent field (CASSCF) energy but is able to provide good results when variationally selecting between wave functions using the energy of the translated PBE (tPBE) functional from multiconfiguration pair-density functional theory (MC-PDFT). Applying this DVS-tPBE approach to selection among state-averaged DMRG wave functions, we obtain a mean unsigned error of only 0.17 eV using hybrid MC-PDFT. This result matches that of our previous benchmark without the need to filter out poor active spaces and with no further orbital optimization following active space selection of the SA-DMRG wave functions. Furthermore, we find that DVS-tPBE is able to robustly and effectively select between the new SA-DMRG wave functions and our previous SA-CASSCF results.
Collapse
Affiliation(s)
- Daniel S King
- Department of Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Group, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Laura Gagliardi
- Department of Chemistry, Pritzker School of Molecular Engineering, James Franck Institute, Chicago Center for Theoretical Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
47
|
Melosso M, Alessandrini S, Spada L, Melli A, Wang X, Zheng Y, Duan C, Li J, Du W, Gou Q, Bizzocchi L, Dore L, Barone V, Puzzarini C. Rotational spectra and semi-experimental structures of furonitrile and its water cluster. Phys Chem Chem Phys 2023; 25:31281-31291. [PMID: 37955344 DOI: 10.1039/d3cp03984f] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2023]
Abstract
Rotational spectroscopy represents an invaluable tool for several applications: from the identification of new molecules in interstellar objects to the characterization of van der Waals complexes, but also for the determination of very accurate molecular structures and for conformational analyses. In this work, we used high-resolution rotational spectroscopic techniques in combination with high-level quantum-chemical calculations to address all these aspects for two isomers of cyanofuran, namely 2-furonitrile and 3-furonitrile. In particular, we have recorded and analyzed the rotational spectra of both of them from 6 to 320 GHz; rotational transitions belonging to several singly-substituted isotopologues have been identified as well. The rotational constants derived in this way have been used in conjunction with computed rotation-vibration interaction constants in order to derive a semi-experimental equilibrium structure for both isomers. Moreover, we observed the rotational spectra of four different intermolecular adducts formed by furonitrile and water, whose identification has been supported by a conformational analysis and a theoretical spectroscopic characterization. A semi-experimental determination of the intermolecular parameters has been achieved for all of them and the results have been compared with those obtained for the analogous system formed by benzonitrile and water.
Collapse
Affiliation(s)
- Mattia Melosso
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy.
| | - Silvia Alessandrini
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy.
| | - Lorenzo Spada
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Alessio Melli
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy.
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Xiujuan Wang
- Department of Chemistry, School of Chemistry and Chemical Engineering, Chongqing University, Daxuecheng South Rd. 55, 401331, Chongqing, China.
| | - Yang Zheng
- Department of Chemistry, School of Chemistry and Chemical Engineering, Chongqing University, Daxuecheng South Rd. 55, 401331, Chongqing, China.
| | - Chunguo Duan
- Department of Chemistry, School of Chemistry and Chemical Engineering, Chongqing University, Daxuecheng South Rd. 55, 401331, Chongqing, China.
| | - Jiayi Li
- Department of Chemistry, School of Chemistry and Chemical Engineering, Chongqing University, Daxuecheng South Rd. 55, 401331, Chongqing, China.
| | - Weiping Du
- Department of Chemistry, School of Chemistry and Chemical Engineering, Chongqing University, Daxuecheng South Rd. 55, 401331, Chongqing, China.
| | - Qian Gou
- Department of Chemistry, School of Chemistry and Chemical Engineering, Chongqing University, Daxuecheng South Rd. 55, 401331, Chongqing, China.
| | - Luca Bizzocchi
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy.
| | - Luca Dore
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy.
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Cristina Puzzarini
- Dipartimento di Chimica "Giacomo Ciamician", Università di Bologna, Via F. Selmi 2, 40126 Bologna, Italy.
| |
Collapse
|
48
|
Crisci L, Di Grande S, Cavallotti C, Barone V. Toward an Accurate Black-Box Tool for the Kinetics of Gas-Phase Reactions Involving Barrier-less Elementary Steps. J Chem Theory Comput 2023; 19:7626-7639. [PMID: 37880932 PMCID: PMC10653117 DOI: 10.1021/acs.jctc.3c00857] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2023] [Revised: 10/06/2023] [Accepted: 10/13/2023] [Indexed: 10/27/2023]
Abstract
An enhanced computational protocol has been devised for the accurate characterization of gas-phase barrier-less reactions in the framework of the reaction-path (RP) and variable reaction coordinate variational transition-state theory. In particular, the synergistic combination of density functional theory and Monte Carlo sampling to optimize reactive fluxes led to a reliable yet effective computational workflow. A black-box strategy has been developed for selecting the most suited density functional with reference to a high-level one-dimensional reference potential. At the same time, different descriptions of hindered rotations are automatically selected, depending on the corresponding harmonic frequencies along the RP. The performance of the new tool is investigated by means of two prototypical reactions involving different degrees of static and dynamic correlation, namely, H2S + Cl and CH3 + CH3. The remarkable agreement of the computed kinetic parameters with the available experimental data confirms the accuracy and robustness of the proposed approach. Together with their intrinsic interest, these results also pave the way toward systematic investigations of gas-phase reactions involving barrier-less elementary steps by a reliable, user-friendly tool, which can be confidently used also by nonspecialists.
Collapse
Affiliation(s)
- Luigi Crisci
- Scuola
Normale Superiore di Pisa, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
- Department
of Chemical Sciences, University of Napoli
Federico II, Complesso Universitario di M.S. Angelo, via Cintia 21, I-80126 Napoli, Italy
| | - Silvia Di Grande
- Scuola
Normale Superiore di Pisa, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
- Scuola
Superiore Meridionale, Largo San Marcellino 10, I-80138 Napoli, Italy
| | - Carlo Cavallotti
- Department
of Chemistry, Materials and Chemical Engineering “G. Natta”, Politecnico di Milano, I-20131 Milano, Italy
| | - Vincenzo Barone
- Scuola
Normale Superiore di Pisa, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| |
Collapse
|
49
|
Hazrah AS, Insausti A, Ma J, Al-Jabiri MH, Jäger W, Xu Y. Wetting vs Droplet Aggregation: A Broadband Rotational Spectroscopic Study of 3-Methylcatechol⋅⋅⋅Water Clusters. Angew Chem Int Ed Engl 2023; 62:e202310610. [PMID: 37697450 DOI: 10.1002/anie.202310610] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2023] [Revised: 09/07/2023] [Accepted: 09/11/2023] [Indexed: 09/13/2023]
Abstract
Two competing solvation pathways of 3-methylcatechol (MC), an atmospherically relevant aromatic molecule, with up to five water molecules were explored in detail by using a combination of broadband rotational spectroscopy and computational chemistry. Theoretically, two different pathways of solvation emerge: the commonly observed droplet pathway which involves preferential binding among the water molecules while the solute serves as an anchor point for the formation of a water cluster, and an unexpected wetting pathway which involves interactions between the water molecules and the aromatic face of MC, i.e., a wetting of the π-surface. Conclusive identification of the MC hydrate structures, and therefore the wetting pathway, was facilitated by rotational spectra of the parent MC hydrates and several H2 18 O and 13 C isotopologues which exhibit splittings associated with methyl internal rotation and/or water tunneling motions. Theoretical modelling and analyses offer insights into the tunneling and conversion barriers associated with the observed hydrate conformers and the nature of the non-covalent interactions involved in choosing the unusual wetting pathway.
Collapse
Affiliation(s)
- Arsh S Hazrah
- Department of Chemistry, University of Alberta, Edmonton, AB T6G 2G2, Canada
- Current Address: Department of Molecular Spectroscopy, Max Planck Institute for Polymer Research, Ackermannweg 10, 55128, Mainz, Germany
| | - Aran Insausti
- Department of Chemistry, University of Alberta, Edmonton, AB T6G 2G2, Canada
- Departamento de Química Física, Facultad de Ciencia y Tecnología, Universidad del País Vasco (UPV-EHU), 48080, Bilbao, Spain
- Biofisika Institute (CSIC, UPV/EHU), University of the Basque Country (UPV/EHU), 48940, Leioa, Spain
| | - Jiarui Ma
- Department of Chemistry, University of Alberta, Edmonton, AB T6G 2G2, Canada
| | - Mohamad H Al-Jabiri
- Department of Chemistry, University of Alberta, Edmonton, AB T6G 2G2, Canada
| | - Wolfgang Jäger
- Department of Chemistry, University of Alberta, Edmonton, AB T6G 2G2, Canada
| | - Yunjie Xu
- Department of Chemistry, University of Alberta, Edmonton, AB T6G 2G2, Canada
| |
Collapse
|
50
|
Olive LN, Dornshuld EV, Schaefer HF, Tschumper GS. Competition between Solvent···Solvent and Solvent···Solute Interactions in the Microhydration of the Tetrafluoroborate Anion, BF 4-(H 2O) n=1,2,3,4. J Phys Chem A 2023; 127:8806-8820. [PMID: 37774368 DOI: 10.1021/acs.jpca.3c04014] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/01/2023]
Abstract
This study systematically examines the interactions of the tetrafluoroborate anion (BF4-) with up to four water molecules (BF4-(H2O)n=1,2,3,4). Full geometry optimizations and subsequent harmonic vibrational frequency computations are performed using a variety of density functional theory (DFT) methods (B3LYP, B3LYP-D3BJ, and M06-2X) and the MP2 ab initio method with a triple-ζ correlation consistent basis set augmented with diffuse functions on all non-hydrogen atoms (cc-pVTZ for H and aug-cc-pVTZ for B, O, and F; denoted as haTZ). Optimized structures and harmonic vibrational frequencies were also obtained with the CCSD(T) ab initio method and the haTZ basis set for the mono- and dihydrate (n = 1, 2) structures. The 2-body:Many-body (2b:Mb) technique, in which CCSD(T) computations capture the 1- and 2-body contributions to the interactions and MP2 computations recover all higher-order contributions, was used to extend these demanding computations to the tri- and tetrahydrate (n = 3, 4) systems. Four, five, and eight new stationary points have been identified for the di-, tri-, and tetrahydrate systems, respectively. The global minimum of the monohydrate adopts a symmetric double ionic hydrogen bond motif with C2v symmetry and an electronic dissociation energy of 13.17 kcal mol-1 at the CCSD(T)/haTZ level of theory. This strong solvent···solute interaction, however, competes with solute···solute interactions in the lowest-energy BF4-(H2O)n=2,3,4 minima that are not seen in the other di-, tri-, or tetrahydrate minima. The latter interactions help increase the 2b:Mb dissociation energies to more than 26, 41, and 51 kcal mol-1 for n = 2, 3, and 4, respectively. Structures that form hydrogen bonds between the solvating water molecules also exhibit the largest shifts in the harmonic OH stretching frequencies for the waters of hydration. These shifts can exceed -280 cm-1 relative to an isolated H2O molecule at the 2b:Mb/haTZ level of theory.
Collapse
Affiliation(s)
- Laura N Olive
- Center for Computational Quantum Chemistry, University of Georgia, Athens, Georgia 30602, United States
| | - Eric V Dornshuld
- Department of Chemistry, Mississippi State University, Mississippi State, Mississippi 39762, United States
| | - Henry F Schaefer
- Center for Computational Quantum Chemistry, University of Georgia, Athens, Georgia 30602, United States
| | - Gregory S Tschumper
- Department of Chemistry and Biochemistry, University of Mississippi, University, Mississippi 38677, United States
| |
Collapse
|