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:

[ERROR] STOP signal received while in :[07] Bethe-Salpeter Kernel

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

I guess I didn't calculate quasiparticle corrections for all the states for which they are needed.

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

The input for the GW calculations is the following:

gw0 # [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|

#%

whereas the input of the BSE calculation which is reading from the SAVE directory where the

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"