We have carried out a study of the X-ray properties of the supernova remnant (SNR) population in M33 with XMM-Newton, comprising deep observations of 8 fields in M33 covering all of the area within the D$_{25}$ contours, and with a typical luminosity of 7.1$\times$10$^{34}$ erg s$^{-1}$ (0.2-2.0 keV) . Here we report our work to characterize the X-ray properties of the previously identified SNRs in M33, as well as our search for new X-ray detected SNRs. With our deep observations and large field of view we have detected 105 SNRs at the 3$\sigma$ level, of which 54 SNRs are newly detected in X-rays, and three are newly discovered SNRs. Combining XMM-Newton data with deep Chandra survey data allows detailed spectral fitting of 15 SNRs, for which we have measured temperatures, ionization timescales, and individual abundances. This large sample of SNRs allows us to construct an X-ray luminosity function, and compare its shape to luminosity functions from host galaxies of differing metallicities and star formation rates to look for environmental effects on SNR properties. We conclude that while metallicity may play a role in SNR population characteristics, differing star formation histories on short timescales, and small-scale environmental effects appear to cause more significant differences between X-ray luminosity distributions. In addition, we analyze the X-ray detectability of SNRs, and find that in M33 SNRs with higher [SII]/H$\alpha$ ratios, as well as those with smaller galactocentric distances, are more detectable in X-rays.

We have obtained a deep 8-field XMM-Newton mosaic of M33 covering the galaxy out to the D$_{25}$ isophote and beyond to a limiting 0.2--4.5 keV unabsorbed flux of 5$\times$10$^{-16}$ erg cm$^{-2}$ s$^{-1}$ (L${>}$4$\times$10$^{34}$ erg s$^{-1}$ at the distance of M33). These data allow complete coverage of the galaxy with high sensitivity to soft sources such as diffuse hot gas and supernova remnants. Here we describe the methods we used to identify and characterize 1296 point sources in the 8 fields. We compare our resulting source catalog to the literature, note variable sources, construct hardness ratios, classify soft sources, analyze the source density profile, and measure the X-ray luminosity function. As a result of the large effective area of XMM-Newton below 1 keV, the survey contains many new soft X-ray sources. The radial source density profile and X-ray luminosity function for the sources suggests that only $\sim$15% of the 391 bright sources with L${>}$3.6$\times$10$^{35}$ erg s$^{-1}$ are likely to be associated with M33, and more than a third of these are known supernova remnants. The log(N)--log(S) distribution, when corrected for background contamination, is a relatively flat power-law with a differential index of 1.5, which suggests many of the other M33 sources may be high-mass X-ray binaries. Finally, we note the discovery of an interesting new transient X-ray source, which we are unable to classify.