Why haven't the observed high-frequency QPOs been produced by (MHD) computer simulations of black hole accretion disks?

We compare some predictions of Wagoner&Tandon (2021) with the results of the hydrodynamic and magnetohydrodynamic (MHD) simulations of Reynolds&Miller (2009). It appears that the MHD simulations were not run for long enough and the numerical damping was not small enough to produce the observed high-frequency QPOs (and the g-mode seen in the hydro simulations).


INTRODUCTION
We do not consider the low-frequency QPOs, since their frequencies are not stable. As an example of a possible high-frequency QPO (HFQPO) in a black hole accretion disk, we choose the fundamental (axisymmetric) g-mode seen in the hydro simulations of RM. We shall employ the dimensionless time (g-mode phase) τ = 2πf 0 t ≈ (2.8 × 10 −2 c 3 /GM)t, where f 0 ≈ 4.5 × 10 −3 c 3 /GM is the frequency of the g-mode. The run time of the canonical MHD simulation MHD 1 of RM was τ (run) = 1.1 × 10 3 .
The quantity plotted in Fig. 2 of RM is of order the turbulent stresses (typically greater than the magnetic stresses within the eigenfunction of the mode). We have assumed that the divergence of the turbulent stress tensor is the main driver of the mode. From that Figure, it appears that the damping time of the hydro runs was τ (damp) ≈ 11. We (and RM) assume that this is due to numerical dissipation, in which case τ (damp) ≈ (grid size) −1/2 . (WT have shown that 'turbulent viscosity' does not operate on modes of oscillation.) The damping time of the MHD run was about 10 times longer. This is consistent with the fact that the grid sizes of runs HD3d 1 and MHD 1 were similar. Goldreich & Kumar (1988) showed that convective turbulence excites (drives) the p-modes in the sun. However, the origin of the damping of these modes is still unknown (references in WT). The possibilities quoted do not include turbulent viscosity.

RESULTS
We take the damping rate parameter D * introduced by WT to be equal to 1/τ (damp). The amplitude of the mode should grow until a time of order 1/(damping rate), as WT found. WT also found that if equilibrium (roughly constant energy E of the mildly nonlinear oscillations) is reached, the short time average is The (dimensionless) oscillation amplitude A(τ ) is approximately equal to the fluid displacement divided by the size of the smallest dimension of the eigenfunction (of order the thickness h of the accretion disk). The parameter η is the volume of the maximally coupled edd(ies) divided by the volume of the eigenfunction, and M is the Mach number of those eddies. Like the MHD (but not the hydro) runs of RM, our evolutions were begun with A = dA/dτ = 0. Since equilibrium was achieved by RM [τ (run) ≫ τ (damp)], we find from equation (1) [with D * = 1/τ (damp) ∼ 0.01 for the MHD run] that for representative values η ∼ M ∼ 0.1, A 2 ∼ 10 −5 . This amplitude would presumably be undetectable in the MHD simulations.
Employing the observed PSDs of black-hole binaries (BHBs) and narrow-line Seyfert 1 AGNs (NLS1s), WT found that the visibility of a mode requires that where 1/Q is the fractional width of the QPO. Most observed HFQPOs have Q 4. Also, WT verified that nonlinearity reduced the value of Q. We have assumed an axisymmetric mode of extent δr ∼ δz ∼ h(r) ∼ 10 −2 r. We note that leakage of a mode into a traveling wave (O'Neill et al. 2009) can reduce the required amplitude A. WT found that a 'blowout' (onset of the effects of the lowest order nonlinearity of the restoring force, terminating the run) occurred if the amplitude reached a value |A| = 1 (at E = 1/6), which from equation (1) occurred if D * 0.1(ηM 2 ) 2 . Therefore, in order to reach a (more observable) blowout during a run, τ (damp) = 1/D * 10 A 2 (ηM 2 ) −2 ∼ 10 7 .
3. CONCLUSIONS From these results, it would appear that in order for simulations of black hole accretion disks to produce the observed high-frequency QPOs, much smaller grid sizes and consequently much longer runs must be employed. This result applies to QPOs produced by normal modes of oscillation whose minimum extent is of order the local thickness of the accretion disk. Topics for future research then include: (1) How can this finite grid size source of damping be reduced? (2) Why is our predicted damping within black hole accretion disks with observed HFQPOs so small, and what is its origin?
We thank Marek Abramowicz and Cole Miller for comments. CT acknowledges financial support from the Stanford Physics Department Summer Undergraduate Research Program. We would like to thank Stanford University and the Stanford Research Computing Center for providing computational resources on the Sherlock cluster and support.