Comprehensive Benchmarking of a Density-Dependent Dispersion Correction

Uncovering the Activity of Alkaline Earth Metal Hydrogenation Catalysis Through Molecular Volcano Plots

Recent advances in alkaline earth (Ae) metal hydrogenation catalysis have broadened the spectrum of potential catalysts to include candidates from the main group, providing a sustainable alternative to the commonly used transition metals. Although Ae-amides have already been demonstrated to catalyze hydrogenation of imines and alkenes, a lucid understanding of how different metal/ligand combinations influence the catalytic activity is yet to be established. In this article, we use linear scaling relationships and molecular volcano plots to assess the potential of the Ae metal-based catalysts for the hydrogenation of alkenes. By analyzing combinations of eight metals (mono-, bi-, tri-, and tetravalent) and seven ligands, we delineate the impact of metal-ligand interplay on the hydrogenation activity. Our findings highlight that the catalytic activity is majorly determined by the charge and the size of the metal ions. While bivalent Ae metal cations delicately regulate the binding and the release of the reactants and the products, respectively, providing the right balance for this reaction, ligands play only a minor role in determining their catalytic activity. We show how volcano plots can be utilized for the rapid screening of prospective Ae catalysts to establish a guideline to achieve maximum activity in facilitating the hydrogenation process.

SUPPLEMENTARY INFORMATION

The online version of this article at 10.1007/s11244-021-01480-7.

Identification of active catalysts for the acceptorless dehydrogenation of alcohols to carbonyls

Acceptorless dehydrogenation into carbonyls and molecular hydrogen is an attractive strategy to valorize (biobased) alcohols. Using 2-octanol dehydrogenation as benchmark reaction in a continuous reactor, a library of metal-supported catalysts is tested to validate the predictive level of catalytic activity for combined DFT and micro-kinetic modeling. Based on a series of transition metals, scaling relations are determined as a function of two descriptors, i.e. the surface binding energies of atomic carbon and oxygen. Then, a volcano-shape relation based on both descriptors is derived, paving the way to further optimization of active catalysts. Evaluation of 294 diluted alloys but also a series of carbides and nitrides with the volcano map identified 12 promising candidates with potentially improved activity for alcohol dehydrogenation, which provides useful guidance for experimental catalyst design. Further screening identifies β-Mo2N and γ-Mo2N exposing mostly (001) and (100) facets as potential candidates for alcohol dehydrogenation.

Computational Prediction of Superlubric Layered Heterojunctions

Structural superlubricity has attracted increasing interest in modern tribology. However, experimental identification of superlubric interfaces among the vast number of heterojunctions is a trial-and-error and time-consuming approach. In this work, based on the requirements on the in-plane stiffnesses of layered materials and the interfacial interactions at the sliding incommensurate interfaces of heterojunctions for structural superlubricity, we propose criteria for predicting structural superlubricity between heterojunctions. Based on these criteria, we identify 61 heterojunctions with potential superlubricity features from 208 candidates by screening the data of first-principles calculations. This work provides a universal route for accelerating the discovery of new superlubric heterojunctions.

Requirements for an accurate dispersion-corrected density functional

Post-self-consistent dispersion corrections are now the norm when applying density-functional theory to systems where non-covalent interactions play an important role. However, there is a wide range of base functionals and dispersion corrections available from which to choose. In this work, we opine on the most desirable requirements to ensure that both the base functional and dispersion correction, individually, are as accurate as possible for non-bonded repulsion and dispersion attraction. The base functional should be dispersionless, numerically stable, and involve minimal delocalization error. Simultaneously, the dispersion correction should include finite damping, higher-order pairwise dispersion terms, and electronic many-body effects. These criteria are essential for avoiding reliance on error cancellation and obtaining correct results from correct physics.

Formation of a Ti-Cu(111) single atom alloy: Structure and CO binding

A single atom Ti-Cu(111) surface alloy can be generated by depositing small amounts of Ti onto Cu(111) at slightly elevated surface temperatures (∼500 to 600 K). Scanning tunneling microscopy shows that small Ti-rich islands covered by a Cu single layer form preferentially on ascending step edges of Cu(111) during Ti deposition below about 400 K but that a Ti-Cu(111) alloy replaces these small islands during deposition between 500 and 600 K, producing an alloy in the brims of the steps. Larger partially Cu-covered Ti-containing islands also form on the Cu(111) terraces at temperatures between 300 and 700 K. After surface exposure to CO at low temperatures, reflection absorption infrared spectroscopy (RAIRS) reveals distinct C-O stretch bands at 2102 and 2050 cm^-1 attributed to CO adsorbed on Cu-covered Ti-containing domains vs sites in the Ti-Cu(111) surface alloy. Calculations using density functional theory (DFT) suggest that the lower frequency C-O stretch band originates specifically from CO adsorbed on isolated Ti atoms in the Ti-Cu(111) surface alloy and predicts a higher C-O stretch frequency for CO adsorbed on Cu above subsurface Ti ensembles. DFT further predicts that CO preferentially adsorbs in flat-lying configurations on contiguous Ti surface structures with more than one Ti atom and thus that CO adsorbed on such structures should not be observed with RAIRS. The ability to generate a single atom Ti-Cu(111) alloy will provide future opportunities to investigate the surface chemistry promoted by a representative early transition metal dopant on a Cu(111) host surface.

Judicious design functionalized 3D-COF to enhance CO2 adsorption and separation

The effects of functional groups (including OH, OCH₃ , NH₂ , CH₂ NH₂ , COOH, SO₃ H, OCO(CH₂ )₂ COOH(E-COOH), and (CH₂ )₄ COOH(c-COOH)) in 3D covalent organic frameworks (3D-COFs) on CO₂ adsorption and separation are investigated by grand canonical Monte Carlo (GCMC) simulations and density functional theory calculations. The results indicate that interaction between CO₂ and the framework is the main factor for determining CO₂ uptakes at low pressure, while pore size becomes the decisive factor at high pressure. The binding energy of CO₂ with functionalized linker is correlated to CO₂ uptake at 0.3 bar and 298 K on 3D-COF-1, suggesting functional groups play a key role in CO₂ capture in microporous 3D-COFs. Moreover, CO₂ selectivity over CH₄ , N₂ , and H₂ can be significantly enhanced by functionalization, where CH₂ NH₂ , COOH, SO₃ H, and E-COOH enhance CO₂ adsorption more effectively at 1 bar. Among them, SO₃ H is the most promising functional group in 3D-COFs for CO₂ separation.

FeCoP2 Nanoparticles Embedded in N and P Co-doped Hierarchically Porous Carbon for Efficient Electrocatalytic Water Splitting

The design and synthesis of low-cost and efficient bifunctional electrocatalysts for water splitting are critical and challenging. Hereby, a bimetallic phosphide embedded in a N and P co-doped porous carbon (FeCoP₂@NPPC) material was synthesized by using sustainable biomass-derived N- and P-containing carbohydrates and non-noble metal salts as precursors. The obtained material exhibits good catalytic activities in hydrogen evolution reaction (HER), oxygen evolution reaction (OER), and overall water splitting. The bimetallic alloy phosphide (FeCoP₂) is the active site for electrocatalysis. Theoretical calculation indicates that the sub-layer Fe atoms and top-layer Co atoms in FeCoP₂ exhibit a synergistic effect for enhanced electrocatalytic performance. The carbon matrix around the FeCoP₂ can prevent the corrosion during the catalytic reactions. The hierarchically porous structure of the FeCoP₂@NPPC material can promote the transfer of electrons and electrolyte, and increase the contact area of the active sites and electrolytes. N- and P-containing functionalities improve the wetting and conductivity properties of the porous carbon. Due to the synergistic effects, FeCoP₂@NPPC requires a low overpotential of 114 and 150 mV at the current density of 10 mA cm^-2 for HER in 0.5 M H₂SO₄ and 1.0 M KOH, and an overpotential of 236 mV for OER in 1.0 M KOH solution, which are much lower than those of FeP@NPPC and CoP@NPPC. Based on the density functional theory calculation, FeCoP₂ yields the smallest Gibbs free energy change of rate-determining step among the samples, which leads to better electrochemical performances. In addition, when FeCoP₂@NPPC was used as a bifunctional catalyst in water splitting, the electrolyzer needed a low voltage of 1.60 V to deliver the current density of 10 mA cm^-2. Furthermore, this work provides a strategy for preparing sustainable, stable, and highly active electrocatalysts for water splitting.

What does graphitic carbon nitride really look like?

Graphitic carbon nitrides (g-CNs) have become popular light absorbers in photocatalytic water splitting cells. Early theoretical work on these structures focused on fully polymerized g-C3N4. Experimentally, it is known that the typically employed melamine polycondensation does not go toward completion, yielding structures with ∼15 at% hydrogen. Here, we study the conformational stability of "melon", with the [C6N9H3]n structural formula using DFT. Referencing to a 2D melon sheet, B3LYP-dDsC and PBE-MBD computations revealed the same qualitative trend in stability of the 3D structures, with several of them within 5 kJ mol-1 per tecton. Fina's orthorhombic melon is the most stable of the studied conformers, with Lotsch' monoclinic melon taking an intermediate value. Invoking a simple Wannier-Mott-type approach, Fina's and Lotsch' structures exhibited the lowest optical gaps (2.8 eV), within the error margin of the experimental value (2.7 eV). All conformers yielded gaps below that of the monolayer's (3.2 eV), suggesting Jelley-type ("J") aggregation effects.

DFT benchmark studies on representative species and poisons of methane steam reforming on Ni(111)

Ni catalysts used in methane steam reforming (MSR) are highly susceptible to poisoning by carbon-based species, which poses a major impediment to the productivity of industrial operations. These species encompass graphitic carbon-like formations that are typically modelled as graphene. First principles-based approaches, such as density functional theory (DFT) calculations, can provide valuable insight into the mechanism of graphene growth in the MSR reaction. It is, however, critical that a DFT model of this reaction can accurately describe the interactions of Ni(111) with the MSR intermediates as well as graphene. In this work, a systematic benchmark study has been carried out to identify a suitable DFT functional for the graphene-MSR system. The binding energies of graphene and important MSR species, as well as the reaction energies of methane dissociation and carbon oxidation, were computed on Ni(111) using GGA functionals, DFT-D and van der Waals density functionals (vdW-DF). It is well-established that the GGA functionals are inadequate for describing graphene-Ni(111) interactions. In the case of vdW-DF, the optPBE-vdW functional predicts the binding energies of graphene and several important MSR species with reasonable accuracy; however, it provides poor estimates of CO and O binding energies. Among the DFT-D functionals, PBE-D3 and PBE-dDsC have been found to exhibit acceptable accuracy for graphene and most MSR species (excluding adsorbed CO), and therefore, both functionals are promising for elucidating carbon-based catalytic poisoning in the MSR reaction. Overall, no single DFT functional could estimate the binding energies of all the species with equally high accuracy.

Computational Methods in Heterogeneous Catalysis

The unprecedented ability of computations to probe atomic-level details of catalytic systems holds immense promise for the fundamentals-based bottom-up design of novel heterogeneous catalysts, which are at the heart of the chemical and energy sectors of industry. Here, we critically analyze recent advances in computational heterogeneous catalysis. First, we will survey the progress in electronic structure methods and atomistic catalyst models employed, which have enabled the catalysis community to build increasingly intricate, realistic, and accurate models of the active sites of supported transition-metal catalysts. We then review developments in microkinetic modeling, specifically mean-field microkinetic models and kinetic Monte Carlo simulations, which bridge the gap between nanoscale computational insights and macroscale experimental kinetics data with increasing fidelity. We finally review the advancements in theoretical methods for accelerating catalyst design and discovery. Throughout the review, we provide ample examples of applications, discuss remaining challenges, and provide our outlook for the near future.

A Comprehensive Study of the Bridge Site and Substrate Relaxation Asymmetry for Methanethiol Adsorption on Au(111) at Low Coverage

We use dispersion-corrected density functional theory to explore the bridge-site asymmetry for methanethiol adsorbed on Au(111) with two different S-C bond orientations. We attribute the asymmetry to the intrinsic character of the Au(111) surface rather than the adsorbate. The preference for bridge-fcc versus bridge-hcp SCH3 adsorption sites is controlled by the S-C bond orientation. The system energy difference favors the bridge-fcc site by 8.1 meV on the unrelaxed Au(111) surface. Relaxing the Au substrate increased this energy difference to 26.1 meV. This asymmetry is also reflected in the atomic displacement of the relaxed Au surface. Although in both cases, the bridge-site Au atoms shift away from the fcc 3-fold hollow site, the motion is greater for the bridge-fcc allowing a more favorable geometry for the sulfur atom to bond to the bridging atoms. We confirm that the adsorption energy is strongly dependent on the S-C bond orientation and position, which can be understood in terms of a simple coordination geometry model. This work has important implications for alkanethiol surface diffusion and the structure of their self-assembled monolayers.

Application of XDM to ionic solids: The importance of dispersion for bulk moduli and crystal geometries

Dispersion corrections are essential in the description of intermolecular interactions; however, dispersion-corrected functionals must also be transferrable to hard solids. The exchange-hole dipole moment (XDM) model has demonstrated excellent performance for non-covalent interactions. In this article, we examine its ability to describe the relative stability, geometry, and compressibility of simple ionic solids. For the specific cases of the cesium halides, XDM-corrected functionals correctly predict the energy ranking of the B1 and B2 forms, and a dispersion contribution is required to obtain this result. Furthermore, for the lattice constants of the 20 alkali halides, the performance of XDM-corrected functionals is excellent, provided that the base functional's exchange enhancement factor properly captures non-bonded repulsion. The mean absolute errors in lattice constants obtained with B86bPBE-XDM and B86bPBE-25X-XDM are 0.060 Å and 0.039 Å, respectively, suggesting that delocalization error also plays a minor role in these systems. Finally, we considered the calculation of bulk moduli for alkali halides and alkaline-earth oxides. Previous claims in the literature that simple generalized gradient approximations, such as PBE, can reliably predict experimental bulk moduli have benefited from large error cancellations between neglecting both dispersion and vibrational effects. If vibrational effects are taken into account, dispersion-corrected functionals are quite accurate (4 GPa-5 GPa average error), again, if non-bonded repulsion is correctly represented. Careful comparisons of the calculated bulk moduli with experimental data are needed to avoid systematic biases and misleading conclusions.

Phase-controllable growth of ultrathin 2D magnetic FeTe crystals

Two-dimensional (2D) magnets with intrinsic ferromagnetic/antiferromagnetic (FM/AFM) ordering are highly desirable for future spintronic devices. However, the direct growth of their crystals is in its infancy. Here we report a chemical vapor deposition approach to controllably grow layered tetragonal and non-layered hexagonal FeTe nanoplates with their thicknesses down to 3.6 and 2.8 nm, respectively. Moreover, transport measurements reveal these obtained FeTe nanoflakes show a thickness-dependent magnetic transition. Antiferromagnetic tetragonal FeTe with the Néel temperature (TN) gradually decreases from 70 to 45 K as the thickness declines from 32 to 5 nm. And ferromagnetic hexagonal FeTe is accompanied by a drop of the Curie temperature (TC) from 220 K (30 nm) to 170 K (4 nm). Theoretical calculations indicate that the ferromagnetic order in hexagonal FeTe is originated from its concomitant lattice distortion and Stoner instability. This study highlights its potential applications in future spintronic devices.

Revisiting the Active Sites at the MoS2/H2O Interface via Grand-Canonical DFT: The Role of Water Dissociation

MoS₂ is a promising low-cost catalyst for the hydrogen evolution reaction (HER). However, the nature of the active sites remains a subject of debate. By taking the electrochemcal potential explicitly into account using grand-canonical density functional theory (DFT) in combination with the linearized Poisson-Boltzmann equation, we herein revisit the active sites of 2H-MoS₂. In addition to the well-known catalytically active edge sites, also specific point defects on the otherwise inert basal plane provide highly active sites for HER. Given that HER takes place in water, we also assess the reactivity of these active sites with respect to H₂O. The thermodynamics of proton reduction as a function of the electrochemical potential reveals that four edge sites and three basal plane defects feature thermodynamic overpotentials below 0.2 V. In contrast to current proposals, many of these active sites involve adsorbed OH. The results demonstrate that even though H₂O and OH block "active" sites, HER can also occur on these "blocked" sites, reducing protons on surface OH/H₂O entities. As a consequence, our results revise the active sites, highlighting the so far overlooked need to take the liquid component (H₂O) of the functional interface into account when considering the stability and activity of the various active sites.

Origin of ferromagnetism and the effect of doping on Fe3GeTe2

Recent experimental findings of two dimensional ferromagnetism in Fe3GeTe2, whose critical temperature can reach room temperature by gating, has attracted great research interest. Here we performed elaborate ab initio studies using density functional theory, dynamical mean-field theory and magnetic force response theory. In contrast to the conventional wisdom, it is unambiguously shown that Fe3GeTe2 is not ferromagnetic but is antiferromagnetic, carrying zero net moment in its stoichiometric phase. Fe defect and hole doping are the keys to make this material ferromagnetic as supported by previously disregarded experiments. Furthermore, we found that electron doping also induces the antiferro- to ferro-magnetic transition. It is crucial to understand the notable recent experiments on gate-controlled ferromagnetism. Our results not only reveal the origin of ferromagnetism of this material but also show how it can be manipulated with defects and doping.

Understanding Solid-Gas Reaction Mechanisms by Operando Soft X-Ray Absorption Spectroscopy at Ambient Pressure

Ambient-pressure operando soft X-ray absorption spectroscopy (soft-XAS) was applied to study the reactivity of hydroxylated SnO₂ nanoparticles toward reducing gases. H₂ was first used as a test case, showing that the gas phase and surface states can be simultaneously probed: Soft-XAS at the O K-edge gains sensitivity toward the gas phase, while at the Sn M_4,5-edges, tin surface states are explicitly probed. Results obtained by flowing hydrocarbons (CH₄ and CH₃CHCH₂) unequivocally show that these gases react with surface hydroxyl groups to produce water without producing carbon oxides and release electrons that localize on Sn to eventually form SnO. The partially reduced SnO_2 - x layer at the surface of SnO₂ is readily reoxidized to SnO₂ by treating the sample with O₂ at mild temperatures (>200 °C), revealing the nature of "electron sponge" of tin oxide. The experiments, combined with DFT calculations, allowed devising of a mechanism for dissociative hydrocarbon adsorption on SnO₂, involving direct reduction of Sn sites at the surface via cleavage of C-H bonds and the formation of methoxy- and/or methyl-tin species at the surface.

Acceptorless dehydrogenation and hydrogenation of N- and O-containing compounds on Pd3Au1(111) facets

Catalytic dehydrogenation and hydrogenation of amines and alcohols are important in the synthesis of fine chemicals. Despite several efficient homogeneous catalysts having been identified, highly active heterogeneous catalysts remain elusive, although they would meet an unmet need. Here, we show that bimetallic Pd-Au nanoparticles with Pd-to-Au molar ratios of 3:1 immobilized on multiwall carbon nanotubes (Pd3Au1/CNT) display high catalytic activity in the oxidant-free and acceptorless dehydrogenation and hydrogenation of N- and O-containing heterocyclic compounds, amines/imines, and alcohols/ketones. Transmission electron microscopy analysis demonstrates the preferential exposure of Pd3Au1(111) facets on the Pd3Au1/CNT catalyst. Mechanistic insights combining experimental data with density functional theory calculations reveal that the Pd3Au1(111) surface enhances both dehydrogenation and hydrogenation reactions and provides a rationale for the observed enhancements.

Ten Facets, One Force Field: The GAL19 Force Field for Water-Noble Metal Interfaces

Understanding the structure of the water/metal interfaces plays an important role in many areas ranging from surface chemistry to environmental processes. The size, required phase-space sampling, and the slow diffusion of molecules at the water/metal interfaces motivate the development of accurate force fields. We develop and parametrize GAL19, a novel force field, to describe the interaction of water with two facets (111 and 100) of five metals (Pt, Pd, Au, Ag, Cu). To increase transferability compared to its predecessor GAL17, the water-metal interaction is described as a sum of pairwise terms. The interaction energy has three contributions: (i) physisorption is described via a Tang and Toennies potential, (ii) chemisorption and surface corrugation rely on an attractive Gaussian term, and (iii) the angular dependence is explicitly included as a truncated Fourier series. Thirteen parameters are used for each metal surface and were fitted on 250 water adsorption energies computed at the PBE+dDsC level. The performance of GAL19 was evaluated on a set of more than 600 DFT adsorption energies for each surface, leading to an average root-mean-square deviation of only 1 kcal/mol, correctly reproducing the adsorption trends: strong on Pt and Pd but weaker on Ag, Au, and Cu. This force field was then used to simulate the water/metal interface for all ten surfaces for 1 ns. Structural analyses reveal similar tendencies for all surfaces: a first, dense water layer that is mostly adsorbed on the metal top sites and a second layer up to around 6 Å, which is less structured. On Pt and Pd, the first layer is strongly organized with water lying flat on the surface. The pairwise additive functional form allows one to simulate the water adsorption on alloys, which is demonstrated at the example of Ag/Cu and Au/Pt alloys. The water/Ag-Cu interface is predicted to be disordered with water mostly adsorbed on Cu which should exacerbate the Ag reactivity. On the contrary, incorporating Pt into Au materials leads to a structuring of the water interface. Our promising results make GAL19 an ideal candidate to get representative sampling of complex metal/water interfaces as a first step toward accurate estimation of free energies of reactions in solution at the metal interface.

Dioxygen dissociation over man-made system at room temperature to form the active α-oxygen for methane oxidation

Activation of dioxygen attracts enormous attention due to its potential for utilization of methane and applications in other selective oxidation reactions. We report a cleavage of dioxygen at room temperature over distant binuclear Fe(II) species stabilized in an aluminosilicate matrix. A pair of formed distant α-oxygen species [i.e., (Fe(IV)═O)²⁺] exhibits unique oxidation properties reflected in an outstanding activity in the oxidation of methane to methanol at room temperature. Designing a man-made system that mimicks the enzyme functionality in the dioxygen activation using both a different mechanism and structure of the active site represents a breakthrough in catalysis. Our system has an enormous practical importance as a potential industrial catalyst for methane utilization because (i) the Fe(II)/Fe(IV) cycle is reversible, (ii) the active Fe centers are stable under the reaction conditions, and (iii) methanol can be released to gas phase without the necessity of water or water-organic medium extraction.

Two-sites are better than one: revisiting the OER mechanism on CoOOH by DFT with electrode polarization

We uncover the existence of several competitive mechanisms of water oxidation on the β-CoOOH (10-14) surface by going beyond the classical 4-step mechanism frequently used to study this reaction at the DFT level. Our results demonstrate the importance of two-site reactivity and of purely chemical steps with the associated activation energies. Taking the electrochemical potential explicitly into account leads to modifications of the reaction energy profiles finally leading to the proposition of a new family of mechanisms involving tetraoxidane intermediates. The two-site mechanisms revealed in this work are of key importance to rationalize and predict the impact of dopants in the design of future catalysts.

Parameter-free coordination numbers for solutions and interfaces

Coordination numbers are among the central quantities to describe the local environment of atoms and are thus used in various applications such as structure analysis, fingerprints, and parameters. Yet, there is no consensus regarding a practical algorithm, and many proposed methods are designed for specific systems. In this work, we propose a scale-free and parameter-free algorithm for nearest neighbor identification. This algorithm extends the powerful Solid-Angle based Nearest-Neighbor (SANN) framework to explicitly include local anisotropy. As such, our Anisotropically corrected SANN (ASANN) algorithm provides with a fast, robust, and adaptive method for computing coordination numbers. The ASANN algorithm is applied to flat and corrugated metallic surfaces to demonstrate that the expected coordination numbers are retrieved without the need for any system-specific adjustments. The same applies to the description of the coordination numbers of metal atoms in AuCu nanoparticles, and we show that ASANN based coordination numbers are well adapted for automatically counting neighbors and the establishment of cluster expansions. Analysis of classical molecular dynamics simulations of an electrified graphite electrode reveals a strong link between the coordination number of Cs⁺ ions and their position within the double layer, a relation that is absent for Na⁺, which keeps its first solvation shell even close to the electrode.

Including dispersion in density functional theory for adsorption on flat oxide surfaces, in metal-organic frameworks and in acidic zeolites

We examine the performance of nine commonly used methods for including dispersion interactions in density functional theory (DFT): three different parametrizations of damped 1/Rn terms (n = 6, 8, …) added to the DFT energy (Grimme's D2 and D3 parameterizations as well as that of Tkatchenko and Scheffler), three different implementations of the many-body dispersion approach (MBD, MBD/HI and MBD/FI), the density-dependent energy correction, called dDsC, and two "first generation" van der Waals density functionals, revPBE-vdW and optB86b-vdW. As test set we use eight molecule-surface systems for which agreement has been reached between experiment and hybrid QM:QM calculations within chemical accuracy limits (±4.2 kJ mol-1). It includes adsorption of carbon monoxide and dioxide in the Mg2(2,5-dioxido-1,4-benzenedicarboxylate) metal-organic framework (Mg-MOF-74, CPO-27-Mg), adsorption of carbon monoxide as well as of monolayers of methane and ethane on the MgO(001) surface, as well as adsorption of methane, ethane and propane in H-chabazite (H-CHA). D2 with Ne parameters for Mg2+, D2(Ne), MBD/HI and MBD/FI perform best. With the PBE functional, the mean unsigned errors are 6.1, 5.6 and 5.4 kJ mol-1, respectively.

Impact of linker functionalization on the adsorption of nitrogen-containing compounds in HKUST-1

Functionalization of metal-organic framework (MOF) ligands can tune the adsorption properties of MOFs. The adsorptions of NO, NO2, NH3, C5H5N, C4H5N, and C4H4O on pristine and five X-functionalized HKUST-1, i.e. Cu3(BTC)2 (BTC = 1,3,5-benzenetricarboxylate) (X = CH3, CH3O, NH2, NO2, and Br) are evaluated by van der Waals corrected density functional theory calculations. Despite the fact that the open metal center is the energetically preferred adsorption site for most of them, the functional group site can yield a comparable adsorption ability with the open metal center. This is particularly true for pyrrole C4H5N adsorption on CH3O-functionalized HKUST-1 where the functional group site shows stronger adsorption stability than the open metal center site, probably due to the formed hydrogen bond between pyrrole and the CH3O functional group. While the CH3- or CH3O-functionalized organic linker in these MOFs strengthens the adsorption of all the considered species, that of NO2-, Br-, or NH2-functional groups reduces, which is associated with their topologies. Among them, only CH3- or CH3O-functionalized HKUST-1 presents the fmj (orthorhombic crystal system) topology while all the others are isostructural to the pristine HKUST-1 with the tbo (twisted boracite-type, cubic) topological structure. Among six adsorbates, two basic adsorbates, C5H5N and NH3, always yield the strongest bonding strength upon adsorption on the pristine and five functionalized HKUST-1. Electronic properties including the Bader charges, electron density differences, and electron localization function were investigated to comprehend their adsorption behaviors. This work provides guidance for the proper functionalization of HKUST-1 with improved adsorption properties for specific adsorbates.

Mononuclear Fe in N-doped carbon: computational elucidation of active sites for electrochemical oxygen reduction and oxygen evolution reactions

First principles simulations show that in Fe and N co-doped carbon, Fe coordination controls the activity for oxygen reduction and oxygen evolution reactions, and that including the electrostatic potential has a major influence at high potential.

Importance of the decoration in shaped cobalt nanoparticles in the acceptor-less secondary alcohol dehydrogenation

Ligands matter for shaped decorated Co nanoparticles, at the frontier between homogeneous and heterogeneous catalysis.

Functional Models of the Nickel Pincer Nucleotide Cofactor of Lactate Racemase

A novel nickel pincer cofactor was recently discovered in lactate racemase. Reported here are three synthetic nickel pincer complexes that are both structural and functional models of the pincer cofactor in lactate racemase. DFT computations suggest the ipso-carbon atom of the pyridinium pincer ligands act as a hydride acceptor for lactate isomerization, whereas an organometallic pathway involving nickel-mediated β-hydride elimination is less favored.

Mechanistic Study on the Photogeneration of Hydrogen by Decamethylruthenocene

Detailed studies on hydrogen evolution by decamethylruthenocene ([Cp*₂ Ru^II ]) highlighted that metallocenes are capable of photoreducing hydrogen without the need for an additional sensitizer. Electrochemical, gas chromatographic, and spectroscopic (UV/Vis, ¹ H and ¹³ C NMR) measurements corroborated by DFT calculations indicated that the production of hydrogen occurs by a two-step process. First, decamethylruthenocene hydride [Cp*₂ Ru^IV (H)]⁺ is formed in the presence of an organic acid. Subsequently, [Cp*₂ Ru^IV (H)]⁺ is reversibly reduced in a heterolytic reaction with one-photon excitation leading to a first release of hydrogen. Thereafter, the resultant decamethylruthenocenium ion [Cp*₂ Ru^III ]⁺ is further reduced with a second release of hydrogen by deprotonation of a methyl group of [Cp*₂ Ru^III ]⁺ . Experimental and computational data show spontaneous conversion of [Cp*₂ Ru^II ] to [Cp*₂ Ru^IV (H)]⁺ in the presence of protons. Calculations highlight that the first reduction is endergonic (ΔG⁰ =108 kJ mol^-1 ) and needs an input of energy by light for the reaction to occur. The hydricity of the methyl protons of [Cp*₂ Ru^II ] was also considered.

DFT and spatial confinement: a benchmark study on the structural and electrical properties of hydrogen bonded complexes

An extended set of 37 exchange correlation functionals, representing different DFT approximations, has been evaluated on a difficult playground represented by the dipole moment (μ_z), polarizability (α_zz), first hyperpolarizability (β_zzz), and the corresponding interaction-induced electrical properties (Δμ_z, Δα_zz, Δβ_zzz) of spatially confined hydrogen bonded (HB) dimers. A two-dimensional harmonic oscillator potential was used to exert the effect of spatial restriction. The performance of DFT methods in predicting hydrogen bond lengths in the studied molecular complexes upon confinement has also been examined. The data determined using a high-level CCSD(T) method serve as a reference. The conducted analyses allow us to conclude that methods rooted in DFT constitute a precise tool for the calculation of μ_z and α_zz as well as Δμ_z and Δα_zz, as most of the tested functionals provide results affected by rather small relative errors. On the other hand, an accurate description of the nonlinear optical response of the studied HB systems remains a great challenge for most of the analyzed DFT functionals, both in vacuum and in the presence of an analytical confining potential. Some of the tested DFT methods are found to be prone to catastrophic failure in the prediction of β_zzz as well as Δβ_zzz. The obtained results indicate that there is no great chasm in performance between functionals belonging to different DFT approximations or functionals including different amount of Hartree-Fock exchange when the values of dipole moment and first hyperpolarizability as well as the corresponding interaction-induced electrical properties are considered. However, a higher fraction of Hartree-Fock exchange improves the quality of predictions of α_zz and Δα_zz. Additionally, it has been shown that only three functionals from the examined set, namely B2PLYP, B3LYP and ωB97X-D, provide highly accurate structural parameters for the investigated systems. Of significant importance is the conclusion that the ωB97X-D functional, representing a modern and highly parametrized range-separated hybrid, demonstrates the most coherent behavior, showing rather small deviations from the reference data in the case of μ_z, α_zz, Δμ_z and Δα_zz as well as the structural parameters of the studied HB dimers. Moreover, our results indicate that the presence of spatial confinement has a rather small effect on the performance of DFT methods.