I was trying to calculate the intrinsic resistivity of a metal (as an extension of fcc Pb in tutorial) and some puzzles appeared in these days.
In my understanding, the total (isotropic) resistivity is calculated by averaging the resistivity from all q points defined by nqf1*nqf2*nqf3. Thus,
1, The very fine q mesh (nqf1*nqf2*nqf3) gives us all positive values in all q points, why we don not need to consider the q points with negative values like (-0.5,0.5,0)?
2, If I want to calculate the resistivity along x axis (Suppose the lattice is orthogonal), I revised the a2f.90 in src
Code: Select all
DO iq = 1, nqtotf ! loop over q-points
Thank you in advance!
Jin Zhang