MN15-L: A New Local Exchange-Correlation Functional for Kohn-Sham Density Functional Theory with Broad Accuracy for Atoms, Molecules, and Solids

Ruthenium complexes bearing nile red chromophore and one of its derivative: Theoretical evaluation of PDT-related properties

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.

Harmonic and anharmonic vibrational computations for biomolecular building blocks: Benchmarking DFT and basis sets by theoretical and experimental IR spectrum of glycine conformers

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.

Lone-Pair-π Bond Strength Unveiled by a Combined Experimental and Computational Study

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.

Machine learning for accuracy in density functional approximations

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.

Small Representative Databases for Testing and Validating Density Functionals and Other Electronic Structure Methods

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.

Benchmarking Functionals for Strong-Field Light-Matter Interactions in Adiabatic Time-Dependent Density Functional Theory

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.

Oxidative Addition of (Hetero)aryl (Pseudo)halides at Palladium(0): Origin and Significance of Divergent Mechanisms

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.

Efficiency Conceptualization Model: A Theoretical Method for Predicting the Turnover of Catalysts

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.

Electronic structure analysis and DFT benchmarking of Rydberg-type alkali-metal-crown ether, -cryptand, and -adamanzane complexes

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.

Global Potential Energy Surfaces by Compressed-State Multistate Pair-Density Functional Theory for Hyperthermal Collisions in the O2+O2 System

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.

Origins of the unphysical noncovalent interaction energy curves obtained with the 2011 and 2012 Minnesota density functionals

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.

Comprehensive Multipath Variational Kinetics Study on Hydrogen Abstraction Reactions from Three Typical Dimethylcyclohexane Isomers by Hydroxyl Radicals: from the Electronic Structure to Model Applications

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.

Quantitative kinetics reveal that reactions of HO2 are a significant sink for aldehydes in the atmosphere and may initiate the formation of highly oxygenated molecules via autoxidation

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.

Development of a multi-step screening procedure for redox active molecules in organic radical polymer anodes and as redox flow anolytes

Benzo[d]-X-zolyl-pyridinyl (XO, 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.

Why does the orientation of azulene affect the two-photon activity of a porphyrinoid-azulene system?

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.

DFT Struggles to Predict the Energy Landscape for Iron Pyridine Diimine-Catalyzed [2 + 2] Cycloaddition of Alkenes: Insights into the Problem and Alternative Solutions

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.

A Theoretical Benchmark of the Geometric and Optical Properties for 3d Transition Metal Nanoclusters via Density Functional Theory

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.

Pt-S Bond-Enabled Temperature-Dependent Phosphorescence in S-Heteroaryl Tetradentate Pt(S^C^N^O) Complexes

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.

The p-block challenge: assessing quantum chemistry methods for inorganic heterocycle dimerizations

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.

Mechanistic Origin of Ligand Effects on Exhaustive Functionalization During Pd-Catalyzed Cross-Coupling of Dihaloarenes

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.

Reversible Coupling of Germylone with Isocyanates

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.

Spurious Oscillations Caused by Density Functional Approximations: Who is to Blame? Exchange or Correlation?

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.

IR spectroscopic characterization of 3d transition metal carbene cations, FeCH2+ and CoCH2+: periodic trends and a challenge for DFT approaches

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.

Structure and dynamics of liquid water from ab initio simulations: adding Minnesota density functionals to Jacob's ladder

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.

DELFI: a computer oracle for recommending density functionals for excited states calculations

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.

Conformational energies of biomolecules in solution: Extending the MPCONF196 benchmark with explicit water molecules

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.

Modeling Multi-Step Organic Reactions: Can Density Functional Theory Deliver Misleading Chemistry?

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.

Spectroscopic characterization of carbon monoxide activation by neutral chromium carbides

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.

Rapid Atmospheric Reactions between Criegee Intermediates and Hypochlorous Acid

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.

Regulating the Flexibility to Assemble Porous Single-Atom Fe-Coordinated Metallopolymers for Efficient Heterogeneous Catalytic Oxidations

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.

Paramagnetic Nuclear Magnetic Resonance Shifts for Triplet Systems and Beyond with Modern Relativistic Density Functional Methods

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.

Exploring CPS-Extrapolated DLPNO-CCSD(T1) Reference Values for Benchmarking DFT Methods on Enzymatically Catalyzed Reactions

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.

Revisiting Gauge-Independent Kinetic Energy Densities in Meta-GGAs and Local Hybrid Calculations of Magnetizabilities

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.

M11pz: A Nonlocal Meta Functional with Zero Hartree-Fock Exchange and with Broad Accuracy for Chemical Energies and Structures

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.

DAPD Set of Pd-Containing Diatomic Molecules: Accurate Molecular Properties and the Great Lengths to Obtain Them

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.

Development of molecular cluster models to probe pyrite surface reactivity

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.

A step toward density benchmarking-The energy-relevant "mean field error"

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.

Benchmarking Density Functional Theory Methods for Metalloenzyme Reactions: The Introduction of the MME55 Set

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.

Comparing Density Functional Theory Metal-Ligand Bond Dissociation Enthalpies with Experimental Solution-Phase Enthalpies of Activation for Bond Dissociation

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.

Synthesis and Characterization of U≡C Triple Bonds in Fullerene Compounds

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.

Unraveling actinide-actinide bonding in fullerene cages: a DFT versus ab initio methodological study

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.

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

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.

A Guanidine-Supported π-Complex of Germanium Amenable to Intramolecular C-C Cleavage in Arene and Ge Atom Transfer

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.

Short-range screened density matrix functional for proper descriptions of thermochemistry, thermochemical kinetics, nonbonded interactions, and singlet diradicals

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.

Scandium-Catalyzed Phenol-Directed Construction of 5-Carbonyl-4-hydroxybenzofurans via Intramolecular Friedel-Crafts Reaction

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.

Theoretical Kinetics studies of isoprene peroxy radical chemistry: The fate of Z-δ-(4-OH, 1-OO)-ISOPOO radical

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.

CRYSTAL23: A Program for Computational Solid State Physics and Chemistry

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.

Stereo-electronic factors influencing the stability of hydroperoxyalkyl radicals: transferability of chemical trends across hydrocarbons and ab initio methods

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.

An Accurate Density Coherence Functional for Hybrid Multiconfiguration Density Coherence Functional Theory

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.

Meta-generalized gradient approximations in time dependent generalized Kohn-Sham theory: Importance of the current density correction

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.