Sincere apologies Unberto, some how the Density Mesh Cutoff was not set to 425 Hartree. So I was still using the default 75 Hartree which is way too small for NiGe in HGHbasis. With 425 Hartree, PDOS on Ni are resembling to the reference. But PDOS on Ge are the same as in the image. I ll share the script and the results, need some time for that.