Within density functional theory (DFT), adding a Hubbard U correction can mitigate some of the deficiencies of local and semi-local exchange-correlation functionals, while maintaining computational efficiency. However, the accuracy of DFT+U largely depends on the chosen Hubbard U values. We propose an approach to determining the optimal U parameters for a given material by machine learning. The Bayesian optimization (BO) algorithm is used with an objective function formulated to reproduce the band structures produced by more accurate hybrid functionals. This approach is demonstrated for transition metal oxides, europium chalcogenides, and narrow-gap semiconductors. The band structures obtained using the BO U values are in agreement with hybrid functional results. Additionally, comparison to the linear response (LR) approach to determining U demonstrates that the BO method is superior.

We present Ogre, an open-source code for generating surface slab models from bulk molecular crystal structures. Ogre is written in Python and interfaces with the FHI-aims code to calculate surface energies at the level of density functional theory (DFT). The input of Ogre is the geometry of the bulk molecular crystal. The surface is cleaved from the bulk structure with the molecules on the surface kept intact. A slab model is constructed according to the user specifications for the number of molecular layers and the length of the vacuum region. Ogre automatically identifies all symmetrically unique surfaces for the user-specified Miller indices and detects all possible surface terminations. Ogre includes utilities to analyze the surface energy convergence and Wulff shape of the molecular crystal. We present the application of Ogre to three representative molecular crystals: the pharmaceutical aspirin, the organic semiconductor tetracene, and the energetic material HMX. The equilibrium crystal shapes predicted by Ogre are in agreement with experimentally grown crystals, demonstrating that DFT produces satisfactory predictions of the crystal habit for diverse classes of molecular crystals.

No abstract

We use density functional theory calculations to study the electronic structure of epitaxial (111) interfaces of the topological crystalline insulators SnSe and SnTe with the magnetic insulator EuS and the non-magnetic insulator CaTe, respectively. We consider both interface slab models with a vacuum region and periodic heterostructures without vacuum. We find that gaps of 21 meV at the Γ point and 9 meV at the M point arise in the topological state at the SnSe/EuS interface, due to the magnetic proximity effect, which breaks the time reversal symmetry. The surface state at Γ is shifted below the Fermi level by 88 meV and the surface state at M is shifted above the Fermi level by 47 meV, owing to band bending at the interface. By comparison, the topological state at the interface of SnTe/CaTe is unperturbed by the presence of non-magnetic CaTe.

We present a new version of the Ogre open source Python package with the capability to perform structure prediction of epitaxial inorganic interfaces by lattice and surface matching. In the lattice matching step a scan over combinations of substrate and film Miller indices is performed to identify the domain-matched interfaces with the lowest mismatch. Subsequently, surface matching is conducted by Bayesian optimization to find the optimal interfacial distance and inplane registry between the substrate and film. For the objective function, a geometric score function is proposed, based on the overlap and empty space between atomic spheres at the interface. The score function reproduces the results of density functional theory (DFT) at a fraction of the computational cost. The optimized interfaces are pre-ranked using a score function based on the similarity of the atomic environment at the interface to the bulk environment. Final ranking of the top candidate structures is performed with DFT. Ogre streamlines DFT calculations of interface energies and electronic properties by automating the construction of interface models. The application of Ogre is demonstrated for two interfaces of interest for quantum computing and spintronics, Al/InAs and Fe/InSb.

No abstract