The propagators of massless Bosons are only zero inside the lightcone if the number of spatial dimensions is odd and larger than 1.
You can determine the exact form of the propagator on the light-cone, even though it is divergent, and do so for any number of spatial dimensions, in the following way:
Use this expression to derive the propagator in d spatial dimensions from the 1 dimensional case with,
$P_d{(t,r)}\ =\ \frac{1}{2\pi^a}\ \frac{\partial^a
}{\partial (s^2)^a}\ P_1{(t,r)} $
where $~P{(t,r)}~$ is the propagator, with $~~a=(d-1)/2~$ and $~~s^2=t^2-r^2$
For the real valued Klein Gordon propagator this becomes:
$P_d^{KG}{(t,r)}\ =\ \frac{1}{2\pi^a}\
\frac{\partial^a }{\partial (s^2)^a} \left\{\ {\cal H}(s^2) J_o(ms)\
\right\}$
Where ${\cal H}(s^2)$ is the Heaviside step function which is 1 inside the light-cone and 0 outside the light-cone. In the mass-less case this becomes:
$P_d^{KG}{(t,r)}\ =\ \frac{1}{2\pi^a}\
\frac{\partial^a }{\partial (s^2)^a} \left\{\ {\cal H}(s^2)\
\right\}$
The following graphs showing the mass-less cases are from numerical simulations:
In the case with 1 spatial dimension you see the Heaviside step function. The 3d case is the first order derivative of the step function. The 5d case is the 2nd order derivative of the step-function and so on.
The even dimension cases are non-zero inside the light-cone because of the 1/2 order derivatives.
The operator which derives the d-dimensional propagator from the 1-dimensional propagator is derived in my paper here: http://www.physics-quest.org/Higher_dimensional_EM_radiation.pdf in section V.