The absence of a thin disc in M81
A. J. Young I. McHardy D. Emmanoulopoulos S. Connolly


firstpage–lastpage 2016

Two-field Warm Inflation and Its Scalar Perturbations on Large Scales
    Jian-Yang Zhu
    December 30, 2023
=====================================================================






We present the results of simultaneous Suzaku and NuSTAR observations of the nearest Low-Luminosity Active Galactic Nucleus (LLAGN), M81^*. The spectrum is well described by a cut-off power law plus narrow emission lines from Fe Kα, Fexxv and Fexxvi. There is no evidence of Compton reflection from an optically thick disc, and we obtain the strongest constraint on the reflection fraction in M81^* to date, with a best-fit value of R = 0.0 with an upper limit of R < 0.1. The Fe Kα line may be produced in optically thin, N_H = 1 × 10^23 cm^-2, gas located in the equatorial plane that could be the broad line region. The ionized iron lines may originate in the hot, inner accretion flow. The X-ray continuum shows significant variability on ∼ 40 ks timescales suggesting that the primary X-ray source is ∼ 100s of gravitational radii in size. If this X-ray source illuminates any putative optically thick disc, the weakness of reflection implies that such a disc lies outside a few× 10^3 gravitational radii. An optically thin accretion flow inside a truncated optically thick disc appears to be a common feature of LLAGN that are accreting at only a tiny fraction of the Eddington limit.



accretion, accretion discs – galaxies: active – galaxies: individual: M81 – X-rays: galaxies




§ INTRODUCTION


M81, also known as NGC 3031, is a spiral galaxy hosting the closest Low-Luminosity Active Galactic Nucleus (LLAGN), M81^*, at a distance of 3.63 Mpc <cit.>. Its central black hole has a mass of M = 7 × 10^7 M_⊙ <cit.>. The bolometric luminosity of M81^* is 9.3 × 10^40 erg s^-1 <cit.>, which corresponds to 1.1 × 10^-5 L_Edd, where L_Edd is the Eddington luminosity. Thus, logarithmically speaking, M81^* is midway between ultra-low Eddington fraction accreters such as Sgr A* at the centre of the Galaxy with L_Sgr A*∼ 3 × 10^-10 L_Edd <cit.> and a typical high-efficiency AGN with L_AGN∼ 0.1 - 1 L_Edd. M81^* is particularly interesting to study because it is nearby and bright enough to allow detailed spectroscopy, while being in the sub-Eddington regime in which the accretion flow is thought to be qualitatively different to that of an efficiently radiating, high accretion rate AGN.

Orbiting the black hole in M81^* at large, ∼ 10 pc, scales is a disc inclined at 14^∘ to the line of sight (i.e., the angle between the normal to the disc and the observer is 14^∘) as measured from its narrow optical emission lines <cit.>. On smaller scales, there is a 0.8 pc radius (2.5 × 10^5 r_g, where r_g = GM/c^2) disc inclined at 50^∘ to the line of sight, the presence of which is inferred from its broad Hα line <cit.>. The inclination angle of this inner disc is consistent with that of the short, 0.017 pc long, radio jet from the nucleus which is also inclined at 50^∘ <cit.>. These observations suggest that there is a warp between the pc and 10 pc-scale discs.

The thin disc does not extend all the way into the event horizon of the black hole. While X-ray observations with ASCA <cit.> revealed a possibly broadened Fe K line with an Equivalent Width (EW) of 170 eV, no strong Compton reflection component was detected. Subsequent observations with BeppoSAX detected a broad complex of Fe lines with the same EW, 170 eV, but again no Compton reflection component. The BeppoSAX spectrum extended to over 100 keV allowing the reflection fraction R = Ω / 2π,  where Ω is the solid angle subtended by the optically thick gas from the point of view of the X-ray continuum source, to be constrained to be R < 0.3 <cit.>. More recent observations with XMM-Newton were able to resolve the broad Fe K complex into three separate lines of almost equal strength due to Fe Kα, Fexxv and Fexxvi with EWs of 39, 47 and 37 eV, respectively <cit.>. High spectral resolution observations with the Chandra High Energy Transmission Gratings (HETGS) spectrometer showed that the Fexxvi line is redshifted by ∼ 2500  km s^-1 while the neutral Fe Kα line, with an EW of 47 eV, is consistent with being at rest <cit.>. The Fe Kα line is unresolved, with a Full Width at Half Maximum (FWHM) of < 4670  km s^-1. For a disc inclined at 50^∘ this requires r_in≳ 10^4 r_g. For a face-on disc, however, the constraint is much weaker and formally fitting a diskline model <cit.> only constrains r_in > 55 r_g <cit.>. The Si Kα line appears to be broadened, with a FWHM ∼ 1433  km s^-1 which corresponds to r_in≃ 10^5 r_g in a disc inclined at 50^∘. This radius is comparable to that of the broad Hα emitting gas. In this truncated disc model the inner accretion flow is optically thin, radiatively inefficient and may produce the ionized Fexxv and Fexxvi emission lines, characteristic of emission from plasma at temperatures of 10^7 - 10^8  K close to the truncation radius <cit.>. The origin of the Fe Kα line is less straightforward.

X-ray monitoring of M81^* has shown the continuum power-law photon index to be anti-correlated with X-ray flux, so the source is harder when it is brighter <cit.>. This is consistent with models in which the continuum is produced by the inner accretion flow which Comptonises seed photons from the inner edge of truncated disc, with changes in photon index associated with changes of the truncation radius <cit.>.

In this paper we address the questions of where the Fe Kα line originates, what is the geometry of the inner accretion flow, and what is the transition radius between the outer thin disc and inner thick disc. To answer these questions we make use of simultaneous Suzaku and NuSTAR observations to measure the strength of the Fe Kα line, Compton reflection hump and hence reflection fraction. In this way we hope to disentangle whether the Fe Kα line is produced by scattering in optically thin gas, or is produced by reflection from an optically thick disc.

M81 has a redshift of z = -0.0001. The Galactic column density in this direction is N_H (Gal) = 4.22 × 10^20 cm^-2.



§ OBSERVATIONS


The X-ray data presented here come from a coordinated NuSTAR and Suzaku observing programme in May 2015. The observing times are simultaneous for about two days, while the NuSTAR observations extend for a further two days, as illustrated in Fig. <ref> and summarised in Table <ref>. In addition, we have contemporaneous Swift X-ray and AMI radio observations that will be presented elsewhere.





The X-ray data were extracted and analysed using HEAsoft version 6.18, using standard processing techniques as described below.



 §.§ NuSTAR


The NuSTAR data, from ObsID 60101049002, were processed using CALDB version 20150702. The standard NuSTAR pipeline processing tools, nupipeline and nuproducts, were used to produce spectra and light curves. The on-source time was 343 ks which resulted in a total good exposure time of 209 ks, where the difference between these two times is mostly due to the lost exposure time caused by Earth occultations and passages through the South Atlantic Anomaly (SAA) each orbit. Source counts were extracted from a circular aperture of radius 100^'' centred on M81^*, with background counts taken from a source-free circular aperture of the same radius to the north of the nucleus. The spectra from modules A and B were fit jointly, rather than being combined.

The NuSTAR full-band, 3–79 keV, background subtracted light curve in 1000 sec bins was produced using the pipleline, and is shown in Fig. <ref>. The error in the source and background counts was given by √(N), where N is the number of counts in a bin, and these errors were propagated through the calculation to give the error bars shown in the figure.





 §.§ Suzaku


The Suzaku data, from ObsID 710017010, were processed using the CALDB with XIS calibration products dated 2015-01-05 and XRT mirror calibration products dated 2011-06-30. The HXD was not used during this observation. The on-source time was 217 ks which resulted in a total good exposure time of 98 ks for XIS 0 and 100 ks for XIS 1 and 3. The good time is significantly shorter than the observation duration due to the low Earth orbit of the satellite resulting in Earth occultations and passage through the SAA. Source counts were extracted from a circular aperture of radius 260^'' centred on M81^*, excluding the region within 100^'' of the Ultra-Luminous X-ray Source (ULX) M81 X-6 that is ∼ 200^'' south of the nucleus. Background counts were taken from a source-free circular aperture of radius 150^'' to the northeast of the nucleus. The spectra from XIS 0, 1 and 3 were not combined, but fit jointly.

Note that during these observations there was a problem with charge leakage in XIS 0 that resulted in telemetry saturation and therefore some data drop-outs that were taken care of with the pipeline processing. This is likely related to the expansion of the dead area of XIS 0 as described in the ISAS memo JX-ISIS-SUZAKU-MEMO-2015-05[<http://www.astro.isas.jaxa.jp/suzaku/doc/suzakumemo/suzaku_memo_2015-05.pdf>.].



 §.§ Spectral fitting


We simultaneously fit the three Suzaku XIS spectra and the two NuSTAR spectra. To account for calibration differences between instruments and detectors the spectral model is multiplied by a (different) constant for each detector. We find that the two NuSTAR telescope modules agree very well. If the normalisation of FPMA is fixed at unity, the normalisation of FPMB is ∼ 1.02, well within the anticipated calibration uncertainty[See, e.g., the analysis caveats at <http://heasarc.gsfc.nasa.gov/docs/nustar/analysis/>.]. With Suzaku the XIS 1 and XIS 3 detectors are in good agreement with one another, and within about 4% of the NuSTAR detector normalisation. The constant used for XIS 1 is ∼ 0.94 and the constant for XIS  3 is ∼ 0.98. The constant for XIS 0, however, is ∼ 0.85, which is lower than expected but likely linked to problems with the XIS 0 charge leakage. However, there is no apparent difference in the shape of the spectrum from XIS 0 when we compare the residuals of fit models to individual detectors. The precise values of these normalisations change slightly from fit to fit. Unless otherwise stated, the data are binned such that each point in the spectrum from an individual satellite has a signal to noise ratio of at least 20; the background is included in this calculation.



§ RESULTS




 §.§ Light curve
 

The full-band, 3–79 keV, background subtracted NuSTAR light curve is shown in Fig. <ref>, and shows that there is significant, ∼ 30%, variability on timescales of ∼ 40  ks. For a black hole mass of 7 × 10^7 M_⊙, 40 ks corresponds to the light crossing time of 116 r_g. This suggests that a significant fraction of the X-ray continuum originates in a very compact region within a few hundred gravitational radii of the black hole.



 §.§ Simple spectral models


If we model the 0.5–30 keV continuum as an absorbed power law we find significant residuals in the 6–7 keV band due to iron. Excluding the 5–8 keV band and refitting we obtain a good description of the continuum with Γ = 1.89^+0.01_-0.01 and N_H = 6.9^+0.2_-0.2× 10^20cm^-2, having χ^2 / d.o.f. = 2545 / 2193 = 1.16 (where d.o.f. means degrees of freedom). The error bars on parameter values give the 90% statistical confidence interval for one interesting parameter (i.e., corresponding to Δχ^2 = 2.71).

We now include the Fe K band. The iron lines are described by narrow Gaussian lines fixed at 6.40 keV, 6.68 keV and 6.96 keV, representing neutral Fe Kα, He-like Fexxv and H-like Fexxvi, respectively. This fit is shown in Fig. <ref>, and the model parameters for all fits are given in Table <ref>. The normalisations of each of the lines are free parameters. The spectrum continues to be well described by a power law of photon index Γ = 1.89^+0.00_-0.01, column density N_H = 6.9^+0.2_-0.3× 10^20cm^-2,  with line strengths and equivalent widths as given in Table <ref>. This is a good fit, with χ^2 / d.o.f = 2933 / 2550  = 1.15. The Suzaku data allow us to accurately measure the strength of the neutral Fe Kα line while cleanly separating the contribution of Fexxv and Fexxvi in this band. This is important because the neutral reflection we wish to constrain is only expected to contribute to the neutral Fe Kα line. The 0.5–10 keV model flux and luminosity are 3.2 × 10^-11erg cm^-2 s^-1 and 5.1 × 10^40erg s^-1, respectively.





While the signal to noise ratio is dominated by the lowest energies, there is no strong evidence of a systematic deviation from the power law model at high energies, other than a slight hint of spectral softening above ∼ 15  keV. This softening can be modelled by a cut-off power law with an e-folding energy of 267^+133_-67 keV, with this one additional free parameter improving the fit by Δχ^2 = 25, although the reduced χ^2 value is essentially unchanged. There is, however, no clear evidence of any upturn in the spectrum at higher energies as one might expect if there were a strong Compton hump. At soft X-ray energies, a small amount of absorption is required in addition to the Galactic column, and this intrinsic absorption has been reported previously <cit.>. The best-fit cut-off power law plus iron line model is shown in Fig. <ref>.





 §.§ Optically thick reflection modelling of the continuum and iron lines


The continuum modelling described above does not suggest that there is any strong Compton reflection hump in the spectrum of M81^*. We can measure the strength of any putative reflection continuum from optically thick material by including a Compton reflection component to the spectral model. We replace the cut-off power-law model with the pexrav model in which the normalisation, photon index, e-fold energy and reflection fraction are free parameters. The disc inclination is fixed at 50^∘, although allowing this to be a free parameter does not change our conclusions. We find that the best fit value for R is consistent with zero, i.e., no reflection, with a 90% upper confidence limit of R < 0.11, where R = 1 corresponds to an isotropic point source above an infinite disc covering 2 π  sr of the sky. The power-law cut-off energy is constrained to be E_cut = 220^+173_-86 keV. Contours of the constraints on these two parameters are shown in Fig. <ref>. Thus the Kα line cannot be produced by reprocessing in optically thick material.





 §.§ Optically thin reprocessing modelling of the continuum and iron lines


An alternative scenario to reflection from Compton thick material for producing the neutral Fe Kα line is that the line is produced by reprocessing in optically thin material which does not produce a strong Compton hump. We model the continuum, scattered and iron line components from a centrally illuminated torus using the mytorus model applied to a power-law continuum with a 200 keV high energy cut-off <cit.>. The geometry is that of a torus with half-opening angle fixed at 60^∘, and we find an acceptable fit for viewing angle i < 60^∘, i.e., the view of the nucleus is unobscured by the torus. Note that acceptable fits might be possible with different opening angles for the torus but the version of mytorus we used has a fixed opening angle and an acceptable description could be found with this value. The iron Kα line is then produced by scattering of the Γ = 1.9 power-law continuum by the optically-thin[We consider N_H < 10^24 cm^-2 to be optically thin.] torus with a column density N_H = 1.1^+0.3_-0.2× 10^23 cm^-2. This provides a good description of the spectrum, with χ^2 / d.o.f. = 2837 / 2529 = 1.12.

The “torus” might correspond to the Broad Line Region (BLR). In M81^*, 1  pc = 3 × 10^5 r_g, so the size of the torus may be comparable to the size of the BLR which is particularly large in M81 with a radius of maybe ∼ 1  pc <cit.>. So this hypothesis is compatible with the optical properties of the galaxy. The BLR might very well have a “torus”-like geometry because both disc and spherical geometries are able to explain the optical lines <cit.>. The column density N_H ≃ 1 × 10^23 cm^-2 is a reasonable column densify for a BLR cloud <cit.>. The Fe Kα line and velocity resolved Si Kα line <cit.> could both originate in the BLR. In this model the highly ionised Fexxv and Fexxvi lines would be associated with the inner hot, inefficient flow within a few thousand r_g. <cit.> reached a similar conclusion for NGC 7213. In some AGN, such as NGC 5548, the BLR can contribute significantly to the X-ray / optical lags <cit.>.





§ DISCUSSION


We have considered two different scenarios to explain the overall Suzaku and NuSTAR observations of M81^* and, in neither case, is there a detectable Compton hump, with a constraint on the reflection fraction R < 0.11. If the Fe Kα line originates in reflection off Compton-thick material, the ∼ 40 eV EW of this line suggests that, for an illuminating continuum of Γ = 1.9 and solar abundances, the reflection fraction should be R ≃ 0.3-0.4 <cit.>. An optically-thick disc reflection origin for the Fe Kα line is clearly ruled out. Furthermore, the observed limit on reflection fraction implies that the EW of the iron line should be ≲ 12  eV, indicating that the majority of the Fe Kα line flux we observe does not originate from disc reflection.

We can approximate the X-ray continuum source size as h ∼ 100 r_g based on the variability timescale (see Section <ref>). The weak reflection then suggests that if the disc is centrally illuminated by a source at height h, the inner radius r_in∼ 20 h ∼ 2 × 10^3 r_g <cit.>. This all suggests that the inner accretion flow within at least a few× 10^3 r_g is optically thin, and we have a very “clean” view of the nucleus. The X-ray continuum source size we have assumed here is much larger than ≲ 4 r_g typical of more Eddington luminous AGN <cit.>.

Instead of disc reflection, the Fe Kα line is well explained by the mytorus model in which the continuum, scattered and iron line components originate in an optically thin “torus”. The precise geometry of this model has not been fully explored, and is unlikely to be unique, but what we see is consistent with a torus of half-opening angle larger than the inclination angle of the inner disc, with a column density of ∼ 10^23 cm ^-2. It isn't clear at what distance this material lies from the nucleus, as that is not directly constrained by the model, but the torus might correspond to the BLR. A possible geometry for M81^* is shown in Fig. <ref>. The BLR in M81 is unlikely to be an outflow, as a radiatively driven outflow cannot be supported by such a weak nucleus, so a disc or spherical inflow is a more plausible geometry <cit.>.



Our results are consistent with the hypothesis that the inner accretion environment of LLAGN is very clean, and the optically thick accretion disc is truncated at a large radius. Other similar galaxies include the Seyfert 2 galaxy NGC 2110 with R < 0.14 <cit.> and the low accretion rate LINER NGC 7213 with R < 0.13 <cit.>. These galaxies are all radiating significantly below the Eddington limit, with L / L_Edd = 0.25 - 3.7 × 10^-2 for NGC 2110 <cit.> and L / L_Edd = 1.4 × 10^-3 for NGC 7213 <cit.>. What is not clear, however, is whether there are low Eddington fraction sources that have high reflection fractions. While low Eddington fraction galaxies can have fairly large EW iron lines <cit.> this does not necessarily mean their reflection fractions are large – this is not the case in NGC 2110, for example. For the AGN in the XXL survey <cit.> with L / L_Edd≳ 10^-2 the reflection fraction is typically R > 0.3 <cit.>, although since these authors have to perform a stacking analysis it isn't clear what R is for the very lowest Eddington fraction galaxies. An illustration of the reflection fraction vs. Eddington fraction is shown in Fig. <ref>, which suggests that the transition between high (R > 0.3) and low (R ≃ 0) reflection fraction may be quite abrupt. Higher signal-to-noise ratio observations of low Eddington fraction AGN are necessary to clarify this picture, to fill in or exclude galaxies from the upper left quadrant of this plot.



The Fexxv and Fexxvi lines may originate in hot, optically thin gas within a few× 10^3 r_g. Highly ionized iron lines are also seen in NGC 7213 <cit.>, a galaxy with very similar spectral properties to M81^*. Observations of “bare” Seyfert galaxies also show Fexxv and Fexxvi lines <cit.>, suggesting that these hot iron lines may be a common component of AGN spectra. The LLAGN NGC 2110, however, does not show evidence of ionized iron lines <cit.>, but NGC 2110 is a Seyfert 2 galaxy so does not have a “bare” nucleus.



§ CONCLUSIONS


In this paper we present simultaneous Suzaku and NuSTAR observations of the LLAGN M81^*. We find the following.


    
  * The NuSTAR light curve of shows significant, ∼ 30%, variability on timescales of ∼ 40 ks.
    
  * The continuum spectrum is well described by a cut-off power law with photon index Γ = 1.9 and cut-off energy of E_c ≃ 220-270 keV.
    
  * There are weak Fe Kα, Fexxiv and Fexxvi lines with EWs of 40, 38 and 33 eV, respectively.
    
  * There is evidence of weak intrinsic absorption, with the absorbing column density N_H = 6.3^+0.3_-0.3× 10^20 cm^-2 exceeding the Galactic value N_H(Gal) = 4.2 × 10^20 cm^-2.
    
  * There is no Compton reflection hump, with the reflection fraction measured to be R < 0.11.
    
  * We model the Fe Kα line as being produced in an optically thin, N_H = 1.1 × 10^23 cm^-2 torus using the mytorus model <cit.>.
    
  * We hypothesise that the primary X-ray emission originates in a hot, optically thin inner accretion flow within ∼ 10^2 r_g. Surrounding this is the “torus”, with a half-opening angle larger than 50^∘, that extends from ∼ 10^2 - 10^5 r_g. At larger radii, ≳ 10^5 r_g, is a warped, optically thick, geometrically thin disc.




§ ACKNOWLEDGEMENTS


This research has made use of the NuSTAR Data Analysis Software (NuSTARDAS) jointly developed by the ASI Science Data Center (ASDC, Italy) and the California Institute of Technology (Caltech, USA). We are grateful to Katja Pottschmidt and Kenji Hamaguchi for helpful discussion regarding the Suzaku calibration and the problem we had with XIS 0.

mnras



