Mass of the active galactic nucleus black hole XMMUJ134736.6+173403

TL;DR

This study analyzes X-ray observations of a quasar, revealing quasiperiodic oscillations that suggest the presence of a supermassive black hole with a mass between 10^7 and 10^9 solar masses, and discusses implications for resonance models.

Contribution

It provides the first mass estimates of an active galactic nucleus black hole based on observed QPOs and tests resonance model predictions against observational data.

Findings

Detected two QPOs with a 3:1 ratio in X-ray flux.

Estimated black hole mass between 10^7 and 10^9 solar masses.

Supports the 3:2 epicyclic resonance model for QPOs.

Abstract

A recent study of the X-ray source XMMUJ134736.6+173403 has revealed a strong quasi-periodic modulation of the X-ray flux. The observation of two quasiperiodic oscillations (QPOs) that occur on a daily timescale and exhibit a 3:1 frequency ratio strongly supports the evidence for the presence of an active galactic nucleus black hole (AGN BH). Assuming the orbital origin of QPOs we calculate the upper and lower limit on AGN BH mass M arriving at $M\approx 10^7-10^9M_{\odot}$. Comparing that to mass estimates of other sources, MMUJ134736.6+173403 appears to be the most massive source with commensurable QPO frequencies, and its mass represents the current observational upper limit on AGN BH mass based on QPO observations. We note it will be crucial for the falsification of particular resonance models of QPOs whether only single QPO with frequency that completes the harmonic sequence 3:2:1 is found in this source, or whether there is a new different pair of QPOs with frequencies in the 3:2 ratio. The former case would agree with the prediction of the 3:2 epicyclic resonance model and black hole mass $M\approx(5a^2 + 8a + 8)\times 10^{7}M_{\odot}$, where a is a dimensionless BH spin.

Tables

Table 1. Table 1: Frequency relations corresponding to individual QPO models investigated in Török et al. ( 2011 ) . The relations are expressed in terms of three fundamental frequencies of perturbed circular geodesic motion.

Model	Relations
$\mathrm{𝐑𝐏}$	${\nu }_L={\nu }_K-{\nu }_r$	${\nu }_U={\nu }_K$
${\mathrm{𝐓𝐃}}^{\ast }$	${\nu }_L={\nu }_K$	${\nu }_U={\nu }_K+{\nu }_r$
$\mathrm{𝐖𝐃}$	${\nu }_L=2\left({\nu }_K-{\nu }_r\right)$	${\nu }_U=2{\nu }_K-{\nu }_r$
${\mathrm{𝐄𝐑}}^{\nabla }$	${\nu }_L={\nu }_r$	${\nu }_U={\nu }_{\theta }$
${\mathrm{𝐊𝐞𝐏}}^{\nabla }$	${\nu }_L={\nu }_r$	${\nu }_U={\nu }_K$
$\mathrm{𝐑𝐏𝟏}$	${\nu }_L={\nu }_K-{\nu }_r$	${\nu }_U={\nu }_{\theta }$
$\mathrm{𝐑𝐏𝟐}$	${\nu }_L={\nu }_K-{\nu }_r$	${\nu }_U=2{\nu }_K-{\nu }_{\theta }$

Equations

$${\nu }_U=\frac{2.8\mathrm{kHz}}{M^{\ast }},$$

$${\nu }_{ISCO}=\frac{2.20\mathrm{kHz}}{M^{\ast }},{\nu }_{RISCO}=\frac{4.04\mathrm{kHz}}{M^{\ast }}.$$

$${\nu }_{ISCO}={\nu }_{RISCO}=\frac{16.2\mathrm{kHz}}{M^{\ast }}.$$

$${\nu }_{\mathrm{i}}=\frac{1}{T_{\mathrm{i}}}=\frac{c^3}{2\pi G}\frac{1}{M}{\mathcal{F}}_{\mathrm{i}}(x,a)\le {\nu }_{ISCO}(M,a),$$

$$M\approx (5_{-3}^{+11}{a}^2+8_{-4}^{+17}a+8_{-4}^{+12})\times 10^7{M}_{\text{sun}},$$

Full text

1 Introduction

The X–ray power density spectra (PDS) of several Galactic black hole (BH) systems displays so-called high-frequency quasi-periodic oscillations (HF QPOs) within the range of $40-450$Hz. It has been noted that their frequencies lie within the range corresponding to timescales of orbital motion in the vicinity of a BH. This coincidence is thought to represent a strong indication that the observed signal originates in the innermost parts of an accretion disk. This hypothesis also finds support in the field of Fourier-resolved spectroscopy (e.g., Gilfanov et al., 2000; McClintock & Remillard, 2006; van der Klis, 2006).

Detections of elusive HF QPO peaks in Galactic microquasars are often reported at rather constant frequencies that are characteristic for a given source. It has been found that they usually appear in ratios of small natural numbers (Abramowicz & Kluźniak, 2001; Remillard et al., 2002; McClintock & Remillard, 2006), typically in a $3\!:\!2$ ratio (see, however, Belloni et al., 2012; Belloni & Altamirano, 2013; Varniere & Rodriguez, 2018). The evidence for the rational frequency ratios was also discussed in the context of neutron star (NS) QPOs. In the NS sources, clustering of twin-peak QPO detections most frequently arises as a result of the weakness of (one or both) QPOs outside the limited range of QPO frequencies (frequency ratio; see Abramowicz et al., 2003; Belloni et al., 2005, 2007; Török et al., 2008a, b; Barret & Boutelier, 2008; Boutelier et al., 2010).

The $3\!:\!2$ frequencies observed in the Galactic microquasars are well matched by the relation (McClintock & Remillard, 2006)

[TABLE]

where $\nu_{\scriptscriptstyle\!U}$ is the higher of the two frequencies that form the $3\!:\!2$ ratio, $R=\nu_{\scriptscriptstyle\!U}/\nu_{\scriptscriptstyle\!L}=3/2$, and the BH gravitational mass is given as $M^{*}=M/M_{\odot}$.

Abramowicz et al. (2004) suggested that detections of ${3\!:\!2}$ QPOs and scaling (1) could provide the basis for a method intended to determine the mass of BH sources such as active galactic nuclei (AGN) and ultraluminuos X-ray (ULX) sources. It has been argued that a confirmation of ${3\!:\!2}$ QPOs in other than stellar mass BH sources could be of fundamental importance for the BH accretion theory Abramowicz et al. (2004); Török (2005a, b). It was furthermore discussed that within a wide range of BH masses $M,$ both the rotational parameter $a\equiv\mathrm{c}J/\mathrm{G}M^{2}$ and specific details of a given orbital QPO model are of secondary importance, and that the observed frequencies can be used for to estimate $M$.111In this paper we assume Kerr BH space times.

A full decade after these results have been published, Zhou et al. (2015) completed a survey pointing to a further observational evidence for the large-scale validity of the $1/M$ scaling of the BH QPO frequencies (1). Even more important findings that support this evidence have been accumulated recently. At present, several sources that span the range of a few orders of magnitude seem to follow the relation. This is illustrated in Figure 1, and further details along with several references can be found in Zhou et al. (2015), Pasham et al. (2014, 2015), Czerny et al. (2016), Zhang et al. (2017, 2018), Gupta et al. (2018), and also Smith et al. (2018). We here discuss the implications of the latest QPO observations and mainly focus on the importance of QPOs that have recently been found in the X-ray source XMMUJ134736.6+173403  by Carpano & Jin (2018).

2 Upper bounds on AGN mass based on HF QPOs

Carpano & Jin (2018) carried out a systematic study of Chandra, Swift, and XMM-Newton observations of the X-ray source XMMUJ134736.6+173403,  which has been previously found to be spatially coincident with a pair of galaxies including the Seyfert 2 galaxy SDSS J134736.39+173404.6. Using the Chandra observation from 2008, Carpano & Jin (2018) accurately evaluated the position of this X-ray source and showed that its coordinates coincide with the position of the Seyfert 2 galaxy. Within the analysis of a set of 29 Swift observations conducted between February 6 and May 23, 2008, they discovered strong twin-peak quasi-periodic oscillations with periods of 23.82$\pm$0.07 h and 71.44$\pm$0.57 h.

This measurement clearly indicates that $R=3.00\pm{0.04}$, and it provides very strong evidence for frequency-commensurable QPOs in an AGN BH. The frequency commensurability is crucial for the nonlinear resonance models discussed by Abramowicz, Kluzniak and collaborators (e.g., Abramowicz & Kluźniak, 2001; Rebusco, 2004; Horák, 2008; Török et al., 2005). Motivated by the observed $3:1$ frequency ratio, Carpano & Jin (2018) investigated the mass-spin relations expected for this source in detail from various resonances between disk-oscillation modes that have previously been discussed in the context of QPOs in Galactic microquasars. We make several remarks regarding these models in Sections 3 and 4, where other QPO models are considered as well. In what follows, we explore general implications of the observed rapid variability that are valid for a large class of orbital QPO models.

The Keplerian frequency of matter orbiting a BH monotonically increases as the orbital radius $r$ decreases to the inner edge of the accretion disk. The location of the inner edge depends on the radiative efficiency of the disk. For a very high efficiency, the disk terminates at the marginally stable circular orbit (thin disks), $r_{{\scriptscriptstyle\!I}{\scriptscriptstyle\!S}{\scriptscriptstyle\!C}{\scriptscriptstyle\!O}}$ (often called innermost stable circular orbit, or ISCO), while for a very low efficiency, it terminates at the marginally bound orbit (ion tori, thin disks, and advection dominated accretion flows - ADAFs), $r_{{\scriptscriptstyle\!R}{\scriptscriptstyle\!I}{\scriptscriptstyle\!S}{\scriptscriptstyle\!C}{\scriptscriptstyle\!O}}$. The Keplerian orbital frequency at these orbits for a Schwarzschild BH ($a=0$) scales with BH mass as (e.g., Bardeen et al., 1972)

[TABLE]

For rotating BHs (corrotating disks), these frequencies are higher, and for maximally rotating Kerr BH ($a=1$), we may write

[TABLE]

In Figure 1 we present relations that determine the highest orbital frequencies. We also mark here the higher QPO frequency observed in XMMUJ134736.6+173403. Postulating that this frequency corresponds to the Keplerian frequency inside the disk, we can find that the mass of the source should not be higher than $M\doteq 1.1\times 10^{9}M_{\sun}$.

3 Application of particular orbital models of QPOs and lower bounds on AGN mass

A large group of models, to which we refer as standard geodesic orbital (SGO) models, assume that the observed frequencies, $\nu_{\scriptscriptstyle\!L}$ and $\nu_{\scriptscriptstyle\!U}$, are equal to fundamental frequencies of (geodesic) orbital motion, that is, the Keplerian frequency defined at a circular orbit above ISCO, $r\geq r_{{\scriptscriptstyle\!I}{\scriptscriptstyle\!S}{\scriptscriptstyle\!C}{\scriptscriptstyle\!O}}$, and the radial and vertical epicyclic frequency, or to their linear combinations (including the periastron and Lense-Thirring precession frequency). The three frequencies $\nu_{\scriptscriptstyle\!K}$, $\nu_{\mathrm{r}}$ , and $\nu_{\theta}$ can be written in a general form

[TABLE]

where $x=r/r_{G}\geq r_{{\scriptscriptstyle\!I}{\scriptscriptstyle\!S}{\scriptscriptstyle\!C}{\scriptscriptstyle\!O}}/r_{G}$, $r_{G}=GM/c^{2}$, and $\mathcal{F}_{\mathrm{i}}$ is a specific function that for a fixed BH spin depends only on $x$. The $1/M$ term here provides a physical explanation of the empiric scaling (1). We illustrate the behavior of Keplerian and epicyclic frequencies for the simple case of a Schwarzschild BH (a=0) in Figure 2a. We also include in this figure an illustration of the behavior of the normalized disk flux $f=f/f_{\mathrm{max}}$ calculated for relativistic thin disks (Page & Thorne, 1974).

The family of SGO models includes various physical concepts that hold the assumption that the QPO excitation radii are located within the most luminuos accretion region ($r\!\ll\!r_{90\%}$, see Figure 2a), usually below $20r_{\mathrm{G}}$. Several models assume, for instance, that QPOs are produced by a local motion of accreted inhomogeneities such as blobs or vortices. This subset of SGO models includes the so-called relativistic precession model (RP) or tidal disruption (TD) model (Abramowicz et al., 1992; Stella & Vietri, 1998; Stella et al., 1999; Čadež et al., 2008; Kostić et al., 2009; Bakala et al., 2014; Karssen et al., 2017; Germanà, 2017). Another possibility is to relate the QPOs to a collective motion of the accreted matter, in particular, to some accretion disk oscillatory modes Wagoner et al. (2001); Rezzolla et al. (2003); Abramowicz et al. (2006); Ingram & Done (2010); Fragile et al. (2016). Specific models of this type often deal with oscillations in a slender accretion torus and assume some kind of resonance between these oscillation modes.

A general idea considering such resonances has originated in Abramowicz & Kluźniak (2001) (see also Aliev & Galtsov, 1981) and was later continued in the context of BH QPOs and spin extensively discussed in several subsequent works covering both parametric and forced resonances (Török et al., 2005; Török, 2005a, b; Török et al., 2011; Stuchlík et al., 2013). Various combinations of both axisymmetric and non-axisymmetric epicyclic modes (EP class of models), as well as combinations of these modes and the Keplerian oscillation (KP class of models), have been considered. An often-quoted model is represented by the $3\!:\!2$ epicyclic resonance model, or by a similar model that consideres the Keplerian frequency instead of the vertical epicyclic frequency (ER and KeP models).

The mass-spin relations predicted for Galactic microquasars by several QPO models have been investigated in a series of works (Török et al., 2011; Kotrlová et al., 2014, 2017). Along with the TD, RP, ER, and KeP model, these include the specific diskoseismic model of Kato (2001, 2007, 2008) and two models introduced by Bursa (2005) and Török et al. (2007, 2010) that were motivated by the observations of NS QPOs (the RP1 and RP2 models). We recall in Table1 how the QPO frequencies are identified with orbital frequencies of geodesic motion.

We calculated the mass-spin relations for XMMUJ134736.6+173403  based on the models listed in Table 1. The obtained results are presented in Figure 2b along with all mass-spin relations relevant to the models considered by Carpano & Jin (2018). The mass predicted by these models is always below the ISCO limit $M\doteq 1.1\times 10^{9}M_{\sun}$. It is nevertheless close to this limit (see the horizontal range of Figure 1), and there is $10^{7}M_{\sun}<M<10^{9}M_{\sun}$, which is generic for the proposed SGO models when $a\geq 0$ because they involve QPO excitation radii $r<r_{90\%}$. In Figure 2b we illustrate these findings and plot functions $f_{90\%}\equiv M(r_{90\%})$ and $f_{{\scriptscriptstyle\!I}{\scriptscriptstyle\!S}{\scriptscriptstyle\!C}{\scriptscriptstyle\!O}}\equiv M(r_{{\scriptscriptstyle\!I}{\scriptscriptstyle\!S}{\scriptscriptstyle\!C}{\scriptscriptstyle\!O}})$ that correspond to an identification of the observed QPO frequency with the appropriate Keplerian frequency.

4 Implications for models assuming parametric resonances

Kluźniak & Abramowicz (2002) have suggested that the lower and upper kHz QPOs may be identified with the radial and vertical axisymmetric epicyclic oscillations of an accretion disk. In their scenario, the resonance arises as a result of nonlinear coupling between the oscillation modes. The radial mode supplies energy to the vertical mode by means of the parametric resonance. As the parametric resonance occurs when it is $\nu_{r}/\nu_{\theta}=2/n$ (with $n$ being an integer number) and $\nu_{r}\leq\nu_{\theta}$, their model naturally explains the $3\!:\!2$ frequency ratio that is often observed in the QPO sources. A quantitative analysis of this model has been presented by Rebusco (2004) and Horák (2008). The latter work is devoted to a detailed study of the nature of nonlinear coupling between disk oscillation modes and extends the previous analysis to non-axisymmetric epicyclic modes. Consideration of these modes implies frequency identifications of the lower and upper QPO that involve combinations of both Keplerian and epicyclic frequencies because they oscillate with frequencies $m\nu_{\mathrm{K}}\pm\nu_{\mathrm{r},\theta}$, where $m$ is an integer azimuthal wavenumber.

4.1 Restrictions on resonant interactions

Formulae of the Ep, RP1, and RP2 models from Table 1 combine various combinations of the radial and vertical epicylic modes with azimuthal wavenumbers $(m_{\mathrm{r}},m_{\theta})=(0,0)$, $(1,0),$ and $(1,2)$, respectively. The symmetry properties of the accretion flow strongly limit the resonant interaction between the modes. In particular, when the unperturbed flow is both axially symmetric and symmetric with respect to the equatorial plane, possible resonances up to the fourth order occur only when $\nu_{\theta}/\nu_{r}=1/2,1,1/4,$ and $3/2$, and the same condition has also to be satisfied by the azimuthal wavenumbers, that is, $m_{\theta}/m_{r}=1/2$, 1, $1/4$ or $3/2$. In this context, in the case of the $3\!:\!2$ observed QPO frequency ratio, the idea of parametric resonance is consistent only with the frequency identification provided by the ER model (Horák, 2008). Even this model nevertheless fails in explaining the 3:1 frequency ratio observed in XMMUJ134736.6+173403.

4.2 $3\!:\!2$ resonance

A resolution of the problem may be connected to a nonlinear nature of both the flow oscillations and the modulation mechanism through which the oscillations manifest themselves in the observed light curves. The observed power spectra often contain higher harmonics and subharmonics of the fundamental mode frequencies. The 3:1 frequency ratio then still may be explained by the $3\!:\!2$ resonance provided that the first-order harmonics of the upper fundamental mode, or a subharmonic of the lower fundamental mode, is observed. We illustrate the underlying mass-spin relation in Figure 2. It is considered that $\nu_{\scriptscriptstyle\!U}$ equals $\nu_{\theta}$ and $\nu_{\scriptscriptstyle\!L}$ equals $\nu_{\theta}-\nu_{\mathrm{r}}$ (ER3:2 model).

5 Discussion and conclusions

If the observed BH HF QPOs can be described by a model similar to some of the SGO models, XMMUJ134736.6+173403  likely represents the most massive BH source with commensurable QPO frequencies. Based on the mass estimation implied by the ISCO limit, we can establish $M\approx 10^{9}M_{\sun}$ as the current observational upper limit on an AGN BH mass inferred from QPOs. This limit is valid as long as the real oscillatory frequencies of the disk are not much higher than the Keplerian frequency. This case is nevertheless not very likely to occur (e.g., Straub & Šrámková, 2009). We note that consideration of miscellaneous combinational frequencies within the QPO models may lead to obtaining observable frequencies that are (although not much) higher than the Keplerian limit.

5.1 Local versus global models, resonances

The problem of determination of the right qpo model exceeds the scope of this Letter. It is nevertheless worthwhile to note that HF QPOs commonly found in Seyefert 1 and Seyefert 2 AGNs pose a challenge for physical concepts that have been proposed to explain modulation of the observed flux within the QPO models.

Seyfert 1 AGNs are expected to have accretion disks with small inclination angles (i.e., disks that are nearly parallel to the observer’s sky). Models that include orbiting accreted inhomogeneities assume a rather small extent of the modulated region of the disk, which implies a low ratio of the modulated to total disk flux. This ratio often approaches zero for the nearly face-on view (e.g, Schnittman & Bertschinger, 2004). Either somewhat non-equatorial orbits or a very powerful enhancement of the inhomogeneities are therefore required in order to explain the observed source variability (e.g., tidal disruption events). In principle, models that introduce a collective motion of the accreted fluid allow for a high spatial extent of the modulated accretion region (e.g., Wagoner et al., 2001; Bursa et al., 2004; Schnittman & Rezzolla, 2006; Ingram & Done, 2010; Bakala et al., 2014). In this sense, this group of models that include the diskoseismic models and concepts that incorporate oscillating tori or hot inner accretion flow are favored.

Regarding the latter class of models, we suggest that it will be crucial for the falsification of models that are based on parametric resonances between disk oscillation modes whether only single QPO with a frequency that completes the harmonic sequence $3\!:\!2\!:\!1$ is found in this source, or if a new different pair of QPOs with frequencies in the $3:\!2$ ratio is found. The former case would agree with the prediction of the $3\!:\!2$ epicyclic resonance (ER3:2) model, which explains the observed 3:1 frequency ratio by means of the combinational frequencies. For both corotating and counterrotating disks, the corresponding mass-spin relation shown in Figure 2 (the dashed red-blue line) can be matched by a simple approximative quadratic relation, ${M\approx(5a^{2}+8a+8)\times 10^{7}M_{\sun}}$.

5.2 Quantification of results common to SGO models

Figure 2b shows that the ER3:2 model prediction coincides for $a=0$ with the predictions of several other models such as the RP model. On a large scale of $M,$ all the considered models provide similar mass-spin relations. Following the quadratic approximation, we may write

[TABLE]

where the upper limit corresponds to models that require a high BH mass (e.g., the ER model) and the lower limit corresponds to models that require a low BH mass (e.g., the WD model). The boundaries of relation (5) for $a\in[-1,\leavevmode\nobreak\ 1]$ are denoted by spots in Figure 2b. This figure shows that relation (5) describes the predictions of the whole group of the considered SGO models well.

Acknowledgments

We would like to acknowledge the Czech Science Foundation grant No. 17-16287S, the INTER-EXCELLENCE project No. LTI17018 that supports the collaboration between the Silesian University in Opava and the Astronomical Institute in Prague, and internal grants of the Silesian University, SGS/14,15/2016. We thank the anonymous referee for their comments and suggestions that have greatly helped us to improve the paper.