Your pulse width fall in an uncomfortable range, for
analysis. It's not short enough to be "adiabatic" (though
this analysis produces the most uncomfortable
numbers) yet it's not wide enough to be fully settled
against the likely thermal time constant of the mass.
A 1/1000 duty cycle would likely make the DC equiv
dissipation look fine.
Questions include "where in the device, is the dissipation?"
(collector is a good answer because it's large, the dissipation
likely is diffuse and it's well bonded to the paddle for a high
heat capacity; a small region surrounding the base, might
put that 0.5W*100us=50uJ into a much smaller volume
that doesn't respond well to (say) time at 800C. Not that
this is the case, necessarily, but you'd need to know device
geometry and intra-device power dissipation distribution to
say what is.
A rule of thumb I have used, and gotten away with, is if
peak power is less than 1/10 the DC power / current rating,
and duty cycle averaged power / current is less than DC
rating, you're OK.
But maybe I just haven't found the challenge that breaks
that thumb, yet.