Well, the log file says MemoryError, so yeah, that's the problem.
You have a lot of k-points and a lot of atoms. Try setting
save_partial = False
inside DensityOfStates call
Note that you do not need to rerun your entire script. Just read the converged calculation from the NC file and just run
bulk_configuration = nlread('GaN25-Al75.nc', BulkConfiguration)[-1]
# -------------------------------------------------------------
# Density Of States
# -------------------------------------------------------------
kpoint_grid = MonkhorstPackGrid(
na=15,
nb=15,
nc=15,
)
density_of_states = DensityOfStates(
configuration=bulk_configuration,
kpoints=kpoint_grid,
energy_zero_parameter=FermiLevel,
bands_above_fermi_level=None,
save_partial = False,
)
nlsave('GaN25-Al75.nc', density_of_states)
nlprint(density_of_states)