1
|
Barretta P, Scoditti S, Belletto D, Ponte F, Vigna V, Mazzone G, Sicilia E. Ruthenium complexes bearing nile red chromophore and one of its derivative: Theoretical evaluation of PDT-related properties. J Comput Chem 2024; 45:2034-2041. [PMID: 38733370 DOI: 10.1002/jcc.27392] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2024] [Revised: 03/25/2024] [Accepted: 04/24/2024] [Indexed: 05/13/2024]
Abstract
The outcomes of DFT-based calculations are here reported to assess the applicability of two synthesized polypyridyl Ru(II) complexes, bearing ethynyl nile red (NR) on a bpy ligand, and two analogues, bearing modified-NR, in photodynamic therapy. The absorption spectra, together with the non-radiative rate constants for the S1 - Tn intersystem crossing transitions, have been computed for this purpose. Calculations evidence that the structural modification on the chromophore destabilizes the HOMO of the complexes thus reducing the H-L gap and, consequently, red shifting the maximum absorption wavelength within the therapeutic window, up to 620 nm. Moreover, the favored ISC process from the bright state involves the triplet state closest in energy, which is also characterized by the highest SOC value and by the involvement of the whole bpy ligand bearing the chromophore in delocalising the unpaired electrons. These outcomes show that the photophysical behavior of the complexes is dominated by the chromophore.
Collapse
Affiliation(s)
- Pierraffaele Barretta
- Department of Chemistry and Chemical Technologies, Università della Calabria, Arcavacata di Rende (CS), Italy
| | - Stefano Scoditti
- Department of Chemistry and Chemical Technologies, Università della Calabria, Arcavacata di Rende (CS), Italy
| | - Daniele Belletto
- Department of Chemistry and Chemical Technologies, Università della Calabria, Arcavacata di Rende (CS), Italy
| | - Fortuna Ponte
- Department of Chemistry and Chemical Technologies, Università della Calabria, Arcavacata di Rende (CS), Italy
| | - Vincenzo Vigna
- Department of Chemistry and Chemical Technologies, Università della Calabria, Arcavacata di Rende (CS), Italy
| | - Gloria Mazzone
- Department of Chemistry and Chemical Technologies, Università della Calabria, Arcavacata di Rende (CS), Italy
| | - Emilia Sicilia
- Department of Chemistry and Chemical Technologies, Università della Calabria, Arcavacata di Rende (CS), Italy
| |
Collapse
|
2
|
Xu R, Jiang Z, Yang Q, Bloino J, Biczysko M. Harmonic and anharmonic vibrational computations for biomolecular building blocks: Benchmarking DFT and basis sets by theoretical and experimental IR spectrum of glycine conformers. J Comput Chem 2024; 45:1846-1869. [PMID: 38682874 DOI: 10.1002/jcc.27377] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2024] [Revised: 04/01/2024] [Accepted: 04/02/2024] [Indexed: 05/01/2024]
Abstract
Advanced vibrational spectroscopic experiments have reached a level of sophistication that can only be matched by numerical simulations in order to provide an unequivocal analysis, a crucial step to understand the structure-function relationship of biomolecules. While density functional theory (DFT) has become the standard method when targeting medium-size or larger systems, the problem of its reliability and accuracy are well-known and have been abundantly documented. To establish a reliable computational protocol, especially when accuracy is critical, a tailored benchmark is usually required. This is generally done over a short list of known candidates, with the basis set often fixed a priori. In this work, we present a systematic study of the performance of DFT-based hybrid and double-hybrid functionals in the prediction of vibrational energies and infrared intensities at the harmonic level and beyond, considering anharmonic effects through vibrational perturbation theory at the second order. The study is performed for the six-lowest energy glycine conformers, utilizing available "state-of-the-art" accurate theoretical and experimental data as reference. Focusing on the most intense fundamental vibrations in the mid-infrared range of glycine conformers, the role of the basis sets is also investigated considering the balance between computational cost and accuracy. Targeting larger systems, a broad range of hybrid schemes with different computational costs is also tested.
Collapse
Affiliation(s)
- Ruiqin Xu
- Department of Physics, College of Sciences, Shanghai University, Shanghai, China
| | | | - Qin Yang
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Science, Prague, Czechia
| | - Julien Bloino
- Classe di Scienze, Scuola Normale Superiore, Pisa, Italy
| | - Malgorzata Biczysko
- Department of Physics, College of Sciences, Shanghai University, Shanghai, China
| |
Collapse
|
3
|
Gutiérrez-Peña CL, Gutiérrez-Blanco A, Gusev DG, Poyatos M, Peris E. Lone-Pair-π Bond Strength Unveiled by a Combined Experimental and Computational Study. Angew Chem Int Ed Engl 2024; 63:e202407817. [PMID: 38748473 DOI: 10.1002/anie.202407817] [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: 04/24/2024] [Indexed: 07/02/2024]
Abstract
A series of naphthalene-diimide (NDI) and perylene-diimide (PDI) connected bis-N-heterocyclic carbene complexes of iridium(III) have been prepared and fully characterized. The analysis of their NMR spectroscopic features, together with their molecular structures show that these species display lone-pair-π interactions between the chloride ligands of the Ir(III) complex and the heterocycles of the NDI/PDI moieties. The detection of this type of interaction in solution is due to the formation of two atropisomers, which are formed as a result of the restricted rotation about the Ir-Ccarbene bond imposed by the (Cl)lp⋅⋅⋅π interaction. Variable-temperature 1H NMR analysis allowed the determination of the strength of this non-covalent interaction, which lies between ΔH=6.6 and 10 kcal/mol. The computational studies performed fully support the experimental findings.
Collapse
Affiliation(s)
- Cristian L Gutiérrez-Peña
- Institute of Advanced Materials (INAM)., Universitat Jaume I, Av. Vicente Sos Baynat s/n., Castellón., E-12071., Spain
| | - Ana Gutiérrez-Blanco
- Institute of Advanced Materials (INAM)., Universitat Jaume I, Av. Vicente Sos Baynat s/n., Castellón., E-12071., Spain
| | - Dmitry G Gusev
- Department of Chemistry and Biochemistry, Wilfrid Laurier University, 75 University Avenue West, Waterloo, Ontario, N2 L3 C5, Canada
| | - Macarena Poyatos
- Institute of Advanced Materials (INAM)., Universitat Jaume I, Av. Vicente Sos Baynat s/n., Castellón., E-12071., Spain
| | - Eduardo Peris
- Institute of Advanced Materials (INAM)., Universitat Jaume I, Av. Vicente Sos Baynat s/n., Castellón., E-12071., Spain
| |
Collapse
|
4
|
Voss J. Machine learning for accuracy in density functional approximations. J Comput Chem 2024; 45:1829-1845. [PMID: 38668453 DOI: 10.1002/jcc.27366] [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: 10/30/2023] [Revised: 02/16/2024] [Accepted: 03/25/2024] [Indexed: 07/21/2024]
Abstract
Machine learning techniques have found their way into computational chemistry as indispensable tools to accelerate atomistic simulations and materials design. In addition, machine learning approaches hold the potential to boost the predictive power of computationally efficient electronic structure methods, such as density functional theory, to chemical accuracy and to correct for fundamental errors in density functional approaches. Here, recent progress in applying machine learning to improve the accuracy of density functional and related approximations is reviewed. Promises and challenges in devising machine learning models transferable between different chemistries and materials classes are discussed with the help of examples applying promising models to systems far outside their training sets.
Collapse
Affiliation(s)
- Johannes Voss
- SUNCAT Center for Interface Science and Catalysis, SLAC National Accelerator Laboratory, Menlo Park, California, USA
| |
Collapse
|
5
|
Shu Y, Zhu Z, Kanchanakungwankul S, Truhlar DG. Small Representative Databases for Testing and Validating Density Functionals and Other Electronic Structure Methods. J Phys Chem A 2024. [PMID: 39046238 DOI: 10.1021/acs.jpca.4c03137] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/25/2024]
Abstract
Broad and diverse sets of accurate data provide useful metrics for assessing the performance of new theoretical methods. However, assessing methods against large databases can be an arduous task. Here, we present 17 representative energetic databases, defined as small databases whose errors and error spreads are representative of larger databases and which therefore can serve as efficient benchmarks for developing and testing electronic structure methods and density functionals. In 15 cases, the representative databases have 6 entries while being representative of larger databases with 14-107 entries, and in the other two cases, they have 14 entries while being representative of larger databases with 418-455 entries. The mean unsigned error (MUE) of 100 electronic structure methods on a given representative database is typically within about 8% of the MUE on its parent database, and the root-mean-square error (RMSE) is typically within about 11% of the RMSE on the parent database. Thus, the representative databases are quite successful in indicating accuracy while maintaining good diversity. The databases include both main-group and transition-metal compounds and reactions, and they include bond energies, reaction energies, barrier heights, noncovalent interactions, ionization potentials, and absolute energies.
Collapse
Affiliation(s)
- Yinan Shu
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Zhaohan Zhu
- Department of Biostatistics, University of Florida, Gainesville, Florida 32611, United States
| | - Siriluk Kanchanakungwankul
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
6
|
Neufeld O, Tancogne-Dejean N, Rubio A. Benchmarking Functionals for Strong-Field Light-Matter Interactions in Adiabatic Time-Dependent Density Functional Theory. J Phys Chem Lett 2024; 15:7254-7264. [PMID: 38976844 PMCID: PMC11261632 DOI: 10.1021/acs.jpclett.4c01383] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2024] [Revised: 06/29/2024] [Accepted: 07/03/2024] [Indexed: 07/10/2024]
Abstract
In recent years, time-dependent density functional theory (TDDFT) has been extensively employed for highly nonlinear optics in molecules and solids, including high harmonic generation (HHG), photoemission, and more. TDDFT exhibits a relatively low numerical cost while still describing both light-matter and electron-electron interactions ab initio, making it highly appealing. However, the majority of implementations of the theory utilize the simplest possible approximations for the exchange-correlation (XC) functional-either the local density or generalized gradient approximations, which are traditionally considered to have rather poor chemical accuracy. We present the first systematic study of the XC functional effect on molecular HHG, testing various levels of theory. Our numerical results suggest justification for using simpler approximations for the XC functional, showing that hybrid and meta functionals (as well as Hartree-Fock) can, at times, lead to poor and unphysical results. The specific source of the failure in more elaborate functionals should be topic of future work, but we hypothesize that its origin might be connected to the adiabatic approximation of TDDFT.
Collapse
Affiliation(s)
- Ofer Neufeld
- Max
Planck Institute for the Structure and Dynamics of Matter and Center
for Free-Electron Laser Science, Hamburg 22761, Germany
| | - Nicolas Tancogne-Dejean
- Max
Planck Institute for the Structure and Dynamics of Matter and Center
for Free-Electron Laser Science, Hamburg 22761, Germany
| | - Angel Rubio
- Max
Planck Institute for the Structure and Dynamics of Matter and Center
for Free-Electron Laser Science, Hamburg 22761, Germany
- Center
for Computational Quantum Physics (CCQ), The Flatiron Institute, New York, New York 10010, United States
| |
Collapse
|
7
|
Kania MJ, Reyes A, Neufeldt SR. Oxidative Addition of (Hetero)aryl (Pseudo)halides at Palladium(0): Origin and Significance of Divergent Mechanisms. J Am Chem Soc 2024; 146:19249-19260. [PMID: 38959060 DOI: 10.1021/jacs.4c04496] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/04/2024]
Abstract
Two limiting mechanisms are possible for oxidative addition of (hetero)aryl (pseudo)halides at Pd(0): a 3-centered concerted and a nucleophilic displacement mechanism. Until now, there has been little understanding about when each mechanism is relevant. Prior investigations to distinguish between these pathways were limited to a few specific combinations of the substrate and ligand. Here, we computationally evaluated over 180 transition structures for oxidative addition in order to determine mechanistic trends based on substrate, ligand(s), and coordination number. Natural abundance 13C kinetic isotope effects provide experimental results consistent with computational predictions. Key findings include that (1) differences in highest occupied molecular orbital (HOMO) symmetries dictate that, although 12e- PdL is strongly biased toward a 3-centered concerted mechanism, 14e- PdL2 often prefers a nucleophilic displacement mechanism; (2) ligand electronics and sterics, including ligand bite angle, influence the preferred mechanism of the reaction at PdL2; (3) phenyl triflate always reacts through a displacement mechanism regardless of the catalyst structure due to the stability of a triflate anion and the inability of oxygen to effectively donate electron density to Pd; and (4) the high reactivity of C-X bonds adjacent to nitrogen in pyridine substrates relates to stereoelectronic stabilization of a nucleophilic displacement transition state. This work has implications for controlling rate and selectivity in catalytic couplings, and we demonstrate application of the mechanistic insight toward chemodivergent cross-couplings of bromochloroheteroarenes.
Collapse
Affiliation(s)
- Matthew J Kania
- Department of Chemistry and Biochemistry, Montana State University, Bozeman, Montana 59717, United States
| | - Albert Reyes
- Department of Chemistry and Biochemistry, Montana State University, Bozeman, Montana 59717, United States
| | - Sharon R Neufeldt
- Department of Chemistry and Biochemistry, Montana State University, Bozeman, Montana 59717, United States
| |
Collapse
|
8
|
Bhattacharyya HP, Sarma M. Efficiency Conceptualization Model: A Theoretical Method for Predicting the Turnover of Catalysts. Chemphyschem 2024; 25:e202400004. [PMID: 38619023 DOI: 10.1002/cphc.202400004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/02/2024] [Revised: 04/11/2024] [Accepted: 04/15/2024] [Indexed: 04/16/2024]
Abstract
In recent times, the theoretical prediction of catalytic efficiency is of utmost urgency. With the advent of density functional theory (DFT), reliable computations can delineate a quantitative aspect of the study. To this state-of-the-art approach, valuable incorporation would be a tool that can acknowledge the efficiency of a catalyst. In the current work, we developed the efficiency conceptualization model (ECM) that utilizes the quantum mechanical tool to achieve efficiency in terms of turnover frequency (TOF). Twenty-six experimentally designed transition metal (TM) water oxidation catalysts were chosen under similar experimental conditions of temperature, pressure, and pH to execute the same. The computations conclude that the Fe-based [Fe(OTf)2(Me2Pytacn)] (MWOC-17) is a highly active catalyst and, therefore, can endure for more time in the catalytic cycle. Our results conclude that the Ir-based catalysts [Cp*Ir(κ2-N,O)X] with MWOC-23: X=Cl; and MWOC-24: X=NO3 report the highest computed turnover numbers (TONs),τ c o m p u t e d T O N 0 ${\tau _{computed\;TON}^0 }$ of 406 and 490 against the highest experimental TONs,τ e x p e r i m e n t a l T O N ${\tau _{experimental\;TON} }$ of 1200 and 2000 respectively, whereas the Co-based [Co(12-TMC)]2+ (MWOC-19) has the lowest TONs (τ c o m p u t e d T O N 0 ${\tau _{computed\;TON}^0 }$ =19, τexperimental TON=16) among the chosen catalysts and thereby successful in corroborating the previous experimental results.
Collapse
Affiliation(s)
| | - Manabendra Sarma
- Department of Chemistry, Indian Institute of Technology Guwahati, Guwahati, Assam, 781039, India
| |
Collapse
|
9
|
Ariyarathna IR. Electronic structure analysis and DFT benchmarking of Rydberg-type alkali-metal-crown ether, -cryptand, and -adamanzane complexes. Phys Chem Chem Phys 2024; 26:16989-16997. [PMID: 38666396 DOI: 10.1039/d4cp00723a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/20/2024]
Abstract
Density functional theory (DFT) and electron propagator theory (EPT) calculations were performed to study ground and excited electronic structures of alkali-metal (M) coordinated 9-crown-3, 24-crown-8, [2.1.1]cryptand, o-Me2-1.1.1, and 36Adamanzane complexes. Each complex bears an expanded electron in the periphery and occupies diffuse 1p-, 1d-, 1f-type molecular orbitals (or superatomic 1P, 1D, 1F orbitals) in excited electronic states. The calculated superatomic shell model of the M(9-crown-3)2 is 1S, 1P, 1D, 1F, 2S, 2P, 2D, 1G and it is held by all other complexes up to the studied 1F level. Due to the highly diffuse nature of the electron, the ionization energies of these complexes are significantly lower (1.6-2.0 eV) and hence these complexes belong to the superalkali category. The ab initio EPT ionization energy and the excitation energies of the Li(9-crown-3)2 were used to evaluate DFT errors associated with a series of exchange correlation functionals that span multiple rungs of Jacob's ladder (i.e., GGA, meta-GGA, global GGA hybrid, meta-GGA hybrid, range-separated hybrid, double-hybrid). Among these, the best performing functional is the range-separated hybrid CAM-B3LYP and the errors are within 6% of high-level ab initio EPT results. The accuracy of CAM-B3LYP is indeed transferable to similar complexes and hence the findings are expected to accelerate the progression of studies of Rydberg-type systems.
Collapse
Affiliation(s)
- Isuru R Ariyarathna
- Physics and Chemistry of Materials (T-1), Los Alamos National Laboratory, Los Alamos, NM 87545, USA.
| |
Collapse
|
10
|
Jiang J, Yang J, Hong Q, Sun Q, Li J. Global Potential Energy Surfaces by Compressed-State Multistate Pair-Density Functional Theory for Hyperthermal Collisions in the O 2+O 2 System. Chemphyschem 2024; 25:e202400078. [PMID: 38526528 DOI: 10.1002/cphc.202400078] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2024] [Revised: 03/02/2024] [Accepted: 03/25/2024] [Indexed: 03/26/2024]
Abstract
Interactions between oxygen molecules play an important role in atmospheric chemistry and hypersonic flow chemistry in atmospheric entries. Recently, high-quality ab initio potential energy surface (PES) of the quintet O4 was reported by Paukku et al. [J. Chem. Phys. 147, 034301 (2017)]. 10543 configurations were sampled and calculated at the level of MS-CASPT2/maug-cc-pVTZ with scaled external correlation. The PES was fitted to a many-body (MB) form with the many-body part described by the permutationally invariant polynomial approach (MB-PIP). In this work, the PIP-Neural Network (PIP-NN) and MB-PIP-NN methods were used to refit the PES based on the same data by Paukku et al. Three PESs were compared. It was found that the performances differ significantly in the O+O3 region as well as in the long-range region. Therefore, additional 1300 points were sampled, and the efficient compressed-state multistate pair-density functional theory (CMS-PDFT) was used to calculate the electronic structure of these 1300 points and 10543 points by Paukku et al. Then, a completely new quintet PES was fitted using the MB-PIP-NN method. Based on this PES, the quasi-classical trajectory (QCT) approach was used to reveal all possible reaction channels for hyperthermal O2-O2 collisions.
Collapse
Affiliation(s)
- Jie Jiang
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Chemical Theory and Mechanism, Chongqing University, Chongqing, 401331, China
| | - Jiawei Yang
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Chemical Theory and Mechanism, Chongqing University, Chongqing, 401331, China
| | - Qizhen Hong
- State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy of Sciences, 100190, Beijing, China
| | - Quanhua Sun
- State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy of Sciences, 100190, Beijing, China
- School of Engineering Science, University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Jun Li
- School of Chemistry and Chemical Engineering & Chongqing Key Laboratory of Chemical Theory and Mechanism, Chongqing University, Chongqing, 401331, China
| |
Collapse
|
11
|
Yu F. Origins of the unphysical noncovalent interaction energy curves obtained with the 2011 and 2012 Minnesota density functionals. J Chem Phys 2024; 160:214120. [PMID: 38836783 DOI: 10.1063/5.0212534] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2024] [Accepted: 05/22/2024] [Indexed: 06/06/2024] Open
Abstract
With the noncovalent interaction energy curves of the methane dimer [(CH4)2], we have clarified two different origins of the unphysical noncovalent interaction energy curves obtained with the Minnesota density functionals of M11-L, MN12-L, and MN12-SX. For the M11-L functional, the unphysical inflection point on the (CH4)2 interaction energy curve originates from the inclusion of the long-range exchange. As to the MN12-L and MN12-SX functionals, the lack of smoothness restraints results in unphysical inflection points on the corresponding (CH4)2 interaction energy curves. As a result, exchange functionals are as important as dispersion corrections for density functionals to map noncovalent interaction energy surfaces reasonably. Moreover, very highly parameterized functionals with smoothness restraints are suggested for investigating noncovalent interaction energy surfaces.
Collapse
Affiliation(s)
- Feng Yu
- Department of Physics, School of Freshmen, Xi'an Technological University, No. 4 Jinhua North Road, Xi'an, Shaanxi 710032, China
| |
Collapse
|
12
|
Yang M, Wang J. Comprehensive Multipath Variational Kinetics Study on Hydrogen Abstraction Reactions from Three Typical Dimethylcyclohexane Isomers by Hydroxyl Radicals: from the Electronic Structure to Model Applications. J Phys Chem A 2024; 128:4517-4531. [PMID: 38804972 DOI: 10.1021/acs.jpca.4c00480] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/29/2024]
Abstract
Cycloalkanes serve as an important class of chemical components in both fossil and alternative transportation fuels and have attracted considerable attention from the combustion community. Hydrogen abstractions from cycloalkanes by hydroxyl radicals initiate the fuel decomposition process and trigger off the subsequent chain reactions and thus play an important role in both combustion and atmospheric chemistry. The target of this study is to fill the vacancy in kinetics data toward the H-abstraction reactions by hydroxyl radical from three typical dimethylcyclohexane isomers through first-principles and direct dynamics. The rate constants involving 18 elementary reactions in total were accurately determined by the multipath canonical variational transition state theory with the multidimensional small-curvature correction for tunneling (MP-CVT/SCT), over a broad temperature range of 200-2000 K. The significant roles of multistructural torsional anharmonicity and recrossing effects were stressed per abstraction site, while the quantum tunneling effect was found to be slight at temperatures of interest in combustion. The discrepancies observed among different reaction systems at a similar abstraction site highlight the fuel molecular effects on site-specific rate constants. The comparison results of total rate constants given by different dynamics approaches prove the importance of considering the torsional anharmonicity, recrossing, and tunneling effects, and the robust feature of the simplified MS-CVT/SCT. The calculated total constants for dimethylcyclohexane isomers by OH are consistent with those measured for methylcyclohexane and 1,4-dimethylcyclohexane at low temperatures. The branching ratio analysis confirms the predominant role of the tertiary abstraction at low-to-intermediate temperatures and its growing competition with distinct secondary abstractions as temperature increases. The calculated rate constants were eventually fitted into the analytical expressions and incorporated into the kinetic models to learn about the influences on modeling performance.
Collapse
Affiliation(s)
- Mo Yang
- National Key Laboratory of Science and Technology on Aero-Engine Aero-Thermodynamics, School of Energy and Power Engineering, Beihang University, Beijing 100191, PR China
| | - Juan Wang
- National Key Laboratory of Science and Technology on Aero-Engine Aero-Thermodynamics, School of Energy and Power Engineering, Beihang University, Beijing 100191, PR China
| |
Collapse
|
13
|
Gao Q, Shen C, Zhang H, Long B, Truhlar DG. Quantitative kinetics reveal that reactions of HO 2 are a significant sink for aldehydes in the atmosphere and may initiate the formation of highly oxygenated molecules via autoxidation. Phys Chem Chem Phys 2024; 26:16160-16174. [PMID: 38787752 DOI: 10.1039/d4cp00693c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/26/2024]
Abstract
Large aldehydes are widespread in the atmosphere and their oxidation leads to secondary organic aerosols. The current understanding of their chemical transformation processes is limited to hydroxyl radical (OH) oxidation during daytime and nitrate radical (NO3) oxidation during nighttime. Here, we report quantitative kinetics calculations of the reactions of hexanal (C5H11CHO), pentanal (C4H9CHO), and butanal (C3H7CHO) with hydroperoxyl radical (HO2) at atmospheric temperatures and pressures. We find that neither tunneling nor multistructural torsion anharmonicity should be neglected in computing these rate constants; strong anharmonicity at the transition states is also important. We find rate constants for the three reactions in the range 3.2-7.7 × 10-14 cm3 molecule-1 s-1 at 298 K and 1 atm, showing that the HO2 reactions can be competitive with OH and NO3 oxidation under some conditions relevant to the atmosphere. Our findings reveal that HO2-initiated oxidation of large aldehydes may be responsible for the formation of highly oxygenated molecules via autoxidation. We augment the theoretic studies with laboratory flow-tube experiments using an iodide-adduct time-of-flight chemical ionization mass spectrometer to confirm the theoretical predictions of peroxy radicals and the autoxidation pathway. We find that the adduct from HO2 + C5H11CHO undergoes a fast unimolecular 1,7-hydrogen shift with a rate constant of 0.45 s-1. We suggest that the HO2 reactions make significant contributions to the sink of aldehydes.
Collapse
Affiliation(s)
- Qiao Gao
- School of Physics and Mechatronic Engineering, Guizhou Minzu University, Guiyang 550025, China.
| | - Chuanyang Shen
- Department of Chemistry, University of California, Riverside, California, 92507, USA.
| | - Haofei Zhang
- Department of Chemistry, University of California, Riverside, California, 92507, USA.
| | - Bo Long
- School of Physics and Mechatronic Engineering, Guizhou Minzu 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 Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, USA.
| |
Collapse
|
14
|
Achazi AJ, Fataj X, Rohland P, Hager MD, Schubert US, Mollenhauer D. Development of a multi-step screening procedure for redox active molecules in organic radical polymer anodes and as redox flow anolytes. J Comput Chem 2024; 45:1112-1129. [PMID: 38258532 DOI: 10.1002/jcc.27299] [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: 09/12/2023] [Revised: 12/13/2023] [Accepted: 12/15/2023] [Indexed: 01/24/2024]
Abstract
Benzo[d]-X-zolyl-pyridinyl (XO, S, NH) radicals represent a promising class of redox-active molecules for organic batteries. We present a multistep screening procedure to identify the most promising radical candidates. Experimental investigations and highly correlated wave function-based calculations are performed to determine benchmark redox potentials. Based on these, the accuracies of different methods (semi-empirical, density functional theory, wave function-based), solvent models, dispersion corrections, and basis sets are evaluated. The developed screening procedure consists of three steps: First, a conformer search is performed with CREST. The molecules are selected based on the redox potentials calculated using GFN2-xTB. Second, HOMO energies calculated with reparametrized B3LYP-D3(BJ) and the def2-SVP basis set are used as selection criteria. The final molecules are selected based on the redox potentials calculated from Gibbs energies using BP86-D3(BJ)/def2-TZVP. With this multistep screening approach, promising molecules can be suggested for synthesis, and structure-property relationships can be derived.
Collapse
Affiliation(s)
- Andreas J Achazi
- Physikalisch-Chemisches Institut, Justus-Liebig-Universität Gießen, Gießen, Germany
- Zentrum für Materialforschung, Justus-Liebig-Universität Gießen, Gießen, Germany
| | - Xhesilda Fataj
- Laboratory of Organic and Macromolecular Chemistry (IOMC), Friedrich Schiller University Jena, Jena, Germany
- Center for Energy and Environmental Chemistry Jena (CEEC Jena), Friedrich Schiller University Jena, Jena, Germany
| | - Philip Rohland
- Laboratory of Organic and Macromolecular Chemistry (IOMC), Friedrich Schiller University Jena, Jena, Germany
- Center for Energy and Environmental Chemistry Jena (CEEC Jena), Friedrich Schiller University Jena, Jena, Germany
| | - Martin D Hager
- Laboratory of Organic and Macromolecular Chemistry (IOMC), Friedrich Schiller University Jena, Jena, Germany
- Center for Energy and Environmental Chemistry Jena (CEEC Jena), Friedrich Schiller University Jena, Jena, Germany
| | - Ulrich S Schubert
- Laboratory of Organic and Macromolecular Chemistry (IOMC), Friedrich Schiller University Jena, Jena, Germany
- Center for Energy and Environmental Chemistry Jena (CEEC Jena), Friedrich Schiller University Jena, Jena, Germany
| | - Doreen Mollenhauer
- Physikalisch-Chemisches Institut, Justus-Liebig-Universität Gießen, Gießen, Germany
- Zentrum für Materialforschung, Justus-Liebig-Universität Gießen, Gießen, Germany
| |
Collapse
|
15
|
Rajput SS, Raghuvanshi N, Banana T, Yadav P, Alam MM. Why does the orientation of azulene affect the two-photon activity of a porphyrinoid-azulene system? Phys Chem Chem Phys 2024; 26:15611-15619. [PMID: 38758026 DOI: 10.1039/d4cp00438h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/18/2024]
Abstract
Attaching a dipolar molecule in a symmetric system induces a major change in the electronic structure, which may be reflected as the enhancement of the optical and charge-transfer properties of the combined system as compared to the pristine ones. Furthermore, the orientation of the dipolar molecule may also affect the said properties. This idea is explored in this work by taking porphyrinoid molecules as the pristine systems. We attached azulene, a dipolar molecule, at various positions of five porphyrinoid cores and studied the effect on charge-transfer and one- and two-photon absorption properties using the state-of-the-art RICC2 method. The attachment of azulene produces two major effects - firstly it introduces asymmetry in the system and, secondly, being dipolar, it makes the resultant molecule dipolar/quadrupolar. Porphyrin, N-confused porphyrin, sub-porphyrin, sapphyrin, and hexaphyrin are used as core porphyrinoid systems. The change in charge-transfer has been studied using the orbital analysis and charge-transfer distance parameter for the first five singlet states of the systems. The effect of orientation of azulene on the said properties is also explored. The insights gained from our observations are explored further at the dipole and transition dipole moment levels using a three-state model.
Collapse
Affiliation(s)
- Swati Singh Rajput
- Department of Chemistry, Indian Institute of Technology Bhilai, Durg, Chhattisgarh-491001, India.
| | - Nikita Raghuvanshi
- Centre for Basic Sciences, Pt. Ravishankar Shukla University, Raipur, Chhattisgarh, India
- School of Studies in Chemistry, Pt. Ravishankar Shukla University, Raipur, India
| | - Tejendra Banana
- Department of Chemistry, Indian Institute of Technology Bhilai, Durg, Chhattisgarh-491001, India.
| | - Pooja Yadav
- Department of Chemistry, Indian Institute of Technology Bhilai, Durg, Chhattisgarh-491001, India.
| | - Md Mehboob Alam
- Department of Chemistry, Indian Institute of Technology Bhilai, Durg, Chhattisgarh-491001, India.
| |
Collapse
|
16
|
Kumar N, Gupta P. DFT Struggles to Predict the Energy Landscape for Iron Pyridine Diimine-Catalyzed [2 + 2] Cycloaddition of Alkenes: Insights into the Problem and Alternative Solutions. J Phys Chem A 2024; 128:4114-4127. [PMID: 38659086 DOI: 10.1021/acs.jpca.3c08325] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/26/2024]
Abstract
In recent years, noninnocent pyridine diimine (PDI) complexes featuring first-row transition metals have emerged as prominent catalysts, demonstrating efficacy in a diverse range of vital organometallic transformations. However, the inherent complexity of the fundamental reactivity paradigm in these systems arises from the presence of a noninnocent ligand and the multispin feasibility of 3d metals. While density functional theory (DFT) has been widely used to unravel mechanistic insights, its limitations as a single-reference method can potentially misrepresent spin-state energetics, compromising our understanding of these intricate systems. In this study, we employ extensive high-level ab initio state averaged-complete active space self-consistent field/N-electron valence state perturbation theory (SA-CASSCF/NEVPT2) calculations in combination with DFT to investigate an iron-PDI-catalyzed [2 + 2] cycloaddition reaction of alkenes. The transformation proceeds through two major steps: oxidative cyclization and reductive elimination. Contrary to the predictions of DFT calculations, which suggest two-state reactivity in the reaction and identify reductive elimination as the turnover-limiting step, SA-CASSCF/NEVPT2-corrected results unequivocally establish a single-state reactivity scenario with oxidative cyclization as the turnover-limiting step. SA-CASSCF/NEVPT2-based insights into electronic ground states and electron distribution elucidate the intriguing interactions between the PDI ligand and the iron center, revealing the highly multiconfigurational nature of these species and providing a precise depiction of metal-ligand cooperativity throughout the transformation. A comparative assessment of several widely recognized DFT functionals against SA-CASSCF/NEVPT2-corrected data indicates that single-point energy calculations using the modern density functional MN15 on TPSSh geometries offer the most reliable density functional methodology, in scenarios where SA-CASSCF/NEVPT2 computational cost is a consideration.
Collapse
Affiliation(s)
- Nikunj Kumar
- Computational Catalysis Center, Department of Chemistry, Indian Institute of Technology Roorkee, Roorkee, Uttarakhand 247667, India
| | - Puneet Gupta
- Computational Catalysis Center, Department of Chemistry, Indian Institute of Technology Roorkee, Roorkee, Uttarakhand 247667, India
- Center for Sustainable Energy, Indian Institute of Technology Roorkee, Roorkee, Uttarakhand 247667, India
| |
Collapse
|
17
|
Havenridge S, Liu C. A Theoretical Benchmark of the Geometric and Optical Properties for 3d Transition Metal Nanoclusters via Density Functional Theory. J Phys Chem A 2024; 128:3947-3956. [PMID: 38729915 DOI: 10.1021/acs.jpca.4c00408] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/12/2024]
Abstract
Understanding structure-property relationships in atomically precise metal nanoclusters is vital in finding selective and tunable catalysts. In this study, density functional theory (DFT) was used to benchmark seven exchange correlation functionals at different basis sets for 17 atomically precise nanoclusters against experimentally determined geometries, band gaps, and optical gaps. The set contains both monometallic and bimetallic clusters that possess at least two types of 3d transition metals (specifically, Cu, Ni, Fe, or Co). The benchmark highlights that PBE0 is a good functional to use regardless of the basis set, and Minnesota functionals do well with respect to specific metals. Further, while long-range corrected functionals overestimate band and optical gaps, they model absorption features better than the other considered functionals. The study additionally looks at the photoinduced hydrogen evolution reaction (HER) and the CO2 reduction mechanism on nanoclusters reported from the literature.
Collapse
Affiliation(s)
- Shana Havenridge
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Cong Liu
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| |
Collapse
|
18
|
Zhan F, Lin GL, Zhang TS, Xu K, Yang YF, Li G, She Y. Pt-S Bond-Enabled Temperature-Dependent Phosphorescence in S-Heteroaryl Tetradentate Pt(S^C^N^O) Complexes. Inorg Chem 2024; 63:8822-8831. [PMID: 38696545 DOI: 10.1021/acs.inorgchem.4c00601] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/04/2024]
Abstract
This study presents the rare examples of S-heteroaryl tetradentate Pt(S^C^N^O) luminescent complexes (PtSZ and PtSZtBu) containing a Pt-S bond. The presence of the Pt-S bond allows the novel Pt(S^C^N^O) complexes to exhibit temperature-dependent phosphorescent emission behavior. The PtSZtBu exhibits dual-emission phenomena and biexponential transient decay spectra above 250 K, indicating the presence of two minimal excited states in the potential energy surface (PES) of the T1 state. Through complementary experimental and computational studies, we have identified changes in orbital composition between Pt(dxy)-S(px) and Pt(dyz)-S(pz) in excited states with increasing temperature. This results in two energy minima, enabling the excited states to decay selectively and radiatively at different temperatures. Consequently, this leads to remarkable steady-state and transient emission spectra changes. Our work not only provides valuable insights for the development of novel Pt-S bond-based tetradentate Pt(II) complexes but also enhances our understanding of the distinctive properties governed by the Pt-S bond.
Collapse
Affiliation(s)
- Feng Zhan
- College of Chemical Engineering, State Key Laboratory Breeding Base of Green-Chemical Synthesis Technology, Zhejiang University of Technology, Hangzhou, Zhejiang 310014, PR China
| | - Guo-Liang Lin
- College of Chemical Engineering, State Key Laboratory Breeding Base of Green-Chemical Synthesis Technology, Zhejiang University of Technology, Hangzhou, Zhejiang 310014, PR China
| | - Teng-Shuo Zhang
- College of Chemical Engineering, State Key Laboratory Breeding Base of Green-Chemical Synthesis Technology, Zhejiang University of Technology, Hangzhou, Zhejiang 310014, PR China
| | - Kewei Xu
- College of Chemical Engineering, State Key Laboratory Breeding Base of Green-Chemical Synthesis Technology, Zhejiang University of Technology, Hangzhou, Zhejiang 310014, PR China
| | - Yun-Fang Yang
- College of Chemical Engineering, State Key Laboratory Breeding Base of Green-Chemical Synthesis Technology, Zhejiang University of Technology, Hangzhou, Zhejiang 310014, PR China
| | - Guijie Li
- College of Chemical Engineering, State Key Laboratory Breeding Base of Green-Chemical Synthesis Technology, Zhejiang University of Technology, Hangzhou, Zhejiang 310014, PR China
| | - Yuanbin She
- College of Chemical Engineering, State Key Laboratory Breeding Base of Green-Chemical Synthesis Technology, Zhejiang University of Technology, Hangzhou, Zhejiang 310014, PR China
| |
Collapse
|
19
|
Gasevic T, Bursch M, Ma Q, Grimme S, Werner HJ, Hansen A. The p-block challenge: assessing quantum chemistry methods for inorganic heterocycle dimerizations. Phys Chem Chem Phys 2024; 26:13884-13908. [PMID: 38661329 DOI: 10.1039/d3cp06217a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/26/2024]
Abstract
The elements of the p-block of the periodic table are of high interest in various chemical and technical applications like frustrated Lewis-pairs (FLP) or opto-electronics. However, high-quality benchmark data to assess approximate density functional theory (DFT) for their theoretical description are sparse. In this work, we present a benchmark set of 604 dimerization energies of 302 "inorganic benzenes" composed of all non-carbon p-block elements of main groups III to VI up to polonium. This so-called IHD302 test set comprises two classes of structures formed by covalent bonding and by weaker donor-acceptor (WDA) interactions, respectively. Generating reliable reference data with ab initio methods is challenging due to large electron correlation contributions, core-valence correlation effects, and especially the slow basis set convergence. To compute reference values for these dimerization reactions, after thorough testing, we applied a computational protocol using state-of-the-art explicitly correlated local coupled cluster theory termed PNO-LCCSD(T)-F12/cc-VTZ-PP-F12(corr.). It includes a basis set correction at the PNO-LMP2-F12/aug-cc-pwCVTZ level. Based on these reference data, we assess 26 DFT methods in combination with three different dispersion corrections and the def2-QZVPP basis set, five composite DFT approaches, and five semi-empirical quantum mechanical methods. For the covalent dimerizations, the r2SCAN-D4 meta-GGA, the r2SCAN0-D4 and ωB97M-V hybrids, and the revDSD-PBEP86-D4 double-hybrid functional are found to be the best-performing methods among the evaluated functionals of the respective class. However, since def2 basis sets for the 4th period are not associated to relativistic pseudo-potentials, we obtained significant errors in the covalent dimerization energies (up to 6 kcal mol-1) for molecules containing p-block elements of the 4th period. Significant improvements were achieved for systems containing 4th row elements by using ECP10MDF pseudopotentials along with re-contracted aug-cc-pVQZ-PP-KS basis sets introduced in this work with the contraction coefficients taken from atomic DFT (PBE0) calculations. Overall, the IHD302 set represents a challenge to contemporary quantum chemical methods. This is due to a large number of spatially close p-element bonds which are underrepresented in other benchmark sets, and the partial covalent bonding character for the WDA interactions. The IHD302 set may be helpful to develop more robust and transferable approximate quantum chemical methods in the future.
Collapse
Affiliation(s)
- Thomas Gasevic
- Mulliken Center for Theoretical Chemistry, Rheinische Friedrich-Wilhelms-Universität Bonn, Beringstr. 4, 53115 Bonn, Germany.
| | - Markus Bursch
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany.
- FACCTs GmbH, 50677, Koeln, Germany
| | - Qianli Ma
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany.
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Rheinische Friedrich-Wilhelms-Universität Bonn, Beringstr. 4, 53115 Bonn, Germany.
| | - Hans-Joachim Werner
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany.
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Rheinische Friedrich-Wilhelms-Universität Bonn, Beringstr. 4, 53115 Bonn, Germany.
| |
Collapse
|
20
|
Larson NG, Norman JP, Neufeldt SR. Mechanistic Origin of Ligand Effects on Exhaustive Functionalization During Pd-Catalyzed Cross-Coupling of Dihaloarenes. ACS Catal 2024; 14:7127-7135. [PMID: 38911468 PMCID: PMC11192547 DOI: 10.1021/acscatal.4c00646] [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] [Indexed: 06/25/2024]
Abstract
We describe a detailed investigation into why bulky ligands-those that enable catalysis at "12e -" Pd0-tend to promote overfunctionalization during Pd-catalyzed cross-couplings of dihalogenated substrates. After one cross-coupling event takes place, PdL initially remains coordinated to the π system of the nascent product. Selectivity for mono- vs. difunctionalization arises from the relative rates of π-decomplexation versus a second oxidative addition. Under the Suzuki coupling conditions in this work, direct dissociation of 12e - PdL from the π-complex cannot outcompete oxidative addition. Instead, Pd must be displaced from the π-complex as 14e - PdL(L') by a second incoming ligand L'. The incoming ligand is another molecule of dichloroarene if the reaction conditions do not include π-coordinating solvents or additives. More overfunctionalization tends to result when increased ligand or substrate sterics raises the energy of the bimolecular transition state for separating 14e - PdL(L') from the mono-cross-coupled product. This work has practical implications for optimizing selectivity in cross-couplings involving multiple halogens. For example, we demonstrate that small coordinating additives like DMSO can largely suppress overfunctionalization and that precatalyst structure can also impact selectivity.
Collapse
Affiliation(s)
- Nathaniel G. Larson
- Department of Chemistry and Biochemistry, Montana State University, Bozeman, MT 59717, United States
| | - Jacob P. Norman
- Department of Chemistry and Biochemistry, Montana State University, Bozeman, MT 59717, United States
| | - Sharon R. Neufeldt
- Department of Chemistry and Biochemistry, Montana State University, Bozeman, MT 59717, United States
| |
Collapse
|
21
|
Tho Nguyen M, Gusev DG, Dmitrienko A, Pilkington M, Nikonov GI. Reversible Coupling of Germylone with Isocyanates. Chemistry 2024; 30:e202400613. [PMID: 38379193 DOI: 10.1002/chem.202400613] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2024] [Accepted: 02/20/2024] [Indexed: 02/22/2024]
Abstract
The germylone dimNHCGe (5, dimNHC=diimino N-heterocyclic carbene) undergoes a [2+2] cycloaddition with isocyanates RNCO (R=4-tolyl or 3,5-xylyl) to furnish novel alkyl carboxamido germylenes 7 (R=4-tolyl) and 8 (R=3,5-xylyl), featuring a C-C bond between the former carbene carbon and the isocyanate moiety. Heating a mixture of 8 with 4-tolyl isocyanate to 100 °C results in isocyanate metathesis, demonstrating reversible C-C bond formation on the reduced germanium compound. DFT calculations suggest that this process occurs via the reductive dissociation of isocyanate from 8 that regenerates the parent Ge(0) compound 5.
Collapse
Affiliation(s)
- Minh Tho Nguyen
- Department of Chemistry, Brock University, 1812 Sir Isaac Brock Way, St. Catharines, Ontario L2S 3 A1, Canada
| | - Dmitry G Gusev
- Department of Chemistry and Biochemistry, Wilfrid Laurier University, 75 University Ave W, Waterloo, Ontario N2 L 3 C5, Canada
| | - Anton Dmitrienko
- Department of Chemistry, Brock University, 1812 Sir Isaac Brock Way, St. Catharines, Ontario L2S 3 A1, Canada
| | - Melanie Pilkington
- Department of Chemistry, Brock University, 1812 Sir Isaac Brock Way, St. Catharines, Ontario L2S 3 A1, Canada
| | - Georgii I Nikonov
- Department of Chemistry, Brock University, 1812 Sir Isaac Brock Way, St. Catharines, Ontario L2S 3 A1, Canada
| |
Collapse
|
22
|
Sitkiewicz SP, Ferradás RR, Ramos-Cordoba E, Zaleśny R, Matito E, Luis JM. Spurious Oscillations Caused by Density Functional Approximations: Who is to Blame? Exchange or Correlation? J Chem Theory Comput 2024; 20:3144-3153. [PMID: 38570186 PMCID: PMC11044272 DOI: 10.1021/acs.jctc.3c01339] [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/07/2023] [Revised: 03/18/2024] [Accepted: 03/19/2024] [Indexed: 04/05/2024]
Abstract
We analyze the varying susceptibilities of different density functional approximations (DFAs) to present spurious oscillations on the profiles of several vibrational properties. Among other problems, these spurious oscillations cause significant errors in harmonic and anharmonic IR and Raman frequencies and intensities. This work hinges on a judicious strategy to dissect the exchange and correlation components of DFAs and pinpoint the origins of these oscillations. We identify spurious oscillations in derivatives of all energy components with respect to nuclear displacements, including those energy terms that do not involve numerical integrations. These indirect spurious oscillations are attributed to suboptimal electron densities resulting from a self-consistent field procedure using a DFA that exhibits direct spurious oscillations. Direct oscillations stem from inaccurate numerical integration of the exchange and correlation energy density functionals. A thorough analysis of direct spurious oscillations reveals that only a handful of exchange and correlation components are insensitive to spurious oscillations, giving rise to three families of functionals, BH&H, LSDA, and BLYP. Among the functionals in these families, we encounter four widespread DFAs: BLYP, B3LYP, LC-BLYP, and CAM-B3LYP. Certain DFAs like PBE appear less sensitive to spurious oscillations due to compensatory cancellations between their energy components. Additionally, we found non-negligible but small oscillations in PBE and TPSS, which could be safely employed provided a sufficiently large integration grid is used in the calculations. These findings hint at the key components of current approximations to be improved and emphasize the necessity to develop accurate DFAs suitable for studying molecular spectroscopies.
Collapse
Affiliation(s)
- Sebastian P. Sitkiewicz
- Donostia
International Physics Center (DIPC), Donostia 20018, Euskadi, Spain
- Wrocław
Centre for Networking and Supercomputing, Wrocław University of Science and Technology, Wyb. Wyspiańskiego 27, Wrocław PL-50370, Poland
| | - Rubén R. Ferradás
- Donostia
International Physics Center (DIPC), Donostia 20018, Euskadi, Spain
| | - Eloy Ramos-Cordoba
- Donostia
International Physics Center (DIPC), Donostia 20018, Euskadi, Spain
- Polimero
eta Material Aurreratuak: Fisika, Kimika eta Teknologia, Kimika Fakultatea,
Euskal Herriko Unibertsitatea UPV/EHU, P.K. 1072, Donostia 20080, Euskadi, Spain
- Ikerbasque
Foundation for Science, Plaza Euskadi 5, Bilbao 48009, Euskadi, Spain
- Institute
for Advanced Chemistry of Catalonia (IQAC), CSIC, Jordi Girona 18-26, Barcelona 08034, Spain
| | - Robert Zaleśny
- Faculty
of
Chemistry, Wrocław University of Science
and Technology, Wyb.
Wyspiańskiego 27, Wrocław PL-50370, Poland
| | - Eduard Matito
- Donostia
International Physics Center (DIPC), Donostia 20018, Euskadi, Spain
- Ikerbasque
Foundation for Science, Plaza Euskadi 5, Bilbao 48009, Euskadi, Spain
| | - Josep M. Luis
- Institut
de Química Computacional i Catàlisi (IQCC) and Departament
de Química, Universitat de Girona, Girona 17003, Catalonia, Spain
| |
Collapse
|
23
|
Wensink FJ, Smink CE, Armentrout PB, Bakker JM. IR spectroscopic characterization of 3d transition metal carbene cations, FeCH 2+ and CoCH 2+: periodic trends and a challenge for DFT approaches. Phys Chem Chem Phys 2024; 26:9948-9962. [PMID: 38497938 DOI: 10.1039/d4cp00026a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/19/2024]
Abstract
A combination of IR multiple-photon dissociation (IRMPD) action spectroscopy and quantum chemical calculations was employed to investigate the [M,C,2H]+ (M = Fe and Co) species. These were formed by reacting laser ablated M+ ions with oxirane (ethylene oxide, c-C2H4O) in a room temperature ion trap. IRMPD spectra for the Fe and Co species are very similar and exhibit one major band. Comparison with density functional theory (DFT) and coupled cluster with single and double excitations (CCSD) calculations allows assignment of the spectra to MCH2+ carbene structures. For these 3d transition metal systems, experimental IRMPD spectra compare relatively poorly with DFT calculated IR spectra, but CCSD calculated spectra are a much better match primarily because the M-C stretch gains significant intensity. The origins of this behavior are explored in some detail. The present results are also compared to previous results for the 4d and 5d congeners and the periodic trends in these structures are evaluated.
Collapse
Affiliation(s)
- Frank J Wensink
- Institute for Molecules and Materials, FELIX Laboratory, Radboud University, Toernooiveld 7, 6525 ED Nijmegen, The Netherlands.
| | - Corry E Smink
- Institute for Molecules and Materials, FELIX Laboratory, Radboud University, Toernooiveld 7, 6525 ED Nijmegen, The Netherlands.
| | - P B Armentrout
- Department of Chemistry, University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112, USA.
| | - Joost M Bakker
- Institute for Molecules and Materials, FELIX Laboratory, Radboud University, Toernooiveld 7, 6525 ED Nijmegen, The Netherlands.
| |
Collapse
|
24
|
Villard J, Bircher MP, Rothlisberger U. Structure and dynamics of liquid water from ab initio simulations: adding Minnesota density functionals to Jacob's ladder. Chem Sci 2024; 15:4434-4451. [PMID: 38516095 PMCID: PMC10952088 DOI: 10.1039/d3sc05828j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2023] [Accepted: 02/12/2024] [Indexed: 03/23/2024] Open
Abstract
The accurate representation of the structural and dynamical properties of water is essential for simulating the unique behavior of this ubiquitous solvent. Here we assess the current status of describing liquid water using ab initio molecular dynamics, with a special focus on the performance of all the later generation Minnesota functionals. Findings are contextualized within the current knowledge on DFT for describing bulk water under ambient conditions and compared to experimental data. We find that, contrary to the prevalent idea that local and semilocal functionals overstructure water and underestimate dynamical properties, M06-L, revM06-L, and M11-L understructure water, while MN12-L and MN15-L overdistance water molecules due to weak cohesive effects. This can be attributed to a weakening of the hydrogen bond network, which leads to dynamical fingerprints that are over fast. While most of the hybrid Minnesota functionals (M06, M08-HX, M08-SO, M11, MN12-SX, and MN15) also yield understructured water, their dynamical properties generally improve over their semilocal counterparts. It emerges that exact exchange is a crucial component for accurately describing hydrogen bonds, which ultimately leads to corrections in both the dynamical and structural properties. However, an excessive amount of exact exchange strengthens hydrogen bonds and causes overstructuring and slow dynamics (M06-HF). As a compromise, M06-2X is the best performing Minnesota functional for water, and its D3 corrected variant shows very good structural agreement. From previous studies considering nuclear quantum effects (NQEs), the hybrid revPBE0-D3, and the rung-5 RPA (RPA@PBE) have been identified as the only two approximations that closely agree with experiments. Our results suggest that the M06-2X(-D3) functionals have the potential to further improve the reproduction of experimental properties when incorporating NQEs through path integral approaches. This work provides further proof that accurate modeling of water interactions requires the inclusion of both exact exchange and balanced (non-local) correlation, highlighting the need for higher rungs on Jacob's ladder to achieve predictive simulations of complex biological systems in aqueous environments.
Collapse
Affiliation(s)
- Justin Villard
- Laboratory of Computational Chemistry and Biochemistry, Institute of Chemical Sciences and Engineering, École Polytechnique Fédérale de Lausanne (EPFL) Lausanne CH-1015 Switzerland
| | - Martin P Bircher
- Computational and Soft Matter Physics, Universität Wien Wien A-1090 Austria
| | - Ursula Rothlisberger
- Laboratory of Computational Chemistry and Biochemistry, Institute of Chemical Sciences and Engineering, École Polytechnique Fédérale de Lausanne (EPFL) Lausanne CH-1015 Switzerland
| |
Collapse
|
25
|
Avagliano D, Skreta M, Arellano-Rubach S, Aspuru-Guzik A. DELFI: a computer oracle for recommending density functionals for excited states calculations. Chem Sci 2024; 15:4489-4503. [PMID: 38516092 PMCID: PMC10952086 DOI: 10.1039/d3sc06440a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2023] [Accepted: 02/05/2024] [Indexed: 03/23/2024] Open
Abstract
Density functional theory (DFT) is the workhorse of computational quantum chemistry. One of its main limitations is that choosing the right functional is a non-trivial task left for human experts. The choice is particularly hard for excited state calculations when using its time-dependent formulation (TD-DFT). This is due to the approximations of the method, but also because the photophysical properties of a molecule are defined by a manifold of states that all need to be properly described. This includes not only the relative energy of the states, but also capturing the correct character, order, and intensity of the transitions. In this work, we developed a neural network to recommend functionals to be used on molecules for TD-DFT calculations, by simultaneously considering all these properties for a manifold of states. This was possible by developing a scoring system to define the accuracy of an excited state's calculation against a higher-accuracy reference. The scoring system is generalizable to any level of theory; we here applied it to evaluate the performance of common functionals of different rungs against a higher accuracy method on a large set of organic molecules. The results are collected in a database that we released and made open, providing four million data points to the community for future applications. The scoring system assigns a value between zero and one hundred to each functional for each molecule, transforming the complicated task of learning photophysical properties into a simpler regression task. We used the dataset to train a graph attention neural network to predict the scores for unseen molecules. We call this oracle DELFI (Data-driven EvaLuation of Functionals by Inference), which can be used to quickly screen and predict the ranking of functionals to calculate the optical properties of organic molecules. We validated DELFI in two in silico experiments: choosing a common functional for a series of spiropyran-merocyanine isomers and a unique functional to screen a large dataset of over 50 000 organic photovoltaic molecules, for which an extensive benchmark would be unfeasible. A corresponding web application allows DELFI to be easily run and the results to be analyzed, alleviating the hurdle of choosing the right functional for TD-DFT calculations.
Collapse
Affiliation(s)
- Davide Avagliano
- Department of Chemistry, University of Toronto 80 St. George Street Toronto ON M5S 3H6 Canada
- Department of Computer Science, University of Toronto 40 St. George Street Toronto ON M5S 2E4 Canada
| | - Marta Skreta
- Department of Computer Science, University of Toronto 40 St. George Street Toronto ON M5S 2E4 Canada
- Vector Institute for Artificial Intelligence 661 University Ave. Suite 710 ON M5G 1M1 Toronto Canada
| | | | - Alán Aspuru-Guzik
- Department of Chemistry, University of Toronto 80 St. George Street Toronto ON M5S 3H6 Canada
- Department of Computer Science, University of Toronto 40 St. George Street Toronto ON M5S 2E4 Canada
- Vector Institute for Artificial Intelligence 661 University Ave. Suite 710 ON M5G 1M1 Toronto Canada
- Department of Materials Science & Engineering, University of Toronto 184 College St Toronto M5S 3E4 Canada
- Department of Chemical Engineering & Applied Chemistry, University of Toronto 200 College St ON M5S 3E5 Toronto Canada
- Lebovic Fellow, Canadian Institute for Advanced Research (CIFAR) 66118 University Ave. M5G 1M1 Toronto Canada
- Acceleration Consortium 80 St George St M5S 3H6 Toronto Canada
| |
Collapse
|
26
|
Plett C, Grimme S, Hansen A. Conformational energies of biomolecules in solution: Extending the MPCONF196 benchmark with explicit water molecules. J Comput Chem 2024; 45:419-429. [PMID: 37982322 DOI: 10.1002/jcc.27248] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2023] [Revised: 10/16/2023] [Accepted: 10/17/2023] [Indexed: 11/21/2023]
Abstract
A prerequisite for the computational prediction of molecular properties like conformational energies of biomolecules is a reliable, robust, and computationally affordable method usually selected according to its performance for relevant benchmark sets. However, most of these sets comprise molecules in the gas phase and do not cover interactions with a solvent, even though biomolecules typically occur in aqueous solution. To address this issue, we introduce a with explicit water molecules solvated version of a gas-phase benchmark set containing 196 conformers of 13 peptides and other relevant macrocycles, namely MPCONF196 [J. Řezáč et al., JCTC 2018, 14, 1254-1266], and provide very accurate PNO-LCCSD(T)-F12b/AVQZ' reference values. The novel solvMPCONF196 benchmark set features two additional challenges beyond the description of conformers in the gas phase: conformer-water and water-water interactions. The overall best performing method for this set is the double hybrid revDSDPBEP86-D4/def2-QZVPP yielding conformational energies of almost coupled cluster quality. Furthermore, some (meta-)GGAs and hybrid functionals like B97M-V and ω B97M-D with a large basis set reproduce the coupled cluster reference with an MAD below 1 kcal mol- 1 . If more efficient methods are required, the composite DFT-method r2 SCAN-3c (MAD of 1.2 kcal mol- 1 ) is a good alternative, and when conformational energies of polypeptides or macrocycles with more than 500-1000 atoms are in the focus, the semi-empirical GFN2-xTB or the MMFF94 force field (for very large systems) are recommended.
Collapse
Affiliation(s)
- Christoph Plett
- Mulliken Center for Theoretical Chemistry, Clausius-Institut für Physikalische und Theoretische Chemie, Universität Bonn, Bonn, Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Clausius-Institut für Physikalische und Theoretische Chemie, Universität Bonn, Bonn, Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Clausius-Institut für Physikalische und Theoretische Chemie, Universität Bonn, Bonn, Germany
| |
Collapse
|
27
|
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
|
28
|
Wang T, Zhang Z, Jiang S, Yan W, Li S, Zhuang J, Xie H, Li G, Jiang L. Spectroscopic characterization of carbon monoxide activation by neutral chromium carbides. Phys Chem Chem Phys 2024; 26:5962-5968. [PMID: 38293768 DOI: 10.1039/d4cp00011k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2024]
Abstract
Spectroscopic characterization of carbon monoxide activation by neutral metal carbides is of essential importance for understanding the structure-reactivity relationships of catalytic sites, but has been proven to be very challenging owing to the difficulty in size selection. Here, we report a size-specific infrared-vacuum ultraviolet spectroscopic study of the reactions between carbon monoxide with neutral chromium carbides. Quantum chemical calculations were carried out to identify the low-lying structures and to interpret the experimental features. The results reveal that the most stable structure of CrC3(CO)2 consists of a CCO ketenylidene unit and that of CrC4(CO)2 has a semi-bridging CO with a very low CO stretching vibrational frequency at 1821 cm-1. The electron structure analyses show that this semi-bridging CO is highly activated through the delocalized Cr-C-C three-center two-electron (3c-2e) interaction between the antibonding orbitals of CO and the metal carbide skeleton. The formation of these metal carbide carbonyls is found to be both thermodynamically exothermic and kinetically facile in the gas phase. The present findings have important implications for the mechanical understanding of the catalytic processes with isolated metal atoms/clusters dispersed on supports.
Collapse
Affiliation(s)
- Tiantong Wang
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China.
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Zhaoyan Zhang
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China.
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Shuai Jiang
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China.
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Wenhui Yan
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China.
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Shangdong Li
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China.
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Jianxing Zhuang
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China.
- University of Chinese Academy of Sciences, Beijing 100049, China
| | - Hua Xie
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China.
| | - Gang Li
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China.
| | - Ling Jiang
- State Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China.
- Hefei National Laboratory, Hefei 230088, China
| |
Collapse
|
29
|
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
|
30
|
Fang L, Gou G, Wang M, Fan T, Yin Y, Li L. Regulating the Flexibility to Assemble Porous Single-Atom Fe-Coordinated Metallopolymers for Efficient Heterogeneous Catalytic Oxidations. ACS APPLIED MATERIALS & INTERFACES 2024; 16:5823-5833. [PMID: 38285621 DOI: 10.1021/acsami.3c15958] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/31/2024]
Abstract
Metallopolymers as organic-inorganic hybrid materials formulated by metal embedding organic polymers show great potential for novel heterogeneous catalysis, in terms of the facile structural design and tunability. Herein, the disadvantage of nonporous stacking of one-dimensional (1D) structures has been suppressed by chain modulation of the 1D metallopolymers, allowing for the convenient construction of porous assemblies with single-atom dispersion and accessible active sites. By postmodification, the Fe/CM-1 catalyst readily synthesized by coordinating the Fe(II) to the twisted chain of 1D Schiff-base polymer possesses expedient flexibility, showing the highest porosity, remarkable heterogeneous recyclability, and thus prominent catalytic activity for the selective oxidation of benzylamine and alcohols. Moreover, control experiments supported by computational studies demonstrated that the unique pincer structure of Fe/CM-1 effectively maintains the valence state of the anchored single-atom iron, facilitating single-electron transfer and promoting efficient iron redox cycling during the catalytic process. Notably, these 1D metallopolymers have the advantage of cost-effectiveness, easy preparation in gram-scale, and utilization in continuous reaction, providing inspirations for facile synthesis of efficient heterogeneous catalysts from the well-developed 1D metallopolymers.
Collapse
Affiliation(s)
- Lei Fang
- Shanghai Key Lab of Chemical Assessment and Sustainability, School of Chemical Science and Engineering, Tongji University, 1239 Siping Road, Shanghai 200092, China
| | - Gaozhang Gou
- Shanghai Key Lab of Chemical Assessment and Sustainability, School of Chemical Science and Engineering, Tongji University, 1239 Siping Road, Shanghai 200092, China
| | - Man Wang
- Shanghai Key Lab of Chemical Assessment and Sustainability, School of Chemical Science and Engineering, Tongji University, 1239 Siping Road, Shanghai 200092, China
| | - Tao Fan
- Shanghai Key Lab of Chemical Assessment and Sustainability, School of Chemical Science and Engineering, Tongji University, 1239 Siping Road, Shanghai 200092, China
| | - Ying Yin
- Shanghai Key Lab of Chemical Assessment and Sustainability, School of Chemical Science and Engineering, Tongji University, 1239 Siping Road, Shanghai 200092, China
| | - Liangchun Li
- Shanghai Key Lab of Chemical Assessment and Sustainability, School of Chemical Science and Engineering, Tongji University, 1239 Siping Road, Shanghai 200092, China
| |
Collapse
|
31
|
Franzke YJ, Bruder F, Gillhuber S, Holzer C, Weigend F. Paramagnetic Nuclear Magnetic Resonance Shifts for Triplet Systems and Beyond with Modern Relativistic Density Functional Methods. J Phys Chem A 2024; 128:670-686. [PMID: 38195394 DOI: 10.1021/acs.jpca.3c07093] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2024]
Abstract
An efficient framework for the calculation of paramagnetic NMR (pNMR) shifts within exact two-component (X2C) theory and (current-dependent) density functional theory (DFT) up to the class of local hybrid functionals (LHFs) is presented. Generally, pNMR shifts for systems with more than one unpaired electron depend on the orbital shielding contribution and a temperature-dependent term. The latter includes zero-field splitting (ZFS), hyperfine coupling (HFC), and the g-tensor. For consistency, we calculate these three tensors at the same level of theory, i.e., using scalar-relativistic X2C augmented with spin-orbit perturbation theory. Results for pNMR chemical shifts of transition-metal complexes reveal that this X2C-DFT framework can yield good results for both the shifts and the individual tensor contributions of metallocenes and related systems, especially if the HFC constant is large. For small HFC constants, the relative error is often large, and sometimes the sign may be off. 4d and 5d complexes with more complicated structures demonstrate the limitations of a fully DFT-based approach. Additionally, a Co-based complex with a very large ZFS and pronounced multireference character is not well described. Here, a hybrid DFT-multireference framework is necessary for accurate results. Our results show that X2C is sufficient to describe relativistic effects and computationally cheaper than a fully relativistic approach. Thus, it allows use of large basis sets for converged HFCs. Overall, current-dependent meta-generalized gradient approximations and LHFs show some potential; however, the currently available functionals leave a lot to be desired, and the predictive power is limited.
Collapse
Affiliation(s)
- Yannick J Franzke
- Fachbereich Chemie, Philipps-Universität Marburg, Hans-Meerwein-Straße 4, 35032 Marburg, Germany
| | - Florian Bruder
- Fachbereich Chemie, Philipps-Universität Marburg, Hans-Meerwein-Straße 4, 35032 Marburg, Germany
| | - Sebastian Gillhuber
- Institute of Inorganic Chemistry, Karlsruhe Institute of Technology (KIT), Engesserstraße 15, 76131 Karlsruhe, Germany
| | - Christof Holzer
- Institute of Theoretical Solid State Physics, Karlsruhe Institute of Technology (KIT), Wolfgang-Gaede-Straße 1, 76131 Karlsruhe, Germany
| | - Florian Weigend
- Fachbereich Chemie, Philipps-Universität Marburg, Hans-Meerwein-Straße 4, 35032 Marburg, Germany
| |
Collapse
|
32
|
Wappett DA, Goerigk L. Exploring CPS-Extrapolated DLPNO-CCSD(T 1) Reference Values for Benchmarking DFT Methods on Enzymatically Catalyzed Reactions. J Phys Chem A 2024; 128:62-72. [PMID: 38124376 DOI: 10.1021/acs.jpca.3c05086] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2023]
Abstract
Domain-based local pair natural orbital coupled-cluster singles doubles with perturbative triples [DLPNO-CCSD(T)] is regularly used to calculate reliable benchmark reference values at a computational cost significantly lower than that of canonical CCSD(T). Recent work has shown that even greater accuracy can be obtained at only a small additional cost through extrapolation to the complete PNO space (CPS) limit. Herein, we test two levels of CPS extrapolation, CPS(5,6), which approximates the accuracy of standard TightPNO, and CPS(6,7), which surpasses it, as benchmark values to test density functional approximations (DFAs) on a small set of organic and transition-metal-dependent enzyme active site models. Between the different reference levels of theory, there are changes in the magnitudes of the absolute deviations for all functionals, but these are small and there is minimal impact on the relative rankings of the tested DFAs. The differences are more significant for the metalloenzymes than the organic enzymes, so we repeat the tests on our entire ENZYMES22 set of organic enzyme active site models [Wappett, D. A.; Goerigk, L. J. Phys. Chem. A 2019, 123, 7057-7074] to confirm that using the CPS extrapolations for the reference values has negligible impact on the benchmarking outcomes. This means that we can particularly recommend CPS(5,6) as an alternative to standard TightPNO settings for calculating reference values, increasing the applicability of DLPNO-CCSD(T) in benchmarking reaction energies and barrier heights of larger models of organic enzymes. DLPNO-CCSD(T1)/CPS(6,7) energies for ENZYMES22 are finally presented as updated reference values for the set, reflecting the recent improvements in the method.
Collapse
Affiliation(s)
- Dominique A Wappett
- School of Chemistry, The University of Melbourne, Parkville, Victoria 3010, Australia
| | - Lars Goerigk
- School of Chemistry, The University of Melbourne, Parkville, Victoria 3010, Australia
| |
Collapse
|
33
|
Schattenberg C, Wodyński A, Åström H, Sundholm D, Kaupp M, Lehtola S. Revisiting Gauge-Independent Kinetic Energy Densities in Meta-GGAs and Local Hybrid Calculations of Magnetizabilities. J Phys Chem A 2023; 127:10896-10907. [PMID: 38100678 PMCID: PMC10758120 DOI: 10.1021/acs.jpca.3c06244] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2023] [Revised: 10/27/2023] [Accepted: 11/08/2023] [Indexed: 12/17/2023]
Abstract
In a recent study [J. Chem. Theory Comput. 2021, 17, 1457-1468], some of us examined the accuracy of magnetizabilities calculated with density functionals representing the local density approximation (LDA), generalized gradient approximation (GGA), meta-GGA (mGGA), as well as global hybrid (GH) and range-separated (RS) hybrid functionals by assessment against accurate reference values obtained with coupled-cluster theory with singles, doubles, and perturbative triples [CCSD(T)]. Our study was later extended to local hybrid (LH) functionals by Holzer et al. [J. Chem. Theory Comput. 2021, 17, 2928-2947]; in this work, we examine a larger selection of LH functionals, also including range-separated LH (RSLH) functionals and strong-correlation LH (scLH) functionals. Holzer et al. also studied the importance of the physically correct handling of the magnetic gauge dependence of the kinetic energy density (τ) in mGGA calculations by comparing the Maximoff-Scuseria formulation of τ used in our aforementioned study to the more physical current-density extension derived by Dobson. In this work, we also revisit this comparison with a larger selection of mGGA functionals. We find that the newly tested LH, RSLH, and scLH functionals outperform all of the functionals considered in the previous studies. The various LH functionals afford the seven lowest mean absolute errors while also showing remarkably small standard deviations and mean errors. Most strikingly, the best two functionals are scLHs that also perform remarkably well in cases with significant multiconfigurational character, such as the ozone molecule, which is traditionally excluded from statistical error evaluations due to its large errors with common density functionals.
Collapse
Affiliation(s)
- Caspar
J. Schattenberg
- Institut
für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, D-10623 Berlin, Germany
| | - Artur Wodyński
- Institut
für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, D-10623 Berlin, Germany
| | - Hugo Åström
- Department
of Chemistry, Faculty of Science, University
of Helsinki, P.O. Box 55
(A.I. Virtanens plats 1), University of Helsinki FI-00014, Finland
| | - Dage Sundholm
- Department
of Chemistry, Faculty of Science, University
of Helsinki, P.O. Box 55
(A.I. Virtanens plats 1), University of Helsinki FI-00014, Finland
| | - Martin Kaupp
- Institut
für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, D-10623 Berlin, Germany
| | - Susi Lehtola
- Department
of Chemistry, Faculty of Science, University
of Helsinki, P.O. Box 55
(A.I. Virtanens plats 1), University of Helsinki FI-00014, Finland
- Molecular
Sciences Software Institute, Blacksburg, Virginia 24061, United States
| |
Collapse
|
34
|
Kanchanakungwankul S, Verma P, Janesko BG, Scalmani G, Frisch MJ, Truhlar DG. M11pz: A Nonlocal Meta Functional with Zero Hartree-Fock Exchange and with Broad Accuracy for Chemical Energies and Structures. J Chem Theory Comput 2023; 19:9102-9117. [PMID: 38096578 DOI: 10.1021/acs.jctc.2c01315] [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
The accuracy of Kohn-Sham density functional theory depends strongly on the approximation to the exchange-correlation functional. In this work, we present a new exchange-correlation functional called M11pz (M11 plus rung-3.5 terms with zero Hartree-Fock exchange) that is built on the M11plus functional with the goal of using its rung-3.5 terms without a Hartree-Fock exchange term, especially to improve the accuracy for strongly correlated systems. The M11pz functional is optimized with the same local and rung-3.5 ingredients that are used in M11plus but without any percentage of Hartree-Fock exchange. The performance of M11pz is compared with eight local functionals, and M11pz is found to be in top three when the errors or ranks are averaged over eight grouped and partially overlapping databases: AME418/22, atomic and molecular energies; MGBE172, main-group bond energies; TMBE40, transition-metal bond energies; SR309, single-reference systems; MR54, multireference systems; BH192, barrier heights; NC579, noncovalent interaction energies; and MS20, molecular structures. For calculations of band gaps of solids, M11pz is the second best of the nine tested functionals that have zero Hartree-Fock exchange.
Collapse
Affiliation(s)
- Siriluk Kanchanakungwankul
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Pragya Verma
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Benjamin G Janesko
- Department of Chemistry and Biochemistry, Texas Christian University, Fort Worth, Texas 76110, United States
| | | | | | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
35
|
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
|
36
|
Kour M, Taborosi A, Boyd ES, Szilagyi RK. Development of molecular cluster models to probe pyrite surface reactivity. J Comput Chem 2023; 44:2486-2500. [PMID: 37650712 DOI: 10.1002/jcc.27213] [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: 04/02/2023] [Revised: 07/28/2023] [Accepted: 08/09/2023] [Indexed: 09/01/2023]
Abstract
The recent discovery that anaerobic methanogens can reductively dissolve pyrite and utilize dissolution products as a source of iron and sulfur to meet their biosynthetic demands for these elements prompted the development of atomic-scale nanoparticle models, as maquettes of reactive surface sites, for describing the fundamental redox steps that take place at the mineral surface during reduction. The given report describes our computational approach for modeling n(FeS2 ) nanoparticles originated from mineral bulk structure. These maquettes contain a comprehensive set of coordinatively unsaturated Fe(II) sites that are connected via a range of persulfide (S2 2- ) ligation. In addition to the specific maquettes with n = 8, 18, and 32 FeS2 units, we established guidelines for obtaining low-energy structures by considering the pattern of ionic, covalent, and magnetic interactions among the metal and ligand sites. The developed models serve as computational nano-reactors that can be used to describe the reductive dissolution mechanism of pyrite to better understand the reactive sites on the mineral, where microbial extracellular electron-transfer reactions can occur.
Collapse
Affiliation(s)
- Manjinder Kour
- Department of Microbiology and Cell Biology, Montana State University, Bozeman, Montana, USA
| | - Attila Taborosi
- Research Initiative for Supra-Materials, Faculty of Engineering, Shinshu University, Nagano, Japan
| | - Eric S Boyd
- Department of Microbiology and Cell Biology, Montana State University, Bozeman, Montana, USA
| | - Robert K Szilagyi
- Department of Chemistry, The University of British Columbia, Okanagan, Kelowna, British Columbia, Canada
| |
Collapse
|
37
|
Gould T. A step toward density benchmarking-The energy-relevant "mean field error". J Chem Phys 2023; 159:204111. [PMID: 38018751 DOI: 10.1063/5.0175925] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2023] [Accepted: 11/05/2023] [Indexed: 11/30/2023] Open
Abstract
Since the development of generalized gradient approximations in the 1990s, approximations based on density functional theory have dominated electronic structure theory calculations. Modern approximations can yield energy differences that are precise enough to be predictive in many instances, as validated by large- and small-scale benchmarking efforts. However, assessing the quality of densities has been the subject of far less attention, in part because reliable error measures are difficult to define. To this end, this work introduces the mean-field error, which directly assesses the quality of densities from approximations. The mean-field error is contextualized within existing frameworks of density functional error analysis and understanding and shown to be part of the density-driven error. It is demonstrated in several illustrative examples. Its potential use in future benchmarking protocols is discussed, and some conclusions are drawn.
Collapse
Affiliation(s)
- Tim Gould
- Qld Micro- and Nanotechnology Centre, Griffith University, Nathan, Qld 4111, Australia
| |
Collapse
|
38
|
Wappett D, Goerigk L. Benchmarking Density Functional Theory Methods for Metalloenzyme Reactions: The Introduction of the MME55 Set. J Chem Theory Comput 2023; 19:8365-8383. [PMID: 37943578 PMCID: PMC10688432 DOI: 10.1021/acs.jctc.3c00558] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2023] [Revised: 10/13/2023] [Accepted: 10/16/2023] [Indexed: 11/10/2023]
Abstract
We present a new benchmark set of metalloenzyme model reaction energies and barrier heights that we call MME55. The set contains 10 different enzymes, representing eight transition metals, both open and closed shell systems, and system sizes of up to 116 atoms. We use four DLPNO-CCSD(T)-based approaches to calculate reference values against which we then benchmark the performance of a range of density functional approximations with and without dispersion corrections. Dispersion corrections improve the results across the board, and triple-ζ basis sets provide the best balance of efficiency and accuracy. Jacob's ladder is reproduced for the whole set based on averaged mean absolute (percent) deviations, with the double hybrids SOS0-PBE0-2-D3(BJ) and revDOD-PBEP86-D4 standing out as the most accurate methods for the MME55 set. The range-separated hybrids ωB97M-V and ωB97X-V also perform well here and can be recommended as a reliable compromise between accuracy and efficiency; they have already been shown to be robust across many other types of chemical problems, as well. Despite the popularity of B3LYP in computational enzymology, it is not a strong performer on our benchmark set, and we discourage its use for enzyme energetics.
Collapse
Affiliation(s)
- Dominique
A. Wappett
- School of Chemistry, The University
of Melbourne, Melbourne, Victoria 3010, Australia
| | - Lars Goerigk
- School of Chemistry, The University
of Melbourne, Melbourne, Victoria 3010, Australia
| |
Collapse
|
39
|
Brothers EN, Bengali AA, Scalmani G, Janesko BG, Verma P, Truhlar DG, Frisch MJ. Comparing Density Functional Theory Metal-Ligand Bond Dissociation Enthalpies with Experimental Solution-Phase Enthalpies of Activation for Bond Dissociation. J Phys Chem A 2023; 127:9695-9704. [PMID: 37939355 DOI: 10.1021/acs.jpca.3c04838] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2023]
Abstract
The predictive ability of density functional theory is fundamental to its usefulness in chemical applications. Recent work has compared solution-phase enthalpies of activation for metal-ligand bond dissociation to enthalpies of reaction for bond dissociation, and the present work continues those comparisons for 43 density functional methods. The results for ligand dissociation enthalpies of 30 metal-ligand complexes tested in this work reveal significant inadequacies of some functionals as well as challenges from the dispersion corrections to some functionals. The analysis presented here demonstrates the excellent performance of a recent density functional, M11plus, which contains nonlocal rung-3.5 correlation. We also find a good agreement between theory and experiment for some functionals without empirical dispersion corrections such as M06, r2SCAN, M06-L, and revM11, as well as good performance for some functionals with added dispersion corrections such as ωB97X-D (which always has a correction) and BLYP, B3LYP, CAM-B3LYP, and PBE0 when the optional dispersion corrections are added.
Collapse
Affiliation(s)
- Edward N Brothers
- Gaussian, Inc., 340 Quinnipiac Street, Wallingford, Connecticut 06492, United States
| | - Ashfaq A Bengali
- Division of Arts and Sciences, Texas A&M University at Qatar, Doha, Qatar
| | - Giovanni Scalmani
- Gaussian, Inc., 340 Quinnipiac Street, Wallingford, Connecticut 06492, United States
| | - Benjamin G Janesko
- Department of Chemistry and Biochemistry, Texas Christian University, Fort Worth, Texas 76110, United States
| | - Pragya Verma
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Michael J Frisch
- Gaussian, Inc., 340 Quinnipiac Street, Wallingford, Connecticut 06492, United States
| |
Collapse
|
40
|
Yao YR, Zhao J, Meng Q, Hu HS, Guo M, Yan Y, Zhuang J, Yang S, Fortier S, Echegoyen L, Schwarz WHE, Li J, Chen N. Synthesis and Characterization of U≡C Triple Bonds in Fullerene Compounds. J Am Chem Soc 2023; 145:25440-25449. [PMID: 37955678 DOI: 10.1021/jacs.3c10042] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2023]
Abstract
Despite decades of efforts, the actinide-carbon triple bond has remained an elusive target, defying synthesis in any isolable compound. Herein, we report the successful synthesis of uranium-carbon triple bonds in carbide-bridged bimetallic [U≡C-Ce] units encapsulated inside the fullerene cages of C72 and C78. The molecular structures of UCCe@C2n and the nature of the U≡C triple bond were characterized through X-ray crystallography and various spectroscopic analyses, revealing very short uranium-carbon bonds of 1.921(6) and 1.930(6) Å, with the metals existing in their highest oxidation states of +6 and +4 for uranium and cerium, respectively. Quantum-chemical studies further demonstrate that the C2n cages are crucial for stabilizing the [UVI≡C-CeIV] units through covalent and coordinative interactions. This work offers a new fundamental understanding of the elusive uranium-carbon triple bond and informs the design of complexes with similar bonding motifs, opening up new possibilities for creating distinctive molecular compounds and materials.
Collapse
Affiliation(s)
- Yang-Rong Yao
- College of Chemistry, Chemical Engineering and Materials Science & State Key Laboratory of Radiation Medicine and Protection, Soochow University, Suzhou 215123, China
| | - Jing Zhao
- Department of Chemistry and Engineering Research Center of Advanced Rare-Earth Materials of the Ministry of Education, Tsinghua University, Beijing 100084, China
| | - Qingyu Meng
- College of Chemistry, Chemical Engineering and Materials Science & State Key Laboratory of Radiation Medicine and Protection, Soochow University, Suzhou 215123, China
| | - Han-Shi Hu
- Department of Chemistry and Engineering Research Center of Advanced Rare-Earth Materials of the Ministry of Education, Tsinghua University, Beijing 100084, China
| | - Min Guo
- College of Chemistry, Chemical Engineering and Materials Science & State Key Laboratory of Radiation Medicine and Protection, Soochow University, Suzhou 215123, China
| | - Yingjing Yan
- College of Chemistry, Chemical Engineering and Materials Science & State Key Laboratory of Radiation Medicine and Protection, Soochow University, Suzhou 215123, China
| | - Jiaxin Zhuang
- College of Chemistry, Chemical Engineering and Materials Science & State Key Laboratory of Radiation Medicine and Protection, Soochow University, Suzhou 215123, China
| | - Shangfeng Yang
- Department of Materials Science and Engineering, CAS Key Laboratory of Materials for Energy Conversion, Anhui Laboratory of Advanced Photon Science and Technology, University of Science and Technology of China, Hefei 230026, China
| | - Skye Fortier
- Department of Chemistry and Biochemistry, University of Texas at El Paso, El Paso, Texas 79968, United States
| | - Luis Echegoyen
- Institut Catalá d'Investigació Química, Ave. Països Catalans 16, 43007 Tarragona, Spain
- Department of Chemistry and Biochemistry, University of Texas at El Paso, El Paso, Texas 79968, United States
| | - W H Eugen Schwarz
- Department of Chemistry and Engineering Research Center of Advanced Rare-Earth Materials of the Ministry of Education, Tsinghua University, Beijing 100084, China
- Physikalische und Theoretische Chemie, Universität Siegen, Siegen 57068, Germany
| | - Jun Li
- Department of Chemistry and Engineering Research Center of Advanced Rare-Earth Materials of the Ministry of Education, Tsinghua University, Beijing 100084, China
- Department of Chemistry and Guangdong Provincial Key Laboratory of Catalytic Chemistry, Southern University of Science and Technology, Shenzhen 518055, China
| | - Ning Chen
- College of Chemistry, Chemical Engineering and Materials Science & State Key Laboratory of Radiation Medicine and Protection, Soochow University, Suzhou 215123, China
| |
Collapse
|
41
|
Jaroš A, Straka M. Unraveling actinide-actinide bonding in fullerene cages: a DFT versus ab initio methodological study. Phys Chem Chem Phys 2023; 25:31500-31513. [PMID: 37962545 DOI: 10.1039/d3cp03606e] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2023]
Abstract
Actinide-actinide bonding poses a challenge for both experimental and theoretical chemists because of both the scarcity of experimental data and the exotic nature of actinide bonding due to the involvement and mixing of actinide 7s-, 6p-, 6d-, and particularly 5f-orbitals. Only a few experimental examples of An-An bonding have been reported so far. Here, we perform a methodological study of actinide-actinide bonding on experimentally known Th2@C80 and U2@C80 systems. We compared selected GGA, meta-GGA, hybrid-GGA and range-separated hybrid-GGA functionals with the results obtained using a multireference CASPT2 method, which we consider as a reference point. We show that functionals such as BP86, PBE or TPSS perform well for predicting geometries, while range-separated hybrids are superior in the description of the chemical bonding. None of the tested functionals were deemed reliable regarding the correct electronic spin ground state. Based on the results of this methodological study, we re-evaluate selected previously studied diactinide fullerene systems using more reliable protocol.
Collapse
Affiliation(s)
- Adam Jaroš
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo nám. 2, CZ-16610, Prague, Czech Republic.
- Faculty of Science, Charles University, Albertov 2038/6, Prague 2, 128 43, Czech Republic
| | - Michal Straka
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo nám. 2, CZ-16610, Prague, Czech Republic.
| |
Collapse
|
42
|
Boychuk BTA, Meyer SP, Wetmore SD. Generation of an accurate CCSD(T)/CBS data set and assessment of DFT methods for the binding strengths of group I metal-nucleic acid complexes. Front Chem 2023; 11:1296787. [PMID: 38053674 PMCID: PMC10694745 DOI: 10.3389/fchem.2023.1296787] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2023] [Accepted: 10/26/2023] [Indexed: 12/07/2023] Open
Abstract
Accurate information about interactions between group I metals and nucleic acids is required to understand the roles these metals play in basic cellular functions, disease progression, and pharmaceuticals, as well as to aid the design of new energy storage materials and nucleic acid sensors that target metal contaminants, among other applications. From this perspective, this work generates a complete CCSD(T)/CBS data set of the binding energies for 64 complexes involving each group I metal (Li+, Na+, K+, Rb+, or Cs+) directly coordinated to various sites in each nucleic acid component (A, C, G, T, U, or dimethylphosphate). This data have otherwise been challenging to determine experimentally, with highly accurate information missing for many group I metal-nucleic acid combinations and no data available for the (charged) phosphate moiety. Subsequently, the performance of 61 DFT methods in combination with def2-TZVPP is tested against the newly generated CCSD(T)/CBS reference values. Detailed analysis of the results reveals that functional performance is dependent on the identity of the metal (with increased errors as group I is descended) and nucleic acid binding site (with larger errors for select purine coordination sites). Over all complexes considered, the best methods include the mPW2-PLYP double-hybrid and ωB97M-V RSH functionals (≤1.6% MPE; <1.0 kcal/mol MUE). If more computationally efficient approaches are required, the TPSS and revTPSS local meta-GGA functionals are reasonable alternatives (≤2.0% MPE; <1.0 kcal/mol MUE). Inclusion of counterpoise corrections to account for basis set superposition error only marginally improves the computed binding energies, suggesting that these corrections can be neglected with little loss in accuracy when using larger models that are necessary for describing biosystems and biomaterials. Overall, the most accurate functionals identified in this study will permit future works geared towards uncovering the impact of group I metals on the environment and human biology, designing new ways to selectively sense harmful metals, engineering modern biomaterials, and developing improved computational methods to more broadly study group I metal-nucleic acid interactions.
Collapse
Affiliation(s)
| | | | - Stacey D. Wetmore
- Department of Chemistry and Biochemistry, University of Lethbridge, Lethbridge, AB, Canada
| |
Collapse
|
43
|
Segizbayev M, Tho Nguyen M, Gusev DG, Dmitrienko A, Pilkington M, van der Est A, Nikonov GI. A Guanidine-Supported π-Complex of Germanium Amenable to Intramolecular C-C Cleavage in Arene and Ge Atom Transfer. Chemistry 2023; 29:e202301981. [PMID: 37732936 DOI: 10.1002/chem.202301981] [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: 08/17/2023] [Revised: 09/20/2023] [Accepted: 09/21/2023] [Indexed: 09/22/2023]
Abstract
The germylone dimNHCGe (dimNHC=diimino N-heterocyclic carbene) reacts with azides N3 R (R=SiMe3 or p-tolyl) to furnish the first examples of germanium π-complexes, i. e. guanidine-ligated compounds (dimNHI-SiMe3 )Ge (NHI=N-heterocyclic imine, R=SiMe3 ) and (dimNHI-Tol)Ge (R=p-tolyl). DFT calculations suggest that these species are formed by a Staudinger type replacement of dinitrogen in the azide by a nucleophilic germylone, leading to a transient carbene adduct of iminogermylidene. Heating a solution of compound (dimNHI-SiMe3 )Ge to 70 °C results in extrusion of the iminogermylidene that further aggregates to produce the known [Me3 SiNGe]4 tetramer, whereas the imidazolylidene fragment transforms into an unusual heptatriene species that can be considered as a product of carbene insertion into the C-C bond of a pendant Ar substituent at the imidazolylidene nitrogen of the dimNHC. Reaction of (dimNHI-SiMe3 )Ge with tetrachloro-o-benzoquinone results in the net transfer of a germanium atom and formation of the free diimino-guanidine ligand. This ligand also forms when (dimNHI-SiMe3 )Ge is treated with azide N3 (p-Tol), with the germanium product being [(p-Tol)NGe]n.
Collapse
Affiliation(s)
- Medet Segizbayev
- Department of Chemistry, Brock University, 1812 Sir Isaac Brock Way, St. Catharines, Ontario L2S 3A1, Canada
| | - Minh Tho Nguyen
- Department of Chemistry, Brock University, 1812 Sir Isaac Brock Way, St. Catharines, Ontario L2S 3A1, Canada
| | - Dmitry G Gusev
- Department of Chemistry and Biochemistry, Wilfrid Laurier University, 75 University Ave W, Waterloo, Ontario N2 L 3 C5, Canada
| | - Anton Dmitrienko
- Department of Chemistry, Brock University, 1812 Sir Isaac Brock Way, St. Catharines, Ontario L2S 3A1, Canada
| | - Melanie Pilkington
- Department of Chemistry, Brock University, 1812 Sir Isaac Brock Way, St. Catharines, Ontario L2S 3A1, Canada
| | - Art van der Est
- Department of Chemistry, Brock University, 1812 Sir Isaac Brock Way, St. Catharines, Ontario L2S 3A1, Canada
| | - Georgii I Nikonov
- Department of Chemistry, Brock University, 1812 Sir Isaac Brock Way, St. Catharines, Ontario L2S 3A1, Canada
| |
Collapse
|
44
|
Ai W, Su NQ, Fang WH. Short-range screened density matrix functional for proper descriptions of thermochemistry, thermochemical kinetics, nonbonded interactions, and singlet diradicals. J Chem Phys 2023; 159:174110. [PMID: 37933778 DOI: 10.1063/5.0169234] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2023] [Accepted: 10/09/2023] [Indexed: 11/08/2023] Open
Abstract
Common one-electron reduced density matrix (1-RDM) functionals that depend on Coulomb and exchange-only integrals tend to underestimate dynamic correlation, preventing reduced density matrix functional theory (RDMFT) from achieving comparable accuracy to density functional theory in main-group thermochemistry and thermochemical kinetics. The recently developed ωP22 functional introduces a semi-local density functional to screen the erroneous short-range portion of 1-RDM functionals without double-counting correlation, potentially providing a better treatment of dynamic correlation around equilibrium geometries. Herein, we systematically evaluate the performance of this functional model, which consists of two parameters, on main-group thermochemistry, thermochemical kinetics, nonbonded interactions, and more. Tests on atomization energies, vibrational frequencies, and reaction barriers reveal that the ωP22 functional model can reliably predict properties at equilibrium and slightly away from equilibrium geometries. In particular, it outperforms commonly used density functionals in the prediction of reaction barriers, nonbonded interactions, and singlet diradicals, thus enhancing the predictive power of RDMFT for routine calculations of thermochemistry and thermochemical kinetics around equilibrium geometries. Further development is needed in the future to refine short- and long-range approximations in the functional model in order to achieve an excellent description of properties both near and far from equilibrium geometries.
Collapse
Affiliation(s)
- Wenna Ai
- Department of Chemistry, Key Laboratory of Advanced Energy Materials Chemistry (Ministry of Education), Frontiers Science Center for New Organic Matter, Nankai University, Tianjin 300071, China
| | - Neil Qiang Su
- Department of Chemistry, Key Laboratory of Advanced Energy Materials Chemistry (Ministry of Education), Frontiers Science Center for New Organic Matter, Nankai University, Tianjin 300071, China
| | - Wei-Hai Fang
- Department of Chemistry, Key Laboratory of Advanced Energy Materials Chemistry (Ministry of Education), Frontiers Science Center for New Organic Matter, Nankai University, Tianjin 300071, China
| |
Collapse
|
45
|
Kuribara T, Kaneki A, Mihara Y, Yuba H, Nemoto T. Scandium-Catalyzed Phenol-Directed Construction of 5-Carbonyl-4-hydroxybenzofurans via Intramolecular Friedel-Crafts Reaction. Org Lett 2023; 25:7890-7894. [PMID: 37882510 DOI: 10.1021/acs.orglett.3c03123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2023]
Abstract
Herein we report a scandium-catalyzed regioselective synthesis of 5-carbonyl-4-hydroxybenzofurans via a phenol-directed intramolecular Friedel-Crafts reaction. This synthetic method was applied for the total synthesis of furanoflavones. Experimental studies and density functional theory calculations suggest that hydrogen bond interactions between the phenolic hydroxy group and the scandium complex realize regioselective intramolecular cyclization.
Collapse
Affiliation(s)
- Takahito Kuribara
- Institute for Advanced Academic Research, Chiba University, 1-33, Yayoi-cho, Chiba 263-8522, Japan
- Graduate School of Pharmaceutical Sciences, Chiba University, 1-8-1, Inohana, Chuo-ku, Chiba 260-8675, Japan
| | - Ayahito Kaneki
- Graduate School of Pharmaceutical Sciences, Chiba University, 1-8-1, Inohana, Chuo-ku, Chiba 260-8675, Japan
| | - Yuma Mihara
- Graduate School of Pharmaceutical Sciences, Chiba University, 1-8-1, Inohana, Chuo-ku, Chiba 260-8675, Japan
| | - Honoka Yuba
- Graduate School of Pharmaceutical Sciences, Chiba University, 1-8-1, Inohana, Chuo-ku, Chiba 260-8675, Japan
| | - Tetsuhiro Nemoto
- Graduate School of Pharmaceutical Sciences, Chiba University, 1-8-1, Inohana, Chuo-ku, Chiba 260-8675, Japan
| |
Collapse
|
46
|
Li Y, Zhang RM, Xu X. Theoretical Kinetics studies of isoprene peroxy radical chemistry: The fate of Z-δ-(4-OH, 1-OO)-ISOPOO radical. ECOTOXICOLOGY AND ENVIRONMENTAL SAFETY 2023; 266:115553. [PMID: 37839188 DOI: 10.1016/j.ecoenv.2023.115553] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/13/2023] [Revised: 09/19/2023] [Accepted: 10/05/2023] [Indexed: 10/17/2023]
Abstract
The OH radical recycling mechanism in isoprene oxidation is one of the most exciting topics in atmospheric chemistry, and the corresponding studies expand our understanding of oxidation mechanisms of volatile organic compounds in the troposphere and provide reliable evidence to improve and develop conventional atmospheric models. In this work, we performed a detailed theoretical kinetics study on the Z-δ-(4-OH, 1-OO)-ISOPOO radical chemistry, which is proposed as the heart of OH recycling in isoprene oxidation. With the full consideration of its accumulation and consumption channels, we studied and discussed the fate of Z-δ-(4-OH, 1-OO)-ISOPOO radical by solving the energy-resolved master equation over a broad range of conditions, including not only room temperatures but also high temperatures of a forest fire or low temperatures and pressures of the upper troposphere. We found non-negligible pressure dependence of its fate at combustion temperatures (up to two orders of magnitude) and demonstrated the significance of both the multi-structural torsional anharmonicity and tunneling for accurately calculating kinetics of the studied system. More interestingly, the tunneling effect on the phenomenological rate constants of the H-shift reaction channel is also found to be pressure-dependent due to the competition with the O2 loss reaction. In addition, our time evolution calculations revealed a two-stage behavior of critical species in this reaction system and estimated the shortest half-lives for the Z-δ-(4-OH, 1-OO)-ISOPOO radical at various temperatures, pressures and altitudes. This detailed kinetics study of Z-δ-(4-OH, 1-OO)-ISOPOO radical chemistry offers a typical example to deeply understand the core mechanism of OH recycling pathways in isoprene oxidation, and provides valuable insights for promoting the development of relevant atmospheric models.
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
| | - Rui Ming Zhang
- 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
| | - Xuefei Xu
- 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
|
47
|
Erba A, Desmarais JK, Casassa S, Civalleri B, Donà L, Bush IJ, Searle B, Maschio L, Edith-Daga L, Cossard A, Ribaldone C, Ascrizzi E, Marana NL, Flament JP, Kirtman B. CRYSTAL23: A Program for Computational Solid State Physics and Chemistry. J Chem Theory Comput 2023; 19:6891-6932. [PMID: 36502394 PMCID: PMC10601489 DOI: 10.1021/acs.jctc.2c00958] [Citation(s) in RCA: 20] [Impact Index Per Article: 20.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2022] [Indexed: 12/14/2022]
Abstract
The Crystal program for quantum-mechanical simulations of materials has been bridging the realm of molecular quantum chemistry to the realm of solid state physics for many years, since its first public version released back in 1988. This peculiarity stems from the use of atom-centered basis functions within a linear combination of atomic orbitals (LCAO) approach and from the corresponding efficiency in the evaluation of the exact Fock exchange series. In particular, this has led to the implementation of a rich variety of hybrid density functional approximations since 1998. Nowadays, it is acknowledged by a broad community of solid state chemists and physicists that the inclusion of a fraction of Fock exchange in the exchange-correlation potential of the density functional theory is key to a better description of many properties of materials (electronic, magnetic, mechanical, spintronic, lattice-dynamical, etc.). Here, the main developments made to the program in the last five years (i.e., since the previous release, Crystal17) are presented and some of their most noteworthy applications reviewed.
Collapse
Affiliation(s)
- Alessandro Erba
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Jacques K. Desmarais
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Silvia Casassa
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Bartolomeo Civalleri
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Lorenzo Donà
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Ian J. Bush
- STFC
Rutherford Appleton Laboratory, Chilton Didcot, Oxfordshire OX11 0QX, United Kingdom
| | - Barry Searle
- SFTC
Daresbury Laboratory, Daresbury, Cheshire WA4 4AD, United Kingdom
| | - Lorenzo Maschio
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Loredana Edith-Daga
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Alessandro Cossard
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Chiara Ribaldone
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Eleonora Ascrizzi
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Naiara L. Marana
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Jean-Pierre Flament
- Université
de Lille, CNRS, UMR 8523 — PhLAM — Physique des Lasers, Atomes et Molécules, 59000 Lille, France
| | - Bernard Kirtman
- Department
of Chemistry and Biochemistry, University
of California, Santa
Barbara, California 93106, United States
| |
Collapse
|
48
|
Kandpal SC, Otukile KP, Jindal S, Senthil S, Matthews C, Chakraborty S, Moskaleva LV, Ramakrishnan R. Stereo-electronic factors influencing the stability of hydroperoxyalkyl radicals: transferability of chemical trends across hydrocarbons and ab initio methods. Phys Chem Chem Phys 2023; 25:27302-27320. [PMID: 37791466 DOI: 10.1039/d3cp03598k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/05/2023]
Abstract
The hydroperoxyalkyl radicals (˙QOOH) are known to play a significant role in combustion and tropospheric processes, yet their direct spectroscopic detection remains challenging. In this study, we investigate molecular stereo-electronic effects influencing the kinetic and thermodynamic stability of a ˙QOOH along its formation path from the precursor, alkylperoxyl radical (ROO˙), and the depletion path resulting in the formation of cyclic ether + ˙OH. We focus on reactive intermediates encountered in the oxidation of acyclic hydrocarbon radicals: ethyl, isopropyl, isobutyl, tert-butyl, neopentyl, and their alicyclic counterparts: cyclohexyl, cyclohexenyl, and cyclohexadienyl. We report reaction energies and barriers calculated with the highly accurate method Weizmann-1 (W1) for the channels: ROO˙ ⇌ ˙QOOH, ROO˙ ⇌ alkene + ˙OOH, ˙QOOH ⇌ alkene + ˙OOH, and ˙QOOH ⇌ cyclic ether + ˙OH. Using W1 results as a reference, we have systematically benchmarked the accuracy of popular density functional theory (DFT), composite thermochemistry methods, and an explicitly correlated coupled-cluster method. We ascertain inductive, resonance, and steric effects on the overall stability of ˙QOOH and computationally investigate the possibility of forming more stable species. With new reactions as test cases, we probe the capacity of various ab initio methods to yield quantitative insights on the elementary steps of combustion.
Collapse
Affiliation(s)
| | - Kgalaletso P Otukile
- Department of Chemistry, University of the Free State, PO Box 339, Bloemfontein 9300, South Africa.
| | - Shweta Jindal
- Tata Institute of Fundamental Research, Hyderabad 500046, India.
| | - Salini Senthil
- Tata Institute of Fundamental Research, Hyderabad 500046, India.
| | - Cameron Matthews
- Department of Chemistry, University of the Free State, PO Box 339, Bloemfontein 9300, South Africa.
| | | | - Lyudmila V Moskaleva
- Department of Chemistry, University of the Free State, PO Box 339, Bloemfontein 9300, South Africa.
| | | |
Collapse
|
49
|
Zhang D, Truhlar DG. An Accurate Density Coherence Functional for Hybrid Multiconfiguration Density Coherence Functional Theory. J Chem Theory Comput 2023; 19:6551-6556. [PMID: 37708640 DOI: 10.1021/acs.jctc.3c00741] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/16/2023]
Abstract
We present hybrid multiconfiguration density coherence functional theory (HMC-DCFT), and we optimize a density coherence functional by parametrization against a diverse data set of 59 bond energies and 60 barrier heights. We compare the results to calculations on the same data set by CASSCF, CASPT2, six Kohn-Sham and hybrid Kohn-Sham exchange-correlation functionals, and three on-top functionals for pair-density functional theory (PDFT) and hybrid PDFT. The new functional has better accuracy than all compared methods.
Collapse
Affiliation(s)
- Dayou Zhang
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
50
|
Richter R, Aschebrock T, Schelter I, Kümmel S. Meta-generalized gradient approximations in time dependent generalized Kohn-Sham theory: Importance of the current density correction. J Chem Phys 2023; 159:124117. [PMID: 38127400 DOI: 10.1063/5.0167972] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2023] [Accepted: 09/07/2023] [Indexed: 12/23/2023] Open
Abstract
We revisit the use of Meta-Generalized Gradient Approximations (mGGAs) in time-dependent density functional theory, reviewing conceptual questions and solving the generalized Kohn-Sham equations by real-time propagation. After discussing the technical aspects of using mGGAs in combination with pseudopotentials and comparing real-space and basis set results, we focus on investigating the importance of the current-density based gauge invariance correction. For the two modern mGGAs that we investigate in this work, TASK and r2SCAN, we observe that for some systems, the current density correction leads to negligible changes, but for others, it changes excitation energies by up to 40% and more than 0.8 eV. In the cases that we study, the agreement with the reference data is improved by the current density correction.
Collapse
Affiliation(s)
- Rian Richter
- Theoretical Physics IV, University of Bayreuth, 95440 Bayreuth, Germany
| | - Thilo Aschebrock
- Theoretical Physics IV, University of Bayreuth, 95440 Bayreuth, Germany
| | - Ingo Schelter
- Theoretical Physics IV, University of Bayreuth, 95440 Bayreuth, Germany
| | - Stephan Kümmel
- Theoretical Physics IV, University of Bayreuth, 95440 Bayreuth, Germany
| |
Collapse
|