Stochastic Variability of Binary Accretion

技术摘要：双星吸积的随机变率

问题陈述
超大质量黑洞（SMBH）双星是等级化星系形成的预测结果，但目前在观测上仍难以捉摸。虽然研究者们一直试图通过接近双星轨道频率的准周期通量调制来寻找其主要标志，但这类搜索受到活跃星系核（AGN）本质上的随机性性质以及观测基线限制的干扰，后者排除了绝大多数具有长轨道周期的亚秒差距双星。此外，由于流体动力学模拟所预测的脉冲形状是非正弦型的，标准的周期图恢复率会显著下降。因此，有必要识别出在连续谱功率谱密度（PSD）中不同于单星 AGN 典型红噪声变率的时间特征，特别是在短于双星轨道周期的尺度上的随机变率。

方法论
作者使用 Sailfish 代码进行高分辨率 2D 流体动力学模拟，以模拟质量不等（$q = M_2/M_1 = 0.2$）且处于圆轨道上的双星及其周围环双星吸积盘中的气体吸积。模拟采用局部等温状态方程，轨道马赫数 $M=10$，并使用 Shakura-Sunyaev $\alpha$ 处方描述粘性（$\alpha = 0.01$）。

方法论的一个关键方面是系统地改变“汇”（sink）的大小参数 $s$，该参数决定了组分微型吸积盘（minidisks）的内截断半径（$r_{sink} \propto s \cdot a$）。参数 $s$ 从 0.05 变化到 0.40。作者区分了进入微型吸积盘的质量通量（“驱动器”）与向黑洞汇注入的质量吸积率（“滤波器”）。模拟经过一个烧入阶段（2900 个轨道周期）后，进入采样阶段（100 个轨道周期），并采用高分辨率（$\Delta \ln r = 0.0025$）进行采样，其时间序列采样频率比轨道频率高出约 $10^3$ 倍，以解析高于轨道频率的频率成分。

关键结果

1. 连续谱 PSD 形态： 在小汇尺寸极限下（$s \lesssim 0.1$），吸积率的功率谱密度收敛为一种普适的折断幂律形状。PSD 在低频处为平坦（白噪声），在窄带内过渡到 $-2$的斜率，并在高频处陡峭化为$-4$ 的斜率。
2. 折断频率： 向陡峭 $-4$斜率的转变发生在折断频率$f_d \approx 5 f_{bin}$处（$f_{bin}$为双星轨道频率）。该频率对应于次级黑洞微型吸积盘外边缘的局部开普勒频率（$r_{out} \approx 0.18\text{--}0.20 a$）。
3. 对汇尺寸的依赖性： 对于大尺寸汇（$s \gtrsim 0.2$），PSD 并不收敛于上述普适形状。相反，微型吸积盘起到的缓冲作用不足，导致过剩的高频功率到达“汇”，表现为吸积率尖峰。作者认为，当汇半径超过黑洞物理上的最内稳定圆轨道（ISCO）时，这些尖峰属于数值伪影。
4. 驱动器-滤波器机制： 这种变率可以用两阶段随机过程来解释：
   • 驱动器： 环双星流以类似于阻尼随机游走（Ornstein-Uhlenbeck 过程）的方式将气体输送到微型吸积盘的外边缘，其相关时间与双星轨道周期相当。
   • 滤波器： 微型吸积盘充当低通滤波器，阻尼掉快于其外边缘动力学时间尺度的波动。这种滤波产生了额外的 $f^{-2}$ 陡峭化，从而导致在高频处总斜率为 $f^{-4}$。
5. 粘性无关性： PSD 的形状（包括折断频率和高频斜率）在跨越一个数量级的粘性参数 $\alpha$ 变化时保持不变，这证实了该折断是由动力学时间尺度而非粘性扩散决定的。

意义与主张
本文声称，具有高频 $-4$ 斜率的折断幂律 PSD 可作为亚秒差距 SMBH 双星的一个潜在连续谱标志。这种形态不同于通常应用于单星 AGN 的标准阻尼随机游走模型（后者通常预测接近 $-2$ 的斜率）。

作者提出，如果一个 AGN 已被怀疑存在双星（例如通过周期性候选者或宽发射线速度偏移），那么在 $5 f_{bin}$ 附近检测到连续谱折断及 $-4$ 的陡峭斜率将支持双星假设。反之，若不存在此类折断，则会削弱双星假设。

然而，本文对于“所有 AGN 的变率都源于广泛的双星系统”这一更广泛的推测持谨慎态度。虽然模拟的 PSD 与许多 AGN 中观察到的随机变率相似，但作者指出，在 AGN 中观察到的折断频率（对应于数年到数十年的轨道周期）会将这些系统置于脉冲星计时阵列可探测的引力波频率段内。鉴于当前对随机引力波背景的约束，不太可能存在大量发光 AGN 在这些特定周期下作为双星存在，因为这会过度产生观测到的背景信号。因此，尽管该机制是针对特定候选者的有效标志，但它可能无法解释整个 AGN 群体的变率。

最后，本研究为模拟方法建立了一个数值准则：为了准确捕捉双星吸积的物理变率，模拟必须采用足够小的汇尺寸（$s \lesssim 0.1$），以解析微型吸积盘的缓冲效应，避免由过大尺寸汇产生的伪高频噪声。