We use an n-dimensional projection method to calculate helicity amplitudes for photon-photon scattering up to the two-loop level, for massless particles inside the loop. We confirm the recent standard model results found by Bern et al using another calculational method. We further apply the method to N = 1 and N = 2 supersymmetric QED. As predicted by the supersymmetric Ward identity, two of the independent amplitudes vanish. The remaining amplitude simplifies as the number of supersymmetries increases. For N = 2 SUSY QED, only contributions of weight 4 remain.