Concentration of vacancies at metal-oxide surfaces: case study of MgO(100)

Fabrication, structural, and enhanced mechanical behavior of MgO substituted PMMA composites for dental applications

The most common denture material used for dentistry is poly-methyl-methacrylate (PMMA). Usually, the polymeric PMMA material has numerous biological, mechanical and cost-effective shortcomings. Hence, to resolve such types of drawbacks, attempts have been made to investigate fillers of the PMMA like alumina (Al2O3), silica (SiO2), zirconia (ZrO2) etc. For the enhancement of the PMMA properties a suitable additive is required for its orthopedic applications. Herein, the main motive of this study was to synthesize a magnesium oxide (MgO) reinforced polymer-based hybrid nano-composites by using heat cure method with superior optical, biological and mechanical characteristics. For the structural and vibrational studies of the composites, XRD and FT-IR were carried out. Herein, the percentage of crystallinity for all the fabricated composites were also calculated and found to be 14.79-30.31. Various physical and optical parameters such as density, band gap, Urbach energy, cutoff energy, cutoff wavelength, steepness parameter, electron-phonon interaction, refractive index, and optical dielectric constant were also studied and their values are found to be in the range of 1.21-1.394 g/cm3, 5.44-5.48 eV, 0.167-0.027 eV, 5.68 eV, 218 nm, 0.156-0.962, 4.273-0.693, 1.937-1.932, and 3.752-3.731 respectively. To evaluate the mechanical properties like compressive strength, flexural strength, and fracture toughness of the composites a Universal Testing Machine (UTM) was used and their values were 60.3 and 101 MPa, 78 and 40.3 MPa, 5.85 and 9.8 MPa-m1/2 respectively. Tribological tests of the composites were also carried out. In order to check the toxicity, MTT assay was also carried out for the PM0 and PM15 [(x)MgO + (100 - x) (C5O2H8)n] (x = 0 and 15) composites. This study provides a comprehensive insight into the structural, physical, optical, and biological features of the fabricated PMMA-MgO composites, highlighting the potential of the PM15 composite with its enhanced density, mechanical strength, and excellent biocompatibility for denture applications.

CO adsorption on MgO thin-films: formation and interaction of surface charged defects

Two-dimensional (2D) materials formed by thin-films of metal oxides that grow on metal supports are commonly used in heterogeneous catalysis and multilayer electronic devices. Despite extensive research on these systems, the effects of charged defects at supported oxides on surface processes are still not clear. In this work, we perform spin-polarized density-functional theory (DFT) calculations to investigate formation and interaction of charged magnesium and oxygen vacancies, and Al dopants on MgO(001)/Ag(001) surface. The results show a sizable interface compressive effect that decreases the metal work function as electrons are added on the MgO surface with a magnesium vacancy. This surface displays a larger formation energy in a water environment (O-rich condition) even with additional Al-doping. Under these conditions, we found that a polar molecule such as CO is more strongly adsorbed on the low-coordination oxygen sites due to a larger contribution of the channeled electronic transport with the silver interface regardless of the surface charge. Therefore, these findings elucidate how surface intrinsic vacancies can influence or contribute to charge transfer, which allows one to explore more specific reactions at different surface topologies for more efficient catalysts for CO2 conversion.

Electronic and magnetic properties of charged point defects in monolayer CrI3

The two-dimensional magnetic material CrI₃ has gained considerable attention owing to its promising applications in photoelectric and spin-related devices. Recently, various structural defects in CrI₃ have been identified; however, the charge states of these defects have been mainly ignored. Here, we report on an investigation of the charged defects in monolayer CrI₃, focused on the electronic and magnetic properties of the five most stable point defects using first-principles calculations. For positively charged I vacancies and negatively charged Cr vacancies, a blue- and red-shift of defect states near the Fermi level can be observed because of the atom relaxation. Our results also indicate that, among the five defects, the Cr interstitial defect has the smallest ionization energy of 0.34 eV, which makes its ionization easiest. Furthermore, a 0.2 μ_B enhancement of the magnetic moment on the nearest Cr atom can be found for the I vacancy and Cr interstitial defect. The investigation contributes to the atomic-scale comparison and understanding of the charged defects of monolayer CrI₃.

The application of QM/MM simulations in heterogeneous catalysis

The QM/MM simulation method is provenly efficient for the simulation of biological systems, where an interplay of extensive environment and delicate local interactions drives a process of interest through a funnel on a complex energy landscape. Recent advances in quantum chemistry and force-field methods present opportunities for the adoption of QM/MM to simulate heterogeneous catalytic processes, and their related systems, where similar intricacies exist on the energy landscape. Herein, the fundamental theoretical considerations for performing QM/MM simulations, and the practical considerations for setting up QM/MM simulations of catalytic systems, are introduced; then, areas of heterogeneous catalysis are explored where QM/MM methods have been most fruitfully applied. The discussion includes simulations performed for adsorption processes in solvent at metallic interfaces, reaction mechanisms within zeolitic systems, nanoparticles, and defect chemistry within ionic solids. We conclude with a perspective on the current state of the field and areas where future opportunities for development and application exist.

Toward a Consistent Prediction of Defect Chemistry in CeO2

Polarizable shell-model potentials are widely used for atomic-scale modeling of charged defects in solids using the Mott-Littleton approach and hybrid Quantum Mechanical/Molecular Mechanical (QM/MM) embedded-cluster techniques. However, at the pure MM level of theory, the calculated defect energetics may not satisfy the requirement of quantitative predictions and are limited to only certain charged states. Here, we proposed a novel interatomic potential development scheme that unifies the predictions of all relevant charged defects in CeO₂ based on the Mott-Littleton approach and QM/MM electronic-structure calculations. The predicted formation energies of oxygen vacancies accompanied by different excess electron localization patterns at the MM level of theory reach the accuracy of density functional theory (DFT) calculations using hybrid functionals. The new potential also accurately reproduces a wide range of physical properties of CeO₂, showing excellent agreement with experimental and other computational studies. These findings provide opportunities for accurate large-scale modeling of the partial reduction and nonstoichiometry in CeO₂, as well as a prototype for developing robust interatomic potentials for other defective crystals.

Counting Point Defects at Nanoparticle Surfaces by Electron Holography

Metal oxide nanoparticles exhibit outstanding catalytic properties, believed to be related to the presence of oxygen vacancies at the particle's surface. However, little quantitative information is known about concentrations of point defects inside and at surfaces of these nanoparticles, due to the challenges in achieving an atomically resolved experimental access. By employing off-axis electron holography, we demonstrate, using MgO nanoparticles as an example, a methodology that discriminates between mobile charge induced by electron beam irradiation and immobile charge associated with deep traps induced by point defects as well as distinguishes between bulk and surface point defects. Counting the immobile charge provides a quantification of the concentration of F²⁺ centers induced by oxygen vacancies at the MgO nanocube surfaces.

Favorable Role of the Metal-Support Perimeter Region in Electrochemical NH3 Synthesis: A Density Functional Theory Study on Ru/BaCeO3

The catalytic electrochemical synthesis of NH₃ on Ru/BaCeO₃ was investigated using density functional theory. The competition between NH₃ formation and the hydrogen evolution reaction (HER) is a key for a high NH₃ formation rate. Our calculations show that H adsorbs more strongly than N₂ at the Ru particle moiety, while the adsorption of N₂ is stronger than the H adsorption at the Ru/BaCeO₃ perimeter, a model for the triple-phase boundary that is proposed to be an active site by experimental studies. This indicates that, while the HER is more favorable at the Ru particle moiety, it should be suppressed at the Ru/BaCeO₃ perimeter. We also calculated the Gibbs free energy changes along the NH₃ formation and found that the N₂H formation, the NHNH₂ formation, and the NH₃ formation steps have a relatively large Gibbs energy change. Therefore, these are possible candidates for the potential-determining step. The calculated equilibrium potential (U = -0.70 V, vs RHE) is in reasonable agreement with experiments. We also evaluated the reaction energy (ΔE) and the activation barrier (E _a) of the N₂H formation at several sites. ΔE and E _a were high at the Ru particle moiety (ΔE = 1.18 eV and E _a = 1.38 eV) but became low (ΔE = 0.32 eV and E _a = 1.31 eV) at the Ru/BaCeO₃ perimeter. These provide the atomic-scale mechanism how the proton conduction in BaCeO₃ assists the electrochemical NH₃ synthesis.

Defect engineering of oxide surfaces: dream or reality?

In this brief perspective we analyze the present status of the field of defect engineering of oxide surfaces. In particular we discuss the tools and techniques available to generate, identify, quantify, and characterize point defects at oxide surfaces and the main areas where these centers play a role in practical applications.

Discrete modeling of ionic space charge zones in solids

The discrete model of space charge zones in solids reveals and remedies a variety of problems with the classic continuous Gouy-Chapman solution that occur for pronounced space charge potentials. Besides inherent problems of internal consistency, it is essentially the extremely steep profile close to the interface which makes this continuum approach questionable. Not only is quasi-1D discrete modeling a sensible approach for large space charge effects, it can also favorably be combined with the continuum description. A particularly useful application is the explicit implementation of crystallographic details and non-idealities close to the interface. This enables us to consider elastic, structural or saturation effects as well as permittivity variations in a simple but realistic way. We address details of the charge carrier profiles, but also overall properties such as space charge capacitance and space charge resistance. In the latter case the difference in the total charge (at identical concentration) is of importance, in the first case it is the inherent difference in the centroid of charge (at identical total charge) that is remarkable. The model is equally applicable for ionic charge carriers and small polarons.

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.

Artificial-intelligence-driven discovery of catalyst genes with application to CO2 activation on semiconductor oxides

Catalytic-materials design requires predictive modeling of the interaction between catalyst and reactants. This is challenging due to the complexity and diversity of structure-property relationships across the chemical space. Here, we report a strategy for a rational design of catalytic materials using the artificial intelligence approach (AI) subgroup discovery. We identify catalyst genes (features) that correlate with mechanisms that trigger, facilitate, or hinder the activation of carbon dioxide (CO2) towards a chemical conversion. The AI model is trained on first-principles data for a broad family of oxides. We demonstrate that surfaces of experimentally identified good catalysts consistently exhibit combinations of genes resulting in a strong elongation of a C-O bond. The same combinations of genes also minimize the OCO-angle, the previously proposed indicator of activation, albeit under the constraint that the Sabatier principle is satisfied. Based on these findings, we propose a set of new promising catalyst materials for CO2 conversion.

Excited States of Crystalline Point Defects with Multireference Density Matrix Embedding Theory

Accurate and affordable methods to characterize the electronic structure of solids are important for targeted materials design. Embedding-based methods provide an appealing balance in the trade-off between cost and accuracy─particularly when studying localized phenomena. Here, we use the density matrix embedding theory (DMET) algorithm to study the electronic excitations in solid-state defects with a restricted open-shell Hartree-Fock (ROHF) bath and multireference impurity solvers, specifically, complete active space self-consistent field (CASSCF) and n-electron valence state second-order perturbation theory (NEVPT2). We apply the method to investigate the electronic excitations in an oxygen vacancy (OV) on a MgO(100) surface and find absolute deviations within 0.05 eV between DMET using the CASSCF/NEVPT2 solver, denoted as CAS-DMET/NEVPT2-DMET, and the nonembedded CASSCF/NEVPT2 approach. Next, we establish the practicality of DMET by extending it to larger supercells for the OV defect and a neutral silicon vacancy in diamond where the use of nonembedded CASSCF/NEVPT2 is extremely expensive.

First-principles calculations of hybrid inorganic-organic interfaces: from state-of-the-art to best practice

The computational characterization of inorganic-organic hybrid interfaces is arguably one of the technically most challenging applications of density functional theory. Due to the fundamentally different electronic properties of the inorganic and the organic components of a hybrid interface, the proper choice of the electronic structure method, of the algorithms to solve these methods, and of the parameters that enter these algorithms is highly non-trivial. In fact, computational choices that work well for one of the components often perform poorly for the other. As a consequence, default settings for one materials class are typically inadequate for the hybrid system, which makes calculations employing such settings inefficient and sometimes even prone to erroneous results. To address this issue, we discuss how to choose appropriate atomistic representations for the system under investigation, we highlight the role of the exchange-correlation functional and the van der Waals correction employed in the calculation and we provide tips and tricks how to efficiently converge the self-consistent field cycle and to obtain accurate geometries. We particularly focus on potentially unexpected pitfalls and the errors they incur. As a summary, we provide a list of best practice rules for interface simulations that should especially serve as a useful starting point for less experienced users and newcomers to the field.

Self-Consistent Potential Correction for Charged Periodic Systems

Supercell models are often used to calculate the electronic structure of local deviations from the ideal periodicity in the bulk or on the surface of a crystal or in wires. When the defect or adsorbent is charged, a jellium counter charge is applied to maintain overall neutrality, but the interaction of the artificially repeated charges has to be corrected, both in the total energy and in the one-electron eigenvalues and eigenstates. This becomes paramount in slab or wire calculations, where the jellium counter charge may induce spurious states in the vacuum. We present here a self-consistent potential correction scheme and provide successful tests of it for bulk and slab calculations.

First-Principles Atomistic Thermodynamics and Configurational Entropy

In most applications, functional materials operate at finite temperatures and are in contact with a reservoir of atoms or molecules (gas, liquid, or solid). In order to understand the properties of materials at realistic conditions, statistical effects associated with configurational sampling and particle exchange at finite temperatures must consequently be taken into account. In this contribution, we discuss the main concepts behind equilibrium statistical mechanics. We demonstrate how these concepts can be used to predict the behavior of materials at realistic temperatures and pressures within the framework of atomistic thermodynamics. We also introduce and discuss methods for calculating phase diagrams of bulk materials and surfaces as well as point defect concentrations. In particular, we describe approaches for calculating the configurational density of states, which requires the evaluation of the energies of a large number of configurations. The cluster expansion method is therefore also discussed as a numerically efficient approach for evaluating these energies.

Iron and oxygen vacancies at the hematite surface: pristine case and with a chlorine adatom

Defect complexes play critical roles in the dynamics of water molecules in photoelectrochemical cell devices. For the specific case of hematite (α-Fe₂O₃), iron and oxygen vacancies are said to mediate the water splitting process through the localization of optically-derived charges. Using first-principles methods based on density-functional theory we show that both iron and oxygen vacancies can be observed at the surface. For an oxygen-rich environment, usually under wet conditions, the charged iron vacancies should be more frequent. As sea water would be an ideal electrolyte for this kind of device, we have also analyzed the effect of additional chlorine adsorption on this surface. While the chlorine adatom kills the charged oxygen vacancies, entering the void sites, it will not react with the iron vacancies, keeping them active during water splitting processes.

Final-State Simulations of Core-Level Binding Energies at Metal-Organic Hybrid Interfaces: Artifacts Caused by Spurious Collective Electrostatic Effects

Core-level energies are frequently calculated to explain the X-ray photoelectron spectra of metal-organic hybrid interfaces. The current paper describes how such simulations can be flawed when modeling interfaces between physisorbed organic molecules and metals. The problem occurs when applying periodic boundary conditions to correctly describe extended interfaces and simultaneously considering core hole excitations in the framework of a final-state approach to account for screening effects. Since the core hole is generated in every unit cell, an artificial dipole layer is formed. In this work, we study methane on an Al(100) surface as a deliberately chosen model system for hybrid interfaces to evaluate the impact of this computational artifact. We show that changing the supercell size leads to artificial shifts in the calculated core-level energies that can be well beyond 1 eV for small cells. The same applies to atoms at comparably large distances from the substrate, encountered, for example, in extended, upright-standing adsorbate molecules. We also argue that the calculated work function change due to a core-level excitation can serve as an indication for the occurrence of such an artifact and discuss possible remedies for the problem.

Prediction of Site Preference of Implanted Transition Metal Dopants in Rock-salt Oxides

Transition metals (TMs) implanted in oxides with rock-salt crystal structures (for example MgO and BaO) are assumed to substitute cations (Mg in case of MgO) from the lattice sites. We show that not all implanted TMs substitute cations but can be stable in interstitial sites as well. Stability of TM (Sc-Zn) dopants in various charge states in MgO and BaO has been investigated in the framework of density functional theory. We propose an effective way to calculate stability of implanted metals that let us predict site preference (interstitial or substitution) of the dopant in the host. We find that two factors govern the preference for an interstitial site: (i) relative ionic radius and (ii) relative oxygen affinity of cation and the TM dopants. If the radius of the cation is much larger than TM dopant, as in BaO, TM atoms always sit at interstitial sites. On the other hand, if the radius of the cation is comparable to that of the dopant TM, as in case of MgO, the transition of the preferred defect site, from substituting lattice Mg atom (Sc to Mn) to occupying interstitial site (Fe to Zn) is observed. This transition can be attributed to the change in the oxygen affinity of the TM atoms from Sc to Zn. Our results also explain experiments on Ni and Fe atoms implanted in MgO. TM dopants at interstitial sites could show substantially different and new properties from substitutionally doped stable compounds.

Interface enhanced superconductivity in monolayer FeSe films on MgO(001): charge transfer with atomic substitution

Interface enhanced superconductivity over 50 K has been discovered in monolayer FeSe films grown on several TiO₂-terminated oxide substrates. Whether such phenomenon exists in other oxide substrates remains an extremely interesting topic. Here we report enhanced superconductivity with an onset transition temperature of 18 K in monolayer FeSe on MgO(001) substrate by transport measurement. Scanning transmission electron microscopy investigation on the interface structure indicates that FeSe films grow epitaxially on MgO(001) and that overlayer Fe atoms diffuse into the top two layers of MgO and substitute Mg atoms. Our density functional theory calculations reveal that this substitution promotes the charge transfer from the MgO substrate to the FeSe films, an essential process that also occurs in monolayer FeSe on TiO₂-terminated oxides and contributes to the enhanced superconductivity therein. Our finding suggests that superconductivity enhancement in monolayer FeSe films on oxides substrates is rather general as long as charge transfer is allowed at the interface, thus pointing out an explicit direction for searching for new high temperature superconductivity by interface engineering.

The Impact of Atmosphere on the Local Luminescence Properties of Metal Halide Perovskite Grains

Metal halide perovskites are exceptional candidates for inexpensive yet high-performing optoelectronic devices. Nevertheless, polycrystalline perovskite films are still limited by nonradiative losses due to charge carrier trap states that can be affected by illumination. Here, in situ microphotoluminescence measurements are used to elucidate the impact of light-soaking individual methylammonium lead iodide grains in high-quality polycrystalline films while immersing them with different atmospheric environments. It is shown that emission from each grain depends sensitively on both the environment and the nature of the specific grain, i.e., whether it shows good (bright grain) or poor (dark grain) luminescence properties. It is found that the dark grains show substantial rises in emission, while the bright grain emission is steady when illuminated in the presence of oxygen and/or water molecules. The results are explained using density functional theory calculations, which reveal strong adsorption energies of the molecules to the perovskite surfaces. It is also found that oxygen molecules bind particularly strongly to surface iodide vacancies which, in the presence of photoexcited electrons, lead to efficient passivation of the carrier trap states that arise from these vacancies. The work reveals a unique insight into the nature of nonradiative decay and the impact of atmospheric passivation on the microscale properties of perovskite films.

Computational methods for 2D materials: discovery, property characterization, and application design

The discovery of two-dimensional (2D) materials comes at a time when computational methods are mature and can predict novel 2D materials, characterize their properties, and guide the design of 2D materials for applications. This article reviews the recent progress in computational approaches for 2D materials research. We discuss the computational techniques and provide an overview of the ongoing research in the field. We begin with an overview of known 2D materials, common computational methods, and available cyber infrastructures. We then move onto the discovery of novel 2D materials, discussing the stability criteria for 2D materials, computational methods for structure prediction, and interactions of monolayers with electrochemical and gaseous environments. Next, we describe the computational characterization of the 2D materials' electronic, optical, magnetic, and superconducting properties and the response of the properties under applied mechanical strain and electrical fields. From there, we move on to discuss the structure and properties of defects in 2D materials, and describe methods for 2D materials device simulations. We conclude by providing an outlook on the needs and challenges for future developments in the field of computational research for 2D materials.

Understanding band gaps of solids in generalized Kohn-Sham theory

The fundamental energy gap of a periodic solid distinguishes insulators from metals and characterizes low-energy single-electron excitations. However, the gap in the band structure of the exact multiplicative Kohn-Sham (KS) potential substantially underestimates the fundamental gap, a major limitation of KS density-functional theory. Here, we give a simple proof of a theorem: In generalized KS theory (GKS), the band gap of an extended system equals the fundamental gap for the approximate functional if the GKS potential operator is continuous and the density change is delocalized when an electron or hole is added. Our theorem explains how GKS band gaps from metageneralized gradient approximations (meta-GGAs) and hybrid functionals can be more realistic than those from GGAs or even from the exact KS potential. The theorem also follows from earlier work. The band edges in the GKS one-electron spectrum are also related to measurable energies. A linear chain of hydrogen molecules, solid aluminum arsenide, and solid argon provide numerical illustrations.

Trends in Adhesion Energies of Metal Nanoparticles on Oxide Surfaces: Understanding Support Effects in Catalysis and Nanotechnology

Nanoparticles on surfaces are ubiquitous in nanotechnologies, especially in catalysis, where metal nanoparticles anchored to oxide supports are widely used to produce and use fuels and chemicals, and in pollution abatement. We show that for hemispherical metal particles of the same diameter, D, the chemical potentials of the metal atoms in the particles (μ_M) differ between two supports by approximately -2(E_adh,A - E_adh,B)V_m/D, where E_ad,i is the adhesion energy between the metal and support i, and V_m is the molar volume of the bulk metal. This is consistent with calorimetric measurements of metal vapor adsorption energies onto clean oxide surfaces where the metal grows as 3D particles, which proved that μ_M increases with decreasing particle size below 6 nm and, for a given size, decreases with E_adh. Since catalytic activity and sintering rates correlate with metal chemical potential, it is thus crucial to understand what properties of catalyst materials control metal/oxide adhesion energies. Trends in how E_adh varies with the metal and the support oxide are presented. For a given oxide, E_adh increases linearly from metal to metal with increasing heat of formation of the most stable oxide of the metal (per mole metal), or metal oxophilicity, suggesting that metal-oxygen bonds dominate interfacial bonding. For the two different stoichiometric oxide surfaces that have been studied on multiple metals (MgO(100) and CeO₂(111), the slopes of these lines are the same, but their offset is large (∼2 J/m²). Adhesion energies increase as MgO(100) ≈ TiO₂(110) < α-Al₂O₃(0001) < CeO₂(111) ≈ Fe₃O₄(111).

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

Coupled-cluster theory with single, double, and perturbative triple excitations (CCSD(T)) is widely considered to be the "gold standard" of ab initio quantum chemistry. Using the domain-based pair natural orbital local correlation concept (DLPNO-CCSD(T)), these calculations can be performed on systems with hundreds of atoms at an accuracy of ∼99.9% of the canonical CCSD(T) method. This allows for ab initio calculations providing reference adsorption energetics at solid surfaces with an accuracy approaching 1 kcal/mol. This is an invaluable asset, not least for the assessment of density functional theory (DFT) as the prevalent approach for large-scale production calculations in energy or catalysis applications. Here we use DLPNO-CCSD(T) with embedded cluster models to compute entire adsorbate potential energy surfaces for the binding of a set of prototypical closed-shell molecules (H₂O, NH₃, CH₄, CH₃OH, CO₂) to the rutile TiO₂(110) surface. The DLPNO-CCSD(T) calculations show excellent agreement with available experimental data, even for the "infamous" challenge of correctly predicting the CO₂ adsorption geometry. The numerical efficiency of the approach is within 1 order of magnitude of hybrid-level DFT calculations, hence blurring the borders between reference and production technique.

Origin of Symmetric Dimer Images of Si(001) Observed by Low-Temperature Scanning Tunneling Microscopy

It has been a long-standing puzzle why buckled dimers of the Si(001) surface appeared symmetric below ~20 K in scanning tunneling microscopy (STM) experiments. Although such symmetric dimer images were concluded to be due to an artifact induced by STM measurements, its underlying mechanism is still veiled. Here, we demonstrate, based on a first-principles density-functional theory calculation, that the symmetric dimer images are originated from the flip-flop motion of buckled dimers, driven by quantum tunneling (QT). It is revealed that at low temperature the tunneling-induced surface charging with holes reduces the energy barrier for the flipping of buckled dimers, thereby giving rise to a sizable QT-driven frequency of the flip-flop motion. However, such a QT phenomenon becomes marginal in the tunneling-induced surface charging with electrons. Our findings provide an explanation for low-temperature STM data that exhibits apparent symmetric (buckled) dimer structure in the filled-state (empty-state) images.

Modelling the chemistry of Mn-doped MgO for bulk and (100) surfaces

We have investigated the energetic properties of Mn-doped MgO bulk and (100) surfaces using a QM/MM embedding computational method, calculating the formation energy for doped systems, as well as for surface defects, and the subsequent effect on chemical reactivity.

Equilibrium Vacancy Concentration Driven by Undetectable Impurities

By means of ab initio calculations combined to statistical mechanics, we provide new evidence that an experimentally undetectable tiny amount of impurities can be responsible for drastic changes in vacancy concentrations ([V]), inducing large deviations from an Arrhenius law even at low temperature. It is the case of O and N in α-Fe. The present finding is fully compatible with existing experiments, and changes the previous common vision that C has the dominant effect. This study provides a route for bridging the longstanding theoretical-experimental gap on the prediction of [V] in metals.

Determination of formation and ionization energies of charged defects in two-dimensional materials

We present a simple and efficient approach to evaluate the formation energy and, in particular, the ionization energy (IE) of charged defects in two-dimensional (2D) systems using the supercell approximation. So far, first-principles results for such systems can scatter widely due to the divergence of the Coulomb energy with vacuum dimension, denoted here as L_{z}. Numerous attempts have been made in the past to fix the problem under various approximations. Here, we show that the problem can be resolved without any such assumption, and a converged IE can be obtained by an extrapolation of the asymptotic IE expression at large L_{z} (with a fixed lateral area S) back to the value at L_{z}=0. Application to defects in monolayer boron nitride reveal that defects in 2D systems can be unexpectedly deep, much deeper than the bulk.

Heterogeneous catalysis

A heterogeneous catalyst is a functional material that continually creates active sites with its reactants under reaction conditions. These sites change the rates of chemical reactions of the reactants localized on them without changing the thermodynamic equilibrium between the materials.

Adsorption and reactions of O2 on anatase TiO2

CONSPECTUS: The interaction of molecular oxygen with titanium dioxide (TiO2) surfaces plays a key role in many technologically important processes such as catalytic oxidation reactions, chemical sensing, and photocatalysis. While O2 interacts weakly with fully oxidized TiO2, excess electrons are often present in TiO2 samples. These excess electrons originate from intrinsic reducing defects (oxygen vacancies and titanium interstitials), doping, or photoexcitation and form polaronic Ti(3+) states in the band gap near the bottom of the conduction band. Oxygen adsorption involves the transfer of one or more of these excess electrons to an O2 molecule at the TiO2 surface. This results in an adsorbed superoxo (O2(-)) or peroxo (O2(2-)) species or in molecular dissociation and formation of two oxygen adatoms (2 × O(2-)). Oxygen adsorption is also the first step toward oxygen incorporation, a fundamental reaction that strongly affects the chemical properties and charge-carrier densities; for instance, it can transform the material from an n-type semiconductor to a poor electronic conductor. In this Account, we present an overview of recent theoretical work on O2 adsorption and reactions on the reduced anatase (101) surface. Anatase is the TiO2 polymorph that is generally considered most active in photocatalysis. Experiments on anatase powders have shown that the properties of photoexcited electrons are similar to those of excess electrons from reducing defects, and therefore, oxygen on reduced anatase is also a model system for studying the role of O2 in photocatalysis. Experimentally, the characteristic Ti(3+) defect states disappear after adsorption of molecular oxygen, which indicates that the excess electrons are indeed trapped by O2. Moreover, superoxide surface species associated with two different cation surface sites, possibly a regular cation site and a cation close to an anion vacancy, were identified by electron paramagnetic resonance spectroscopy. On the theoretical side, however, density functional theory studies have consistently found that it is energetically more favorable for O2 to adsorb in the peroxo form rather than the superoxo form. As a result, obtaining a detailed understanding of the nature of the observed superoxide species has proven difficult for many years. On reduced anatase (101), both oxygen vacancies and Ti interstitials have been shown to reside exclusively in the susbsurface. We discuss how reaction of O2 with a subsurface O vacancy heals the vacancy while leading to the formation of a surface bridging dimer defect. Similarly, the interaction of O2 with a Ti interstitial causes migration of this defect to the surface and the formation of a surface TiO2 cluster. Finally, we analyze the peroxo and superoxo states of the adsorbed molecule. On the basis of periodic hybrid functional calculations of interfacial electron transfer between reduced anatase and O2, we show that the peroxide form, while energetically more stable, is kinetically less favorable than the superoxide form. The existence of a kinetic barrier between the superoxo and peroxo states is essential for explaining a variety of experimental observations.

Defect-driven interfacial electronic structures at an organic/metal-oxide semiconductor heterojunction

The electronic structure of the hybrid interface between ZnO and the prototypical organic semiconductor PTCDI is investigated via a combination of ultraviolet and X-ray photoelectron spectroscopy (UPS/XPS) and density functional theory (DFT) calculations. The interfacial electronic interactions lead to a large interface dipole due to substantial charge transfer from ZnO to 3,4,9,10-perylenetetracarboxylicdiimide (PTCDI), which can be properly described only when accounting for surface defects that confer ZnO its n-type properties.

Toward a Quantitative Understanding of the Electric Field in Thermal Metal Oxidation and a Self-Consistent Wagner Theory

The electric field in the growing oxide film is important to the kinetics and mechanism of metal oxidation. However, understanding of the essential characteristics of the electric field during oxidation remains insufficient. A special-case analytical model is presented that provides a unified understanding for the electric field from the viewpoints of kinetics and thermodynamics. More general cases are studied by computer simulations that show similar characteristics in the electric field. In particular, simulations indicate that in many situations, the electrostatic potential drop across the bulk oxide is limited to ∼kBT/e, which means that the total electrostatic potential drop across the oxide film, if on the order of 1 V by rough estimation, should have contributions mostly from the electrified interfaces. Finally, regarding the Gibbs-Duhem relation, the commonly used isobaric assumption for the diffusing species is refuted. The results contained herein also provide a self-consistent understanding of Wagner's oxidation theory.

Electronic and chemical properties of a surface-terminated screw dislocation in MgO

Dislocations represent an important and ubiquitous class of topological defect found at the surfaces of metal oxide materials. They are thought to influence processes as diverse as crystal growth, corrosion, charge trapping, luminescence, molecular adsorption, and catalytic activity; however, their electronic and chemical properties remain poorly understood. Here, through a detailed first-principles investigation into the properties of a surface-terminated screw dislocation in MgO we provide atomistic insight into these issues. We show that surface dislocations can exhibit intriguing electron trapping properties which are important for understanding the chemical and electronic characteristics of oxide surfaces. The results presented in this article taken together with recent experimental reports show that surface dislocations can be equally as important as more commonly considered surface defects, such as steps, kinks, and vacancies, but are now just beginning to be understood.

Space-charge transfer in hybrid inorganic-organic systems

We discuss density functional theory calculations of hybrid inorganic-organic systems that explicitly include the global effects of doping (i.e., position of the Fermi level) and the formation of a space-charge layer. For the example of tetrafluoro-tetracyanoquinodimethane on the ZnO(0001[over ¯]) surface we show that the adsorption energy and electron transfer depend strongly on the ZnO doping. The associated work function changes are large, for which the formation of space-charge layers is the main driving force. The prominent doping effects are expected to be quite general for charge-transfer interfaces in hybrid inorganic-organic systems and important for device design.

Stability and metastability of clusters in a reactive atmosphere: theoretical evidence for unexpected stoichiometries of MgMOx

By applying a genetic algorithm and ab initio atomistic thermodynamics, we identify the stable and metastable compositions and structures of MgMOx clusters at realistic temperatures and oxygen pressures. We find that small clusters (M≲5) are in thermodynamic equilibrium when x>M. The nonstoichiometric clusters exhibit peculiar magnetic behavior, suggesting the possibility of tuning magnetic properties by changing environmental pressure and temperature conditions. Furthermore, we show that density-functional theory with a hybrid exchange-correlation functional is needed for predicting accurate phase diagrams of metal-oxide clusters. Neither a (sophisticated) force field nor density-functional theory with (semi)local exchange-correlation functionals is sufficient for even a qualitative prediction.