Mechanistic understanding of N2 activation: a comparison of unsupported and supported Ru catalysts.

N2 dissociative adsorption is commonly the rate-determining step in thermal ammonia synthesis. Herein, we performed density functional theory (DFT) calculations to understand the N2 dissociation mechanism on models of unsupported Ru(0001) terraces, Ru B5 sites, and polar MgO(111)-supported Ru8 cluster mimicking a B5 site geometry, denoted (Ru8(B5-like)/MgO(111)). The activation energy of N2 dissociative adsorption on the Ru8(B5-like)/MgO(111) model (Ea = 0.33 eV) is much lower than that on the unsupported Ru(0001) terrace (Ea = 1.74 eV) and Ru B5 (Ea = 0.62 eV) models. The lower N2 dissociation barrier on Ru B5 sites is facilitated by the enhanced σ donation and π* back-donation between N2(σ, π*) and Ru(d) orbitals resulting in the stronger activation of the molecular side-on N2* dissociation precursor. The Ru8(B5-like)/MgO(111) also exhibits enhanced σ donation because of the B5-like cluster geometry. Furthermore, the Ru cluster of the bare Ru8(B5-like)/MgO(111) model is positively charged. This induced an unusual π donation from N2(π) to Ru(d) orbitals as revealed by analyses of the density of states and partial charge densities. The combined σ and π donation resulted in an increased synergistic π* back-donation. The total interactions between N2(σ, π, π*) and Ru(d) resulted in an overall electron transfer to the adsorbed N2 from the Ru atoms in the B5-like site with no direct involvement of the MgO(111) substrate. Analyses of bond stretching vibrations and bond lengths show that the N2(σ, π, π*) and Ru(d) interactions lead to a weaker N-N bond and stronger Ru-N bonds. These correspond to a lower barrier of N2 dissociation on the Ru8(B5-like)/MgO(111) model, where the highest red-shift of N-N vibration and the longest N-N bond length were observed after side-on N2* adsorption. These results demonstrate that an electron-deficient Ru catalyst are not always inhibited from donating electrons to adsorbed N2. Rather, this study shows that the electron deficiency of Ru can promote π* back-donation and N2 activation. These new insights may therefore open new avenues to design supported Ru catalysts for nitrogen activation.

Cluster Nanozymes with Optimized Reactivity and Utilization of Active Sites for Effective Peroxidase (and Oxidase) Mimicking.

Single-atom catalysts have attracted attention in the past decade since they maximize the utilization of active sites and facilitate the understanding of product distribution in some catalytic reactions. Recently, this idea has been extended to single-atom nanozymes (SAzymes) for the mimicking of natural enzymes such as horseradish peroxidase (HRP) often used in bioanalytical applications. Herein, it is demonstrated that those SAzymes without constructing the reaction pocket of HRP still undergo the OH radical-mediated pathway like most of the reported nanozymes. Their positively charged single-atom centers resulting from support electronegative oxygen/nitrogen hinder the reductive conversion of H2 O2 to OH radicals and hence display low activity per site. In contrast, it is found that this step can be facilitated over their metallic counterparts on cluster nanozymes with much higher site activity and atom efficiency (cf. SAzymes with 100% atom utilization). Besides the mimicking of HRP in glucose detection, cluster nanozymes are also demonstrated as a better oxidase mimetic for glutathione detection.

Structural evolution of atomically dispersed Pt catalysts dictates reactivity.

The use of oxide-supported isolated Pt-group metal atoms as catalytic active sites is of interest due to their unique reactivity and efficient metal utilization. However, relationships between the structure of these active sites, their dynamic response to environments and catalytic functionality have proved difficult to experimentally establish. Here, sinter-resistant catalysts where Pt was deposited uniformly as isolated atoms in well-defined locations on anatase TiO2 nanoparticle supports were used to develop such relationships. Through a combination of in situ atomic-resolution microscopy- and spectroscopy-based characterization supported by first-principles calculations it was demonstrated that isolated Pt species can adopt a range of local coordination environments and oxidation states, which evolve in response to varied environmental conditions. The variation in local coordination showed a strong influence on the chemical reactivity and could be exploited to control the catalytic performance.

Electronic structure of Al, Ga, In and Cu doped ZnO/Cu(111) bilayer films.

The effect of doping with group-III metals (Al, Ga and In) and Cu free standing and Cu(111) supported ZnO bilayer films has been investigated computationally by using the DFT+U method including dispersion contributions. The changes in the electronic properties of doped ZnO and ZnO/Cu(111) films have been tested by adsorbing CO probe molecules. The replacement of a lattice Zn ion in a free-standing ZnO bilayer by a group-III element generates an extra electron whose distribution depends on the dopant. In particular, while the excess electron is delocalized over the conduction band for Al or Ga doping, it is localized on the dopant in the case of In. The situation is different on the supported ZnO/Cu(111) film, where the extra electron is transferred to the underlying Cu support. While the CO adsorption energy at the doped sites in the ZnO bilayer is the same as in the ZnO/Cu(111) ultrathin films, CO exhibits a larger red-shift in the unsupported ZnO bilayer. The oxidation state of Cu replacing Zn in the unsupported ZnO films is 2+, Cu(3d9) state, while it is 1+, Cu(3d10) state, in the ZnO/Cu(111) supported films where a charge transfer from the supporting Cu metal to the Cu impurity occurs. Cu doping results in a stronger interaction with CO and a large red-shift of the CO stretching frequency. In this respect, Cu doping of ZnO/Cu(111) bilayer films could have interesting consequences in gas adsorption while doping with group-III elements does not lead to major changes of the adsorption properties when the free-standing ZnO films are compared to the supported ZnO/Cu(111) counterparts.

The effect of the zeolite pore size on the Lewis acid strength of extra-framework cations.

The catalytic activity and the adsorption properties of zeolites depend on their topology and composition. For a better understanding of the structure-activity relationship it is advantageous to focus just on one of these parameters. Zeolites synthesized recently by the ADOR protocol offer a new possibility to investigate the effect of the channel diameter on the adsorption and catalytic properties of zeolites: UTL, OKO, and PCR zeolites consist of the same dense 2D layers (IPC-1P) that are connected with different linkers (D4R, S4R, O-atom, respectively) resulting in the channel systems of different sizes (14R × 12R, 12R × 10R, 10R × 8R, respectively). Consequently, extra-framework cation sites compensating charge of framework Al located in these dense 2D layers (channel-wall sites) are the same in all three zeolites. Therefore, the effect of the zeolite channel size on the Lewis properties of the cationic sites can be investigated independent of other factors determining the quality of Lewis sites. UTL, OKO, and PCR and pillared 2D IPC-1PI materials were prepared in Li-form and their properties were studied by a combination of experimental and theoretical methods. Qualitatively different conclusions are drawn for Li(+) located at the channel-wall sites and at the intersection sites (Li(+) located at the intersection of two zeolite channels): the Lewis acid strength of Li(+) at intersection sites is larger than that at channel-wall sites. The Lewis acid strength of Li(+) at channel-wall sites increases with decreasing channel size. When intersecting channels are small (10R × 8R in PCR) the intersection Li(+) sites are no longer stable and Li(+) is preferentially located at the channel-wall sites. Last but not least, the increase in adsorption heats with the decreasing channel size (due to enlarged dispersion contribution) is clearly demonstrated.

Measuring the Brønsted acid strength of zeolites--does it correlate with the O-H frequency shift probed by a weak base?

Brønsted-acid zeolites are currently being used as catalysts in a wide range of technological processes, spanning from the petrochemical industry to biomass upgrade, methanol to olefin conversion and the production of fine chemicals. For most of the involved chemical processes, acid strength is a key factor determining catalytic performance, and hence there is a need to evaluate it correctly. Based on simplicity, the magnitude of the red shift of the O-H stretching frequency, Δν(OH), when the Brønsted-acid hydroxyl group of protonic zeolites interacts with an adsorbed weak base (such as carbon monoxide or dinitrogen) is frequently used for ranking acid strength. Nevertheless, the enthalpy change, ΔH(0), involved in that hydrogen-bonding interaction should be a better indicator; and in fact Δν(OH) and ΔH(0) are often found to correlate among themselves, but, as shown herein, that is not always the case. We report on experimental determination of the interaction (at a low temperature) of carbon monoxide and dinitrogen with the protonic zeolites H-MCM-22 and H-MCM-56 (which have the MWW structure type) showing that the standard enthalpy of formation of OH···CO and OH···NN hydrogen-bonded complexes is distinctively smaller than the corresponding values reported in the literature for H-ZSM-5 and H-FER, and yet the corresponding Δν(OH) values are significantly larger for the zeolites having the MWW structure type (DFT calculations are also shown for H-MCM-22). These rather unexpected results should alert the reader to the risk of using the O-H frequency shift probed by an adsorbed weak base as a general indicator for ranking zeolite Brønsted acidity.

Computational investigation of the Lewis acidity in three-dimensional and corresponding two-dimensional zeolites: UTL vs IPC-1P.

The adsorption and catalytic properties of three-dimensional zeolite UTL were investigated computationally along with properties of its two-dimensional analogue IPC-1P that can be obtained from UTL by a removal of D4R units. Adsorption properties and Lewis acidity of extra-framework Li(+) sites were investigated for both two- and three-dimensional forms of UTL using the carbon monoxide as a probe molecule. The CO adsorption enthalpies, calculated with various dispersion-corrected DFT methods, including DFT/CC, DFT-D2, and vdW-DF2, and the CO stretching frequencies obtained with the νCO/rCO correlation method are compared for corresponding Li(+) sites in 3D and 2D UTL zeolite. For the majority of framework Al positions the Li(+) cation is preferably located in one of the channel wall sites and such sites remains unchanged upon the 3D → 2D UTL transformation; consequently, the adsorption enthalpies become only slightly smaller in 2D UTL (less than 3 kJ mol(-1)) due to the missing part of dispersion interactions and νCO becomes also only up to 5 cm(-1) smaller in 2D UTL due to the missing repulsion with framework oxygen atoms from the opposite site of the zeolite channel (effect from the top). However, when Li(+) is located in the intersection site in 3D UTL (about 20% probability), its coordination with the framework is significantly increased in 2D UTL and that is accompanied by significant decrease of both νCO (about 20 cm(-1)) and adsorption enthalpy (about 20 kJ mol(-1)). Because the intersection sites in 3D UTL are the most active adsorption and catalytic Lewis sites, the results reported herein suggest that the 3D → 2D transformation of UTL zeolite is connected with partial decrease of zeolite activity in processes driven by Lewis acid sites.

Dispersed surface Ru ensembles on MgO(111) for catalytic ammonia decomposition.

Ammonia is regarded as an energy vector for hydrogen storage, transport and utilization, which links to usage of renewable energies. However, efficient catalysts for ammonia decomposition and their underlying mechanism yet remain obscure. Here we report that atomically-dispersed Ru atoms on MgO support on its polar (111) facets {denoted as MgO(111)} show the highest rate of ammonia decomposition, as far as we are aware, than all catalysts reported in literature due to the strong metal-support interaction and efficient surface coupling reaction. We have carefully investigated the loading effect of Ru from atomic form to cluster/nanoparticle on MgO(111). Progressive increase of surface Ru concentration, correlated with increase in specific activity per metal site, clearly indicates synergistic metal sites in close proximity, akin to those bimetallic N2 complexes in solution are required for the stepwise dehydrogenation of ammonia to N2/H2, as also supported by DFT modelling. Whereas, beyond surface doping, the specific activity drops substantially upon the formation of Ru cluster/nanoparticle, which challenges the classical view of allegorically higher activity of coordinated Ru atoms in cluster form (B5 sites) than isolated sites.

Do stoichiometric or nonstoichiometric models of a polar surface affect their structural, energetic and electronic structure properties? A DFT case study of Ru/MgO(111).

The structural, energetic and electronic structure properties of stoichiometric and nonstoichiometric slab models of bare MgO(111) and Ru/MgO(111) with different coverages of 1 monolayer (ML), 1/4 ML and 1/16 ML have been investigated using spin-polarized density functional theory. Calculated results show that the structural, energetic properties and charge transfer of both bare MgO(111) and Ru/MgO(111) are independent of the stoichiometric and nonstoichiometric models. In contrast, their density of state (DOS) profiles demonstrate metal and half-metal characters for the stoichiometric and nonstoichiometric bare MgO(111) surfaces, respectively. The Ru-O orbital coupling characters of these two types of Ru/MgO(111) models are also different. This work indicates that for a polar surface model, the calculated features and trends of the structural and energetic properties, charge distributions and magnetic structures might not be affected by their stoichiometric and nonstoichiometric models; however, the detailed features of their DOS features would strongly depend on the models constructed.

Operando X-ray Studies of Ni-Containing Heteropolyvanadate Electrode for High-Energy Lithium-Ion Storage Applications.

Ni-containing heteropolyvanadate, Na6[NiV14O40], was synthesized for the first time to be applied in high-energy lithium storage applications as a negative electrode material. Na6[NiV14O40] can be prepared via a facile solution process that is suitable for low-cost mass production. The as-prepared electrode provided a high capacity of approximately 700 mAh g-1 without degradation for 400 cycles, indicating excellent cycling stability. The mechanism of charge storage was investigated using operando X-ray absorption spectroscopy (XAS), X-ray diffraction (XRD), transition X-ray microscopy (TXM), and density functional theory (DFT) calculations. The results showed that V5+ was reduced to V2+ during lithiation, indicating that Na6[NiV14O40] is an insertion-type material. In addition, Na6[NiV14O40] maintained its amorphous structure with negligible volume expansion/contraction during cycling. Employed as the negative electrode in a lithium-ion battery (LIB), the Na6[NiV14O40]//LiFePO4 full cell had a high energy density of 300 W h kg-1. When applied in a lithium-ion capacitor, the Na6[NiV14O40]//expanded mesocarbon microbead full cell displayed energy densities of 218.5 and 47.9 W h kg-1 at power densities of 175.7 and 7774.2 W kg-1, respectively. These findings reveal that the negative electrode material Na6[NiV14O40] is a promising candidate for Li-ion storage applications.

Electronic structure of CuWO4: dielectric-dependent, self-consistent hybrid functional study of a Mott-Hubbard type insulator.

CuWO4 is a semiconducting oxide with interesting applications in photocatalysis. In this paper we present an accurate study of the electronic properties of stoichiometric and oxygen deficient CuWO4 based on a dielectric dependent hybrid density functional. In CuWO4 the Cu ions (Cu2+) are in a 3d9 configuration, so that the material must be classified as a magnetic insulator. Various magnetic configurations of CuWO4 have been considered, the most stable configuration being anti-ferromagnetic. The band structure, described in terms of density of states (DOS), exhibit the presence of a wide band dominated by W 5d states, separated by about 5 eV from the top of the valence band (VB), consisting of O 2p states partly mixed with Cu 3d states. The empty component of the Cu 3d orbitals forms a narrow band 3.6 eV above the VB maximum. The electronic structure emerging from the DOS curves and the Kohn-Sham energies is hard to reconcile with an experimental band gap of 2.1-2.3 eV. This gap can be rationalized within the Mott-Hubbard model of magnetic insulators, and has been computed from the total energies of the system with one electron removed from the O 2p band and one electron added to the Cu 3d states. Computing the charge transition levels for CuWO4, we come to a theoretical band gap of 2.1 eV, in excellent agreement with the experimental observations. We also studied the nature of the oxygen vacancy in CuWO4 with particular attention to the electron redistribution following the oxygen removal. The excess electrons, in fact, can occupy the localized 3d states of Cu or the localized 5d states of W. The resulting solution depends on various factors, including the concentration of oxygen vacancies.

Theoretical investigation of layered zeolites with MWW topology: MCM-22P vs. MCM-56.

The inter-layer interactions and the possible arrangements of MWW-type layers were investigated computationally at the non-local density functional theory level. Powder XRD patterns were simulated for structures obtained computationally and compared with experimental data. The MCM-22P material corresponds to the layers bound with relatively strong hydrogen bonds between surface silanol groups that is an energetically preferred structure in the presence of a structure directing agent (hexamethyleneimine). The powder XRD pattern of MCM-56 is best matched for relatively disordered (in the ab plane) MWW layers that are partially condensed. The appearance of the powder XRD pattern in the 2θ range of 7.5-10° depends on the extent of interlayer condensation. The combination of density functional investigation of interactions between MWW layers together with simulation of powder XRD patterns brings atomistic insight into the inter-layer arrangement and better understanding of the effects responsible for the differences between various layered materials of the MWW family.