evGW
Posted: Sun Feb 08, 2026 9:28 am
Dear all,
When I was running a G1W0 calculation, I copied into the G1W0 directory all the files generated during my G0W0 run (ndb.HF_and_locXC, ndb.dipoles, ndb.pp*, ndb.cutoff, ndb.RIM, ndb.RIM_W). In the G1W0 input, I used the parameter `GfnQPdb` to read the quasiparticle energies produced by the G0W0 calculation, i.e., `GfnQPdb = "E < ./360b_15Ry/ndb.QP"`.
However, the results I obtained are very strange: compared with G0W0, the G1W0 result not only fails to open the band gap, but instead closes it.
I would like to know whether this is due to an operational mistake in my self-consistent evGW0 procedure, or because the number of bands stored in `ndb.QP` from G0W0 and used for G1W0 is too small (only a dozen or so bands around the Fermi level, though this is sufficient for BSE calculations).
Below I attach the G1W0 input file and the corresponding .r file. Thank you in advance for your help!
When I was running a G1W0 calculation, I copied into the G1W0 directory all the files generated during my G0W0 run (ndb.HF_and_locXC, ndb.dipoles, ndb.pp*, ndb.cutoff, ndb.RIM, ndb.RIM_W). In the G1W0 input, I used the parameter `GfnQPdb` to read the quasiparticle energies produced by the G0W0 calculation, i.e., `GfnQPdb = "E < ./360b_15Ry/ndb.QP"`.
However, the results I obtained are very strange: compared with G0W0, the G1W0 result not only fails to open the band gap, but instead closes it.
I would like to know whether this is due to an operational mistake in my self-consistent evGW0 procedure, or because the number of bands stored in `ndb.QP` from G0W0 and used for G1W0 is too small (only a dozen or so bands around the Fermi level, though this is sufficient for BSE calculations).
Below I attach the G1W0 input file and the corresponding .r file. Thank you in advance for your help!