Magnetized Accretion Disks with Outflows for Changing-look AGNs

Changing-look active galactic nuclei (CL-AGNs) challenge the standard accretion theory, owing to their rapid variability. Recent numerical simulations have shown that, for the sub-Eddington accretion case, the disk is magnetic pressure–dominated, thermally stable, and geometrically thicker than the standard disk. In addition, outflows were found in the simulations. Observationally, high blueshifted velocities absorption lines indicate that outflows exist in AGNs. In this work, based on the simulation results, we investigate the magnetic pressure–dominated disk, and find that the accretion timescale is significantly shorter than that of the standard thin disk. However, such a timescale is still longer than that of the CL-AGNs. Moreover, if the role of outflows is taken into account, then the accretion timescale can be shortened even further. By the detailed comparison of the theoretical accretion timescale with the observations, we propose that the magnetic pressure–dominated disk incorporating outflows can be responsible for the rapid variability of CL-AGNs.


INTRODUCTION
Active galactic nuclei (AGNs) have strong flux variations from infrared to X-rays and are typically classified into two types: type 1, which shows broad emission lines, and type 2, which does not.According to the AGNs unified model, these two types differ in their orientations with respect to the observer's line of sight to the dust torus surrounding the AGNs (e.g., Antonucci 1993).However, an increasing number of multi-wavelength observations (including X-ray (e.g., Krumpe et al. 2017;Noda & Done 2018;Parker et al. 2019;Jana et al. 2021), infrared (e.g., Sheng et al. 2017;Yang et al. 2018;Kokubo & Minezaki 2020), and radio (e.g., Koay et al. 2016;Panessa et al. 2019;Yang et al. 2021;Lyu et al. 2021)) have identified a sub-class of AGNs referred to as changing-look AGNs (CL-AGNs), which are capable of transforming from type 1 to type 2 or vice versa over relatively short timescales.Due to the difficulties in explaining the timescale of CL-AGNs using the standard model (Shakura & Sunyaev 1973, hereafter SSD), it has become an increasingly interesting topic.
Several schemes have been proposed to explain the contradictions between theory and observation.The obscuration mechanism, based on the AGNs unified model (Goodrich & Miller 1989) where obstructions enter or leave our line of sight, can cause the appearance or disappearance of broad emission lines.However, this scheme is ruled out due to significant variation of CL-AGNs' mid-infrared emissions and the absence of polarization or weak polarization (e.g., MacLeod et al. 2016;Sheng et al. 2017;Hutsemékers et al. 2019).Another candidate is tidal disruption events (e.g., Merloni et al. 2015;Runnoe et al. 2016;Ricci et al. 2020), which predicts a change in the light curve with t −5/3 .Although it explains some cases, most observations do not support this model (Sheng et al. 2020).It is generally accepted that these intrinsic changes from the accretion disk and variations in the accretion rate dominate this phenomenon.However, the viscosity timescale of the standard disk, which varies with the accretion rate, is much larger than that of the CL-AGNs (e.g., Stern et al. 2018;Ricci & Trakhtenbrot 2022).Therefore, new models or modifications to the standard model have been proposed to address the issue.The role of the magnetic field in accretion disks has been studied and applied to the CL-AGNs (e.g., Dexter & Begelman 2019;Feng et al. 2021;Scepi et al. 2021Scepi et al. , 2023)).The limit cycle of magnetically dominated disks or disk instabilities located in the narrow region between the inner ADAF and the outer standard disk is used to explain the repeating CL-AGNs (e.g., Sniegowska et Wu & Gu al. 2020;Pan et al. 2021;Kaur et al. 2022).Wang & Bon (2020) demonstrated that tidal interaction between disks in supermassive black hole binaries can also trigger CL-AGNs.
The standard thin disk model has served as the foundation for understanding various astronomical phenomena since its proposal.However, it still faces challenges such as observations of cataclysmic variables and low-mass Xrays binaries, where the disk appears to be geometrically thicker than expected (see Begelman & Pringle (2007) for details).Most recently, numerical simulations (Huang et al. 2023) have shown that, for the sub-Eddington accretion case, the disk is dominated by the magnetic pressure and exhibits distinct properties from the standard thin disk.The magnetic pressure-dominated accretion disk is thermally stable, geometrically thicker than the standard thin disk, and has flatter effective temperature profiles.In addition, outflows were found in the simulations.These features are consistent with Begelman & Pringle (2007)'s magnetically dominated viscous disks and some previous numerical simulations (e.g., Salvesen et al. 2016;Sadowski 2016;Jiang et al. 2019) except for outflow.The twodimensional radiation-magnetohydrodynamic numerical simulation by Ohsuga & Mineshige (2011) likewise discovers outflows driven by magnetic pressure on standard-type thin disks.
Highly blueshifted absorption lines are frequently detected in the UV and X-ray spectra of AGNs (e.g., Tombesi et al. 2010;Serafinelli et al. 2019), indicating the presence of high-velocity winds/outflows in these systems.These outflows, particularly ultrafast ones, are presumed to originate in the inner regions of accretion disks and could be powered by the magnetic field (e.g., Yang et al. 2021;Wang et al. 2022).Recent observations have also validated the existence of outflows in CL-AGNs, such as MrK590 (Gupta et al. 2015) and NGC1566 (Jana et al. 2021).In this study, we utilized the results of numerical simulations to elucidate the CL-AGNs phenomenon by considering the impact of outflows.The remainder is organized as fellows.In Section 2, we provide basic physics and equations for our model.Section 3 shows numerical results and analyses.Section 4 presents the conclusions and discussion.

BASIC EQUATIONS
For simplicity, a steady-state axisymmetric magnetized accretion disk with outflows is considered under the wellknown pseudo-Newtonian potential Φ = −GM BH /(R − R g ) (Paczyńsky & Wiita 1980), where M BH is the mass of the black hole and R g = 2GM BH /c 2 is the Schwarzschild radius.Therefore, we assume the angular velocity is Keplerian, i.e., Ω = Ω In this work, the magnetic pressure P m = B 2 /8π, B is the mean mid-plane magnetic field strength.The gas pressure P g = ρk B T c /µm p , the radiation pressure P r = aT 4 c /3, where T c is the mid-plane temperature, ρ is the midplane density, µ = 0.617 is the mean molecular weight, k B is the Boltzmann constant, and a is the radiation constant.We adopt Begelman & Pringle (2007) in assuming that the magneto-rotational instability-riven dynamo could amplify the turbulent (toroidal) field until saturation.The saturation is determined by the Alfvén speed V A ≡ P m /ρ and roughly equals the geometric means of the Keplerian velocity V K and the gas sound speed c g = P g /ρ, i.e., where we have introduced a parameter λ, which represents the strength of the turbulent (toroidal) magnetic field, and λ is in the range [0, 1].The condition λ = 1 corresponds to the saturation case of the turbulent (toroidal) magnetic field and has been observed by numerical simulations (e.g., Huang et al. 2023).Due to the additional magnetic pressure, the total pressure P = P m + P g + P r , the sound speed is written as Numerical simulations have demonstrated the above magnetically dominated viscous disks (e.g., Salvesen et al. 2016;Sadowski 2016;Jiang et al. 2019;Huang et al. 2023).The model shows a geometrically thicker disk and higher color temperatures, which would be stable against thermal and viscous instabilities.In addition, outflows were found in the simulations (Huang et al. 2023).The two-dimensional radiation-magnetohydrodynamic numerical simulation of Ohsuga & Mineshige (2011) also found outflows driven by magnetic pressure on standard-type thin disks (in their numerical simulation, the toroidal component of the magnetic field is dominant, i.e., |B φ |/B p > 1).Taking inspiration from the simulations above, we developed a new outflow model.The fundamental equations describing accretion disks with outflows employed in this study can be found in our previous work (Wu et al. 2022).The continuity equation is where V R is the inflow speed, which is defined to be negative when the flow is inward, and Σ is the surface density The scaleheight can be redefined based on the above analysis Similarly, we redefine the kinematic viscosity as where α is the viscosity parameter.In principle, α is related to the magnetization of the disk (Salvesen et al. 2016).Thus, for high viscosity, such as α = 0.1, the parameter λ may not be negligible.The relationship between α and the magnetization of the disk is beyond the scope of this paper.For simplicity, we take both α and λ as free parameters in this work.
Due to the effect of outflow, we assume that the mass accretion rate varies with the radius (Ferreira & Pelletier 1995;Blandford & Begelman 1999) Ṁ = Ṁout (R/R out ) p , where Ṁout is the mass accretion rate at the outer boundary R out , in the present investigation, we set p as a free parameter.
The azimuthal momentum equation is reduced to the algebraic form (Wu et al. 2022): where . l is the parameter that characterizes the angular momentum taken away by the outflow (Knigge 1999).In this study, we set l > 1 due to the presence of a strong magnetic field.
The energy equation is written as where Q vis and Q rad are the viscous heating rate and the radiative cooling rate, respectively.Their expressions are as follows (Wu et al. 2022), where τ = κ es ρH represents the optical depth, with κ es = 0.34 cm 2 g −1 being the electric scattering opacity.In the inner region of the accretion disk, electric scattering opacity dominates over free-free absorption, leading the latter to be neglected when calculating the disk's structure.
The equation of state concerning the sound speed c s is Considering the viscosity prescription, we have By solving the ten equations, Equations (1-8), and (11-12), for the ten variables V A , P , c g , c s , H, Σ, V R , ν, ρ and T c with given parameters M BH , α, Ṁ , λ, and l.We obtain the accretion timescale t acc = −R/V R .In the following calculations, we define as ṀEdd = 16L Edd /c 2 = 64πGM BH /κ es c, and fix R out = 50R g , where a typical region of the optical/UV emission.

NUMERICAL RESULTS
Before detailed numerical calculations, we do some analysis of the basic characteristic timescales of the accretion disk, including the dynamical, thermal, and viscous timescales.Their expressions are as follows (Frank et al. 1985;Kato et al. 2008), Regarding the observed timescales of CL-AGNs, t dyn and t th are too short, while t vis is too long (e.g., Stern et al. 2018;Ricci & Trakhtenbrot 2022).From Equations ( 14) and ( 15), it is seen that t th ∝ (H SSD /R) 0 and t vis ∝ (H SSD /R) −2 .Some researchers (e.g., Stern et al. 2018;Ricci & Trakhtenbrot 2022) have suggested that setting t ∝ (H SSD /R) −1 may help address these inconsistencies.When outflows are not considered (i.e., l = 0, p = 0), this corresponds to the case of magnetically dominated viscous disks.By using Equations ( 3) and ( 6), Equation ( 7) can reduce to then, the accretion timescale For λ = 1, V A ≫ c g , using Equations (1-2) and ( 5), we have c s ≈ V A , then According to the above analysis, we have shown that t acc ∝ (H/R) −2 ∼ (H SSD /R) −1 .However, based on the following numerical calculations, we find that t acc ∝ (H SSD /R) −1 cannot fully explain the rapid variability of CL-AGNs.On the other hand, we propose that outflows play an important role in explaining the phenomenon of CL-AGNs.
Figure 1 shows the relationship of outflow parameter p and the strength of the turbulent (toroidal) magnetic field λ for a given timescale when the L bol /L Edd = 0.05, where M BH = 10 8 M ⊙ , α = 0.1, and L bol (= 2πRQ rad dR) is the bolometric luminosity.The dashed line represents (l = √ 2) and the solid line corresponds to (l = 2), We would like to stress that, the scaleheight of the disk in our model (magnetically dominated viscous disks with outflows), is significantly larger than that of a standard thin disk (see Appendix A for the derivation of scaleheight).Therefore, it is possible to further shorten the accretion timescale.The blue, red, and black lines represent 20, 10, and 5 years, respectively.For comparison, we calculated the accretion timescale of the standard thin disks (λ = 0, l = 0, p = 0, note that for standard thin disks: λ ≪ 1.In this work, for simplicity, we followed the conventional approach λ = 0 for standard thin disks, as described in some classic literature, e.g., Eq.(5.36) of Frank et al. (2002) and Eq.(3.41) of Kato et al. (2008)) and the magnetically dominated viscous disks (0 < λ ≤ 1, l = 0, p = 0) under the same conditions (α = 0.1, L bol /L Edd = 0.05): t acc ∼ 3.35 × 10 3 year and t acc ≥ 63.36 year.Our calculations demonstrate that if the role of outflows is considered in the magnetic pressure-dominated disk, the accretion timescale can be shortened to several decades or a few years.
To further understand the impact of the outflow, we adopt M BH = 10 8 M ⊙ and fix λ = 1 based on simulation results (Huang et al. 2023).We set l = 2, which corresponds to the case that significant angular momentum is carried away by outflows.The range of the parameter p is constrained by the given values of λ and l.In our case, the range 0 < p ≤ 0.16 corresponds to λ = 1 and l = 2.In addition, by varying the viscosity parameter α (0.01 ∼ 0.3) and the range of L bol (0.001 ≤ L bol /L Edd ≤ 0.1), we therefore obtain a decade range in accretion timescale in Figure 2. The blue region is the accretion timescale of standard thin disks, the green region is the accretion timescale of magnetically dominated viscous disks, and the yellow region takes into account the role of outflows.Black dots denote the high state (turn-on) of CL-AGNs discovered by Jin et al. (2022) via SDSS, with physical quantities listed in Table 1. Figure 2 shows that the timescale of magnetically dominated viscous disks is significantly shorter than that of standard thin disks, but still longer than the observed timescale.However, when the angular momentum is carried away by the outflow, the radial velocity of the accretion disk can be significantly increased, as shown in Equation (A2).Thus, magnetically dominated viscous disks with outflows can further shorten the accretion timescale and satisfy the rapid variability observed in CL-AGNs, as shown by the yellow region in Figure 2.
Based on the preceding analysis, our model predicts a timescale that can be as short as several years, which aligns with the observed timescale of CL-AGNs.However, this result is based on a fixed mass assumption and may not accurately reflect timescales for CL-AGNs with different masses.Figure 3 illustrates the variations in accretion time with mass.Similar to Figure 2, we choose λ = 1, l = 2, and 0 < p ≤ 0.16.By varying α (0.01 ∼ 0.3) and L bol (0.001 ≤ L bol /L Edd ≤ 0.1), as well as a given range for M BH (10 6 M ⊙ ≤ M BH ≤ 10 9 M ⊙ ), we therefore obtain a decade range in accretion timescale.As depicted in Figure 3, while the magnetically dominated viscous disks model can satisfy CL-AGNs with masses less than 10 7.5 M ⊙ , it is ineffective for higher masses.In contrast, our model can accommodate all sources' timescales.Therefore, the magnetic pressure-dominated disk model, with outflows taken into account, appears to be particularly promising for resolving the inconsistencies between observed timescales and disk theory.

CONCLUSIONS AND DISCUSSION
In this work, following the recent simulation results of Huang et al. (2023), we have investigated the accretion timescale of the magnetic pressure-dominated accretion disk.Our results show that, even the timescale is significantly shorter than that of the standard thin disk, it is still beyond the rapid variability of CL-AGNs.Furthermore, by taking into account the role of outflows, which were generally found in simulations and observations, we have shown that the accretion timescale can be even shortened.By the detailed comparison of the theoretical accretion timescale with the observations, we propose that the magnetic pressure-dominated disk with outflows can be responsible for the rapid variability of CL-AGNs.
The effects of magnetic fields on accretion disks have been extensively studied, and some early research (Ferreira & Pelletier 1995;Ferreira 1997;Jacquemin-Ide et al. 2019) found that magnetic winds can significantly decrease the accretion timescale.This finding has been confirmed by numerical simulations (Jacquemin-Ide et al. 2021;Scepi et al. 2023) and has been used to explain observed phenomena such as YSOs (Combet & Ferreira 2008;Martel & Lesur 2022) and X-ray binaries (Ferreira et al. 2006;Marcel et al. 2018).Full MHD models (Ferreira 1997) relate various parameters, such as p with disk magnetization, α with disk magnetization, and l with p.However, these models do not fully consider the role of magnetic (turbulence) pressure.We would point out that our model includes magnetic pressure and utilizes a simple parameterization, and it is worth further investigation on the relationship among these parameters.
Quasi-periodic eruptions (QPEs) are a newly discovered type of changing-look phenomenon characterized by rapid and extreme X-ray bursts with mass M BH ∼ 10 5−6 M ⊙ , and periods as short as a few hours (e.g., Miniutti et al. 2019;Giustini et al. 2020;Chakraborty et al. 2021;Arcodia et al. 2021).However, these short timescales cannot be explained even at R = 50R g and α = 0.3 using our present one-zone model.In order to interpret QPEs, a thicker disk may be necessary.Begelman et al. (2015) proposed a vertically stratified model with a larger characteristic scale height H ∼ R which agrees with Bai & Stone (2013)'s vertically stratified, local shearing box ideal MHD.Subsequent numerical simulations (Salvesen et al. 2016) have also confirmed this model.In addition, Bai & Stone (2013)   shows that a magnetocentrifugal mechanism could excite a strong outflow.Therefore, our future work will focus on a new outflow model based on these findings to explain the QPE phenomenon.
We thank Zhi-Xiang Zhang for helpful discussions, and the anonymous referee for constructive suggestions that improved the paper.This work was supported by the National Natural Science Foundation of China under grants 11925301, 12033006, and 12221003.

APPENDIX
A. EXPRESSION OF SCALEHEIGHT In this Appendix, we will derive the detailed expression of scaleheight, which is mentioned in Section 3.For λ = 1, V A ≫ c g , according to Equations (1-2) and ( 5), the scaleheight takes the form: The meaning of the colors and symbols is the same as in Figure 2.For the standard thin disks, the accretion timescale is much larger than the observed timescale.Magnetically dominated viscous disks only explain part of the CL-AGNs.However, our model (magnetically dominated viscous disks with outflows) can naturally explain all the observed sources.

Figure 1 .
Figure 1.Timescales are shown in the {p, λ}-plane for a given L bol /L Edd = 0.05.The dashed line represents (l =√ 2) and the solid line corresponds to (l = 2).Each of the color lines corresponds to a specific value on the timescale.

Figure 3 .
Figure 3.The variations of the accretion timescale with mass under the λ = 1, l = 2, α = 0.01 ∼ 0.3, L bol /L Edd = 0.001 ∼ 0.1.The meaning of the colors and symbols is the same as in Figure2.For the standard thin disks, the accretion timescale is much larger than the observed timescale.Magnetically dominated viscous disks only explain part of the CL-AGNs.However, our model (magnetically dominated viscous disks with outflows) can naturally explain all the observed sources.
The accretion timescale of accretion disks at MBH = 10 8 M⊙ for different α and L bol .The regions of color correspond to the cases for the standard thin disks (blue), magnetically dominated viscous disks (green), and our model (magnetically dominated viscous disks with outflows, yellow).The black dots are the high state (turn-on) of CL-AGNs in Table1.

Table 1 .
Information about sources.