my issue is in running a full GW+BSE calculation. I want the BSE calculation to reade the quasiparticle corrections

computed within GW approximation, to check the validity of the "scissor operator" shift of the BSE spectrum for

a system under study.

The error I get is as follows:

I guess I didn't calculate quasiparticle corrections for all the states for which they are needed.[ERROR] STOP signal received while in :[07] Bethe-Salpeter Kernel

[ERROR] Null BSE kernel block dimension(s) found. Increase e/h range

Can you give me some advice on how to perform this kind of test?

The input for the GW calculations is the following:

whereas the input of the BSE calculation which is reading from the SAVE directory where thegw0 # [R GW] GoWo Quasiparticle energy levels

ppa # [R Xp] Plasmon Pole Approximation

HF_and_locXC # [R XX] Hartree-Fock Self-energy and Vxc

em1d # [R Xd] Dynamical Inverse Dielectric Matrix

EXXRLvcs= 2999 RL # [XX] Exchange RL components

% QpntsRXp

1 | 36 | # [Xp] Transferred momenta

%

% BndsRnXp

1 | 348 | # [Xp] Polarization function bands

%

NGsBlkXp= 250 RL # [Xp] Response block size

% LongDrXp

1 | 0 | 0 | # [Xp] [cc] Electric Field

%

PPAPntXp= 27.21138 eV # [Xp] PPA imaginary energy

% GbndRnge

1 | 348 | # [GW] G[W] bands range

%

GDamping= 0.1 eV # [GW] G[W] damping

dScStep= 0.10000 eV # [GW] Energy step to evalute Z factors

DysSolver= "n" # [GW] Dyson Equation solver (`n`,`s`,`g`)

%QPkrange # [GW] QP generalized Kpoint/Band indices

1 | 36 | 160 | 189 |

%

#%QPerange # [GW] QP generalized Kpoint/Energy indices

# 1| 36| 0.0|-1.0|

#%

GW run wrote is:

optics # [R OPT] Optics

bse # [R BSK] Bethe Salpeter Equation.

em1s # [R Xs] Static Inverse Dielectric Matrix

bss # [R BSS] Bethe Salpeter Equation solver

em1d # [R Xd] Dynamical Inverse Dielectric Matrix

ppa # [R Xp] Plasmon Pole Approximation

BSresKmod= "xc" # [BSK] Resonant Kernel mode. (`x`;`c`;`d`)

BScplKmod= "none" # [BSK] Coupling Kernel mode. (`x`;`c`;`d`;`u`)

% BSEBands

160 | 189 | # [BSK] Bands range

%

BSENGBlk= 250 RL # [BSK] Screened interaction block size

BSENGexx= 2999 RL # [BSK] Exchange components

% QpntsRXs

1 | 36 | # [Xs] Transferred momenta

%

% BndsRnXs

1 | 348 | # [Xs] Polarization function bands

%

NGsBlkXs= 250 RL # [Xs] Response block size

% LongDrXs

1.000000 | 0.000000 | 0.000000 | # [Xs] [cc] Electric Field

%

% QpntsRXp

1 | 36 | # [Xp] Transferred momenta

%

% BndsRnXp

1 | 348 | # [Xp] Polarization function bands

%

NGsBlkXp= 253 RL # [Xp] Response block size

% LongDrXp

1.000000 | 0.000000 | 0.000000 | # [Xp] [cc] Electric Field

%

PPAPntXp= 27.21138 eV # [Xp] PPA imaginary energy

BSSmod= "d" # [BSS] Solvers `h/d/i/t`

% BEnRange

0.00000 | 8.00000 | eV # [BSS] Energy range

%

% BDmRange

0.10000 | 0.10000 | eV # [BSS] Damping range

%

BEnSteps= 500 # [BSS] Energy steps

% BLongDir

1.000000 | 0.000000 | 0.000000 | # [BSS] [cc] Electric Field

%

KfnQPdb= "E < SAVE/ndb.QP"