Ab initio calculations of photocurrent in solids based on Wannier function

We present a general ab initio method based on Wannier function 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 relaxation time approximation), that is to say, it includes all intra-band, inter-band and their cross terms. Specifically for the second-order DC photocurrent, it includes all of the following contributions - shift current, (magnetic) injection current, Berry curvature dipole, gyration current and Fermi surface ones, instead of only a part of them as in most previous ab initio methods. 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 WS₂, a 2D ferroelectric material - monolayer GeS, bilayer anti-ferromagnetic MnBi₂Te₄ and topological Weyl semimetal RhSi to simulate their charge and/or spin, DC and/or AC photocurrent. Our theoretical results are in good agreement with previous theoretical works. Our method provides a universal computational tool for reliable and accurate predictions of abundant weak-field photocurrent phenomena in disparate materials.