PAllen said:
But is there some reasonably simple reason it doesn't happen according to QED (which is a theory in flat Minkowski space)?
I would find it extremely unsatisfying to imagine that the only answer in the context of QED was something technical about field theory. After all, QED is only one tiny piece of the standard model.
Classical kinematics says that the lifetime has to go like \tau=(\hbar/c^4m^2)E, where m is some constant with units of mass. Since we observe quasars at cosmological distances at wavelengths on the order of 10 m, m \lesssim 10^{-56} kg. For a point source, the squared wavefunction of the photon would be multiplied by an exponential factor e^{-r/c\tau}.
But this would be exactly equivalent to the result of the Proca equation for a massive, spin-1 field, which has an exponential factor of e^{-2Mcr/\hbar} for the squared field, where M is the mass of the particle.
There have been various experiments that have tested Maxwell's equations to high precision. The exponent in Coulomb's law has been constrained to be within 10^-16 of 2. The mass of the photon has been constrained to be less than about 10^{-51} kg. These experiments all use very different techniques, but suppose that you're doing such an experiment, and your technique isn't sensitive to the decay products of a \gamma\rightarrow 3\gamma decay, only to the probability of receiving the parent photon. Then you may consider yourself to be putting an upper limit on the photon mass M, but this is exactly equivalent to putting an upper limit on 2m^2/E (in units with c=1), where m is the constant with units of mass that characterizes the photon's decay rate.
Now a massive photon breaks gauge invariance. QED can accommodate it, but it would be a horrible problem for QFT:
http://optica.machorro.net/Lecturas/PhotonMass_rpp5_1_RO2.pdf This seems to me like relatively nontechnical fact about QFT, too. So this suggests to me that if there is a relatively nontechnical way of prohibiting a nonzero value of M, then there might also be a relatively nontechnical way of prohibiting a nonzero m.
One interesting thing about the relationship between M and m is that if M has any nonzero value, no matter how ridiculously small, then \gamma\rightarrow 3\gamma decay is prohibited. That is, although a nonzero value of M and a nonzero value of m would act similarly in certain experiments, a nonzero M actually forces m=0.
If \gamma\rightarrow 3\gamma was possible, then the lifetimes of the products would be, on average, 1/3 the lifetime \tau_o of the parent. You would have an accelerating process of decay, and it would accelerate geometrically. Summing a geometric series, basically you expect that within a few times \tau_o, the decay would run to completion, meaning that you would produce a perfectly focused jet consisting of infinitely many photons, each with infinitesimal energy. Because of their infinitesimal energies, they would have to have infinite wavelengths. Their fields would add incoherently, so the average field would be infinitesimal, and therefore undetectable. The only way they would be detectable would be through their gravitational fields.
For the particle-in-a-box version, with the box being at a low enough temperature, the result would be that within a few times \tau_o, any photon that you initially put in the box would end up as a Bose-Einstein condensate.