A comprehensive DFT investigation of bulk and low-index surfaces of ZrO2polymorphs

Low-nuclearity CuZn ensembles on ZnZrOx catalyze methanol synthesis from CO2

Metal promotion could unlock high performance in zinc-zirconium catalysts, ZnZrOx, for CO2 hydrogenation to methanol. Still, with most efforts devoted to costly palladium, the optimal metal choice and necessary atomic-level architecture remain unclear. Herein, we investigate the promotion of ZnZrOx catalysts with small amounts (0.5 mol%) of diverse hydrogenation metals (Re, Co, Au, Ni, Rh, Ag, Ir, Ru, Pt, Pd, and Cu) prepared via a standardized flame spray pyrolysis approach. Cu emerges as the most effective promoter, doubling methanol productivity. Operando X-ray absorption, infrared, and electron paramagnetic resonance spectroscopic analyses and density functional theory simulations reveal that Cu0 species form Zn-rich low-nuclearity CuZn clusters on the ZrO2 surface during reaction, which correlates with the generation of oxygen vacancies in their vicinity. Mechanistic studies demonstrate that this catalytic ensemble promotes the rapid hydrogenation of intermediate formate into methanol while effectively suppressing CO production, showcasing the potential of low-nuclearity metal ensembles in CO2-based methanol synthesis.

MC design and FIB preparation of a YSZ biochemical material microstructure

Aiming at the difficulty of traditional machining of Y2O3-ZrO2 (YSZ) inert ceramic materials, a different method using focused ion beam to selectively create nanoscale microscale structures on the surface of materials was proposed. The sputtering yield, surface damage, and the energy loss of YSZ materials was investigated using the SRIM software using the Monte Carlo method. It is shown that the sputtering yield increases with ion energy in the range 0-30 keV, reaching a maximum of 9.4 atoms/ion at 30 keV. At an ion beam voltage of 30 keV, the most severe damage to the material is 8 mm on the surface. At the same time, the main forms of energy loss in the treatment are phonon energy loss and ionization energy loss, of which phonon energy loss due to the recoil atoms is the largest. In addition, we continue to perform focused ion beam processing experiments on YSZ materials, combining previous MC modeling to optimize different operating conditions such as ion beam, voltage and processing mode. The optimized processing parameters are 30 keV and 2.5 nA. It is shown that the quality of the deep grooves gradually improves with decreasing ion beam current at the same ion beam voltage. However, an excessively small ion beam current leads to an excessively large depth of the deep grooves and lengthy processing times.

Doping induced phase stabilization and electronic properties of alkaline earth metal doped zirconium (IV) oxide: A first principles study

The role of divalent dopant cations such as Ca and Mg in phase stabilization of ZrO2 has been demonstrated experimentally, with Mg emerging as a crucial dopant ion because of its ability to enhance the photocatalytic properties of ZrO2. To provide a theoretical basis for these experimental observations, the modifications of the crystal and electronic structure of the monoclinic phase of zirconia, m-ZrO2, upon doping with Mg have been studied at the atomic level using Density Functional Theory method. Additionally, the effect of dopant ionic radius on the electronic properties has been demonstrated by doping with Ca, which is isoelectronic with Mg. On 6.25 % doping, a structural distortion of the monoclinic crystal structure towards a tetragonal structure is observed. Additionally, the Density of States of doped m-ZrO2 exhibits the characteristics of t-ZrO2 in the Zr d orbitals in the unoccupied states and O unoccupied states emerge upon creation of an O vacancy in Mg/Ca doped m-ZrO2. The calculated band gap of m-ZrO2 is 3.6 eV. Upon doping there is a shift of the Fermi energy towards the valence band maximum.

Control of monomeric Vo's versus Vo clusters in ZrO2-x for solar-light H2 production from H2O at high-yield (millimoles gr-1 h-1)

Pristine zirconia, ZrO₂, possesses high premise as photocatalyst due to its conduction band energy edge. However, its high energy-gap is prohibitive for photoactivation by solar-light. Currently, it is unclear how solar-active zirconia can be designed to meet the requirements for high photocatalytic performance. Moreover, transferring this design to an industrial-scale process is a forward-looking route. Herein, we have developed a novel Flame Spray Pyrolysis process for generating solar-light active nano-ZrO_2−x via engineering of lattice vacancies, Vo. Using solar photons, our optimal nano-ZrO_2−x can achieve milestone H₂-production yield, > 2400 μmolg⁻¹ h⁻¹ (closest thus, so far, to high photocatalytic water splitting performance benchmarks). Visible light can be also exploited by nano-ZrO_2−x at a high yield via a two-photon process. Control of monomeric Vo versus clusters of Vo’s is the key parameter toward Highly-Performing-Photocatalytic ZrO_2−x. Thus, the reusable and sustainable ZrO_2−x catalyst achieves so far unattainable solar activated photocatalysis, under large scale production.

Bishop's hat silicene: a planar square silicon bilayer decorated with adatoms

We investigate a family of free-standing quasi-two-dimensional silicon structures based on a planar square bilayer with adatom decorations. When attached to the bilayer, these adatoms form local reconstructions which resemble either a bishop's hat or elongated square bipyramids. We systematically constructed members of this family via exhaustive enumeration and then studied them using tight-binding and density-functional theory. We find that this geometry contributes significantly to the stability of the resulting structures, with some squared bilayers energetically more stable than the honeycomb bilayers. The most interesting phases were then characterized in more detail, and they all turned out metallic. Finally, we propose the [100] surface of ZrO2 as the most suitable substrate for the synthesis of these two-dimensional phases.

Development of Density-Functional Tight-Binding Parameters for the Molecular Dynamics Simulation of Zirconia, Yttria, and Yttria-Stabilized Zirconia

In this work, a set of density-functional tight-binding (DFTB) parameters for the Zr-Zr, Zr-O, Y-Y, Y-O, and Zr-Y interactions was developed for bulk and surface simulations of ZrO₂ (zirconia), Y₂O₃ (yttria), and yttria-stabilized zirconia (YSZ) materials. The parameterization lays the ground work for realistic simulations of zirconia-, yttria-, and YSZ-based electrolytes in solid oxide fuel cells and YSZ-based catalysts on long timescales and relevant size scales. The parameterization was validated for the zirconia and yttria polymorphs observed under standard conditions based on density functional theory calculations and experimental data. Additionally, we performed DFTB-based molecular dynamics (MD) simulations to compute structural and vibrational properties of these materials. The results show that the parameters can give a qualitatively correct phase ordering of zirconia, where the tetragonal phase is more stable than the cubic phase at a lower temperature. The lattice parameters are only slightly overestimated by 0.05-0.1 Å (2% error), still within the typical accuracy of first-principles methods. Additionally, the MD results confirm that zirconia and yttria phases are stable against transformations under standard conditions. The parameterization also predicts that vibrational spectra are within the range of 100-1000 cm^-1 for zirconia and 100-800 cm^-1 for yttria, which is in good agreement with predictions both from full quantum mechanics and a recently developed classical force field. To further demonstrate the advantage of the developed DFTB parameters in terms of computational resources, we conducted DFTB/MD simulations of the YSZ4 and YS12 models containing approximately 750 atoms.

Modelling the bulk properties of ambient pressure polymorphs of zirconia

We report a detailed survey of the calculated bulk properties of zirconia using GGA and meta-GGA functionals (PBE, PBEsol, RPBE, and TPSS), dispersion (Grimme's D2 and D3 approach), and on-site Coulomb repulsion correction (U = 2-8 eV). Structural, elastic, mechanical, and dielectric properties, as well as energetics, electronic structure, and phonon dispersion curves were computed and compared to previous investigations to identify the best DFT approach for a consistent in silico description of zirconia polymorphs. In general, inclusion of dispersion corrections led to only small changes in the calculated properties, whereas DFT+U (U = 2 or 4 eV) reduced the deviations of calculated properties from the experimental results, although deterioration of the structure and relative stabilities may be observed in some cases. Standard PBEsol, RPBE+U, and PBE+U were the best methodologies for a simultaneous description of the three polymorphs of ZrO2. RPBE+U, however, was the only functional to conserve the distinct structures and stabilities of c-, t-, and m-ZrO2 when U = 4 eV was used, resulting in the best in silico replication of the band gaps of ZrO2, whilst outperforming the other methodologies in the description of elastic, mechanical, and dielectric properties of this material. Overall, these results provide insight into the most appropriate DFT methodology for in silico investigations of ZrO2, and show that simultaneous description of all three ambient pressure zirconia polymorphs by DFT techniques with acceptable levels of accuracy can be achieved only when the correct choice of methodology is applied.

Theoretical insights into inorganic-organic intercalation products of the layered perovskite HLaNb2O7: perspectives for hybrid proton conductors

The modification of metal oxide surfaces with organic moieties has been widely studied as a method of preparing organic-inorganic hybrid materials for various applications. Among the inorganic oxides, ion-exchangeable layered perovskites are particularly interesting, because of their appealing electronic and reactive properties. In particular, their protonated interlayer surface can be easily functionalized with organic groups allowing the production of stable hybrid materials. As a further step in the design of new inorganic-organic hybrid proton conductors, a combined experimental and theoretical study of two intercalated compounds (propanol and imidazole) in HLaNb2O7 is presented here. A generally very good agreement with the available experimental data is found in reproducing both structural features and 13C-NMR chemical shifts, and marked differences between the two considered intercalated compounds are evidenced, with possible important outcomes for proton conduction. Notably, the free imidazole molecules are easily protonated by the acidic protons present in the interlayer spacing, thus inhibiting an efficient charge transport mechanism. In order to overcome this problem, a model system has been considered, where the imidazoles are bound to the end of a butyl chain, the whole being intercalated between two perovskite layers. The obtained theoretical data suggest that, in such a system, proton transfer between two adjacent imidazoles is a barrierless process. These results could then open new perspectives for such hybrid proton conductors.

Computational Design of Novel Hydrogen-Doped, Oxygen-Deficient Monoclinic Zirconia with Excellent Optical Absorption and Electronic Properties

Monoclinic ZrO2 has recently emerged as a new highly efficient material for the photovoltaic and photocatalytic applications. Herein, first-principles calculations were carried out to understand how Hydrogen doping can affect the electronic structure and optical properties of the material. The effects of Hydrogen interstitial and substitutional doping at different sites and concentrations in m-ZrO2 were examined by an extensive model study to predict the best structure with the optimal properties for use in solar energy conversion devices. Hydrogen interstitials (Hi) in pristine m-ZrO2 were found to lower the formation energy but without useful effects on the electronic or optical properties. Hydrogen mono- and co-occupying oxygen vacancy (Ov) were also investigated. At low concentration of Hydrogen mono-occupying oxygen vacancy (HOv), Hydrogen atoms introduced shallow states below the conduction band minimum (CBM) and increase the dielectric constant, which could be very useful for gate dielectric application. The number and position of such defect states strongly depend on the doping sites and concentration. At high oxygen vacancy concentration, the modeled HOv-Ov structure shows the formation of shallow and localized states that are only 1.1 eV below the CBM with significantly high dielectric constant and extended optical absorption to the infrared region. This strong absorption with the high permittivity and low exciton binding energies make the material an ideal candidate for use in solar energy harvesting devices. Finally, the band edge positions of pristine and doped structures with respect to the redox potentials of water splitting indicated that Hydrogen occupying oxygen vacancies can increase the photocatalytic activity of the material for hydrogen generation due the extremely improved optical absorption and the band gap states.

Ab initio insights into the structural, energetic, electronic, and stability properties of mixed CenZr15−nO30 nanoclusters

In this work, we report a DFT investigation combined with Spearman correlation analysis of the energetic, structural and electronic properties of mixed Ce_nZr_15−nO₃₀ nanoclusters as a function of the composition (n = 0, 1,…,14, 15).

Predicting Monolayer Oxide Stability over Low-Index Surfaces of TiO2 Polymorphs Using ab Initio Thermodynamics

Monolayer metal oxide coatings on metal oxide supports have the possibility of tuning the surface chemical properties of the coated systems. However, the (meta)stability of these structures makes experimental discovery challenging. A computational approach can help to determine properties that make a coating/substrate system stable and evaluate the stability of a variety of combinations. Herein, we use density functional theory (DFT) to study the stability of monolayer transitional metal oxides over different facets of anatase, brookite, and rutile phase of TiO₂. We find that coatings that have a stable polymorph matching that of the support, as well as substrates with higher surface energies, are more likely to form monolayer-coated systems. DFT calculations recommend a number of coating/TiO₂ surface facet combinations that may be stable. Despite these predictive observations, we did not find a significant correlation between monolayer stability and a single atomic, surface, or structural property of the coating/support metal/metal oxide and coating oxide monolayer stability. More complex predictive relationships need future study.

Revealing the properties of the cubic ZrO2 (111) surface by periodic DFT calculations: reducibility and stabilization through doping with aliovalent Y2O3

The reducibility of the clean cubic ZrO₂ (111) surface, as well as its stabilization through doping, have been investigated by hybrid DFT calculations within a periodic approach and localized basis sets.

The growth of Nin clusters and their interaction with cubic, monoclinic, and tetragonal ZrO2 surfaces–a theoretical and experimental study

The interaction and growth of Ni clusters on three phases of ZrO₂ surfaces were investigated through theoretical and experimental methods.