Search

In Dou et al. (2021), we introduced the Fundamental Formation Relation (FFR), a tight relation between specific SFR (sSFR), H$_2$ star formation efficiency (SFE$_{\rm H_2}$), and the ratio of H$_2$ to stellar mass. Here we show that atomic gas HI does not follow a similar FFR as H$_2$. The relation between SFE$_{\rm HI}$ and sSFR shows significant scatter and strong systematic dependence on all of the key galaxy properties that we have explored. The dramatic difference between HI and H$_2$ indicates that different processes (e.g., quenching by different mechanisms) may have very different effects on the HI in different galaxies and hence produce different SFE$_{\rm HI}$-sSFR relations, while the SFE$_{\rm H_2}$-sSFR relation remains unaffected. The facts that SFE$_{\rm H_2}$-sSFR relation is independent of other key galaxy properties, and that sSFR is directly related to the cosmic time and acts as the cosmic clock, make it natural and very simple to study how different galaxy populations (with different properties and undergoing different processes) evolve on the same SFE$_{\rm H_2}$-sSFR $\sim t$ relation. In the gas regulator model (GRM), the evolution of a galaxy on the SFE$_{\rm H_2}$-sSFR($t$) relation is uniquely set by a single mass-loading parameter $\lambda_{\rm net,H_2}$. This simplicity allows us to accurately derive the H$_2$ supply and removal rates of the local galaxy populations with different stellar masses, from star-forming galaxies to the galaxies in the process of being quenched. This combination of FFR and GRM, together with the stellar metallicity requirement, provide a new powerful tool to study galaxy formation and evolution., Comment: 16 pages, 7 figures, accepted for publication in ApJ

The massive galaxy population above the characteristic Schechter mass M*~10^10.6 Msun contributes to about half of the total stellar mass in the local Universe. These massive galaxies usually reside in hot dark matter haloes above the critical shock-heating mass ~10^12 Msun, where the external cold gas supply to these galaxies is expected to be suppressed. When galaxies run out of their cold gas reservoir, they become dead and quiescent. Therefore, massive quiescent galaxies living in hot haloes are commonly believed to be gas-poor. Based on the data from SDSS, ALFALFA, GASS and COLD GASS surveys, here we show that the vast majority of the massive, quiescent, central disk galaxies in the nearby Universe have a remarkably large amount of cold atomic hydrogen gas, surprisingly similar to star-forming galaxies. Both star-forming and quiescent disk galaxies show identical symmetric double-horn HI spectra, indicating similar regularly rotating HI disks. Relative to their star-forming counterparts, massive quiescent central disk galaxies are quenched because of their significantly reduced molecular gas content, lower dust content, and lower star formation efficiency. Our findings reveal a new picture, which clearly demonstrates the detailed star-formation quenching process in massive galaxies and provides a stringent constraint on the physical mechanism of quenching., 10 pages, 4 figures (+2 figures in the appendix); accepted for publication in The Astrophysical Journal Letters

The growth of solids from sub-micron to millimeter and centimeter sizes is the early step toward the formation of planets inside protoplanetary disks (PPDs). However, such processes and their potential impact on the later stages of solid growth are still poorly understood. In this work, we test the hypothesis that most disks contain at least one ringed structure with a relative small radius. We have carried out a large family of 1D two-fluid (gas+dust) hydrodynamical simulations by evolving the gas and dust motion self-consistently while allowing dust size to evolve via coagulation and fragmentation. We investigate the joint effects of ringed structures and dust size growth on the overall sub-millimeter and millimeter (mm) flux and spectral index of PPDs. Ringed structures slow down the dust radial drift and speed up the dust growth. In particular, we find that those unresolved disks with a high fragmentation velocity ($\sim10\ {\rm m~s^{-1}}$) and a high dust surface density ($\sim10\ {\rm g\ cm^{-2}}$ in the ring) can have mm spectral indices as low as $\sim2.0$, consistent with mm observations of faint disks in nearby star forming regions. Furthermore, disks with more than one ringed structure can potentially reproduce brighter disks with spectral indices lower than $\sim2.5$. Future multi-wavelength high-resolution observations of these low spectral index sources can be used to test the existence of the ringed structures in the unresolved disks and differentiate the effects of dust size growth from optical depth., Comment: 18 pages, accepted for publication in ApJ

Based on our previous modeling of AGN feedback in isolated elliptical galaxies (Gan et al. 2014) using the MACER (Massive AGN Controlled Ellipticals Resolved) code, we extend and improve the model to include rotation, to limit star formation to regions of high density and low temperature, to facilitate angular momentum transfer via the Toomre instability in gaseous disks, and to improve the treatment of hot mode (low accretion rate) AGN feedback. The model galaxy now has an extended dark matter profile that matches with standard observations, but has a resolution of parsecs in the inner region and resolves the Bondi radius. We find that the results agree reasonably well with a panoply of observations: (1) both AGN activity and star formation are primarily in central cold gaseous disks, are bursty and mainly driven by the Toomre instability; (2)AGN duty cycle agrees well with the Soltan argument, i.e., the AGN spends most of its lifetime when it is in low luminosity (half of time with $L/L_{Edd}0.06$); (3) the total star formation is $\sim$ few percents of the initial stellar mass, occurring in the bursts that would be associated with the observed E+A phenomenon. Most of the star formation occurs in the circumnuclear disk of a size, 23 pages, 18 figures; accepted by ApJ

We report the analysis of simultaneous XMM-Newton+NuSTAR observations of two low-luminosity Active Galactic Nuclei (LLAGN), NGC 3998 and NGC 4579. We do not detect any significant variability in either source over the ~3 day length of the NuSTAR observations. The broad-band 0.5-60 keV spectrum of NGC 3998 is best fit with a cutoff power-law, while the one for NGC 4579 is best fit with a combination of a hot thermal plasma model, a power-law, and a blend of Gaussians to fit an Fe complex observed between 6 and 7 keV. Our main spectral results are the following: (1) neither source shows any reflection hump with a $3\sigma$ reflection fraction upper-limits $R, Comment: 14 pages, 11 figures, 4 tables, accepted for publication in ApJ

Simultaneous multiwavelength observations were recently performed for three black hole candidates --- SWIFT J1753.5-0127, GRO J1655-40 and XTE J1720-318. In this paper we test the accretion-jet model originally proposed to XTE J1118+480 by investigating the hard state of these three sources using this model. The accretion flow in the model is composed of an inner hot accretion flow and an outer truncated thin disk. We find that the model satisfactorily explains the spectrum ranging from radio to X-rays, with the radio and X-ray spectra dominated by the synchrotron and thermal Comptonization emissions in the jet and the hot accretion flow respectively, while the infrared and optical being the sum of the emissions from the jet, hot accretion flow, and the truncated thin disk. Similar to the case of XTE J1118+480, the model can also explain, although only qualitatively in some cases, the observed timing features including QPO, and positive and negative time lags between the optical and X-ray emissions detected in SWIFT J1753.5-0127. The origin of the ejection events detected in XTE J1720-318 is also briefly discussed., 7 pages, 6 figures accepted for publication in ApJ, minor modification after proof

We investigate the origin of the X-ray emission in low-luminosity AGNs (LLAGNs). Yuan & Cui (2005) predicted that the X-ray emission should originate from jets rather than from an advection-dominated accretion flow (ADAF) when the X-ray luminosity $L_{\rm X}$ of the source is below a critical value of $L_{\rm X,crit} \approx 10^{-6}L_{\rm Edd}$. This prediction implies that the X-ray spectrum in such sources should be fitted by jets rather than ADAFs. Furthermore, below $L_{\rm X,crit}$ the correlation between radio ($L_{\rm R}$) and X-ray ($L_{\rm X}$) luminosities and the black hole mass ($M$)--the so-called fundamental plane of black hole activity--should deviate from the general correlation obtained by Merloni, Heinz & Di Matteo (2003) and become steeper. The Merloni et al. correlation is described by ${\rm log}L_{\rm R} =0.6{\rm log}L_{\rm X}+0.78{\rm log}M+7.33$, while the predicted correlation is ${\rm log}L_{\rm R}=1.23{\rm log}L_{\rm X} +0.25{\rm log}M-13.45$. We collect data from the literature to check the validity of these two expectations. We find that among the 16 LLAGNs with good X-ray and radio spectra, 13 are consistent with the Yuan & Cui prediction. For the 22 LLAGNs with $L_{\rm X} < L_{\rm X,crit}$, the fundamental plane correlation is described by ${\rm log}L_{\rm R}=1.22{\rm log}L_{\rm X}+0.23{\rm log}M-12.46 $, also in excellent agreement with the prediction., 24 pages, 2 tables, 3 figures; accepted by ApJ

We investigate the X-ray origin in FR Is using the radio, submillimetre, optical, and Chandra X-ray data of a small sample consisting of eight FR I sources. These sources are very dim, with X-ray luminosities LX/LEdd ∼ 10 −4 − 10 −8 (LX is the X-ray luminosity between 2-10 keV). We try to fit the multiwaveband spectrum using a coupled accretion-jet model. In this model, the accretion flow is described by an advection-dominated accretion flow (ADAF) while in the innermost region of ADAF a fraction of accretion flow is transferred into the vertical direction and forms a jet. We find that X-ray emission in the source with the highest LX (∼ 1.8×10 −4 LEdd) is from the ADAF. The results for the four sources with moderate LX (∼ several ×10 −6 LEdd) are complicated. Two are mainly from the ADAFs, one from the jet, and the other from the sum of the jet and ADAF. The X-ray emission in the three least luminous sources (LX . 1.0 × 10 −6 LEdd) is dominated by the jet although for one source it can also be interpreted by the ADAF since the quality of X-ray data is low. We conclude that these results roughly support the predictions of Yuan & Cui (2005) where they predict that when the X-ray luminosity of the system is below a critical value, the X-radiation will not be dominated by the emission from the ADAF any longer, but by the jet. We also investigate the fuel supply in these sources. We find that the accretion rates in four sources among the five in which we can have good constraints to their accretion rates must be higher than the Bondi rates. This implies that other fuel supply, such as the gas released by the stellar population inside the Bondi radius, should be important.

Recent radio observations by the VLBA at 7 and 3.5 mm produced the high-resolution images of the compact radio source located at the center of our Galaxy--Sgr A*, and detected its wavelength-dependent intrinsic sizes at the two wavelengths. This provides us with a good chance of testing previously-proposed theoretical models for Sgr A*. In this {\em Letter}, we calculate the size based on the radiatively inefficient accretion flow (RIAF) model proposed by Yuan, Quataert & Narayan (2003). We find that the predicted sizes after taking into account the scattering of the interstellar electrons are consistent with the observations. We further predict an image of Sgr A* at 1.3 mm which can be tested by future observations.

Recently a correlation between the radio and X-ray luminosities is found, $L_{\rm R}\propto L_{\rm X}^{0.7}$, in black hole sources including black hole candidates in our galaxy and active galactic nuclei. We first show that the correlation can be understood in the context of an accretion-jet model developed for explaining the spectral and timing properties of XTE J1118+480. More importantly, we show that when the X-ray luminosity is below a critical value, $\la (10^{-5}$--$10^{-6}) L_{\rm Edd}$, if the jet persists, the correlation should turn and become steeper, $L_{\rm R}\propto L_{\rm X}^{1.23}$, and the X-ray radiation of the system should be dominated by the emission from the jet, rather than by the accretion flow. Possible observational evidence for our predictions is presented and future observations to further test our predictions are proposed., 16 pages, 3 figures, the final version accepted by ApJ

Multi-wavelength observations of the black hole X-ray binary XTE J1118+480 have offered abundant spectral and timing information about the source, and have thus provided serious challenges to theoretical models. We propose a coupled accretion-jet model to interpret the observations. We model the accretion flow as an outer standard thin accretion disk truncated at a transition radius by an inner hot accretion flow. The accretion flow accounts for the observed UV and X-ray emission, but it substantially under-predicts the radio and infrared fluxes, even after we allow for nonthermal electrons in the hot flow. We attribute the latter components to a jet. We model the jet emission by means of the internal shock scenario which is widely employed for gamma-ray bursts. In our accretion-jet model of XTE J1118+480, the jet dominates the radio and infrared emission, the thin disk dominates the UV emission, and the hot flow produces most of the X-ray emission. The optical emission has contributions from all three components: jet, thin disk, and hot flow. The model qualitatively accounts for timing features, such as the intriguing positive and negative time lags between the optical and X-ray emission, and the wavelength-dependent variability amplitude., Comment: 27 pages, 4 figures (one in color); to appear in ApJ in Feb. 2005

Recent infrared (IR) observations of the center of our Galaxy indicate that the supermassive black hole source Sgr A* is strongly variable in the IR. The timescale for the variability, $\sim 30$ min, is comparable to that of the X-ray flares observed by {\em Chandra} and {\em XMM}, suggesting a common physical origin. In this paper, we investigate the nature of the IR emission in the context of models recently proposed to interpret the X-ray flares. We show that the IR emission in Sgr A* can be well explained by nonthermal synchrotron emission if a small fraction of the electrons in the innermost region of the accretion flow around the black hole are accelerated {into a broken power-law distribution}, perhaps due to transient events such as magnetic reconnection. The model predicts differences in the variability amplitudes of flares in the IR and X-rays, in general agreement with observations. It also predicts that the IR emission should be linearly polarized, as has indeed been observed during one epoch. IR and X-ray flares analogous to those observed in Sgr A* may be detectable from other accreting SMBHs, provided $L \lsim 10^{-8} L_{EDD}$; at higher luminosities the flaring emission is dominated by thermal synchrotron-self Compton emission, which is likely to be less variable.

We investigate radiatively inefficient accretion flow models for Sgr A*, the supermassive black hole in our Galactic Center, in light of new observational constraints. Confirmation of linear polarization in the submm emission argues for accretion rates much less than the canonical Bondi rate. We consider models with low accretion rates and calculate the spectra produced by a hybrid electron population, consisting of both thermal and nonthermal particles. The thermal electrons produce the submm emission and can account for its linear polarization properties. As noted in previous work, the observed low-frequency radio spectrum can be explained if a small fraction ($\approx 1.5$%) of the electron thermal energy resides in a soft power-law tail. In the innermost region of the accretion flow, turbulence and/or magnetic reconnection events may occasionally accelerate a fraction of the electrons into a harder power-law tail. We show that the synchrotron emission from these electrons, or the Compton up-scattering of synchrotron photons by the same electrons, may account for the X-ray flares observed by {\it Chandra}., Comment: 32 pages, 7 figures, the final version accepted by ApJ

Global evolution and dispersal of protoplanetary disks (PPDs) is governed by disk angular momentum transport and mass-loss processes. Recent numerical studies suggest that angular momentum transport in the inner region of PPDs is largely driven by magnetized disk wind, yet the wind mass-loss rate remains unconstrained. On the other hand, disk mass loss has conventionally been attributed to photoevaporation, where external heating on the disk surface drives a thermal wind. We unify the two scenarios by developing a 1D model of magnetized disk winds with a simple treatment of thermodynamics as a proxy for external heating. The wind properties largely depend on 1) the magnetic field strength at the wind base, characterized by the poloidal Alfv\'en speed $v_{Ap}$, 2) the sound speed $c_s$ near the wind base, and 3) how rapidly poloidal field lines diverge (achieve $R^{-2}$ scaling). When $v_{Ap}\gg c_s$, corotation is enforced near the wind base, resulting in centrifugal acceleration. Otherwise, the wind is accelerated mainly by the pressure of the toroidal magnetic field. In both cases, the dominant role played by magnetic forces likely yields wind outflow rates that well exceed purely hydrodynamical mechanisms. For typical PPD accretion-rate and wind-launching conditions, we expect $v_{Ap}$ to be comparable to $c_s$ at the wind base. The resulting wind is heavily loaded, with total wind mass loss rate likely reaching a considerable fraction of wind-driven accretion rate. Implications for modeling global disk evolution and planet formation are also discussed., Comment: 21 pages, 12 figures, submitted to ApJ

Previous works show strong winds exist in hot accretion flows around black holes. Those works focus only on the region close to the black hole thus it is unknown whether or where the wind production stops at large radii. In this paper, we investigate this problem by hydrodynamical simulations. We take into account gravities of both the black hole and the nuclear star clusters. For the latter, we assume that the velocity dispersion of stars is a constant and its gravitational potential $\propto \sigma^2 \ln (r)$, where $\sigma$ is the velocity dispersion of stars and $r$ is the distance from the center of the galaxy. We focus on the region where the gravitational potential is dominated by the star cluster. We find, same as the accretion flow at small radii, the mass inflow rate decreases inward and the flow is convectively unstable. However, trajectory analysis shows that there is very few wind launched from the flow. Our result, combined with the results of Yuan et al. (2015), indicates that the mass flux of wind launched from hot accretion flow $\dot{M}_{\rm wind}=\dot{M}_{\rm BH}(r/20r_s)$, with $r\la R_A\equiv GM_{\rm BH}/\sigma^2$. Here $\dot{M}_{\rm BH}$ is accretion rate at black hole horizon. $R_A$ is similar to Bondi radius. We argue that the inward decrease of inflow rate is not because of mass loss via wind, but because of convective motion. The disappearance of wind outside of $R_A$ must be because of the change of the gravitational potential, but the exact reason remains to be probed., Comment: 9 pages, 11 figures, accepted by ApJ

In a previous work, we have shown that the formation of Fermi bubbles can be due to the interaction between winds launched from the hot accretion flow in Sgr A* and the interstellar medium (ISM). In that work, we focus only on the morphology. In this paper we continue our study by calculating the gamma-ray radiation. Some cosmic-ray protons (CRp) and electrons (CRe) must be contained in the winds, which are likely formed by physical processes such as magnetic reconnection. We have performed MHD simulations to study the spatial distribution of CRp, considering the advection and diffusion of CRp in the presence of magnetic field. We find that a permeated zone is formed just outside of the contact discontinuity between winds and the ISM, where the collisions between CRp and thermal nuclei mainly occur. The decay of neutral pions generated in the collisions, combined with the inverse Compton scattering of background soft photons by the secondary leptons generated in the collisions and primary CRe, can well explain the observed gamma-ray spectral energy distribution. Other features such as the uniform surface brightness along the latitude and the boundary width of the bubbles are also explained. The advantage of this "accretion wind" model is that the adopted wind properties come from the detailed small-scale MHD numerical simulation of accretion flows and the value of mass accretion rate has independent observational evidences. The success of the model suggests that we may seriously consider the possibility that cavities and bubbles observed in other contexts such as galaxy clusters may be formed by winds rather than jets.

The Local Group compact elliptical galaxy M32 hosts one of the nearest candidate supermassive black holes (SMBHs), which has a previously suggested X-ray counterpart. Based on sensitive observations taken with the Karl G. Jansky Very Large Array (VLA), we detect for the first time a compact radio source coincident with the nucleus of M32, which exhibits an integrated flux density of similar to 47.3. +/- 6.1 mu Jy at 6.6 GHz. We discuss several possibilities for the nature of this source, favoring an origin of the long-sought radio emission from the central SMBH, for which we also revisit the X-ray properties based on recently acquired Chandra and XMM-Newton data. Our VLA observations also discover radio emission from three previously known optical planetary nebulae in the inner region of M32.

We study the dynamics of super-Eddington accretion flows by performing two-dimensional radiation-hydrodynamic simulations. Compared with previous works, in this paper we include the $T_{\theta\phi}$ component of the viscous stress and consider various values of the viscous parameter $\alpha$. We find that when $T_{\theta\phi}$ is included, the rotational speed of the high-latitude flow decreases, while the density increases and decreases at the high and low latitudes, respectively. We calculate the radial profiles of inflow and outflow rates. We find that the inflow rate decreases inward, following a power law form of $\dot{M}_{\rm in}\propto r^s$. The value of $s$ depends on the magnitude of $\alpha$ and is within the range of $\sim 0.4-1.0$. Correspondingly, the radial profile of density becomes flatter compared with the case of a constant $\dot{M}(r)$. We find that the density profile can be described by $\rho(r)\propto r^{-p}$, and the value of $p$ is almost same for a wide range of $\alpha$ ranging from $\alpha=0.1$ to $0.005$. The inward decrease of inflow accretion rate is very similar to hot accretion flows, which is attributed to the mass loss in outflows. To study the origin of outflow, we analyze the convective stability of slim disk. We find that depending on the value of $\alpha$, the flow is marginally stable (when $\alpha$ is small) or unstable (when $\alpha$ is large). This is different from the case of hydrodynamical hot accretion flow where radiation is dynamically unimportant and the flow is always convectively unstable. We speculate that the reason for the difference is because radiation can stabilize convection. The origin of outflow is thus likely because of the joint function of convection and radiation, but further investigation is required., Comment: 16 pages, 13 figures, accepted for publication in ApJ

Previous hydrodynamical (HD) and magnetohydrodynamical (MHD) numerical simulations of hot accretion flows have shown that the mass accretion rate decreases with decreasing radius. Two models have been proposed to explain this result. In the ADIOS model, the inward decrease of accretion rate is because of the loss of gas in the outflow. In the CDAF model, the gas is assumed to be locked in convective eddies, which results in the inward decrease of the accretion rate. We investigate the nature of inward decrease of accretion rate using HD and MHD simulations. We calculate various properties of inflow and outflow, including the mass flux, radial and rotational velocities, temperature, and the Bernoulli parameter ($Be$). Systematic and significant differences between inflow and outflow are found. These results suggest that the inflow and outflow are not dominated by convective turbulence, but are systematic inward and outward motion. We have also analyzed the convective stability of MHD accretion flow and found that they are convectively stable. These results indicate that the ADIOS scenario is favored. The different properties of inflow and outflow also suggest that the mechanisms of producing outflow in HD and MHD flows are buoyancy and centrifugal force associated with the magnetic field, respectively. The latter mechanism is similar to the Blandford & Payne mechanism. We also study the effect of initial conditions in the simulations. We find that the value of $Be$ is mainly determined by the value of $Be$ of the initial condition. We discuss some possible observational applications of our outflow model. These observations include the Fermi bubble observed in the Galaxy center, and winds widely observed in AGNs and black hole X-ray binaries., 24 emulateapj pages, 15 figures; one figure and some references added, conclusion unchanged; ApJ in press

In a hot accretion flow, the radiation from the innermost region of the flow propagates outward and heats the electrons at large radii via Compton scattering. It has been shown in previous works that if the radiation is strong enough, L greater than or similar to 2% L-Edd, the electrons at the Bondi radius (r(B) similar to 10(5)r(s)) will be heated to above the virial temperature; thus, the accretion will be stopped. The accretion will recover after the gas cools down. This results in the oscillation of the black hole activity. In this paper, we show that this mechanism is the origin of the intermittent activity of some compact young radio sources. Such intermittency is required to explain the population of these sources. We calculate the timescales of the black hole oscillation and find that the durations of active and inactive phases are 3 x 10(4)(0.1/alpha)(M/10(8) M-circle dot)(L/2% L-Edd)(-1/2) yr and 10(5)(alpha/0.1)(M/10(8) M-circle dot) yr, respectively, consistent with those required to explain observations. Such feedback occurring at the parsec scale should be common in low-luminosity active galactic nuclei and should be considered when we consider their matter and energy output.

Low-luminosity active galactic nuclei (LLAGNs) are generally believed to be powered by an inner radiatively inefficient, advection-dominated accretion flow (ADAF), an outer truncated thin disk, and a jet. Maoz recently challenged this picture based on the observation that the strength of ultraviolet emission relative to the X-ray and radio bands does not depart from empirical trends defined by more luminous sources. He advocates that AGNs across all luminosities have essentially the same accretion and radiative processes, which in luminous sources are described by a standard optically thick, geometrically thin disk. We calculate ADAF models and demonstrate that they can successfully fit the observed spectral energy distributions of the LLAGNs in Maoz's sample. Our model naturally accommodates the radio and X-ray emission, and the ultraviolet flux is well explained by a combination of the first-order Compton scattering in the ADAF, synchrotron emission in the jet, and black body emission in the truncated thin disk. It is premature to dismiss the ADAF model for LLAGNs. The UV data can be fitted equally well using a standard thin disk, but an additional corona and jet would be required to account for the X-ray and radio emission. We argue that there are strong theoretical reasons to prefer the ADAF model over the thin disk scenario. We discuss testable predictions that can potentially discriminate between the two accretion models.