Surface Adsorption Energetics Studied with "Gold Standard" Wave-Function-Based Ab Initio Methods: Small-Molecule Binding to TiO2(110)

State-of-the-art local correlation methods enable affordable gold standard quantum chemistry for up to hundreds of atoms

In this feature, we review the current capabilities of local electron correlation methods up to the coupled cluster model with single, double, and perturbative triple excitations [CCSD(T)], which is a gold standard in quantum chemistry. The main computational aspects of the local method types are assessed from the perspective of applications, but the focus is kept on how to achieve chemical accuracy (i.e., <1 kcal mol-1 uncertainty), as well as on the broad scope of chemical problems made accessible. The performance of state-of-the-art methods is also compared, including the most employed DLPNO and, in particular, our local natural orbital (LNO) CCSD(T) approach. The high accuracy and efficiency of the LNO method makes chemically accurate CCSD(T) computations accessible for molecules of hundreds of atoms with resources affordable to a broad computational community (days on a single CPU and 10-100 GB of memory). Recent developments in LNO-CCSD(T) enable systematic convergence and robust error estimates even for systems of complicated electronic structure or larger size (up to 1000 atoms). The predictive power of current local CCSD(T) methods, usually at about 1-2 order of magnitude higher cost than hybrid density functional theory (DFT), has become outstanding on the palette of computational chemistry applicable for molecules of practical interest. We also review more than 50 LNO-based and other advanced local-CCSD(T) applications for realistic, large systems across molecular interactions as well as main group, transition metal, bio-, and surface chemistry. The examples show that properly executed local-CCSD(T) can contribute to binding, reaction equilibrium, rate constants, etc. which are able to match measurements within the error estimates. These applications demonstrate that modern, open-access, and broadly affordable local methods, such as LNO-CCSD(T), already enable predictive computations and atomistic insight for complicated, real-life molecular processes in realistic environments.

Investigating the Basis Set Convergence of Diagrammatically Decomposed Coupled-Cluster Correlation Energy Contributions for the Uniform Electron Gas

We investigate the convergence of coupled-cluster (CC) correlation energies and related quantities with respect to the employed basis set size for the uniform electron gas (UEG) to gain a better understanding of the basis set incompleteness error (BSIE). To this end, coupled-cluster doubles (CCD) theory is applied to the three-dimensional UEG for a range of densities, basis set sizes, and electron numbers. We present a detailed analysis of individual diagrammatically decomposed contributions to the amplitudes at the level of CCD theory. In particular, we show that only two terms from the amplitude equations contribute to the asymptotic large-momentum behavior of the transition structure factor, corresponding to the cusp region at short interelectronic distances. However, due to the coupling present in the amplitude equations, all decomposed correlation energy contributions show the same asymptotic convergence behavior to the complete basis set limit. These findings provide an additional rationale for the success of a recently proposed correction to the BSIE of CC theory. Lastly, we examine the BSIE in the CCD plus perturbative triples [CCD(T)] method, as well as in the newly proposed CCD plus complete perturbative triples [CCD(cT)] method.

Going for Gold(-Standard): Attaining Coupled Cluster Accuracy in Oxide-Supported Nanoclusters

The structure of oxide-supported metal nanoclusters plays an essential role in their sharply enhanced catalytic activity over that of bulk metals. Simulations provide the atomic-scale resolution needed to understand these systems. However, the sensitive mix of metal-metal and metal-support interactions, which govern their structure, puts stringent requirements on the method used, requiring calculations beyond standard density functional theory (DFT). The method of choice is coupled cluster theory [specifically CCSD(T)], but its computational cost has so far prevented its application to these systems. In this work, we showcase two approaches to make CCSD(T) accuracy readily achievable in oxide-supported nanoclusters. First, we leverage the SKZCAM protocol to provide the first benchmarks of oxide-supported nanoclusters, revealing that it is specifically metal-metal interactions that are challenging to capture with DFT. Second, we propose a CCSD(T) correction (ΔCC) to the metal-metal interaction errors in DFT, reaching accuracy comparable to that of the SKZCAM protocol at significantly lower cost. This approach forges a path toward studying larger systems at reliable accuracy, which we highlight by identifying a ground-state structure in agreement with experiments for Au20 on MgO, a challenging system where DFT models have yielded conflicting predictions.

A theoretical spectroscopy study of the photoluminescence properties of narrow band Eu2+-doped phosphors containing multiple candidate doping centers. Prediction of an unprecedented narrow band red phosphor

We have previously presented a computational protocol that is based on an embedded cluster model and operates in the framework of TD-DFT in conjunction with the excited state dynamics (ESD) approach. The protocol is able to predict the experimental absorption and emission spectral shapes of Eu2+-doped phosphors. In this work, the applicability domain of the above protocol is expanded to Eu2+-doped phosphors bearing multiple candidate Eu doping centers. It will be demonstrated that this protocol provides full control of the parameter space that describes the emission process. The stability of Eu doping at various centers is explored through local energy decomposition (LED) analysis of DLPNO-CCSD(T) energies. This enables further development of the understanding of the electronic structure of the targeted phosphors, the diverse interactions between Eu and the local environment, and their impact on Eu doping probability, and control of the emission properties. Hence, it can be employed to systematically improve deficiencies of existing phosphor materials, defined by the presence of various intensity emission bands at undesired frequencies, towards classes of candidate Eu2+-doped phosphors with desired narrow band red emission. For this purpose, the chosen study set consists of three UCr4C4-based narrow-band phosphors, namely the known alkali lithosilicates RbNa[Li3SiO4]2:Eu2+ (RNLSO2), RbNa3[Li3SiO4]4:Eu2+ (RNLSO) and their isotypic nitridolithoaluminate phosphors consisting of CaBa[LiAl3N4]2:Eu2+ (CBLA2) and the proposed Ca3Ba[LiAl3N4]4:Eu2+ (CBLA), respectively. The theoretical analysis presented in this work led us to propose a modification of the CBLA2 phosphor that should have improved and unprecedented narrow band red emission properties. Finally, we believe that the analysis presented here is important for the future rational design of novel Eu2+-doped phosphor materials, with a wide range of applications in science and technology.

Level-Shifted Embedded Cluster Method for Modeling the Chemistry of Metal Oxides

The embedded cluster method has been used extensively in the study of the chemical and physical properties of metal oxides. This method has been a popular tool due to its relatively high accuracy and low computational cost. An even more promising option may entail integrating the embedded cluster method with the combined quantum mechanical and molecular mechanical (QM/MM) approach, thereby enabling further consideration of interactions within the entire system for superior results. We aim to accurately model the chemistry of metal oxides using this combined scheme. Here, using the prototypical MgO(100) surface as a test system, with Mg9O14 as the cluster in the quantum mechanical region, we show that the embedded cluster with untailored boundary effective core potentials (ECPs) can have frontier orbital energy levels that substantially deviate from the quantum mechanical reference results. This occurs even when Mg9O9, which retains the stoichiometry of MgO, is used as the cluster in the quantum mechanical region. As a result, the chemical properties of the embedded cluster models differ from those of the quantum mechanical reference model. To address this issue, we propose a new variant of the embedded cluster method called the level-shifted embedded cluster (LSEC) method, which allows the energy levels to be shifted to match the reference levels by tuning the boundary ECPs. Our validation calculations on the adsorption of various adsorbates with different properties on the MgO(100) surface show that the overall performance of QM/MM with the LSEC method is excellent for the adsorption energies, geometries, and charge properties. The excellent performance holds for both the nonstoichiometric and stoichiometric clusters (i.e., Mg9O14 and Mg9O9, respectively), demonstrating the robustness of the LSEC method. We expect that the LSEC method can be combined with QM/MM or used separately for future chemical studies of metal oxides and other ionically bonded systems.

Many-Body Methods for Surface Chemistry Come of Age: Achieving Consensus with Experiments

The adsorption energy of a molecule onto the surface of a material underpins a wide array of applications, spanning heterogeneous catalysis, gas storage, and many more. It is the key quantity where experimental measurements and theoretical calculations meet, with agreement being necessary for reliable predictions of chemical reaction rates and mechanisms. The prototypical molecule-surface system is CO adsorbed on MgO, but despite intense scrutiny from theory and experiment, there is still no consensus on its adsorption energy. In particular, the large cost of accurate many-body methods makes reaching converged theoretical estimates difficult, generating a wide range of values. In this work, we address this challenge, leveraging the latest advances in diffusion Monte Carlo (DMC) and coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] to obtain accurate predictions for CO on MgO. These reliable theoretical estimates allow us to evaluate the inconsistencies in published temperature-programed desorption experiments, revealing that they arise from variations in employed pre-exponential factors. Utilizing this insight, we derive new experimental estimates of the (electronic) adsorption energy with a (more) precise pre-exponential factor. As a culmination of all of this effort, we are able to reach a consensus between multiple theoretical calculations and multiple experiments for the first time. In addition, we show that our recently developed cluster-based CCSD(T) approach provides a low-cost route toward achieving accurate adsorption energies. This sets the stage for affordable and reliable theoretical predictions of chemical reactions on surfaces to guide the realization of new catalysts and gas storage materials.

High-level ab initio evidence of bipyramidal Cu5 clusters as fluxional Jahn-Teller molecules

Novel highly selective synthesis techniques have enable the production of atomically precise monodisperse metal clusters (AMCs) of subnanometer size. These AMCs exhibit 'molecule-like' structures that have distinct physical and chemical properties, significantly different from those of nanoparticles and bulk material. In this work, we study copper pentamer Cu5 clusters as model AMCs by applying both density functional theory (DFT) and high-level (wave-function-based) ab initio methods, including those which are capable of accounting for the multi-state multi-reference character of the wavefunction at the conical intersection (CI) between different electronic states and augmenting the electronic basis set till achieving well-converged energy values and structures. After assessing the accuracy of a high-level multi-multireference ab initio protocol for the well-known Cu3 case, we apply it to demonstrate that bypiramidal Cu5 clusters are distorted Jahn-Teller (JT) molecules. The method is further used to evaluate the accuracy of single-reference approaches, finding that the coupled cluster singles and doubles and perturbative triples CCSD(T) method delivers the results closer to our ab initio predictions and that dispersion-corrected DFT can outperform the CCSD method. Finally, we discuss how JT effects and, more generally, conical intersections, are intimately connected to the fluxionality of AMCs, giving them a 'floppy' character that ultimately facilitates their interaction with environmental molecules and thus enhances their functioning as catalysts.

New Insights on Singlet Oxygen Release from Li-Air Battery Cathode: Periodic DFT Versus CASPT2 Embedded Cluster Calculations

Li-air batteries are a promising energy storage technology for large-scale applications, but the release of highly reactive singlet oxygen (1O2) during battery operation represents a main concern that sensibly limits their effective deployment. An in-depth understanding of the reaction mechanisms underlying the 1O2 formation is crucial to prevent its detrimental reactions with the electrolyte species. However, describing the elusive chemistry of highly correlated species such as singlet oxygen represents a challenging task for state-of-the-art theoretical tools based on density functional theory. Thus, in this study, we apply an embedded cluster approach, based on CASPT2 and effective point charges, to address the evolution of 1O2 at the Li2O2 surface during oxidation, i.e., the battery charging process. Based on recent hypothesis, we depict a feasible O22-/O2-/O2 mechanisms occurring from the (112̅0)-Li2O2 surface termination. Our highly accurate calculations allow for the identification of a stable superoxide as local minimum along the potential energy surface (PES) for 1O2 release, which is not detected by periodic DFT. We find that 1O2 release proceeds via a superoxide intermediate in a two-step one-electron process or another still accessible pathway featuring a one-step two-electron mechanism. In both cases, it represents a feasible product of Li2O2 oxidation upon battery charging. Thus, tuning the relative stability of the intermediate superoxide species can enable key strategies aiming at controlling the detrimental development of 1O2 for new and highly performing Li-air batteries.

Polyhedral Dicobaltadithiaboranes and Dicobaltdiselenaboranes as Examples of Bimetallic Nido Structures without Bridging Hydrogens

The geometries and energetics of the n-vertex polyhedral dicobaltadithiaboranes and dicobaltadiselenaboranes Cp2Co2E2Bn−4Hn−4 (E = S, Se; n = 8 to 12) have been investigated via the density functional theory. Most of the lowest-energy structures in these systems are generated from the (n + 1)-vertex most spherical closo deltahedra by removal of a single vertex, leading to a tetragonal, pentagonal, or hexagonal face depending on the degree of the vertex removed. In all of these low-energy structures, the chalcogen atoms are located at the vertices of the non-triangular face. Alternatively, the central polyhedron in most of the 12-vertex structures can be derived from a Co2E2B8 icosahedron with adjacent chalcogen (E) vertices by breaking the E–E edge and 1 or more E–B edges to create a hexagonal face. Examples of the arachno polyhedra with two tetragonal and/or pentagonal faces derived from the removal of two vertices from isocloso deltahedra were found among the set of lowest-energy Cp2Co2E2Bn−4Hn−4 (E = S, Se; n = 8 and 12) structures.

Aromaticity in P8 allotropes and (CH)8 analogues: significance of their 40 valence electrons?

The currently unknown phosphorus allotrope P₈ is of interest since its 40 total valence electrons is a "magic number" corresponding to a filled 1S²1P⁶1D¹⁰1S²1F¹⁴2P⁶ shell such as found in the relatively stable main group element clusters Al₁₃^- and Ge₉^4-. However, P₈ still remains as an elusive structure not realized experimentally. The lowest energy P₈ structure by a margin of ∼9 kcal mol^-1 is shown by density functional theory to be a cuneane analogue with no PP double bonds and two each of P₅, P₄, and P₃ rings. Higher energy P₈ structures are polycyclic systems having at most a single PP double bond. These P₈ systems are not "carbon copies" of the corresponding (CH)₈ hydrocarbons with exactly one hydrogen atom bonded to each carbon atom. Thus the lowest energy (CH)₈ structure is cyclooctatetraene with four CC bonds followed by benzocyclobutene with three CC bonds. The cuneane (CH)₈ structure is a relatively high energy isomer lying ∼36 kcal mol^-1 above cyclooctatetraene. The cubane P₈ and (CH)₈ structures are even higher energy structures, lying ∼37 and ∼74 kcal mol^-1 in energy above the corresponding global minima. Our results demonstrate differences in medium sized aggregates of elemental phosphorus and isolobal hydrocarbon species.

Electronic and Optical Properties of Eu2+-Activated Narrow-Band Phosphors for Phosphor-Converted Light-Emitting Diode Applications: Insights from a Theoretical Spectroscopy Perspective

In this work, we present a computational protocol that is able to predict the experimental absorption and emission spectral shapes of Eu²⁺-doped phosphors. The protocol is based on time-dependent density functional theory and operates in conjunction with an excited-state dynamics approach. It is demonstrated that across the study set consisting of representative examples of nitride, oxo-nitride, and oxide Eu²⁺-doped phosphors, the energy distribution and the band shape of the emission spectrum are related to the nature of the 4f–5d transitions that are probed in the absorption process. Since the 4f orbitals are very nearly nonbonding, the decisive quantity is the covalency of the 5d acceptor orbitals that become populated in the electronically excited state that leads to emission. The stronger the (anti) bonding interaction between the lanthanide and the ligands is in the excited state, the larger will be the excited state distortion. Consequently, the corresponding emission will get broader due to the vibronic progression that is induced by the structural distortion. In addition, the energy separation of the absorption bands that are dominated by states with valence 4f–5d and a metal to ligand charge transfer character defines a measure for the thermal quenching of the studied Eu²⁺-doped phosphors. Based on this analysis, simple descriptors are identified that show a strong correlation with the energy position and bandwidth of the experimental emission bands without the need for elaborate calculations. Overall, we believe that this study serves as an important reference for designing new Eu²⁺-doped phosphors with desired photoluminescence properties.

First-principles modelling of the new generation of subnanometric metal clusters: Recent case studies

The very recent development of highly selective techniques making possible the synthesis and experimental characterization of subnanometric (subnanometer-sized) metal clusters (even single atoms) is pushing our understanding far beyond the present knowledge in materials science, driving these clusters as a new generation of quantum materials at the lower bounds of nanotechnology. When the size of the metal cluster is reduced to a small number of atoms, the d-band of the metal splits into a subnanometric d-type molecular orbitals network in which all metal atoms are inter-connected, with the inter-connections having the length of a chemical bond (1-2 Å). These molecular characteristics are at the very core of the high stability and novel properties of the smallest metal clusters, with their integration into colloidal materials interacting with the environment having the potential to further boost their performance in applications such as luminescence, sensing, bioimaging, theranostics, energy conversion, catalysis, and photocatalysis. Through the presentation of very recent case studies, this Feature Article is aimed to illustrate how first-principles modelling, including methods beyond the state-of-the-art and an interplay with cutting-edge experiments, is helping to understand the special properties of these clusters at the most fundamental level. Moreover, it will be discussed how superfluid helium droplets can act both as nano-reactors and carriers to achieve the synthesis and surface deposition of metal clusters. This concept will be illustrated with the quantum simulation of the helium droplet-assisted soft-landing of a single Au atom onto a titanium dioxide (TiO₂) surface. Next, it will be shown how the application of first-principles methods have disclosed the fundamental reasons why subnanometric Cu₅ clusters are resistant to irreversible oxidation, and capable of increasing and extending into the visible region the solar absorption of TiO₂, of augmenting its efficiency for photo-catalysis beyond a factor of four, also considering the decomposition and photo-activation of CO₂ as a prototypical (photo-) catalytic reaction. Finally, I will discuss how the modification of the same material with subnanometric Ag₅ clusters has converted it into a "reporter" of a surface polaron property as well as a novel two-dimensional polaronic material.

General embedded cluster protocol for accurate modeling of oxygen vacancies in metal-oxides

The O vacancy (Ov) formation energy, EOv, is an important property of a metal-oxide, governing its performance in applications such as fuel cells or heterogeneous catalysis. These defects are routinely studied with density functional theory (DFT). However, it is well-recognized that standard DFT formulations (e.g., the generalized gradient approximation) are insufficient for modeling the Ov, requiring higher levels of theory. The embedded cluster method offers a promising approach to compute EOv accurately, giving access to all electronic structure methods. Central to this approach is the construction of quantum(-mechanically treated) clusters placed within suitable embedding environments. Unfortunately, current approaches to constructing the quantum clusters either require large system sizes, preventing application of high-level methods, or require significant manual input, preventing investigations of multiple systems simultaneously. In this work, we present a systematic and general quantum cluster design protocol that can determine small converged quantum clusters for studying the Ov in metal-oxides with accurate methods, such as local coupled cluster with single, double, and perturbative triple excitations. We apply this protocol to study the Ov in the bulk and surface planes of rutile TiO2 and rock salt MgO, producing the first accurate and well-converged determinations of EOv with this method. These reference values are used to benchmark exchange–correlation functionals in DFT, and we find that all the studied functionals underestimate EOv, with the average error decreasing along the rungs of Jacob’s ladder. This protocol is automatable for high-throughput calculations and can be generalized to study other point defects or adsorbates.

Theoretical investigation of CH-bond activation by photocatalytic excited SO2 and the effects of C-, N-, S-, and Se-doped TiO2

The photocatalytic sulfoxidation on TiO₂ discovered by Parrino et al. represents a new, interesting and lower energy route for the synthesis of sulfonic acids. Sulfonic acids are important precursors for dyes, detergents and drugs. In the commonly known industrial process, SO₂ and a specific hydrocarbon are converted into sulfonic acids using high-energy UV light. In this reaction, SO₂ is excited into a metastable triplet state (³SO₂), which has the potential to activate a CH-bond of hydrocarbons and start a radical reaction cycle. By introducing TiO₂ as a photocatalyst, it has been shown that visible light can be used for the synthesis. This offers the potential to be a cost-effective reaction approach for industrial use. However, experimental studies indicate that the initial excitation mechanism of SO₂ on TiO₂ is significantly different from the catalyst-free mechanism. Parrino et al. were able to reveal first evidence for the existence of a charge-transfer process from SO₂ to the TiO₂ surface by means of electrochemical experiments. First theoretical investigations from first principles were able to further substantiate the existence of a charge-transfer. However, to fully understand this mechanism, more accurate methods such as Time Dependent Density Functional Theory (TD-DFT) or ab initio multireference methods such as the Complete Active Space Self Consistent Field (CASSCF) method are required. Furthermore, after understanding the charge-transfer mechanism, the introduction of dopants into TiO₂ can be investigated in order to possibly redshift the excitation energy. This might open the route to using lower energy light for the sulfoxidation of hydrocarbons on TiO₂ as a new potential industrial reaction for the synthesis of sulfonic acids. In this work, we will study the initial step of the photocatalytic sulfoxidation of hydrocarbons using the TD-DFT and CASSCF methods by using a combined approach consisting of calculations with periodic boundary conditions and a newly constructed embedded cluster model. Furthermore, we will explore the effects of doping by introducing four heteroatoms (C, N, S, and Se) into the TiO₂ surfaces anatase[101] and rutile[110] to find a possible enhancement of the photocatalytic reactivity by lowering the electronic excitation energy.

Mechanism for the Reaction of White Phosphorus with Cp2Cr2(CO)6 Leading Ultimately to the Triple-Decker Sandwich Cp2Cr2(μ-η5,η5-P5): A Theoretical Study

The experimentally known reaction of Cp₂Cr₂(CO)₆ with white phosphorus (P₄) to give CpCr(CO)₂(η³-P₃), Cp₂Cr₂(CO)₄(μ-η,²η²-P₂), and the triple-decker sandwich Cp₂Cr₂(μ-η,⁵η⁵-P₅) is of interest since the P₄ reactant having a tetrahedral cluster of four phosphorus atoms is converted to products having P₂, P₃, and P₅ ligands. The mechanism of this obviously complicated reaction can be dissected into three stages using a coupled cluster theoretical method that has been benchmarked with the P₂, Mn(CO)₅, and CpCr(CO)₃ dimerization processes. The first stage of the Cp₂Cr₂(CO)₆/P₄ reaction mechanism generates the unsaturated singlet intermediate Cp₂Cr₂(CO)₅ that combines with the P₄ reactant. Decarbonylation of the resulting Cp₂Cr₂(CO)₅(P₄) complex provides a singlet tetracarbonyl readily fragmenting into the stable triphosphacyclopropenyl complex CpCr(CO)₂(η³-P₃) and the chromium phosphide CpCr(CO)₂(P). The isomeric triplet tetracarbonyl Cp₂Cr₂(CO)₄(P₄), readily fragments into CpCr(CO)₂(η²-P₂), which can generate the stable diphosphaacetylene complex Cp₂Cr₂(CO)₄(η,²η²-P₂) as well as the pentamer [CpCr(CO)₂]₅(P₁₀). Combination of the coordinately unsaturated CpCr(CO)(η³-P₃) with CpCr(CO)₂(η²-P₂) can lead to a ring expansion. This generates the P₅ pentagonal ligand in a Cp₂Cr₂(CO)₃(P₅) precursor to the experimentally observed carbonyl-free triple-decker sandwich Cp₂Cr₂(μ-η,⁵η⁵-P₅) after three successive decarbonylations.

Selective Oxidation of 5-Hydroxymethylfurfural to 2,5-Diformylfuran by Visible Light-Driven Photocatalysis over In Situ Substrate-Sensitized Titania

Solar energy-driven processes for biomass valorization are priority for the growing industrialized society. To address this challenge, efficient visible light-active photocatalyst for the selective oxidation of biomass-derived platform chemical is highly desirable. Herein, selective oxidation of 5-hydroxymethylfurfural (HMF) to 2,5-diformylfuran (DFF) was achieved by visible light-driven photocatalysis over titania. Pristine titania is photocatalytically inactive under visible light, so an unconventional approach was employed for the visible light (λ=515 nm) sensitization of titania via a formation of a visible light-absorbing complex of HMF (substrate) on the titania surface. Surface-complexation of HMF on titania mediated ligand-to-metal charge transfer (LMCT) under visible light, which efficiently catalyzed the oxidation of HMF to DFF. A high DFF selectivity of 87 % was achieved with 59 % HMF conversion after 4 h of illumination. The apparent quantum yield obtained for DFF production was calculated to be 6.3 %. It was proposed that the dissociative interaction of hydroxyl groups of HMF and the titania surface is responsible for the surface-complex formation. When the hydroxyl groups of titania were modified via surface-fluorination or calcination the oxidation of HMF was inhibited under visible light, signifying that hydroxyl groups are decisive for photocatalytic activity.

Computation of NMR Shielding Constants for Solids Using an Embedded Cluster Approach with DFT, Double-Hybrid DFT, and MP2

In this work, we explore the accuracy of post-Hartree–Fock (HF) methods and double-hybrid density functional theory (DFT) for the computation of solid-state NMR chemical shifts. We apply an embedded cluster approach and investigate the convergence with cluster size and embedding for a series of inorganic solids with long-range electrostatic interactions. In a systematic study, we discuss the cluster design, the embedding procedure, and basis set convergence using gauge-including atomic orbital (GIAO) NMR calculations at the DFT and MP2 levels of theory. We demonstrate that the accuracy obtained for the prediction of NMR chemical shifts, which can be achieved for molecular systems, can be carried over to solid systems. An appropriate embedded cluster approach allows one to apply methods beyond standard DFT even for systems for which long-range electrostatic effects are important. We find that an embedded cluster should include at least one sphere of explicit neighbors around the nuclei of interest, given that a sufficiently large point charge and boundary effective potential embedding is applied. Using the pcSseg-3 basis set and GIAOs for the computation of nuclear shielding constants, accuracies of 1.6 ppm for ⁷Li, 1.5 ppm for ²³Na, and 5.1 ppm for ³⁹K as well as 9.3 ppm for ¹⁹F, 6.5 ppm for ³⁵Cl, 7.4 ppm for ⁷⁹Br, and 7.5 ppm for ²⁵Mg as well as 3.8 ppm for ⁶⁷Zn can be achieved with MP2. Comparing various DFT functionals with HF and MP2, we report the superior quality of results for methods that include post-HF correlation like MP2 and double-hybrid DFT.

Photodesorption mechanism of water on WO3(001) - a combined embedded cluster, computational intelligence and wave packet approach

In this work we investigate the mechanism of photodesorption of water from a WO3(001) surface by theoretical calculations, applying an embedded cluster model. Using the CASSCF method, we have calculated both the ground state as well as the energetically preferred charge-transfer state in three degrees of freedom of the water molecule on the surface. The calculated potential energy surfaces were afterwards fitted with a neural network optimized by a genetic algorithm. A final quantum dynamic wave packet study provided insight into the photodesorption mechanism.

Efficient and accurate description of adsorption in zeolites

Accurate theoretical methods are needed to correctly describe adsorption on solid surfaces or in porous materials. The random phase approximation (RPA) with singles corrections scheme and the second order Møller-Plesset perturbation theory (MP2) are two schemes, which offer high accuracy at affordable computational cost. However, there is little knowledge about their applicability and reliability for different adsorbates and surfaces. Here, we calculate adsorption energies of seven different molecules in zeolite chabazite to show that RPA with singles corrections is superior to MP2, not only in terms of accuracy but also in terms of computer time. Therefore, RPA with singles is a suitable scheme for obtaining highly accurate adsorption energies in porous materials and similar systems.

Ab Initio Calculations for Molecule-Surface Interactions with Chemical Accuracy

Atomistic understanding of complex surface phenomena such as heterogeneous catalysis or storage and separation of energy-relevant gases in nanoporous materials (zeolites; metal-organic frameworks, MOFs) requires knowledge about reaction energies and energy barriers for elementary steps. This is difficult to obtain from experiment since the number of possible chemical, adsorption/desorption, and diffusion steps coupled to complex reaction networks is large, and so is the number of possible surface sites. Here is an important role of quantum chemistry which can provide rate and equilibrium constants for elementary steps "ab initio." To be useful, the predictions have to reach chemical accuracy (4 kJ/mol) which is difficult to achieve because realistic models of the surface systems may comprise of the order of a thousand atoms. While density functional theory (DFT) as a rule cannot be trusted to yield results within chemical accuracy limits, methods that are accurate enough (Coupled Cluster with Single, Double, and perturbative Triple Substitution, CCSD(T)) cannot be applied because of their exponential scaling with system size. This Account presents a hybrid high-level-low-level quantum method that combines DFT with dispersion for the full periodic system with second order Møller-Plesset perturbation theory (MP2) for the reaction site within a mechanical embedding scheme. In addition, to check if MP2 is accurate enough, we calculate Coupled Cluster (CC) corrections with Single, Double, and perturbatively treated Triple substitutions (CCSD(T)) for sufficiently small models of the reaction site. This multilevel hybrid MP2:DFT-D+ΔCC method is shown to yield chemical accuracy for a set of 12 molecule-surface interaction systems for which reliable experimental data are available. For CO/MgO(001), the history of the experiment-theory comparison illustrates two problems: (i) Do experiment and theory look at the same surface site? (ii) Does theory calculate the same quantity as derived from experiment? The hybrid MP2:DFT-D+ΔCC data set generated includes the MgO(001) surface, the Mg₂(dobdc) metal-organic framework, and the proton forms of the CHA and MFI zeolites interacting with the H₂, N₂, CO, CO₂, CH₄, and C₂H₆ molecules. It serves two purposes. First, it will be useful for testing density functionals with respect to their performance for molecule-surface interactions. Second, it establishes the hybrid MP2:DFT-D+ΔCC method as a reliable and powerful tool for ab initio predictions of adsorption and reaction energies as well as energy barriers when testing reaction mechanisms. For adsorption of small molecules in MOFs, isotherm predictions have reached a level of accuracy that deviations between theoretical predictions and experiments indicate sample imperfections. For elementary steps of the industrially important methanol-to-olefin process, our hybrid MP2:PBE+D+ΔCC calculations yield rate constants in agreement with experiment within chemical accuracy limits, finally achieving for molecule-surface reactions which was possible so hitherto only for gas phase reactions involving not more than 10 atoms.

Duality of Ring and Ladder Diagrams and Its Importance for Many-Electron Perturbation Theories

We present a diagrammatic decomposition of the transition pair correlation function for the uniform electron gas. We demonstrate explicitly that ring and ladder diagrams are dual counterparts that capture significant long- and short-ranged interelectronic correlation effects, respectively. Our findings help to guide the further development of approximate many-electron theories and reveal that the contribution of the ladder diagrams to the electronic correlation energy can be approximated in an effective manner using second-order perturbation theory. We employ the latter approximation to reduce the computational cost of coupled cluster theory calculations for insulators and semiconductors by 2 orders of magnitude without compromising accuracy.

Water Adsorption on Ideal Anatase-TiO2(101) – An Embedded Cluster Model for Accurate Adsorption Energetics and Excited State Properties

A combined theoretical approach towards the accurate description of water on anatase-TiO2(101) was pursued in this study. Firstly, periodic slab calculations on the basis of density hybrid functionals (PBE0, HSE06) were performed in order to gain insight into the adsorption sites and geometric structure of the surface. For submonolayer coverage of H2O, the molecular adsorption of water is found to be the most stable one with quite similar energetics in PBE0 and HSE06. Moreover, the transition states towards the less preferred dissociative adsorption forms are predicted to be greater than 0.7 eV. Thus, water will not spontaneously dissociate and based on the Computational Hydrogen Electrode model an overpotential of about 1.71 V is needed to drive the overall oxidation. In addition, to validate our results for molecular adsorption of H2O, an embedded cluster model is carefully evaluated for the a-TiO2(101) surface based on the periodic slab calculations. Subsequent high-level DLPNO-CCSD(T) results are in close agreement with our periodic slab calculations since the interaction is found to mainly consist of electrostatic contributions which are captured by hybrid functionals. Finally, first results on optimized geometries in the excited state based on the photogenerated charge-transfer state are presented.

Catalytic synthesis of 9-cis-retinoids: mechanistic insights

The regioselective Z-isomerization of thermodynamically stable all-trans retinoids remains challenging, and ultimately limits the availability of much needed therapeutics for the treatment of human diseases. We present here a novel, straightforward approach for the catalytic Z-isomerization of retinoids using conventional heat treatment or microwave irradiation. A screen of 20 transition metal-based catalysts identified an optimal approach for the regioselective production of Z-retinoids. The most effective catalytic system was comprised of a palladium complex with labile ligands. Several mechanistic studies, including isotopic H/D exchange and state-of-the-art quantum chemical calculations using coupled cluster methods indicate that the isomerization is initiated by catalyst dimerization followed by the formation of a cyclic, six-membered chloropalladate catalyst-substrate adduct, which eventually opens to produce the desired Z-isomer. The synthetic development described here, combined with thorough mechanistic analysis of the underlying chemistry, highlights the use of readily available transition metal-based catalysts in straightforward formats for gram-scale drug synthesis.

Accurate Band Gap Predictions of Semiconductors in the Framework of the Similarity Transformed Equation of Motion Coupled Cluster Theory

In this work, we present a detailed comparison between wave-function-based and particle/hole techniques for the prediction of band gap energies of semiconductors. We focus on the comparison of the back-transformed Pair Natural Orbital Similarity Transformed Equation of Motion Coupled-Cluster (bt-PNO-STEOM-CCSD) method with Time Dependent Density Functional Theory (TD-DFT) and Delta Self Consistent Field/DFT (Δ-SCF/DFT) that are employed to calculate the band gap energies in a test set of organic and inorganic semiconductors. Throughout, we have used cluster models for the calculations that were calibrated by comparing the results of the cluster calculations to periodic DFT calculations with the same functional. These calibrations were run with cluster models of increasing size until the results agreed closely with the periodic calculation. It is demonstrated that bt-PNO-STEOM-CC yields accurate results that are in better than 0.2 eV agreement with the experiment. This holds for both organic and inorganic semiconductors. The efficiency of the employed computational protocols is thoroughly discussed. Overall, we believe that this study is an important contribution that can aid future developments and applications of excited state coupled cluster methods in the field of solid-state chemistry and heterogeneous catalysis.

In this work, it is shown that a combination of the embedded cluster approach with wave-function-based ab initio methods in the framework of the Similarity Transformed Equation of Motion Coupled Cluster (bt-PNO STEOM-CC) provides an accurate protocol for band gap energy predictions in classes of organic and inorganic semiconductors.

Chemistry and Quantum Mechanics in 2019: Give Us Insight and Numbers

This Perspective revisits Charles Coulson’s famous statement from 1959 “give us insight not numbers” in which he pointed out that accurate computations and chemical understanding often do not go hand in hand. We argue that today, accurate wave function based first-principle calculations can be performed on large molecular systems, while tools are available to interpret the results of these calculations in chemical language. This leads us to modify Coulson’s statement to “give us insight and numbers”. Examples from organic, inorganic, organometallic and surface chemistry as well as molecular magnetism illustrate the points made.

Reversible complexation of ammonia by breaking a manganese–manganese bond in a manganese carbonyl ethylenedithiolate complex: a theoretical study of an unusual type of Lewis acid

The addition of bases such as ammonia and trimethylphosphine to H₂C₂S₂Mn₂(CO)₆ to give yellow 1 : 1 adducts is shown to break the metal–metal bond rather than displace the coordinated double bond.

Ab initio design of light absorption through silver atomic cluster decoration of TiO2

A first-principles study of the stability and optical response of subnanometer silver clusters Agn (n ≤ 5) on a TiO2(110) surface is presented. First, the adequacy of the vdW-corrected DFT-D3 approach is assessed using the domain-based pair natural orbital correlation DLPNO-CCSD(T) calculations along with the Symmetry-Adapted Perturbation Theory [SAPT(DFT)] applied to a cluster model. Next, using the DFT-D3 treatment with a periodic slab model, we analyze the interaction energies of the atomic silver clusters with the TiO2(110) surface. Finally, the hybrid HSE06 functional and a reduced density matrix treatment are applied to obtain the projected electronic density of states and photo-absorption spectra of the TiO2(110) surface, with and without adsorbed silver clusters. Our results show the stability of the supported clusters, the enhanced light absorbance intensity of the material upon their deposition, and the appearance of intense secondary broad peaks in the near-infrared and the visible regions of the spectrum, with positions depending on the size and shape of the supported clusters. The secondary peaks arise from the photo-induced transfer of electrons from intra-band valence 5s orbitals of the noble-metal cluster to 3d Ti band states of the supporting material.

Addressing global uncertainty and sensitivity in first-principles based microkinetic models by an adaptive sparse grid approach

In the last decade, first-principles-based microkinetic modeling has been developed into an important tool for a mechanistic understanding of heterogeneous catalysis. A commonly known, but hitherto barely analyzed issue in this kind of modeling is the presence of sizable errors from the use of approximate Density Functional Theory (DFT). We here address the propagation of these errors to the catalytic turnover frequency (TOF) by global sensitivity and uncertainty analysis. Both analyses require the numerical quadrature of high-dimensional integrals. To achieve this efficiently, we utilize and extend an adaptive sparse grid approach and exploit the confinement of the strongly non-linear behavior of the TOF to local regions of the parameter space. We demonstrate the methodology on a model of the oxygen evolution reaction at the Co₃O₄ (110)-A surface, using a maximum entropy error model that imposes nothing but reasonable bounds on the errors. For this setting, the DFT errors lead to an absolute uncertainty of several orders of magnitude in the TOF. We nevertheless find that it is still possible to draw conclusions from such uncertain models about the atomistic aspects controlling the reactivity. A comparison with derivative-based local sensitivity analysis instead reveals that this more established approach provides incomplete information. Since the adaptive sparse grids allow for the evaluation of the integrals with only a modest number of function evaluations, this approach opens the way for a global sensitivity analysis of more complex models, for instance, models based on kinetic Monte Carlo simulations.