Source author record

Ye-Fei Yuan

Ye-Fei Yuan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

31works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

31 published item(s)

preprint2020arXiv

Broad-line Region of the Quasar PG 2130+099 from a Two-Year Reverberation Mapping Campaign with High Cadence

As one of the most interesting Seyfert 1 galaxies, PG 2130+099 has been the target of several reverberation mapping (RM) campaigns over the years. However, its measured broad H$β$ line responses have been inconsistent, with time lags of $\sim$200 days, $\sim$25 days, and $\sim$10 days being reported for different epochs while its optical luminosity changed no more than 40%. To investigate this issue, we conducted a new RM-campaign with homogenous and high cadence (about $\sim$3 days) for two years during 2017--2019 to measure the kinematics and structure of the ionized gas. We successfully detected time lags of broad H$β$, He II, He I, and Fe II lines with respect to the varying 5100Åcontinuum, revealing a stratified structure that is likely virialized with Keplerian kinematics in the first year of observations, but an inflow kinematics of the broad-line region from the second year. With a central black hole mass of $0.97_{-0.18}^{+0.15}\times 10^7~M_{\odot}$, PG 2130+099 has an accretion rate of $10^{2.1\pm0.5}L_{\rm Edd}c^{-2}$, where $L_{\rm Edd}$ is the Eddington luminosity and $c$ is speed of light, implying that it is a super-Eddington accretor and likely possesses a slim, rather than thin, accretion disk. The fast changes of the ionization structures of the three broad lines remain puzzling.

preprint2020arXiv

On the Cooling Flow Problem in the Gaseous Halo of the Milky Way

Theoretical and observational arguments suggest that there is a large amount of hot ($\sim 10^6$ K), diffuse gas residing in the Milky Way's halo, while its total mass and spatial distribution are still unclear. In this work, we present a general model for the gas density distribution in the Galactic halo, and investigate the gas evolution under radiative cooling with a series of 2D hydrodynamic simulations. We find that the mass inflow rate in the developed cooling flow increases with gas metallicity and the total gas mass in the halo. For a fixed halo gas mass, the spatial gas distribution affects the onset time of the cooling catastrophe, which starts earlier when the gas distribution is more centrally-peaked, but does not substantially affect the final mass inflow rate. The gravity from the Galactic bulge and disk affects gas properties in inner regions, but has little effect on the final inflow rate either. We confirm our results by investigating cooling flows in several density models adopted from the literature, including the Navarro-Frenk-White (NFW) model, the cored-NFW model, the Maller & Bullock model, and the $β$ model. Typical mass inflow rates in our simulations range from $\sim 5 M_{\odot}$ yr$^{-1}$ to $\sim 60 M_{\odot}$ yr$^{-1}$, and are much higher than the observed star formation rate in our Galaxy, suggesting that stellar and active galactic nucleus feedback processes may play important roles in the evolution of the Milky Way (MW) and MW-type galaxies.

preprint2020arXiv

Revisiting the Relationship between the Long GRB Rate and Cosmic Star Formation History Based on a Large Swift Sample

The exact relationship between the long gamma-ray burst (LGRB) rate and the cosmic star formation rate (CSFR) is essential for using LGRBs as cosmological probes. In this work, we collect a large sample composed of 371 Swift LGRBs with known redshifts and prompt emission properties. We first compare the rest-frame prompt properties of these bursts in different redshift bins, finding negligible redshift evolution of the luminosity of LGRBs with $L_{\mathrm{iso}}\gtrsim10^{51}\,\mathrm{erg\, s^{-1}}$ between $z\sim1$ and $z\sim4$. Then, by utilizing the CSFR obtained from the large-scale cosmological hydrodynamical simulation, the Illustris simulation, we calculate the cumulative redshift distribution of LGRBs under different metallicity thresholds. After comparing with our sample, we find that the predictions with a moderate threshold between $0.3\,Z_{\odot}\leqslant Z_{\mathrm{th}}\leqslant1.0\,Z_{\odot}$ are consistent with the sample between redshift $0<z<3$, while at higher redshifts between $3<z<5$, all metallicity thresholds fit the data well. When changing to an empirical model based on observations, the predictions show similar results as well. After comparing with the metallicity distribution of the observed LGRB host galaxies between $0<z<1$, we confirm that the production of LGRBs in galaxies with super-solar metallicity is suppressed. Nevertheless, considering that a significant fraction of stars are born in sub-solar metallicity environments at $z\gtrsim3$, we suggest that, as a first approximation, LGRBs can be used as direct tracers of the CSFR in this redshift range.

preprint2020arXiv

Supermassive black holes with high accretion rates in active galactic nuclei. XI. Accretion disk reverberation mapping of Mrk 142

We performed an intensive accretion disk reverberation mapping campaign on the high accretion rate active galactic nucleus Mrk 142 in early 2019. Mrk 142 was monitored with the Neil Gehrels Swift Observatory for 4 months in X-rays and 6 UV/optical filters. Ground-based photometric monitoring was obtained from the Las Cumbres Observatory, Liverpool Telescope and Dan Zowada Memorial Observatory in ugriz filters and the Yunnan Astronomical Observatory in V. Mrk 142 was highly variable throughout, displaying correlated variability across all wavelengths. We measure significant time lags between the different wavelength light curves, finding that through the UV and optical the wavelength-dependent lags, $τ(λ)$, generally follow the relation $τ(λ) \propto λ^{4/3}$, as expected for the $T\propto R^{-3/4}$ profile of a steady-state optically-thick, geometrically-thin accretion disk, though can also be fit by $τ(λ) \propto λ^{2}$, as expected for a slim disk. The exceptions are the u and U band, where an excess lag is observed, as has been observed in other AGN and attributed to continuum emission arising in the broad-line region. Furthermore, we perform a flux-flux analysis to separate the constant and variable components of the spectral energy distribution, finding that the flux-dependence of the variable component is consistent with the $f_ν\proptoν^{1/3}$ spectrum expected for a geometrically-thin accretion disk. Moreover, the X-ray to UV lag is significantly offset from an extrapolation of the UV/optical trend, with the X-rays showing a poorer correlation with the UV than the UV does with the optical. The magnitude of the UV/optical lags is consistent with a highly super-Eddington accretion rate.

preprint2020arXiv

Testing the equivalence principle via the shadow of black holes

We study the equivalence principle, regarded as the cornerstone of general relativity, by analyzing the deformation observable of black hole shadows. Such deformation can arise from new physics and may be expressed as a phenomenological violation of the equivalence principle. Specifically, we assume that there is an additional background vector field that couples to the photons around the supermassive black hole. This type of coupling yields impact on the way the system depends on initial conditions, and affects the black hole shadow at different wavelengths by a different amount, and therefore observations of the shadow in different wavelengths could constrain such couplings. This can be tested by future multi-band observations. Adopting a specific form of the vector field, we obtain constraints on model parameters from Event Horizon Telescope observations and measurements of gas/stellar orbits.

preprint2019arXiv

Active Galactic Nuclei with Ultra-fast Outflows Monitoring Project: The Broad-line Region of Mrk 79 as a Disk Wind

We developed a spectroscopic monitoring project to investigate the kinematics of the broad-line region (BLR) in active galactic nuclei (AGN) with ultra-fast outflows (UFOs). Mrk~79 is a radio-quiet AGN with UFOs and warm absorbers, had been monitored by three reverberation mapping (RM) campaigns, but its BLR kinematics is not understood yet. In this paper, we report the results from a new RM-campaign of Mrk~79, which was undertaken by Lijiang 2.4-m telescope. Mrk~79 is seeming to come out the faint state, the mean flux approximates a magnitude fainter than historical record. We successfully measured the lags of the broad emission lines including H$β~\lambda4861$, H$γ~\lambda4340$, He II $\lambda4686$ and He I $\lambda5876$ with respect to the varying AGN continuum. Based on the broad H$β~\lambda4861$ line, we measured black hole (BH) mass of $M_{\bullet}=5.13^{+1.57}_{-1.55}\times10^{7}M_{\odot}$, estimated accretion rates of ${\dot{M}_{\bullet}}=(0.05\pm0.02)~L_{\rm Edd}~c^{-2}$, indicating that Mrk~79 is a sub-Eddington accretor. We found that Mrk~79 deviates from the canonical Radius$-$Luminosity relationship. The marginal blueshift of the broad He II $\lambda4686$ line detected from rms spectrum indicates outflow of high-ionization gas. The velocity-resolved lag profiles of the broad H$γ~\lambda4340$, H$β~\lambda4861$, and He I $\lambda5876$ lines show similar signatures that the largest lag occurs in the red wing of the lines then the lag decreases to both sides. These signatures should suggest that the BLR of Keplerian motion probably exists the outflow gas motion. All findings including UFOs, warm absorbers, and the kinematics of high- and low-ionization BLR, may provide an indirect evidence that the BLR of Mrk~79 probably originates from disk wind.

preprint2019arXiv

Two-dimensional Inflow-outflow Solution of Supercritical Accretion Flow

We present the two-dimensional inflow-outflow solutions of radiation hydrodynamic (RHD) equations of supercritical accretion flows. Compared with prior studies, we include all components of the viscous stress tensor. We assume steady state flow and use self-similar solutions in the radial direction to solve the equations in $ r-θ$ domain of the spherical coordinates. The set of differential equations have been integrated from the rotation axis to the equatorial plane. We find that the self-similarity assumption requires that the radial profile of density is described by $ ρ(r) \propto r^{-0.5} $. Correspondingly, the radial profile of the mass inflow rate decreases with decreasing radii as $ \dot{M}_\mathrm{in} \propto r $. Inflow-outflow structure has been found in our solution. In the region $ θ> 65^{\circ} $ there exist inflow while above that flow moves outward and outflow could launch. The driving forces of the outflow are analyzed and found that the radiation force is dominant and push the gas particles outwards with poloidal velocity $ \sim 0.25 c $. The properties of outflow are also studied. The results show that the mass flux weighted angular momentum of the inflow is lower than that of outflow, thus the angular momentum of the flow can be transported by the outflow. We also analyze the convective stability of the supercritical disk and find that in the absence of the magnetic field, the flow is convectively unstable. Our analytical results are fully consistent with the previous numerical simulations of the supercritical accretion flow.

preprint2016arXiv

Effects of anisotropic thermal conduction on wind properties in hot accretion flow

Previous works have clearly shown the existence of winds from black hole hot accretion flow and investigated their detailed properties. In extremely low accretion rate systems, the collisional mean-free path of electrons is large compared with the length-scale of the system, thus thermal conduction is dynamically important. When the magnetic field is present, the thermal conduction is anisotropic and energy transport is along magnetic field lines. In this paper, we study the effects of anisotropic thermal conduction on the wind production in hot accretion flows by performing two-dimensional magnetohydrodynamic simulations. We find that thermal conduction has only moderate effects on the mass flux of wind. But the energy flux of wind can be increased by a factor of $\sim 10$ due to the increase of wind velocity when thermal conduction is included. The increase of wind velocity is because of the increase of driving forces (e.g. gas pressure gradient force and centrifugal force) when thermal conduction is included. This result demonstrates that thermal conduction plays an important role in determining the properties of wind.

preprint2016arXiv

Hot accretion flow with radiative cooling: state transitions in black hole X-ray binaries

We investigate state transitions in black hole X-ray binaries through different parameters by using two-dimensional axisymmetric hydrodynamical simulation method. For radiative cooling in hot accretion flow, we take into account the bremsstrahlung, synchrotron and synchrotron-self Comptonization self-consistently in the dynamics. Our main result is that the state transitions occur when the accretion rate reaches a critical value $\dot M \sim 3α \dot M_{\rm Edd}$, above which cold and dense clumpy/filamentary structures are formed, embedded within the hot gas. We argued this mode likely corresponds to the proposed two-phase accretion model, which may be responsible for the intermediate state of black hole X-ray binaries. When the accretion rate becomes sufficiently high, the clumpy/filamentary structures gradually merge and settle down onto the mid-plane. Eventually the accretion geometry transforms to a disc-corona configuration. In summary our results are consistent with the truncated accretion scenario for the state transition.

preprint2016arXiv

Supermassive Black Holes with High Accretion Rates in Active Galactic Nuclei. V. A New Size-Luminosity Scaling Relation for the Broad-Line Region

This paper reports results of the third-year campaign of monitoring super-Eddington accreting massive black holes (SEAMBHs) in active galactic nuclei (AGNs) between 2014-2015. Ten new targets were selected from quasar sample of Sloan Digital Sky Survey (SDSS), which are generally more luminous than the SEAMBH candidates in last two years. H$β$ lags ($τ_{_{\rm Hβ}}$) in five of the 10 quasars have been successfully measured in this monitoring season. We find that the lags are generally shorter, by large factors, than those of objects with same optical luminosity, in light of the well-known $R_{_{\rm Hβ}}-L_{5100}$ relation. The five quasars have dimensionless accretion rates of $\dot{\mathscr{M}}=10-10^3$. Combining measurements of the previous SEAMBHs, we find that the reduction of H$β$ lags tightly depends on accretion rates, $τ_{_{\rm Hβ}}/τ_{_{R-L}}\propto\dot{\mathscr{M}}^{-0.42}$, where $τ_{_{R-L}}$ is the H$β$ lag from the normal $R_{_{\rm Hβ}}-L_{5100}$ relation. Fitting 63 mapped AGNs, we present a new scaling relation for the broad-line region: $R_{_{\rm Hβ}}=α_1\ell_{44}^{β_1}\,\min\left[1,\left(\dot{\mathscr{M}}/\dot{\mathscr{M}}_c\right)^{-γ_1}\right]$, where $\ell_{44}=L_{5100}/10^{44}\,\rm erg~s^{-1}$ is 5100 Å continuum luminosity, and coefficients of $α_1=(29.6_{-2.8}^{+2.7})$ lt-d, $β_1=0.56_{-0.03}^{+0.03}$, $γ_1=0.52_{-0.16}^{+0.33}$ and $\dot{\mathscr{M}}_c=11.19_{-6.22}^{+2.29}$. This relation is applicable to AGNs over a wide range of accretion rates, from $10^{-3}$ to $10^3$. Implications of this new relation are briefly discussed.

preprint2016arXiv

Supermassive Black Holes with High Accretion Rates in Active Galactic Nuclei. VI. Velocity-resolved Reverberation Mapping of H$β$ Line

In the sixth of the series of papers reporting on a large reverberation mapping (RM) campaign of active galactic nuclei (AGNs) with high accretion rates, we present velocity-resolved time lags of H$β$ emission lines for nine objects observed in the campaign during 2012$-$2013. In order to correct the line-broadening caused by seeing and instruments before the analysis of velocity-resolved RM, we adopt Richardson-Lucy deconvolution to reconstruct their H$β$ profiles. The validity and effectiveness of the deconvolution are checked out by Monte Carlo simulation. Five among the nine objects show clear dependence of time delay on velocity. Mrk 335 and Mrk 486 show signatures of gas inflow whereas the clouds in the broad-line regions (BLRs) of Mrk 142 and MCG +06-26-012 tend to be radial outflowing. Mrk 1044 is consistent with the case of virialized motions. The lags of the rest four are not velocity-resolvable. The velocity-resolved RM of super-Eddington accreting massive black holes (SEAMBHs) shows that they have diversity of the kinematics in their BLRs. Comparing with the AGNs with sub-Eddington accretion rates, we do not find significant differences in the BLR kinematics of SEAMBHs.

preprint2015arXiv

Constraints on black hole spins with a general relativistic accretion disk corona model

The peaks of the spectra of the accretion disks surrounding massive black holes in quasars are in the far-UV or soft X-ray band, which are usually not observed. However, in the disk corona model, the soft photons from the disk are Comptonized to high energy in the hot corona, and the hard X-ray spectra (luminosity and spectral shape) contain the information of the incident spectra from the disk. The values of black hole spin parameter $a_{\ast}$ are inferred from the spectral fitting, which spread over a large range, $\sim -0.94$ to $0.998$. We find that the inclination angles and mass accretion rates are well determined by the spectral fitting, while the results are sensitive to the accuracy of black hole mass estimates. No tight constraints on the black hole spins are achieved, if the uncertainties of black hole mass measurements are a factor of four, which are typical for the single-epoch reverberation mapping method. Recently, the accuracy of black hole mass measurement has been significantly improved to $0.2-0.4$~dex with velocity resolved reverberation mapping method (Pancoast et al., 2014). The black hole spin can be well constrained if the mass measurement accuracy is $\la 50$\%. In the accretion disk corona scenario, a fraction of power dissipated in the disk is transported into the corona, and therefore the accretion disk is thinner than a bare disk for the same mass accretion rate, because the radiation pressure in the disk is reduced. We find that the thin disk approximation, $H/R\la 0.1$, is still valid if $0.3<\dot{m}<0.5$, provided a half of the dissipated power is radiated in the corona above the disk.

preprint2015arXiv

Cosmic reionization of hydrogen and helium: contribution from both mini-quasars and stars

Observations on the high-redshift galaxies at $z>6$ imply that their ionizing emissivity is unable to fully reionize the Universe at $z\sim 6$. Either a high escape fraction of ionizing photons from these galaxies or a large population of faint galaxies below the detection limit are required. However, these requirements are somewhat in tension with present observations. In this work, we explored the combined contribution of mini-quasars and stars to the reionization of cosmic hydrogen and helium. Our model is roughly consistent with: (1) the low escape fractions of ionizing photons from the observed galaxies, (2) the optical depth of Cosmic Microwave Background (CMB) measured by the WMAP-7, and (3) the redshift of the end of hydrogen and helium reionization at $z\approx 6$ and $z\approx 3$, respectively. Neither an extremely high escape fraction nor a large population of fainter galaxies is required in this scenario. In our most optimistic model, more than $\sim20\%$ of the cosmic helium is reionized by $z\sim6$, and the ionized fraction of cosmic helium rapidly climbs to more than $50\%$ by $z\sim5$. These results may imply that better measurements of helium reionization, especially at high redshifts, could be helpful in constraining the growth of intermediate-mass black holes (IMBHs) in the early Universe, which would shed some light on the puzzles concerning the formation of supermassive black holes (SMBHs).

preprint2015arXiv

Long GRBs as a Tool to Investigate Star Formation in Dark Matter Halos

First stars can only form in structures that are suitably dense, which can be parametrized by the minimum dark matter halo mass $M_{\rm min}$. $M_{\rm min}$ must plays an important role in star formation. The connection of long gamma-ray bursts (LGRBs) with the collapse of massive stars has provided a good opportunity for probing star formation in dark matter halos. We place some constraints on $M_{\rm min}$ using the latest $Swift$ LGRB data. We conservatively consider that LGRB rate is proportional to the cosmic star formation rate (CSFR) and an additional evolution parametrized as $(1+z)^α$, where the CSFR model as a function of $M_{\rm min}$. Using the $χ^{2}$ statistic, the contour constraints on the $M_{\rm min}$--$α$ plane show that at the $1σ$ confidence level, we have $M_{\rm min}<10^{10.5}$ $\rm M_{\odot}$ from 118 LGRBs with redshift $z<4$ and luminosity $L_{\rm iso}>1.8\times10^{51}$ erg $\rm s^{-1}$. We also find that adding 12 high-\emph{z} $(4<z<5)$ LGRBs (consisting of 104 LGRBs with $z<5$ and $L_{\rm iso}>3.1\times10^{51}$ erg $\rm s^{-1}$) could result in much tighter constraints on $M_{\rm min}$, for which, $10^{7.7}\rm M_{\odot}<M_{\rm min}<10^{11.6}\rm M_{\odot}$ ($1σ$). Through Monte Carlo simulations, we estimate that future five years of Sino-French spacebased multiband astronomical variable objects monitor (\emph{SVOM}) observations would tighten these constraints to $10^{9.7}\rm M_{\odot}<M_{\rm min}<10^{11.3}\rm M_{\odot}$. The strong constraints on $M_{\rm min}$ indicate that LGRBs are a new promising tool for investigating star formation in dark matter halos.

preprint2015arXiv

Supermassive Black Holes with High Accretion Rates in Active Galactic Nuclei. IV. H$β$ Time Lags and Implications for Super-Eddington Accretion

We have completed two years of photometric and spectroscopic monitoring of a large number of active galactic nuclei (AGNs) with very high accretion rates. In this paper, we report on the result of the second phase of the campaign, during 2013--2014, and the measurements of five new H$β$ time lags out of eight monitored AGNs. All five objects were identified as super-Eddington accreting massive black holes (SEAMBHs). The highest measured accretion rates for the objects in this campaign are $\dot{\mathscr{M}}\gtrsim 200$, where $\dot{\mathscr{M}}= \dot{M}_{\bullet}/L_{\rm Edd}c^{-2}$, $\dot{M}_{\bullet}$ is the mass accretion rates, $L_{\rm Edd}$ is the Eddington luminosity and $c$ is the speed of light. We find that the H$β$ time lags in SEAMBHs are significantly shorter than those measured in sub-Eddington AGNs, and the deviations increase with increasing accretion rates. Thus, the relationship between broad-line region size ($R_{_{\rm Hβ}}$) and optical luminosity at 5100Å, $R_{_{\rm Hβ}}-L_{5100}$, requires accretion rate as an additional parameter. We propose that much of the effect may be due to the strong anisotropy of the emitted slim-disk radiation. Scaling $R_{_{\rm Hβ}}$ by the gravitational radius of the black hole, we define a new radius-mass parameter ($Y$) and show that it saturates at a critical accretion rate of $\dot{\mathscr{M}}_c=6\sim 30$, indicating a transition from thin to slim accretion disk and a saturated luminosity of the slim disks. The parameter $Y$ is a very useful probe for understanding the various types of accretion onto massive black holes. We briefly comment on implications to the general population of super-Eddington AGNs in the universe and applications to cosmology.

preprint2014arXiv

Episodic jet power extracted from a spinning black hole surrounded by a neutrino-dominated accretion flow in gamma-ray bursts

It was suggested that the relativistic jets in gamma-ray bursts (GRBs) are powered via the Blandford-Znajek (BZ) mechanism or the annihilation of neutrinos and anti-neutrinos from a neutrino cooling-dominated accretion flow (NDAF). The advection and diffusion of the large-scale magnetic field of a NDAF is calculated, and the external magnetic field is found to be dragged inward efficiently by the accretion flow for a typical magnetic Prandtl number P_m=1. The maximal BZ jet power can be ~10^53-10^54 erg/sec for an extreme Kerr black hole, if an external magnetic field with 10^14 Gauss is advected by the NDAF. This is roughly consistent with the field strength of the disk formed after a tidal disrupted magnetar. The accretion flow near the black hole horizon is arrested by the magnetic field if the accretion rate is below than a critical value for a given external field. The arrested accretion flow fails to drag the field inward and the field strength decays, and then the accretion re-starts, which leads to oscillating accretion. The typical timescale of such episodic accretion is in an order of one second. This can qualitatively explain the observed oscillation in the soft extend emission of short-type GRBs.

preprint2013arXiv

Is the metallicity of the progenitor of long gamma-ray bursts really low?

Observations of long gamma-ray bursts (LGRBs) offer a unique opportunity for probing the cosmic star formation history, although whether or not LGRB rates are biased tracers of star formation rate history is highly debated. Based on an extensive sample of LGRBs compiled by Robertson & Ellis (2012), we analyze various models of star formation rate and the possible effect of the evolution of cosmic metallicity under the assumption that LGRBs tend to occur in low-metallicity galaxies. The models of star formation rate tested in this work include empirical fits from observational data as well as a self-consistent model calculated in the framework of the hierarchical structure formation. Comparing with the observational data, we find a relatively higher metallicity cut of $Z\gtrsim0.6Z_{\odot}$ for the empirical fits and no metallicity cut for the self-consistent model. These results imply that there is no strong metallicity preference for the host galaxy of LGRBs, in contrast to previous work which suggest a cut of $Z\sim0.1-0.3Z_{\odot}$, and that the inferred dependencies of LGRBs on their host galaxy properties are strongly related to the specific models of star formation rate. Furthermore, a significant fraction of LGRBs occur in small dark matter halos down to $3\times10^{8}\,\mathrm{M_{\odot}}$ can provide an alternative explanation for the discrepancy between the star formation rate history and LGRB rate history.

preprint2013arXiv

Mass flow in the circumbinary disk with gap around supermassive binary black holes

In this paper, we study the interaction between the supermassive binary black holes in elliptical orbit and their surrounding disk with a gap. The gap in the disk is a low density region formed due to the tidal effects of the less massive black hole. The binary we have investigated has a sub-parsec separation and is coplanar with the disk. We find that the maximum variation of the surface density in the gap reaches 50% during an orbital period. However, in other regions of the disk, the density variation is much less than 1%. Furthermore, we calculate the corresponding variation of spectral energy distribution within a period, but little variation is found. The reason for these results is that the viscosity timescale of the disk at the binary radius is much longer than the orbital period of the binary.

preprint2013arXiv

Progenitor delay-time distribution of short gamma-ray bursts: Constraints from observations

Context. The progenitors of short gamma-ray bursts (SGRBs) have not yet been well identified. The most popular model is the merger of compact object binaries (NS-NS/NS-BH). However, other progenitor models cannot be ruled out. The delay-time distribution of SGRB progenitors, which is an important property to constrain progenitor models, is still poorly understood. Aims. We aim to better constrain the luminosity function of SGRBs and the delay-time distribution of their progenitors with newly discovered SGRBs. Methods. We present a low-contamination sample of 16 Swift SGRBs that is better defined by a duration shorter than 0.8 s. By using this robust sample and by combining a self-consistent star formation model with various models for the distribution of time delays, the redshift distribution of SGRBs is calculated and then compared to the observational data. Results. We find that the power-law delay distribution model is disfavored and that only the lognormal delay distribution model with the typical delay tau >= 3 Gyr is consistent with the data. Comparing Swift SGRBs with T90 > 0.8 s to our robust sample (T90 < 0.8 s), we find a significant difference in the time delays between these two samples. Conclusions. Our results show that the progenitors of SGRBs are dominated by relatively long-lived systems (tau >= 3 Gyr), which contrasts the results found for Type Ia supernovae. We therefore conclude that primordial NS-NS systems are not favored as the dominant SGRB progenitors. Alternatively, dynamically formed NS-NS/BH and primordial NS-BH systems with average delays longer than 5 Gyr may contribute a significant fraction to the overall SGRB progenitors.

preprint2012arXiv

A general relativistic model of accretion disks with coronae surrounding Kerr black holes

We calculate the structure of a standard accretion disk with corona surrounding a massive Kerr black hole in general relativistic frame, in which the corona is assumed to be heated by the reconnection of the strongly buoyant magnetic fields generated in the cold accretion disk. The emergent spectra of the accretion disk-corona systems are calculated by using the relativistic ray-tracing method. We propose a new method to calculate the emergent Comptonized spectra from the coronae. The spectra of the disk-corona systems with a modified $α$-magnetic stress show that both the hard X-ray spectral index and the hard X-ray bolometric correction factor $L_{\rm bol}/L_{\rm X,2-10keV}$ increase with the dimensionless mass accretion rate, which are qualitatively consistent with the observations of active galactic nuclei (AGNs). The fraction of the power dissipated in the corona decreases with increasing black hole spin parameter $a$, which leads to lower electron temperatures of the coronas for rapidly spinning black holes. The X-ray emission from the coronas surrounding rapidly spinning black holes becomes weak and soft. The ratio of the X-ray luminosity to the optical/UV luminosity increases with the viewing angle, while the spectral shape in the X-ray band is insensitive with the viewing angle. We find that the spectral index in the infrared waveband depends on the mass accretion rate and the black hole spin $a$, which deviates from $f_ν\proptoν^{1/3}$ expected by the standard thin disk model.

preprint2012arXiv

Neutrino and anti-neutrino transport in accretion disks

We numerically solve the one dimensional Boltzmann equation of the neutrino and anti-neutrino transport in accretion disks and obtain the fully energy dependent and direction dependent neutrino and anti-neutrino emitting spectra, under condition that the distribution of the mass density,temperature and chemical components are given. Then, we apply the resulting neutrino and anti-neutrino emitting spectra to calculate the corresponding annihilation rate of neutrino pairs above the neutrino dominated accretion disk and find that the released energy resulting from the annihilation of neutrino pairs can not provide sufficient energy for the most energetic short gamma ray bursts whose isotropic luminosity can be as high as $10^{52}$ ergs/s unless the high temperature zone where the temperature is beyond 10 MeV can stretch over 200 km in the disk. We also compare the resulting luminosity of neutrinos and anti-neutrinos with the results from the two commonly used approximate treatment of the neutrino and anti-neutrino luminosity: the Fermi-Dirac black body limit and a simplified model of neutrino transport, i.e., the gray body model, and find that both of them overestimate the neutrino/anti-neutrino luminosity and their annihilation rate greatly. Additionally, as did in Sawyer (2003), we also check the validity of the two stream approximation, and find that it is a good approximation to high accuracy.

preprint2012arXiv

Vertical Structure of Neutrino Dominated Accretion Disks and Neutrino Transport in the disks

We investigate the vertical structure of neutrino dominated accretion disks by self-consistently considering the detailed microphysics, such as the neutrino transport, vertical hydrostatic equilibrium, the conservation of lepton number, as well as the balance between neutrino cooling, advection cooling and viscosity heating. After obtaining the emitting spectra of neutrinos and antineutrinos by solving the one dimensional Boltzmann equation of neutrino and antineutrino transport in the disk, we calculate the neutrino/antineutrino luminosity and their annihilation luminosity. We find that the total neutrino and antineutrino luminosity is about $10^{54}$ ergs/s and their annihilation luminosity is about $5\times10^{51}$ ergs/s with an extreme accretion rate $10 M_{\rm {sun}}$/s and an alpha viscosity $α=0.1$. In addition, we find that the annihilation luminosity is sensitive to the accretion rate and will not exceed $10^{50}$ ergs/s which is not sufficient to power the most fireball of GRBs, if the accretion rate is lower than $1 M_{\rm {sun}}$/s. Therefore, the effects of the spin of black hole or/and the magnetic field in the accretion flow might be introduced to power the central engine of GRBs.

preprint2011arXiv

A general relativistic external Compton-scattering model for TeV emission from M87

M87 is the first detected non-blazar extragalactic Tera-Electron-Volt (TeV) source with rapid variation and very flat spectrum in the TeV band. To explain the two-peaks in the spectral energy distribution (SED) of the nucleus of M87 which is similar to those of blazars, the most commonly adopted models are the synchrotron self-Compton scattering (SSC) models and the external inverse Compton scattering (EIC) models. Considering that there is no correlated variation in the soft band (from radio to X-ray) matching the TeV variation, and the TeV sources should not suffer from the gamma-gamma absorption due to the flat TeV spectrum, the EIC models are advantageous in modeling the TeV emission from M87. In this paper, we propose a self-consistent EIC model to explain the flat TeV spectrum of M87 within the framework of fully general relativity, where the background soft photons are from the advection-dominated accretion flow (ADAF) around the central black hole, and the high energy electrons are from the mini-jets which are powered by the magnetic reconnection in the main jet (Giannios et al. 2010). In our model, both the TeV flares observed in the years of 2005 and 2008 could be well explained: the gamma-gamma absorption for TeV photons is very low, even inside the region very close to the black hole 20Rg~50Rg; at the same region, the average EIC cooling time (~ 10^2-10^3s) is short, which is consistent with the observed time scale of TeV variation. Furthermore, we also discuss the possibility that the accompanying X-ray flare in 2008 is due to the direct synchrotron radiation of the mini-jets.

preprint2010arXiv

Emergent Spectra From Disks Surrounding Kerr Black Holes: Effect of Photon Trapping and Disk Self-Shadowing

Based on a new estimation of their thickness, the global properties of relativistic slim accretion disks are investigated in this work. The resulting emergent spectra are calculated using the relativistic ray-tracing method, in which we neglect the self-irradiation of the accretion disk. The angular dependence of the disk luminosity, the effects of the heat advection and the disk thickness on the estimation of the black hole spin are discussed. Compare to the previous works, our improvements are that we use the self-consistent disk equations and we consider the disk self-shadowing effect. We find that at the moderate accretion rate, the radiation trapped in the outer region of the accretion disks will escape in the inner region of the accretion disk and contribute to the emergent spectra. At the high accretion rate, for the large inclination and large black hole spin, both the disk thickness and the heat advection have significant influence on the emergent spectra. Consequently, these effects will influence the measurement of the black hole spin based on the spectra fitting and influence the angular dependence of the luminosity. For the disks around Kerr black holes with $a=0.98$, if the disk inclination is greater than $60^\circ$, and their luminosity is beyond 0.2 Eddington luminosity, the spectral model which is based on the relativistic standard accretion disk is no longer applicable for the spectra fitting. We also confirm that the effect of the self-shadowing is significantly enhanced by the light-bending, which implies that the non-relativistic treatment of the self-shadowing is inaccurate. According to our results, the observed luminosity dependence of the measured spin suggests that the disk self-shadowing significantly shapes the spectra of GRS 1915+105, which might lead to the underestimation of the black hole spin for the high luminosity states.

preprint2010arXiv

Implication of the Steady State Equilibrium Condition for Electron-Positron Gas in the Neutrino-driven Wind from Proto-Neutron Star

Based on the steady state equilibrium condition for neutron-proton-electron-positron gas in the neutrino-driven wind from protoneutron star, we estimate the initial electron fraction in the wind in a simple and effective way. We find that the condition in the wind might be propriate for the r-process nucleosynthesis.

preprint2010arXiv

Instabilities in the Gamma Ray Burst central engine. What makes the jet variable?

Both types of long and short gamma ray bursts involve a stage of a hyper-Eddington accretion of hot and dense plasma torus onto a newly born black hole. The prompt gamma ray emission originates in jets at some distance from this 'central engine' and in most events is rapidly variable, having a form of spikes and subpulses. This indicates at the variable nature of the engine itself, for which a plausible mechanism is an internal instability in the accreting flow. We solve numerically the structure and evolution of the neutrino-cooled torus. We take into account the detailed treatment of the microphysics in the nuclear equation of state that includes the neutrino trapping effect. The models are calculated for both Schwarzschild and Kerr black holes. We find that for sufficiently large accretion rates (> 10 Msun/s for non-rotating black hole, and >1 Msun/s for rotating black hole, depending on its spin), the inner regions of the disk become opaque, while the helium nuclei are being photodissociated. The sudden change of pressure in this region leads to the development of a viscous and thermal instability, and the neutrino pressure acts similarly to the radiation pressure in sub-Eddington disks. In the case of rapidly rotating black holes, the instability is enhanced and appears for much lower accretion rates. We also find the important and possibly further destabilizing role of the energy transfer from the rotating black hole to the torus via the magnetic coupling.

preprint2009arXiv

Constraining spins of supermassive black holes from TeV variability. II. fully general relativistic calculations

The fast variability of energetic TeV photons from the center of M87 has been detected, offering a new clue to estimate spins of supermassive black holes (SMBHs). We extend the study of Wang et al. (2008) by including all of general relativistic effects. We numerically solve the full set of relativistic hydrodynamical equations of the radiatively inefficient accretion flows (RIAFs) and then obtain the radiation fields around the black hole. The optical depth of the radiation fields to TeV photons due to pair productions are calculated in the Kerr metric. We find that the optical depth strongly depends on: (1) accretion rates as $\tautev\propto \dot{M}^{2.5-5.0}$; (2) black hole spins; and (3) location of the TeV source. Jointly considering the optical depth and the spectral energy distribution radiated from the RIAFs, the strong degeneration of the spin with the other free parameters in the RIAF model can be largely relaxed. We apply the present model to M87, wherein the RIAFs are expected to be at work, and find that the minimum specific angular momentum of the hole is $a\sim0.8$. The present methodology is applicable to M87-like sources with future detection of TeV emissions to constrain the spins of SMBHs.

preprint2009arXiv

Effect of Electron Screening on the Collapsing Process of Core-Collapse Supernovae

By using an average heavy nuclei model, the effects of the electron screening on electron capture (EC) in core-collapse supernovae are investigated. A one-dimension code based on the Ws15$M_{\sun}$ progenitor model is utilized to test the effects of electron screening during the collapsing process. The results show that, at high densities, the effects of EC on electron capture becomes significant. During the collapsing stage, the EC rate is decreased, the collapse timescale is prolonged and the leakage of the neutrino energy is increased. These effects lead to an appreciable decrease in the initial energy of the bounce shock wave. The effects of electron screeening in the other progenitor models are also discussed.

preprint2009arXiv

Images of the radiatively inefficient accretion flow surrounding a Kerr black hole: application in Sgr A*

In fully general relativity, we calculate the images of the radiatively inefficient accretion flow (RIAF) surrounding a Kerr black hole with arbitrary spins, inclination angles, and observational wavelengths. For the same initial conditions, such as the fixed accretion rate, it is found that the intrinsic size and radiation intensity of the images become larger, but the images become more compact in the inner region, while the size of the black hole shadow decreases with the increase of the black hole spin. With the increase of the inclination angles, the shapes of the black hole shadows change and become smaller, even disappear at all due to the obscuration by the thick disks. For median inclination angles, the radial velocity observed at infinity is larger because of both the rotation and radial motion of the fluid in the disk, which results in the luminous part of the images is much brighter. For larger inclination angles, such as the disk is edge on, the emission becomes dimmer at longer observational wavelengths (such as at 7.0mm and 3.5mm wavelengths), or brighter at shorter observational wavelengths (such as at 1.3 mm wavelength) than that of the face on case, except for the high spin and high inclination images. These complex behaviors are due to the combination of the Lorentz boosting effect and the radiative absorption in the disk. We hope our results are helpful to determine the spin parameter of the black hole in low luminosity sources, such as the Galactic center. A primary analysis by comparison with the observed sizes of Sgr A* at millimeters strongly suggests that the disk around the central black hole at Sgr A* is highly inclined or the central black hole is rotating fast.

preprint2009arXiv

Polarized Emission of Sagittarius A*

We explore the parameter space of the two temperature pseudo-Newtonian Keplerian accretion flow model for the millimeter and shorter wavelength emission from Sagittarius A*. A general relativistic ray-tracing code is used to treat the radiative transfer of polarized synchrotron emission from the flow. The synchrotron self-Comptonization and bremsstrahlung emission components are also included. It is shown that the model can readily account for the millimeter to sub-millimeter emission characteristics with an accretion rate of ~6x10^17g.s^-1 and an inclination angle of ~40 deg. However, the corresponding model predicted near-infrared and X-ray fluxes are more than one order of magnitude lower than the observed 'quiescent' state values. While the extended quiescent-state X-ray emission has been attributed to thermal emission from the large-scale accretion flow, the NIR emission and flares are likely dominated by emission regions either within the last stable orbit of a Schwarzschild black hole or associated with outflows. With the viscous parameter derived from numerical simulations, there is still a degeneracy between the electron heating rate and the magnetic parameter. A fully general relativistic treatment with the black hole spin incorporated will resolve these issues.

preprint2009arXiv

The role of black hole spin and magnetic field threading the unstable neutrino disk in Gamma Ray Bursts

We report on the third phase of our study of the neutrino-cooled hyperaccreting torus around a black hole that powers the jet in Gamma Ray Bursts. We focus on the influence of the black hole spin on the properties of the torus. The structure of a stationary torus around the Kerr black hole is solved numerically. We take into account the detailed treatment of the microphysics in the nuclear equation of state that includes the neutrino trapping effect. We find, that in the case of rapidly rotating black holes, the thermal instability discussed in our previous work is enhanced and develops for much lower accretion rates. We also find the important role of the energy transfer from the rotating black hole to the torus, via the magnetic coupling.