hlee wrote: ↑Wed Feb 08, 2023 4:33 pm

Dear Chendix:

.., I'm replacing the hr.dat file and .eig file with my modified Hamiltonians and eigenvalues. However, after I fed the modified hr.dat and .eig files into EPW ...

How did you feed them into EPW?

The correct way is to (1) perform scf, nscf, phonon calculations with external fields, (2) starting from the obtained results to carry out EPW calculations.

It is not the correct way just to use the Hamiltonian and the eigenvalues under external fields without going through the steps above.

Sincerely,

H. Lee

Hi Dr. Lee,

Thank you for the response.

Actually the external fields are regarded as ultrafast pump fields. Considering the retardedness of phonon in ultrafast timescale, we assume there's no change on relaxation, scf, nscf and phonon. Based on this idea, we ran the Wannierization step of EPW to generate .ukk, .mmn and .bvec files and then we replaced the Hamiltonian and eigenvalues by our modified version on el-ph coupling step to calculate lambda and a2f. (Please correct me if I'm wrong, but I don't think it'll read H_R, since H_R (variable chw(nbndsub, nbndsub, nrr)) is directly delivered from subroutine 'hambloch2wan' to 'hamwan2bloch')

However, with the lack of correct eigenvectors, the band structure shows great difference. To solve this problem once and for all, I'm now learning the EPW source code to find the real space vectors required for my calculation, which can be derived from lattice and orbital indices (first five columns in hr.dat). I'm wondering if you are familiar with subroutine 'hambloch2wan' in bloch2wan.f90, specifically for the variable wslen(nrr), which represents Wigner-Seitz vector length?

Thank you again.

Best regards,

Chendix