Showing 1–50 of 104 results for author: Papadakis, I E

Broadband X-ray/UV/optical time-resolved spectroscopy of NGC 5548: The origin of the UV/optical variability in active galactic nuclei

Authors: E. Kammoun, I. E. Papadakis, M. Dovčiak, C. Panagiotou

Abstract: Recently, nearby active galactic nuclei (AGN) have been subject to long X-rays/UV/optical monitoring campaigns. These campaigns reveal a strong correlation between the various UV and optical bands, with time lags increasing with wavelength. In a series of papers, we demonstrated that a scenario in which a central X-ray source illuminates the accretion disc explains the observed correlations. Howev… ▽ More Recently, nearby active galactic nuclei (AGN) have been subject to long X-rays/UV/optical monitoring campaigns. These campaigns reveal a strong correlation between the various UV and optical bands, with time lags increasing with wavelength. In a series of papers, we demonstrated that a scenario in which a central X-ray source illuminates the accretion disc explains the observed correlations. However, some of the monitored AGN show low/moderate X-rays-UV correlations, which could challenge this scenario. In this paper, we study the broadband X-ray/UV/optical spectral energy distributions (SEDs) of NGC 5548, one of the most intensively monitored AGN. We aim to test if the X-ray illumination model explains the broadband spectral behaviour of the source, despite the moderate X-ray-UV/optical correlation. We model the broadband time-averaged SED, from the STORM monitoring campaign of the source, using the KYNSED model which assumes an X-ray illuminated disc. We assume that the accretion process powers the X-ray corona. We also model 15 time-resolved SEDs from the same campaign to check whether this scenario can account for the observed spectral variability. The proposed model describes well the time-averaged and the time-resolved SEDs of NGC 5548. In this scenario, the corona height, the X-ray photon index, and the power transferred to the corona all vary. This explains the variability behaviour at different wavelengths. The best-fit model is obtained for a non-rotating black hole accreting at a constant rate of 5% its Eddington limit. Since each of the variable parameters affects the observed flux in a particular way, the combined variability of all parameters explains the moderate X-ray-UV/optical correlation. The X-ray illuminated disc model provides a complete description of the behaviour of NGC 5548, explaining its broadband SEDs, time-lag spectrum, and its power spectral distribution. △ Less

Submitted 18 March, 2024; originally announced March 2024.

Comments: Accepted for publication in A&A

The X-ray variability of AGN: power-spectrum and variance analysis of the Swift/BAT light curves

Authors: I. E. Papadakis, V. Binas-Valavanis

Abstract: We study the X-ray power spectrum of Active Galactic Nuclei (AGN) to investigate whether Seyfert I and II power spectra are similar or not, whether the AGN variability depends on black hole mass and accretion rate, and to compare the AGN power spectra with the Galactic X-ray black hole binaries power-spectra. We used 14-195 keV band light curves from the 157th SWIFT/BAT hard X-ray survey and we co… ▽ More We study the X-ray power spectrum of Active Galactic Nuclei (AGN) to investigate whether Seyfert I and II power spectra are similar or not, whether the AGN variability depends on black hole mass and accretion rate, and to compare the AGN power spectra with the Galactic X-ray black hole binaries power-spectra. We used 14-195 keV band light curves from the 157th SWIFT/BAT hard X-ray survey and we computed the mean power spectrum and excess variance of AGN in narrow black-hole mass/luminosity bins. We fitted a power-law model to the AGN power spectra, and we investigated whether the power spectrum parameters and the excess variance depend on black hole mass, luminosity and accretion rate. The Seyfert I and Seyfert II power spectra are identical, in agreement with AGN unification models. The mean AGN X-ray power spectrum has the same, power-law like shape with a slope of -1 in all AGN, irrespective of their luminosity and BH mass. We do not detect any flattening to a slope of zero at frequencies as low as 10e-9 Hz. We detect an anti-correlation between the PSD amplitude and the accretion rate, similar to what has been seen in the past in the 2-10 keV band. This implies that the variability amplitude in AGN decreases with increasing accretion rate. The universal AGN power-spectrum is consistent with the mean, 2-9 keV band Cyg X-1 power spectrum in its soft state. The mean, low frequency AGN X-ray power spectrum is consistent with the extension of the mean, 0.01-25 Hz Cyg X-1 power spectrum to lower frequencies. The agreement between the AGN and the Cyg X-1 power spectrum (either in the soft or the hard state) over many decades in frequency indicates that the X-ray variability process is probably the same in all accreting objects, irrespective of the mass of the compact object. We plan to investigate this issue further in the near future. △ Less

Submitted 23 February, 2024; originally announced February 2024.

Comments: 13 pages, 10 figures. Accepted for publication in A&A

Journal ref: A&A 685, A50 (2024)

Comparing indirect methods for black hole masses in AGN: the good, the bad, and the ugly

Authors: M. Gliozzi, J. K. Williams, A. Akylas, I. E. Papadakis, O. I. Shuvo, A. Halavatkar, A. Alt

Abstract: The black hole mass MBH is crucial in constraining the growth of supermassive BHs within their host galaxies. Since direct measurements of MBH with dynamical methods are restricted to a limited number of nearly quiescent nearby galaxies and a small minority of active galactic nuclei (AGN), we must rely on indirect methods. In this work, we utilize an unbiased, volume-limited, hard X-ray selected s… ▽ More The black hole mass MBH is crucial in constraining the growth of supermassive BHs within their host galaxies. Since direct measurements of MBH with dynamical methods are restricted to a limited number of nearly quiescent nearby galaxies and a small minority of active galactic nuclei (AGN), we must rely on indirect methods. In this work, we utilize an unbiased, volume-limited, hard X-ray selected sample of AGN to compare the reliability of some commonly used indirect methods, emphasising those that can be applied to obscured AGN. Based on a subsample of AGN with MBH determined via dynamical methods, our study suggests that X-ray based techniques, such as the scaling method and the one based on the variability measured through the excess variance, are in good agreement with the dynamical methods. On the other hand, the M-sigma correlation based on inactive galaxies tends to systematically overestimate MBH, regardless of the level of obscuration. We provide a correcting factor that produces an acceptable agreement with dynamical values and can be used to quickly correct the MBH computed with this method. We also derive an alternative M-sigma correlation based on this unbiased sample of AGN with a slope considerably shallower than the ones obtained using inactive galaxies, suggesting that the latter correlation may not be appropriate to compute the MBH in AGN. Finally, we find that no quick fix can be applied to correct the MBH obtained from the fundamental plane of black hole activity, casting doubts on the reliability of this method. △ Less

Submitted 21 December, 2023; originally announced December 2023.

Comments: 15 pages, 9 figures, 3 tables, Accepted for publication in MNRAS

Detection of X-ray/UV delay in NGC 4051 using AstroSat observations

Authors: Kavita Kumari, G. C. Dewangan, I. E. Papadakis, K. P. Singh

Abstract: We study the connection between the variations in the far ultra-violet (FUV), near ultra-violet (NUV) and X-ray band emission from NGC 4051 using 4-days long AstroSat observations performed during 5-9 June 2016. NGC 4051 showed rapid variability in all three bands with the strongest variability amplitude in the X-ray band ($F_{var} \sim 37\%$) and much weaker variability in the UV bands (… ▽ More We study the connection between the variations in the far ultra-violet (FUV), near ultra-violet (NUV) and X-ray band emission from NGC 4051 using 4-days long AstroSat observations performed during 5-9 June 2016. NGC 4051 showed rapid variability in all three bands with the strongest variability amplitude in the X-ray band ($F_{var} \sim 37\%$) and much weaker variability in the UV bands ($F_{var} \sim 3 - 5\%$). Cross-correlation analysis performed using Interpolated cross-correlation Functions (ICCF) and Discrete cross-correlation Functions (DCF) revealed a robust correlation ($\sim 0.75$) between the UV and X-ray light curves. The variations in the X-ray band are found to lead those in the FUV and NUV bands by $\sim 7.4{\rm~ks}$ and $\sim 24.2{\rm~ks}$, respectively. The UV lags favour the thermal disc reprocessing model. The FUV and NUV bands are strongly correlated ($\sim 0.9$) and the variations in the FUV band lead those in the NUV band by $\sim 13{\rm~ks}$. Comparison of the UV lags found using the AstroSat observations with those reported earlier and the theoretical model for thermal reverberation time-lag suggests a possible change in either the geometry of the accretion disc/corona or the height of the corona. △ Less

Submitted 5 November, 2023; originally announced November 2023.

Comments: 8 pages, 5 figures, Accepted for publication in MNRAS

Revisiting UV/optical continuum time lags in AGN

Authors: E. S. Kammoun, L. Robin, I. E. Papadakis, M. Dovčiak, C. Panagiotou

Abstract: In this paper, we present an updated version of our model (KYNXiltr) which considers thermal reverberation of a standard Novikov-Thorne accretion disc illuminated by an X-ray point-like source. Previously, the model considered only two cases of black hole spins, and assumed a colour correction factor $f_{\rm col} = 2.4$. Now, we extend the model to any spin value and colour correction. In addition… ▽ More In this paper, we present an updated version of our model (KYNXiltr) which considers thermal reverberation of a standard Novikov-Thorne accretion disc illuminated by an X-ray point-like source. Previously, the model considered only two cases of black hole spins, and assumed a colour correction factor $f_{\rm col} = 2.4$. Now, we extend the model to any spin value and colour correction. In addition, we consider two scenarios of powering the X-ray corona, either via accretion, or external to the accretion disc. We use KYNXiltr to fit the observed time lags obtained from intense monitoring of four local Seyfert galaxies (NGC 5548, NGC 4395, Mrk 817, and Fairall 9). We consider various combinations of black hole spin, colour correction, corona height, and fraction of accretion power transferred to the corona. The model fits well the overall time-lags spectrum in these sources (for a large parameter space). For NGC 4593 only, we detect a significant excess of delays in the U-band. The contribution of the diffuse BLR emission in the time-lags spectrum of this source is significant. It is possible to reduce the large best-fitting parameter space by combining the results with additional information, such as the observed Eddington ratio and average X-ray luminosity. We also provide an update to the analytic expression provided by Kammoun et al., for an X-ray source that is not powered by the accretion process, which can be used for any value of colour correction, and for two values of the black hole spin (0 and 0.998). △ Less

Submitted 13 September, 2023; v1 submitted 11 September, 2023; originally announced September 2023.

Comments: Accepted for publication in MNRAS

Far-Ultraviolet Spectroscopy Of Active Galactic Nuclei With AstroSat/UVIT

Authors: Shrabani Kumar, G. C. Dewangan, K. P. Singh, P. Gandhi, I. E. Papadakis, P. Tripathi, L. Mallick

Abstract: $\require{mediawiki-texvc}$ We study accretion disk emission from eight Seyfert $1 - 1.5$ active galactic nuclei (AGN) using far ultra-violet ($1300-1800$ $Å… ▽ More $\require{mediawiki-texvc}$ We study accretion disk emission from eight Seyfert $1 - 1.5$ active galactic nuclei (AGN) using far ultra-violet ($1300-1800$ $Å$) slit-less grating spectra acquired with AstroSat/UVIT. We correct for the Galactic and intrinsic extinction, contamination from the host galaxies, narrow and broad-line regions, Fe II emission and Balmer continuum, and derive the intrinsic continua. We use HST COS/FOS spectra to account for the emission/absorption lines in the low-resolution UVIT spectra. We find generally redder power-law ($f_ν\propto ν^α$) slopes ($α\sim -1.1 - 0.3$) in the far UV band than predicted by the standard accretion disk model in the optical/UV band. We fit accretion disk models such as the multi-temperature disk blackbody ($\texttt{DISKBB}$) and relativistic disk ($\texttt{ZKERRBB}$, $\texttt{OPTXAGNF}$) models to the observed intrinsic continuum emission. We measure the inner disk temperatures using the $\texttt{DISKBB}$ model for seven AGN. These temperatures in the range $\sim 3.6-5.8$ eV are lower than the peak temperatures predicted for standard disks around maximally spinning super-massive black holes accreting at Eddington rates. The inner disks in two AGN, NGC 7469 and Mrk 352, appear to be truncated at $\sim 35-125r_{g}$ and $50-135r_{g}$, respectively. While our results show that the intrinsic FUV emission from the AGN are consistent with the standard disks, it is possible that UV continua may be affected by the presence of soft X-ray excess emission, X-ray reprocessing, and thermal Comptonisation in the hot corona. Joint spectral modeling of simultaneously acquired UV/X-ray data may be necessary to further investigate the nature of accretion disks in AGN. △ Less

Submitted 21 March, 2023; originally announced March 2023.

Comments: 29 pages, 12 figures, accepted for publication in ApJ

Contrasting X-ray/UV time-lags in Seyfert 1 galaxies NGC 4593 and NGC 7469 using AstroSat observations

Authors: Kavita Kumari, G. C. Dewangan, I. E. Papadakis, Max W. J. Beard, I. M. McHardy, K. P. Singh, D. Bhattacharya, S. Bhattacharyya, S. Chandra

Abstract: We study accretion disk-corona connection in Seyfert 1 galaxies using simultaneous UV/X-ray observations of NGC 4593 (July 14-18, 2016) and NGC 7469 (October 15-19, 2017) performed with AstroSat. We use the X-ray (0.5-7.0 keV) data acquired with the Soft X-ray Telescope (SXT) and the UV (FUV: 130-180 nm, NUV: 200-300 nm) data obtained with the Ultra-Violet Imaging Telescope (UVIT). We also use the… ▽ More We study accretion disk-corona connection in Seyfert 1 galaxies using simultaneous UV/X-ray observations of NGC 4593 (July 14-18, 2016) and NGC 7469 (October 15-19, 2017) performed with AstroSat. We use the X-ray (0.5-7.0 keV) data acquired with the Soft X-ray Telescope (SXT) and the UV (FUV: 130-180 nm, NUV: 200-300 nm) data obtained with the Ultra-Violet Imaging Telescope (UVIT). We also use the contemporaneous Swift observations of NGC 4593 and demonstrate AstroSat's capability for X-ray/UV correlation studies. We performed UV/X-ray cross-correlation analysis using the Interpolated and the Discrete Cross-Correlation Functions and found similar results. In the case of NGC 4593, we found that the variations in the X-rays lead to those in the FUV and NUV bands by ~ 38 ks and ~ 44 ks, respectively. These UV lags favour the disk reprocessing model, they are consistent with the previous results within uncertainties. In contrast, we found an opposite trend in NGC 7469 where the soft X-ray variations lag those in the FUV and NUV bands by ~ 41 ks and ~ 49 ks, respectively. The hard lags in NGC 7469 favour the Thermal Comptonization model. Our results may provide direct observational evidence for the variable intrinsic UV emission from the accretion disk which acts as the seed for thermal Comptonization in a hot corona in a lamp-post like geometry. The non-detection of disk reverberation photons in NGC 7469, using AstroSat data, is most likely due to a high accretion rate resulting in a hot accretion disk and large intrinsic emission. △ Less

Submitted 2 March, 2023; originally announced March 2023.

Comments: 14 pages, 20 figures, Accepted for publication in MNRAS

The universal shape of the X-ray variability power spectrum of AGN up to $z\sim 3$

Authors: M. Paolillo, I. E. Papadakis, W. N. Brandt, F. E. Bauer, G. Lanzuisi, V. Allevato, O. Shemmer, X. C. Zheng, D. De Cicco, R. Gilli, B. Luo, M. Thomas, P. Tozzi, F. Vito, Y. Q. Xue

Abstract: We study the ensemble X-ray variability properties of Active Galactic Nuclei (AGN) over a large range of timescales (20 ks $\leq T\leq$ 14 yrs), redshift ($0\leq z \lesssim 3$), luminosities ($10^{40}$ erg s$^{-1}\leq L_X\leq 10^{46}$ erg s$^{-1}$) and black hole (BH) masses ($10^6 \leq $M$_\odot \leq 10^9$). We propose the use of the variance-frequency diagram, as a viable alternative to the stud… ▽ More We study the ensemble X-ray variability properties of Active Galactic Nuclei (AGN) over a large range of timescales (20 ks $\leq T\leq$ 14 yrs), redshift ($0\leq z \lesssim 3$), luminosities ($10^{40}$ erg s$^{-1}\leq L_X\leq 10^{46}$ erg s$^{-1}$) and black hole (BH) masses ($10^6 \leq $M$_\odot \leq 10^9$). We propose the use of the variance-frequency diagram, as a viable alternative to the study of the power spectral density (PSD), which is not yet accessible for distant, faint and/or sparsely sampled AGN. We show that the data collected from archival observations and previous literature studies are fully consistent with a universal PSD form which does not show any evidence for systematic evolution of shape or amplitude with redshift or luminosity, even if there may be differences between individual AGN at a given redshift or luminosity. We find new evidence that the PSD bend frequency depends on BH mass and, possibly, on accretion rate. We finally discuss the implications for current and future AGN population and cosmological studies. △ Less

Submitted 16 February, 2023; originally announced February 2023.

Comments: In press on Astronomy and Astrophysics

Journal ref: A&A 673, A68 (2023)

A theoretical study of the time-lags due to Comptonization and the constraints on the X-ray corona in AGN

Authors: W. Zhang, I. E. Papadakis, M. Dovčiak, M. Bursa, V. Karas

Abstract: We study the Fourier time-lags due to the Comptonization of disc-emitted photons in a spherical, uniform, and stationary X-ray corona, which located on the rotational axis of the black hole. We use Monk, a general relativistic Monte-Carlo radiative transfer code, to calculate Compton scattering of photons emitted by a thin disc with a Novikov-Thorne temperature profile. We find that the model time… ▽ More We study the Fourier time-lags due to the Comptonization of disc-emitted photons in a spherical, uniform, and stationary X-ray corona, which located on the rotational axis of the black hole. We use Monk, a general relativistic Monte-Carlo radiative transfer code, to calculate Compton scattering of photons emitted by a thin disc with a Novikov-Thorne temperature profile. We find that the model time-lags due to Comptonization remain constant up to a characteristic frequency and then rapidly decrease to zero at higher frequencies. We provide equations which can be used to determine the time-lags and cross spectra for a wide range of values for the corona radius, temperature, optical depth, height, and for various accretion rates and black hole masses. We also provide an equation for the X-ray luminosity of a single corona, as a function of the its characteristics and location above the disc. Remarkably, the observed X-ray time-lags of nearby, bright active galaxies can be successfully reproduced by inverse Comptonization process of multiple dynamic coronae. △ Less

Submitted 5 December, 2022; originally announced December 2022.

Comments: 16 pages, 17 figures; accepted for publication in MNRAS

X-ray illuminated accretion discs and quasar microlensing disc sizes

Authors: I. E. Papadakis, M. Dovciak, E. Kammoun

Abstract: We study the half-light radius versus black hole mass as well as the luminosity versus black hole mass relations in active galactic nuclei (AGN) when the disc is illuminated by the X-ray corona. We used KYNSED, a recently developed spectral model for studying broadband spectral energy distribution in AGN. We considered non-illuminated Novikov-Thorne discs and X-ray illuminated discs based on a Nov… ▽ More We study the half-light radius versus black hole mass as well as the luminosity versus black hole mass relations in active galactic nuclei (AGN) when the disc is illuminated by the X-ray corona. We used KYNSED, a recently developed spectral model for studying broadband spectral energy distribution in AGN. We considered non-illuminated Novikov-Thorne discs and X-ray illuminated discs based on a Novikov-Thorne temperature radial profile. We also considered the case where the temperature profile is modified by a colour-correction factor. We assumed that the X-ray luminosity is equal to the accretion power that is dissipated to the disc below a transition radius. The half-light radius of X-ray illuminated radii can be up to some three to four times greater than the radius of a standard disc, even for a non-spinning black hole, due to the fact that the absorbed X-rays act as a secondary source of energy, increasing the disc temperature. Non-illuminated discs are consistent with observations, but only at the 2.5sigma level. On the other hand, X-ray illuminated discs can explain both the half-light radius-black hole mass as well as the luminosity-black hole mass relation in AGN, for a wide range of physical parameters. In addition, we show that the observed X-ray luminosity of the gravitationally lensed quasars is fully consistent with the X-ray luminosity that is necessary for heating the disc. X-ray disc illumination was proposed many years ago to explain various features that are commonly observed in the X-ray spectra of AGN. Recently, we showed that X-ray illumination of accretion disc can also explain the observed UV/optical time-lags in AGN, while in this work, we show that the same model can also account for the quasar micro-lensing disc size problem. These results support the hypothesis of the disc X-ray illumination in AGN. △ Less

Submitted 14 September, 2022; v1 submitted 25 July, 2022; originally announced July 2022.

Comments: Paper replaced with a new version to correct for minor mistakes (this is the version that will appear in A&A)

Journal ref: A&A 666, A11 (2022)

A physical model for the broadband energy spectrum of X-ray illuminated accretion discs: fitting the spectral energy distribution of NGC 5548

Authors: M. Dovciak, I. E. Papadakis, E. S. Kammoun, W. Zhang

Abstract: We develop a new physical model for the broadband spectral energy distribution (SED) of X-ray illuminated accretion discs, that takes into account the mutual interaction of the accretion disc and the X-ray corona, including all relativistic effects. We assume a Keplerian, optically thick and geometrically thin accretion disc and an X-ray source in the lamp-post geometry that emits an isotropic pow… ▽ More We develop a new physical model for the broadband spectral energy distribution (SED) of X-ray illuminated accretion discs, that takes into account the mutual interaction of the accretion disc and the X-ray corona, including all relativistic effects. We assume a Keplerian, optically thick and geometrically thin accretion disc and an X-ray source in the lamp-post geometry that emits an isotropic power-law spectrum with a high-energy cut-off. We assume that all the energy that would be released by thermal radiation in the standard disc model in its innermost part, is transported to the corona, effectively cooling the disc in this region. We include the disc heating due to thermalisation of the absorbed part of the disc illumination by X-ray corona. The X-ray reflection from the disc is also included. We compute the X-ray luminosity and the low-energy X-ray cut-off through an iterative process, taking full account of the interplay between the X-ray illumination of the disc and the resulting accretion disc spectrum which enters the corona so that the energy balance is preserved. The corona radius is also computed from the conservation of the photon's number during Comptonization. We discuss the model SEDs and their dependence on system parameters. The disc-corona interaction has profound effects - it constrains the X-ray luminosity and changes the shape and normalisation of the UV/optical blue bump. We use the new code to fit the broad-band SED of a typical Seyfert 1 galaxy, NGC 5548. We infer a high black-hole spin, an intermediate system inclination, and an accretion rate below 10% of Eddington. The X-ray luminosity in this source could be supported by 45-70% of the accretion energy dissipated in the disc. The new model, named KYNSED, is publicly available to be used for fitting AGN SEDs inside the XSPEC spectral analysis tool. △ Less

Submitted 21 March, 2022; v1 submitted 4 October, 2021; originally announced October 2021.

Comments: 20 pages, 14 figures, accepted by Astronomy & Astrophysics

Journal ref: A&A 661, A135 (2022)

Revealing Thermal Comptonization of accretion-disk photons in IC4329A with AstroSat

Authors: Prakash Tripathi, Gulab Chand Dewangan, I. E. Papadakis, K. P. Singh

Abstract: We present five simultaneous UV/X-ray observations of IC4329A by AstroSat, performed over {a five-month} period. We utilize the excellent spatial resolution of the Ultra-Violet Imaging Telescope (UVIT) onboard AstroSat to reliably separate the intrinsic AGN flux from the host galaxy emission and to correct for the Galactic and internal reddening, as well as the contribution from the narrow and bro… ▽ More We present five simultaneous UV/X-ray observations of IC4329A by AstroSat, performed over {a five-month} period. We utilize the excellent spatial resolution of the Ultra-Violet Imaging Telescope (UVIT) onboard AstroSat to reliably separate the intrinsic AGN flux from the host galaxy emission and to correct for the Galactic and internal reddening, as well as the contribution from the narrow and broad-line regions. We detect large-amplitude UV variability, which is unusual for a large black hole mass AGN, like IC4329A, over such a small period. In fact, the fractional variability amplitude is larger in the UV band than in the X-ray band. This demonstrates that the observed UV variability is intrinsic to the disk, and is not due to X-ray illumination. The joint X-ray spectral analyses of five SXT and LAXPC spectral data reveal a soft-X-ray excess component, a narrow iron-line (with no indication of a significant Compton hump), and a steepening power-law ($ΔΓ\sim 0.21$) with increasing X-ray flux. The soft excess component could arise due to thermal Comptonization of the inner disk photons in a warm corona with $kT_e\sim 0.26$ keV. The UV emission we detect acts as the primary seed photons for the hot corona, which produces the broadband X-ray continuum. The X-ray spectral variability is well described by the cooling of this corona from $kT_e\sim42$ keV to $\sim 32$ keV with increasing UV flux, while the optical depth remains constant at $τ\sim 2.3$. △ Less

Submitted 27 April, 2021; originally announced April 2021.

Comments: Accepted for publication in ApJ, 18 pages, 9 figures, 4 Tables

AstroSat/UVIT observations of IC4329A: Constraining the accretion disc inner radius

Authors: G. C. Dewangan, P. Tripathi, I. E. Papadakis, K. P. Singh

Abstract: We present a study of far and near-ultraviolet emission from the accretion disk in a powerful Seyfert 1 galaxy IC4329A using observations performed with the Ultraviolet Imaging Telescope (UVIT) onboard AstroSat. These data provide the highest spatial resolution and deepest images of IC4329A in the far and near UV bands acquired to date. The excellent spatial resolution of the UVIT data has allowed… ▽ More We present a study of far and near-ultraviolet emission from the accretion disk in a powerful Seyfert 1 galaxy IC4329A using observations performed with the Ultraviolet Imaging Telescope (UVIT) onboard AstroSat. These data provide the highest spatial resolution and deepest images of IC4329A in the far and near UV bands acquired to date. The excellent spatial resolution of the UVIT data has allowed us to accurately separate the extended emission from the host galaxy and the AGN emission in the far and near UV bands. We derive the intrinsic AGN flux after correcting for the Galactic and internal reddening, as well as for the contribution of emission lines from the broad and narrow-line regions. The intrinsic UV continuum emission shows a marked deficit compared to that expected from the "standard" models of the accretion disk around an estimated black hole mass of 1-2x10^8Msun when the disk extends to the innermost stable circular orbit. We find that the intrinsic UV continuum is fully consistent with the standard disk models, but only if the disk emits from distances larger than 80-150 gravitational radii. △ Less

Submitted 16 April, 2021; originally announced April 2021.

Comments: Accepted for publication in MNRAS, 10 pages, 5 figures

Modeling the UV/optical continuum time-lags in AGN

Authors: E. S. Kammoun, I. E. Papadakis, M. Dovciak

Abstract: Thermal reverberation in accretion discs of active galactic nuclei is thought to be the reason of the continuum UV/optical time lags seen in these sources. Recently, we studied thermal reverberation of a standard Novikov-Thorne accretion disc illuminated by an X-ray point-like source, and we derived an analytic prescription for the time lags as function of wavelength. In this work, we use this ana… ▽ More Thermal reverberation in accretion discs of active galactic nuclei is thought to be the reason of the continuum UV/optical time lags seen in these sources. Recently, we studied thermal reverberation of a standard Novikov-Thorne accretion disc illuminated by an X-ray point-like source, and we derived an analytic prescription for the time lags as function of wavelength. In this work, we use this analytic function to fit the time-lags spectra of seven Seyferts, that have been intensively monitored, in many wave-bands, in the last few years. We find that thermal reverberation can explain the observed UV/optical time lags in all these sources. Contrary to previous claims, the magnitude of the observed UV/optical time-lags is exactly as expected in the case of a standard accretion disc in the lamp-post geometry, given the black hole mass and the accretion rate estimates for the objects we study. We derive estimates of the disc accretion rates and corona height for a non-spinning and a maximally spinning black hole scenarios. We also find that the modelling of the continuum optical/UV time-lags can be used to estimate the black hole spin, when combined with additional information. We also find that the model under-predicts the observed X-ray to UV time-lags, but this difference is probably due to the broad X-ray auto-correlation function of these sources. △ Less

Submitted 8 March, 2021; originally announced March 2021.

Comments: Accepted for publication in MNRAS

UV/optical disk thermal reverberation in AGN: an in-depth study with an analytic prescription for the time-lag spectra

Authors: E. S. Kammoun, M. Dovciak, I. E. Papadakis, M. D. Caballero-Garcia, V. Karas

Abstract: Several active galactic nuclei show correlated variations in the ultraviolet/optical range, with time delays increasing at longer wavelengths. Thermal reprocessing of the X-rays illuminating the accretion disk has been proposed as a viable explanation. In this scenario, the variable X-ray flux irradiating the accretion disk is partially reflected in X-rays, and partially absorbed, thermalized and… ▽ More Several active galactic nuclei show correlated variations in the ultraviolet/optical range, with time delays increasing at longer wavelengths. Thermal reprocessing of the X-rays illuminating the accretion disk has been proposed as a viable explanation. In this scenario, the variable X-ray flux irradiating the accretion disk is partially reflected in X-rays, and partially absorbed, thermalized and re-emitted with some delay by the accretion disk at longer wavelengths. We investigate this scenario assuming an X-ray point-like source illuminating a standard Novikov-Thorne accretion disk, around a rotating black hole. We consider all special and general relativistic effects to determine the incident X-ray flux on the disk and in propagating light from the source to the disk and to the observer. We also compute the disk reflection flux taking into consideration the disk ionization. We investigate the dependence of the disk response function and time lags on various physical parameters, such as black hole mass and spin, X-ray corona height, luminosity, and photon index, accretion rate, inclination, and inner/outer disk radii. We found it is important to consider relativistic effects and the disk ionization in estimating the disk response. We also found a strong non-linearity between the X-ray luminosity and the disk response. We present an analytic function for the time-lags dependence on wavelength, which can be used to fit observed time-lag spectra. We also estimate the fraction of the reverberation signal with respect to the total flux and we suggest possible explanation for the lack of X-ray-ultraviolet/optical correlated variations in a few sources. △ Less

Submitted 17 November, 2020; originally announced November 2020.

Comments: Accepted for publication in ApJ

A combined timing/spectral study of IRAS 13224-3809 using XMM-Newton data

Authors: M. D. Caballero-Garcia, I. E. Papadakis, M. Dovciak, M. Bursa, J. Svoboda, V. Karas

Abstract: We present the results from an X-ray variability study of IRAS 13224-3809. This is probably the best source for X-ray reverberation studies since it is X-ray bright, extremely variable, and it has been extensively observed with XMM-Newton. We used all the archival XMM-Newton data from the three EPIC cameras (to increase the signal-to-noise) and, given the many observations of the source, we were a… ▽ More We present the results from an X-ray variability study of IRAS 13224-3809. This is probably the best source for X-ray reverberation studies since it is X-ray bright, extremely variable, and it has been extensively observed with XMM-Newton. We used all the archival XMM-Newton data from the three EPIC cameras (to increase the signal-to-noise) and, given the many observations of the source, we were able to compute the time-lags spectra in three different flux levels/periods. We fitted the time-lags and energy spectra, simultaneously, using a new X-ray reverberation code which computes the time dependent reflection spectra of the disc as a response to an X-ray flash from a point source located on the axis of the black-hole (BH) accretion disc (lamp-post geometry). To the best of our knowledge, this is the first time for an AGN that both time-lags and energy spectra are fitted by a model simultaneously in different flux periods. The model fits in the case when the BH is rapidly rotating are significantly better than the model fits in the case of a Schwarzschild BH. This result strongly favours the hypothesis of a rotating central BH in this source. We also detect significant variations in the height of the X-ray corona. The X-ray height appears to increase from 3-5 gravitational radii when the X-ray luminosity is of the order of 1.5-3 percent of the Eddington limit, up to 10 gravitational radii, when the luminosity doubles. △ Less

Submitted 18 August, 2020; v1 submitted 1 July, 2020; originally announced July 2020.

Comments: Accepted for publication in MNRAS (18th Aug. 2020)

A hard look at thermal reverberation and optical/UV lags in NGC 5548

Authors: E. S. Kammoun, I. E. Papadakis, M. Dovciak

Abstract: The UV/optical variations in many AGN are very well correlated, showing delays which increase with increasing wavelength. It is thought that this is due to thermal reprocessing of the X-ray emission by the accretion disk. In this scenario, the variable X-ray flux from the corona illuminates the accretion disk where it is partially reflected, and partially absorbed and thermalized in the disk produ… ▽ More The UV/optical variations in many AGN are very well correlated, showing delays which increase with increasing wavelength. It is thought that this is due to thermal reprocessing of the X-ray emission by the accretion disk. In this scenario, the variable X-ray flux from the corona illuminates the accretion disk where it is partially reflected, and partially absorbed and thermalized in the disk producing a UV/optical reverberation signal. This will lead to a time lag increasing with wavelength. However, although the shape of the observed time-lags as a function of wavelength is consistent with the model predictions, their amplitude suggested a disk which is significantly hotter than expected. In this work, we estimate the response functions and the corresponding time lags assuming a standard Novikov-Thorne accretion disk illuminated by a point-like X-ray source. We take into account all relativistic effects in the light propagation from the X-ray source to the disk then to the observer. We also compute the disk reflection, accounting for its ionization profile. Our results show that thermal reverberation effects are stronger in sources with large X-ray source height and low accretion rate. We also found that the time lags increase with height and accretion rate. We apply our model to NGC 5548 and we show that the observed lags in this source can be explained by the model, for a source height of $\sim 60 ~\rm r_g$ and an accretion rate of a few percent of the Eddington limit for a maximally-spinning black hole. △ Less

Submitted 19 June, 2019; v1 submitted 18 June, 2019; originally announced June 2019.

Comments: Accepted for publication in ApJL

The long term X-ray time-lags of NGC 4051

Authors: I. E. Papadakis, A. Rigas, A. Markowitz, I. M. McHardy

Abstract: We present the long term, frequency-dependent, X-ray time lags of the Seyfert galaxy NGC 4051. We used Rossi X-ray Timing Explorer light curves in the 2-4, 4-7 and 7-10 keV bands and we measured the time-lags at $10^{-7}-10^{-6}$ Hz. This is the lowest frequency range that AGN time-lags have ever been measured. The variations in the higher energy bands are delayed with respect to the variations we… ▽ More We present the long term, frequency-dependent, X-ray time lags of the Seyfert galaxy NGC 4051. We used Rossi X-ray Timing Explorer light curves in the 2-4, 4-7 and 7-10 keV bands and we measured the time-lags at $10^{-7}-10^{-6}$ Hz. This is the lowest frequency range that AGN time-lags have ever been measured. The variations in the higher energy bands are delayed with respect to the variations we observe at 2-4 keV, in agreement with the time-lags at high frequencies. When we combine our results with those from the model fitting of the time lags at higher frequencies we find that that the X-ray hard lags in NGC 4051 follow a power law of slope $\sim -1$ over a broad frequency range, from $\sim 10^{-7}$ to $\sim 10^{-3}$ Hz. This is consistent with the time-lags of Cyg X-1, a result which further supports the analogy between active galaxies and Galactic X-ray black hole binaries. △ Less

Submitted 14 March, 2019; v1 submitted 12 March, 2019; originally announced March 2019.

Comments: Accepted for publication in MNRAS

Observatory science with eXTP

Authors: Jean J. M. in 't Zand, Enrico Bozzo, Jinlu Qu, Xiang-Dong Li, Lorenzo Amati, Yang Chen, Immacolata Donnarumma, Victor Doroshenko, Stephen A. Drake, Margarita Hernanz, Peter A. Jenke, Thomas J. Maccarone, Simin Mahmoodifar, Domitilla de Martino, Alessandra De Rosa, Elena M. Rossi, Antonia Rowlinson, Gloria Sala, Giulia Stratta, Thomas M. Tauris, Joern Wilms, Xuefeng Wu, Ping Zhou, Iván Agudo, Diego Altamirano , et al. (159 additional authors not shown)

Abstract: In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to stu… ▽ More In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to study one common aspect of these objects: their often transient nature. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s. △ Less

Submitted 10 December, 2018; originally announced December 2018.

Comments: Accepted for publication on Sci. China Phys. Mech. Astron. (2019)

The nature of X-ray spectral variability in SWIFT J2127.4+5654

Authors: E. S. Kammoun, I. E. Papadakis

Abstract: We study the flux-flux plots (FFPs) of the Seyfert galaxy SWIFT J2127.4+5654 which was observed simultaneously by XMM-Newton and NuSTAR. The 0.7-2 keV FFPs show a non-linear behaviour, while they are well fitted by a straight line, in the $\sim 2-40$ keV range. Without any additional modelling, this result strongly suggests that neither absorption nor spectral slope variations can contribute signi… ▽ More We study the flux-flux plots (FFPs) of the Seyfert galaxy SWIFT J2127.4+5654 which was observed simultaneously by XMM-Newton and NuSTAR. The 0.7-2 keV FFPs show a non-linear behaviour, while they are well fitted by a straight line, in the $\sim 2-40$ keV range. Without any additional modelling, this result strongly suggests that neither absorption nor spectral slope variations can contribute significantly to the observed variability above 2 keV (on time scales $\lesssim$5 days, at least). The FFP method is capable of separating the variable and the constant X-ray spectral components in this source. We found that the low-, average- and high-flux spectra of the variable components are consistent with a power law varying in normalisation only plus a relativistic reflection component varying simultaneously with the continuum. At low energies, we identify an additional variable component, consistent with a neutral absorber with a variable covering fraction. We also detect two spectral components whose fluxes remain constant over the duration of the observations. These components can be well described by a blackbody emission of $kT_{\rm BB}\sim 0.18\, \rm keV$, dominating at soft X-rays, and a reflection component from neutral material. The spectral model resulting from the FFP analysis can fit well the energy spectra of the source. We suggest that a combined study of FFPs and energy spectra in the case of X-ray bright and variable Seyferts can identify, unambiguously, the various spectral components, and can lift the degeneracy that is often observed within different models fitting equally well time-averaged spectra of these objects. △ Less

Submitted 30 July, 2018; originally announced July 2018.

Comments: Accepted for publication in MNRAS

Testing the X-ray reverberation model KYNREFREV in a sample of Seyfert 1 Active Galactic Nuclei

Authors: M. D. Caballero-Garcia, I. E. Papadakis, M. Dovciak, M. Bursa, A. Epitropakis, V. Karas, J. Svoboda

Abstract: We present the first results obtained by the application of the KYNREFREV-reverberation model, which is ready for its use in XSPEC. This model computes the time dependent reflection spectra of the disc as a response to a flash of primary power-law radiation from a point source corona located on the axis of the black hole accretion disc (lamp-post geometry). Full relativistic effects are taken into… ▽ More We present the first results obtained by the application of the KYNREFREV-reverberation model, which is ready for its use in XSPEC. This model computes the time dependent reflection spectra of the disc as a response to a flash of primary power-law radiation from a point source corona located on the axis of the black hole accretion disc (lamp-post geometry). Full relativistic effects are taken into account. The ionisation of the disc is set for each radius according to the amount of the incident primary flux and the density of the accretion disc. We tested the model by fitting model predictions to the observed time-lag spectra of three Narrow-Line Seyfert 1 galaxies (ARK 564, MCG-6-30-15 and 1H 0707-495), assuming either a rapidly or zero spinning black hole (BH). The time-lags strongly suggest a compact X-ray source, located close to the BH, at a height of approx. 4 gravitational radii. This result does not depend either on the BH spin or the disc ionization. There is no significant statistical difference between the quality of the best-fits in the rapidly and zero spinning BH scenarios in Ark 564 and MCG-6-30-15. But there is an indication that the hypothesis of a non-rotating BH in 1H 0707-495 is not consistent with its time-lag spectrum. Finally, the best-fits to the Ark 564 and 1H 0707-495 data are of rather low quality. We detect wavy-residuals around the best-fit reverberation model time-lags at high frequencies. This result suggests that the simple lamp-post geometry does not fully explain the X-ray source/disc configuration in Active Galactic Nuclei. △ Less

Submitted 23 July, 2018; v1 submitted 10 April, 2018; originally announced April 2018.

Comments: Accepted in MNRAS after minor revision (23/Jul./2018). 10 pages

X-ray/UV/optical variability of NGC 4593 with Swift: Reprocessing of X-rays by an extended reprocessor

Authors: I M McHardy, S D Connolly, K Horne E M Cackett, J Gelbord, B M Peterson, M Pahari, N Gehrels, R Edelson, M Goad, P Lira, P Arevalo, R D Baldi, N Brandt, E Breedt, H Chand, G Dewangan, C Done, M Elvis, D Emmanoulopoulos, M M Fausnaugh, S Kaspi, C S Kochanek, K Korista, I E Papadakis, A R Rao , et al. (3 additional authors not shown)

Abstract: We report the results of intensive X-ray, UV and optical monitoring of the Seyfert 1 galaxy NGC 4593 with Swift. There is no intrinsic flux-related spectral change in the the variable components in any band with small apparent variations due only to contamination by a second constant component, possibly a (hard) reflection component in the X-rays and the (red) host galaxy in the UV/optical bands.… ▽ More We report the results of intensive X-ray, UV and optical monitoring of the Seyfert 1 galaxy NGC 4593 with Swift. There is no intrinsic flux-related spectral change in the the variable components in any band with small apparent variations due only to contamination by a second constant component, possibly a (hard) reflection component in the X-rays and the (red) host galaxy in the UV/optical bands. Relative to the shortest wavelength band, UVW2, the lags of the other UV and optical bands are mostly in agreement with the predictions of reprocessing of high energy emission from an accretion disc. The U-band lag is, however, far larger than expected, almost certainly because of reprocessed Balmer continuum emission from the more distant broad line region gas. The UVW2 band is well correlated with the X-rays but lags by ~6x more than expected if the UVW2 results from reprocessing of X-rays on the accretion disc. However, if the lightcurves are filtered to remove variations on timescales >5d, the lag approaches the expectation from disc reprocessing. MEMEcho analysis shows that direct X-rays can be the driver of most of the variations in the UV/optical bands as long as the response functions for those bands all have long tails (up to 10d) in addition to a strong peak (from disc reprocessing) at short lag (<1d). We interpret the tails as due to reprocessing from the surrounding gas. Comparison of X-ray to UVW2 and UVW2 to V-band lags for 4 AGN, including NGC 4593, shows that all have UVW2 to V-band lags which exceed the expectations from disc resprocessing by factor < 2. However the X-ray to UVW2 lags are, mostly, in greater excess from the expectations from disc reprocessing and differ between AGN. The largest excess is in NGC 4151. Absorption and scattering may be affecting X-ray to UV lags. △ Less

Submitted 13 December, 2017; originally announced December 2017.

Comments: 18 pages, 17 figures. Submitted to Monthly Notices of the Royal Astronomical Society

The photon-index - time-lag correlation in black-hole X-ray binaries

Authors: Pablo Reig, Nikolaos D. Kylafis, Iossif E. Papadakis, Maria Teresa Costado

Abstract: We have performed a timing and spectral analysis of a set of black-hole binaries to study the correlation between the photon index and the time lag of the hard photons with respect to the soft ones. We provide further evidence that the timing and spectral properties in black-hole X-ray binaries are coupled. In particular, we find that the average time lag increases as the X-ray emission becomes so… ▽ More We have performed a timing and spectral analysis of a set of black-hole binaries to study the correlation between the photon index and the time lag of the hard photons with respect to the soft ones. We provide further evidence that the timing and spectral properties in black-hole X-ray binaries are coupled. In particular, we find that the average time lag increases as the X-ray emission becomes softer. Although a correlation between the hardness of the X-ray spectrum and the time (or phase) lag has been reported in the past for a handful of systems, our study confirms that this correlated behaviour is a global property of black-hole X-ray binaries. We also demonstrate that the photon-index - time-lag correlation can be explained as a result of inverse Comptonization in a jet. △ Less

Submitted 13 October, 2017; originally announced October 2017.

Comments: 10 pages, 7 figures, accepted for publications in MNRAS

The nature of X-ray spectral variability in MCG-6-30-15

Authors: E. S. Kammoun, I. E. Papadakis

Abstract: The flux-flux plot (FFP) method can provide model-independent clues regarding the X-ray variability of active galactic nuclei. To use it properly, the bin size of the light curves should be as short as possible, provided the average counts in the light curve bins are larger than $\sim 200$. We apply the FFP method to the 2013, simultaneous XMM-Newton and NuSTAR observations of the Seyfert galaxy M… ▽ More The flux-flux plot (FFP) method can provide model-independent clues regarding the X-ray variability of active galactic nuclei. To use it properly, the bin size of the light curves should be as short as possible, provided the average counts in the light curve bins are larger than $\sim 200$. We apply the FFP method to the 2013, simultaneous XMM-Newton and NuSTAR observations of the Seyfert galaxy MCG$-$6-30-15, in the 0.3-40 keV range. The FFPs above $\sim 1.6$ keV are well-described by a straight line. This result rules out spectral slope variations and the hypothesis of absorption driven variability. Our results are fully consistent with a power-law component varying in normalization only, with a spectral slope of $\sim 2$, plus a variable, relativistic reflection arising from the inner accretion disc around a rotating black hole. We also detect spectral components which remain constant over $\sim 4.5$ days (at least). At energies above $\sim 1.5$ keV, the stable component is consistent with reflection from distant, neutral material. The constant component at low energies is consistent with a blackbody spectrum of $kT_{\rm BB} \sim 100$ eV. The fluxes of these components are $\sim 10-20\%$ of the average continuum flux (in the respective bands). They should always be included in the models that are used to fit the spectrum of the source. The FFPs below 1.6 keV are non-linear, which could be due to the variable warm absorber in this source. △ Less

Submitted 14 September, 2017; v1 submitted 4 September, 2017; originally announced September 2017.

Comments: 18 pages, 12 figures, 6 tables. Accepted for publication in MNRAS

Complex UV/X-ray variability of 1H~0707-495

Authors: P. K. Pawar, G. C. Dewangan, I. E. Papadakis, M. K. Patil, Main Pal, A. K. Kembhavi

Abstract: We study the relationship between the UV and X-ray variability of the narrow-line Seyfert 1 galaxy 1H~0707-495. Using a year long {\it Swift} monitoring and four long {\it XMM-Newton} observations, we perform cross-correlation analyses of the UV and X-ray light curves, on both long and short time scales. We also perform time-resolved X-ray spectroscopy on 1-2 ks scale, and study the relationship b… ▽ More We study the relationship between the UV and X-ray variability of the narrow-line Seyfert 1 galaxy 1H~0707-495. Using a year long {\it Swift} monitoring and four long {\it XMM-Newton} observations, we perform cross-correlation analyses of the UV and X-ray light curves, on both long and short time scales. We also perform time-resolved X-ray spectroscopy on 1-2 ks scale, and study the relationship between the UV emission and the X-ray spectral components - soft X-ray excess and a power-law. We find that the UV and X-ray variations anti-correlate on short, and possibly on long time scales as well. Our results rule out reprocessing as the dominant mechanism for the UV variability, as well as the inward propagating fluctuations in the accretion rate. Absence of a positive correlation between the photon index and the UV flux suggests that the observed UV emission is unlikely to be the seed photons for the thermal Comptonisation. We find a strong correlation between the continuum flux and the soft-excess temperature which implies that the soft excess is most likely the reprocessed X-ray emission in the inner accretion disc. Strong X-ray heating of the innermost regions in the disc, due to gravitational light bending, appears to be an important effect in 1H~0707-495, giving rise to a significant fraction of the soft excess as reprocessed thermal emission. We also find indications for a non-static, dynamic X-ray corona, where either the size or height (or both) vary with time. △ Less

Submitted 29 August, 2017; originally announced August 2017.

Comments: 10 pages, 7 figures, Accepted to MNRAS

The X-ray continuum time-lags and intrinsic coherence in AGN

Authors: A. Epitropakis, I. E. Papadakis

Abstract: We present the results from a systematic analysis of the X-ray continuum (`hard') time-lags and intrinsic coherence between the $2-4\,\mathrm{keV}$ and various energy bands in the $0.3-10\,\mathrm{keV}$ range, for ten X-ray bright and highly variable active galactic nuclei (AGN). We used all available archival \textit{XMM-Newton} data, and estimated the time-lags following Epitropakis \& Papadakis… ▽ More We present the results from a systematic analysis of the X-ray continuum (`hard') time-lags and intrinsic coherence between the $2-4\,\mathrm{keV}$ and various energy bands in the $0.3-10\,\mathrm{keV}$ range, for ten X-ray bright and highly variable active galactic nuclei (AGN). We used all available archival \textit{XMM-Newton} data, and estimated the time-lags following Epitropakis \& Papadakis (2016). By performing extensive numerical simulations, we arrived at useful guidelines for computing intrinsic coherence estimates that are minimally biased, have known errors, and are (approximately) Gaussian distributed. Owing to the way we estimated the time-lags and intrinsic coherence, we were able to do a proper model fitting to the data. Regarding the continuum time-lags, we are able to demonstrate that they have a power-law dependence on frequency, with a slope of $-1$, and that their amplitude scales with the logarithm of the light-curve mean-energy ratio. We also find that their amplitude increases with the square root of the X-ray Eddington ratio. Regarding the intrinsic coherence, we found that it is approximately constant at low frequencies. It then decreases exponentially at frequencies higher than a characteristic `break frequency.' Both the low-frequency constant intrinsic-coherence value and the break frequency have a logarithmic dependence on the light-curve mean-energy ratio. Neither the low-frequency constant intrinsic-coherence value, nor the break frequency exhibit a universal scaling with either the central black hole mass, or the the X-ray Eddington ratio. Our results could constrain various theoretical models of AGN X-ray variability. △ Less

Submitted 13 March, 2017; originally announced March 2017.

Comments: 36 pages, 6 tables, 42 figures (accepted for publication in MNRAS)

Long-term monitoring of Ark 120 with Swift

Authors: M. Gliozzi, I. E. Papadakis, D. Grupe, W. P. Brinkmann, C. Raeth

Abstract: We report the results of a six-month Swift monitoring campaign of Ark120, a prototypical bare Seyfert1 galaxy. The lack of intrinsic absorption combined with the nearly contemporaneous coverage of the UV and X-ray bands makes it possible to investigate the link between the accretion disk and the Comptonization corona. Our observations confirm the presence of substantial temporal variability, with… ▽ More We report the results of a six-month Swift monitoring campaign of Ark120, a prototypical bare Seyfert1 galaxy. The lack of intrinsic absorption combined with the nearly contemporaneous coverage of the UV and X-ray bands makes it possible to investigate the link between the accretion disk and the Comptonization corona. Our observations confirm the presence of substantial temporal variability, with the X-rays characterized by large-amplitude flux changes on timescales of few days, while the variations in the UV bands are smoother on timescales of several weeks. The source also shows spectral variability with the X-ray spectrum steepening when the source is brighter. We do not detect any correlation between the UV flux and the X-ray spectral slope. A cross correlation analysis suggests positive delays between X-rays and the UV emission, favoring a scenario of disk reprocessing. Although the strength of the correlation is moderate with a delay not well constrained (7.5+-7 days), it is nevertheless indicative of a large disk reprocessing region, with a separation between the X-ray and the UV emitting regions of the order of 1000 rG. The Ark120 correlation results are in agreement with those obtained in similar monitoring studies. When combined together, the observations can be well described by a linear relation between the X-ray/UV delay and the black hole mass. Within the context of the simplest scenario where these delays correspond to light-travel times, the implied distance between the X-ray source and the UV disk reprocessing region is of the order of many hundreds of gravitational radii. △ Less

Submitted 13 October, 2016; originally announced October 2016.

Comments: 11 pages, 11 figures, 4 tables, accepted for publication in MNRAS

Theoretical modelling of the AGN iron-line vs continuum time-lags in the lamp-post geometry

Authors: A. Epitropakis, I. E. Papadakis, M. Dovčiak, T. Pecháček, D. Emmanoulopoulos, V. Karas, I. M. McHardy

Abstract: Context: Theoretical modelling of time-lags between variations in the Fe K$α$ emission and the X-ray continuum might shed light on the physics and geometry of the X-ray emitting region in active galaxies (AGN) and X-ray binaries. We here present the results from a systematic analysis of time-lags between variations in two energy bands ($5-7$ vs $2-4\,\mathrm{keV}$) for seven X-ray bright and varia… ▽ More Context: Theoretical modelling of time-lags between variations in the Fe K$α$ emission and the X-ray continuum might shed light on the physics and geometry of the X-ray emitting region in active galaxies (AGN) and X-ray binaries. We here present the results from a systematic analysis of time-lags between variations in two energy bands ($5-7$ vs $2-4\,\mathrm{keV}$) for seven X-ray bright and variable AGN. Aims: We estimate time-lags as accurately as possible and fit them with theoretical models in the context of the lamp-post geometry. We also constrain the geometry of the X-ray emitting region in AGN. Methods: We used all available archival \textit{XMM-Newton} data for the sources in our sample and extracted light curves in the $5-7$ and $2-4\,\mathrm{keV}$ energy bands. We used these light curves and applied a thoroughly tested (through extensive numerical simulations) recipe to estimate time-lags that have minimal bias, approximately follow a Gaussian distribution, and have known errors. Using traditional $χ^2$ minimisation techniques, we then fitted the observed time-lags with two different models: a phenomenological model where the time-lags have a power-law dependence on frequency, and a physical model, using the reverberation time-lags expected in the lamp-post geometry. The latter were computed assuming a point-like primary X-ray source above a black hole surrounded by a neutral and prograde accretion disc with solar iron abundance. We took all relativistic effects into account for various X-ray source heights, inclination angles, and black hole spin values. △ Less

Submitted 4 August, 2016; v1 submitted 9 July, 2016; originally announced July 2016.

Comments: Accepted for publication in A&A

Journal ref: A&A 594, A71 (2016)

A search for X-ray reprocessing echoes in the power spectral density functions of AGN

Authors: D. Emmanoulopoulos, I. E. Papadakis, A. Epitropakis, T. Pecháček, M. Dovčiak, I. M. McHardy

Abstract: We present the results of a detailed study of the X-ray power spectra density (PSD) functions of twelve X-ray bright AGN, using almost all the archival XMM-Newton data. The total net exposure of the EPIC-pn light curves is larger than 350 ks in all cases (and exceeds 1 Ms in the case of 1H 0707-497). In a physical scenario in which X-ray reflection occurs in the inner part of the accretion disc of… ▽ More We present the results of a detailed study of the X-ray power spectra density (PSD) functions of twelve X-ray bright AGN, using almost all the archival XMM-Newton data. The total net exposure of the EPIC-pn light curves is larger than 350 ks in all cases (and exceeds 1 Ms in the case of 1H 0707-497). In a physical scenario in which X-ray reflection occurs in the inner part of the accretion disc of AGN, the X-ray reflection component should be a filtered echo of the X-ray continuum signal and should be equal to the convolution of the primary emission with the response function of the disc. Our primary objective is to search for these reflection features in the 5-7 keV (iron line) and 0.5-1 keV (soft) bands, where the X-ray reflection fraction is expected to be dominant. We fit to the observed periodograms two models: a simple bending power law model (BPL) and a BPL model convolved with the transfer function of the accretion disc assuming the lamp-post geometry and X-ray reflection from a homogeneous disc. We do not find any significant features in the best-fitting BPL model residuals either in individual PSDs in the iron band, soft and full band (0.3-10 keV) or in the average PSD residuals of the brightest and more variable sources (with similar black hole mass estimates). The typical amplitude of the soft and full-band residuals is around 3-5 per cent. It is possible that the expected general relativistic effects are not detected because they are intrinsically lower than the uncertainty of the current PSDs, even in the strong relativistic case in which X-ray reflection occurs on a disc around a fast rotating black hole having an X-ray source very close above it. However, we could place strong constrains to the X-ray reflection geometry with the current data sets if we knew in advance the intrinsic shape of the X-ray PSDs, particularly its high frequency slope. △ Less

Submitted 2 June, 2016; originally announced June 2016.

Comments: Accepted for publication in MNRAS. The paper is 15 pages long and contains 7 figures and 6 tables

MSC Class: 85-02 ACM Class: J.2

Modelling the variable broad-band optical/UV/X-ray spectrum of PG1211+143: Implications for the ionized outflow

Authors: I. E. Papadakis, F. Nicastro, C. Panagiotou

Abstract: We present the results from a detailed analysis of the 2007 Swift monitoring campaign of the quasar PG1211+143. We constructed broad-band, optical/UV/X-ray spectral energy distributions over three X-ray flux intervals. We fitted them with a model which accounts for the disc and the X-ray coronal emission and the warm absorber (well established in this source). The three flux spectra are well fitte… ▽ More We present the results from a detailed analysis of the 2007 Swift monitoring campaign of the quasar PG1211+143. We constructed broad-band, optical/UV/X-ray spectral energy distributions over three X-ray flux intervals. We fitted them with a model which accounts for the disc and the X-ray coronal emission and the warm absorber (well established in this source). The three flux spectra are well fitted by the model we considered. The disc inner temperature remains constant at ~2 eV, while X-rays are variable both in spectral slope and normalization. The absorber covers almost 90% of the central source. It is outflowing with a velocity less than 2.3*10^4 km/s (3sigma upper limit), and has a column density of ~10^23.2. Its ionization parameter varies by a factor of 1.6, and it is in photo-ionizing equilibrium with the ionizing flux. It is located at a distance of less than 0.35 pc from the central source and its relative thickness, DR/R is less than 0.1. The absorber' s ionization parameter variations can explain the larger than average amplitude of the X-ray variations. The absence of optical/UV variations (consistent with the high black hole mass estimate) argues against the presence of inward propagating disc fluctuations and strong X-ray illumination of the disc (in agreement with the low ratio of X-ray over the bolometric luminosity of ~20-35). We estimate an upper limit for the mass outflow of ~5 solar masses per year (~2.3 times the Eddington mass accretion rate). If the outflow rate is indeed that high, then it must be a short-lived episode in the quasar's life time. Finally, we estimate an upper limit for the kinetic power of the outflow of ~1.4*10^43 ergs/s. This outflow cannot deploy significant mechanical energy to the surrounding ISM of the quasar's host galaxy, but is sufficient to heat the ISM to 10^7 K and to produce a fast decline to the star formation rate of the galaxy. △ Less

Submitted 18 May, 2016; originally announced May 2016.

Comments: Accepted for publication by A&A

Journal ref: A&A 591, A102 (2016)

Extensive X-ray variability studies of NGC 7314 using long XMM-Newton observations

Authors: D. Emmanoulopoulos, I. M. McHardy, S. Vaughan, I. E. Papadakis

Abstract: We present a detailed X-ray variability study of the low mass Active Galactic Nuclei (AGN) NGC 7314 using the two newly obtained XMM-Newton observations ($140$ and $130$ ks), together with two archival data sets of shorter duration ($45$ and $84$ ks). The relationship between the X-ray variability characteristics and other physical source properties (such as the black hole mass) are still relative… ▽ More We present a detailed X-ray variability study of the low mass Active Galactic Nuclei (AGN) NGC 7314 using the two newly obtained XMM-Newton observations ($140$ and $130$ ks), together with two archival data sets of shorter duration ($45$ and $84$ ks). The relationship between the X-ray variability characteristics and other physical source properties (such as the black hole mass) are still relatively poorly defined, especially for low-mass AGN. We perform a new, fully analytical, power spectral density (PSD) model analysis method, which will be described in detail in a forthcoming paper, that takes into consideration the spectral distortions, caused by red-noise leak. We find that the PSD in the $0.5-10$ keV energy range, can be represented by a bending power-law with a bend around $6.7\times10^{-5}$ Hz, having a slope of $0.51$ and $1.99$ below and above the bend, respectively. Adding our bend time-scale estimate, to an already published ensemble of estimates from several AGN, supports the idea that the bend time-scale depends linearly only on the black hole mass and not on the bolometric luminosity. Moreover, we find that as the energy range increases, the PSD normalization increases and there is a hint that simultaneously the high frequency slope becomes steeper. Finally, the X-ray time-lag spectrum of NGC 7314 shows some very weak signatures of relativistic reflection, and the energy resolved time-lag spectrum, for frequencies around $3\times10^{-4}$ Hz, shows no signatures of X-ray reverberation. We show that the previous claim about ks time-delays in this source, is simply an artefact induced by the minuscule number of points entering during the time-lag estimation in the low frequency part of the time-lag spectrum (i.e. below $10^{-4}$ Hz). △ Less

Submitted 10 May, 2016; originally announced May 2016.

Comments: Accepted for publication in MNRAS. The paper is 21 pages long and contains 15 figures and 3 tables

MSC Class: 85-02 ACM Class: J.2

Statistical properties of Fourier-based time-lag estimates

Authors: A. Epitropakis, I. E. Papadakis

Abstract: The study of X-ray time-lag spectra in active galactic nuclei (AGN) is currently an active research area, since it has the potential to illuminate the physics and geometry of the innermost region (i.e. close to the putative super-massive black hole) in these objects. To obtain reliable information from these studies, the statistical properties of time-lags estimated from data must be known as accu… ▽ More The study of X-ray time-lag spectra in active galactic nuclei (AGN) is currently an active research area, since it has the potential to illuminate the physics and geometry of the innermost region (i.e. close to the putative super-massive black hole) in these objects. To obtain reliable information from these studies, the statistical properties of time-lags estimated from data must be known as accurately as possible. Aims: We investigated the statistical properties of Fourier-based time-lag estimates (i.e. based on the cross-periodogram), using evenly sampled time series with no missing points. Our aim is to provide practical `guidelines' on estimating time-lags that are minimally biased (i.e. whose mean is close to their intrinsic value) and have known errors.} Methods: Our investigation is based on both analytical work and extensive numerical simulations. The latter consisted of generating artificial time series with various signal-to-noise ratios and sampling patterns/durations similar to those offered by AGN observations with present and past X-ray satellites. We also considered a range of different model time-lag spectra commonly assumed in X-ray analyses of compact accreting systems. Results: Discrete sampling, binning and finite light curve duration cause the mean of the time-lag estimates to have a smaller magnitude than their intrinsic values. Smoothing (i.e. binning over consecutive frequencies) of the cross-periodogram can add extra bias at low frequencies. The use of light curves with low signal-to-noise ratio reduces the intrinsic coherence, and can introduce a bias to the sample coherence, time-lag estimates, and their predicted error. △ Less

Submitted 23 April, 2016; v1 submitted 6 April, 2016; originally announced April 2016.

Comments: Accepted for publication in A&A

Journal ref: A&A 591, A113 (2016)

A comprehensive analysis of the hard X-ray spectra of bright Seyfert galaxies

Authors: P. Lubinski, V. Beckmann, L. Gibaud, S. Paltani, I. E. Papadakis, C. Ricci, S. Soldi, M. Türler, R. Walter, A. A. Zdziarski

Abstract: Hard X-ray spectra of 28 bright Seyfert galaxies observed with INTEGRAL were analyzed together with the X-ray spectra from XMM-Newton, Suzaku and RXTE. These broad-band data were fitted with a model assuming a thermal Comptonization as a primary continuum component. We tested several model options through a fitting of the Comptonized continuum accompanied by a complex absorption and a Compton refl… ▽ More Hard X-ray spectra of 28 bright Seyfert galaxies observed with INTEGRAL were analyzed together with the X-ray spectra from XMM-Newton, Suzaku and RXTE. These broad-band data were fitted with a model assuming a thermal Comptonization as a primary continuum component. We tested several model options through a fitting of the Comptonized continuum accompanied by a complex absorption and a Compton reflection. Both the large data set used and the model space explored allowed us to accurately determine a mean temperature kTe of the electron plasma, the Compton parameter y and the Compton reflection strength R for the majority of objects in the sample. Our main finding is that a vast majority of the sample (20 objects) is characterized by kTe < 100 keV, and only for two objects we found kTe > 200 keV. The median kTe for entire sample is 48(-14,+57) keV. The distribution of the y parameter is bimodal, with a broad component centered at ~0.8 and a narrow peak at ~1.1. A complex, dual absorber model improved the fit for all data sets, compared to a simple absorption model, reducing the fitted strength of Compton reflection by a factor of about 2. Modest reflection (median R ~0.32) together with a high ratio of Comptonized to seed photon fluxes point towards a geometry with a compact hard X-ray emitting region well separated from the accretion disc. Our results imply that the template Seyferts spectra used in AGN population synthesis models should be revised. △ Less

Submitted 26 February, 2016; originally announced February 2016.

Comments: 26 pages, 12 figures, accepted for publication in MNRAS

RoboPol: optical polarization-plane rotations and flaring activity in blazars

Authors: D. Blinov, V. Pavlidou, I. E. Papadakis, T. Hovatta, T. J. Pearson, I. Liodakis, G. V. Panopoulou, E. Angelakis, M. Baloković, H. Das, P. Khodade, S. Kiehlmann, O. G. King, A. Kus, N. Kylafis, A. Mahabal, A. Marecki, D. Modi, I. Myserlis, E. Paleologou, I. Papamastorakis, B. Pazderska, E. Pazderski, C. Rajarshi, A. Ramaprakash , et al. (4 additional authors not shown)

Abstract: We present measurements of rotations of the optical polarization of blazars during the second year of operation of RoboPol, a monitoring programme of an unbiased sample of gamma-ray bright blazars specially designed for effective detection of such events, and we analyse the large set of rotation events discovered in two years of observation. We investigate patterns of variability in the polarizati… ▽ More We present measurements of rotations of the optical polarization of blazars during the second year of operation of RoboPol, a monitoring programme of an unbiased sample of gamma-ray bright blazars specially designed for effective detection of such events, and we analyse the large set of rotation events discovered in two years of observation. We investigate patterns of variability in the polarization parameters and total flux density during the rotation events and compare them to the behaviour in a non-rotating state. We have searched for possible correlations between average parameters of the polarization-plane rotations and average parameters of polarization, with the following results: (1) there is no statistical association of the rotations with contemporaneous optical flares; (2) the average fractional polarization during the rotations tends to be lower than that in a non-rotating state; (3) the average fractional polarization during rotations is correlated with the rotation rate of the polarization plane in the jet rest frame; (4) it is likely that distributions of amplitudes and durations of the rotations have physical upper bounds, so arbitrarily long rotations are not realised in nature. △ Less

Submitted 15 January, 2016; v1 submitted 13 January, 2016; originally announced January 2016.

Comments: 12 pages, 12 figures, accepted to MNRAS

Flux and spectral variability of the blazar PKS 2155-304 with XMM-Newton: Evidence of Particle Acceleration and Synchrotron Cooling

Authors: Jai Bhagwana, A. C. Gupta, I. E. Papadakis, Paul J. Wiita

Abstract: We have analyzed XMM-Newton observations of the high energy peaked blazar, PKS 2155-304, made on 24 May 2002 in the 0.3 - 10 keV X-ray band. These observations display a mini-flare, a nearly constant flux period and a strong flux increase. We performed a time-resolved spectral study of the data, by dividing the data into eight segments. We fitted the data with a power-law and a broken power-law mo… ▽ More We have analyzed XMM-Newton observations of the high energy peaked blazar, PKS 2155-304, made on 24 May 2002 in the 0.3 - 10 keV X-ray band. These observations display a mini-flare, a nearly constant flux period and a strong flux increase. We performed a time-resolved spectral study of the data, by dividing the data into eight segments. We fitted the data with a power-law and a broken power-law model, and in some of the segments we found a noticeable spectral flattening of the source's spectrum below 10 keV. We also performed time-resolved cross-correlation analyses and detected significant hard and soft lags (for the first time in a single observation of this source) during the first and last parts of the observation, respectively. Our analysis of the spectra, the variations of photon-index with flux as well as the correlation and lags between the harder and softer X-ray bands indicate that both the particle acceleration and synchrotron cooling processes make an important contribution to the emission from this blazar. The hard lags indicate a variable acceleration process. We also estimated the magnetic field value using the soft lags. The value of the magnetic field is consistent with the values derived from the broad-band SED modeling of this source. △ Less

Submitted 7 October, 2015; originally announced October 2015.

Comments: 30 pages, 4 figures, 1 tables. Accepted for publication in New Astropnomy (in press)

Journal ref: New Astronomy, 44, 21-28 (2016)

Variability of the soft X-ray excess in IRAS 13224-3809

Authors: E. S. Kammoun, I. E. Papadakis, B. M. Sabra

Abstract: We study the soft excess variability of the narrow line Seyfert 1 galaxy IRAS 13224-3809. We considered all five archival XMM-Newton observations, and we applied the 'flux-flux plot' (FFP) method. We found that the flux-flux plots were highly affected by the choice of the light curves' time bin size, most probably because of the fast and large amplitude variations, and the intrinsic non-linear flu… ▽ More We study the soft excess variability of the narrow line Seyfert 1 galaxy IRAS 13224-3809. We considered all five archival XMM-Newton observations, and we applied the 'flux-flux plot' (FFP) method. We found that the flux-flux plots were highly affected by the choice of the light curves' time bin size, most probably because of the fast and large amplitude variations, and the intrinsic non-linear flux--flux relations in this source. Therefore, we recommend that the smallest bin-size should be used in such cases. Hence, We constructed FFPs in 11 energy bands below 1.7 keV, and we considered the 1.7-3 keV band, as being representative of the primary emission. The FFPs are reasonably well fitted by a 'power-law plus a constant' model. We detected significant positive constants in three out of five observations. The best-fit slopes are flatter than unity at energies below $\sim 0.9$ keV, where the soft excess is strongest. This suggests the presence of intrinsic spectral variability. A power-law-like primary component, which is variable in flux and spectral slope (as $Γ\propto N_{\rm PL}^{0.1}$) and a soft-excess component, which varies with the primary continuum (as $F_{\rm excess}\propto F_{\rm primary}^{0.46}$), can broadly explain the FFPs. In fact, this can create positive `constants', even when a stable spectral component does not exist. Nevertheless, the possibility of a stable, soft--band constant component cannot be ruled out, but its contribution to the observed 0.2-1 keV band flux should be less than $\sim 15$ %. The model constants in the FFPs were consistent with zero in one observation, and negative at energies below 1 keV in another. It is hard to explain these results in the context of any spectral variability scenario, but they may signify the presence of a variable, warm absorber in the source. △ Less

Submitted 4 November, 2015; v1 submitted 11 August, 2015; originally announced August 2015.

Comments: Accepted for publication in A&A, 10 pages, 7 figures

Journal ref: A&A 582, A40 (2015)

Spectral energy distributions of the BL Lac PKS 2155-304 from XMM-Newton

Authors: Jai Bhagwan, Alok C. Gupta, I. E. Papadakis, Paul J. Wiita

Abstract: We have used all 20 archival XMM-Newton observations of PKS 2155-304 with simultaneous X-ray and UV/optical data to study its long term flux and spectral variability. We find significant variations, in all bands, on time-scales of years with an rms amplitude of ~ 35-45 per cent, through the optical/UV variations are not correlated with those in the X-ray. We constructed SEDs that span more than th… ▽ More We have used all 20 archival XMM-Newton observations of PKS 2155-304 with simultaneous X-ray and UV/optical data to study its long term flux and spectral variability. We find significant variations, in all bands, on time-scales of years with an rms amplitude of ~ 35-45 per cent, through the optical/UV variations are not correlated with those in the X-ray. We constructed SEDs that span more than three orders of magnitude in frequency and we first fitted them with a log-parabolic model; such models have been applied many times in the past for this, and other, blazars. These fits were poor, so we then examined combined power-law and log-parabolic fits that are improvements. These models indicate that the optical/UV and X-ray flux variations are mainly driven by model normalization variations, but the X-ray band flux is also affected by spectral variations, as parametrized with the model 'curvature' parameter, b. Overall, the energy at which the emitted power is maximum correlates positively with the total flux. As the spectrum shifts to higher frequencies, the spectral 'curvature' increase, in contrast to what is expected if a single log-parabolic model were an acceptable representation of the broad band SEDs. Our results suggest that the optical/UV and X-ray emissions in this source may arise from different lepton populations. △ Less

Submitted 10 September, 2014; originally announced September 2014.

Comments: MNRAS in press (11 pages, 6 figures, 3 tables)

Journal ref: 2014, MNRAS, 444, 3647

Active galactic nuclei synapses: X-ray versus optical classifications using artificial neural networks

Authors: O. Gonzalez-Martin, D. Diaz-Gonzalez, J. A. Acosta-Pulido, J. Masegosa, I. E. Papadakis, J. M. Rodriguez-Espinosa, I. Marquez, L. Hernandez-Garcia

Abstract: (Abridged) Many classes of active galactic nuclei (AGN) have been defined entirely throughout optical wavelengths while the X-ray spectra have been very useful to investigate their inner regions. However, optical and X-ray results show many discrepancies that have not been fully understood yet. The aim of this paper is to study the "synapses" between the X-ray and optical classifications. For th… ▽ More (Abridged) Many classes of active galactic nuclei (AGN) have been defined entirely throughout optical wavelengths while the X-ray spectra have been very useful to investigate their inner regions. However, optical and X-ray results show many discrepancies that have not been fully understood yet. The aim of this paper is to study the "synapses" between the X-ray and optical classifications. For the first time, the new EFLUXER task allowed us to analyse broad band X-ray spectra of emission line nuclei (ELN) without any prior spectral fitting using artificial neural networks (ANNs). Our sample comprises 162 XMM-Newton/pn spectra of 90 local ELN in the Palomar sample. It includes starbursts (SB), transition objects (T2), LINERs (L1.8 and L2), and Seyferts (S1, S1.8, and S2). The ANNs are 90% efficient at classifying the trained classes S1, S1.8, and SB. The S1 and S1.8 classes show a wide range of S1- and S1.8-like components. We suggest that this is related to a large degree of obscuration at X-rays. The S1, S1.8, S2, L1.8, L2/T2/SB-AGN (SB with indications of AGN), and SB classes have similar average X-ray spectra within each class, but these average spectra can be distinguished from class to class. The S2 (L1.8) class is linked to the S1.8 (S1) class with larger SB-like component than the S1.8 (S1) class. The L2, T2, and SB-AGN classes conform a class in the X-rays similar to the S2 class albeit with larger fractions of SB-like component. This SB-like component is the contribution of the star-formation in the host galaxy, which is large when the AGN is weak. An AGN-like component seems to be present in the vast majority of the ELN, attending to the non-negligible fraction of S1-like or S1.8-like component. This trained ANN could be used to infer optical properties from X-ray spectra in surveys like eRosita. △ Less

Submitted 3 April, 2014; originally announced April 2014.

Comments: 15 pages, 7 figures, accepted for publication in A&A. Appendix B only in the full version of the paper here: https://dl.dropboxusercontent.com/u/3484086/AGNSynapsis_OGM_online.pdf

Journal ref: A&A 567, A92 (2014)

General relativistic modelling of the negative reverberation X-ray time delays in AGN

Authors: D. Emmanoulopoulos, I. E. Papadakis, M. Dovciak, I. M. McHardy

Abstract: We present the first systematic physical modelling of the time-lag spectra between the soft (0.3-1 keV) and the hard (1.5-4 keV) X-ray energy bands, as a function of Fourier frequency, in a sample of 12 active galactic nuclei which have been observed by XMM-Newton. We concentrate particularly on the negative X-ray time-lags (typically seen above $10^{-4}$ Hz) i.e. soft band variations lag the hard… ▽ More We present the first systematic physical modelling of the time-lag spectra between the soft (0.3-1 keV) and the hard (1.5-4 keV) X-ray energy bands, as a function of Fourier frequency, in a sample of 12 active galactic nuclei which have been observed by XMM-Newton. We concentrate particularly on the negative X-ray time-lags (typically seen above $10^{-4}$ Hz) i.e. soft band variations lag the hard band variations, and we assume that they are produced by reprocessing and reflection by the accretion disc within a lamp-post X-ray source geometry. We also assume that the response of the accretion disc, in the soft X-ray bands, is adequately described by the response in the neutral iron line (Fe k$α$) at 6.4 keV for which we use fully general relativistic ray-tracing simulations to determine its time evolution. These response functions, and thus the corresponding time-lag spectra, yield much more realistic results than the commonly-used, but erroneous, top-hat models. Additionally we parametrize the positive part of the time-lag spectra (typically seen below $10^{-4}$ Hz) by a power-law. We find that the best-fitting BH masses, M, agree quite well with those derived by other methods, thus providing us with a new tool for BH mass determination. We find no evidence for any correlation between M and the BH spin parameter, $α$, the viewing angle, $θ$, or the height of the X-ray source above the disc, $h$. Also on average, the X-ray source lies only around 3.7 gravitational radii above the accretion disc and the viewing angles are distributed uniformly between 20 and 60 degrees. Finally, there is a tentative indication that the distribution of spin parameters may be bimodal above and below 0.62. △ Less

Submitted 4 February, 2014; originally announced February 2014.

Comments: Accepted for publication in MNRAS. The paper is 22 pages long and contains 19 figures and 2 tables

MSC Class: 85-02 ACM Class: J.2

Evidence for a change in the radiation mechanism in the hard state of GRO J1655-40. Hysteresis in the broad-band noise components

Authors: P. Reig, I. E. Papadakis, M. A. Sobolewska, J. Malzac

Abstract: We have analysed RXTE archival data to study the aperiodic variability of the black-hole binary GRO J1655-40 during the hard state of the 2005 outburst. We computed the 0.008-64 Hz power spectral density during the rise and decay of the 2005 outburst, and we found that they were reasonably well modelled by the sum of two, broad Lorenztian functions in most cases (plus a narrow QPO) which correspon… ▽ More We have analysed RXTE archival data to study the aperiodic variability of the black-hole binary GRO J1655-40 during the hard state of the 2005 outburst. We computed the 0.008-64 Hz power spectral density during the rise and decay of the 2005 outburst, and we found that they were reasonably well modelled by the sum of two, broad Lorenztian functions in most cases (plus a narrow QPO) which correspond to three different variability components. Our aim is to study the evolution of the timing properties of the source during the outburst, by studying the correlation between the characteristics of the broad-band noise components in the power spectra and the source luminosity. Our results suggest that the whole power spectrum shifts to high (low) frequencies as the source luminosity increases (decreases), in agreement with previous studies of other black hole binaries. However, we also detect a strong "hysteresis" pattern in the "frequency-luminosity" plots, and show that the "critical" luminosity limit, above which the timing properties of the source change, is different during the rise and the decay phase of the outburst. We discuss the general implications of these results in the context of the truncated disc model. △ Less

Submitted 2 September, 2013; originally announced September 2013.

Comments: 11 pages, 7 figures, 3 tables

Long-term monitoring of PKS0558-504 with Swift: the disc-corona link

Authors: M. Gliozzi, I. E. Papadakis, D. Grupe, W. P. Brinkmann, C. Raeth

Abstract: PKS 0558-504 is a highly variable, X-ray bright, radio-loud, Narrow-Line Seyfert 1 galaxy with super-Eddington accretion rate and extended jets that do not dominate the emission beyond the radio band. Therefore, this source represents an ideal laboratory to shed some light on the central engine in highly accreting systems and specifically on the link between accretion disc and corona. Here we pres… ▽ More PKS 0558-504 is a highly variable, X-ray bright, radio-loud, Narrow-Line Seyfert 1 galaxy with super-Eddington accretion rate and extended jets that do not dominate the emission beyond the radio band. Therefore, this source represents an ideal laboratory to shed some light on the central engine in highly accreting systems and specifically on the link between accretion disc and corona. Here we present the results from a 1.5 year monitoring with Swift XRT and UVOT. The simultaneous coverage at several wavelengths confirms that PKS 0558-504 is highly variable in any band from optical, to UV and X-rays, with the latter showing the largest amplitude changes but with the UV emission dominating the radiative output. A cross-correlation analysis reveals a tight link between the emission in the optical and UV bands and provides suggestive evidence in favor of a scenario where the variability originates in the outer part of the accretion flow and propagates inwards before triggering the activity of the X-ray emitting corona. Finally, a positive correlation between the soft X-ray flux and the hard photon index suggests that in PKS 0558-504 the seed photons are provided to the corona by the soft excess component. △ Less

Submitted 15 May, 2013; originally announced May 2013.

Comments: 10 pages, 9 figures, 3 tables, accepted for publication in MNRAS

Generating artificial light curves: Revisited and updated

Authors: D. Emmanoulopoulos, I. M. McHardy, I. E. Papadakis

Abstract: The production of artificial light curves with known statistical and variability properties is of great importance in astrophysics. Consolidating the confidence levels during cross-correlation studies, understanding the artefacts induced by sampling irregularities, establishing detection limits for future observatories are just some of the applications of simulated data sets. Currently, the widely… ▽ More The production of artificial light curves with known statistical and variability properties is of great importance in astrophysics. Consolidating the confidence levels during cross-correlation studies, understanding the artefacts induced by sampling irregularities, establishing detection limits for future observatories are just some of the applications of simulated data sets. Currently, the widely used methodology of amplitude and phase randomisation is able to produce artificial light curves which have a given underlying power spectral density (PSD) but which are strictly Gaussian distributed. This restriction is a significant limitation, since the majority of the light curves e.g. active galactic nuclei, X-ray binaries, gamma-ray bursts show strong deviations from Gaussianity exhibiting `burst-like' events in their light curves yielding long-tailed probability distribution functions (PDFs). In this study we propose a simple method which is able to precisely reproduce light curves which match both the PSD and the PDF of either an observed light curve or a theoretical model. The PDF can be representative of either the parent distribution or the actual distribution of the observed data, depending on the study to be conducted for a given source. The final artificial light curves contain all of the statistical and variability properties of the observed source or theoretical model i.e. same PDF and PSD, respectively. Within the framework of Reproducible Research, the code, together with the illustrative example used in this manuscript, are both made publicly available in the form of an interactive Mathematica notebook. △ Less

Submitted 1 May, 2013; originally announced May 2013.

Comments: Accepted for publication in MNRAS. The paper is 23 pages long and contains 21 figures and 2 tables. The Mathematica notebook can be found in the web as part of this paper (Online Material) or at http://www.astro.soton.ac.uk/~de1e08/ArtificialLightCurves/

MSC Class: 85-04; 85-08; 62-04 ACM Class: J.2; G.3; G.4; I.6.8

X-ray spectral analysis of the low-luminosity active galactic nucleus NGC 7213 using long XMM-Newton observations

Authors: D. Emmanoulopoulos, I. E. Papadakis, F. Nicastro, I. M. McHardy

Abstract: We present the X-ray spectral results from the longest X-ray multi-mirror mission-Newton observation, 133 ks, of the low luminosity active galactic nucleus NGC 7213. The hardness ratio analysis of the X-ray light curves discloses a rather constant X-ray spectral shape, at least for the observed exposure time, enabling us to perform X-ray spectral studies using the total observed spectrum. Apart fr… ▽ More We present the X-ray spectral results from the longest X-ray multi-mirror mission-Newton observation, 133 ks, of the low luminosity active galactic nucleus NGC 7213. The hardness ratio analysis of the X-ray light curves discloses a rather constant X-ray spectral shape, at least for the observed exposure time, enabling us to perform X-ray spectral studies using the total observed spectrum. Apart from a neutral Fe Kαemission line, we also detect narrow emission lines from the ionised iron species, Fe xxv and Fe xxvi. Our analysis suggests that the neutral Fe Kαoriginates from a Compton-thin reflector, while the gas responsible for the high ionisation lines is collisionally excited. The overall spectrum, in the 0.3-10 keV energy band, registered by the European Photon Imaging Camera, can be modelled by a power-law component (with a slope of Γ\simeq1.9) plus two thermal components at 0.36 and 8.84 keV. The low-energy thermal component is entirely consistent with the X-ray spectral data obtained by the Reflection Grating Spectrometer between 0.35-1.8 keV. △ Less

Submitted 12 December, 2012; originally announced December 2012.

Comments: Accepted for publication in MNRAS. The paper is 12 pages long and contains 10 figures and 2 tables. Analysis of different source from arXiv:1104.0044

MSC Class: 85-02 ACM Class: J.2

The `harder when brighter' X-ray behaviour of the low luminosity active galactic nucleus NGC 7213

Authors: D. Emmanoulopoulos, I. E. Papadakis, I. M. McHardy, P. Arevalo, D. E. Calvelo, P. Uttley

Abstract: We present the first robust evidence of an anti-correlation between the X-ray photon index, Γ, and the X-ray luminosity in a single low luminosity active galactic nuclei (LLAGN), NGC 7213. Up to now, such anti-correlation trends have been seen only in large samples of LLAGN that span a wide range of X-ray fluxes, although the opposite behaviour (i.e. a positive correlation between Γand X-ray lumin… ▽ More We present the first robust evidence of an anti-correlation between the X-ray photon index, Γ, and the X-ray luminosity in a single low luminosity active galactic nuclei (LLAGN), NGC 7213. Up to now, such anti-correlation trends have been seen only in large samples of LLAGN that span a wide range of X-ray fluxes, although the opposite behaviour (i.e. a positive correlation between Γand X-ray luminosity) has been extensively studied for individual X-ray bright active galactic nuclei. For NGC 7213, we use the long-term X-ray monitoring data of Rossi X-ray Timing Explorer (RXTE), regularly obtained on average every two days from March 2006 to December 2009. Based on our X-ray data, we derive the Γversus flux and the hardness ratio versus flux relations, indicating clearly that NGC 7213 follows a `harder when brighter' spectral behaviour. Additionally, by analysing radio and optical data, and combining data from the literature, we form the most complete spectral energy distribution (SED) of the source across the electromagnetic spectrum yielding a bolometric luminosity of 1.7*10^43 erg s^-1. Phenomenologically, the SED of NGC 7213 is similar to that of low-ionization nuclear emission-line region. The robust anti-correlation trend that we find between Γand X-ray luminosity together with the low accretion rate of the source, 0.14 per cent that of Eddington limit, make NGC 7213 the first LLAGN exhibiting a similar spectral behaviour with that of black hole X-ray binaries in `hard state'. △ Less

Submitted 15 May, 2012; originally announced May 2012.

Comments: Accepted for publication in MNRAS. The paper is 9 pages long and contains 5 figures and 1 table

Revisiting algorithms for generating surrogate time series

Authors: C. Raeth, M. Gliozzi, I. E. Papadakis, W. Brinkmann

Abstract: The method of surrogates is one of the key concepts of nonlinear data analysis. Here, we demonstrate that commonly used algorithms for generating surrogates often fail to generate truly linear time series. Rather, they create surrogate realizations with Fourier phase correlations leading to non-detections of nonlinearities. We argue that reliable surrogates can only be generated, if one tests sepa… ▽ More The method of surrogates is one of the key concepts of nonlinear data analysis. Here, we demonstrate that commonly used algorithms for generating surrogates often fail to generate truly linear time series. Rather, they create surrogate realizations with Fourier phase correlations leading to non-detections of nonlinearities. We argue that reliable surrogates can only be generated, if one tests separately for static and dynamic nonlinearities. △ Less

Submitted 17 August, 2012; v1 submitted 6 November, 2011; originally announced November 2011.

Comments: 5 pages, 4 figures, accepted for publication in PRL

Hard X-ray spectral variability of the brightest Seyfert AGN in the Swift/BAT sample

Authors: M. D. Caballero-Garcia, I. E. Papadakis, F. Nicastro, M. Ajello

Abstract: Aims: We used data from the 58 month long, continuous Swift/Burst Alert Telescope (BAT) observations of the five brightest Seyfert galaxies at hard X-rays, to study their flux and spectral variability in the 20-100 keV energy band. The column density in these objects is less than 10^24 cm-2, which implies that the Swift/BAT data allow us to study the "true" variability of the central source. Resul… ▽ More Aims: We used data from the 58 month long, continuous Swift/Burst Alert Telescope (BAT) observations of the five brightest Seyfert galaxies at hard X-rays, to study their flux and spectral variability in the 20-100 keV energy band. The column density in these objects is less than 10^24 cm-2, which implies that the Swift/BAT data allow us to study the "true" variability of the central source. Results: All objects show significant variations, with an amplitude which is similar to the AGN variability amplitude at energies below 10 keV. We found evidence for an anti-correlation between variability amplitude and black hole mass. The light curves in both bands are well correlated, with no significant delays on time scales as short as 2 days. NGC 4151 and NGC 2110 do not show spectral variability, but we found a significant anti-correlation between hardness ratios and source flux in NGC 4388 (and NGC 4945, IC 4329, to a lesser extent). This "softer when brighter" behaviour is similar to what has been observed at energies below 10 keV, and cannot be explained if the continuum varies only in flux; the intrinsic shape should also steepen with increasing flux. Conclusions: The presence of significant flux variations indicate that the central source in these objects is intrinsically variable on time scales as short as 1-2 days. The intrinsic slope of the continuum varies with the flux (at least in NGC 4388). The positive "spectral slope-flux" correlation can be explained if the temperature of the hot corona decreases with increasing flux. The lack of spectral variations in two objects, could be due to the fact that they may operate in a different "state", as their accretion rate is less than 1% of the Eddington limit (significantly smaller than the rate of the other three objects in the sample). △ Less

Submitted 29 October, 2011; v1 submitted 30 August, 2011; originally announced August 2011.

Comments: Accepted (29/10/11) for publication in A&A (12 pages, containing 14 figures and 2 tables). (Abstract shortened --see link for the complete one)

Negative X-ray reverberation time delays from MCG-6-30-15 and Mrk 766

Authors: D. Emmanoulopoulos, I. M. McHardy, I. E. Papadakis

Abstract: We present an X-ray time lag analysis, as a function of Fourier frequency, for MCG-6-30-15 and Mrk 766 using long term XMM-Newton light curves in the 0.5-1.5 keV and the 2-4 keV energy bands, together with some physical modelling of the corresponding time lag spectra. Both the time lag spectra of MCG-6-30-15 and Mrk 766 show negative values (i.e. soft band variations lag behind the corresponding h… ▽ More We present an X-ray time lag analysis, as a function of Fourier frequency, for MCG-6-30-15 and Mrk 766 using long term XMM-Newton light curves in the 0.5-1.5 keV and the 2-4 keV energy bands, together with some physical modelling of the corresponding time lag spectra. Both the time lag spectra of MCG-6-30-15 and Mrk 766 show negative values (i.e. soft band variations lag behind the corresponding hard band variations) at high frequencies, around 0.001 Hz, similar to those previously observed from 1H 0707-495. The remarkable morphological resemblance between the time lag spectra of MCG-6-30-15 and Mrk 766 indicate that the physical processes responsible for the observed soft time delays is very similar in the two sources, favouring a reflection scenario from material situated very nearby to the central black hole. △ Less

Submitted 29 June, 2011; originally announced June 2011.

Comments: Accepted for publication in MNRAS Letters. The paper (6 pages) contains 2 figures and 1 table

Evidence for a change in the X-ray radiation mechanism in the hard state of Galactic black holes

Authors: M. A. Sobolewska, I. E. Papadakis, C. Done, J. Malzac

Abstract: We present results on spectral variability of two Galactic black hole X-ray binaries, GRO J1655-40 and GX 339-4, in the hard state. We confirm a transition in behaviour of the photon index with luminosity, such that the well known decrease in X-ray photon index with decreasing luminosity only continues down to L_bol ~ 0.01 L_Edd . Below this point the photon index increases again. For Comptonisati… ▽ More We present results on spectral variability of two Galactic black hole X-ray binaries, GRO J1655-40 and GX 339-4, in the hard state. We confirm a transition in behaviour of the photon index with luminosity, such that the well known decrease in X-ray photon index with decreasing luminosity only continues down to L_bol ~ 0.01 L_Edd . Below this point the photon index increases again. For Comptonisation models, this implies that the ratio of the Compton luminosity to seed photon luminosity, lh/ls, changes with bolometric luminosity, consistent with a scenario where seed photons change from cyclo-synchrotron at the lowest luminosities to those from a truncated disc. Alternatively, the transition could mark the point below which the non-thermal jet starts to dominate, or where reprocessed photons replace the viscous ones in an outflowing corona model. △ Less

Submitted 8 June, 2011; originally announced June 2011.

Comments: 10 pages, 5 figures, MNRAS, in press

An XMM-Newton view of the `bare' nucleus of Fairall 9

Authors: D. Emmanoulopoulos, I. E. Papadakis, I. M. McHardy, F. Nicastro, S. Bianchi, P. Arevalo

Abstract: We present the spectral results from a 130 ks observation, obtained from the X-ray Multi-Mirror Mission-Newton (XMM-Newton) observatory, of the type I Seyfert galaxy Fairall 9. An X-ray hardness-ratio analysis of the light-curves, reveals a `softer-when-brighter' behaviour which is typical for radio-quiet type I Seyfert galaxies. Moreover, we analyse the high spectral-resolution data of the reflec… ▽ More We present the spectral results from a 130 ks observation, obtained from the X-ray Multi-Mirror Mission-Newton (XMM-Newton) observatory, of the type I Seyfert galaxy Fairall 9. An X-ray hardness-ratio analysis of the light-curves, reveals a `softer-when-brighter' behaviour which is typical for radio-quiet type I Seyfert galaxies. Moreover, we analyse the high spectral-resolution data of the reflection grating spectrometer and we did not find any significant evidence supporting the presence of warm-absorber in the low X-ray energy part of the source's spectrum. This means that the central nucleus of Fairall 9 is `clean' and thus its X-ray spectral properties probe directly the physical conditions of the central engine. The overall X-ray spectrum in the 0.5-10 keV energy-range, derived from the EPIC data, can be modelled by a relativistically blurred disc-reflection model. This spectral model yields for Fairall 9 an intermediate black-hole best-fit spin parameter of $α=0.39^{+0.48}_{-0.30}$. △ Less

Submitted 31 March, 2011; originally announced April 2011.

Comments: Accepted for publication in MNRAS. The paper contains 11 figures and 1 table

A panchromatic view of PKS 0558-504: an ideal laboratory to study the disk-jet link

Authors: M. Gliozzi, I. E. Papadakis, D. Grupe, W. P. Brinkmann, C. Raeth, L. Kedziora-Chudczer

Abstract: PKS 0558-504 is the brightest radio-loud Narrow-Line Seyfert 1 galaxy at X-ray energies. Here we present results from the radio, optical, UV, and X-ray bands obtained with Swift, XMM, and ATCA during a 10-day monitoring campaign in September 2008. The simultaneous coverage at several wavelengths makes it possible to investigate in detail the broadband spectral energy distribution (SED) and the ene… ▽ More PKS 0558-504 is the brightest radio-loud Narrow-Line Seyfert 1 galaxy at X-ray energies. Here we present results from the radio, optical, UV, and X-ray bands obtained with Swift, XMM, and ATCA during a 10-day monitoring campaign in September 2008. The simultaneous coverage at several wavelengths makes it possible to investigate in detail the broadband spectral energy distribution (SED) and the energetic of this source. The main results can be summarized as follows. The ATCA reveals the presence of an extended radio emission in PKS 0558-504 with two lobe-like structures at ~7" from the bright central source. The extended radio structure and the low value of the radio-loudness similar to radio-quiet Seyfert galaxies coupled with constraints from higher energy bands argue against a jet-dominated emission. The study of the SED, which is dominated by a nearly constant optical-UV emission, supports the conclusion that PKS 0558-504 is accreting at super-Eddington rate. This conclusion was reached assuming M_BH=2.5e8 M_sun, which was obtained with a new scaling method based on X-ray spectral variability results. A comparison between the accretion luminosity and the kinetic power associated with the jet suggests that in this source the accretion power dominates in agreement with the results obtained from Radiation-MHD simulations of Galactic black holes (GBHs) accreting at the Eddington rate. The combined findings from this panchromatic investigation strongly suggest that PKS 0558-504 is a large-scale analog of GBHs in their highly accreting intermediate state. Importantly, PKS 0558-504 may also be the prototype of the parent population of the very radio-loud NLS1s recently detected at gamma-ray energies. △ Less

Submitted 26 May, 2010; originally announced May 2010.

Comments: 11 pages, 5 figures, 3 tables, accepted for publication in ApJ