We review the recent progress in the calculation of the amplitude for γ γ -> π0 π0 to two loops in chiral perturbation theory. We match the low-energy amplitude in chiral perturbation theory with the result of the dispersion theoretic analysis. The neutral pion polarizabilities are also given to two-loop accuracy. Then, the results are compared with the dispersion relation calculation of the pion polarizabilities