Four Years of Extreme Ultraviolet Observations of Markarian 421: II. Temporal Analysis

The Extreme Ultraviolet Explorer (EUVE) satellite accumulated ~one million seconds of public data between 1994 and 1997 for the BL Lacertae object Markarian 421. This is the second of two papers in which we present the results of spectral and temporal analysis of this EUVE data set. We analyze in the present paper the imaging data by means of power spectrum and structure function techniques, while the spectral analysis is presented in a companion paper. We find for MRK 421 a power spectrum with slope -2.14 +- 0.28 with a break at ~3 days. This is the first time that a break in the power spectrum of a BL Lacertae object has been found. We also find evidence of non-stationarity for MRK 421 EUV emissionComment: Accepted for publication in the Astrophysical Journal. 16 pages, 14 Postscript figures, 3 Table

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

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.Comment: 18 pages, 12 figures, 6 tables. Accepted for publication in MNRA

Generating artificial light curves: Revisited and updated

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.Comment: 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

Measuring the broadband power spectra of active galactic nuclei with RXTE

We have developed a Monte Carlo technique to test models for the true power spectra of intermittently sampled lightcurves against the noisy, observed power spectra, and produce a reliable estimate of the goodness of fit of the given model. We apply this technique to constrain the broadband power spectra of a sample of four Seyfert galaxies monitored by the Rossi X-ray Timing Explorer (RXTE) over three years. We show that the power spectra of three of the AGN in our sample (MCG-6-30-15, NGC5506 and NGC3516) flatten significantly towards low frequencies, while the power spectrum of NGC5548 shows no evidence of flattening. We fit two models for the flattening, a `knee' model, analogous to the low-frequency break seen in the power spectra of BHXRBs in the low state (where the power-spectral slope flattens to \alpha=0) and a `high-frequency break' model (where the power-spectral slope flattens to \alpha=1), analogous to the high-frequency break seen in the high and low-state power spectra of the classic BHXRB Cyg X-1. Both models provide good fits to the power spectra of all four AGN. For both models, the characteristic frequency for flattening is significantly higher in MCG-6-30-15 than in NGC 3516 (by factor ~10) although both sources have similar X-ray luminosities, suggesting that MCG-6-30-15 has a lower black hole mass and is accreting at a higher rate than NGC 3516. Assuming linear scaling of characteristic frequencies with black hole mass, the high accretion rate implied for MCG-6-30-15 favours the high-frequency break model for this source and further suggests that MCG-6-30-15 and possibly NGC 5506, may be analogues of Cyg X-1 in the high state [ABRIDGED].Comment: 23 pages, accepted for publication in MNRA

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

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.Comment: Accepted for publication by A&

X-ray and optical counterparts of hard X-ray selected sources from the SHEEP survey: first results

We present followup observations of five hard X-ray sources from the ASCA 5-10 keV SHEEP survey, which has a limiting flux of $\sim 10^{-13}$ erg cm$^{-2}$ s$^{-1}$. Chandra data have been obtained to improve the X-ray positions from a few arcmin to $<1''$, which allows unambiguous optical identification. While the objects almost certainly house AGN based on their X-ray luminosity, optical spectroscopy reveals a variety of properties. The identifications indicate that the SHEEP survey samples the same populations as deeper surveys which probe the origin of the X-ray background, but because the SHEEP sources are far brighter, they are more amenable to detailed followup work. We find a variety of classifications and properties, including a type II QSO, a galaxy undergoing star formation, and a broad-line AGN which has a very hard X-ray spectrum, indicating substantial absorption in the X-ray but none in the optical. Two objects have X-ray/optical flux ratios which, were they at an X-ray flux level typical of objects in Chandra deep surveys, would place them in the ``optically faint'' category. They are both identified with broad line QSOs at z$\sim 1$. Clearly this survey - which is relatively unbiased against obscured objects - is revealing a set of remarkable objects quite different to the familiar classes of AGN found in previous optical and soft X-ray surveys.Comment: 5 pages, 3 figures. MNRAS, in pres

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

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 $\Gamma\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.Comment: Accepted for publication in A&A, 10 pages, 7 figure