The 2-79 keV X-ray Spectrum of the Circinus Galaxy with NuSTAR, XMM-Newton and Chandra: a Fully Compton-Thick AGN

The Circinus galaxy is one of the nearest obscured AGN, making it an ideal target for detailed study. Combining archival Chandra and XMM-Newton data with new NuSTAR observations, we model the 2-79 keV spectrum to constrain the primary AGN continuum and to derive physical parameters for the obscuring material. Chandras high angular resolution allows a separation of nuclear and off-nuclear galactic emission. In the off-nuclear diffuse emission we find signatures of strong cold reflection, including high equivalent-width neutral Fe lines. This Compton-scattered off-nuclear emission amounts to 18% of the nuclear flux in the Fe line region, but becomes comparable to the nuclear emission above 30 keV. The new analysis no longer supports a prominent transmitted AGN component in the observed band. We find that the nuclear spectrum is consistent with Compton-scattering by an optically-thick torus, where the intrinsic spectrum is a powerlaw of photon index $Gamma = 2.2-2.4$, the torus has an equatorial column density of $N_{rm H} = (6-10)times10^{24}$cm$^{-2}$ and the intrinsic AGN $2-10$ keV luminosity is $(2.3-5.1)times 10^{42}$ erg/s. These values place Circinus along the same relations as unobscured AGN in accretion rate-vs-$Gamma$ and $L_X$-vs-$L_{IR}$ phase space. NuSTARs high sensitivity and low background allow us to study the short time-scale variability of Circinus at X-ray energies above 10 keV for the first time. The lack of detected variability favors a Compton-thick absorber, in line with the the spectral fitting results.

A Mexican Hat with holes: calculating low resolution power spectra from data with gaps

A simple method for calculating a low-resolution power spectrum from data with gaps is described. The method is a modification of the $Delta$-variance method previously described by Stutzki and Ossenkopf. A Mexican Hat filter is used to single out fluctuations at a given spatial scale and the variance of the convolved image is calculated. The gaps in the image, defined by the mask, are corrected for by representing the Mexican Hat filter as a difference between two Gaussian filters with slightly different widths, convolving the image and mask with these filters and dividing the results before calculating the final filtered image. This method cleanly compensates for data gaps even if these have complicated shapes and cover a significant fraction of the data. The method was developed to deal with problematic 2D images, where irregular detector edges and masking of contaminating sources compromise the power spectrum estimates, but it can also be straightforwardly applied to 1D timing analysis or 3D data cubes from numerical simulations.

Optical and Near-IR long-term monitoring of NGC3783 and MR2251-178: evidence for variable near-IR emission from thin accretion discs

We present long term near-IR light curves for two nearby AGN: NGC3783 and MR2251-178. The near-IR data are complemented with optical photometry obtained over the same period of time. The light curves in all bands are highly variable and good correlations can be seen between optical and NIR variations. Cross-correlation analysis for NGC 3783 suggests that some disc near-IR emission is present in the J-band flux, while the H and K-bands are dominated by emission from a torus located at the dust sublimation radius. For MR2251-178 the cross-correlation analysis and the optical-near-IR flux-flux plots suggest that the near-IR flux is dominated by disc emission. We model the optical to near-IR Spectral Energy Distributions (SED) of both sources and find that disc flaring might be a necessary modification to the geometry of a thin disc in order to explain the observations. The SED of MR2251-178 gives some indications for the presence of NIR emission from a torus. Finally, we consider the implications of the standard alpha disc model to explain the different origin of the variable NIR emission in these AGN.

Twelve years of X-ray and optical variability in the Seyfert galaxy NGC 4051

We discuss the origin of the optical variations in the Narrow line Seyfert 1 galaxy NGC 4051 and present the results of a cross-correlation study using X-ray and optical light curves spanning more than 12 years. The emission is highly variable in all wavebands, and the amplitude of the optical variations is found to be smaller than that of the X-rays, even after correcting for the contaminating host galaxy flux falling inside the photometric aperture. The optical power spectrum is best described by an unbroken power law model with slope $alpha=1.4^{+0.6}_{-0.2}$ and displays lower variability power than the 2-10 keV X-rays on all time-scales probed. We find the light curves to be significantly correlated at an optical delay of $1.2^{+1.0}_{-0.3}$ days behind the X-rays. This time-scale is consistent with the light travel time to the optical emitting region of the accretion disc, suggesting that the optical variations are driven by X-ray reprocessing. We show, however, that a model whereby the optical variations arise from reprocessing by a flat accretion disc cannot account for all the optical variability. There is also a second significant peak in the cross-correlation function, at an optical delay of $39^{+2.7}_{-8.4}$ days. The lag is consistent with the dust sublimation radius in this source, suggesting that there is a measurable amount of optical flux coming from the dust torus. We discuss the origin of the additional optical flux in terms of reprocessing of X-rays and reflection of optical light by the dust.

A correlation between the spectral and timing properties of AGN

We present the results from a combined study of the average X-ray spectral and timing properties of 14 nearby AGN. For 11 of the sources in the sample, we used all the available data from the RXTE archive, which were taken until the end of 2006. There are 7795 RXTE observations in total for these AGN, obtained over a period of ~7-11 years. We extracted their 3-20 keV spectra and fitted them with a simple power-law model, modified by the presence of a Gaussian line (at 6.4 keV) and cold absorption, when necessary. We used these best-fit slopes to estimate the mean spectral slope for each object, while we used results from the literature to estimate the average spectral slope of the three objects without archival, monitorin RXTE data. Our results show that the AGN average spectral slopes are not correlated either with the black hole mass or the characteristic frequencies that were detected in the power spectra.They are positively correlated, though, with the characteristic frequency when normalised to the sources black hole mass. This is similar to the spectral-timing correlation that has been observed in Cyg X-1, but not the same.The AGN spectral-timing correlation can be explained if we assume that the accretion rate determines both the average spectral slope and the characteristic time scales in AGN. The spectrum should steepen and the characteristic frequency should increase, proportionally, with increasing accretion rate. We also provide a quantitative expression between spectral slope and accretion rate. Thermal Comptonisation models are broadly consistent with our result, but only if the ratio of the soft photons luminosity to the power injected to the hot corona is proportionally related to the accretion rate.

Long term optical and X-ray variability of the Seyfert galaxy Markarian 79

We present the results of concurrent X-ray and optical monitoring of the Seyfert 1 galaxy Mrk 79 over a period of more than five years. We find that on short to medium time-scales (days to a few tens of days) the 2-10 keV X-ray and optical u and V band fluxes are significantly correlated, with a delay between the bands consistent with zero days. We show that most of these variations may be well reproduced by a model where the short-term optical variations originate from reprocessing of X-rays by an optically thick accretion disc. The optical light curves, however, also display long time-scale variations over thousands of days, which are not present in the X-ray light curve. These optical variations must originate from an independent variability mechanism and we show that they can be produced by variations in the (geometrically) thin disc accretion rate as well as by varying reprocessed fractions through changes in the location of the X-ray corona.

Correlated X-ray/Optical Variability in the Quasar MR2251-178

Emission from Active Galactic Nuclei is known to vary strongly over time over a wide energy band, but the origin of the variability and especially of the inter-band correlations is still not well established. Here we present the results of our X-ray and optical monitoring campaign of the quasar MR2251-178, covering a period of 2.5 years. The X-ray 2-10 keV flux is remarkably well correlated with the optical B, V and R bands, their fluctuations are almost simultaneous with a delay consistent with 0 days and not larger than 4 days in either direction. The amplitude of variations shows an intriguing behaviour: rapid, large amplitude fluctuations over tens of days in the X-rays have only small counterparts in the optical bands, while the long-term trends over hundreds of days are stronger in the B band than in X-rays. We show that simple reprocessing models, where all the optical variability arises from the variable X-ray heating, cannot simultaneously explain the discrepant variability amplitudes on different time-scales and the short delays between X-ray and optical bands. We interpret the variability and correlations, in the optically-thick accretion disc plus corona scenario, as the result of intrinsic accretion rate variations modulating both X-ray and optical emission, together with reprocessing of X-rays by the accretion disc.

Fourier-resolved energy spectra of the Narrow Line Seyfert 1 Mkn 766

We compute Fourier-resolved X-ray spectra of the Seyfert 1 Markarian 766 to study the shape of the variable components contributing to the 0.3-10 keV energy spectrum and their time-scale dependence. The fractional variability spectra peak at 1-3 keV, as in other Seyfert 1 galaxies, consistent with either a constant contribution from a soft excess component below 1 keV and Compton reflection component above 2 keV, or variable warm absorption enhancing the variability in the 1-3 keV range. The rms spectra, which shows the shape of the variable components only, is well described by a single power law with an absorption feature around 0.7 keV, which gives it an apparent soft excess. This spectral shape can be produced by a power law varying in normalisation, affected by an approximately constant (within each orbit) warm absorber, with parameters similar to those found by Turner et al. for the warm-absorber layer covering all spectral components in their scattering scenario. The total soft excess in the average spectrum can therefore be produced by a combination of constant warm absorption on the power law plus an additional less variable component. On shorter time-scales, the rms spectrum hardens and this evolution is well described by a change in power law slope, while the absorption parameters remain the same. The frequency dependence of the rms spectra can be interpreted as variability arising from propagating fluctuations through an extended emitting region, whose emitted spectrum is a power law that hardens towards the centre. This scenario reduces the short time-scale variability of lower energy bands making the variable spectrum harder on shorter time-scales and at the same time explains the hard lags found in these data by Markowitz et al.

Discovery of multiple Lorentzian components in the X-ray timing properties of the Narrow Line Seyfert 1 Ark 564

We present a power spectral analysis of a 100 ksec XMM-Newton observation of the narrow line Seyfert 1 galaxy Ark~564. When combined with earlier RXTE and ASCA observations, these data produce a power spectrum covering seven decades of frequency which is well described by a power law with two very clear breaks. This shape is unlike the power spectra of almost all other AGN observed so far, which have only one detected break, and resemble Galactic binary systems in a soft state. The power spectrum can also be well described by the sum of two Lorentzian-shaped components, the one at higher frequencies having a hard spectrum, similar to those seen in Galactic binary systems. Previously we have demonstrated that the lag of the hard band variations relative to the soft band in Ark 564 is dependent on variability time-scale, as seen in Galactic binary sources. Here we show that the time-scale dependence of the lags can be described well using the same two-Lorentzian model which describes the power spectrum, assuming that each Lorentzian component has a distinct time lag. Thus all X-ray timing evidence points strongly to two discrete, localised, regions as the origin of most of the variability. Similar behaviour is seen in Galactic X-ray binary systems in most states other than the soft state, i.e. in the low-hard and intermediate/very high states. Given the very high accretion rate of Ark 564 the closest analogy is with the very high (intermediate) state rather than the low-hard state. We therefore strengthen the comparison between AGN and Galactic binary sources beyond previous studies by extending it to the previously poorly studied very high accretion rate regime.