Your system is periodic in the directions B and C (Y and Z) so please increase ky and kz in SCF and DOS calculation
like
numerical_accuracy_parameters = NumericalAccuracyParameters(
k_point_sampling=(1, 27, 27),
)
(I am not sure k value 27 is the best one)
and
density_of_states = DensityOfStates(
configuration=bulk_configuration,
kpoints=MonkhorstPackGrid(1,101,101),
energy_zero_parameter=FermiLevel,
bands_above_fermi_level=None,
)
Usually we need lager values of ks in DOS calculation than SCF's ones.