I hope you are doing well. Recently, I am running BSE calculation of absorption spectrum of 2D MoS2. The DFT cal is done in the QE (lattice vector a=3.168 A, c=20 A, ecut=60 Ry, nbands=500, k-grid=27*27*1). Then, I have G0W0 and G1W1 calcultions (see the gw_try1.in in the attachment). Later, I followed the tutorial listed in the webcite of Yambo tutorial and three input files are attached. My problem is that I can not reproduce the A or B peak in the absorption spectrum at around 2.0 eV(I know I did not include spin-orbit coupling, there should be split of peaks at 2.0 eV, but there should be a single peak at around 2.0 eV. ). I can reproduce the absorption spectrum at around 3.0 eV. The convergence test of paprameters is the same as the nature paper mentioned. The absortion spectrum is below

Xin