Search

Previous X-ray works on Mrk 1239 have revealed a complex Narrow Line Seyfert 1 (NLS1) that exhibits substantial absorption and strong emission from both collisional (CIE) and photoionized (PIE) plasmas. Here, we report on deep-pointed observations with $XMM{\rm -}Newton$ and $NuSTAR$, along with $Swift$ monitoring, to understand the $0.3-30$ keV continuum emission and the central engine geometry. A strong X-ray flare, where the AGN brightens by a factor of five in $\sim30$ ks, is captured between $4-30$ keV and can be attributed to a brightening of the primary continuum. However, the lack of any variability below $\sim3$ keV on long- or short-time scales requires complete absorption of the AGN continuum with a neutral medium of column density $\sim 10^{23.5}{\rm cm}^{-2}$. The timing and spectral properties are consistent with a blurred reflection interpretation for the primary emission. The variability and presence of a Compton hump disfavours ionized partial covering. The neutral absorber, if outflowing, could be crashing into the surrounding medium and ISM to produce the low-energy continuum and CIE. Scattered emission off the inner torus could produce the PIE. The intricate scenario is demanded by the data and highlights the complexity of the environment that is normally invisible when overwhelmed by the AGN continuum. Objects like Mrk 1239 serve as important sources for unveiling the interface between the AGN and host galaxy environments., Comment: 26 pages, accepted for publication in ApJ

We present an analysis of the HST COS spectrum of IZw1 aiming to probe the absorbing medium associated with the active galactic nucleus (AGN). We fitted the emission spectrum and performed spectral analysis of the identified absorption features to derive the corresponding ionic column densities and covering fractions of the associated outflows. We employed photoionisation modelling to constrain the total column density and the ionisation parameter of four detected kinematic components. By investigating the implications of the results together with the observed kinematic properties of both emission and absorption features, we derived constraints on the structure and geometry of the absorbing medium in the AGN environment. We find and characterise absorption line systems from outflowing ionised gas in four distinct kinematic components, located at -60, -280, -1950, and -2900 km/s with respect to the source rest frame. While the two slower outflows are consistent with a full covering of the underlying radiation source, the well-constrained doublet line ratios of the faster two, higher column density, outflows suggest partial covering, with a covering fraction of C_f~0.4. The faster outflows show also line-locking in the NV doublet, a signature of acceleration via line absorption. This makes IZw1 possibly the closest object that shows evidence for hosting line-driven winds. The observed -1950 km/s absorption is likely due to the same gas as an X-ray warm absorber. Furthermore, the behaviour in UV and X-ray bands implies that this outflow has a clumpy structure. We find that the highly asymmetric broad emission lines in IZw1, indicative of a collimated, outflowing broad line region, are covered by the absorbing gas. Finally, the strongest UV--X-ray absorber may be connected to some of the blueshifted line emission, indicative of a more spatially extended structure of this ionised medium., Comment: Accepted for publication in A&A

We discuss applications of the study of the new and barely explored class of changing-look (CL) narrow-line Seyfert 1 (NLS1) galaxies and comment on their detection with the space mission SVOM (Space Variable Objects Monitor). We highlight the case of NGC 1566, which is outstanding in many respects, for instance as one of the nearest known CL AGN undergoing exceptional outbursts. Its NLS1 nature is discussed, and we take it as a nearby prototype for systems that could be discovered and studied in the near future, including with SVOM. Finally, we briefly examine the broader implications and applications of CL events in NLS1 galaxies and show that such systems, once discovered in larger numbers, will greatly advance our understanding of the physics of the environment of rapidly growing supermassive black holes. This White Paper is part of a sequence of publications which explore aspects of our understanding of (CL) NLS1 galaxy physics with future missions., Comment: 14 pages, 4 figures. Accepted for publication in the Universe Special Issue "A Multimessenger View of Supermassive Black Holes and the Quasar Main Sequence"

In this work, we examine the application of the wavelet transform to the X-ray timing analyses of active galactic nuclei (AGN) and quasi-periodic eruption sources (QPEs). Several scenarios are simulated to test the effectiveness of the wavelet analysis to stationary and non-stationary data. We find that the power spectral density (PSD) slope and the nature of the periodic signal can influence the ability to identify important features in the wavelet power spectrum. In general, weak and transient features can be discerned, which make the wavelet spectrum an important tool in examining AGN light curves. We carried out a wavelet analysis to four unique objects: Ark 120, IRAS 13224-3809, RE J1034+396, and the QPE GSN 069. The well-known quasi-periodic oscillation (QPO) in RE J1034+396 is significantly detected in the wavelet power spectrum. In IRAS 13224-3809, significant transient features appear during a flare at frequencies coincident with previously detected reverberation signals. Finally, the wavelet power spectrum of the QPE GSN 069 significantly reveals four persistent signals that exhibit a 3:2 ratio in oscillation frequencies, consistent with high-frequency QPOs in stellar mass X-ray binaries, but we cannot rule out the possibility this is an artefact of the calculation., Comment: 14 pages, 14 figures, Accepted for publication to MNRAS

Based on our dedicated Swift monitoring program, MOMO, OJ 287 is one of the best-monitored blazars in the X-ray--UV--optical regime. Here, we report results from our accompanying, dense, multi-frequency (1.4--44 GHz) radio monitoring of OJ 287 between 2015 and 2022 covering a broad range of activity states. Fermi gamma-ray observations are added. We characterize the radio flux and spectral variability in detail, including DCF and other variability analyses, and discuss its connection with the multiwavelength emission. Deep fades of radio and optical--UV fluxes are found to occur every 1--2 years. Further, it is shown that a precursor flare of thermal bremsstrahlung predicted by one of the binary supermassive black hole (SMBH) models of OJ 287 was absent. We then focus on the nature of the extraordinary, nonthermal 2016/2017 outburst that we initially discovered with Swift. We interpret it as the latest of the famous optical double-peaked outbursts of OJ 287, favoring binary scenarios that do not require a highly precessing secondary SMBH., Comment: ApJ, in press. The full radio-data tables will be available in the online material of the journal. Note that in a follow-up publication (MNRAS Letters, in press), the black hole mass of OJ 287 will be revised

The project MOMO (Multiwavelength Observations and Modelling of OJ 287) was set up to test predictions of binary supermassive black hole (SMBH) scenarios and to understand disk-jet physics of the blazar OJ 287. After a correction, the precessing binary (PB) SMBH model predicted the next main outburst of OJ 287 in 2022 October, no longer in July, making the outburst well observable and the model testable. We have densely covered this period in our ongoing multi-frequency radio, optical, UV, and X-ray monitoring. The predicted outburst was not detected. Instead, OJ 287 was at low optical-UV emission levels, declining further into November. The predicted thermal bremsstrahlung spectrum was not observed either, at any epoch. Further, applying scaling relations, we estimate a SMBH mass of OJ 287 of 10^8 M_sun. The latest in a sequence of deep low-states that recur every 1-2 yrs is used to determine an upper limit on the Eddington ratio and on the accretion-disk luminosity. This limit is at least a factor of 10 lower than required by the PB model with its massive primary SMBH of >10^{10} M_sun. All these results favor alternative binary SMBH models of OJ 287 that neither require strong orbital precession nor a very large mass of the primary SMBH., Comment: MNRAS Letter, in press

Obscuring winds driven away from active supermassive black holes are rarely seen due to their transient nature. They have been observed with multi-wavelength observations in a few Seyfert 1 galaxies and one broad absorption line radio-quiet quasar so far. An X-ray obscuration event in MR 2251-178 was caught in late 2020, which triggered multi-wavelength (NIR to X-ray) observations targeting this radio-quiet quasar. In the X-ray band, the obscurer leads to a flux drop in the soft X-ray band from late 2020 to early 2021. X-ray obscuration events might have a quasi-period of two decades considering earlier events in 1980 and 1996. In the UV band, a forest of weak blueshifted absorption features emerged in the blue wing of Ly$\alpha$ $\lambda1216$ in late 2020. Our XMM-Newton, NuSTAR, and HST/COS observations are obtained simultaneously, hence, the transient X-ray obscuration event is expected to account for the UV outflow, although they are not necessarily caused by the same part of the wind. Both blueshifted and redshifted absorption features were found for He {\sc i} $\lambda10830$, but no previous NIR spectra are available for comparison. The X-ray observational features of MR 2251-178 shared similarities with some other type 1 AGNs with obscuring wind. However, observational features in the UV to NIR bands are distinctly different from those seen in other AGN with obscuring winds. A general understanding of the observational variety and the nature of obscuring wind is still lacking., Comment: ApJ accepted

We are carrying out the densest and longest multiyear, multiwavelength monitoring project of OJ 287 ever done. The project MOMO (Multiwavelength Observations and Modelling of OJ 287) covers wavelengths from the radio to the high-energy regime. A few selected observations are simultaneous with those of the Event Horizon Telescope (EHT). MOMO aims at understanding disk-jet physics and at testing predictions of the binary black hole scenario of OJ 287. Here, we present a discussion of extreme outburst and minima states in context, and then focus on the recent flux and spectral evolution between 2021 and May 2022, including an ongoing bright radio flare. Further, we show that there is no evidence for precursor flare activity in our optical-UV-X-ray light curves that would be associated with any secondary supermassive black hole (SMBH) disk impact and that was predicted to start as thermal flare on 2021 December 23., Comment: Submitted to AN. Contribution to the XMM-Newton Workshop 2022 (Madrid, Spain). 5 figures

Understanding if and when the accretion disc extends down to the innermost stable circular orbit (ISCO) is important since it is the fundamental assumption behind measuring black hole spin. Here, we examine the 2013 and 2018 NuSTAR and Swift data (0.5 - 50 keV) of the narrow-line Seyfert 1 galaxy, WKK 4438. The X-ray emission can be fitted well with models depicting a corona and blurred reflection originating from a disc around a low spin (a* ~ 0) black hole. However, such models result in unconventional values for some of the parameters (e.g. inverted emissivity profile and high coronal height). Alternatively, equally good fits can be achieved if the disc is truncated at ~10 rg and the black hole is spinning at the Thorne limit (a* = 0.998). In these cases, the model parameters are consistent with the interpretation that the corona is centrally located close to the black hole and illuminating the disc at a larger distance., Comment: 8 pages, 5 figures. Accepted for publication in MNRAS

We report results from our ongoing project MOMO (Multiwavelength Observations and Modelling of OJ 287). In this latest publication of a sequence, we combine our Swift UVOT--XRT and Effelsberg radio data (2.6-44 GHz) between 2019 and 2022.04 with public SMA data and gamma-ray data from the Fermi satellite. The observational epoch covers OJ 287 in a high state of activity from radio to X-rays. The epoch also covers two major events predicted by the binary supermassive black hole (SMBH) model of OJ 287. Spectral and timing analyses clearly establish: a new UV-optical minimum state in 2021 December at an epoch where the secondary SMBH is predicted to cross the disk surrounding the primary SMBH; an overall low level of gamma-ray activity in comparison to pre-2017 epochs; the presence of a remarkable, long-lasting UV--optical flare event of intermediate amplitude in 2020--2021; a high level of activity in the radio band with multiple flares; and particularly a bright, ongoing radio flare peaking in 2021 November that may be associated with a gamma-ray flare, the strongest in 6 years. Several explanations for the UV--optical minimum state are explored, including the possibility that a secondary SMBH launches a temporary jet, but the observations are best explained by variability associated with the main jet., Comment: 16 pages, 11 figures. Published in MNRAS (submitted in January, accepted in March)

We report on X-ray flares that were observed from the active galactic nucleus I Zwicky 1 (I Zw 1) in 2020 January by the NuSTAR and XMM-Newton observatories. The X-ray spectrum is well-described by a model comprised of the continuum emission from the corona and its reflection from the accretion disc around a rapidly spinning (a > 0.94) black hole. In order to model the broadband spectrum, it is necessary to account for the variation in ionisation across the disc. Analysis of the X-ray spectrum in time periods before, during and after the flares reveal the underlying changes to the corona associated with the flaring. During the flares, the reflection fraction drops significantly, consistent with the acceleration of the corona away from the accretion disc. We find the first evidence that during the X-ray flares, the temperature drops from 140(-20,+100)keV before to 45(-9,+40)keV during the flares. The profile of the iron K line reveals the emissivity profile of the accretion disc, showing it to be illuminated by a compact corona extending no more than 7(-2,+4)rg over the disc before the flares, but with tentative evidence that the corona expands as it is accelerated during the flares. Once the flares subsided, the corona had collapsed to a radius of 6(-2,+2)rg. The rapid timescale of the flares suggests that they arise within the black-hole magnetosphere rather than in the accretion disc, and the variation of the corona is consistent with the continuum arising from the Comptonisation of seed photons from the disc., Comment: 15 pages, 10 figures. Accepted for publication in MNRAS

We are carrying out a dense monitoring of the blazar OJ 287 with Swift since late 2015 as part of our project MOMO (Multiwavelength Observations and Modeling of OJ 287). This is the densest existing monitoring of OJ 287 involving X-ray and UV data. In this latest publication of a sequence, we characterize the multiwavelength variability of OJ 287 based on >4000 Swift single-wave-band data sets including archival data since 2005. A structure function analysis reveals a characteristic timescale of ~5 days in the optical-UV at epochs of low-level activity, and larger during outbursts. The discrete correlation function shows zero lag between optical and UV, with tau = 0+-1 days at the epoch of densest cadence. During outbursts (in 2016/17 and 2020) the X-rays follow the UV with near-zero lags. However, during quiescence, the delay is 7-18 days with X-rays leading or lagging, interpreted as due to a different X-ray component dominated by inverse Compton emission. Scaling relations are used to derive the characteristic length scales of broad-line region and torus in OJ 287. A remarkable, symmetric UV--optical deep fade is identified in late 2017, lasting for 2 months. We rule out occultation from the passage of a dusty cloud and a model where the secondary black hole deflects the jet between the primary and observer. We speculate about a temporary dispersion or jet swing event in the core or in a bright quasi-stationary jet feature. The deep fade reveals an additional, spatially distinct X-ray component. The epoch 2020.9-2021.1 was searched for precursor flare activity predicted by the binary black hole model of OJ 287., Comment: Results are part of a sequence of publications, previous and upcoming, from our ongoing project MOMO focussed on OJ 287. Accepted in July 2021, published December 9 in the Astrophysical Journal (ApJ 923, 51, 2021)

The innermost regions of accretion disks around black holes are strongly irradiated by X-rays that are emitted from a highly variable, compact corona, in the immediate vicinity of the black hole. The X-rays that are seen reflected from the disk and the time delays, as variations in the X-ray emission echo or reverberate off the disk provide a view of the environment just outside the event horizon. I Zwicky 1 (I Zw 1), is a nearby narrow line Seyfert 1 galaxy. Previous studies of the reverberation of X-rays from its accretion disk revealed that the corona is composed of two components; an extended, slowly varying component over the surface of the inner accretion disk, and a collimated core, with luminosity fluctuations propagating upwards from its base, which dominates the more rapid variability. Here we report observations of X-ray flares emitted from around the supermassive black hole in I Zw 1. X-ray reflection from the accretion disk is detected through a relativistically broadened iron K line and Compton hump in the X-ray emission spectrum. Analysis of the X-ray flares reveals short flashes of photons consistent with the re-emergence of emission from behind the black hole. The energy shifts of these photons identify their origins from different parts of the disk. These are photons that reverberate off the far side of the disk and bent around the black hole and magnified by the strong gravitational field. Observing photons bent around the black hole confirms a key prediction of General Relativity., Comment: Preprint version. 22 pages, 8 figures. Published in Nature

We present the first high signal-to-noise broadband X-ray spectrum of the radio-quiet type-2 Seyfert ESO 033-G002, combining data from $XMM$-$Newton$ and $NuSTAR$. The nuclear X-ray spectrum is complex, showing evidence for both neutral and ionised absorption, as well as reflection from both the accretion disc and more distant material, but our broadband coverage allows us to disentangle all of these different components. The total neutral column during this epoch is $N_{\rm{H}} \sim 5-6 \times 10^{22}$ cm$^{-2}$, consistent with the optical classification of ESO 033-G002 as a type-2 Seyfert but not so large as to prevent us from robustly determining the properties of the innermost accretion flow. The ionised absorption - dominated by lines from Fe XXV and Fe XXVI - reveals a moderately rapid outflow ($v_{\rm{out}} \sim 5400$ km s$^{-1}$) which has a column comparable to the neutral absorption. We find the disc reflection from the innermost regions to be extreme, with a reflection fraction of $R_{\rm{frac}} \sim 5$. This requires strong gravitational lightbending and, in turn, both an extremely compact corona (within $\sim$2 $R_{\rm{G}}$ of the black hole) and a rapidly rotating black hole ($a^* > 0.96$). Despite this tight size constraint, with a temperature of $kT_{\rm{e}} = 40-70$ keV the X-ray corona in ESO 033-G002 appears similar to other AGN in terms of its placement in the compactness-temperature plane, consistent with sitting close to the limit determined by runaway pair production. Finally, combining X-ray spectroscopy, timing and updated optical spectroscopy, we also estimate the mass of the black hole to be $\log[M_{\rm{BH}} / M_{\odot}] \sim 7.0 - 7.5$., Comment: 15 pages, 8 figures; published in MNRAS

Our project MOMO (Multiwavelength observations and modelling of OJ 287) consists of dedicated, dense, long-term flux and spectroscopic monitoring and deep follow-up observations of the blazar OJ 287 at >13 frequencies from the radio to the X-ray band since late 2015. In particular, we are using Swift to obtain optical-UV-X-ray spectral energy distributions (SEDs) and the Effelsberg telescope to obtain radio measurements between 2 and 40 GHz. MOMO is the densest long-term monitoring of OJ 287 involving X-rays and broad-band SEDs. The theoretical part of the project aims at understanding jet and accretion physics of the blazar central engine in general and the supermassive binary black hole scenario in particular. Results are presented in a sequence of publications and so far included: detection and detailed analysis of the bright 2016/17 and 2020 outbursts and the long-term light curve; Swift, XMM and NuSTAR spectroscopy of the 2020 outburst around maximum; and interpretation of selected events in the context of the binary black hole scenario of OJ 287 (papers I-IV). Here, we provide a description of the project MOMO, a summary of previous results, the latest results, and we discuss future prospects., Comment: Review. 18 pages, 11 figures. Submitted to Universe Special Issue "Panchromatic View of the Life-Cycle of AGN"

Narrow-line Seyfert 1 (NLS1) galaxies are a class of active galactic nuclei (AGN) that, in some cases, can harbor powerful relativistic jets. One of them, PKS 2004-447, shows gamma-ray emission, and underwent its first recorded multifrequency flare in 2019. However, past studies revealed that in radio this source can be classified as a compact steep-spectrum source (CSS), suggesting that, unlike other gamma-ray sources, the relativistic jets of PKS 2004-447 have a large inclination with respect to the line of sight. We present here a set of spectroscopic observations of this object, aimed at carefully measuring its black hole mass and Eddington ratio, determining the properties of its emission lines, and characterizing its long term variability. We find that the black hole mass is $(1.5\pm0.2)\times10^7$ M$_\odot$, and the Eddington ratio is 0.08. Both values are within the typical range of NLS1s. The spectra also suggest that the 2019 flare was caused mainly by the relativistic jet, while the accretion disk played a minor role during the event. In conclusion, we confirm that PKS 2004-447 is one of the rare examples of gamma-ray emitting CSS/NLS1s hybrid, and that these two classes of objects are likely connected in the framework of AGN evolution., Comment: 14 pages, 8 figures, 4 tables, resubmitted to A&A after minor comments from the referee

Supermassive binary black holes (SMBBHs) are laboratories par excellence for relativistic effects, including precession effects in the Kerr metric and the emission of gravitational waves. Binaries form in the course of galaxy mergers, and are a key component in our understanding of galaxy evolution. Dedicated searches for SMBBHs in all stages of their evolution are therefore ongoing and many systems have been discovered in recent years. Here we provide a review of the status of observations with a focus on the multiwavelength detection methods and the underlying physics. Finally, we highlight our ongoing, dedicated multiwavelength program MOMO (for Multiwavelength Observations and Modelling of OJ 287). OJ 287 is one of the best candidates to date for hosting a sub-parsec SMBBH. The MOMO program carries out a dense monitoring at >13 frequencies from radio to X-rays and especially with Swift since 2015. Results so far included: (1) The detection of two major UV-X-ray outbursts with Swift in 2016/17 and 2020; exhibiting softer-when-brighter behaviour. The non-thermal nature of the outbursts was clearly established and shown to be synchrotron radiation. (2) Swift multi-band dense coverage and XMM-Newton spectroscopy during EHT campaigns caught OJ 287 at an intermediate flux level with synchrotron and IC spectral components. (3) Discovery of a remarkable, giant soft X-ray excess with XMM and NuSTAR during the 2020 outburst. (4) Spectral evidence (at 2sigma) for a relativistically shifted iron absorption line in 2020. (5) The non-thermal 2020 outburst is consistent with an after-flare predicted by the SMBBH model of OJ 287., Comment: Invited review at the COST Special Session on "The gravitational-wave Universe" during the XIX Serbian Astronomical Conference (Oct. 2020). 14 pages incl. 7 figures

The narrow-line Seyfert 1 galaxy Mrk 335 was one of the X-ray brightest AGN, but it has systematically faded since 2007. Here, we report the discovery with Swift of a sequence of bright and rapid X-ray flare events that reveal the emergence of Mrk 335 from its ultra-deep multiyear low state. Results are based on our dedicated multiyear monitoring with Swift. Unlike other bright AGN, the optical-UV is generally not correlated with the X-rays in Mrk 335 on a timescale of days to months. This fact either implies the absence of a direct link between the two emission components; or else implies that the observed X-rays are significantly affected by (dust-free) absorption along our l.o.s. The UV and optical, however, are closely correlated at the 99.99 percent confidence level. The UV is leading the optical by delta t = 1.5+/-1.5 d. The Swift X-ray spectrum shows strong deviations from a single power law in all brightness states of the outbursts, indicating that significant absorption or reprocessing is taking place. Mrk 335 displays a softer-when-brighter variability pattern at intermediate X-ray levels, which has been seen in our Swift data since 2007 (based on a total of 590 observations). This pattern breaks down at the highest and lowest count rates. We interpret the 2020 brightening of Mrk 335 as a decrease in column density and covering factor of a partial-covering absorber along our line of sight in form of a clumpy accretion-disk wind that reveals an increasing portion of the intrinsic emission of Mrk 335 from the disk and/or corona region, while the optical emission-line regions receive a less variable spectral energy distribution. This then also explains why Mrk 335 was never seen to change its optical Seyfert type (not 'changing look') despite its factor 50 X-ray variability with Swift., Comment: A&A Letters, appeared in Vol. 643, L7

We present new broadband X-ray observations of the type-I Seyfert galaxy IRAS 09149-6206, taken in 2018 with $XMM$-$Newton$, $NuSTAR$ and $Swift$. The source is highly complex, showing a classic 'warm' X-ray absorber, additional absorption from highly ionised iron, strong relativistic reflection from the innermost accretion disc and further reprocessing by more distant material. By combining X-ray timing and spectroscopy, we have been able to fully characterise the supermassive black hole in this system, constraining both its mass and - for the first time - its spin. The mass is primarily determined by X-ray timing constraints on the break frequency seen in the power spectrum, and is found to be $\log[M_{\rm{BH}}/M_{\odot}] = 8.0 \pm 0.6$ (1$\sigma$ uncertainties). This is in good agreement with previous estimates based on the H$\alpha$ and H$\beta$ line widths, and implies that IRAS 09149-6206 is radiating at close to (but still below) its Eddington luminosity. The spin is constrained via detailed modelling of the relativistic reflection, and is found to be $a^* = 0.94^{+0.02}_{-0.07}$ (90% confidence), adding IRAS 09149-6206 to the growing list of radio-quiet AGN that host rapidly rotating black holes. The outflow velocities of the various absorption components are all relatively modest ($v_{\rm{out}} \lesssim 0.03c$), implying these are unlikely to drive significant galaxy-scale AGN feedback., Comment: 19 pages, 15 figures, accepted for publication in MNRAS

Multi-wavelength monitoring of Mrk 335 with Swift between 2007-2019 are used to construct annual spectral energy distributions and track year-to-year changes. Non-contemporaneous archival data prior to 2007 are used to build a bright state SED. In this work, the changes are examined and quantified to build the foundation for future SED modelling. The yearly SEDs trace a downward trend on the average, with the X-ray portion varying significantly and acquiring further lower values in the past two years when compared to the optical/UV portion of SED. The bolometric Eddington ratios derived using optical/UV to X-ray SEDs and the calculated X-ray luminosities show a gradual decrease over the monitoring period. Changes in the parameters over time are examined. Principal component analysis suggests that the primary variability is in the X-ray properties of Mrk 335. When looking at the broader picture of Mrk 335 and its behaviour, the X-rays, accounting most of the variability in the 13-year data, are possibly driven by physical processes related to the corona or absorption whereas the modest optical-UV variations suggest their origin within the accretion disc. These results are consistent with the previous interpretation of Mrk 335 using the timing analyses on the monitoring data and spectral modelling of deep observations., Comment: Accepted for publication in MNRAS. Long tables 3 and 4 are available as supplementary material of MNRAS

We present results of temporal and spectral analyses on four XMM-Newton EPIC pn observations of IRAS 17020+4544, a narrow-line Seyfert 1 galaxy with evidence of a radio jet. Analysis of the light curves reveals that this radio-loud source does not behave like the bulk population of its radio-quiet counterparts. A trend of spectral hardening with increased flux is found. Variability is found to increase with energy, though it decreases as the spectrum hardens. The first 40 ks of the most recent observation behave uniquely among the epochs, exhibiting a softer spectral state than at any other time. Possible non-stationarity at low energies is found, with no such effect present at higher energies, suggesting at least two distinct spectral components. A reverberation signature is confirmed, with the lag-frequency, lag-energy, and covariance spectra changing significantly during the soft-state epoch. The temporal analysis suggests a variable power-law in the presence of a reflection component, thus motivating such a fit for the 0.3-10 keV EPIC pn spectra from all epochs. We find an acceptable spectral fit using the timing-motivated parameters and report the detection of a broad Fe K emission line, requiring an additional model component beyond the reflection spectrum. We discuss links between this source and other narrow-line Seyfert 1 sources that show evidence of jet activity, finding similarities among this currently very limited sample of interesting objects., Comment: Accepted for publication in MNRAS. 19 pages, 14 figures

All interesting and fascinating collective properties of a complex system arise from the intricate way in which its components interact. Various systems in physics, biology, social sciences and engineering have been successfully modelled as networks of coupled dynamical systems, where the graph links describe pairwise interactions. This is, however, too strong a limitation, as recent studies have revealed that higher-order many-body interactions are present in social groups, ecosystems and in the human brain, and they actually affect the emergent dynamics of all these systems. Here, we introduce a general framework that allows to study coupled dynamical systems accounting for the precise microscopic structure of their interactions at any possible order. We consider the most general ensemble of identical dynamical systems, organized on the nodes of a simplicial complex, and interacting through synchronization-non-invasive coupling function. The simplicial complex can be of any dimension, meaning that it can account, at the same time, for pairwise interactions, three-body interactions and so on. In such a broad context, we show that complete synchronization exists as an invariant solution, and we give the necessary condition for it to be observed as a stable state in terms of a Master Stability Function. This generalizes the existing results valid for pairwise interactions (i.e. graphs) to the case of complex systems with the most general possible architecture. Moreover, we show how the approach can be simplified for specific, yet frequently occurring, instances, and we verify all our theoretical predictions in synthetic and real-world systems. Given the completely general character of the method proposed, our results contribute to the theory of dynamical systems with many-body interactions and can find applications in an extremely wide range of practical cases.

We explore a disc origin for the highly-blueshifted, variable absorption lines seen in the X-ray spectrum of the Narrow Line Seyfert 1 galaxy IRAS13224-3809. The blueshift corresponds to a velocity of about 0.25c. Such features in other Active Galactic Nuclei are often interpreted as UltraFast Outflows (UFOs). The velocity is of course present in the orbital motions of the inner disk. The absorption lines in IRAS13224-3809 are best seen when the flux is low and the reflection component of the disk is strong relative to the power-law continuum. The spectra are consistent with a model in which the reflection component passes through a thin, highly-ionized absorbing layer at the surface of the inner disc, the blue-shifted side of which dominates the flux due to relativistic aberration (the disc inclination is about 70 deg). No fast outflow need occur beyond the disc., Comment: 6 pages, 5 figures, accepted for publication in MNRAS

The next generation of electromagnetic and gravitational wave observatories will open unprecedented windows to the birth of the first supermassive black holes. This has the potential to reveal their origin and growth in the first billion years, as well as the signatures of their formation history in the local Universe. With this in mind, we outline three key focus areas which will shape research in the next decade and beyond: (1) What were the "seeds" of the first quasars; how did some reach a billion solar masses before z$\sim7$? (2) How does black hole growth change over cosmic time, and how did the early growth of black holes shape their host galaxies? What can we learn from intermediate mass black holes (IMBHs) and dwarf galaxies today? (3) Can we unravel the physics of black hole accretion, understanding both inflows and outflows (jets and winds) in the context of the theory of general relativity? Is it valid to use these insights to scale between stellar and supermassive BHs, i.e., is black hole accretion really scale invariant? In the following, we identify opportunities for the Canadian astronomical community to play a leading role in addressing these issues, in particular by leveraging our strong involvement in the Event Horizon Telescope, the {\it James Webb Space Telescope} (JWST), Euclid, the Maunakea Spectroscopic Explorer (MSE), the Thirty Meter Telescope (TMT), the Square Kilometer Array (SKA), the Cosmological Advanced Survey Telescope for Optical and ultraviolet Research (CASTOR), and more. We also discuss synergies with future space-based gravitational wave (LISA) and X-ray (e.g., Athena, Lynx) observatories, as well as the necessity for collaboration with the stellar and galactic evolution communities to build a complete picture of the birth of supermassive black holes, and their growth and their influence over the history of the Universe., Comment: White paper submitted to Canadian 2020 Long Range Plan committee. 11 pages, 3 figures

We present XMM-Newton, NuSTAR, Swift and Hubble Space Telescope observations of the Narrow-line Seyfert 1 galaxy Mrk 335 in a protracted low state in 2018 and 2019. The X-ray flux is at the lowest level so far observed, and the extremely low continuum flux reveals a host of soft X-ray emission lines from photoionised gas. The simultaneous UV flux drop suggests that the variability is intrinsic to the source, and we confirm this with broad-band X-ray spectroscopy. The dominance of the soft X-ray lines at low energies and distant reflection at high energies is therefore due to the respective emission regions being located far enough from the X-ray source that they have not yet seen the flux drop. Between the two XMM-Newton spectra, taken 6 months apart, the emission line ratio in the Ovii triplet changes drastically. We attribute this change to a drop in the ionisation of intervening warm absorption, which means that the absorber must cover a large fraction of the line emitting region, and extend much further from the black hole than previously assumed. The HST spectrum, taken in 2018, shows that new absorption features have appeared on the blue wings of Ciii*, Ly{\alpha}, Nv, Siiv and Civ, likely due to absorbing gas cooling in response to the low flux state., Comment: 15 pages, 11 figures, accepted for publication in MNRAS

A multi-epoch X-ray spectral and variability analysis is conducted for the narrow-line Seyfert 1 (NLS1) active galactic nucleus (AGN) Mrk 478. All available X-ray data from XMM-Newton and Suzaku satellites, spanning from 2001 to 2017, are modelled with a variety of physical models including partial covering, soft-Comptonisation, and blurred reflection, to explain the observed spectral shape and variability over the 16 years. All models are a similar statistical fit to the data sets, though the analysis of the variability between data sets favours the blurred reflection model. In particular, the variability can be attributed to changes in flux of the primary coronal emission. Different reflection models fit the data equally well, but differ in interpretation. The use of reflionx predicts a low disc ionisation and power law dominated spectrum, while relxill predicts a highly ionised and blurred reflection dominated spectrum. A power law dominated spectrum might be more consistent with the normal X-ray-to-UV spectral shape (aox). Both blurred reflection models suggest a rapidly spinning black hole seen at a low inclination angle, and both require a sub-solar (~0.5) abundance of iron. All physical models require a narrow emission feature at 6.7 keV likely attributable to Fe xxv emission, while no evidence for a narrow 6.4 keV line from neutral iron is detected., Comment: 15 pages. 9 figures. Accepted for publication in MNRAS

We have carried out a systematic analysis of the nearby (z=0.0279) active galaxy Zw 229.015 using multi-epoch, multi-instrument and deep pointed observations with XMM-Newton, Suzaku, Swift and NuSTAR. Spectral and temporal variability are examined in detail on both the long (weeks-to-years) and short (hours) timescales. A deep Suzaku observation of the source shows two distinct spectral states; a bright-soft state and a dim-hard state in which changes in the power law component account for the differences. Partial covering, blurred reflection and soft Comptonisation models describe the X-ray spectra comparably well, but the smooth, rather featureless, spectrum may be favouring the soft Comptonisation scenario. Moreover, independent of the spectral model, the observed spectral variability is ascribed to the changes in the power law continuum only and do not require changes in the properties of the absorber or blurred reflector incorporated in the other scenarios. The multi-epoch observations between 2009 and 2018 can be described in similar fashion. This could be understood if the primary emission is originating at a large distance from a standard accretion disc or if the disc is optically thin and geometrically thick as recently proposed for Zw 229.015. Our investigation shows that Zw 229.015 behaves similar to sources like Akn 120 and Mrk 530, that exhibit a strong soft-excess, but weak Compton hump and Fe K${\alpha}$ emission., Comment: 14 pages, 7 figures. Accepted for publication in MNRAS

A number of works point to the presence of narrow emission features at unusual energies in the X-ray spectra of active galactic nuclei (AGN) or to potentially low iron abundances as possible evidence for the spallation of iron. With the imminence of high-resolution calorimeter spectroscopy, the potential to test spallation models in astrophysical sources will soon be possible. Previously determined nuclear spallation reactions of Fe are recalculated making use of improved total inelastic and partial reaction cross sections that result in different absolute and relative abundances of the main spallation elements Mn, Cr, V, and Ti. The effects of ionisation and dynamics near the black hole are examined in simulated spectra with CCD and calorimeter (i.e. Hitomi) resolution. With high-resolution, differences in relative abundances and ionisation should be detectable if spallation is originating at large distances from the black hole (e.g. torus or disc wind) where blurring is not significant. If spallation were occurring in the inner accretion disc, it would likely be undetected as blurring effects would cause significant blending of spectral features., Comment: 6 pages, 8 figures, accepted for publication in MNRAS

A triggered 140 ks XMM-Newton observation of the narrow-line Seyfert 1 (NLS1) Mrk 335 in December 2015 caught the active galaxy at its lowest X-ray flux since 2007. The NLS1 is relatively quiescent for the first ~120 ks of the observation before it flares in brightness by a factor of about five in the last 20 ks. Although only part of the flare is captured before the observation is terminated, the data reveal significant differences between the flare and quiescent phases. During the low-flux state, Mrk 335 demonstrates a reflection-dominated spectrum that results from a compact corona around a Kerr black hole. In addition to the rapid brightening, the flare is further described by spectral softening and a falling reflection fraction that are consistent with previous observations advocating at least part of the corona in Mrk 335 could be the base of an aborted jet. The spectrum during the flaring interval reveals several residuals between the 2-3 sigma level that could be attributed to absorption lines from a highly ionised plasma that is moving outward at v~0.12c. It could be that the increased luminosity during the flare enhances the radiation pressure sufficiently to launch a possible wind. If the wind is indeed responding to the change in corona luminosity then it must be located within ~80 Rg. The escape velocity at this distance is comparable to the estimated wind velocity. If confirmed, this is the first example of a radio-quiet AGN exhibiting behaviour consistent with both diffuse and collimated outflow., Comment: 11 pages, 10 figures, accepted for publication in MNRAS

We present results from a deep, coordinated $XMM$-$Newton$+$NuSTAR$ observation of the Seyfert 2 galaxy IRAS 00521-7054. The $NuSTAR$ data provide the first detection of this source in high-energy X-rays ($E > 10$ keV), and the broadband data show this to be a highly complex source which exhibits relativistic reflection from the inner accretion disc, further reprocessing by more distant material, neutral absorption, and evidence for ionised absorption in an extreme, ultrafast outflow ($v_{\rm{out}} \sim 0.4c$). Based on lamppost disc reflection models, we find evidence that the central supermassive black hole is rapidly rotating ($a > 0.77$), consistent with previous estimates from the profile of the relativistic iron line, and that the accretion disc is viewed at a fairly high inclination ($i \sim 59^{\circ}$). Based on extensive simulations, we find the ultrafast outflow is detected at $\sim$4$\sigma$ significance (or greater). We also estimate that the extreme outflow should be sufficient to power galaxy-scale feedback, and may even dominate the energetics of the total output from the system., Comment: 12 pages, 8 figures; accepted for publication in MNRAS

It is arguably in the X-ray regime that Narrow-line Seyfert 1 galaxies (NLS1s) exhibit the most extreme behaviour. Spectral complexity, rapid and large amplitude flux variations, and exceptional spectral variability are well known characteristics. However, NLS1s are not eccentric, but form a continuous sequence with typical Seyfert 1 galaxies. Understanding the extreme behaviour displayed by NLS1s will provide insight to the general AGN phenomenon. In this review, I will examine some of the important NLS1 X-ray discoveries over the past twenty years. I will then explore recent work that looks at the nature of the primary X-ray source (i.e. the corona) in NLS1s, demonstrating how the corona can be compact, dynamic, and in some cases consistent with collimated outflow. X-ray observations of NLS1s will be key in determining the nature of the corona, resolving the disc-jet connection, and determining the origin of the radio loud/quiet dichotomy in AGN., Comment: Invited review article for the conference "Revisiting narrow-line Seyfert 1 galaxies and their place in the Universe" (Padova, April 2018). Accepted for publication in Proceedings of Science, PoS (NLS1-2018) 034

We present new XMM-Newton observations of the intriguing warm absorber in I Zwicky 1. This luminous and nearby narrow-line Seyfert 1 galaxy shows ionized absorption by two components of outflowing gas; a low and a high-ionization phase with log $\xi$~0 and log $\xi$~2 respectively. Detailed modelling of these data reveal a complex and variable multi-phase warm absorber. However, we find the changes in the ionization state of the gas not to be straightforwardly correlated with the variability of the intrinsic continuum source, in apparent contrast with photoionization equilibrium. The observed variability hints instead at a close connection between the two gas components, possibly both directly connected to the accretion disc activity. We thus suggest a phenomenological model capable of explaining these observations, consisting of a clumpy outflow where the high and the low-ionization components are closely linked. Changes in ionization over the years are mainly driven by the different densities of the clumps crossing the observer's line-of-sight, in which the `skin' layer facing the source accounts for the more ionized component., Comment: 10 pages, 7 figures. Accepted for publication in MNRAS

We collate active galactic nuclei (AGN) with reported detections of both relativistic reflection and ultra-fast outflows. By comparing the inclination of the inner disc from reflection with the line-of-sight velocity of the outflow, we show that it is possible to meaningfully constrain the geometry of the absorbing material. We find a clear relation between the velocity and inclination, and demonstrate that it can potentially be explained either by simple wind geometries or by absorption from the disc surface. Due to systematic errors and a shortage of high- quality simultaneous measurements our conclusions are tentative, but this study represents a proof-of-concept that has great potential., Comment: 5 pages, 3 figures, accepted to MNRAS letters

A 2012 Suzaku observation of the Seyfert 1.5 galaxy Markarian 530 was analysed and found to exhibit two distinct modes of variability, which were found to be independent from one another. Firstly, the spectrum undergoes a smooth transition from a soft to a hard spectrum. Secondly, the spectrum displays more rapid variability seemingly confined to a very narrow energy band (~ 1-3 keV). Three physical models (blurred reflection, partial covering, and soft Comptonisation) were explored to characterise the average spectrum of the observation as well as the spectral state change. All three models were found to fit the average spectrum and the spectral changes equally well. The more rapid variability appears as two cycles of a sinusoidal function, but we cannot attribute this to periodic variability. The Fe K-alpha band exhibits a narrow 6.4 keV emission line consistent with an origin from the distant torus. In addition, features blueward of the neutral iron line are consistent with emission from He-like and H-like iron that could be originating from the highly ionised layer of the torus, but a broad Gaussian profile at ~ 6.7 keV also fits the spectrum well., Comment: 12 pages, 6 figures, 2 tables. Accepted for publication in MNRAS

The narrow-line Seyfert 1 galaxy (NLS1) Mrk 335 has been continuously monitored with Swift since May 2007 when it fell into a long-lasting, X-ray low-flux interval. Results from the nearly 11 years of monitoring are presented here. Structure functions are used to measure the UV-optical and X-ray power spectra. The X-ray structure function measured between 10-100 days is consistent with the flat, low-frequency part of the power spectrum measured previously in Mrk 335. The UV-optical structure functions of Mrk 335 are comparable with those of other Seyfert 1 galaxies and of Mrk 335 itself when it was in a normal bright state. There is no indication that the current X-ray low-flux state is attributed to changes in the accretion disc structure of Mrk 335. The characteristic timescales measured in the structure functions can be attributed to thermal (for the UV) and dynamic (for the optical) timescales in a standard accretion disc. The high-quality UVW2 (~1800 A in the source frame) structure function appears to have two breaks and two different slopes between 10-160 days. Correlations between the X-ray and other bands are not highly significant when considering the entire 11-year light curves, but more significant behaviour is present when considering segments of the light curves. A correlation between the X-ray and UVW2 in 2014 (Year-8) may be predominately caused by an giant X-ray flare that was interpreted as jet-like emission. In 2008 (Year-2), possible lags between the UVW2 emission and other UV-optical waveband may be consistent with reprocessing of X-ray or UV emission in the accretion disc., Comment: Figure 8b has been corrected. 12 pages. Accepted for publication in MNRAS

We present a detailed spectral analysis of the recent 1.5\,Ms XMM-Newton observing campaign on the narrow line Seyfert 1 galaxy IRAS~13224$-$3809, taken simultaneously with 500\,ks of NuSTAR data. The X-ray lightcurve shows three flux peaks, registering at about 100 times the minimum flux seen during the campaign, and rapid variability with a time scale of kiloseconds. The spectra are well fit with a primary powerlaw continuum, two relativistic-blurred reflection components from the inner accretion disk with very high iron abundance, and a simple blackbody-shaped model for the remaining soft excess. The spectral variability is dominated by the power law continuum from a corona region within a few gravitational radii from the black hole. Additionally, blueshifted Ne \textsc{x}, Mg \textsc{xii}, Si \textsc{xiv} and S \textsc{xvi} absorption lines are identified in the stacked low-flux spectrum, confirming the presence of a highly ionized outflow with velocity up to $v= 0.267$ and $0.225$\,c. We fit the absorption features with \texttt{xstar} models and find a relatively constant velocity outflow through the whole observation. Finally, we replace the \texttt{bbody} and supersolar abundance reflection models by fitting the soft excess successfully with the extended reflection model \texttt{relxillD}, which allows for higher densities than the standard \texttt{relxill} model. This returns a disk electron density $n_{\rm e}>10^{18.7}$\,cm$^{-3}$ and lowers the iron abundance from $Z_{\rm Fe}=24^{+3}_{-4}Z_\odot$ with $n_{\rm e}\equiv10^{15}$\,cm$^{-3}$ to $Z_{\rm Fe}=6.6^{+0.8}_{-2.1}Z_\odot$., Comment: 18 pages, 14 figures. Accepted by MNRAS

We present a detailed X-ray timing analysis of the highly variable NLS1 galaxy, IRAS 13224-3809. The source was recently monitored for 1.5 Ms with XMM-Newton which, combined with 500 ks archival data, makes this the best studied NLS1 galaxy in X-rays to date. We apply standard time- and Fourier-domain in order to understand the underlying variability process. The source flux is not distributed lognormally, as would be expected for accreting sources. The first non-linear rms-flux relation for any accreting source in any waveband is found, with $\mathrm{rms} \propto \mathrm{flux}^{2/3}$. The light curves exhibit significant strong non-stationarity, in addition to that caused by the rms-flux relation, and are fractionally more variable at lower source flux. The power spectrum is estimated down to $\sim 10^{-7}$ Hz and consists of multiple peaked components: a low-frequency break at $\sim 10^{-5}$ Hz, with slope $\alpha < 1$ down to low frequencies; an additional component breaking at $\sim 10^{-3}$ Hz. Using the high-frequency break we estimate the black hole mass $M_\mathrm{BH} = [0.5-2] \times 10^{6} M_{\odot}$, and mass accretion rate in Eddington units, $\dot m_{\rm Edd} \gtrsim 1$. The non-stationarity is manifest in the PSD with the normalisation of the peaked components increasing with decreasing source flux, as well as the low-frequency peak moving to higher frequencies. We also detect a narrow coherent feature in the soft band PSD at $0.7$ mHz, modelled with a Lorentzian the feature has $Q \sim 8$ and an $\mathrm{rms} \sim 3$ %. We discuss the implication of these results for accretion of matter onto black holes., Comment: Accepted to MNRAS. 19 pages, 17 figures

We examine a 200 ks XMM-Newton observation of the narrow-line Seyfert 1 galaxy Mrk 493. The active galaxy was half as bright as in a previous 2003 snapshot observation and the current lower flux enables a study of the putative reflection component in detail. We determine the characteristics of the 2015 X-ray continuum by first analyzing the short-term variability using model-independent techniques. We then continue with a time-resolve analysis including spectral fitting and modelling the fractional variability. We determine that the variability arises from changes in the amount of primary flux striking the accretion disk, which induces changes in the ionization parameter and flux of the blurred reflection component. The observations seem consistent with the picture that the primary source is of roughly constant brightness and that variations arise from changes in the degree of light bending happening in the vicinity of the supermassive black hole., Comment: 10 pages. Accepted for publication in MNRAS

We present results from the optical, ultraviolet and X-ray monitoring of the NLS1 galaxy IRAS 13224-3809 taken with Swift and XMM-Newton during 2016. IRAS 13224-3809 is the most variable bright AGN in the X-ray sky and shows strong X-ray reflection, implying that the X-rays strongly illuminate the inner disc. Therefore, it is a good candidate to study the relationship between coronal X-ray and disc UV emission. However, we find no correlation between the X-ray and UV flux over the available ~40 day monitoring, despite the presence of strong X-ray variability and the variable part of the UV spectrum being consistent with irradiation of a standard thin disc. This means either that the X-ray flux which irradiates the UV emitting outer disc does not correlate with the X-ray flux in our line of sight and/or that another process drives the majority of the UV variability. The former case may be due to changes in coronal geometry, absorption or scattering between the corona and the disc., Comment: 9 pages, 8 figures, MNRAS submitted: this version in response to reviewers comments; in original form 2017 October 2