Engineering surface dipoles on mixed conducting oxides with ultra-thin oxide decoration layers.

Improving materials for energy conversion and storage devices is deeply connected with an optimization of their surfaces and surface modification is a promising strategy on the way to enhance modern energy technologies. This study shows that surface modification with ultra-thin oxide layers allows for a systematic tailoring of the surface dipole and the work function of mixed ionic and electronic conducting oxides, and it introduces the ionic potential of surface cations as a readily accessible descriptor for these effects. The combination of X-ray photoelectron spectroscopy (XPS) and density functional theory (DFT) illustrates that basic oxides with a lower ionic potential than the host material induce a positive surface charge and reduce the work function of the host material and vice versa. As a proof of concept that this strategy is widely applicable to tailor surface properties, we examined the effect of ultra-thin decoration layers on the oxygen exchange kinetics of pristine mixed conducting oxide thin films in very clean conditions by means of in-situ impedance spectroscopy during pulsed laser deposition (i-PLD). The study shows that basic decorations with a reduced surface work function lead to a substantial acceleration of the oxygen exchange on the surfaces of diverse materials.

Nonzero spontaneous electric polarization in metals: novel predictive methods and applications.

Ferroelectricity in metals has advanced since the initial discovery of nonmagnetic ferroelectric-like metal LiOsO[Formula: see text], anchored in the Anderson and Blount prediction. However, evaluating the spontaneous electric polarization (SEP) of this metal has been hindered by experimental and theoretical obstacles. The experimental challenge arises from difficulties in switching polarization using an external electric field, while the theoretical limitation lies in existing methods applicable only to nonmetals. Zabalo and Stengel (Phys Rev Lett 126:127601, 2021, https://doi.org/10.1103/PhysRevLett.126.127601 ) addressed the experimental obstacle by proposing flexoelectricity as an alternative for practical polarization switching in LiOsO[Formula: see text], which requires a critical bending radius similar to BaTiO[Formula: see text]. In this study, we focus on resolving the theoretical obstacle by modifying the Berry phase and Wannier functions approaches within density functional theory plus modern theory of polarization. By employing these modifications, we calculate the SEP of LiOsO[Formula: see text], comparable to the polarization of BaTiO[Formula: see text]. We validate our predictions using various ways. This study confirms the coexistence of ferroelectricity and metallicity in this new class of ferroelectric-like metals. Moreover, by addressing the theoretical limitation and providing new insights into polarization properties, our study complements the experimental flexoelectricity proposal and opens avenues for further exploration and manipulation of polarization characteristics. The developed approaches, incorporating modified Berry phase and Wannier function techniques, offer promising opportunities for studying and designing novel materials, including bio- and nano-ferroelectric-like metals. This study contributes to the advancement of ferroelectricity in metals and provides a foundation for future research in this exciting field.

Origin of the success of mGGAs for bandgaps.

One of the well-known limitations of Kohn-Sham density functional theory is the tendency to strongly underestimate bandgaps. Meta-generalized gradient approximations (mGGAs), which include the kinetic energy density in the functional form, have been shown to significantly alleviate this deficiency. In this study, we explore the mechanisms responsible for this improvement from the angle of the underlying local densities. We find that the highest occupied and lowest unoccupied states are distinct in the space of the underlying descriptors. The gap opening is compared to a simple scaling of the local density approximation, and two mechanisms responsible for opening the mGGA gaps are identified. First of all, the relatively large negative derivative of the functional form with respect to reduced kinetic energy tends to elevate the lowest unoccupied state. Second, the curvature of functional, which ensures that it is bounded, tends to lower the highest occupied state. Remarkably, these two mechanisms are found to be transferable over a large and diverse database of compounds.

DFT + U Simulation of the X-ray Absorption Near-Edge Structure of Bulk UO2 and PuO2.

Hubbard U-corrected density functional theory within the periodic boundary condition model in the WIEN2k code is used to simulate the actinide LIII and O K edge X-ray absorption near-edge structure (XANES) for UO2 and PuO2. Spin-orbit coupling effects are included, as are possible excitonic effects using supercells with a core hole on one of the atoms. Our calculations yield spectra in excellent agreement with previous experiments and superior to previous simulations. Density of states analysis reveals the mechanism behind the XANES peaks: the main contribution to the U/Pu LIII edges comes from the U/Pu d states hybridized with O p states, while as expected, the O p states primarily determine the O K edges of both UO2 and PuO2. The O K edges also feature O p hybridizing with U/Pu d and f states in the low-energy region and with U/Pu s and p states for the higher-energy peaks.

Electronic and ionic effects of sulphur and other acidic adsorbates on the surface of an SOFC cathode material.

The effects of sulphur adsorbates and other typical solid oxide fuel cell (SOFC) poisons on the electronic and ionic properties of an SrO-terminated (La,Sr)CoO3 (LSC) surface and on its oxygen exchange kinetics have been investigated experimentally with near ambient pressure X-ray photoelectron spectroscopy (NAP-XPS), low energy ion scattering (LEIS) and impedance spectroscopy as well as computationally with density functional theory (DFT). The experiment shows that trace amounts of sulphur in measurement atmospheres form SO2- 4 adsorbates and strongly deactivate a pristine LSC surface. They induce a work function increase, indicating a changing surface potential and a surface dipole. DFT calculations reveal that the main participants in these charge transfer processes are not sub-surface transition metals, but surface oxygen atoms. The study further shows that sulphate adsorbates strongly affect oxygen vacancy formation energies in the LSC (sub-)surface, thus affecting defect concentrations and oxygen transport properties. To generalize these results, the investigation was extended to other acidic oxides which are technologically relevant as SOFC cathode poisons, such as CO2 and CrO3. The results unveil a clear correlation of work function changes and redistributed charge with the Smith acidity of the adsorbed oxide and clarify fundamental mechanistic details of atomic surface modifications. The impact of acidic adsorbates on various aspects of the oxygen exchange reaction rate is discussed in detail.

What is the optimal mGGA exchange functional for solids?

The space of generalized gradient approximation (GGA) and meta-GGA (mGGA) exchange approximations is systematically explored by training 25 new functionals to produce accurate lattice parameter, cohesive energy, and bandgap predictions. The trained functionals are used to reproduce exact constraints in a data-driven way and to understand the accuracy trade-off between the mentioned properties. The functionals are compared to notable mGGA functionals to analyze how changes in the enhancement factor maps influence the accuracy of predictions. Some of the trained functionals are found to perform on par with specialized functionals for bandgaps, while outperforming them on the other two properties. The error surface of our trained functionals can serve as a soft-limit of what mGGA functionals can achieve.

CO Adsorption and Disproportionation on Smooth and Defect-Rich Ir(111).

CO adsorption and dissociation on "perfect" and "defect-rich" Ir(111) surfaces were studied by a combination of surface-analytical techniques, including polarization-dependent (PPP and SSP) sum frequency generation (SFG) vibrational spectroscopy, low-energy electron diffraction (LEED), Auger electron spectroscopy, X-ray photoelectron spectroscopy (XPS), and density functional theory (DFT) calculations. CO was found to be ordered and tilted from the surface normal at high coverage on the "perfect" surface (e.g., θ = 30° at 0.70 ML), whereas it was less ordered and preferentially upright (θ = 4-10°) on the "defect-rich" surface for coverages of 0.55-0.70 ML. SFG, LEED, and XPS revealed that CO adsorption at low pressure/high temperature and high pressure/low temperature was reversible. In contrast, upon heating to ∼600 K in near mbar CO pressure, "perfect" and even more "defect-rich" Ir(111) surfaces were irreversibly modified by carbon deposits, which, according to DFT, result from CO disproportionation.

Elucidating the formation and active state of Cu co-catalysts for photocatalytic hydrogen evolution.

The design of active and selective co-catalysts constitutes one of the major challenges in developing heterogeneous photocatalysts for energy conversion applications. This work provides a comprehensive insight into thermally induced bottom-up generation and transformation of a series of promising Cu-based co-catalysts. We demonstrate that the volcano-type HER profile as a function of calcination temperature is independent of the type of the Cu precursor but is affected by changes in oxidation state and location of the copper species. Supported by DFT modeling, our data suggest that low temperature (<200 °C) treatments facilitate electronic communication between the Cu species and TiO2, which allows for a more efficient charge utilization and maximum HER rates. In contrast, higher temperatures (>200 °C) do not affect the Cu oxidation state, but induce a gradual, temperature-dependent surface-to-bulk diffusion of Cu, which results in interstitial, tetra-coordinated Cu+ species. The disappearance of Cu from the surface and the introduction of new defect states is associated with a drop in HER performance. This work examines electronic and structural effects that are in control of the photocatalytic activity and can be transferred to other systems for further advancing photocatalysis.

Bandgap of two-dimensional materials: Thorough assessment of modern exchange-correlation functionals.

The density-functional theory (DFT) approximations that are the most accurate for the calculation of bandgap of bulk materials are hybrid functionals, such as HSE06, the modified Becke-Johnson (MBJ) potential, and the GLLB-SC potential. More recently, generalized gradient approximations (GGAs), such as HLE16, or meta-GGAs, such as (m)TASK, have also proven to be quite accurate for the bandgap. Here, the focus is on two-dimensional (2D) materials and the goal is to provide a broad overview of the performance of DFT functionals by considering a large test set of 298 2D systems. The present work is an extension of our recent studies [T. Rauch, M. A. L. Marques, and S. Botti, Phys. Rev. B 101, 245163 (2020); Patra et al., J. Phys. Chem. C 125, 11206 (2021)]. Due to the lack of experimental results for the bandgap of 2D systems, G0W0 results were taken as reference. It is shown that the GLLB-SC potential and mTASK functional provide the bandgaps that are the closest to G0W0. Following closely, the local MBJ potential has a pretty good accuracy that is similar to the accuracy of the more expensive hybrid functional HSE06.

Efficient Band Structure Calculation of Two-Dimensional Materials from Semilocal Density Functionals.

The experimental and theoretical realization of two-dimensional (2D) materials is of utmost importance in semiconducting applications. Computational modeling of these systems with satisfactory accuracy and computational efficiency is only feasible with semilocal density functional theory methods. In the search for the most useful method in predicting the band gap of 2D materials, we assess the accuracy of recently developed semilocal exchange-correlation (XC) energy functionals and potentials. Though the explicit forms of exchange-correlation (XC) potentials are very effective against XC energy functionals for the band gap of bulk solids, their performance needs to be investigated for 2D materials. In particular, the LMBJ [J. Chem. Theory Comput.2020, 16, 2654] and GLLB-SC [Phys. Rev. B82, 2010, 115106] potentials are considered for their dominance in bulk band gap calculation. The performance of recently developed meta generalized gradient approximations, like TASK [Phys. Rev. Res.1, 2019, 033082] and MGGAC [Phys. Rev. B. 100, 2019, 155140], is also assessed. We find that the LMBJ potential constructed for 2D materials is not as successful as its parent functional, i.e., MBJ [Phys. Rev. Lett.102, 2009, 226401] in bulk solids. Due to a contribution from the derivative discontinuity, the band gaps obtained with GLLB-SC are in a certain number of cases, albeit not systematically, larger than those obtained with other methods, which leads to better agreement with the quasi-particle band gap obtained from the GW method. The band gaps obtained with TASK and MGGAC can also be quite accurate.

Pristine quantum criticality in a Kondo semimetal.

The observation of quantum criticality in diverse classes of strongly correlated electron systems has been instrumental in establishing ordering principles, discovering new phases, and identifying the relevant degrees of freedom and interactions. At focus so far have been insulators and metals. Semimetals, which are of great current interest as candidate phases with nontrivial topology, are much less explored in experiments. Here, we study the Kondo semimetal CeRu4Sn6 by magnetic susceptibility, specific heat, and inelastic neutron scattering experiments. The power-law divergence of the magnetic Grünesien ratio reveals that, unexpectedly, this compound is quantum critical without tuning. The dynamical energy over temperature scaling in the neutron response throughout the Brillouin zone and the temperature dependence of the static uniform susceptibility, indicate that temperature is the only energy scale in the criticality. Such behavior, which has been associated with Kondo destruction quantum criticality in metallic systems, could be generic in the semimetal setting.

Giant spontaneous Hall effect in a nonmagnetic Weyl-Kondo semimetal.

Nontrivial topology in condensed-matter systems enriches quantum states of matter to go beyond either the classification into metals and insulators in terms of conventional band theory or that of symmetry-broken phases by Landau's order parameter framework. So far, focus has been on weakly interacting systems, and little is known about the limit of strong electron correlations. Heavy fermion systems are a highly versatile platform to explore this regime. Here we report the discovery of a giant spontaneous Hall effect in the Kondo semimetal [Formula: see text] that is noncentrosymmetric but preserves time-reversal symmetry. We attribute this finding to Weyl nodes-singularities of the Berry curvature-that emerge in the immediate vicinity of the Fermi level due to the Kondo interaction. We stress that this phenomenon is distinct from the previously detected anomalous Hall effect in materials with broken time-reversal symmetry; instead, it manifests an extreme topological response that requires a beyond-perturbation-theory description of the previously proposed nonlinear Hall effect. The large magnitude of the effect in even tiny electric and zero magnetic fields as well as its robust bulk nature may aid the exploitation in topological quantum devices.

Ca-doped rare earth perovskite materials for tailored exsolution of metal nanoparticles.

Perovskite-type oxide materials (nominal composition ABO3) are a very versatile class of materials, and their properties are tuneable by varying and doping A- and B-site cations. When the B-site contains easily reducible cations (e.g. Fe, Co or Ni), these can exsolve under reducing conditions and form metallic nanoparticles on the surface. This process is very interesting as a novel route for the preparation of catalysts, since oxide surfaces decorated with finely dispersed catalytically active (often metallic) nanoparticles are a key requirement for excellent catalyst performance. Five doped perovskites, namely, La0.9Ca0.1FeO3-δ, La0.6Ca0.4FeO3-δ, Nd0.9Ca0.1FeO3-δ, Nd0.6Ca0.4FeO3-δ and Nd0.6Ca0.4Fe0.9Co0.1O3-δ, have been synthesized and characterized by experimental and theoretical methods with respect to their crystal structures, electronic properties, morphology and exsolution behaviour. All are capable of exsolving Fe and/or Co. Special emphasis has been placed on the influence of the A-site elemental composition on structure and exsolution capability. Using Nd instead of La increased structural distortions and, at the same time, hindered exsolution. Increasing the amount of Ca doping also increased distortions and additionally changed the Fe oxidation states, resulting in exsolution being shifted to higher temperatures as well. Using the easily reducible element Co as the B-site dopant significantly facilitated the exsolution process and led to much smaller and homogeneously distributed exsolved particles. Therefore, the Co-doped perovskite is a promising material for applications in catalysis, even more so as Co is catalytically a highly active element. The results show that fine-tuning of the perovskite composition will allow tailored exsolution of nanoparticles, which can be used for highly sophisticated catalyst design.

Coverage-Induced Orientation Change: CO on Ir(111) Monitored by Polarization-Dependent Sum Frequency Generation Spectroscopy and Density Functional Theory.

Polarization-dependent sum frequency generation (SFG) spectroscopy was applied to study the adsorption of carbon monoxide (CO) on the well-ordered (annealed) Ir(111) single-crystal surface at various CO coverages. Coverage was adjusted by varying the substrate temperature (300-575 K) and/or gas pressure (10-7 to 1.0 mbar). Under all conditions investigated, only a single absorption band at 2038-2094 cm-1 was observed, characteristic of linearly bonded (on-top) CO. Using different polarizations, PPP and SSP spectra were acquired with a high signal-to-noise ratio, whereby tilt angles of CO on Ir(111) could be determined for the first time by SFG. It was found that not only the vibrational frequency of on-top CO but also the tilt angle was strongly coverage-dependent. The higher the coverage was, the larger the vibrational frequency and the tilt angle were. At about 0.7 ML coverage, a CO tilt angle of at least 20° was observed, which is in good agreement with density functional theory (DFT) calculations. In addition, the molecular hyperpolarizability ratio (R) of CO (at 0.13 ML in UHV) was determined to be 0.08. Based on the combined SFG/DFT results, it may change to 0.29 at 0.77 ML coverage.

WIEN2k: An APW+lo program for calculating the properties of solids.

The WIEN2k program is based on the augmented plane wave plus local orbitals (APW+lo) method to solve the Kohn-Sham equations of density functional theory. The APW+lo method, which considers all electrons (core and valence) self-consistently in a full-potential treatment, is implemented very efficiently in WIEN2k, since various types of parallelization are available and many optimized numerical libraries can be used. Many properties can be calculated, ranging from the basic ones, such as the electronic band structure or the optimized atomic structure, to more specialized ones such as the nuclear magnetic resonance shielding tensor or the electric polarization. After a brief presentation of the APW+lo method, we review the usage, capabilities, and features of WIEN2k (version 19) in detail. The various options, properties, and available approximations for the exchange-correlation functional, as well as the external libraries or programs that can be used with WIEN2k, are mentioned. References to relevant applications and some examples are also given.

Local geometry around B atoms in B/Si(1 1 1) from polarized x-ray absorption spectroscopy.

The arrangement of B atoms in a doped Si(1 1 1)-[Formula: see text]:B system was studied using a near-edge x-ray absorption fine structure (NEXAFS). Boron atoms were deposited via segregation from the bulk by flashing the sample repeatedly. The positions of B atoms are determined by comparing measured polarized (angle-dependent) NEXAFS spectra with spectra calculated for various structural models based on ab initio total energy calculations. It is found that most of boron atoms are located in sub-surface L[Formula: see text] positions, beneath a Si atom. However, depending on the preparation method a significant portion of B atoms may be located elsewhere. A possible location of these non-L[Formula: see text]-atoms is at the surface, next to those Si atoms which form the [Formula: see text] reconstruction.

On the calculation of the bandgap of periodic solids with MGGA functionals using the total energy.

During the last few years, it has become more and more clear that functionals of the meta generalized gradient approximation (MGGA) are more accurate than GGA functionals for the geometry and energetics of electronic systems. However, MGGA functionals are also potentially more interesting for the electronic structure, in particular, when the potential is nonmultiplicative (i.e., when MGGAs are implemented in the generalized Kohn-Sham framework), which may help to get more accurate bandgaps. Here, we show that the calculation of bandgap of solids with MGGA functionals can also be done very accurately in a non-self-consistent manner. This scheme uses only the total energy and can, therefore, be very useful when the self-consistent implementation of a particular MGGA functional is not available. Since self-consistent MGGA calculations may be difficult to converge, the non-self-consistent scheme may also help to speed up the calculations. Furthermore, it can be applied to any other types of functionals, for which the implementation of the corresponding potential is not trivial.

Comparative study of the PBE and SCAN functionals: The particular case of alkali metals.

The SCAN meta-generalized gradient approximation (GGA) functional is known to describe multiple properties of various materials with different types of bonds with greater accuracy, compared to the widely used PBE GGA functional. Yet, for alkali metals, SCAN shows worse agreement with experimental results than PBE despite using more information about the system. In the current study, this behavior for alkali metals is explained by identifying an inner semicore region which, within SCAN, contributes to an underbinding. The inner semicore push toward larger lattice constants is a general feature but is particularly important for very soft materials, such as the alkali metals, while for harder materials the valence region dominates.

Orbital-free approximations to the kinetic-energy density in exchange-correlation MGGA functionals: Tests on solids.

A recent study of Mejia-Rodriguez and Trickey [Phys. Rev. A 96, 052512 (2017)] showed that the deorbitalization procedure (replacing the exact Kohn-Sham kinetic-energy density by an approximate orbital-free expression) applied to exchange-correlation functionals of the meta-generalized gradient approximation (MGGA) can lead to important changes in the results for molecular properties. For the present work, the deorbitalization of MGGA functionals is further investigated by considering various properties of solids. It is shown that depending on the MGGA, common orbital-free approximations to the kinetic-energy density can be sufficiently accurate for the lattice constant, bulk modulus, and cohesive energy. For the bandgap, calculated with the modified Becke-Johnson MGGA potential, the deorbitalization has a larger impact on the results.