Reading electron-phonon matrix elements file by EPW
Posted: Thu Feb 20, 2025 6:27 pm
Hi,
I am trying to read the ephmat files generated by EPW in binary format using Python. I was reading the source code write_ephmat routine within the io_eliashberg.f90 file. I am doing for GaN with 1 band and there are 12 phonon modes. The dense k and q grid is 40^3. I am using 4 nodes and 16 cores so there are 64 processes. Thus I got 64 ephmat files with 1000 kpts in each of them. I wanted to verify the format of the ephmat file generated by EPW. I am using EPW 5.7 and QE7.2. I am also using a large fsthick to ensure all the k and q points are taken into the calculation
So the first line contains the pool id and the number of k points within the pool right? And after that the matrix elements are written. As far as I understood, the first matrix element is the one corresponding to the 1st kpoint 1st qpoint , 1st band and the 1st mode[/b] right? Then the next elements is 1st kpoint, 1st qpoint,1st band and 2nd mode right? In this way it goes upto 12 modes and then is it is for 1st kpt, 2nd qpt, 1st band and 1st modes right?
It would be really great if someone could verify whether I am understanding it correctly.
I am trying to read the ephmat files generated by EPW in binary format using Python. I was reading the source code write_ephmat routine within the io_eliashberg.f90 file. I am doing for GaN with 1 band and there are 12 phonon modes. The dense k and q grid is 40^3. I am using 4 nodes and 16 cores so there are 64 processes. Thus I got 64 ephmat files with 1000 kpts in each of them. I wanted to verify the format of the ephmat file generated by EPW. I am using EPW 5.7 and QE7.2. I am also using a large fsthick to ensure all the k and q points are taken into the calculation
So the first line contains the pool id and the number of k points within the pool right? And after that the matrix elements are written. As far as I understood, the first matrix element is the one corresponding to the 1st kpoint 1st qpoint , 1st band and the 1st mode[/b] right? Then the next elements is 1st kpoint, 1st qpoint,1st band and 2nd mode right? In this way it goes upto 12 modes and then is it is for 1st kpt, 2nd qpt, 1st band and 1st modes right?
It would be really great if someone could verify whether I am understanding it correctly.