finite Q exciton

[Bill Ash]

Dear BGW developers, I want to compute the exciton band structures. As shown in the figure below. 

So far I have added the flag of exciton_Q_shift in kernel.inp and absorption.inp and computed it, may I know in which file I should find the result?

Here is my input file：

epsilon.inp

epsilon_cutoff           8.0
number_bands             178
cell_slab_truncation
frequency_dependence 0
degeneracy_check_override
begin qpoints
0.001000000 0.000000000 0.000000000 1.0 1
0.000000000 0.333333333 0.000000000 1.0 0
0.000000000 0.666666667 0.000000000 1.0 0
0.333333333 0.000000000 0.000000000 1.0 0
0.333333333 0.333333333 0.000000000 1.0 0
0.333333333 0.666666667 0.000000000 1.0 0
0.666666667 0.000000000 0.000000000 1.0 0
0.666666667 0.333333333 0.000000000 1.0 0
0.666666667 0.666666667 0.000000000 1.0 0
end

sigma.inp

screened_coulomb_cutoff  8.0
number_bands            178
band_index_min   35
band_index_max   44
screening_semiconductor
degeneracy_check_override
dont_use_vxcdat
cell_slab_truncation
exact_static_ch 1
begin kpoints
0.000000000 0.000000000 0.000000000 1.0
0.000000000 0.333333333 0.000000000 1.0
0.000000000 0.666666667 0.000000000 1.0
0.333333333 0.000000000 0.000000000 1.0
0.333333333 0.333333333 0.000000000 1.0
0.333333333 0.666666667 0.000000000 1.0
0.666666667 0.000000000 0.000000000 1.0
0.666666667 0.333333333 0.000000000 1.0
0.666666667 0.666666667 0.000000000 1.0
end

kernel.inp

number_val_bands 4
number_cond_bands 4
cell_slab_truncation
screening_semiconductor
exciton_Q_shift 1 0.333333333 0.333333333 0.000000000

absorption.inp

diagonalization
number_val_bands_coarse 4
number_val_bands_fine 2
number_cond_bands_coarse 4
number_cond_bands_fine 2
no_symmetries_coarse_grid
no_symmetries_fine_grid
no_symmetries_shifted_grid
screening_semiconductor
cell_slab_truncation
use_velocity
gaussian_broadening
energy_resolution 0.02
degeneracy_check_override
eqp_co_corrections
eqp_co_q_corrections
write_eigenvectors 10
exciton_Q_shift 1 0.333333333 0.333333333 0.000000000