The transition metal dipnictides TaAs2, TaSb2, NbAs2 and NbSb2 have recently sparked interest for exhibiting giant magnetoresistance. While the exact nature of the magnetoresistance in these materials is still under active investigation, there are experimental results indicating that it is of the anisotropic negative variety. We study the effect of magnetic fields on the band structure topology of these materials by applying Zeeman splitting. In the absence of a magnetic field, we find that the materials are weak topological insulators, which is in agreement with previous studies. When the magnetic field is applied, we find that type-II Weyl points form. This result is found first from a symmetry argument, and then numerically for a model of TaAs2 and a tight-binding model of NbSb2. This effect could be of help in the search for an explanation of the anomalous magnetoresistance in these materials.
We present a general technique for capturing various non-trivial topologies in the band structure of materials, which often arise from spin-orbit coupling. The technique is aimed at insulators and semimetals. Of insulators, Chern, Z2, and crystalline topological insulators can be identified. Of semimetals, the technique captures non-trivial topologies associated with the presence of Weyl and Dirac points in the spectrum. A public software package -- Z2Pack -- based on this technique will be presented. Z2Pack is an easy-to-use, well documented Python package that computes topological invariants and illustrates non-trivial features of Berry curvature. It works as a post-processing tool with all major first-principles codes, as well as with tight-binding models. As such, it can be used to investigate materials with strong spin-orbit coupling.
The electronic structure of surfaces plays a key role in the properties of quantum devices. However, surfaces are also the most challenging to simulate and engineer. Here, we study the electronic structure of InAs(001), InAs(111), and InSb(110) surfaces using a combination of density functional theory (DFT) and angle-resolved photoemission spectroscopy (ARPES). We were able to perform large-scale first principles simulations and capture effects of different surface reconstructions by using DFT calculations with a machine-learned Hubbard U correction [npj Comput. Mater. 6, 180 (2020)]. To facilitate direct comparison with ARPES results, we implemented a "bulk unfolding" scheme by projecting the calculated band structure of a supercell surface slab model onto the bulk primitive cell. For all three surfaces, we find a good agreement between DFT calculations and ARPES. For InAs(001), the simulations clarify the effect of the surface reconstruction. Different reconstructions are found to produce distinctive surface states. For InAs(111) and InSb(110), the simulations help elucidate the effect of oxidation. Owing to larger charge transfer from As to O than from Sb to O, oxidation of InAs(111) leads to significant band bending and produces an electron pocket, whereas oxidation of InSb(110) does not. Our combined theoretical and experimental results may inform the design of quantum devices based on InAs and InSb semiconductors, e.g., topological qubits utilizing the Majorana zero modes.
The prediction of material properties through electronic-structure simulations based on density-functional theory has become routinely common, thanks, in part, to the steady increase in the number and robustness of available simulation packages. This plurality of codes and methods aiming to solve similar problems is both a boon and a burden. While providing great opportunities for cross-verification, these packages adopt different methods, algorithms, and paradigms, making it challenging to choose, master, and efficiently use any one for a given task. Leveraging recent advances in managing reproducible scientific workflows, we demonstrate how developing common interfaces for workflows that automatically compute material properties can tackle the challenge mentioned above, greatly simplifying interoperability and cross-verification. We introduce design rules for reproducible and reusable code-agnostic workflow interfaces to compute well-defined material properties, which we implement for eleven different quantum engines and use to compute three different material properties. Each implementation encodes carefully selected simulation parameters and workflow logic, making the implementer's expertise of the quantum engine directly available to non-experts. Full provenance and reproducibility of the workflows is guaranteed through the use of the AiiDA infrastructure. All workflows are made available as open-source and come pre-installed with the Quantum Mobile virtual machine, making their use straightforward.
Fermions--elementary particles such as electrons--are classified as Dirac, Majorana or Weyl. Majorana and Weyl fermions had not been observed experimentally until the recent discovery of condensed matter systems such as topological superconductors and semimetals, in which they arise as low-energy excitations. Here we propose the existence of a previously overlooked type of Weyl fermion that emerges at the boundary between electron and hole pockets in a new phase of matter. This particle was missed by Weyl because it breaks the stringent Lorentz symmetry in high-energy physics. Lorentz invariance, however, is not present in condensed matter physics, and by generalizing the Dirac equation, we find the new type of Weyl fermion. In particular, whereas Weyl semimetals--materials hosting Weyl fermions--were previously thought to have standard Weyl points with a point-like Fermi surface (which we refer to as type-I), we discover a type-II Weyl point, which is still a protected crossing, but appears at the contact of electron and hole pockets in type-II Weyl semimetals. We predict that WTe2 is an example of a topological semimetal hosting the new particle as a low-energy excitation around such a type-II Weyl point. The existence of type-II Weyl points in WTe2 means that many of its physical properties are very different to those of standard Weyl semimetals with point-like Fermi surfaces.
Abstract W annier90 is an open-source computer program for calculating maximally-localised Wannier functions (MLWFs) from a set of Bloch states. It is interfaced to many widely used electronic-structure codes thanks to its independence from the basis sets representing these Bloch states. In the past few years the development of W annier90 has transitioned to a community-driven model; this has resulted in a number of new developments that have been recently released in W annier90 v3.0. In this article we describe these new functionalities, that include the implementation of new features for wannierisation and disentanglement (symmetry-adapted Wannier functions, selectively-localised Wannier functions, selected columns of the density matrix) and the ability to calculate new properties (shift currents and Berry-curvature dipole, and a new interface to many-body perturbation theory); performance improvements, including parallelisation of the core code; enhancements in functionality (support for spinor-valued Wannier functions, more accurate methods to interpolate quantities in the Brillouin zone); improved usability (improved plotting routines, integration with high-throughput automation frameworks), as well as the implementation of modern software engineering practices (unit testing, continuous integration, and automatic source-code documentation). These new features, capabilities, and code development model aim to further sustain and expand the community uptake and range of applicability, that nowadays spans complex and accurate dielectric, electronic, magnetic, optical, topological and transport properties of materials.
