Dear developer:
I am doing a calculation with both extended_kernel and exciton_Q_shift turned on. I found some indications that BerkeleyGW may not did it correctly.
When comparing besmat.h5 generated for Q=0 and a finite Q1. Something strange happened when I check the results in bsemat.h5 ,specifically the kernel matrix.
I found when the four band indices are all conduction bands. The head matrix ,wing matrix and body matrix elements of Q=0 calculation are exactly the same as those of a Q=Q1 calculation. Otherwise, those matrix are different as expected(when at least one band indices is valance bands, or for exchange matrix)
I suspect, when writing the code, finite Q shift is only applied to valance band state. Which works fine without exciton_Q_shift but fails in finite Q calculation .
I wonder whether BerkeleyGW correctly treat extended kernel with finite Q, especially the extended kernel involving four conduction bands states