Abstract:We present a general ab initio method based on Wannier functions using the covariant derivative for simulating the photocurrent in solids. The method is widely applicable to charge/spin DC and AC photocurrent at any perturbation levels in both semiconductors and metals for both linearly and circularly polarized light. This is because the method is theoretically complete (within the relaxation time approximation), that is to say, it includes all intraband, interband and their cross terms. It is also free from the degeneracy issue, i.e., applicable to arbitrary band structures with arbitrary numbers of degenerate bands. We apply the method to various semiconductors and metals, including GaAs, graphene-hBN heterostructure, monolayer WS2, a 2D ferroelectric material - monolayer GeS, bilayer anti-ferromagnetic MnBi2Te4 and topological Weyl semimetal RhSi, to simulate their charge and/or spin, DC and/or AC photocurrent. Our theoretical results are in agreement with previous theoretical works. Our numerical tests of GaAs, WS2 and GeS suggest setting the degeneracy threshold in the conventional method as \hbar\Gamma^{2}, with \Gamma^{2} the relaxation rate of the off-diagonal elements of the density matrix between two states with close energies. We find that compared with the conventional Wannier-function-based method using non-dgenerate perturbation theory, the numerical errors of optical susceptibilities of bilayer anti-ferromagnetic MnBi2Te4 with the PT symmetry can be reduced by 1-2 orders of magnitude by our method for circularly polarized light. Our method provides a universal computational tool for reliable and accurate predictions of abundant weak-field photocurrent phenomena in disparate materials.

Wannier Interpolation with Nonorthogonal Localized Orbitals: Application to ab initio Calculation of Nonlinear Optical Responses.

First-principles calculation of optical responses based on nonorthogonal localized orbitals

High-order and adaptive optical conductivity calculations using Wannier interpolation

Wannier function localisation using Bloch intrinsic atomic orbitals

Maximally localized Wannier functions: Theory and applications

A highly accurate procedure for computing globally optimal Wannier functions in one-dimensional crystalline insulators

Electronic Structure Interpolation Via Atomic Orbitals

First-principles calculation of orbital Hall effect by Wannier interpolation: Role of orbital dependence of the anomalous position

Wannier interpolation of reciprocal-space periodic and non-periodic matrix elements in the optimally smooth subspace

All-electron periodic G 0 W 0 implementation with numerical atomic orbital basis functions: Algorithm and benchmarks

Electronic and Optical Excitations in van der Waals Materials from a Non-Empirical Wannier-Localized Optimally-Tuned Screened Range-Separated Hybrid Functional

First-Principles Calculations of Conductance Within A Plane Wave Basis Set Via Non-Orthogonal Wannier-Type Atomic Orbitals

\textit{Ab initio} Wannier-representation-based calculations of photocurrent in semiconductors and metals

An unambiguous and robust formulation for Wannier localization

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

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

Automated mixing of maximally localized Wannier functions into target manifolds

Closest Wannier functions to a given set of localized orbitals

Linear Scaling Approach for Optical Excitations Using Maximally Localized Wannier Functions

Excitonic effects in nonlinear optical responses: Exciton-state formalism and first-principles calculations

Obtaining Wannier functions of a crystalline insulator within a Hartree-Fock approach: Applications to LiF and LiCl