Transition metal oxides hold great potential for the development of new device paradigms because of the field-tunable functionalities driven by their strong electronic correlations, combined with their earth abundance and environmental friendliness. Recently, the interfaces between transition-metal oxides have revealed striking phenomena, such as insulator-metal transitions, magnetism, magnetoresistance and superconductivity. Such oxide interfaces are usually produced by sophisticated layer-by-layer growth techniques, which can yield high-quality, epitaxial interfaces with almost monolayer control of atomic positions. The resulting interfaces, however, are fixed in space by the arrangement of the atoms. Here we demonstrate a route to overcoming this geometric limitation. We show that the electrical conductance at the interfacial ferroelectric domain walls in hexagonal ErMnO(3) is a continuous function of the domain wall orientation, with a range of an order of magnitude. We explain the observed behaviour using first-principles density functional and phenomenological theories, and relate it to the unexpected stability of head-to-head and tail-to-tail domain walls in ErMnO(3) and related hexagonal manganites. As the domain wall orientation in ferroelectrics is tunable using modest external electric fields, our finding opens a degree of freedom that is not accessible to spatially fixed interfaces.

Systematic and automatic calculations of the electronic band structure are a crucial component of computationally-driven high-throughput materials screening. An algorithm, for any crystal, to derive a unique description of the crystal structure together with a recommended band path is indispensable for this task. The electronic band structure is typically sampled along a path within the first Brillouin zone including the surface in reciprocal space. Some points in reciprocal space have higher site symmetries and/or have higher constraints than other points regarding the electronic band structure and therefore are likely to be more important than other points. This work categorizes 2 points in reciprocal space according to their symmetry and provides recommended band paths that cover all special wavevector ( k -vector) points and lines necessarily and sufficiently. Points in reciprocal space are labeled such that there is no conflict with the crystallographic convention. The k -vector coefficients of labeled points, which are located at Brillouin zone face and edge centers as well as vertices, are derived based on a primitive cell compatible with the crystallographic convention, including those with axial ratio-dependent coordinates. Furthermore, we provide an open-source implementation of the algorithms within our SeeK-path python code, to allow researchers to obtain k -vector coefficients and recommended band paths in an automated fashion. Finally, we created a free online service to compute and visualise Brillouin Zone, labeled k -points and suggested band paths for any crystal structure, that we made available at http://www.materialscloud.org/tools/seekpath/ .

Finite-size corrections for charged defect supercell calculations typically consist of image-charge and potential alignment corrections. A wide variety of schemes for both corrections have been proposed for decades. Regarding the image-charge correction, Freysoldt, Neugebauer, and Van de Walle (FNV) recently proposed a novel method that enables us to accurately estimate the correction energy a posteriori through alignment of the defect-induced potential to the model charge potential [C. Freysoldt, J. Neugebauer, and C. G. Van de Walle, Phys. Rev. Lett. 102, 016402 (2009).] This method, however, still has two issues in practice. Firstly, it uses planar-averaged electrostatic potential for determining the potential offset, which cannot be readily applied to relaxed atomic structure. Secondly, the long-range Coulomb interaction is assumed to be screened by a macroscopic dielectric constant. This is valid only for cubic systems and can bring forth huge errors for defects in anisotropic materials, particularly with layered and low-dimensional structures. In the present study, we use the atomic site electrostatic potential as a potential marker instead of the planar-averaged potential, and extend the FNV scheme by adopting the point charge model in an anisotropic medium for estimating long-range interactions. We also revisit the conventional potential alignment correction and show that it is fully included in the image-charge correction and therefore unnecessary. In addition, we show that the potential alignment corresponds to a part of first-order and full of third-order image-charge correction; thus the third-order imagecharge contribution is absent after the potential alignment. Finally, a systematic assessment of the accuracy of the extended FNV correction scheme is performed for a wide range of material classes: β-Li 2 TiO 3 , ZnO, MgO, corundum Al 2 O 3 , monoclinic HfO 2 , cubic and hexagonal BN, Si, GaAs, and diamond. The defect formation energies with -6 to +3 charges calculated using around 100-atom supercells are successfully corrected even after atomic relaxation within a few tenths of eV compared to those in the dilute limit.

Aerobic oxidation of 5-hydroxymethylfurfural (HMF) to 2,5-furandicarboxylic acid (FDCA) as a bioplastics monomer is efficiently promoted by a simple system based on a nonprecious-metal catalyst of MnO2 and NaHCO3. Kinetic studies indicate that the oxidation of 5-formyl-2-furancarboxylic acid (FFCA) to FDCA is the slowest step for the aerobic oxidation of HMF to FDCA over activated MnO2. We demonstrate through combined computational and experimental studies that HMF oxidation to FDCA is largely dependent on the MnO2 crystal structure. Density functional theory (DFT) calculations reveal that vacancy formation energies at the planar oxygen sites in α- and γ-MnO2 are higher than those at the bent oxygen sites. β- and λ-MnO2 consist of only planar and bent oxygen sites, respectively, with lower vacancy formation energies. Consequently, β- and λ-MnO2 are likely to be good candidates as oxidation catalysts. On the other hand, experimental studies reveal that the reaction rates per surface area for the slowest step (FFCA oxidation to FDCA) decrease in the order of β-MnO2 > λ-MnO2 > γ-MnO2 ≈ α-MnO2 > δ-MnO2 > ε-MnO2; the catalytic activity of β-MnO2 exceeds that of the previously reported activated MnO2 by three times. The order is in good agreement not only with the DFT calculation results, but also with the reduction rates per surface area determined by the H2-temperature-programmed reduction measurements for MnO2 catalysts. The successful synthesis of high-surface-area β-MnO2 significantly improves the catalytic activity for the aerobic oxidation of HMF to FDCA.