1
|
Li H, Mansoori Kermani M, Ottochian A, Crescenzi O, Janesko BG, Truhlar DG, Scalmani G, Frisch MJ, Ciofini I, Adamo C. Modeling Multi-Step Organic Reactions: Can Density Functional Theory Deliver Misleading Chemistry? J Am Chem Soc 2024; 146:6721-6732. [PMID: 38413362 DOI: 10.1021/jacs.3c12713] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/29/2024]
Abstract
Many organic reactions are characterized by a complex mechanism with a variety of transition states and intermediates of different chemical natures. Their correct and accurate theoretical characterization critically depends on the accuracy of the computational method used. In this work, we study a complex ambimodal cycloaddition with five transition states, two intermediates, and three products, and we ask whether density functional theory (DFT) can provide a correct description of this type of complex and multifaceted reaction. Our work fills a gap in that most systematic benchmarks of DFT for chemical reactions have considered much simpler reactions. Our results show that many density functionals not only lead to seriously large errors but also differ from one another in predicting whether the reaction is ambimodal. Only a few of the available functionals provide a balanced description of the complex and multifaceted reactions. The parameters varied in the tested functionals are the ingredients, the treatment of medium-range and nonlocal correlation energy, and the inclusion of Hartree-Fock exchange. These results show a clear need for more benchmarks on the mechanisms of large molecules in complex reactions.
Collapse
Affiliation(s)
- Hanwei Li
- Chimie ParisTech, PSL Research University, CNRS, Institute of Chemistry for Life and Health Sciences, Paris F-75005, France
| | - Maryam Mansoori Kermani
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Alistar Ottochian
- Chimie ParisTech, PSL Research University, CNRS, Institute of Chemistry for Life and Health Sciences, Paris F-75005, France
| | - Orlando Crescenzi
- Dipartimento di Scienze Chimiche, Università di Napoli Federico II, Complesso Universitario di Monte Sant'Angelo, Via Cinthia, Napoli 80126, Italy
| | - Benjamin G Janesko
- Department of Chemistry & Biochemistry, Texas Christian University, Fort Worth, Texas 76129, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | | | | | - Ilaria Ciofini
- Chimie ParisTech, PSL Research University, CNRS, Institute of Chemistry for Life and Health Sciences, Paris F-75005, France
| | - Carlo Adamo
- Chimie ParisTech, PSL Research University, CNRS, Institute of Chemistry for Life and Health Sciences, Paris F-75005, France
- Institut Universitaire de France, 103 Boulevard Saint Michel, Paris F-75005, France
| |
Collapse
|
2
|
Ram H, Sadej TP, Murphy CC, Mallo TJ, Westmoreland PR. Thermochemistry of Species in Gas-Phase Thermal Oxidation of C 2 to C 8 Perfluorinated Carboxylic Acids. J Phys Chem A 2024; 128:1313-1326. [PMID: 38335280 DOI: 10.1021/acs.jpca.3c06937] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/12/2024]
Abstract
New thermochemical properties, Cp°(T), H°(T), S°(T), and G°(T), are predicted for 123 species involved in the thermal destruction of perfluorinated carboxylic acids (PFCAs) using computational quantum chemistry and ideal-gas statistical mechanics. Relevant species were identified from the development of mechanisms for the pyrolysis and oxidation of PFCAs of C2 to C8 in length. Partition functions were obtained from the results of calculations at the G4 level for species up to C4 in length and M06-2X-D3(0)/def2-QZVPP for species C5 to C8 in length. The 1D hindered-rotor approximation was used to correct for torsional modes in the larger species. Ideal-gas thermochemistry was computed and fitted to 7-parameter NASA polynomials over a 200-2500 K temperature range, and the data are provided in standardized format. To gauge the effects of both method and basis set choice, enthalpies of formation at 0 K are calculated from various other density functionals (including B3LYP and ωB97XD), basis sets, and composite model chemistries (CBS-QB3). They are benchmarked against data from the Active Thermochemical Tables, high-level ANL0 calculations from the literature, and G4 calculations from this work. The effects of internal rotations and other anharmonicities are discussed, and bond dissociation energies and reaction equilibria provide mechanistic insights.
Collapse
Affiliation(s)
- Hrishikesh Ram
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27606, United States
| | - Thomas P Sadej
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27606, United States
| | - C Claire Murphy
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27606, United States
| | - Tim J Mallo
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27606, United States
| | - Phillip R Westmoreland
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27606, United States
| |
Collapse
|
3
|
Le TH, Ferro-Costas D, Fernández-Ramos A, Ortuño MA. Combined DFT and Kinetic Monte Carlo Study of UiO-66 Catalysts for γ-Valerolactone Production. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2024; 128:1049-1057. [PMID: 38293690 PMCID: PMC10823797 DOI: 10.1021/acs.jpcc.3c06053] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/08/2023] [Revised: 12/21/2023] [Accepted: 01/02/2024] [Indexed: 02/01/2024]
Abstract
Zr-based metal-organic frameworks (MOFs) are excellent heterogeneous porous catalysts due to their thermal stability. Their tunability via node and linker modifications makes them amenable for theoretical studies on catalyst design. However, detailed benchmarks on MOF-based reaction mechanisms combined with kinetics analysis are still scarce. Thus, we here evaluate different computational models and density functional theory (DFT) methods followed by kinetic Monte Carlo studies for a case reaction relevant in biomass upgrading, i.e., the conversion of methyl levulinate to γ-valerolactone catalyzed by UiO-66. We show the impact of cluster versus periodic models, the importance of the DF of choice, and the direct comparison to experimental data via simulated kinetics data. Overall, we found that Perdew-Burke-Ernzerhof (PBE), a widely employed method in plane-wave periodic calculations, greatly overestimates reaction rates, while M06 with cluster models better fits the available experimental data and is recommended whenever possible.
Collapse
Affiliation(s)
- Thanh-Hiep
Thi Le
- Centro
Singular de Investigación en Química Biolóxica
e Materiais Moleculares (CIQUS), Universidade
de Santiago de Compostela, 15782 Santiago de Compostela, Spain
| | - David Ferro-Costas
- Centro
Singular de Investigación en Química Biolóxica
e Materiais Moleculares (CIQUS), Universidade
de Santiago de Compostela, 15782 Santiago de Compostela, Spain
- Departamento
de Química Física, Facultade de Química, Universidade de Santiago de Compostela, 15782 Santiago
de Compostela, Spain
| | - Antonio Fernández-Ramos
- Centro
Singular de Investigación en Química Biolóxica
e Materiais Moleculares (CIQUS), Universidade
de Santiago de Compostela, 15782 Santiago de Compostela, Spain
- Departamento
de Química Física, Facultade de Química, Universidade de Santiago de Compostela, 15782 Santiago
de Compostela, Spain
| | - Manuel A. Ortuño
- Centro
Singular de Investigación en Química Biolóxica
e Materiais Moleculares (CIQUS), Universidade
de Santiago de Compostela, 15782 Santiago de Compostela, Spain
| |
Collapse
|
4
|
Pederson MR, Withanage KPK, Hooshmand Z, Johnson AI, Baruah T, Yamamoto Y, Zope RR, Kao DY, Shukla PB, Johnson JK, Peralta JE, Jackson KA. Use of FLOSIC for understanding anion-solvent interactions. J Chem Phys 2023; 159:154112. [PMID: 37861122 DOI: 10.1063/5.0172300] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2023] [Accepted: 09/20/2023] [Indexed: 10/21/2023] Open
Abstract
An Achille's heel of lower-rung density-functional approximations is that the highest-occupied-molecular-orbital energy levels of anions, known to be stable or metastable in nature, are often found to be positive in the worst case or above the lowest-unoccupied-molecular-orbital levels on neighboring complexes that are not expected to accept charge. A trianionic example, [Cr(C2O4)3]3-, is of interest for constraining models linking Cr isotope ratios in rock samples to oxygen levels in Earth's atmosphere over geological timescales. Here we describe how crowd sourcing can be used to carry out self-consistent Fermi-Löwdin-Orbital-Self-Interaction corrected calculations (FLOSIC) on this trianion in solution. The calculations give a physically correct description of the electronic structure of the trianion and water. In contrast, uncorrected local density approximation (LDA) calculations result in approximately half of the anion charge being transferred to the water bath due to the effects of self-interaction error. Use of group-theory and the intrinsic sparsity of the theory enables calculations roughly 125 times faster than our initial implementation in the large N limit reached here. By integrating charge density densities and Coulomb potentials over regions of space and analyzing core-level shifts of the Cr and O atoms as a function of position and functional, we unambiguously show that FLOSIC, relative to LDA, reverses incorrect solute-solvent charge transfer in the trianion-water complex. In comparison to other functionals investigated herein, including Hartree-Fock and the local density approximation, the FLOSIC Cr 1s eigenvalues provide the best agreement with experimental core ionization energies.
Collapse
Affiliation(s)
- Mark R Pederson
- Physics Department, University of Texas at El Paso, El Paso, Texas 79968, USA
| | | | - Zahra Hooshmand
- Physics Department, University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Alex I Johnson
- Physics Department, University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Tunna Baruah
- Physics Department, University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Yoh Yamamoto
- Physics Department, University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Rajendra R Zope
- Physics Department, University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Der-You Kao
- NASA Postdoctoral Program, NASA Goddard Space Flight Center, Greenbelt, Maryland 20771, USA
| | - Priyanka B Shukla
- Department of Chemical and Petroleum Engineering, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, USA
| | - J Karl Johnson
- Department of Chemical and Petroleum Engineering, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, USA
| | - Juan E Peralta
- Physics Department, Central Michigan University, Mt. Pleasant, Michigan 48859, USA
| | - Koblar A Jackson
- Physics Department, Central Michigan University, Mt. Pleasant, Michigan 48859, USA
| |
Collapse
|
5
|
Kee CW. Molecular Understanding and Practical In Silico Catalyst Design in Computational Organocatalysis and Phase Transfer Catalysis-Challenges and Opportunities. Molecules 2023; 28:molecules28041715. [PMID: 36838703 PMCID: PMC9966076 DOI: 10.3390/molecules28041715] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2022] [Revised: 02/03/2023] [Accepted: 02/05/2023] [Indexed: 02/25/2023] Open
Abstract
Through the lens of organocatalysis and phase transfer catalysis, we will examine the key components to calculate or predict catalysis-performance metrics, such as turnover frequency and measurement of stereoselectivity, via computational chemistry. The state-of-the-art tools available to calculate potential energy and, consequently, free energy, together with their caveats, will be discussed via examples from the literature. Through various examples from organocatalysis and phase transfer catalysis, we will highlight the challenges related to the mechanism, transition state theory, and solvation involved in translating calculated barriers to the turnover frequency or a metric of stereoselectivity. Examples in the literature that validated their theoretical models will be showcased. Lastly, the relevance and opportunity afforded by machine learning will be discussed.
Collapse
Affiliation(s)
- Choon Wee Kee
- Institute of Sustainability for Chemicals, Energy and Environment (ISCE2), Agency for Science, Technology and Research (A*STAR), 1 Pesek Road, Jurong Island, Singapore 627833, Republic of Singapore
| |
Collapse
|