On scale functions for Lévy processes with negative phase-type jumps