High performance Wannier interpolation of Berry curvature and related quantities with WannierBerri code

Wannier interpolation is a powerful tool for performing Brillouin zone integrals over dense grids ofkpoints, which are essential to evaluate such quantities as the intrinsic anomalous Hall conductivity or Boltzmann transport coefficients. However, more complex physical problems and materials create harder numerical challenges, and computations with the existing codes become very expensive, which often prevents reaching the desired accuracy. In this article, I present a series of methods that boost the speed of Wannier interpolation by several orders of magnitude. They include a combination of fast and slow Fourier transforms, explicit use of symmetries, and recursive adaptive grid refinement among others. The proposed methodology has been implemented in the python code WannierBerri, which also aims to serve as a convenient platform for the future development of interpolation schemes for other phenomena

Ab initio calculation of the shift photocurrent by Wannier interpolation

We describe and implement a first-principles algorithm based on maximally-localized Wannier functions for calculating the shift-current response of piezoelectric crystals in the independent-particle approximation. The proposed algorithm presents several advantages over existing ones, including full gauge invariance, low computational cost, and a correct treatment of the optical matrix elements with nonlocal pseudopotentials. Band-truncation errors are avoided by a careful formulation of $k\cdot p$ perturbation theory within the subspace of wannierized bands. The needed ingredients are the matrix elements of the Hamiltonian and of the position operator in the Wannier basis, which are readily available at the end of the wannierization step. If the off-diagonal matrix elements of the position operator are discarded, our expressions reduce to the ones that have been used in recent tight-binding calculations of the shift current. We find that this `diagonal' approximation can introduce sizeable errors, highlighting the importance of carefully embedding the tight-binding model in real space for an accurate description of the charge transfer that gives rise to the shift current.Comment: 13 pages, 7 figure

On the separation of Hall and Ohmic nonlinear responses

The symmetric and antisymmetric parts of the linear conductivity describe the dissipative (Ohmic) and nondissipative (Hall) parts of the current. The Hall current is always transverse to the applied electric field regardless of its orientation; the Ohmic current is purely longitudinal in cubic crystals, but in lower-symmetry crystals it has a transverse component whenever the field is not aligned with a principal axis. In this work, we extend that analysis beyond the linear regime. We consider all possible ways of partitioning the current at any order in the electric field without taking symmetry into account, and find that the Hall vs Ohmic decomposition is the only one that satisfies certain basic requirements. A general prescription is given for achieving that decomposition, and the case of the quadratic conductivity is analyzed in detail. By performing a symmetry analysis we find that in five of the 122 magnetic point groups the quadratic dc conductivity is purely Ohmic and even under time reversal, a type of response that is entirely disorder mediated

Unfolding spinor wavefunctions and expectation values of general operators: Introducing the unfolding-density operator

We show that the spectral weights $W_{m\vec K}(\vec k)$ used for the unfolding of two-component spinor eigenstates $| {\psi_{m\vec K}^\mathrm{SC}} > = | \alpha > | {\psi_{m\vec{K}}^\mathrm{SC, \alpha}} > + | \beta > | {\psi_{m\vec{K}}^\mathrm{SC, \beta}} >$ can be decomposed as the sum of the partial spectral weights $W_{m\vec{K}}^{\mu}(\vec k)$ calculated for each component $\mu = \alpha, \beta$ independently, effortlessly turning a possibly complicated problem involving two coupled quantities into two independent problems of easy solution. Furthermore, we define the unfolding-density operator $\hat{\rho}_{\vec{K}}(\vec{k}_{i}; \, \varepsilon)$, which unfolds the primitive cell expectation values $\varphi^{pc}(\vec{k}; \varepsilon)$ of any arbitrary operator $\mathbf{\hat\varphi}$ according to $\varphi^{pc}(\vec{k}_{i}; \varepsilon) = \mathit{Tr}(\hat{\rho}_{\vec{K}}(\vec{k}_{i}; \, \varepsilon)\,\,\hat{\varphi})$. As a proof of concept, we apply the method to obtain the unfolded band structures, as well as the expectation values of the Pauli spin matrices, for prototypical physical systems described by two-component spinor eigenfunctions

Multi-band nodal links in triple-point materials

We study a class of topological materials which in their momentum-space band structure exhibit three-fold degeneracies known as triple points. Specifically, we investigate and classify triple points occurring along high-symmetry lines of $\mathcal{P}\mathcal{T}$-symmetric crystalline solids with negligible spin-orbit coupling. By employing the recently discovered non-Abelian band topology, we argue that a rotation-symmetry-breaking strain transforms a certain class of triple points into multi-band nodal links. Although multi-band nodal-line compositions were previously theoretically conceived, a practical condensed-matter platform for their manipulation and inspection has hitherto been missing. By reviewing the known triple-point materials in the considered symmetry class, and by performing first-principles calculations to predict new ones, we identify suitable candidates for the realization of multi-band nodal links. In particular, we find that Li$_2$NaN is an ideal compound to study this phenomenon, where the band nodes facilitate largely tunable density of states and optical conductivity with doping and strain, respectively. The multi-band linking is expected to equip the nodal rings with monopole charges, making such triple-point materials a versatile platform to probe the non-Abelian band topology.Comment: 4 pages (3 figures, 1 table) + 13 pages of Supplemental Material (13 figures, 3 tables) + reference

Ab initio study of the nonlinear optical properties and d.c. photocurrent of the Weyl semimetal TaIrTe4_4

We present a first principles theoretical study employing nonlinear response theory to investigate the d.c. photocurrent generated by linearly polarized light in the type-II Weyl semimetal TaIrTe4. We report the low energy spectrum of several nonlinear optical effects. At second-order, we consider the shift and injection currents. Assuming the presence of a built-in static electric field, at third-order we study the current-induced shift and injection currents, as well as the jerk current. We discuss our results in the context of a recent experiment measuring an exceptionally large photoconductivity in this material [J. Ma et at., Nat. Mater. 18, 476 (2019)]. According to our results, the jerk current is the most likely origin of the large response. Finally, we propose means to discern the importance of the various mechanisms involved in a time-resolved experiment.Comment: 9 pages, 11 figure

Topological zero-dimensional defect and flux states in three-dimensional insulators

In insulating crystals, it was previously shown that defects with two fewer dimensions than the bulk can bind topological electronic states. We here further extend the classification of topological defect states by demonstrating that the corners of crystalline defects with integer Burgers vectors can bind 0D higher-order end (HEND) states with anomalous charge and spin. We demonstrate that HEND states are intrinsic topological consequences of the bulk electronic structure and introduce new bulk topological invariants that are predictive of HEND dislocation states in solid-state materials. We demonstrate the presence of first-order 0D defect states in PbTe monolayers and HEND states in 3D SnTe crystals. We relate our analysis to magnetic flux insertion in insulating crystals. We find that π-flux tubes in inversion- and time-reversal-symmetric (helical) higher-order topological insulators bind Kramers pairs of spin-charge-separated HEND states, which represent observable signatures of anomalous surface half quantum spin Hall states