1
|
Chan B. Limiting factors in the accuracy of DFT calculation for redox potentials. J Comput Chem 2024; 45:1177-1186. [PMID: 38311976 DOI: 10.1002/jcc.27320] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2023] [Revised: 01/17/2024] [Accepted: 01/18/2024] [Indexed: 02/06/2024]
Abstract
In the present study, we have investigated factors affecting the accuracy of computational chemistry calculation of redox potentials, namely the gas-phase ionization energy (IE) and electron affinity (EA), and the continuum solvation effect. In general, double-hybrid density functional theory methods yield IEs and EAs that are on average within ~0.1 eV of our high-level W3X-L benchmark, with the best performing method being DSD-BLYP/ma-def2-QZVPP. For lower-cost methods, the average errors are ~0.2-0.3 eV, with ωB97X-3c being the most accurate (~0.15 eV). For the solvation component, essentially all methods have an average error of ~0.3 eV, which shows the limitation of the continuum solvation model. Curiously, the directly calculated redox potentials show errors of ~0.3 eV for all methods. These errors are notably smaller than what can be expected from error propagation with the two components (IE and EA, and solvation effect). Such a discrepancy can be attributed to the cancellation of errors, with the lowest-cost GFN2-xTB method benefiting the most, and the most accurate ωB97X-3c method benefiting the least. For organometallic species, the redox potentials show large deviations exceeding ~0.5 eV even for DSD-BLYP. The large errors are attributed to those for the gas-phase IEs and EAs, which represents a major barrier to the accurate calculation of redox potentials for such systems.
Collapse
Affiliation(s)
- Bun Chan
- Graduate School of Engineering, Nagasaki University, Nagasaki, Japan
- RIKEN Center for Computational Science, Kobe, Japan
| |
Collapse
|
2
|
Semidalas E, Karton A, Martin JML. W4Λ: Leveraging Λ Coupled-Cluster for Accurate Computational Thermochemistry Approaches. J Phys Chem A 2024; 128:1715-1724. [PMID: 38400740 PMCID: PMC10926103 DOI: 10.1021/acs.jpca.3c08158] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2023] [Revised: 01/28/2024] [Accepted: 01/31/2024] [Indexed: 02/26/2024]
Abstract
High-accuracy composite wave function methods like Weizmann-4 (W4) theory, high-accuracy extrapolated ab initio thermochemistry (HEAT), and the Feller-Peterson-Dixon (FPD) approach enable sub-kJ/mol accuracy in gas-phase thermochemical properties. Their biggest computational bottleneck is the evaluation of the valence post-CCSD(T) correction term. We demonstrate here, for the W4-17 thermochemistry benchmark and subsets thereof, that the Λ coupled-cluster expansion converges more rapidly and smoothly than the regular coupled-cluster series. By means of CCSDT(Q)Λ and CCSDTQ(5)Λ, we can considerably (up to an order of magnitude) accelerate W4- and W4.3-type calculations without loss in accuracy, leading to the W4Λ and W4.3Λ computational thermochemistry protocols.
Collapse
Affiliation(s)
- Emmanouil Semidalas
- Department
of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| | - Amir Karton
- School
of Science and Technology, University of
New England, Armidale, New South Wales 2351, Australia
| | - Jan M. L. Martin
- Department
of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| |
Collapse
|
3
|
Salo VT, Chen J, Runeberg N, Kjaergaard HG, Kurtén T. Multireference and Coupled-Cluster Study of Dimethyltetroxide (MeO 4Me) Formation and Decomposition. J Phys Chem A 2024. [PMID: 38417845 DOI: 10.1021/acs.jpca.3c08043] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/01/2024]
Abstract
Peroxyl radicals (RO2) are important intermediates in the atmospheric oxidation processes. The RO2 can react with other RO2 to form reactive intermediates known as tetroxides, RO4R. The reaction mechanisms of RO4R formation and its various decomposition channels have been investigated in multiple computational studies, but previous approaches have not been able to provide a unified methodology that is able to connect all relevant reactions together. An apparent difficulty in modeling the RO4R formation and its decomposition is the involvement of open-shell singlet electronic states along the reaction pathway. Modeling such electronic states requires multireference (MR) methods, which we use in the present study to investigate in detail a model reaction of MeO2 + MeO2 → MeO4Me, and its decomposition, MeO4Me → MeO + O2 + MeO, as well as the two-body product complexes MeO···O2 + MeO and MeO···MeO + O2. The used MR methods are benchmarked against relative energies of MeO2 + MeO2, MeO4Me, and MeO + MeO + O2, calculated with CCSD(T)/CBS, W2X, and W3X-L methods. We found that the calculated relative energies of the overall MeO2 + MeO2 → MeO4Me → MeO + O2 + MeO reaction are very sensitive to the chosen MR method and that the CASPT2(22e,14o)-IPEA method is able to reproduce the relative energies obtained by the various coupled-cluster methods. Furthermore, CASPT2(22e,14o)-IPEA and W3X-L results show that the MeO···O2 product complex + MeO is lower in energy than the MeO···MeO complex + O2. The formation of MeO4Me is exothermic, and its decomposition is endothermic, relative to the tetroxide. Both the formation and decomposition reactions appear to follow pathways with no saddle points. According to potential energy surface scans of the decomposition pathway, the concerted cleavage of both MeO···O bonds in MeO4Me is energetically preferred over the corresponding sequential decomposition.
Collapse
Affiliation(s)
- Vili-Taneli Salo
- Department of Chemistry, Faculty of Science, University of Helsinki, Helsinki FI-00014, Finland
| | - Jing Chen
- Department of Chemistry, University of Copenhagen, Copenhagen 2100, Denmark
| | - Nino Runeberg
- Department of Chemistry, Faculty of Science, University of Helsinki, Helsinki FI-00014, Finland
| | | | - Theo Kurtén
- Department of Chemistry, Faculty of Science, University of Helsinki, Helsinki FI-00014, Finland
| |
Collapse
|
4
|
Yang X, Wang H, Lu K, Ma X, Tan Z, Long B, Chen X, Li C, Zhai T, Li Y, Qu K, Xia Y, Zhang Y, Li X, Chen S, Dong H, Zeng L, Zhang Y. Reactive aldehyde chemistry explains the missing source of hydroxyl radicals. Nat Commun 2024; 15:1648. [PMID: 38388476 PMCID: PMC10883920 DOI: 10.1038/s41467-024-45885-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2023] [Accepted: 02/06/2024] [Indexed: 02/24/2024] Open
Abstract
Hydroxyl radicals (OH) determine the tropospheric self-cleansing capacity, thus regulating air quality and climate. However, the state-of-the-art mechanisms still underestimate OH at low nitrogen oxide and high volatile organic compound regimes even considering the latest isoprene chemistry. Here we propose that the reactive aldehyde chemistry, especially the autoxidation of carbonyl organic peroxy radicals (R(CO)O2) derived from higher aldehydes, is a noteworthy OH regeneration mechanism that overwhelms the contribution of the isoprene autoxidation, the latter has been proved to largely contribute to the missing OH source under high isoprene condition. As diagnosed by the quantum chemical calculations, the R(CO)O2 radicals undergo fast H-migration to produce unsaturated hydroperoxyl-carbonyls that generate OH through rapid photolysis. This chemistry could explain almost all unknown OH sources in areas rich in both natural and anthropogenic emissions in the warm seasons, and may increasingly impact the global self-cleansing capacity in a future low nitrogen oxide society under carbon neutrality scenarios.
Collapse
Affiliation(s)
- Xinping Yang
- State Key Joint Laboratory of Environmental Simulation and Pollution Control, State Environmental Protection Key Laboratory of Atmospheric Ozone Pollution Control, College of Environmental Sciences and Engineering, Peking University, Beijing, 100871, China
- State Environmental Protection Key Laboratory of Vehicle Emission Control and Simulation, Vehicle Emission Control Center, Chinese Research Academy of Environmental Sciences, Beijing, 100012, China
| | - Haichao Wang
- School of Atmospheric Sciences, Sun Yat-sen University and Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai), Zhuhai, 519082, China
- Guangdong Provincial Observation and Research Station for Climate Environment and Air Quality Change in the Pearl River Estuary, Key Laboratory of Tropical Atmosphere-Ocean System, Ministry of Education, Zhuhai, 519082, China
| | - Keding Lu
- State Key Joint Laboratory of Environmental Simulation and Pollution Control, State Environmental Protection Key Laboratory of Atmospheric Ozone Pollution Control, College of Environmental Sciences and Engineering, Peking University, Beijing, 100871, China.
| | - Xuefei Ma
- State Key Joint Laboratory of Environmental Simulation and Pollution Control, State Environmental Protection Key Laboratory of Atmospheric Ozone Pollution Control, College of Environmental Sciences and Engineering, Peking University, Beijing, 100871, China
| | - Zhaofeng Tan
- State Key Joint Laboratory of Environmental Simulation and Pollution Control, State Environmental Protection Key Laboratory of Atmospheric Ozone Pollution Control, College of Environmental Sciences and Engineering, Peking University, Beijing, 100871, China
| | - Bo Long
- College of Material Science and Engineering, Guizhou Minzu University, Guizhou, China
| | - Xiaorui Chen
- State Key Joint Laboratory of Environmental Simulation and Pollution Control, State Environmental Protection Key Laboratory of Atmospheric Ozone Pollution Control, College of Environmental Sciences and Engineering, Peking University, Beijing, 100871, China
| | - Chunmeng Li
- State Key Joint Laboratory of Environmental Simulation and Pollution Control, State Environmental Protection Key Laboratory of Atmospheric Ozone Pollution Control, College of Environmental Sciences and Engineering, Peking University, Beijing, 100871, China
| | - Tianyu Zhai
- State Key Joint Laboratory of Environmental Simulation and Pollution Control, State Environmental Protection Key Laboratory of Atmospheric Ozone Pollution Control, College of Environmental Sciences and Engineering, Peking University, Beijing, 100871, China
| | - Yang Li
- State Key Joint Laboratory of Environmental Simulation and Pollution Control, State Environmental Protection Key Laboratory of Atmospheric Ozone Pollution Control, College of Environmental Sciences and Engineering, Peking University, Beijing, 100871, China
| | - Kun Qu
- State Key Joint Laboratory of Environmental Simulation and Pollution Control, State Environmental Protection Key Laboratory of Atmospheric Ozone Pollution Control, College of Environmental Sciences and Engineering, Peking University, Beijing, 100871, China
| | - Yu Xia
- College of Material Science and Engineering, Guizhou Minzu University, Guizhou, China
| | - Yuqiong Zhang
- College of Material Science and Engineering, Guizhou Minzu University, Guizhou, China
| | - Xin Li
- State Key Joint Laboratory of Environmental Simulation and Pollution Control, State Environmental Protection Key Laboratory of Atmospheric Ozone Pollution Control, College of Environmental Sciences and Engineering, Peking University, Beijing, 100871, China
| | - Shiyi Chen
- State Key Joint Laboratory of Environmental Simulation and Pollution Control, State Environmental Protection Key Laboratory of Atmospheric Ozone Pollution Control, College of Environmental Sciences and Engineering, Peking University, Beijing, 100871, China
| | - Huabin Dong
- State Key Joint Laboratory of Environmental Simulation and Pollution Control, State Environmental Protection Key Laboratory of Atmospheric Ozone Pollution Control, College of Environmental Sciences and Engineering, Peking University, Beijing, 100871, China
| | - Limin Zeng
- State Key Joint Laboratory of Environmental Simulation and Pollution Control, State Environmental Protection Key Laboratory of Atmospheric Ozone Pollution Control, College of Environmental Sciences and Engineering, Peking University, Beijing, 100871, China
| | - Yuanhang Zhang
- State Key Joint Laboratory of Environmental Simulation and Pollution Control, State Environmental Protection Key Laboratory of Atmospheric Ozone Pollution Control, College of Environmental Sciences and Engineering, Peking University, Beijing, 100871, China.
| |
Collapse
|
5
|
Zhang YQ, Francisco JS, Long B. Rapid Atmospheric Reactions between Criegee Intermediates and Hypochlorous Acid. J Phys Chem A 2024; 128:909-917. [PMID: 38271208 DOI: 10.1021/acs.jpca.3c06144] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2024]
Abstract
Hypochlorous acid (HOCl) is a paramount compound in the atmosphere due to its significant contribution to both tropospheric oxidation capacity and ozone depletion. The main removal routes for HOCl are photolysis and the reaction with OH during the daytime, while these processes are unimportant during the nighttime. Here, we report the rapid reactions of Criegee intermediates (CH2OO and anti/syn-CH3CHOO) with HOCl by using high-level quantum chemical methods as the benchmark; their accuracy is close to coupled cluster theory with single, double, and triple excitations and quasiperturbative connected quadruple excitations with a complete basis limit by extrapolation [CCSDT(Q)/CBS]. Their rate constants have been calculated by using a dual-level strategy; this combines conventional transition state theory calculated at the benchmark level with variational transition state theory with small-curvature tunneling by a validated density functional method. We find that the introduction of the methyl group into Criegee intermediates not only affects their reactivities but also exerts a remarkable influence on anharmonicity. The calculated results uncover that anharmonicity increases the rate constants of CH2OO + HOCl by a factor of 18-5, while it is of minor importance in the anti/syn-CH3CHOO + HOCl reaction at 190-350 K. The present findings reveal that the loose transition state for anti-CH3CHOO and HOCl is a rate-determining step at 190-350 K. We also find that the reaction of Criegee intermediates with HOCl contributes significantly to the sink of HOCl during the nighttime in the atmosphere.
Collapse
Affiliation(s)
- Yu-Qiong Zhang
- College of Physics and Mechatronic Engineering, Guizhou Minzu University, Guiyang 550025, China
| | - Joseph S Francisco
- Department of Earth and Environmental Sciences and Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - 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
|
6
|
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
|
7
|
Chan B. DAPD Set of Pd-Containing Diatomic Molecules: Accurate Molecular Properties and the Great Lengths to Obtain Them. J Chem Theory Comput 2023; 19:9260-9268. [PMID: 38096563 DOI: 10.1021/acs.jctc.3c01060] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2023]
Abstract
In the present study, we obtained reliable bond energy, bond length, and zero-point vibrational frequency for a set of diatomic Pd species (the DAPD set). It includes PdH, Pd2, and PdX (X = B, C, N, O, F, Al, Si, P, S, and Cl). Our highest-level protocol (W4X-L) represents scalar and spin-orbit relativistic, valence- and inner-valence correlated, extrapolated CCSDTQ(5) energy. The DAPD set of molecules is challenging for computational chemistry methods in different manners; for Pd2, the spin-orbit contribution to the bond energy is fairly large, whereas for PdC and PdSi, the post-CCSD(T) correlation components are considerable. The diverse range of requirements represents a significant challenge for lower-level methods. While density functional theory (DFT) methods generally yield good agreements for bond lengths and vibrational frequencies, large deviations are found for bond energies. In general, hybrid DFT methods are more accurate than nonhybrid functionals, but the agreement in individual cases varies. This illustrates the critical role that new high-quality reference data would play in the continual development of lower-cost methods.
Collapse
Affiliation(s)
- Bun Chan
- Graduate School of Engineering, Nagasaki University, Bunkyo 1-14, Nagasaki 852-8521, Japan
| |
Collapse
|
8
|
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
|
9
|
Zhang Y, Wang Z, Wang H, Cheng Y, Zhang T, Ou T, Wang R. Atmospheric Chemistry of NH 2SO 3H in Polluted Areas: An Unexpected Isomerization of NH 2SO 3H in Acid-Polluted Regions. J Phys Chem A 2023; 127:8935-8942. [PMID: 37844321 DOI: 10.1021/acs.jpca.3c04982] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2023]
Abstract
NH2SO3H is an effective nucleation agent for the formation of atmospheric aerosols and cloud particles. So, the ammonolysis of SO3 to form NH2SO3H without and with neutral (H2O) and basic (NH3) trace gases has been extensively investigated. However, the acidic trace gas X (X = H2SO4 and CH3SO3H)-assisted ammonolysis of SO3 is still up for debate. In this work, a comprehensive theoretical investigation of X-assisted ammonolysis of SO3 and its reverse reaction (the isomerization of NH2SO3H to form SO3-···NH3+) was carried out in the gas phase and at the air-water interface. The gas-phase results show that X-assisted isomerization of NH2SO3H to form SO3-···NH3+ is more energetically and kinetically favorable than its reverse reaction and the isomerization of NH2SO3H in the presence of H2O and NH3. Such unexpected findings revealed that gas-phase NH2SO3H is highly reactive in the presence of acidic trace gas in contrast to the high stability of NH2SO3H in neutral and basic conditions. At the air-water interface, the X-assisted isomerization reaction of NH2SO3H involves multiple water molecules. The loop structure of the reaction center (X···NH2SO3H···3H2O) promotes the transfer of protons in the water molecules to form the SO3-···NH3+ ion pair, which can then interact with several interfacial water molecules to form ammonium bisulfate. These interfacial reaction channels follow a stepwise mechanism and proceed at the picosecond time-scale. The findings of this study will contribute to a better understanding of the atmospheric behavior of NH2SO3H in polluted acidic trace gases.
Collapse
Affiliation(s)
- Yongqi Zhang
- Institute of Theoretical and Computational Chemistry, Shaanxi Key Laboratory of Catalysis, School of Chemical & Environment Science, Shaanxi University of Technology, Hanzhong, Shaanxi 723001, P. R. China
- National and Local Joint Biomedical Engineering Research Center on Photodynamic Technologies, College of Chemistry, Fuzhou University, Fuzhou, Fujian 350116, P. R. China
| | - Zehui Wang
- Institute of Theoretical and Computational Chemistry, Shaanxi Key Laboratory of Catalysis, School of Chemical & Environment Science, Shaanxi University of Technology, Hanzhong, Shaanxi 723001, P. R. China
| | - Hui Wang
- Institute of Theoretical and Computational Chemistry, Shaanxi Key Laboratory of Catalysis, School of Chemical & Environment Science, Shaanxi University of Technology, Hanzhong, Shaanxi 723001, P. R. China
| | - Yang Cheng
- Institute of Theoretical and Computational Chemistry, Shaanxi Key Laboratory of Catalysis, School of Chemical & Environment Science, Shaanxi University of Technology, Hanzhong, Shaanxi 723001, P. R. China
| | - Tianlei Zhang
- Institute of Theoretical and Computational Chemistry, Shaanxi Key Laboratory of Catalysis, School of Chemical & Environment Science, Shaanxi University of Technology, Hanzhong, Shaanxi 723001, P. R. China
| | - Ting Ou
- Institute of Theoretical and Computational Chemistry, Shaanxi Key Laboratory of Catalysis, School of Chemical & Environment Science, Shaanxi University of Technology, Hanzhong, Shaanxi 723001, P. R. China
| | - Rui Wang
- Institute of Theoretical and Computational Chemistry, Shaanxi Key Laboratory of Catalysis, School of Chemical & Environment Science, Shaanxi University of Technology, Hanzhong, Shaanxi 723001, P. R. China
| |
Collapse
|
10
|
Long B, Xia Y, Zhang YQ, Truhlar DG. Kinetics of Sulfur Trioxide Reaction with Water Vapor to Form Atmospheric Sulfuric Acid. J Am Chem Soc 2023; 145:19866-19876. [PMID: 37651227 DOI: 10.1021/jacs.3c06032] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/02/2023]
Abstract
Although experimental methods can be used to obtain the quantitative kinetics of atmospheric reactions, experimental data are often limited to a narrow temperature range. The reaction of SO3 with water vapor is important for elucidating the formation of sulfuric acid in the atmosphere; however, the kinetics is uncertain at low temperatures. Here, we calculate rate constants for reactions of sulfur trioxide with two water molecules. We consider two mechanisms: the SO3···H2O + H2O reaction and the SO3 + (H2O)2 reaction. We find that beyond-CCSD(T) contributions to the barrier heights are very large, and multidimensional tunneling, unusually large anharmonicity of high-frequency modes, and torsional anharmonicity are important for obtaining quantitative kinetics. We find that at lower temperatures, the formation of the termolecular precursor complexes, which is often neglected, is rate-limiting compared to passage through the tight transition states. Our calculations show that the SO3···H2O + H2O mechanism is more important than the SO3 + (H2O)2 mechanism at 5-50 km altitudes. We find that the rate ratio between SO3···H2O + H2O and SO3 + (H2O)2 is greater than 20 at altitudes between 10 and 35 km, where the concentration of SO3 is very high.
Collapse
Affiliation(s)
- Bo Long
- College of Materials Science and Engineering, Guizhou Minzu University, Guiyang 550025, China
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Yu Xia
- College of Materials Science and Engineering, Guizhou Minzu University, Guiyang 550025, China
| | - Yu-Qiong Zhang
- College of Materials Science and Engineering, Guizhou Minzu University, Guiyang 550025, China
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
11
|
Sun Y, Long B, Truhlar DG. Unimolecular Reactions of E-Glycolaldehyde Oxide and Its Reactions with One and Two Water Molecules. RESEARCH (WASHINGTON, D.C.) 2023; 6:0143. [PMID: 37435010 PMCID: PMC10332847 DOI: 10.34133/research.0143] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/05/2023] [Accepted: 04/20/2023] [Indexed: 07/13/2023]
Abstract
The kinetics of Criegee intermediates are important for atmospheric modeling. However, the quantitative kinetics of Criegee intermediates are still very limited, especially for those with hydroxy groups. Here, we calculate rate constants for the unimolecular reaction of E-glycolaldehyde oxide [E-hydroxyethanal oxide, E-(CH2OH)CHOO], for its reactions with H2O and (H2O)2, and for the reaction of the E-(CH2OH)CHOO…H2O complex with H2O. For the highest level of electronic structure, we use W3X-L//CCSD(T)-F12a/cc-pVDZ-F12 for the unimolecular reaction and the reaction with water and W3X-L//DF-CCSD(T)-F12b/jun-cc-pVDZ for the reaction with 2 water molecules. For the dynamics, we use a dual-level strategy that combines conventional transition state theory with the highest level of electronic structure and multistructural canonical variational transition state theory with small-curvature tunneling with a validated density functional for the electronic structure. This dynamical treatment includes high-frequency anharmonicity, torsional anharmonicity, recrossing effects, and tunneling. We find that the unimolecular reaction of E-(CH2OH)CHOO depends on both temperature and pressure. The calculated results show that E-(CH2OH)CHOO…H2O + H2O is the dominant entrance channel, while previous investigations only considered Criegee intermediates + (H2O)2. In addition, we find that the atmospheric lifetime of E-(CH2OH)CHOO with respect to 2 water molecules is particularly short with a value of 1.71 × 10-6 s at 0 km, which is about 2 orders of magnitude shorter than those usually assumed for Criegee intermediate reactions with water dimer. We also find that the OH group in E-(CH2OH)CHOO enhances its reactivity.
Collapse
Affiliation(s)
- Yan Sun
- Department of Physics, Guizhou University, Guiyang 550025, China
| | - Bo Long
- Department of Physics, Guizhou University, Guiyang 550025, China
- College of Materials Science and Engineering, Guizhou Minzu University, Guiyang 550025, China
| | - Donald G. Truhlar
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, MN 55455-0431, USA
| |
Collapse
|
12
|
Molteni G, Ponti A. Is DFT Accurate Enough to Calculate Regioselectivity? The Case of 1,3-Dipolar Cycloaddition of Azide to Alkynes and Alkenes. Chemphyschem 2023; 24:e202300114. [PMID: 36896728 DOI: 10.1002/cphc.202300114] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2023] [Revised: 03/07/2023] [Accepted: 03/09/2023] [Indexed: 03/11/2023]
Abstract
The importance of regioselectivity in 1,3-dipolar cycloadditions (DCs) makes it surprising that no benchmarking study on this problem has appeared. We investigated whether DFT calculations are an accurate tool to predict the regioselectivity of uncatalyzed thermal azide 1,3-DCs. We considered the reaction between HN3 and 12 dipolarophiles, comprising ethynes HC≡C-R and ethenes H2 C=CH-R (R=F, OH, NH2 , Me, CN, CHO), which cover a broad range of electron demand and conjugation ability. We established benchmark data by the W3X protocol [complete-basis-set-extrapolated CCSD(T)-F12 energy with T-(T) and (Q) corrections and MP2-calculated core/valence and relativistic effects] and showed that core/valence effects and high-order excitations are important for accurate regioselectivity. Regioselectivities calculated using an extensive set of density functional approximations (DFAs) were compared with benchmark data. Range-separated and meta-GGA hybrids gave the best results. Good treatment of self-interaction and electron exchange are the key features for accurate regioselectivity. Dispersion correction slightly improves agreement with W3X results. The best DFAs provide the isomeric TS energy difference with an expected error ≈0.7 mh and errors ≈2 mh can occur. The isomer yield provided by the best DFA has an expected error of ±5 %, though errors up to 20 % are not rare. At present, an accuracy of 1-2 % is unfeasible but it seems that we are not far from achieving this goal.
Collapse
Affiliation(s)
- Giorgio Molteni
- Dipartimento di Chimica, Università degli Studi di Milano, Via C. Golgi 19, 20133, Milano, Italy
| | - Alessandro Ponti
- Istituto di Scienze e Tecnologie Chimiche "Giulio Natta", Consiglio Nazionale delle Ricerche, Via C. Golgi 19, 20133, Milano, Italy
| |
Collapse
|
13
|
Reactions with criegee intermediates are the dominant gas-phase sink for formyl fluoride in the atmosphere. FUNDAMENTAL RESEARCH 2023. [DOI: 10.1016/j.fmre.2023.02.012] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/09/2023] Open
|
14
|
Long B, Xia Y, Truhlar DG. Quantitative Kinetics of HO 2 Reactions with Aldehydes in the Atmosphere: High-Order Dynamic Correlation, Anharmonicity, and Falloff Effects Are All Important. J Am Chem Soc 2022; 144:19910-19920. [PMID: 36264240 DOI: 10.1021/jacs.2c07994] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Kinetics provides the fundamental parameters for elucidating sources and sinks of key atmospheric species and for atmospheric modeling more generally. Obtaining quantitative kinetics in the laboratory for the full range of atmospheric temperatures and pressures is quite difficult. Here, we use computational chemistry to obtain quantitative rate constants for the reactions of HO2 with HCHO, CH3CHO, and CF3CHO. First, we calculate the high-pressure-limit rate constants by using a dual-level strategy that combines conventional transition state theory using a high level of electronic structure wave function theory with canonical variational transition state theory including small-curvature tunneling using density functional theory. The wave-function level is beyond-CCSD(T) for HCHO and CCSD(T)-F12a (Level-A) for XCHO (X = CH3, CF3), and the density functional (Level-B) is specifically validated for these reactions. Then, we calculate the pressure-dependent rate constants by using system-specific quantum RRK theory (SS-QRRK) and also by an energy-grained master equation. The two treatments of the pressure dependence agree well. We find that the Level-A//Level-B method gives good agreement with CCSDTQ(P)/CBS. We also find that anharmonicity is an important factor that increases the rate constants of all three reactions. We find that the HO2 + HCHO reaction has a significant dependence on pressure, but the HO2 + CF3CHO reaction is almost independent of pressure. Our findings show that the HO2 + HCHO reaction makes important contribution to the sink for HCHO, and the HO2 + CF3CHO reaction is the dominant sink for CF3CHO in the atmosphere.
Collapse
Affiliation(s)
- Bo Long
- College of Materials Science and Engineering, Guizhou Minzu University, Guiyang 550025, China
| | - Yu Xia
- College of Materials Science and Engineering, Guizhou Minzu University, Guiyang 550025, China
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
15
|
Zhao YC, Long B, Francisco JS. Quantitative Kinetics of the Reaction between CH 2OO and H 2O 2 in the Atmosphere. J Phys Chem A 2022; 126:6742-6750. [DOI: 10.1021/acs.jpca.2c04408] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Yong-Chao Zhao
- College of Mechanical and Electrical Engineering, Guizhou Minzu University, Guiyang 550025, China
- College of Materials Science and Engineering, Guizhou Minzu University, Guiyang 550025, China
| | - Bo Long
- College of Mechanical and Electrical Engineering, Guizhou Minzu University, Guiyang 550025, China
- College of Materials Science and Engineering, Guizhou Minzu University, Guiyang 550025, China
| | - Joseph S. Francisco
- Department of Earth and Environmental Sciences and Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
16
|
Chan B, Karton A. Assessment of DLPNO-CCSD(T)-F12 and its use for the formulation of the low-cost and reliable L-W1X composite method. J Comput Chem 2022; 43:1394-1402. [PMID: 35709311 DOI: 10.1002/jcc.26892] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2022] [Revised: 03/28/2022] [Accepted: 04/28/2022] [Indexed: 02/03/2023]
Abstract
In the present study, we have investigated the performance of RIJCOSX DLPNO-CCSD(T)-F12 methods for a wide range of systems. Calculations with a high-accuracy option ["DefGrid3 RIJCOSX DLPNO-CCSD(T1 )-F12"] extrapolated to the complete-basis-set limit using the maug-cc-pV[D+d,T+d]Z basis sets provides fairly good agreements with the canonical CCSD(T)/CBS reference for a diverse set of thermochemical and kinetic properties [with mean absolute deviations (MADs) of ~1-2 kJ mol-1 except for atomization energies]. On the other hand, the low-cost "RIJCOSX DLPNO-CCSD(T)-F12D" option leads to substantial deviations for certain properties, notably atomization energies (MADs of up to tens of kJ mol-1 ). With the high-accuracy CBS approach, we have formulated the L-W1X method, which further includes a low-cost core-valence plus scalar-relativistic term. It shows generally good accuracy. For improved accuracies in specific cases, we advise replacing maug-cc-pV(n+d)Z with jun-cc-pV(n+d)Z for the calculation of electron affinities, and using well-constructed isodesmic-type reactions to obtain atomization energies. For medium-sized systems, DefGrid3 RIJCOSX DLPNO-CCSD(T1 )-F12 calculations are several times faster than the corresponding canonical computation; the use of the local approximations (RIJCOSX and DLPNO) leads to a better scaling than that for the canonical calculation (from ~6-7 down to ~2-4 for our test systems). Thus, the DefGrid3 RIJCOSX DLPNO-CCSD(T1 )-F12 method, and the L-W1X protocol that based on it, represent a useful means for obtaining accurate thermochemical quantities for larger systems.
Collapse
Affiliation(s)
- Bun Chan
- Graduate School of Engineering, Nagasaki University, Nagasaki, Japan
| | - Amir Karton
- School of Molecular Sciences, University of Western Australia, Perth, WA, Australia
| |
Collapse
|
17
|
Chan B. High-Level Quantum Chemistry Reference Heats of Formation for a Large Set of C, H, N, and O Species in the NIST Chemistry Webbook and the Identification and Validation of Reliable Protocols for Their Rapid Computation. J Phys Chem A 2022; 126:4981-4990. [PMID: 35878062 DOI: 10.1021/acs.jpca.2c03846] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
A recent study has examined the accuracy of NIST heats of formation for a set of C, H, and O-containing species with a proposed low-cost quantum chemistry approach. In the present study, we have used high-level methods such as W1X-2 to obtain these data more rigorously, which we have then used to assess the NIST and the previously computed values. We find that many of these NIST data that are as suggested to be unreliable by the previous study are indeed inconsistent with our high-level reference values. However, we also find substantial deviations for the previously computed values from our benchmark. Thus, we have assessed the performance of alternative low-cost methods. In our assessment, we have additionally examined C, H, N, and O-containing species for which heats of formation are available from the NIST database. We find the ωB97M-V/ma-def2-TZVP, DSD-PBEP86/ma-def2-TZVP, and CCSD(T)-F12b/aug'-cc-pVDZ methods to be adequate for obtaining heats of formation with the atomization approach, once their atomic energies are optimized with our benchmark. Notably, the low-cost ωB97M-V method yields values that agree to be within 10 kJ mol-1 for more than 90% of the (∼1500) species. A higher 20 kJ mol-1 threshold captures 98% of the data. The outlier species typically contain many electron-withdrawing (nitro) groups. In these cases, the use of isodesmic-type reactions rather than the atomization approach is more reliable. Our assessment has also identified significant outliers from the NIST database, for which experimental re-determination of the heats of formation would be desirable.
Collapse
Affiliation(s)
- Bun Chan
- Graduate School of Engineering, Nagasaki University, Bunkyo 1-14, Nagasaki 852-8521, Japan
| |
Collapse
|
18
|
Bakowies D. ATOMIC-2 Protocol for Thermochemistry. J Chem Theory Comput 2022; 18:4142-4163. [PMID: 35658473 DOI: 10.1021/acs.jctc.1c01272] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
ATOMIC is a midlevel thermochemistry protocol that uses Pople's concept of bond separation reactions (BSRs) as a theoretical framework to reduce computational demands in the evaluation of atomization energies and enthalpies of formation. Various composite models are available that approximate bond separation energies at the complete-basis-set limit of all-electron CCSD(T), each balancing computational cost with achievable accuracy. Evaluated energies are then combined with very high-level, precomputed atomization energies of all auxiliary molecules appearing in the BSR to obtain the atomization energy of the molecule under study. ATOMIC-2 is a new version of the protocol that retains the overall concept and all previously defined composite models but improves on ATOMIC-1 in various other ways: Geometry optimization and zero-point-energy evaluation are performed at the density functional level (PBE0-D3/6-311G(d)), which shows significant computational savings and better accuracy than the previously employed RI-MP2/cc-pVTZ. The BSR framework is improved, using more accurate complete-basis-set (CBS) extrapolations toward the Full CI limit for the atomization energies of all auxiliary molecules. Finally, and most importantly, an error and uncertainty model termed ATOMIC-2um is added that estimates average bias and uncertainty for each of the atomization energy contributions that arise from the simplified treatment of some contributions to bond separation energies (CCSD(T)) and the neglect of others (such as higher order, scalar relativistic, or diagonal Born-Oppenheimer corrections) or from residual error in the energies of auxiliary molecules. Large and diverse benchmarks including up to 1179 molecules are used to evaluate necessary reference data and to correlate the observed error for each of the contributions with appropriate proxies that are available without additional quantum-chemical calculations for a particular molecule and represent its size and type. The implementation of ATOMIC-2 considers neutral, closed-shell molecules containing H, C, N, O, and F atoms; compared to ATOMIC-1, the framework has been extended to cover a few challenging but rare bond topologies. In comparison to highly accurate reference data for 184 molecules taken from the ATcT database (V. 1.122r), regular ATOMIC-2 shows noticeable underbinding, but the bias-corrected protocol ATOMIC-2um is found to be more accurate than either ATOMIC-1 or standard Gaussian-4 theory, and the uncertainty model is consistent with statistics of actually observed errors. Problems arising from ambiguous or challenging Lewis-valence structures defining BSRs are discussed, and computational efficiency is demonstrated. Computer code is made available to perform ATOMIC-2um analyses.
Collapse
Affiliation(s)
- Dirk Bakowies
- Institute of Physical Chemistry, Department of Chemistry, University of Basel, Klingelbergstraße 80, CH 4056 Basel, Switzerland
| |
Collapse
|
19
|
Wang PB, Truhlar DG, Xia Y, Long B. Temperature-dependent kinetics of the atmospheric reaction between CH 2OO and acetone. Phys Chem Chem Phys 2022; 24:13066-13073. [PMID: 35583864 DOI: 10.1039/d2cp01118b] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Criegee intermediates are important oxidants produced in the ozonolysis of alkenes in the atmosphere. Quantitative kinetics of the reactions of Criegee intermediates are required for atmospheric modeling. However, the experimental studies do not cover the full relevant range of temperature and pressure. Here we report the quantitative kinetics of CH2OO + CH3C(O)CH3 by using our recently developed dual strategy that combines coupled cluster theory with high excitation levels for conventional transition state theory and well validated levels of density functional theory for direct dynamics calculations using canonical variational transition theory including tunneling. We find that the W3X-L//DF-CCSD(T)-F12b/jun-cc-pVDZ electronic structure method can be used to obtain quantitative kinetics of the CH2OO + CH3C(O)CH3 reaction. Whereas previous investigations considered a one-step mechanistic pathway, we find that the CH2OO + CH3C(O)CH3 reaction occurs in a stepwise manner. This has implications for the modeling of Criegee-intermediate reactions with other ketones and with aldehydes. In the kinetics calculations, we show that recrossing effects of the conventional transition state are negligible for determining the rate constant of CH2OO + CH3C(O)CH3. The present findings reveal that the rate ratio between CH2OO + CH3C(O)CH3 and OH + CH3C(O)CH3 has a significant negative dependence on temperature such that the CH2OO + CH3C(O)CH3 reaction can contribute as a significant sink for atmospheric CH3C(O)CH3 at low temperature. The present findings should have broad implications in understanding the reactions of Criegee intermediates with carbonyl compounds and ketones in the atmosphere.
Collapse
Affiliation(s)
- Peng-Biao Wang
- Department of Physics, Guizhou University, Guiyang, 550025, China.
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA
| | - Yu Xia
- College of Materials Science and Engineering, Guizhou Minzu University, Guiyang 550025, China
| | - Bo Long
- Department of Physics, Guizhou University, Guiyang, 550025, China. .,College of Materials Science and Engineering, Guizhou Minzu University, Guiyang 550025, China
| |
Collapse
|
20
|
Long B, Xia Y, Bao JL, Carmona-García J, Gómez Martín JC, Plane JMC, Saiz-Lopez A, Roca-Sanjuán D, Francisco JS. Reaction of SO 3 with HONO 2 and Implications for Sulfur Partitioning in the Atmosphere. J Am Chem Soc 2022; 144:9172-9177. [PMID: 35576167 DOI: 10.1021/jacs.2c03499] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Sulfur trioxide is a critical intermediate for the sulfur cycle and the formation of sulfuric acid in the atmosphere. The traditional view is that sulfur trioxide is removed by water vapor in the troposphere. However, the concentration of water vapor decreases significantly with increasing altitude, leading to longer atmospheric lifetimes of sulfur trioxide. Here, we utilize a dual-level strategy that combines transition state theory calculated at the W2X//DF-CCSD(T)-F12b/jun'-cc-pVDZ level, with variational transition state theory with small-curvature tunneling from direct dynamics calculations at the M08-HX/MG3S level. We also report the pressure-dependent rate constants calculated using the system-specific quantum Rice-Ramsperger-Kassel (SS-QRRK) theory. The present findings show that falloff effects in the SO3 + HONO2 reaction are pronounced below 1 bar. The SO3 + HONO2 reaction can be a potential removal reaction for SO3 in the stratosphere and for HONO2 in the troposphere, because the reaction can potentially compete well with the SO3 + 2H2O reaction between 25 and 35 km, as well as the OH + HONO2 reaction. The present findings also suggest an unexpected new product from the SO3 + HONO2 reaction, which, although very short-lived, would have broad implications for understanding the partitioning of sulfur in the stratosphere and the potential for the SO3 reaction with organic acids to generate organosulfates without the need for heterogeneous chemistry.
Collapse
Affiliation(s)
- Bo Long
- College of Materials Science and Engineering, Guizhou Minzu University, 550025 Guiyang, China
| | - Yu Xia
- College of Materials Science and Engineering, Guizhou Minzu University, 550025 Guiyang, China
| | - Junwei Lucas Bao
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts 02467, United States
| | - Javier Carmona-García
- Institut de Ciència Molecular, Universitat de València, València 46071, Spain.,Department of Atmospheric Chemistry and Climate, Institute of Physical Chemistry Rocasolano, CSIC, Madrid 28006, Spain
| | | | - John M C Plane
- School of Chemistry, University of Leeds, LS2 9JT Leeds, U.K
| | - Alfonso Saiz-Lopez
- Department of Atmospheric Chemistry and Climate, Institute of Physical Chemistry Rocasolano, CSIC, Madrid 28006, Spain
| | - Daniel Roca-Sanjuán
- Institut de Ciència Molecular, Universitat de València, València 46071, Spain
| | - Joseph S Francisco
- Department of Earth and Environmental Sciences and Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
21
|
Xia Y, Long B, Lin S, Teng C, Bao JL, Truhlar DG. Large Pressure Effects Caused by Internal Rotation in the s-cis-syn-Acrolein Stabilized Criegee Intermediate at Tropospheric Temperature and Pressure. J Am Chem Soc 2022; 144:4828-4838. [PMID: 35262353 DOI: 10.1021/jacs.1c12324] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Criegee intermediates are important atmospheric oxidants, and quantitative kinetics for stabilized Criegee intermediates are key parameters for atmospheric modeling but are still limited. Here we report barriers and rate constants for unimolecular reactions of s-cis-syn-acrolein oxide (scsAO), in which the vinyl group makes it a prototype for Criegee intermediates produced in the ozonolysis of isoprene. We find that the MN15-L and M06-2X density functionals have CCSD(T)/CBS accuracy for the unimolecular cyclization and stereoisomerization of scsAO. We calculated high-pressure-limit rate constants by the dual-level strategy that combines (a) high-level wave function-based conventional transition-state theory (which includes coupled-cluster calculations with quasiperturbative inclusion of quadruple excitations because of the strongly multiconfigurational character of the electronic wave function) and (b) canonical variational transition-state theory with small-curvature tunneling based on a validated density functional. We calculated pressure-dependent rate constants both by system-specific quantum Rice-Ramsperger-Kassel theory and by solving the master equation. We report rate constants for unimolecular reactions of scsAO over the full range of atmospheric temperature and pressure. We found that the unimolecular reaction rates of this larger-than-previously studied Criegee intermediate depend significantly on pressure. Particularly, we found that falloff effects decrease the effective unimolecular cyclization rate constant of scsAO by about a factor of 3, but the unimolecular reaction is still the dominant atmospheric sink for scsAO at low altitudes. The large falloff caused by the inclusion of the stereoisomerization channel in the master equation calculations has broad implications for mechanistic analysis of reactions with competitive internal rotations that can produce stable rotamers.
Collapse
Affiliation(s)
- Yu Xia
- College of Mechanical and Electrical Engineering, Guizhou Minzu University, Guiyang 550025, China.,College of Materials Science and Engineering, Guizhou Minzu University, Guiyang 550025, China
| | - Bo Long
- College of Mechanical and Electrical Engineering, Guizhou Minzu University, Guiyang 550025, China.,College of Materials Science and Engineering, Guizhou Minzu University, Guiyang 550025, China
| | - Shiru Lin
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts 02467, United States
| | - Chong Teng
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts 02467, United States
| | - Junwei Lucas Bao
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts 02467, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
22
|
Simmie JM. C 2H 5NO Isomers: From Acetamide to 1,2-Oxazetidine and Beyond. J Phys Chem A 2022; 126:924-939. [PMID: 35113546 PMCID: PMC8859852 DOI: 10.1021/acs.jpca.1c09984] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
This work documents
the properties of a number of isomers of molecular
formula C2H5NO from the most stable, acetamide,
through 1,2-oxazetidine and including even higher energy species largely
of a dipolar nature. Only two of the isomers have been detected in
emissions from the interstellar medium (ISM); possible further candidates
are identified, and the likelihood of their being detectable is considered.
In general, hardly any of these compounds have been discussed in the
existing chemical literature, so this work represents an important
contribution extending the canon of chemical bonding which can contribute
to machine learning, providing a more exacting test of AI applications.
The presence in the ISM of acetamide, CH3C(O)NH2, is the subject of current debate with no clear and obvious paths
to its formation; it is shown that a 1,3-[H]-transfer from (E,Z)-ethanimidic acid, CH3C(OH)=NH, is
feasible in spite of an energy barrier of 130 kJ mol–1. It is speculated that imidic acid can itself be formed from abundant
precursors, H2O and CH3C≡N, in an acid-induced,
water addition, autocatalytic reaction on water–ice grains.
H3CC≡NH3CC≡NH+ +
H2OH3CC(O+H2)=NHH3CC(OH)=NH
+ H3O+
Collapse
Affiliation(s)
- John M Simmie
- School of Chemistry, National University of Ireland, Galway H91 TK33, Ireland
| |
Collapse
|
23
|
Ji P, Luo YR, Xue XS, Cheng JP. Efficient estimation of bond dissociation energies of organic compounds. ADVANCES IN PHYSICAL ORGANIC CHEMISTRY 2022. [DOI: 10.1016/bs.apoc.2022.10.001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
24
|
Zhang Y, Cheng Y, Zhang T, Wang R, Ji J, Xia Y, Makroni L, Wang Z, M B. A computational study of the HO2 + SO3 → HOSO2 + 3O2 reaction catalyzed by water monomer, water dimer and small clusters of sulfuric acid: kinetics and atmospheric implications. Phys Chem Chem Phys 2022; 24:18205-18216. [DOI: 10.1039/d1cp03318b] [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/2022]
Abstract
Herein, the reaction mechanisms and kinetics for the HO2 + SO3 → HOSO2 + 3O2 reaction catalyzed by water monomer, water dimer and small clusters of sulfuric acid have been...
Collapse
|
25
|
Klippenstein SJ. Spiers Memorial Lecture: theory of unimolecular reactions. Faraday Discuss 2022; 238:11-67. [DOI: 10.1039/d2fd00125j] [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/2022]
Abstract
One hundred years ago, at an earlier Faraday Discussion meeting, Lindemann presented a mechanism that provides the foundation for contemplating the pressure dependence of unimolecular reactions. Since that time, our...
Collapse
|
26
|
Chan B. Accurate Thermochemistry for Main-Group Elements up to Xenon with the W n-P34 Series of Composite Methods. J Chem Theory Comput 2021; 17:5704-5714. [PMID: 34410730 DOI: 10.1021/acs.jctc.1c00598] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
In the present study, we introduce the accurate Wn-P34 quantum chemistry composite methods with applicability to heavy p-block elements up to xenon. For a set of thermochemical properties for prototypical third- and fourth-row species and for a diverse set of small light-main-group species, they show accuracies of ∼3 kJ mol-1 or better. Overall, the Wn-P34 methods are comparable in accuracy to Wn, with a widened applicability to heavier elements. We have used Wn-P34 to compile the P34 set of accurate thermochemical values for heavy p-block species, and we have applied this set to assess a wide range of lower-cost methods. The results of our assessment show that the G4(MP2)-XK composite method provides adequate treatments for these species, but several widely used double-hybrid density functional theory (DH-DFT) methods show uncharacteristically large deviations. In contrast, we find it presently surprising that some pure and hybrid DFT methods such as TPSS and SCANh perform quite well. We hope that our findings and new tools would facilitate the application of computational chemistry for heavy elements, of which the properties are yet to be broadly explored.
Collapse
Affiliation(s)
- Bun Chan
- Graduate School of Engineering, Nagasaki University, Bunkyo 1-14, Nagasaki 852-8521, Japan
| |
Collapse
|
27
|
Chan B, Karton A. Polycyclic aromatic hydrocarbons: from small molecules through nano-sized species towards bulk graphene. Phys Chem Chem Phys 2021; 23:17713-17723. [PMID: 34378574 DOI: 10.1039/d1cp01659h] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We have examined the use of systematic bond-separation reactions and purposely constructed chemistry-preserving isodesmic reactions for the thermochemical calculation of aromatic hydrocarbon species. The bond-separation approach yields somewhat disappointing accuracy even when the reaction energies are obtained with generally robust composite and double-hybrid (DH) density functional theory (DFT) methods. In contrast, for the purposely constructed reactions, we find a dramatic improvement in the accuracy for energies calculated with all methods examined. Notably, for medium-sized aromatic hydrocarbons, we find that an effective approach for formulating a well-balanced reaction is to split the target species into two halves with an aromatic overlapping region. Overall, the G4(MP2)-XK, MPW2PLYP, MN15, PBE, and DC-DFTB3 methods are reasonable within their respective classes of methods for the calculation of bond-separation as well as chemistry-preserving isodesmic reactions. We have further computed per-carbon atomization energy (AE) for a series of D6h benzene-type molecules, and thus obtained a formula for extrapolation to the graphene limit [AEn = 711.5 × (1 - 1/n0.640) kJ mol-1, where n = number of carbons]. It suggests that nano-graphene with a length larger than 10 nm would resemble properties of bulk graphene, and conversely, downsizing a nano-graphene beyond this point may lead to considerably altered properties from the bulk.
Collapse
Affiliation(s)
- Bun Chan
- Graduate School of Engineering, Nagasaki University, Bunkyo 1-14, Nagasaki-Shi, Nagasaki 852-8521, Japan.
| | | |
Collapse
|
28
|
Long B, Wang Y, Xia Y, He X, Bao JL, Truhlar DG. Atmospheric Kinetics: Bimolecular Reactions of Carbonyl Oxide by a Triple-Level Strategy. J Am Chem Soc 2021; 143:8402-8413. [PMID: 34029069 DOI: 10.1021/jacs.1c02029] [Citation(s) in RCA: 26] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Criegee intermediates in the atmosphere serve as oxidizing agents to initiate aerosol formation, which are particularly important for atmospheric modeling, and understanding their kinetics is one of the current outstanding challenges in climate change modeling. Because experimental kinetics are still limited, we must rely on theory for the complete picture, but obtaining absolute rates from theory is a formidable task. Here, we report the bimolecular reaction kinetics of carbonyl oxide with ammonia, hydrogen sulfide, formaldehyde, and water dimer by designing a triple-level strategy that combines (i) benchmark results close to the complete-basis limit of coupled-cluster theory with the single, double, triple, and quadruple excitations (CCSDTQ/CBS), (ii) a new hybrid meta density functional (M06CR) specifically optimized for reactions of Criegee intermediates, and (iii) variational transition-state theory with both variable rection coordinates and optimized reaction paths, with multidimensional tunneling, and with pressure effects. For (i) we have found that quadruple excitations are required to obtain quantitative reaction barriers, and we designed new composite methods and strategies to reach CCSDTQ/CBS accuracy. The present findings show that (i) the CH2OO + HCHO reaction can make an important contribution to the sink of HCHO under wide atmospheric conditions in the gas phase and that (ii) CH2OO + (H2O)2 dominates over the CH2OO + H2O reaction below 10 km.
Collapse
Affiliation(s)
- Bo Long
- College of Materials Science and Engineering, Guizhou Minzu University, Guiyang 550025, China.,Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Ying Wang
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China.,The National and Local Joint Engineering Laboratory of Animal Peptide Drug Development, College of Life Sciences, Hunan Normal University, Changsha 410006, China
| | - Yu Xia
- College of Materials Science and Engineering, Guizhou Minzu University, Guiyang 550025, China
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| | - Junwei Lucas Bao
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts 02467, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
29
|
Li Y, Javoy S, Mevel R, Xu X. A chemically consistent rate constant for the reaction of nitrogen dioxide with the oxygen atom. Phys Chem Chem Phys 2021; 23:585-596. [PMID: 33331363 DOI: 10.1039/d0cp05131d] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
In the present study, a chemically consistent rate constant for the reaction between nitrogen dioxide and the oxygen atom has been obtained by combining low-temperature experimental data from the literature and new high-temperature quantum chemical calculations. The expression for our rate constant is kNO2+O=NO+O2= 2.589 × 1015T-1.035 exp(-226/RT) + 4.242 × 1016T-0.861exp(-50 917/RT) cm3 mol-1 s-1, where R = 8.314 J mol-1 K-1, and is valid over the temperature range T = 221 to 3000 K. The effect of the inclusion of the new rate constant on the prediction of three detailed reaction models from the literature has been studied using (i) new experimental oxygen atom profiles obtained in a shock tube during nitrogen dioxide pyrolysis, and (ii) published shock tube and jet-stirred reactor data for H2-NOx mixtures with and without dioxygen. The impact of the new rate constant on the sensitivity coefficients and reaction pathways has also been analyzed under some conditions. Overall, the predictive capability of the reaction models were improved. The present study suggests that our chemically consistent rate constant should be included in detailed reaction models for combustion applications.
Collapse
Affiliation(s)
- Yan Li
- Center for Combustion Energy, Department of Energy and Power Engineering, and Key Laboratory for Thermal Science and Power Engineering of Ministry of Education, Tsinghua University, Beijing 100084, China.
| | | | | | | |
Collapse
|
30
|
Semidalas E, Martin JML. Canonical and DLPNO-Based Composite Wavefunction Methods Parametrized against Large and Chemically Diverse Training Sets. 2: Correlation-Consistent Basis Sets, Core-Valence Correlation, and F12 Alternatives. J Chem Theory Comput 2020; 16:7507-7524. [PMID: 33200931 PMCID: PMC7735707 DOI: 10.1021/acs.jctc.0c01106] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
A hierarchy
of wavefunction composite methods (cWFT), based on
G4-type cWFT methods available for elements H through Rn, was recently
reported by the present authors [J. Chem. Theor. Comput.2020, 16, 4238]. We extend this hierarchy
by considering the inner-shell correlation energy in the second-order
Møller–Plesset correction and replacing the Weigend–Ahlrichs
def2-mZVPP(D) basis sets used with complete basis
set extrapolation from augmented correlation-consistent core–valence
triple-ζ, aug-cc-pwCVTZ(-PP), and quadruple-ζ, aug-cc-pwCVQZ(-PP),
basis sets, thus creating cc-G4-type methods. For the large and chemically
diverse GMTKN55 benchmark suite, they represent a substantial further
improvement and bring WTMAD2 (weighted mean absolute deviation) down
below 1 kcal/mol. Intriguingly, the lion’s share of the improvement
comes from better capture of valence correlation; the inclusion of
core–valence correlation is almost an order of magnitude less
important. These robust correlation-consistent cWFT methods approach
the CCSD(T) complete basis limit with just one or a few fitted parameters.
Particularly, the DLPNO variants such as cc-G4-T-DLPNO are applicable
to fairly large molecules at a modest computational cost, as is (for
a reduced range of elements) a different variant using MP2-F12/cc-pVTZ-F12
for the MP2 component.
Collapse
Affiliation(s)
- Emmanouil Semidalas
- Department of Organic Chemistry, Weizmann Institute of Science, 7610001 Rehovot, Israel
| | - Jan M L Martin
- Department of Organic Chemistry, Weizmann Institute of Science, 7610001 Rehovot, Israel
| |
Collapse
|
31
|
Zhang T, Wen M, Zhang Y, Chen X, Qiao Z, Su Y, Lily M, Wang Z. Sulfuric acid catalyzed HCl + HO → Cl + H2O reaction in troposphere: A quantum chemical investigation. COMPUT THEOR CHEM 2020. [DOI: 10.1016/j.comptc.2020.112936] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
32
|
Chan B. Fullerene Thermochemical Stability: Accurate Heats of Formation for Small Fullerenes, the Importance of Structural Deformation on Reactivity, and the Special Stability of C 60. J Phys Chem A 2020; 124:6688-6698. [PMID: 32786665 DOI: 10.1021/acs.jpca.0c04732] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
We have used quantum chemistry computations, in conjunction with isodesmic-type reactions, to obtain accurate heats of formation (HoFs) for the small fullerenes C20 (2358.2 ± 8.0 kJ mol-1), C24 (2566.2 ± 7.6), and the lowest-energy isomers of C32 (2461.1 ± 15.4), C42 (2629.0 ± 20.5), and C54 (2686.2 ± 25.3). As part of this endeavor, we have also obtained accurate HoFs for several medium-sized molecules, namely 216.6 ± 1.4 for fulvene, 375.5 ± 1.5 for pentalene, 670.8 ± 2.9 for acepentalene, and 262.7 ± 2.5 for acenaphthylene. We combine the energies of the small fullerenes and previously obtained energies for larger fullerenes (from C60 to C6000) into a full picture of fullerene thermochemical stability. In general, the per-carbon energies can be reasonably approximated by the "R+D" model that we have previously developed [Chan et al. J. Chem. Theory Comput. 2019, 15, 1255-1264], which takes into account Resonance and structural Deformation factors. In a case study on C54, we find that most of the high-deformation-energy atoms correspond to the sites of the C-Cl bond in the experimentally captured C54Cl8. In another case study, we find that C60 has the lowest value for the maximum local-deformation energy when compared with similar-sized fullerenes, which is consistent with its "special stability". These results are indicative of structural deformation playing an important role in the reactivity of fullerenes.
Collapse
Affiliation(s)
- Bun Chan
- Graduate School of Engineering, Nagasaki University, Bunkyo 1-14, Nagasaki-shi, Nagasaki 852-8521, Japan
| |
Collapse
|
33
|
Zhang T, Bi X, Wen M, Liu S, Chai G, Zeng Z, Wang R, Wang W, Long B. The HO 4H → O 3 + H 2O reaction catalysed by acidic, neutral and basic catalysts in the troposphere. Mol Phys 2020. [DOI: 10.1080/00268976.2019.1673912] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Affiliation(s)
- Tianlei Zhang
- School of Chemical & Environment Science, Shaanxi Key Laboratory of Catalysis, Institute of Theoretical and Computational Chemistry, Shaanxi University of Technology, Hanzhong, People’s Republic of China
- Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, Fudan University, Shanghai, People’s Republic of China
| | - Xiujuan Bi
- School of Chemical & Environment Science, Shaanxi Key Laboratory of Catalysis, Institute of Theoretical and Computational Chemistry, Shaanxi University of Technology, Hanzhong, People’s Republic of China
| | - Mingjie Wen
- School of Chemical & Environment Science, Shaanxi Key Laboratory of Catalysis, Institute of Theoretical and Computational Chemistry, Shaanxi University of Technology, Hanzhong, People’s Republic of China
- Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, Fudan University, Shanghai, People’s Republic of China
| | - Shuai Liu
- School of Chemical & Environment Science, Shaanxi Key Laboratory of Catalysis, Institute of Theoretical and Computational Chemistry, Shaanxi University of Technology, Hanzhong, People’s Republic of China
| | - Guang Chai
- School of Chemical & Environment Science, Shaanxi Key Laboratory of Catalysis, Institute of Theoretical and Computational Chemistry, Shaanxi University of Technology, Hanzhong, People’s Republic of China
| | - Zhaopeng Zeng
- School of Chemical & Environment Science, Shaanxi Key Laboratory of Catalysis, Institute of Theoretical and Computational Chemistry, Shaanxi University of Technology, Hanzhong, People’s Republic of China
| | - Rui Wang
- School of Chemical & Environment Science, Shaanxi Key Laboratory of Catalysis, Institute of Theoretical and Computational Chemistry, Shaanxi University of Technology, Hanzhong, People’s Republic of China
| | - Wenliang Wang
- Key Laboratory for Macromolecular Science of Shaanxi Province, School of Chemistry & Chemical Engineering, Shaanxi Normal University, Xi’an, People’s Republic of China
| | - Bo Long
- School of Materials Science and Engineering, Guizhou Minzu University, Guiyang, People’s Republic of China
| |
Collapse
|
34
|
Semidalas E, Martin JML. Canonical and DLPNO-Based G4(MP2)XK-Inspired Composite Wave Function Methods Parametrized against Large and Chemically Diverse Training Sets: Are They More Accurate and/or Robust than Double-Hybrid DFT? J Chem Theory Comput 2020; 16:4238-4255. [PMID: 32456427 PMCID: PMC7366511 DOI: 10.1021/acs.jctc.0c00189] [Citation(s) in RCA: 26] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
The
large and chemically diverse GMTKN55 benchmark was used as
a training set for parametrizing composite wave function thermochemistry
protocols akin to G4(MP2)XK theory (Chan, B.; Karton, A.; Raghavachari,
K. J. Chem. Theory Comput. 2019, 15, 4478–4484). On account of their availability
for elements H through Rn, Karlsruhe def2 basis sets were employed.
Even after reparametrization, the GMTKN55 WTMAD2 (weighted mean absolute
deviation, type 2) for G4(MP2)-XK is actually inferior to that of
the best rung-4 DFT functional, ωB97M-V. By increasing the basis
set for the MP2 part to def2-QZVPPD, we were able to substantially
improve performance at modest cost (if an RI-MP2 approximation is
made), with WTMAD2 for this G4(MP2)-XK-D method now comparable to
the better rung-5 functionals (albeit at greater cost). A three-tier
approach with a scaled MP3/def2-TZVPP intermediate step, however,
leads to a G4(MP3)-D method that is markedly superior to even the
best double hybrids ωB97M(2) and revDSD-PBEP86-D4. Evaluating
the CCSD(T) component with a triple-ζ, rather than split-valence,
basis set yields only a modest further improvement that is incommensurate
with the drastic increase in computational cost. G4(MP3)-D and G4(MP2)-XK-D
have about 40% better WTMAD2, at similar or lower computational cost,
than their counterparts G4 and G4(MP2), respectively: detailed comparison
reveals that the difference lies in larger molecules due to basis
set incompleteness error. An E2/{T,Q} extrapolation and a CCSD(T)/def2-TZVP
step provided the G4-T method of high accuracy and with just three
fitted parameters. Using KS orbitals in MP2 leads to the G4(MP3|KS)-D
method, which entirely eliminates the CCSD(T) step and has no steps
costlier than scaled MP3; this shows a path forward to further improvements
in double-hybrid density functional methods. None of our final selections
require an empirical HLC correction; this cuts the number of empirical
parameters in half and avoids discontinuities on potential energy
surfaces. G4-T-DLPNO, a variant in which post-MP2 corrections are
evaluated at the DLPNO-CCSD(T) level, achieves nearly the accuracy
of G4-T but is applicable to much larger systems.
Collapse
Affiliation(s)
- Emmanouil Semidalas
- Department of Organic Chemistry, Weizmann Institute of Science, 7610001 Reḩovot, Israel
| | - Jan M L Martin
- Department of Organic Chemistry, Weizmann Institute of Science, 7610001 Reḩovot, Israel
| |
Collapse
|
35
|
Zhang T, Zhai K, Zhang Y, Geng L, Geng Z, Zhou M, Lu Y, Shao X, Lily M. Effect of water and ammonia on the HO + NH3 → NH2 + H2O reaction in troposphere: Competition between single and double hydrogen atom transfer pathways. COMPUT THEOR CHEM 2020. [DOI: 10.1016/j.comptc.2020.112747] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
36
|
Zhang W, Kong X, Liu S, Zhao Y. Multi‐coefficients correlation methods. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2020. [DOI: 10.1002/wcms.1474] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Wenna Zhang
- The Institute of Technological Sciences, Wuhan University Wuhan Hubei People's Republic of China
| | - Xirui Kong
- The Institute of Technological Sciences, Wuhan University Wuhan Hubei People's Republic of China
| | - Sheng Liu
- The Institute of Technological Sciences, Wuhan University Wuhan Hubei People's Republic of China
| | - Yan Zhao
- The Institute of Technological Sciences, Wuhan University Wuhan Hubei People's Republic of China
- State Key Laboratory of Advanced Technology for Materials Synthesis and Processing Wuhan University of Technology Wuhan Hubei People's Republic of China
| |
Collapse
|
37
|
Association of Cl with C 2H 2 by unified variable-reaction-coordinate and reaction-path variational transition-state theory. Proc Natl Acad Sci U S A 2020; 117:5610-5616. [PMID: 32123079 DOI: 10.1073/pnas.1920018117] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Barrierless unimolecular association reactions are prominent in atmospheric and combustion mechanisms but are challenging for both experiment and kinetics theory. A key datum for understanding the pressure dependence of association and dissociation reactions is the high-pressure limit, but this is often available experimentally only by extrapolation. Here we calculate the high-pressure limit for the addition of a chlorine atom to acetylene molecule (Cl + C2H2→C2H2Cl). This reaction has outer and inner transition states in series; the outer transition state is barrierless, and it is necessary to use different theoretical frameworks to treat the two kinds of transition state. Here we study the reaction in the high-pressure limit using multifaceted variable-reaction-coordinate variational transition-state theory (VRC-VTST) at the outer transition state and reaction-path variational transition state theory (RP-VTST) at the inner turning point; then we combine the results with the canonical unified statistical (CUS) theory. The calculations are based on a density functional validated against the W3X-L method, which is based on coupled cluster theory with single, double, and triple excitations and a quasiperturbative treatment of connected quadruple excitations [CCSDT(Q)], and the computed rate constants are in good agreement with some of the experimental results. The chlorovinyl (C2H2Cl) adduct has two isomers that are equilibrium structures of a double-well C≡C-H bending potential. Two procedures are used to calculate the vibrational partition function of chlorovinyl; one treats the two isomers separately and the other solves the anharmonic energy levels of the double well. We use these results to calculate the standard-state free energy and equilibrium constant of the reaction.
Collapse
|
38
|
Tan XF, Zhang L, Long B. New mechanistic pathways for the formation of organosulfates catalyzed by ammonia and carbinolamine formation catalyzed by sulfuric acid in the atmosphere. Phys Chem Chem Phys 2020; 22:8800-8807. [DOI: 10.1039/c9cp06297a] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Sulfuric acid exerts a remarkable catalytic role in the H2SO4 + HCHO + NH3 reaction that leads to the formation of carbinolamine.
Collapse
Affiliation(s)
- Xing-Feng Tan
- School of Mechatronics Engineering
- Guizhou Minzu University
- Guiyang
- China
| | - Lin Zhang
- Department of Physics
- Guizhou University
- Guiyang
- China
| | - Bo Long
- School of Materials Science and Engineering, Guizhou Minzu University
- Guiyang
- China
| |
Collapse
|
39
|
Zhang T, Wen M, Zeng Z, Lu Y, Wang Y, Wang W, Shao X, Wang Z, Makroni L. Effect of NH 3 and HCOOH on the H 2O 2 + HO → HO 2 + H 2O reaction in the troposphere: competition between the one-step and stepwise mechanisms. RSC Adv 2020; 10:9093-9102. [PMID: 35496523 PMCID: PMC9050117 DOI: 10.1039/d0ra00024h] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/02/2020] [Accepted: 02/13/2020] [Indexed: 11/21/2022] Open
Abstract
The H2O2 + HO → HO2 + H2O reaction is an important reservoir for both radicals of HO and HO2 catalyzing the destruction of O3. Here, this reaction assisted by NH3 and HCOOH catalysts was explored using the CCSD(T)-F12a/cc-pVDZ-F12//M06-2X/aug-cc-pVTZ method and canonical variational transition state theory with small curvature tunneling. Two possible sets of mechanisms, (i) one-step routes and (ii) stepwise processes, are possible. Our results show that in the presence of both NH3 and HCOOH catalysts under relevant atmospheric temperature, mechanism (i) is favored both energetically and kinetically than the corresponding mechanism (ii). At 298 K, the relative rate for mechanism (i) in the presence of NH3 (10, 2900 ppbv) and HCOOH (10 ppbv) is respectively 3–5 and 2–4 orders of magnitude lower than that of the water-catalyzed reaction. This is due to a comparatively lower concentration of NH3 and HCOOH than H2O which indicates the positive water effect under atmospheric conditions. Although NH3 and HCOOH catalysts play a negligible role in the reservoir for both radicals of HO and HO2 catalyzing the destruction of O3, the current study provides a comprehensive example of how acidic and basic catalysts assisted the gas-phase reactions. The H2O2 + HO → HO2 + H2O reaction is an important reservoir for both radicals of HO and HO2 catalyzing the destruction of O3.![]()
Collapse
Affiliation(s)
- Tianlei Zhang
- Shaanxi Key Laboratory of Catalysis
- School of Chemical & Environment Science
- Shaanxi University of Technology
- Hanzhong
- P. R. China
| | - Mingjie Wen
- Shaanxi Key Laboratory of Catalysis
- School of Chemical & Environment Science
- Shaanxi University of Technology
- Hanzhong
- P. R. China
| | - Zhaopeng Zeng
- Shaanxi Key Laboratory of Catalysis
- School of Chemical & Environment Science
- Shaanxi University of Technology
- Hanzhong
- P. R. China
| | - Yousong Lu
- Key Laboratory for Macromolecular Science of Shaanxi Province
- School of Chemistry & Chemical Engineering
- Shaanxi Normal University
- Xi'an
- P. R. China
| | - Yan Wang
- Shaanxi Key Laboratory of Catalysis
- School of Chemical & Environment Science
- Shaanxi University of Technology
- Hanzhong
- P. R. China
| | - Wei Wang
- Shaanxi Key Laboratory of Catalysis
- School of Chemical & Environment Science
- Shaanxi University of Technology
- Hanzhong
- P. R. China
| | - Xianzhao Shao
- Shaanxi Key Laboratory of Catalysis
- School of Chemical & Environment Science
- Shaanxi University of Technology
- Hanzhong
- P. R. China
| | - Zhiyin Wang
- Shaanxi Key Laboratory of Catalysis
- School of Chemical & Environment Science
- Shaanxi University of Technology
- Hanzhong
- P. R. China
| | - Lily Makroni
- Key Laboratory for Macromolecular Science of Shaanxi Province
- School of Chemistry & Chemical Engineering
- Shaanxi Normal University
- Xi'an
- P. R. China
| |
Collapse
|
40
|
Zhang T, Wen M, Zhang Y, Lan X, Long B, Wang R, Yu X, Zhao C, Wang W. Atmospheric chemistry of the self-reaction of HO 2 radicals: stepwise mechanism versus one-step process in the presence of (H 2O) n (n = 1-3) clusters. Phys Chem Chem Phys 2019; 21:24042-24053. [PMID: 31646308 DOI: 10.1039/c9cp03530c] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The effects of water on radical-radical reactions are of great importance for the elucidation of the atmospheric oxidation process of free radicals. In the present work, the HO2 + HO2 reactions with (H2O)n (n = 1-3) have been investigated using quantum chemical methods and canonical variational transition state theory with small curvature tunneling. We have explored both one-step and stepwise mechanisms, in particular the stepwise mechanism initiated by ring enlargement. The calculated results have revealed that the stepwise mechanism is the dominant one in the HO2 + HO2 reaction that is catalyzed by one water molecule. This is because its pseudo-first-order rate constant (kRWM1') is 3 orders of magnitude larger than that of the corresponding one-step mechanism. Additionally, the value of kRWM1' at 298 K has been found to be 4.3 times larger than that of the rate constant of the HO2 + HO2 reaction (kR1) without catalysts, which is in good agreement with the experimental findings. The calculated results also showed that the stepwise mechanism is still dominant in the (H2O)2 catalyzed reaction due to its higher pseudo-first-order rate constant, which is 3 orders of magnitude larger than that of the corresponding one-step mechanism. On the other hand, the one-step process is much faster than the stepwise mechanism by a factor of 105-106 in the (H2O)3 catalyzed reaction. However, the pseudo-first-order rate constants for the (H2O)2 and (H2O)3-catalyzed reactions are lower than that of the H2O-catalyzed reaction by 3-4 orders of magnitude, which indicates that the water monomer is the most efficient one among all the catalysts of (H2O)n (n = 1-3). The present results have provided a definitive example that water and water clusters have important influences on atmospheric reactions.
Collapse
Affiliation(s)
- Tianlei Zhang
- Shaanxi Key Laboratory of Catalysis, School of Chemical & Environment Science, Shaanxi University of Technology, Hanzhong, Shaanxi 723001, P. R. China.
| | | | | | | | | | | | | | | | | |
Collapse
|
41
|
Bakowies D. Estimating Systematic Error and Uncertainty in Ab Initio Thermochemistry. I. Atomization Energies of Hydrocarbons in the ATOMIC(hc) Protocol. J Chem Theory Comput 2019; 15:5230-5251. [DOI: 10.1021/acs.jctc.9b00343] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Dirk Bakowies
- Institute of Physical Chemistry, Department of Chemistry, University of Basel, Klingelbergstr. 80, CH 4056 Basel, Switzerland
| |
Collapse
|
42
|
Shee J, Arthur EJ, Zhang S, Reichman DR, Friesner RA. Singlet–Triplet Energy Gaps of Organic Biradicals and Polyacenes with Auxiliary-Field Quantum Monte Carlo. J Chem Theory Comput 2019; 15:4924-4932. [DOI: 10.1021/acs.jctc.9b00534] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Affiliation(s)
- James Shee
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Evan J. Arthur
- Schrodinger Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Shiwei Zhang
- Center for Computational Quantum Physics, Flatiron Institute, 162 Fifth Avenue, New York, New York 10010, United States
- Department of Physics, College of William and Mary, Williamsburg, Virginia 23187, United States
| | - David R. Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Richard A. Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| |
Collapse
|
43
|
Chan B. The CUAGAU Set of Coupled-Cluster Reference Data for Small Copper, Silver, and Gold Compounds and Assessment of DFT Methods. J Phys Chem A 2019; 123:5781-5788. [PMID: 31241947 DOI: 10.1021/acs.jpca.9b03976] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We have obtained benchmark data for a set of small molecular systems of Cu, Ag, and Au using coupled-cluster methods. Using this collection of reference data (that we termed the CUAGAU set) for assessing DFT-type methods, we find the MN15-L nonhybrid DFT to be cost-effective for geometry optimization [mean absolute deviation (MAD) in bond length = 0.20 Å], with an accuracy that is comparable to that for the double-hybrid (DH) DFT method DSD-PBEP86 (MAD = 0.19 Å). For the computation of thermochemical properties, among "conventional" (non-MP2-based) DFT methods, the best performance is found for the global-hybrid meta-GGA functional MN15, with an MAD of 11.4 kJ mol-1. We also find the nonhybrid method B97M-rV to have a reasonable performance (MAD = 14.4 kJ mol-1), and it may serve as a cost-effective means for qualitative study. If we look beyond conventional functionals, we find DSD-PBEP86 (MAD = 7.3 kJ mol-1) to be more accurate than even MN15. Nonetheless, this level of accuracy is still not sufficient for quantitative studies. In this regard, high-level wave function methods such as composite procedures that are based on coupled cluster are still indispensable for obtaining reliable reference data for transition-metal species.
Collapse
Affiliation(s)
- Bun Chan
- Graduate School of Engineering , Nagasaki University , Bunkyo 1-14 , Nagasaki 852-8521 , Japan
| |
Collapse
|
44
|
Chan B, Karton A, Raghavachari K. G4(MP2)-XK: A Variant of the G4(MP2)-6X Composite Method with Expanded Applicability for Main-Group Elements up to Radon. J Chem Theory Comput 2019; 15:4478-4484. [DOI: 10.1021/acs.jctc.9b00449] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Bun Chan
- Graduate School of Engineering, Nagasaki University, Bunkyo-machi 1-14, Nagasaki 852-8521, Japan
| | - Amir Karton
- School of Molecular Sciences, The University of Western Australia, Perth, WA 6009, Australia
| | - Krishnan Raghavachari
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States
| |
Collapse
|
45
|
Welch BK, Dawes R, Bross DH, Ruscic B. An Automated Thermochemistry Protocol Based on Explicitly Correlated Coupled-Cluster Theory: The Methyl and Ethyl Peroxy Families. J Phys Chem A 2019; 123:5673-5682. [DOI: 10.1021/acs.jpca.9b04381] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Bradley K. Welch
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65409, United States
| | - Richard Dawes
- Department of Chemistry, Missouri University of Science and Technology, Rolla, Missouri 65409, United States
| | - David H. Bross
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Branko Ruscic
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
- Consortium for Advanced Science and Engineering, University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
46
|
Rapid unimolecular reaction of stabilized Criegee intermediates and implications for atmospheric chemistry. Nat Commun 2019; 10:2003. [PMID: 31043594 PMCID: PMC6494847 DOI: 10.1038/s41467-019-09948-7] [Citation(s) in RCA: 57] [Impact Index Per Article: 11.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2018] [Accepted: 04/03/2019] [Indexed: 11/10/2022] Open
Abstract
Elucidating atmospheric oxidation mechanisms is necessary for estimating the lifetimes of atmospheric species and understanding secondary organic aerosol formation and atmospheric oxidation capacity. We report an unexpectedly fast mechanistic pathway for the unimolecular reactions of large stabilized Criegee intermediates, which involves the formation of bicyclic structures from large Criegee intermediates containing an aldehyde group. The barrier heights of the mechanistic pathways are unexpectedly low – about 2–3 kcal/mol – and are at least 10 kcal/mol lower than those of hydrogen shift processes in large syn Criegee intermediates; and the calculated rate constants show that the mechanistic pathways are 105-109 times faster than those of the corresponding hydrogen shift processes. The present findings indicate that analogous low-energy pathways can now also be expected in other large Criegee intermediates and that oxidative capacity of some Criegee intermediates is smaller than would be predicted by existing models. Criegee intermediates have received much attention in atmospheric chemistry because of their importance in ozonolysis mechanisms. Here, using quantum mechanical kinetics, the authors reveal an unexpectedly fast mechanistic pathway for unimolecular reactions of large stabilized Criegee intermediates.
Collapse
|
47
|
Trogolo D, Arey JS, Tentscher PR. Gas-Phase Ozone Reactions with a Structurally Diverse Set of Molecules: Barrier Heights and Reaction Energies Evaluated by Coupled Cluster and Density Functional Theory Calculations. J Phys Chem A 2019; 123:517-536. [DOI: 10.1021/acs.jpca.8b10323] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Daniela Trogolo
- School of Architecture, Civil and Environmental Engineering, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - J. Samuel Arey
- School of Architecture, Civil and Environmental Engineering, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
- Eawag, Swiss Federal Institute of Aquatic Science and Technology, Überlandstrasse 133, 8600 Dübendorf, Switzerland
| | - Peter R. Tentscher
- Eawag, Swiss Federal Institute of Aquatic Science and Technology, Überlandstrasse 133, 8600 Dübendorf, Switzerland
| |
Collapse
|
48
|
Mai TVT, Nguyen HT, Huynh LK. Ab initio dynamics of hydrogen abstraction from N2H4 by OH radicals: an RRKM-based master equation study. Phys Chem Chem Phys 2019; 21:23733-23741. [DOI: 10.1039/c9cp04585f] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The detailed kinetic mechanism of the N2H4 + OH reaction is comprehensively reported for a wide condition range of conditions (i.e., 200–3000 K & 1–7600 Torr) using the CCSD(T)/CBS//M06-2X/6-311++G(3df,2p) level and the RRKM-based master equation rate model.
Collapse
Affiliation(s)
- Tam V.-T. Mai
- Molecular Science and Nano-Materials Lab
- Institute for Computational Science and Technology
- Ho Chi Minh City
- Vietnam
- University of Science
| | - Hieu T. Nguyen
- Molecular Science and Nano-Materials Lab
- Institute for Computational Science and Technology
- Ho Chi Minh City
- Vietnam
| | - Lam K. Huynh
- International University
- Vietnam National University – HCMC
- Ho Chi Minh City
- Vietnam
| |
Collapse
|
49
|
Long B, Bao JL, Truhlar DG. Kinetics of the Strongly Correlated CH3O + O2 Reaction: The Importance of Quadruple Excitations in Atmospheric and Combustion Chemistry. J Am Chem Soc 2018; 141:611-617. [DOI: 10.1021/jacs.8b11766] [Citation(s) in RCA: 43] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Bo Long
- College of Materials Science and Engineering, Guizhou Minzu University, Guiyang, 550025, China
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Junwei Lucas Bao
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Donald G. Truhlar
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
50
|
Zhao Y, Xia L, Liao X, He Q, Zhao MX, Truhlar DG. Extrapolation of high-order correlation energies: the WMS model. Phys Chem Chem Phys 2018; 20:27375-27384. [PMID: 30357169 DOI: 10.1039/c8cp04973d] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We have developed a new composite model chemistry method called WMS (Wuhan-Minnesota scaling method) with three characteristics: (1) a composite scheme to approximate the complete configuration interaction valence energy with the affordability condition of requiring no calculation more expensive than CCSD(T)/jul-cc-pV(T+d)Z, (2) low-cost methods for the inner-shell correlation contribution and scalar relativistic correction, and (3) accuracy comparable to methods with post-CCSD(T) components. The new method is shown to be accurate for the W4-17 database of 200 atomization energies with an average mean unsigned error (averaged with equal weight over strongly correlated and weakly correlated subsets of the data) of 0.45 kcal mol-1, and the performance/cost ratio of these results compares very favorably to previously available methods. We also assess the WMS method against the DBH24-W4 database of diverse barrier heights and the energetics of the reactions of three strongly correlated Criegee intermediates with water. These results demonstrate that higher-order correlation contributions necessary to obtain high accuracy for molecular thermochemistry may be successfully extrapolated from the lower-order components of CCSD(T) calculations, and chemical accuracy can now be obtained for larger and more complex molecules and reactions.
Collapse
Affiliation(s)
- Yan Zhao
- State Key Laboratory of Silicate Materials for Architectures, International School of Materials Science and Engineering, Wuhan University of Technology, Wuhan 430070, People's Republic of China.
| | | | | | | | | | | |
Collapse
|