Intermediate-Mass Black Holes: The Essential Population to Explore the Unified Model for Accretion and Ejection Processes

We study radio and X-ray emissions from IMBHs and explore the unified model for accretion and ejection processes. The radio band survey of IMBH (candidate) hosted galaxies indicates that only a small fraction ($\sim$0.6\%) of them are radio-band active. In addition, very long baseline interferometry observations reveal parsec-scale radio emission of IMBHs, further resulting in a lower fraction of actively ejecting objects (radio emission is produced by IMBHs other than hosts), which is consistent with a long quiescent state in the evolution cycle of IMBHs. Most (75\%, i.e., 3 out of 4 samples according to a recent mini-survey) of the radio-emitting IMBHs are associated with radio relics and there is also evidence of dual radio blobs from episodic ejecting phases. Taking the radio emission and the corresponding core X-ray emission of IMBH, we confirm a universal fundamental plane relation (FMP) of black hole activity. Furthermore, state transitions can be inferred by comparing a few cases in XRBs and IMBHs in FMP, i.e., both radio luminosity and emission regions evolve along these state transitions. These signatures and evidence suggest an analogy among all kinds of accretion systems which span from stellar mass to supermassive black holes, hinting at unified accretion and ejection physics. To validate the unified model, we explore the correlation between the scale of outflows (corresponding to ejection powers) and the masses of central engines; it shows that the largest scale of outflows $\hat{LS}_\mathrm{out}$ follows a power-law correlation with the masses of accretors $M_\mathrm{core}$, i.e., $\log{\hat{LS}_\mathrm{out}} = (0.73\pm0.01)\log{M_\mathrm{core}} - (3.34\pm0.10)$. In conclusion, this work provides evidence to support the claim that the ejection (and accretion) process behaves as scale-invariant and their power is regulated by the masses of accretors.


Introduction
Two types of astrophysical black holes, stellar-mass black holes (SBHs) and supermassive black holes (SMBHs), are widely identified in the Universe through their observational signatures. The death of massive stars will leave behind SBHs as fossil [1], which have the masses 3-100 M . SBHs are abundant in our Galaxy as the X-ray binary systems (XRBs, which accrete matters from a companion star). On the other hand, SMBHs (≥10 6 M ) are generally identified at the center of massive galaxies (they co-evolve with bulges and hosts) [2]. Remarkably, SMBHs have an essential role in regulating the evolution of their host galaxies, which conduce to the correlation between SMBHs and their host properties, e.g., the correlation between the SMBH mass and the bulge velocity dispersion (M − σ). In particular, accreting SMBHs from active galactic nuclei (AGNs) can maintain strong feedback to the hosts. The black holes with The analogy between XRBs and AGNs implies a scale-invariant accretion and ejection process among accreting black holes. This is motivated by a radio-to-X-ray correlation (i.e., the disk-jet coupling) in the XRB systems [44]; then Merloni et al. [39] shows that the radio-to-X-ray correlation is applicable in both XRBs and AGNs and is regulated by the black hole masses, i.e., the fundamental plane of black hole activity. Unfortunately, the original 'fundamental plane' has a break between SBHs and SMBHs due to the lack of IMBH samples. Subsequent studies, e.g., Greene et al. [41], Gültekin et al. [40] and Bariuan et al. [45] investigated the radio properties of the low-mass RQ-AGNs using the arcsec-scale resolution Very Large Array (VLA) observations and found that they deviate from Merloni's fundamental plane relation. However, our recent studies [46] show that the inconsistency between the fundamental plane and the observations for low-mass AGNs may simply correspond to a deficit of resolution.
Radio emissions from radio-loud AGNs (RL-AGNs; defined as R > 10) are primarily from synchrotron mechanisms, where the electrons are accelerated to relativistic energies in jets (e.g., [47,48]). Jets can exist for a long time (up to mega years) and the physical scales of the relativistic jets can extend from less than a parsec to megaparsec scales. Jets are the dominating radio-emitting source in RL-AGNs; however, most AGNs are radio-quiet [49]. On the weak radio emissions of radio-quiet AGNs, several mechanisms are responsible for the radio origin (see [50] and references therein). The possible radio emissions of RQ-AGNs include lower power jets (e.g., [51][52][53]), free-free emission from photo-ionized gas in the circum-nuclear region [54], star-forming regions (e.g., [52]), accretion disk-corona activity (e.g., [43,[55][56][57]) and sub-relativistic wide-angled winds (e.g., [53,58,59]). Hereinafter, we term the jets and the windlike ejections as outflows. Therefore, in constraining the radio emission of radio-quiet AGNs, high-resolution (e.g., the very long baseline interferometry) radio observations are essential.
As one of the highest-resolution approaches, the very long baseline interferometry (VLBI) method at the radio band provides the milli-arcsec scale resolutions at GHz bands. The highresolution VLBI observations of dwarf galaxies can reveal pc-scale radio-emitting structures of their central engines, which can be cores, core-jets or jet-knots. Therefore, VLBI observations can directly probe the jet powered by potentially accreting IMBHs. At present, the VLBI observations of IMBH candidates in dwarf galaxies are only limited to a few cases: NGC 4293 [13], NGC 4395 [60], Henize 2-10 [61], NGC 404 [62], RGG 9 [63] and the four samples from the mini-survey of IMBHs by Yang et al. [46]. The radio observation with the High Sensitivity Array (HSA) 1.4 GHz of NGC 4395 (R X = −5, ref. [64]) reveals an elongated sub-parsec scale structure, which was explained as a core-jet [60]; however, the recent study with European VLBI Network (EVN) 5 GHz high-resolution observations of NGC 4395 by Yang et al. [65] failed to detect a radio-active nucleus at the optical (Gaia) position. Alternatively, a relatively low-resolution Very Large Array (VLA) 15 GHz observation detected the radio emission that coincides with the Gaia position. The above observations are consistent with the scenario of multiple radio lobes from episodic ejections. On the other hand, the EVN 5 GHz observation of NGC 404 [62] failed to detect radio emission from its nucleus, while the large-scale (parsec-scale) radio emission was detected with the VLA 1.4 GHz observations [66]. Furthermore, the Chandra X-ray observation of NGC 404 indicated a steady source of log R X < −3.8 in the nucleus.
In this work, we explore the universal accretion and ejection processes spanning from stellar mass to supermassive black holes. In particular, this work involves the largest sample (at present) of radio-emitting IMBHs to unveil the vacuum region between SBHs and SMBHs. Throughout this work, we adopt the standard ΛCDM cosmology with a Hubble constant H 0 = 70 km s −1 Mpc −1 and matter density and dark energy density parameters Ω Λ = 0.73, Ω m = 0.27, respectively.

The Universal Picture for Black Hole Accretion and Ejection
Both SBHs and SMBHs can maintain accretion and ejection processes, which form Xray binaries (XRBs) and active galactic nuclei (AGNs), respectively. The accretion discs and associated coronae and jets in both XRBs and AGNs produce multi-band emissions, especially at X-ray and radio bands. Though with different environments and size scales, observations and theoretical studies reveal similar accretion and ejection behaviors, e.g., the evolution of accretion flows, the structures of accretion discs and the features of the ejection process. In individual sources, the accretion and ejection processes evolve as the fluctuation of accretion rates in terms of Eddington ratios (see [67,68]).
On the activity of SBHs, monitoring and simultaneous observations at radio and X-ray bands reveal a clear evolution cycle among them. They can be well characterized as several accretion and ejection states through X-ray and radio properties [69]. The X-ray emissions in accreting systems are thought to be primarily and generally from the accretion disk and corona system (and jets as well), while radio emissions are thought to be dominated by jets. The association between X-ray and radio properties in X-ray binary indicates a strong coupling between accretion discs and jets [67]. The jet-disk coupling in X-ray binaries can be observed with convenient time scales (months to years), which is long expected to be a common feature in both stellar mass and supermassive black holes. Furthermore, observations of X-ray binaries show that the accretion flow and jet production is correlated with the accretion rate/Eddington ratio. The fluctuations of the accretion rate will essentially trigger the accretion state transition [70].
Based on the monitoring observations of X-ray binaries, we now have a picture of state transition and evolution: (a) the accretion with a low Eddington ratio (typically below a few percent of the Eddington luminosity) is characterized by a 'low/hard' X-ray state; it is associated with compact radio emission and identified to be steady and short jets; (b) with the increase of accretion rates, the X-ray spectrum becomes dominated by the soft emission and the radio emission dramatically drops to an undetectable level. The Eddington ratio of the soft state is substantially close to the Eddington limit; (c) the transition from the 'low/hard' state to the 'soft state' is generally associated with an unstable 'very high state', where the Eddington ratio will be up to and over one, i.e., near-and super-Eddington rates. Episodic jet blobs can be triggered as the X-ray binaries in the 'very high state', while mechanisms are relatively unclear due to the short duration. The steady jets in quiescent and low/hard state are generally believed to be launched by extracting the energy of either an accretion disk (the 'BP' jets, [71]) or the spin of the central compact object (the 'BZ' jets, [72]), while episodic jets are alternatively explained by other models (e.g., the instability of accretion disks [73]).
Several works indicate that the accretion in the transition state can be described as a 'slim disk' [74], which supports the near and moderately super-Eddington accretion. As the super-Eddington accretion has an essential role in ultraluminous X-ray sources and high redshift quasars, the transition state is also named an 'ultraluminous state' [75]. Due to the bright episodic jet eruption (which is brighter than steady jets), the radio emissions observed in IMBHs and low-luminosity AGNs likely correspond to the ejecta of the 'ultraluminous state'. However, the understanding of the accretion and ejection properties in the 'ultraluminous state' is limited to only a few XRBs, e.g., Cygnus X-3, GRS 1915 + 105 and SS 433.
It was theoretically studied that the accretion flow in supermassive black holes in active galactic nuclei resembles the stellar mass black holes in X-ray binary systems. Therefore, the two systems, even though they have a mass range of 10 dexes, should share the same correlations. Indeed, there are several correlations that work in both stellar mass and supermassive black holes (e.g., [39,43,76]). The universal correlations among the diverse accretion systems indicate a similar physics in accretion, that naturally induces the scale-invariant accretion and ejection phenomena in both radio and X-ray. Historically, there are several correlations that are successful in unifying XRBs and AGNs: (1) the most applicable correlation among X-ray luminosity, radio luminosity and black hole mass was named the fundamental plane of black hole activity [39,77]. The fundamental plane correlation was initially identified in X-ray binaries when their accretion state is in the low/hard state, where the accretion is in low and moderate rates. As the accretion transits to the 'high/soft' state, the radio emission (thus the jet) dramatically drops while the X-ray emission remains in a high luminosity, which deviates from the fundamental plane relation. The similarity between XRBs and AGNs not only indicated that both of them satisfy the fundamental plane relation of black hole activity but also the similar radio quenching as the increase of accretion rates. At present, the fundamental plane relation has a clear gap between the stellar mass and supermassive black holes. Filling the gap would build a more fundamental and comprehensive understanding of the accretion systems; therefore, IMBHs are the key to linking the stellar mass and supermassive black holes in the fundamental plane relation. The fundamental plane relation is important as it links both accretion and ejection processes and is applicable based on X-ray and radio surveys; (2) because the fundamental plane relation of black hole activity only works well in the low/hard state, a universal correlation developed, i.e., the inverse correlation between the radio loudness and the Eddington ratio [43,[78][79][80]. The inverse correlation indicates the suppression of the jet producing as the increase of Eddington ratios, which is applicable for the radio emission produced in the high and super-Eddington XRBs and AGNs; (3) the correlation among the bolometric luminosity, black hole mass and the characteristic timescales of X-ray variability [81]; (4) the hardness-intensity diagram (HID, [82]) and the accretion rate ratio-intensity diagram (ARRID, [83,84]) for individual XRB and the statistical version, the disk-fraction luminosity diagrams [76,82]; (5) the inverse correlation between QPO frequencies and black hole masses [22,25,85]. However, evidence is still lacking that the above correlations are applicable for all accreting systems, i.e., these correlations need to be explored in the mass scale of 10 3 -10 5 M .

Radio Activities of IMBHs
Given that dwarf galaxies tend to host low-mass AGNs (or IMBHs), Reines et al. [86] conducted a radio census of 111 FIRST detected dwarf galaxies (the stellar mass range from 3 × 10 7 to 3 × 10 9 M ). Among them, only 13 were identified to be powered by AGNs, through radio diagnoses, e.g., the compact point-like morphologies and the excess radio emissions compared with those afforded by star-forming activities. This yields an AGN rate of 11% in radio-active dwarf galaxies. Only one of these 13 was cross identified in optical bands and therefore it results in 0.9% AGN fraction in dwarf galaxies simultaneously. However, the fraction may not imply a physical distribution due to the sample selection effect and the incompleteness of the sample.
Alternatively, we [46] compiled a sample of 598 IMBH candidates which were uniformly selected from the Sloan Digital Sky Survey (SDSS, [18][19][20][21]35]). Through cross-matching with the NRAO VLA Sky Survey (NVSS, [87]) and Faint Images of the Radio Sky at Twenty centimeters (FIRST, [88]) catalogs within 1 arcsec of their optical positions, Yang et al. [46] found that 36 sources (6%) have radio counterparts (with signal-to-noise ratios >9) in the FIRST survey. Furthermore, we selected four IMBH candidates (0.6%) for VLBI observations with the further criteria: (1) the estimated black hole mass should be <10 6 M , (2) the 1.4 GHz radio flux density of the FIRST survey is >2 mJy, which yields a signal-to-noise ratio of >12, and (3) objects associated with the point X-ray emissions. Yang et al. [46] successfully detected parsec-scale radio emissions in these four candidates and confirmed that the IMBH candidates are radioactive, i.e., they maintain ejection processes. In Table 1, we summarize the IMBH candidates with both radio and X-ray observations.
Comparing the optical and radio emission regions of these IMBH candidates, most of the radio-active IMBHs have radio emission regions offset from the optical nuclei. This indicates that the radio emissions in these cases are more likely relics of ejecta than radio-active nuclei. Given the small fraction of radio-active IMBHs and the lesser rate of actively ejecting ones of them, it satisfies the scenario of episodic ejections and a long evolution cycle from a quiescent state to a very high state. The episodic ejection nature in IMBHs has been further identified in a few other candidates with multiple radio-emitting blobs, e.g., NGC 4293 [13], NGC 4395 [65] and RGG9 [63].

The Black Hole Fundamental Plane of Involving IMBHs
Accretion onto compact objects follows the so-called fundamental plane relation (FMP) of black hole activity, which is among nuclear radio and X-ray luminosities and black hole masses of any accretion systems. The fundamental plane relation originates from jet-disk coupling, where the jet produces radio emission and the accretion disk and corona system produce X-ray emission and the mass of the central engine regulates both the accretion and ejection power. The correlation between radio and X-ray luminosity was initially inferred from X-ray binary systems that host stellar mass black holes, which only work in quiescent and low/hard accretion states that are associated with steady ejections [39,44].
In addition, the very high/intermediate state will substantially trigger episodic radio ejection, then subsequently the radio emission will dramatically decrease as the accretion state evolves from the very high state to the soft state. Similarly, we expect that the radio power of blobs may also be correlated with the X-ray luminosity. Indeed, several works (e.g., [39,44]) explored the radio-to-X-ray correlation in transient (very high) states, in X-ray binaries; the results indicate that the correlation may still be maintained but with high dispersion. It was already known that the fundamental plane of black hole activity can be contaminated by several mechanisms, e.g., the enhancement of the radio emission of lobes (i.e., the shock waves) when propagating through a dense medium [13,95]; the radio contamination by AGN hosts; the relativistic beaming effect of jets; the X-ray contamination by the jet [96]. Therefore, the fundamental plane relation works in both low/hard and very high/intermediate states.
The fundamental plane of black hole activity presented by Merloni et al. [39] is There are also recent updates, e.g., Saikia et al. [94] log L R = (0.48 ± 0.04) log L X + (0.79 ± 0.03) log M BH + 11.71, and Gültekin et al. [97] log M BH 10 8 M = (1.09 ± 0.10) log L R 10 38 erg s −1 + −0.59 +0. 16 −0.15 log L X 10 40 erg s −1 + (0.55 ± 0.22). (3) Here we only explore the fundamental plane from Merloni et al. [39] based on the scope of this work, i.e., the radio emission of AGNs (especially low luminosity and low-mass AGNs) should be obtained from the parsec-scale region and the comparison sample used in this work is from the same work [39]. The fundamental plane of black hole activity indicates a scale-free accretion and ejection physics in both SBHs and AGNs. Recently, Yang et al. [46] explored the fundamental plane relation in a sample of IMBH candidates with pc-scale radio emissions obtained from VLBI observations; which implies that IMBHs follow the fundamental plane of black hole activity. In Figure 1, we use the IMBH sample from Yang et al. [46] and Yang et al. [13]. There is one more IMBH candidate (J111552.01 − 000436.1) from the same observation of Yang et al. [46] but it was not detected; here, we take its upper limit. Furthermore, we involve a few IMBH candidates from the literature (see comments in Table 1). In particular, we only use the radio emission from VLBI observations (see Table 1) to reduce the contamination from hosts. However, it is hard to determine whether to take the radio emissions from lobes or core regions. Neither the radio emission from lobes nor core regions indicate the power of ejection by considering that the lobe may receive enhancement by circumnuclear medium and beaming effect and the effect of episodic ejection and state transition. To reduce the deficit (see Figure 2), we use the radio emissions close to the core region, which are accompanied by the core X-ray emissions. Specifically, we use the radio emission from the core region of NGC 4293 as the lobes of this source received a strong enhancement by circumnuclear medium (see [13]). The linear distribution trend presented in Figure 1 indicates that the accretion and ejection processes in IMBHs are similar to those in SBHs and SMBHs, i.e., there exists a unified model for the accretion and ejection process for SBHs, IMBHs and SMBHs. Obviously, involving IMBHs in the fundamental plane fills the gap between the stellar mass and supermassive black holes, i.e., low luminosity AGNs (powered by supermassive black holes) extend to the lower-left in the fundamental plane relation and approximately overlap with the highly accreting IMBHs and low-rate accreting IMBHs touch the most luminous stellar-mass black holes.  [39] and the IMBH candidates are from Yang et al. [46] and Yang et al. [13]. We will emphasize that the fundamental plane relation is devoted to depicting the correlation between the radio and X-ray luminosity in the ejection phase (excepting the canonical 'soft state'), which is relatively difficult for IMBHs and low-mass AGNs. Radio observations of these sources are generally limited to the radio emissions produced from very high states because the radio emissions in quiescent and low/hard states are likely dominated by the compact and steady jet, which is weaker than the episodic radio blobs in very high states. However, we do not track the evolution of each radio blob, but the radio luminosity of the instance ejection and associated X-ray luminosity. Again, low-resolution radio observations may collect the radio emissions from multiple blobs of one or more evolution cycles, which will obviously overestimate the radio emission from one single blob. Due to the long timescale of the evolution cycle, it is generally impossible to record the radio and X-ray emission in the very high state of individual IMBH. Fortunately, the radio blobs in IMBHs and low-mass AGNs will live for a sufficiently longer time than in XRBs and move to a distance from parsec to kilo-parsec scales. Therefore, it is applicable to take the radio emission from a recent blob along with the core X-ray emission.
In Figure 1, we also involve pc and sub-pc-scale radio constraints of two IMBH candidates, NGC 4395 and NGC 404, having low accretion rates. They have the Eddington ratios of λ Edd = 1.5 × 10 −6 [62] and 1.2 × 10 −3 [60], respectively. In this plot, the radio luminosity of NGC 404 was taken from the nuclear 7 parsec scale region through VLA A-array 5 GHz observation [66] and we estimated the 5 GHz radio luminosity of NGC 4395 via VLA A-array 15 GHz observation [94], from the nuclear 4 parsec scale region. Both the radio luminosities of the above two sources are obtained from the parsec-scale region, which is comparable to the VLBI observations of the other IMBHs presented in this plot due to the low redshift of NGC 4395 and NGC 404. Furthermore, both NGC 4395 and NGC 404 received the new measurements of black hole masses; see Woo et al. [15] and Davis et al. [89], respectively. In particular, we take the X-ray luminosity of NGC 404 from Chandra observation [62]. It should be noted here that the two sources have radio emission at parsec scales, but VLBI observations only reveal the upper limit of them. Interestingly, a sample of low luminosity AGNs has a similar radio drop from parsec to sub-parsec scales [98]. Certainly, the non-detection of sub-parsec scale radio emissions in a few low accretion rate AGNs [98] and IMBHs (NGC 4395 and NGC 404) challenges the (Merloni) fundamental plane and deeper VLBI observations are required to unveil the nature of these central engines. Taking the unified scheme for the accretion and ejection process, it was already shown that the episodic radio ejection and the high Eddington ratios in a few IMBHs [13,43,46,63,65] draw a similarity with the canonical very high/intermediate state in XRBs. On the other hand, the low Eddington ratios in a few low luminosity AGNs [98] and IMBHs (NGC 4395 and NGC 404) alternatively indicate quiescent or low/hard accretion states. Studying XRBs shows that radio emissions will drop with the transition from a very high/intermediate state to a soft state and comparing between a very high state and a quiescent state [99]. In Figure 1, we mark the monitoring observations of the microquasar GRS 1915 + 105 to illustrate the similarity and it is likely consistent with the observational signatures in NGC 4395 and NGC 404.

The Correlation between the Size Scale of Outflows and the Masses of Accretors
Several critical factors regarding the size scales of outflows include the power of the central engine and the external environment. Given that external environments vary in each object, the largest scale of outflows should be only related to the power or type of the central engines/accretors for ejecting objects with similar ages, statistically. Again, the distribution of ejecting ages directly results in the positive distribution of outflow sizes and the largest outflow size tends to be associated with the older (over a critical age) ones of them. Naturally, the largest size scale of outflow-dominated feedback relies on the power of the central engine, hence the mass of the central engine according to the unified accretion and ejection process. Therefore, a positive correlation between the largest (linear size) scale of outflowsLS out and the masses of accretors M core should be expected. Now we try to explore the M core − LS out correlation among a sample of accretors with masses ranging from 10 −2 to 10 10 M (the data will be presented in a forthcoming paper, Yang [100], along with additional samples). This sample includes young stellar objects (YSOs), X-ray binary systems (XRBs), intermediate-mass black holes (IMBHs), normal radio galaxies (RGs) and giant radio galaxies (GRGs). Obviously, theLS out − M core correlation can be identified by inspecting the upper envelope of jet sizes (the line 'best fit −1' in Figure 3). However, jets can further extend to form relics and bubbles; on the largest scale of outflows, the jet relics in GRGs and a few YSOs and XRBs may represent the largest structures of the ejected plasma (fossil plasma); in addition, the shocks driven by outflows tend to form a more elongated structure (i.e., shock-inflated bubbles) than fossil plasma (see [100]). Interestingly, the jet relics in GRGs and a few YSOs and XRBs follow the trend identified in (active) jets, which ensures we fit a linear correlation between the size of jet relics and the masses of accretors and gives the following parameters, whereLS out is in kpc and M core is in solar mass. The distribution between the linear size of outflows and the masses of accretors is presented in Figure 3. We further extend the linear correlation by −1 and 2 dex, which likely links the largest size scale of active jets and shock-inflated bubbles, respectively. This implies that the shock-inflated bubbles can reach a 3-dex longer extension than the active jets. Accordingly, a hypothesis is that the jets in IMBH may provide a few hundred kilo-parsec scale feedback, which is beyond the size of their dwarf hosts. Again, based on the size of the shock-inflated bubbles in XRBs, jets produced by XRBs will drive shocks to further blow bubbles of 2-dex larger than the jet relics. If this scaling relation holds for both IMBHs and SMBHs, then we expect to find shock-inflated bubbles of 10-100 kpc and 100 Mpc scale in IMBHs and the most massive SMBHs, respectively. Furthermore, LS out − M core correlation lacks the data of jet relics in the IMBH mass scale, which requires further exploration in observation. Again, a theoretical deduction is required as well (see a forthcoming paper, Yang [100]). The kinetic power is the key parameter to drive jet expansion, which is thought to be correlated with the black hole spin and mass [71,72,101] and the structure of the accretion disk [73,102,103]. The correlation between the jet kinetic power and black hole masses is indicated in several statistical analyses (e.g., [104][105][106]), while the observational correlation between the jet power and black hole spin is still a mystery (e.g., [107]). Except for the kinetic power, the size scale of an outflow in an individual object may be regulated by ages [95], inclination angle and interplay between jets and surrounding media [108]. These parameters insert dispersion in the LS out distribution along the mass of accretors (see, e.g., [101]). However, the largest size scale of outflow (the upper envelope in LS out − M core distribution) is likely free from those parameters due to the selection effect, i.e., the largest scale of outflow tends to be associated with the longest-lasting and the highest spin central engines. Again, a low frustration from the surrounding media is an assistance to allow the jet expansion further. Therefore, the correlation between the size scale of outflows and the masses of accretors can be considered to be intrinsic.