  Using Bogoliubov-de Gennes (BdG) equations we numerically calculate the
disorder averaged density of states of disordered semiconductor nanowires
driven into a putative topological p-wave superconducting phase by spin-orbit
coupling, Zeeman spin splitting and s-wave superconducting proximity effect
induced by a nearby superconductor. Comparing with the corresponding
theoretical self-consistent Born approximation (SCBA) results treating disorder
effects, we comment on the topological phase diagram of the system in the
presence of increasing disorder. Although disorder strongly suppresses the
zero-bias peak (ZBP) associated with the Majorana zero mode, we find some clear
remnant of a ZBP even when the topological gap has essentially vanished in the
SCBA theory because of disorder. We explicitly compare effects of disorder on
the numerical density of states in the topological and trivial phases.
