We calculate the P0 to gamma gamma, V0 to P0 gamma and V0to V'0 gamma gamma
decays in the framework of a U(3)xU(3) linear sigma model which includes
constituent quarks. For the first two decays this approach improves results
based on the anomalous Wess-Zumino term, with contributions due to SU(3)
symmetry breaking and vector mixing. The phi to (omega,rho) gamma gamma decays
are dominated by resonant eta' exchange . Our calculation for the later decays
improves and update similar calculations in the -closely related- framework of
vector meson dominance. We obtain BR(phi to rho gamma gamma)=2.5x10^{-5} and
BR(phi to omega gamma gamma)=2.8x10^{-6} within the scope of the
high-luminosity phi factories.Comment: 8 pages, submitted to Phys. Rev.