In Yambo’s RPA/BSE calculations, I want to use spin-, k- and band-resolved QP (qusi-particle) corrections obtained using another GW code (WEST code). Could anyone tell me how to do that?

Personally, I tried to modify the PWSCF output files eigenval1(2).xml to add state-resolved QP corrections onto the stored eigen energies. However, I found that with the modified eigenval1(2).xml files, the calculated RPA absorption spectra by Yambo has right peak positions but wrong (low) peak heights or amplitudes.

Please find the attached figure to see the problem. In the figure, I compared three Yambo RPA spectra: (i) grey - no QP corrections; (ii) red - QP corrections set by KfnQP_up_E and KfnQP_dn_E; (iii) blue – QP corrections are added to eigenval1(2).xml files of PWSCF.

The studied system is a supercell of h-BN with a defect. The first peak in the spectra corresponds to a transition between two defect states.

Could anyone help to explain why adding QP corrections to PWSCF output files eigenval1(2).xml does not work? Is this trouble related to a kind of sum rule?

I also attached one RPA input file using variables KfnQP_up_E and KfnQP_dn_E.

- Code: Select all
`rim_cut # [R RIM CUT] Coulomb potential`

optics # [R OPT] Optics

bse # [R BSE] Bethe Salpeter Equation.

bsk # [R BSK] Bethe Salpeter Equation kernel

bss # [R BSS] Bethe Salpeter Equation solver

RandQpts= 100000 # [RIM] Number of random q-points in the BZ

RandGvec= 10 RL # [RIM] Coulomb interaction RS components

CUTGeo= "box z" # [CUT] Coulomb Cutoff geometry: box/cylinder/sphere X/Y/Z/XY..

% CUTBox

0.0 | 0.0 | 19.0 | # [CUT] [au] Box sides

%

CUTRadius= 0.000000 # [CUT] [au] Sphere/Cylinder radius

CUTCylLen= 0.000000 # [CUT] [au] Cylinder length

BSEmod= "retarded" # [BSE] resonant/retarded/coupling

BSKmod= "Hartree" # [BSE] IP/Hartree/HF/ALDA/SEX

BSSmod= "d" # [BSS] (h)aydock/(d)iagonalization/(i)nversion/(t)ddft`

BSENGBlk= 1 Ry # [BSK] Screened interaction block size

% KfnQP_up_E

3.200000 | 1.000000 | 1.000000 | # [EXTQP BSK BSS] E parameters UP (c/v) eV|adim|adim

%

% KfnQP_dn_E

3.800000 | 1.000000 | 1.000000 | # [EXTQP BSK BSS] E parameters UP (c/v) eV|adim|adim

%

WRbsWF

% BEnRange

0.00000 | 9.00000 | eV # [BSS] Energy range

%

% BDmRange

0.02000 | 0.02000 | eV # [BSS] Damping range

%

BEnSteps= 9001 # [BSS] Energy steps

% BLongDir

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

%

% BSEBands

131 | 162 | # [BSK] Bands range

%

Junqing