Sky surveys produce enormous quantities of data on extensive regions of the sky. The easiest way to access this information is through catalogues of standardised data products. {\em XMM-Newton} has been surveying the sky in the X-ray, ultra-violet, and optical bands for 20 years. The {\em XMM-Newton} Survey Science Centre has been producing standardised data products and catalogues to facilitate access to the serendipitous X-ray sky. Using improved calibration and enhanced software, we re-reduced all of the 14041 {\em XMM-Newton} X-ray observations, of which 11204 observations contained data with at least one detection and with these we created a new, high quality version of the {\em XMM-Newton} serendipitous source catalogue, 4XMM-DR9. 4XMM-DR9 contains 810795 detections down to a detection significance of 3 $\sigma$, of which 550124 are unique sources, which cover 1152 degrees$^{2}$ (2.85\%) of the sky. Filtering 4XMM-DR9 to retain only the cleanest sources with at least a 5 $\sigma$ detection significance leaves 433612 detections. Of these detections, 99.6\% have no pileup. Furthermore, 336 columns of information on each detection are provided, along with images. The quality of the source detection is shown to have improved significantly with respect to previous versions of the catalogues. Spectra and lightcurves are also made available for more than 288000 of the brightest sources (36\% of all detections).

The XMM-Newton Survey Science Centre Consortium (SSC) develops software in close collaboration with the Science Operations Centre to perform a pipeline analysis of all XMM-Newton observations. In celebration of the 20th launch anniversary, the SSC has compiled the 4th generation of serendipitous source catalogues, 4XMM. The catalogue described here, 4XMM-DR9s, explores sky areas that were observed more than once by XMM-Newton. It was constructed from simultaneous source detection on the overlapping observations, which were bundled in groups ("stacks"). Stacking leads to a higher sensitivity, resulting in newly discovered sources and better constrained source parameters, and unveils long-term brightness variations. As a novel feature, positional rectification was applied beforehand. Observations with all filters and suitable camera settings were included. Exposures with a high background were discarded, which was determined through a statistical analysis of all exposures in each instrument configuration. The X-ray background maps used in source detection were modelled via adaptive smoothing with newly determined parameters. Source fluxes were derived for all contributing observations, irrespective of whether the source would be detectable in an individual observation. From 1,329 stacks with 6,604 contributing observations over repeatedly covered 300 square degrees in the sky, 4XMM-DR9s lists 288,191 sources. 218,283 of them were observed several times. Most stacks are composed of two observations, the largest one comprises 352. The number of observations of a source ranges from 1 to 40. Auxiliary products like X-ray images, long-term light curves, and optical finding charts are published as well. 4XMM-DR9s is considered a prime resource to explore long-term variability of X-ray sources discovered by XMM-Newton. Regular incremental releases including new public observations are planned.

We aimed to identify the variable X-ray source 3XMM J000511.8+634018, which was serendipitously discovered through routine inspections while the 3XMM catalogue was compiled. We analysed the archival XMM-Newton observation of the source, obtained BUSCA photometry in three colours, and performed optical spectroscopy with the LBT. These data were supplemented by archival observations from the Zwicky Transient Facility. Based on its optical and X-ray properties, 3XMM J000511.8+634018 is classified as a magnetic cataclysmic variable, or polar. The flux is modulated with a period of 2.22 hr (8009.1 +- 0.2 s), which we identify with the orbital period. The bright phases are highly variable in X-ray luminosity from one cycle to the next. The source shows a thermal plasma spectrum typical of polars without evidence of a luminous soft blackbody-like component. It is non-eclipsing and displays one-pole accretion. The X-ray and BUSCA light curves show a stream absorption dip, which suggests an inclination 50degr < i < 75degr. The phasing of this feature, which occurs at the end of the bright phase, requires a somewhat special accretion geometry with a stream running far around the white dwarf before it is magnetically channelled. The period of this polar falls within the period gap of the cataclysmic variables (2.15-3.18 h), but appears to fall just below the minimum period when only polars are considered.

Aims: Spectral and temporal analysis of the NuSTAR observation Galactic Be-XRB Swift J1845.7-0037. during its recent outburst. Methods: For the spectral analysis we use both phenomenological and physics-based models. We employ an often used empirical model to identify the main characteristics of the spectral shape in relation to nominal spectral characteristics of X-ray pulsars. Additionally, we used the latest version of Bulk \& Thermal comptonization model (BW), to assess the validity of the spectral components required by the empirical model and to investigate the origin of the hard X-ray emission. We also analyzed the source light-curve, studying the pulse shape at different energy ranges and tracking the spectral evolution with pulse phase by using the model independent hardness ratio (HR). Results: We find that while both the empirical and physical (BW) spectral models can produce good spectral fits, the BW model returns physically plausible best-fit values for the source parameters and does not require any additional spectral components to the non-thermal, accretion column emission. The BW model also yielded an estimation of the neutron star magnetic field placing it in the 10^12G range. Conclusions: Our results, show that the spectral and temporal characteristics of the source emission are consistent with the scattering processes expected for radiation dominated shocks within the accretion column of highly magnetized accreting neutron stars. We further indicate that physically-derived spectral models such as BW, can be used to tentatively infer fundamental source parameters, in the absence of more direct observational signatures.

We have analyzed the X-ray spectra of all known Ultra Compact X-ray Binaries (UCXBs), with the purpose of constraining the chemical composition of their accretion disk and donor star. Our investigation was focused on the presence (or absence) of the Fe K${\alpha}$ emission line, which was used as the probe of chemical composition of the disk, based on previously established theoretical predictions for the reflection of X-ray radiation off the surface of C/O-rich or He-rich accretion disks in UCXBs. We have contrasted the results of our spectral analysis to the history of type I X-ray bursts from these systems, which can also indicate donor star composition. We found that UCXBs with prominent and persistent iron K${\alpha}$ emission also featured repeat bursting activity. On the other hand, the UCXBs for which no iron line was detected, appear to have few or no type I X-ray bursts detected over more than a decade of monitoring. Based on Monte Carlo simulations, demonstrating a strong correlation between the Fe K${\alpha}$ line strength and the abundance of C and O in the accretion disk material and given the expected correlation between the H/He abundance and the recurrence rate of type I X-ray bursts, we propose that there is a considerable likelihood that UCXBs with persistent iron emission have He-rich donors, while those that do not, likely have C/O or O/Ne/Mg-rich donors. Our results strongly advocate for the development of more sophisticated simulations of X-ray reflection from hydrogen-poor accretion disks.

We investigate the long-term variability of the K$\alpha$ line of iron in the spectra of two Ultra Compact X-ray Sources (UCXBs) with C/O-rich donors. We revisit archival observations by five different X-ray telescopes, over a ~twenty year period. Adopting physically motivated models for the spectral continuum, we probe the long-term evolution of the source emission in a self-consistent manner enabling physical interpretation of potential variability in the primary X-ray emission continuum and/or any emission lines from reflection off the accretion disk. We find that the spectral shape and flux of the source emission (for both objects) has remained almost constant throughout all the observations, displaying only minor variability in some spectral parameters and the source flux (largest variation is a ~25% drop in the flux of Swift J1756.9-2508). We note a striking variability of the Fe K$\alpha$ line which fluctuates from a notable equivalent width of ~66-100 eV in 4U 1543-624 and ~170 eV in Swift J1756.9-2508, to non-detections with upper limits of 2-8 eV. We argue that the disappearance of the iron line is due to the screening of the Fe K$\alpha$ line by the overabundant oxygen in the C/O-rich UCXBs. This effect is cancelled when oxygen becomes fully ionized in the inner disk region, resulting in the variability of the Fe K$\alpha$ line in an otherwise unaltered spectral shape. This finding supports earlier predictions on the consequences of H-poor, C/O-rich accretion disk on reflection induced fluorescent lines in the spectra of UCXBs.

In the current work we explore the applicability of standard theoretical models of accretion to the observed properties of M51 ULX-7. The spin-up rate and observed X-ray luminosity are evidence of a neutron star with a surface magnetic field of $2-7\times10^{13}$ G, rotating near equilibrium. Analysis of the X-ray light-curve of the system (Swift/XRT data) reveals the presence of a $\sim$39 d super-orbital period. We argue that the super-orbital periodicity is due to disc precession, and that material is accreted onto the neutron star at a constant rate throughout it. Moreover, by attributing this modulation to the free precession of the neutron star we estimate a surface magnetic field strength of $3-4\times10^{13}$ G. The agreement of these two independent estimates provide strong constraints on the surface polar magnetic field strength of the neutron star.

NGC300 ULX1 is an ultra-luminous X-ray pulsar, showing an unprecedented spin evolution, from about 126 s to less than 20 s in only 4 years, consistent with steady mass accretion rate. Following its discovery we have been monitoring the system with Swift/XRT and NICER to further study its properties. We found that even though the observed flux of the system dropped by a factor of $\gtrsim$20, the spin-up rate remained almost constant. A possible explanation is that the decrease in the observed flux is a result of increased absorption of obscuring material due to outflows or a precessing accretion disk.

In recent years, the discovery of massive quasars at z~7 has provided a striking challenge to our understanding of the origin and growth of supermassive black holes in the early Universe. Mounting observational and theoretical evidence indicates the viability of massive seeds, formed by the collapse of supermassive stars, as a progenitor model for such early, massive accreting black holes. Although considerable progress has been made in our theoretical understanding, many questions remain regarding how (and how often) such objects may form, how they live and die, and how next generation observatories may yield new insight into the origin of these primordial titans. This review focusses on our present understanding of this remarkable formation scenario, based on discussions held at the Monash Prato Centre from November 20--24, 2017, during the workshop "Titans of the Early Universe: The Origin of the First Supermassive Black Holes."

XMM-Newton has observed the X-ray sky since early 2000. The XMM-Newton Survey Science Centre Consortium has published catalogues of X-ray and ultraviolet sources found serendipitously in the individual observations. This series is now augmented by a catalogue dedicated to X-ray sources detected in spatially overlapping XMM-Newton observations. The aim of this catalogue is to explore repeatedly observed sky regions. It thus makes use of the long(er) effective exposure time per sky area and offers the opportunity to investigate long-term flux variability directly through the source detection process. A new standardised strategy for simultaneous source detection on multiple observations is introduced. It is coded as a new task within the XMM-Newton Science Analysis System and used to compile a catalogue of sources from 434 stacks comprising 1,789 overlapping XMM-Newton observations that entered the 3XMM-DR7 catalogue, have a low background and full-frame readout of all EPIC cameras. The first stacked catalogue is called 3XMM-DR7s. It contains 71,951 unique sources with positions and parameters such as fluxes, hardness ratios, quality estimates, and information on inter-observation variability. About 15% of the sources are new with respect to 3XMM-DR7. Through stacked source detection, the parameters of repeatedly observed sources can be determined with higher accuracy than in the individual observations. The method is more sensitive to faint sources and tends to produce fewer spurious detections. With this first stacked catalogue we demonstrate the feasibility and benefit of the approach. It supplements the large data base of XMM-Newton detections by additional, in particular faint, sources and adds variability information. In the future, the catalogue will be expanded to larger samples and continued within the series of serendipitous XMM-Newton source catalogues.

Methods. We analyzed broadband XMM-Newton and NuSTAR observations of NGC 300 ULX1, performing phase-averaged and phase-resolved spectroscopy. We compared two physically motivated models for the source spectrum: Non-thermal accretion column emission modeled by a power law with high-energy exponential roll-off (AC model) vs multicolor thermal emission from an optically thick accretion envelope plus a hard power-law tail (MCAE model). We combine the findings of our Bayesian analysis with qualitative physical considerations to evaluate the suitability of each model. Results. The low-energy part (<2 keV) of the source spectrum is dominated by non-pulsating, multicolor thermal emission. The (pulsating) high energy continuum is more ambiguous. If modelled with the AC model a residual structure is detected, that can be modeled using a broad Gaussian absorption line centered at ~12 keV. However, the same residuals can be successfully modeled using the MCAE model, without the need for the absorption-like feature. Model comparison, using the Bayesian approach strongly indicates that the MCAE model -- without the absorption line -- is the preferred model. Conclusions. All models considered strongly indicate the presence of an accretion disk truncated at a large distance from the central object, as has been recently suggested for a large fraction of both pulsating and non-pulsating ULXs. The hard, pulsed emission is not described by a smooth spectral continuum. If modelled by a broad Gaussian absorption line, the fit residuals can be interpreted as a cyclotron scattering feature (CRSF) compatible with a ~10^12 G magnetic field. However, the MCAE model describes the spectral and temporal characteristics of the source emission, without the need for an absorption feature and yields physically meaningful parameter values. Therefore strong doubts are cast on the presence of a CRSF in NGC 300 ULX1.

We study the temporal and spectral characteristics of SMC X-3 during its 2016 outburst. To probe accretion onto highly magnetized neutron stars (NSs), at the Eddington limit. We obtained XMM-Newton observations of SMC X-3 and combined them with long-term observations by Swift. We studied the temporal and spectral behavior of the source, and its short- and long-term evolution. We constructed a simple toy-model to gain insight into the complex emission pattern of SMC X-3. We confirm the pulse period of the system, derived by previous works and note that the pulse has a complex, three peaked shape. We find that the pulsed emission is dominated by hard photons, while at energies below ~1 keV, the emission is virtually non-pulsating. We further find that the shape of the pulse profile and the short and long-term evolution of the source light-curve can be explained by invoking a combination of a "fan" and a "polar" beam. Our spectroscopic analysis, reveals a two-component emission: a hard power law and a soft thermal component. We find that the latter produces the bulk of the non-pulsating emission and is most likely the result of reprocessing of the primary hard emission, by optically thick material that partially obscures the central source. We also detect strong emission lines from highly ionized metals. The strength of the emission lines are strongly phase-depended. The energy and temporal evolution and the shape of the pulse profile and the long-term spectra evolution of the source are consistent with the expected emission pattern of the accretion column in the super-critical regime, while the presence of a large reprocessing region is consistent with the analysis of previously studied X-ray pulsars observed at high accretion rates. The presence of this region is consistent with recently proposed works that suggest that high-B NSs occupy a considerable fraction of ULXs.

Intermediate mass black holes (IMBHs) are an (as yet) elusive class of black holes that are expected to lie in the $10^2-10^5\,M_{\odot}$ range, between the firmly established stellar-mass black holes and ${\gtrsim}10^6\,M_{\odot}$ supermassive black holes. Predicted by a variety of theoretical models, IMBHs are the potential seeds of supermassive black holes and are expected to power some of the brightest extra-nuclear X-ray sources. This brief review is the result of a presentation and subsequent discussion of IMBHs that took place during the 12th International Frascati workshop on "Multifrequency Behaviour of High Energy Cosmic Sources". The manuscript aims to provide a concise and up-to-date review of the different evolutionary scenarios for the creation of IMBHs. Throughout the text I emphasize the importance of the identification and classification of IMBHSs in our effort to understand the formation of supermassive black holes and their co-evolution with their host galaxies.

Aims: In light of recent discoveries of pulsating ULXs and recently introduced models, that propose neutron stars (NSs) as the central engines of ULXs, we revisit the (XMM-Newton and NuSTAR) spectra of eighteen well known ULXs, in search of indications that favor this new hypothesis. Results: We confirm the, previously noted, presence of the low energy (<6 keV) spectral rollover and argue that it could be interpreted as thermal emission. The spectra are well described by a double thermal model consisting of a "hot" (>1 keV) and a "cool" (<0.7 keV) multicolor black body (MCB). Under the assumption that the "cool" MCD emission originates in a disk truncated at the NS magnetosphere, we find that all ULXs in our sample are consistent with accretion onto a highly magnetized (B>10^12G) NS. We note a strong correlation between the strength of the magnetic field, the temperature of the "hot" thermal component and the total unabsorbed luminosity. Examination of the NuSTAR data supported this interpretation and confirmed the presence of a weak, high energy (>15 keV) tail, most likely the result of modification of the MCB emission by inverse Compton scattering. We also note that the apparent high energy tail, may simply be the result of mis-modeling of MCB emission with an atypical temperature (T) vs radius (r) gradient, using a standard MCD model with a fixed gradient of T~r^-0.75. Conclusions: We have offered a new physical interpretation for the dual-thermal spectra of ULXs. We find that the best-fit derived parameters, are in excellent agreement with recent predictions that favor super-critically accreting NSs as the engines of a large fraction of ULXs. Nevertheless, the considerable degeneracy between models and the lack of unequivocal evidence cannot rule out other equally plausible interpretations. Our results warrant further exploration of this new framework.

Aims: This work is the first stage of a campaign to search for IMBHs, in low luminosity AGN (LLAGN) and dwarf galaxies. An additional, and equally important, aim of this pilot study is to investigate the consistency between the predictions of several BH scaling relations and the fundamental plane of black hole activity (FP-BH). Methods: We use X-ray and radio luminosity relations in accreting BHs, along with the latest scaling relations between the mass of the central black hole (MBH) and the properties of its host spheroid, to predict MBH in seven LLAGN, that were previously reported to be in the IMBH regime. Namely, we use the recently re-evaluated MBH - Msph (Msph : spheroid absolute magnitude at 3.6 microm) scaling relation for spiral galaxies, the MBH - nsph (nsph : major axis Sersic index of the bulge) relation, the MBH - PA (PA: pitch angle) relation and the FP-BH for weakly accreting BHs, to independently estimate MBH in all seven galaxies. Results: We find that all LLAGN in our list have low-mass central black holes - with log MBH/Msol~6.5 on average - but they are, not IMBHs. All 4 methods used predicted consistent BH masses, in the 1sigma range. Furthermore, we report that, in contrast to previous classification, galaxy NGC 4470 is bulge-less, and cast doubts on the AGN classification of NGC 3507. Conclusions: We find that our latest, state-of-the-art techniques for bulge magnitude & Sersic index computations and the most recent updates of the MBH - Lsph, MBH - nsph , MBH - PA relations and the FP-BH produce consistent results in the low mass regime. We establish a multiple method approach for predicting BH masses in the regime where their spheres-of- grav. influence cannot be spatially resolved. Our approach mitigates against outliers from any one relation and provides a more robust average prediction. We will use our new method to revisit more IMBH candidates in LLAGN

Thanks to the large collecting area (3 x ~1500 cm$^2$ at 1.5 keV) and wide field of view (30' across in full field mode) of the X-ray cameras on board the European Space Agency X-ray observatory XMM-Newton, each individual pointing can result in the detection of hundreds of X-ray sources, most of which are newly discovered. Recently, many improvements in the XMM-Newton data reduction algorithms have been made. These include enhanced source characterisation and reduced spurious source detections, refined astrometric precision, greater net sensitivity and the extraction of spectra and time series for fainter sources, with better signal-to-noise. Further, almost 50\% more observations are in the public domain compared to 2XMMi-DR3, allowing the XMM-Newton Survey Science Centre (XMM-SSC) to produce a much larger and better quality X-ray source catalogue. The XMM-SSC has developed a pipeline to reduce the XMM-Newton data automatically and using improved calibration a new catalogue version has been produced from XMM-Newton data made public by 2013 Dec. 31 (13 years of data). Manual screening ensures the highest data quality. This catalogue is known as 3XMM. In the latest release, 3XMM-DR5, there are 565962 X-ray detections comprising 396910 unique X-ray sources. For the 133000 brightest sources, spectra and lightcurves are provided. For all detections, the positions on the sky, a measure of the quality of the detection, and an evaluation of the X-ray variability is provided, along with the fluxes and count rates in 7 X-ray energy bands, the total 0.2-12 keV band counts, and four hardness ratios. To identify the detections, a cross correlation with 228 catalogues is also provided for each X-ray detection. 3XMM-DR5 is the largest X-ray source catalogue ever produced. Thanks to the large array of data products, it is an excellent resource in which to find new and extreme objects.

We detect variability of the Fe Ka emission line in the spectrum of X-ray pulsar 4U 1626-67, correlated with changes in its luminosity and in the shape of its pulse profile. Analysis of archival Chandra and RXTE observations revealed the presence of an intrinsically narrow Fe Ka emission line in the spectrum obtained during the source's current high luminosity period. However, the line was not present during an XMM-Newton observation seven years earlier, when the source was ~three times fainter. The line is resolved by the high energy grating of Chandra at the 98% confidence level, and its small intrinsic width, sigma=36.4 (-/+ 11.3,15.3)eV, suggests reflection off an accretion disk at the radius R~7.5(-/+ 3.8,8.2)x10^8cm assuming a Keplerian disk, viewed at an inclination angle of 20deg. This value is consistent with the radius of the magnetosphere of the pulsar, suggesting that the line originates near the inner edge of a disk that is truncated by the magnetic field of the neutron star. Timing analysis of the XMM-Newton and RXTE data revealed a major change in the pulse profile of the source from a distinct double peaked shape during the high luminosity state when the line was present, to a much more complex multi-peak structure during the low luminosity state. We argue that the appearance of the line and the change in the shape of the pulse profile are correlated and are the result of a major change in the emission diagram of the accretion column, from a pencil-beam pattern at low luminosity, to a fan-beam pattern at high luminosity.

We search for the Fe K${\alpha}$ line in spectra of Ultra Compact X-ray Binaries (UCXBs). For this purpose we have analyzed XMM-Newton observations of five confirmed UCXBs. We find that the object 2S 0918-549 - whose optical spectrum bears tentative signatures of a C/O accretion disk - is devoid of any emission features in the 6-7 keV range, with an upper limit of less than 10 eV for the equivalent width (EW) of the iron line. 4U 1916-05 - whose optical spectrum is consistent with reflection from a He-rich accretion disk - exhibits a bright broad iron emission line. This behavior is in agreement with the theoretical predictions presented in Koliopanos, Gilfanov and Bildsten (2013). Namely, we expect strong suppression of the Fe K${\alpha}$ emission line in spectra originating in moderately bright (LogLx less than $\approx$ 37.5) UCXBs with C/O or O/Ne/Mg-rich donors. On the other hand the EW of the iron line in spectra from UCXBs with He-rich donors is expected to retain its nominal value of $\approx$ 100 eV. Our analysis also reveals a strong Fe K${\alpha}$ line in the spectrum of 4U 0614+091. This detection points towards a He-rich donor and seems to be at odds with the source's classification as C/O-rich. Nevertheless, a He-rich donor would explain the bursting activity reported for this system. Lastly, based on our theoretical predictions, we attribute the lack of a strong iron emission line - in the two remaining UCXB sources in our sample (XTE J1807-294, 4U 0513-40) - as an indication of a C/O or O/Ne/Mg white dwarf donor. From the upper limits of the Fe K${\alpha}$ line EW in 4U 0513-40, 2S 0918-549 and XTE J1807-294 we obtain a lower limit on the oxygen-to-iron ratio,O/Fe$\ge$ 10$\times$[O/Fe]$_{\odot}$.

Non-solar composition of the donor star in ultra-compact X-ray binaries may have a pronounced effect on the fluorescent lines appearing in their spectra due to reprocessing of primary radiation by the accretion disk and the white dwarf surface. We show that the most dramatic and easily observable consequence of the anomalous C/O abundance, is the significant, by more than an order of magnitude, attenuation of the Ka line of iron. It is caused by screening of the presence of iron by oxygen - in the C/O dominated material the main interaction process for a E ~ 7keV photon is absorption by oxygen rather than by iron, contrary to the solar composition case. Ionization of oxygen at high mass accretion rates adds a luminosity dependence to this behavior - the iron line is significantly suppressed only at low luminosity, log(LX) less than 37-37.5, and should recover its nominal strength at higher luminosity. The increase of the EW of the Ka lines of carbon and oxygen, on the other hand, saturates at rather moderate values. Screening by He is less important, due to its low ionization threshold and because in the accretion disk it is mostly ionized. Consequently, in the case of the He-rich donor, the iron line strength remains close to its nominal value, determined by the iron abundance in the accretion disk. This opens the possibility of constraining the nature of donor stars in UCXBs by means of X-ray spectroscopy with moderate energy resolution.