Paramagnetic States in Oxygen-Doped Boron Nitride Extend Light Harvesting and Photochemistry to the Deep Visible Region.

A family of boron nitride (BN)-based photocatalysts for solar fuel syntheses have recently emerged. Studies have shown that oxygen doping, leading to boron oxynitride (BNO), can extend light absorption to the visible range. However, the fundamental question surrounding the origin of enhanced light harvesting and the role of specific chemical states of oxygen in BNO photochemistry remains unanswered. Here, using an integrated experimental and first-principles-based computational approach, we demonstrate that paramagnetic isolated OB3 states are paramount to inducing prominent red-shifted light absorption. Conversely, we highlight the diamagnetic nature of O-B-O states, which are shown to cause undesired larger band gaps and impaired photochemistry. This study elucidates the importance of paramagnetism in BNO semiconductors and provides fundamental insight into its photophysics. The work herein paves the way for tailoring of its optoelectronic and photochemical properties for solar fuel synthesis.

Role of electron localisation in H adsorption and hydride formation in the Mg basal plane under aqueous corrosion: a first-principles study.

Understanding hydrogen-metal interactions is important in various fields of surface science, including the aqueous corrosion of metals. The interaction between atomic H and a Mg surface is a key process for the formation of sub-surface Mg hydride, which may play an important role in Mg aqueous corrosion. In the present work, we performed first-principles Density Functional Theory (DFT) calculations to study the mechanisms for hydrogen adsorption and crystalline Mg hydride formation under aqueous conditions. The Electron Localisation Function (ELF) is found to be a promising indicator for predicting stable H adsorption in the Mg surface. It is found that H adsorption and hydride layer formation is dominated by high ELF adsorption sites. Our calculations suggest that the on-surface adsorption of atomic H, OH radicals and atomic O could enhance the electron localisation at specific sites in the sub-surface region, thus forming effective H traps locally. This is predicted to result in the formation of a thermodynamically stable sub-surface hydride layer, which is a potential precursor of the crucial hydride corrosion product of magnesium.

Mind the Interface Gap: Exposing Hidden Interface Defects at the Epitaxial Heterostructure between CuO and Cu2O.

Well designed and optimized epitaxial heterostructures lie at the foundation of materials development for photovoltaic, photocatalytic, and photoelectrochemistry applications. Heterostructure materials offer tunable control over charge separation and transport at the same time preventing recombination of photogenerated excitations at the interface. Thus, it is of paramount importance that a detailed understanding is developed as the basis for further optimization strategies and design. Oxides of copper are nontoxic, low cost, abundant materials with a straightforward and stable manufacturing process. However, in individual applications, they suffer from inefficient charge transport of photogenerated carriers. Hence, in this work, we investigate the role of the interface between epitaxially aligned CuO and Cu2O to explore the potential benefits of such an architecture for more efficient electron and hole transfer. The CuO/Cu2O heterojunction nature, stability, bonding mechanism, interface dipole, electronic structure, and band bending were rationalized using hybrid density functional theory calculations. New electronic states are identified at the interface itself, which are originating neither from lattice mismatch nor strained Cu-O bonds. They form as a result of a change in coordination environment of CuO surface Cu2+ cations and an electron transfer across the interface Cu1+-O bond. The first process creates occupied defect-like electronic states above the valence band, while the second leaves hole states below the conduction band. These are constitutional to the interface and are highly likely to contribute to recombination effects competing with the improved charged separation from the suitable band bending and alignment and thus would limit the expected output photocurrent and photovoltage. Finally, a favorable effect of interstitial oxygen defects has been shown to allow for band gap tunability at the interface but only to the point of the integral geometrical contact limit of the heterostructure itself.

Strain-Tuneable Magnetism and Spintronics of Distorted Monovacancies in Graphene.

The electronic and spintronic properties of the monovacancies in freestanding and isotopically compressed graphene are investigated using hybrid exchange density functional perturbation theory. When the effects of electronic self-interaction are taken into account, an integer magnetic moment of 2 µB is identified for a Jahn-Teller reconstructed V1(5-9) monovacancy in freestanding graphene. For graphene with stable ripples induced by a compressive strain of 5%, a bond reconstruction produces a V1(55-66) structure for the monovacancy, which is localized at the saddle points of the ripple. The sizeable local distortion induced by reconstruction modifies both the geometric and electronic properties of rippled graphene and quenches the magnetic moment of the vacancy due to the sp3 hybridization of the central atom. The nonmagnetic V1(55-66) structure is found to be stable on rippled structures, with the formation energy ∼2.3 eV lower than that of the metastable distorted V1(5-9) structures localized at sites other than the saddle points. The electronic ground state of distorted V1(5-9) corresponds to a wide range of fractional magnetic moments (0.50-1.25 µB). The computed relative stabilities and the electronic and magnetic properties of the V1(5-9) structures are found to be closely related to their local distortions. This analysis of the fundamental properties of defective graphene under compression suggests a number of strategies for generating regular defect patterns with tuneable magnetic and electronic properties and may, therefore, be used as a novel technique to achieve more precise control of graphene electronic structure for various application scenarios such as transistors, strain sensors, and directed chemisorption.

A transferable prediction model of molecular adsorption on metals based on adsorbate and substrate properties.

Surface adsorption is one of the fundamental processes in numerous fields, including catalysis, the environment, energy and medicine. The development of an adsorption model which provides an effective prediction of binding energy in minutes has been a long term goal in surface and interface science. The solution has been elusive as identifying the intrinsic determinants of the adsorption energy for various compositions, structures and environments is non-trivial. We introduce a new and flexible model for predicting adsorption energies to metal substrates. The model is based on easily computed, intrinsic properties of the substrate and adsorbate, which are the same for all the considered systems. It is parameterised using machine learning based on first-principles calculations of probe molecules (e.g., H2O, CO2, O2, N2) adsorbed to a range of pure metal substrates. The model predicts the computed dissociative adsorption energy to metal surfaces with a correlation coefficient of 0.93 and a mean absolute error of 0.77 eV for the large database of molecular adsorption energies provided by Catalysis-Hub.org which have a range of 15 eV. As the model is based on pre-computed quantities it provides near-instantaneous estimates of adsorption energies and it is sufficiently accurate to eliminate around 90% of candidates in screening study of new adsorbates. The model, therefore, significantly enhances current efforts to identify new molecular coatings in many applied research fields.

Corrosion inhibition in acidic environments: key interfacial insights with photoelectron spectroscopy.

In many engineering scenarios, surface-active organic species are added to acidic solutions to inhibit the corrosion of metallic components. Given suitable selection, such corrosion inhibitors are highly effective, preventing significant degradation even in highly aggressive environments. Nevertheless, there are still considerable gaps in fundamental knowledge of corrosion inhibitor functionality, severely restricting rational development. Here, we demonstrate the capability of X-ray photoelectron spectroscopy (XPS), supported by ab initio modelling, for revealing key details of inhibited substrates. Attention is focussed on the corrosion inhibition of carbon steel through the addition of an exemplar imidazoline-based corrosion inhibitor (OMID) to aqueous solutions of both HCl and H2SO4. Most notably, it is demonstrated that interfacial chemistry varies with the identity of the acid. High resolution Fe 2p, O 1s, N 1s, and Cl 2p XPS spectra, acquired from well-inhibited carbon steel in 1 M HCl, show that there are two different singly protonated OMID species bound directly to the metallic carbon steel substrate. In sharp contrast, in 0.01 M H2SO4, OMID adsorbs onto an ultra-thin surface film, composed primarily of a ferric sulfate (Fe2(SO4)3)-like phase. Such insight is essential to efforts to develop a mechanistic description of corrosion inhibitor functionality, as well as knowledge-based identification of next generation corrosion inhibitors.

Monomolecular Cracking of Propane: Effect of Zeolite Confinement and Acidity.

The effect of zeolite pore geometry and intrinsic acidity on the activation energy of propane monomolecular cracking was investigated for six topologically distinct zeolites with different pore sizes. Periodic density functional theory calculations were used to calculate the activation energy, while cluster models were used to calculate deprotonation energies. The computed intrinsic activation energies showed a smaller variation with topology than the adsorption energies. No correlation was found between the computed deprotonation and ammonia adsorption energies at the acid site and the intrinsic activation energy. Detailed analysis of the computed structures and properties suggests that acid sites with different pore topologies impose geometrical constraints on the ion-pair formed by the ammonium molecule, which differs significantly from those that affect the propane reaction.

The CRYSTAL code, 1976-2020 and beyond, a long story.

CRYSTAL is a periodic ab initio code that uses a Gaussian-type basis set to express crystalline orbitals (i.e., Bloch functions). The use of atom-centered basis functions allows treating 3D (crystals), 2D (slabs), 1D (polymers), and 0D (molecules) systems on the same grounds. In turn, all-electron calculations are inherently permitted along with pseudopotential strategies. A variety of density functionals are implemented, including global and range-separated hybrids of various natures and, as an extreme case, Hartree-Fock (HF). The cost for HF or hybrids is only about 3-5 times higher than when using the local density approximation or the generalized gradient approximation. Symmetry is fully exploited at all steps of the calculation. Many tools are available to modify the structure as given in input and simplify the construction of complicated objects, such as slabs, nanotubes, molecules, and clusters. Many tensorial properties can be evaluated by using a single input keyword: elastic, piezoelectric, photoelastic, dielectric, first and second hyperpolarizabilities, etc. The calculation of infrared and Raman spectra is available, and the intensities are computed analytically. Automated tools are available for the generation of the relevant configurations of solid solutions and/or disordered systems. Three versions of the code exist: serial, parallel, and massive-parallel. In the second one, the most relevant matrices are duplicated on each core, whereas in the third one, the Fock matrix is distributed for diagonalization. All the relevant vectors are dynamically allocated and deallocated after use, making the code very agile. CRYSTAL can be used efficiently on high performance computing machines up to thousands of cores.

Corrosion Protection through Naturally Occurring Films: New Insights from Iron Carbonate.

Despite intensive study over many years, the chemistry and physics of the atomic level mechanisms that govern corrosion are not fully understood. In particular, the occurrence and severity of highly localized metal degradation cannot currently be predicted and often cannot be rationalized in failure analysis. We report a first-principles model of the nature of protective iron carbonate films coupled with a detailed chemical and physical characterization of such a film in a carefully controlled environment. The fundamental building blocks of the protective film, siderite (FeCO3) crystallites, are found to be very sensitive to the growth environment. In iron-rich conditions, cylindrical crystallites form that are highly likely to be more susceptible to chemical attack and dissolution than the rhombohedral crystallites formed in iron-poor conditions. This suggests that local degradation of metal surfaces is influenced by structures that form during early growth and provides new avenues for the prevention, detection, and mitigation of carbon steel corrosion.

First Principles Calculations on the Stoichiometric and Defective (101) Anatase Surface and Upon Hydrogen and H2Pc Adsorption: The Influence of Electronic Exchange and Correlation and of Basis Set Approximations.

Anatase TiO2 provides photoactivity with high chemical stability at a reasonable cost. Different methods have been used to enhance its photocatalytic activity by creating band gap states through the introduction of oxygen vacancies, hydrogen impurities, or the adorption of phthalocyanines, which are usually employed as organic dyes in dye-sensitized solar cells. Predicting how these interactions affect the electronic structure of anatase requires an efficient and robust theory. In order to document the efficiency and accuracy of commonly used approaches we have considered two widely used implementations of density functional theory (DFT), namely the all-electron linear combination of atomic orbitals (AE-LCAO) and the pseudo-potential plane waves (PP-PW) approaches, to calculate the properties of the stoichiometric and defective anatase TiO2 (101) surface. Hybrid functionals, and in particular HSE, lead to a computed band gap in agreement with that measured by using UV adsorption spectroscopy. When using PBE+U, the gap is underestimated by 20 % but the computed position of defect induced gap states relative to the conduction band minimum (CBM) are found to be in good agreement with those calculated using hybrid functionals. These results allow us to conclude that hybrid functionals based on the use of AE-LCAO provide an efficient and robust approach for predicting trends in the band gap and the position of gap states in large model systems. We extend this analysis to surface adsorption and use the AE-LCAO approach with the hybrid functional HSED3 to study the adsorption of the phthalocyanine H2Pc on anatase (101). Our results suggest that H2Pc prefers to be adsorbed on the surface Ti5c rows of anatase (101), in agreement with that seen in recent STM experiments on rutile (110).

Unravelling Some of the Structure-Property Relationships in Graphene Oxide at Low Degree of Oxidation.

Graphene oxide (GO) is a versatile 2D material whose properties can be tuned by changing the type and concentration of oxygen-containing functional groups attached to its surface. However, a detailed knowledge of the dependence of the chemo/physical features of this material on its chemical composition is largely unknown. We combine classical molecular dynamics and density functional theory simulations to predict the structural and electronic properties of GO at low degree of oxidation and suggest a revision of the Lerf-Klinowski model. We find that layer deformation is larger for samples containing high concentrations of epoxy groups and that correspondingly the band gap increases. Targeted chemical modification of the GO surface appears to be an effective route to tailor the electronic properties of the monolayer for given applications. Our simulations also show that the chemical shift of the C-1s XPS peak allows one to unambiguously characterize GO composition, resolving the peak attribution  uncertainty often encountered in experiments.

The atomistic structure of yttria stabilised zirconia at 6.7 mol%: an ab initio study.

Yttria stabilized zirconia (YSZ) is an important oxide ion conductor used in solid oxide fuel cells, oxygen sensing devices, and for oxygen separation. Doping pure zirconia (ZrO2) with yttria (Y2O3) stabilizes the cubic structure against phonon induced distortions and this facilitates high oxide ion conductivity. The local atomic structure of the dopant is, however, not fully understood. X-ray and neutron diffraction experiments have established that, for dopant concentrations below 40 mol% Y2O3, no long range order is established. A variety of local structures have been suggested on the basis of theoretical and computational models of dopant energetics. These studies have been restricted by the difficulty of establishing force field models with predictive accuracy or exploring the large space of dopant configurations with first principles theory. In the current study a comprehensive search for all symmetry independent configurations (2857 candidates) is performed for 6.7 mol% YSZ modelled in a 2 × 2 × 2 periodic supercell using gradient corrected density functional theory. The lowest energy dopant structures are found to have oxygen vacancy pairs preferentially aligned along the 〈210〉 crystallographic direction in contrast to previous results which have suggested that orientation along the 〈111〉 orientation is favourable. Analysis of the defect structures suggests that the Y3+-Ovac interatomic separation is an important parameter for determining the relative configurational energies. Current force field models are found to be poor predictors of the lowest energy structures. It is suggested that the energies from a simple point charge model evaluated at unrelaxed geometries is actually a better descriptor of the energy ordering of dopant structures. Using these observations a pragmatic procedure for identifying low energy structures in more complicated material models is suggested. Calculation of the oxygen vacancy migration activation energies within the lowest energy 〈210〉 oriented structures gives results consistent with experimental observations.

Growth of Epitaxial Oxide Thin Films on Graphene.

The transfer process of graphene onto the surface of oxide substrates is well known. However, for many devices, we require high quality oxide thin films on the surface of graphene. This step is not understood. It is not clear why the oxide should adopt the epitaxy of the underlying oxide layer when it is deposited on graphene where there is no lattice match. To date there has been no explanation or suggestion of mechanisms which clarify this step. Here we show a mechanism, supported by first principles simulation and structural characterisation results, for the growth of oxide thin films on graphene. We describe the growth of epitaxial SrTiO3 (STO) thin films on a graphene and show that local defects in the graphene layer (e.g. grain boundaries) act as bridge-pillar spots that enable the epitaxial growth of STO thin films on the surface of the graphene layer. This study, and in particular the suggestion of a mechanism for epitaxial growth of oxides on graphene, offers new directions to exploit the development of oxide/graphene multilayer structures and devices.

Investigation of the Switching Mechanism in TiO2-Based RRAM: A Two-Dimensional EDX Approach.

The next generation of nonvolatile memory storage may well be based on resistive switching in metal oxides. TiO2 as transition metal oxide has been widely used as active layer for the fabrication of a variety of multistate memory nanostructure devices. However, progress in their technological development has been inhibited by the lack of a thorough understanding of the underlying switching mechanisms. Here, we employed high-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) combined with two-dimensional energy dispersive X-ray spectroscopy (2D EDX) to provide a novel, nanoscale view of the mechanisms involved. Our results suggest that the switching mechanism involves redistribution of both Ti and O ions within the active layer combined with an overall loss of oxygen that effectively render conductive filaments. Our study shows evidence of titanium movement in a 10 nm TiO2 thin-film through direct EDX mapping that provides a viable starting point for the improvement of the robustness and lifetime of TiO2-based resistive random access memory (RRAM).

A hybrid-exchange density functional study of the bonding and electronic structure in bulk CuFeS2.

The geometric, electronic, and magnetic properties of bulk chalcopyrite CuFeS2 have been investigated using hybrid-exchange density functional theory calculations. The results are compared with available theoretical and experimental data. The theoretical description of the bonding and electronic structure in CuFeS2 is analyzed in detail and compared to those computed for chalcocite (CuS2) and greigite (Fe3S4). In particular, the behavior of the 3d electrons of Fe(3+) is discussed in terms of the Hubbard-Anderson model in the strongly correlated regime and found to be similarly described in both materials by an on-site Coulomb repulsion (U) of ∼8.9 eV and a transfer integral (t) of ∼0.3 eV.

Chemical Descriptors of Yttria-Stabilized Zirconia at Low Defect Concentration: An ab Initio Study.

Yttria-stabilized zirconia (YSZ) is an important oxide ion conductor with applications in solid oxide fuel cells (SOFCs) and oxygen sensing devices. Doping the cubic phase of zirconia (c-ZrO2) with yttria (Y2O3) is isoelectronic, as two Zr(4+) ions are replaced by two Y(3+) ions, plus a charge compensating oxygen vacancy (Ovac). Typical doping concentrations include 3, 8, 10, and 12 mol %. For these concentrations, and all below 40 mol %, no phase with long-range order has been observed in either X-ray or neutron diffraction experiments. The prediction of local defect structure and the interaction between defects is therefore of great interest. This has not been possible to date as the number of possible defect topologies is very large and to perform reliable total energy calculations for all of them would be prohibitively expensive. Previous theoretical studies have only considered a selection of representative structures. In this study, a comprehensive search for low-energy defect structures using a combined classical modeling and density functional theory approach is used to identify the low-energy isolated defect structures at the dilute limit, 3.2 mol %. Through analysis of energetics computed using the best available Born-Mayer-Huggins empirical potential model, a point charge model, DFT, and a local strain energy estimated in the harmonic approximation, the main chemical and physical descriptors that correlate to the low-energy DFT structures are discussed. It is found that the empirical potential model reproduces a general trend of increasing DFT energetics across a series of locally strain relaxed structures but is unreliable both in predicting some incorrect low-energy structures and in finding some metastable structures to be unstable. A better predictor of low-energy defect structures is found to be the total electrostatic energy of a simple point charge model calculated at the unrelaxed geometries of the defects. In addition, the strain relaxation energy is estimated effectively in the harmonic approximation to the imaginary phonon modes of undoped c-ZrO2 but is found to be unimportant in determining the low-energy defect structures. These results allow us to propose a set of easily computed descriptors that can be used to identify the low-energy YSZ defect structures, negating the combinatorial complexity and number of defect structures that need to be considered.

Diffraction of helium on MgO(100) surface calculated from first-principles.

In this work we simulate the diffraction peak intensities of He beams scattered on the MgO(100) surface from first principles. It turns out that diffraction peak intensities are extremely sensitive to the quality of the potential describing the He-MgO surface interaction. Achieving the required accuracy in first principles calculations is very challenging indeed. The present work describes a first principles protocol able to achieve very high accuracy for reasonable computational cost. This method is based on periodic local second-order Møller-Plesset perturbation theory where systematic corrections for basis set truncation and for high-order electronic correlation are introduced using coupled cluster calculations on finite model systems mimicking the target system. For the He-MgO system the requirements with respect to the level of theory are very high; it turns out that contributions from connected quadruple excitations are non-negligible. Here we demonstrate that using this protocol, it is possible to reach the accuracy in the He-MgO potential that is required to predict the observed He diffraction peak intensities.

High-temperature antiferromagnetism in molecular semiconductor thin films and nanostructures.

The viability of dilute magnetic semiconductors in applications is linked to the strength of the magnetic couplings, and room temperature operation is still elusive in standard inorganic systems. Molecular semiconductors are emerging as an alternative due to their long spin-relaxation times and ease of processing, but, with the notable exception of vanadium-tetracyanoethylene, magnetic transition temperatures remain well below the boiling point of liquid nitrogen. Here we show that thin films and powders of the molecular semiconductor cobalt phthalocyanine exhibit strong antiferromagnetic coupling, with an exchange energy reaching 100 K. This interaction is up to two orders of magnitude larger than in related phthalocyanines and can be obtained on flexible plastic substrates, under conditions compatible with routine organic electronic device fabrication. Ab initio calculations show that coupling is achieved via superexchange between the singly occupied a1g () orbitals. By reaching the key milestone of magnetic coupling above 77 K, these results establish quantum spin chains as a potentially useable feature of molecular films.

Structural reconstruction of the graphene monovacancy.

Two distinct configurations of the monovacancy in graphene have been observed using aberration-corrected transmission electron microscopy (AC-TEM) at 80 kV. The predicted lower energy asymmetric monovacancy (MV), exhibiting a Jahn-Teller reconstruction (r-MV), has been observed, but in addition, we have imaged instances of a symmetric monovacancy (s-MV). We have used geometric phase analysis (GPA) to quantitatively determine the strain in the lattice surrounding these two defect configurations and show that the Jahn-Teller reconstruction generates significant extra strain compared to the symmetric MV case. Density functional theory calculations demonstrate that our experimental images of the two different monovacancies show good agreement with both the low energy r-MV and the metastable structures.

Ab initio calculation of the MgO(100) interaction with He and Ne: a HF + MP2 and HF + MP2(B3LYP) comparison.

Second order Rayleigh Schrödinger perturbation theory is applied to calculate the correlation energy contribution to the London dispersion interaction to approximate the interaction of the He and Ne with the MgO(100) surface; single particle orbitals using either Hartree-Fock theory or hybrid-exchange density functional theory are used as the reference state.