Source author record

Mouyuan Sun

Mouyuan Sun 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

22works
4topics
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

22 published item(s)

preprint2026arXiv

Light curve modeling of the nearest neutron star candidate LAMOST J235456.73+335625.9

The discovery of heavy radioactive elements (e.g., $^{60}\mathrm{Fe}$) on Earth suggests that supernova explosions may have occurred near our planet within the past million years, potentially having a significant impact on the ecological environment. This finding has motivated the search for nearby neutron stars in the Solar neighborhood. In a recent study, a candidate for one of the closest neutron stars to Earth, LAMOST J235456.73+335625.9 (hereafter J2354), was reported. Based on dynamical mass measurements under different inclination angle assumptions, the inferred mass range for the unseen compact companion in the system is $1.4$--$1.6$ $M_{\odot}$. Hence, the unseen companion in J2354 is either a massive cold white dwarf or a neutron star. Here we model the flux variations of J2354 as a combination of ellipsoidal modulation and surface spots. We test both cold spot and hot spot models, setting the number of spots to two in each case, and constrain the spot properties through light curve fitting. In the cold spot scenario, the spots are mostly visible at phases $0.5$--$0.75$, whereas in the hot spot scenario, the spots appear predominantly at phases $0.25$--$0.5$. The hot spot model shows better agreement with the observed H$α$ phase variation than the cold spot model. Furthermore, the thermal radiation of a massive but cold white dwarf cannot produce the level of localized heating required to explain the hot spot unless additional heating mechanisms are involved; in contrast, a neutron star can naturally provide such heating through energetic winds. Our results are consistent with the neutron star interpretation of the compact object in J2354.

preprint2022arXiv

Evidence for quasar fast outflows being accelerated at the scale of tens of parsecs

Quasar outflows may play a crucial role in regulating the host galaxy, although the spatial scale of quasar outflows remain a major enigma, with their acceleration mechanism poorly understood. The kinematic information of outflow is the key to understanding its origin and acceleration mechanism. Here, we report the galactocentric distances of different outflow components for both a sample and an individual quasar. We find that the outflow distance increases with velocity, with a typical value from several parsecs to more than one hundred parsecs, providing direct evidence for an acceleration happening at a scale of the order of 10 parsecs. These outflows carry ~1% of the total quasar energy, while their kinematics are consistent with a dust driven model with a launching radius comparable to the scale of a dusty torus, indicating that the coupling between dust and quasar radiation may produce powerful feedback that is crucial to galaxy evolution.

preprint2022arXiv

Nucleosynthesis Contribution of Neutrino-dominated Accretion Flows to the Chemical Evolution of Active Galactic Nuclei

Recent observations of quasars show high line-flux ratios in their broad emission lines and the ratios appear to be independent of redshift up to $z \gtrsim 6$, which indicate that the broad-line regions of these early quasars are surprisingly metal-rich. Here, we revisit the chemical evolution of high-redshift quasars by adding a new ingredient, i.e., the neutrino-dominated accretion flows (NDAFs) with outflows, on top of the conventional core-collapse supernovae (CCSNe). In the presence of the chemical contribution from NDAFs with outflows, the total metal mass (i.e., the summation of the conventional CCSN and NDAFs with outflows) per CCSN depends weakly upon the mass of the progenitor star if the mass is in the range of $\sim 25-55~M_{\odot}$. We model the chemical evolution by adopting a improved open-box model with three typical initial mass functions (IMFs). We find that, with the additional chemical contribution from NDAFs with outflows, the quasar metallicity can be enriched more rapidly in the very early Universe ($z \sim 10$) and reaches higher saturation than the no-NDAF case at $z \sim 8$, after which they evolve slowly with redshift. The quasar metallicity can reach $\sim 20~Z_{\odot}$ ($Z_\odot$ denotes the metallicity of the Sun; and $\sim 20\%$ of which is produced by NDAF outflows) at $z \sim 8$ for the ``top-heavy'' IMF model in \citet{Toyouchi2022}, which readily explains the quasar observations on the super-solar metal abundance and redshift-independent evolution.

preprint2022arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: UV-Optical Accretion Disk Measurements with Hubble Space Telescope

We present accretion-disk structure measurements from UV-optical reverberation mapping observations of a sample of eight quasars at 0.24<z<0.85. Ultraviolet photometry comes from two cycles of Hubble Space Telescope monitoring, accompanied by multi-band optical monitoring by the Las Cumbres Observatory network and Liverpool Telescopes. The targets were selected from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project sample with reliable black-hole mass measurements from Hbeta reverberation mapping results. We measure significant lags between the UV and various optical griz bands using JAVELIN and CREAM methods. We use the significant lag results from both methods to fit the accretion-disk structure using a Markov chain Monte Carlo approach. We study the accretion disk as a function of disk normalization, temperature scaling, and efficiency. We find direct evidence for diffuse nebular emission from Balmer and FeII lines over discrete wavelength ranges. We also find that our best-fit disk color profile is broadly consistent with the Shakura \& Sunyaev disk model. We compare our UV-optical lags to the disk sizes inferred from optical-optical lags of the same quasars and find that our results are consistent with these quasars being drawn from a limited high-lag subset of the broader population. Our results are therefore broadly consistent with models that suggest longer disk lags in a subset of quasars, for example, due to a nonzero size of the ionizing corona and/or magnetic heating contributing to the disk response.

preprint2022arXiv

Thermal Equilibrium Solutions of Black Hole Accretion Flows: Outflows VS Advection

Observations and numerical simulations have shown that outflows generally exist in the accretion process. We revisit the thermal equilibrium solutions of black hole accretion flows by including the role of outflows. Our study focuses on the comparison of the cooling rate of outflows with that of advection. Our results show that, except for the inner region, outflows can dominate over advection in a wide range of the flow, which is in good agreement with previous numerical simulations. We argue that an advection-dominated inner region together with an outflow-dominated outer region should be a general radial distribution for both super-Eddington accretion flows and optically thin flows with low accretion rates.

preprint2022arXiv

X-ray emission of contact binary variables within 1 kpc

By assembling the largest sample to date of X-ray emitting EW-type binaries (EWXs), we carried out correlation analyses for the X-ray luminosity log$L_{\textrm{X}}$, and X-ray activity level log($L_{\textrm{X}}$/$L_{\textrm{bol}}$) versus the orbital period $P$ and effective temperature $T_{\rm eff}$. We find strong $P$-log$L_{\textrm{X}}$ and $P$-log($L_{\textrm{X}}$/$L_{\textrm{bol}}$) correlations for EWXs with $P$ < 0.44 days and we provide the linear parametrizations for these relations, on the basis of which the orbital period can be treated as a good predictor for log$L_{\textrm{X}}$ and log($L_{\textrm{X}}$/$L_{\textrm{bol}}$). The aforementioned binary stellar parameters are all correlated with log$L_{\textrm{X}}$, while only $T_{\rm eff}$ exhibits a strong correlation with log($L_{\textrm{X}}$/$L_{\textrm{bol}}$). Then, EWXs with higher temperature show lower X-ray activity level, which could indicate the thinning of the convective area related to the magnetic dynamo mechanism. The total X-ray luminosity of an EWX is essentially consistent with that of an X-ray saturated main sequence star with the same mass as its primary, which may imply that the primary star dominates the X-ray emission. The monotonically decreasing $P$-log($L_{\textrm{X}}$/$L_{\textrm{bol}}$) relation and the short orbital periods indicate that EWXs could all be in the X-ray saturated state, and they may inherit the changing trend of the saturated X-ray luminosities along with the mass shown by single stars. For EWXs, the orbital period, mass, and effective temperature increase in concordance. We demonstrate that the period $P=0.44$ days corresponds to the primary mass of $\sim1.1 \rm M_\odot$, beyond which the saturated X-ray luminosity of single stars will not continue to increase with mass. This explains the break in the positive $P$-log$L_{\textrm{X}}$ relation for EWXs with $P>0.44$ days.

preprint2021arXiv

An Extraordinary Response of Iron Emission to the Central Outburst in a Tidal Disruption Event Candidate

Understanding the origin of \feii\ emission is important because it is crucial to construct the main sequence of Active Galactic Nuclei (AGNs). Despite several decades of observational and theoretical effort, the location of the optical iron emitting region and the mechanism responsible for the positive correlation between the \feii\ strength and the black hole accretion rate remain open questions as yet. In this letter, we report the optical \feii\ response to the central outburst in PS1-10adi, a candidate tidal disruption event (TDE) taking place in an AGN at $z = 0.203$ that has aroused extensive attention. For the first time, we observe that the \feii\ response in the rising phase of its central luminosity is significantly more prominent than that in the decline phase, showing a hysteresis effect. We interpret this hysteresis effect as a consequence of the gradual sublimation of the dust grains situating at the inner surface of the torus into gas when the luminosity of the central engine increases. It is the iron element released from the sublimated dust that contributes evidently to the observed \feii\ emission. This interpretation, together with the weak response of the \hb\ emission as we observe, naturally explains the applicability of relative \feii\ strength as a tracer of the Eddington ratio. In addition, optical iron emission of this origin renders the \feii\ time lag a potential "standard candle" with cosmological implications.

preprint2021arXiv

The Disk Veiling Effect of the Black Hole Low-Mass X-ray Binary A0620-00

The optical light curves of quiescent black hole low-mass X-ray binaries often exhibit significant non-ellipsoidal variabilities, showing the photospheric radiation of the companion star is veiled by other source of optical emission. Assessing this "veiling" effect is critical to the black hole mass measurement. Here in this work, we carry out a strictly simultaneous spectroscopic and photometric campaign on the prototype of black hole low-mass X-ray binary A0620-00. We find that for each observation epoch, the extra optical flux beyond a pure ellipsoidal modulation is positively correlated with the fraction of veiling emission, indicating the accretion disk contributes most of the non-ellipsoidal variations. Meanwhile, we also obtain a K2V spectral classification of the companion, as well as the measurements of the companion's rotational velocity $v \sin i = 83.8\pm1.9$ km s$^{-1}$ and the mass ratio between the companion and the black hole $q=0.063\pm0.004$.

preprint2020arXiv

AGNs are not that cool: revisiting the intrinsic AGN far-infrared spectral energy distribution

We investigate the intrinsic spectral energy distribution (SED) of active galactic nuclei (AGNs) at infrared (IR) bands with 42 $z < 0.5$ optically luminous Palomar Green survey quasars through SED decomposition. We decompose the SEDs of the 42 quasars by combining an AGN IR template library Siebenmorgen2015 that covers a wide range of the AGN parameter space with three commonly used galaxy template libraries. We determine the median AGN SED from the best-fitting results. The far-IR (FIR) contribution of our median AGN SED is significantly smaller than that of Symeonidis et al. 2016, but roughly consistent with that of Lyu et al. 2017. The AGN IR SED becomes cooler with increasing bolometric luminosity, which might be due to that more luminous AGNs might have stronger radiative feedback to change torus structures and/or their tori might have higher metallicities. Our conclusions do not depend on the choice of galaxy template libraries. However, since the predicted polycyclic aromatic hydrocarbon (PAH) emission line flux is galaxy template-dependent, cautions should be taken on deriving galaxy FIR contribution from PAH fluxes.

preprint2020arXiv

Corona-Heated Accretion-disk Reprocessing (CHAR): A Physical Model to Decipher the Melody of AGN UV/optical Twinkling

Active galactic nuclei (AGNs) have long been observed to "twinkle" (i.e., their brightness varies with time) on timescales from days to years in the UV/optical bands. Such AGN UV/optical variability is essential for probing the physics of supermassive black holes (SMBHs), the accretion disk, and the broad-line region. Here we show that the temperature fluctuations of an AGN accretion disk, which is magnetically coupled with the corona, can account for observed high-quality AGN optical light curves. We calculate the temperature fluctuations by considering the gas physics of the accreted matter near the SMBH. We find that the resulting simulated AGN UV/optical light curves share the same statistical properties as the observed ones as long as the dimensionless viscosity parameter $α$, which is widely believed to be controlled by magnetohydrodynamic (MHD) turbulence in the accretion disk, is about $0.01$---$0.2$. Moreover, our model can simultaneously explain the larger-than-expected accretion disk sizes and the dependence of UV/optical variability upon wavelength for NGC 5548. Our model also has the potential to explain some other observational facts of AGN UV/optical variability, including the timescale-dependent bluer-when-brighter color variability and the dependence of UV/optical variability on AGN luminosity and black hole mass. Our results also demonstrate a promising way to infer the black-hole mass, the accretion rate, and the radiative efficiency, thereby facilitating understanding of the gas physics and MHD turbulence near the SMBH and its cosmic mass growth history by fitting the AGN UV/optical light curves in the era of time-domain astronomy.

preprint2020arXiv

EUCLIA. II. On the puzzling large UV to X-ray lags in Seyfert galaxies

Recent intense X-ray and UV monitoring campaigns with Swift have detected clear UV lags behind X-ray in several local active galactic nuclei (AGNs). The UV to X-ray lags are often larger (by a factor up to ~ 20) than expected if the UV variation is simply due to the X-ray reprocessing. We previously developed a model in which the UV/optical variations are attributed to disk turbulences, and the effect of large-scale turbulence is considered. Our model, which overcomes many severe challenges to the reprocessing scheme, can well explain the observed variations in NGC 5548, particularly the correlations and lags among the UV/optical bands. In this work, assuming the corona heating is associated with turbulences in the inner accretion disk, we extend our study to model the correlations and lags between the X-ray and UV/optical bands. We find that our model, without the need of light echoing, can well reproduce the observed UV to X-ray lags and the optical to UV lags simultaneously in four local Seyfert galaxies, including NGC 4151, NGC 4395, NGC 4593, and NGC 5548. In our scenario, relatively larger UV to X-ray lag is expected for AGN with smaller innermost disk radius and thus more compact corona. Interestingly, for these Seyfert galaxies studied in this work, sources with relatively larger UV to X-ray lags do have broader Fe Ka lines, indicative of relativistic broadening due to more compact corona and smaller innermost disk radius. If confirmed with more X-ray and UV monitoring campaigns, this interesting discovery would provide a new probe to the inner disk/corona.

preprint2019arXiv

Mining for Candidates of Galactic Stellar-mass Black Hole Binaries with LAMOST

We study the prospects of searching for black hole (BH) binary systems with a stellar-mass BH and a non-compact visible companion, by utilizing the spectroscopic data of Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST). We simulate the Galactic BH binary population and determine its optical visibility by considering the stellar synthetic population model and the distributions of binary orbital parameters. By convolving the visibility of BH binaries with the LAMOST detection sensitivity, we predict that $\gtrsim$ 400 candidate BH binaries can be found by the low-resolution, non-time-domain survey, and $\sim$ 50-350 candidates by the LAMOST ongoing medium-resolution, time-domain spectroscopic survey. Most of the candidates are short-period (0.2-2 days) binaries with M-, K-, G-, or F-type companions, in which $\sim$ 47% have a mass function (the lower limit of the BH mass) larger than 3 $M_{\odot}$. By complementing the LAMOST spectroscopic data with other photometric/spectroscopic surveys or follow-up observations, these candidates could be confirmed. Therefore, by exploring the LAMOST data, we can enlarge the sample of dynamically confirmed BH binaries significantly, which can improve our understanding of the mass distribution of BHs and the stellar evolution model.

preprint2019arXiv

Understanding Broad Mg II Variability in Quasars with Photoionization: Implications for Reverberation Mapping and Changing-Look Quasars

The broad Mg II line in quasars has distinct variability properties compared with broad Balmer lines: it is less variable, and usually does not display a "breathing" mode, the increase in the average cloud distance when luminosity increases. We demonstrate that these variability properties of Mg II can be reasonably well explained by simple Locally Optimally Emitting Cloud (LOC) photoionization models, confirming earlier photoionization results. In the fiducial LOC model, the Mg II-emitting gas is on average more distant from the ionizing source than the H$α$/H$β$ gas, and responds with a lower amplitude to continuum variations. If the broad-line region (BLR) is truncated at a physical radius of $\sim 0.3$ pc (for a $10^{8.5}M_{\odot}$ BH accreting at Eddington ratio of 0.1), most of the Mg II flux will always be emitted near this outer boundary and hence will not display breathing. These results indicate that reverberation mapping results on broad Mg II, while generally more difficult to obtain due to the lower line responsivity, can still be used to infer the Mg II BLR size and hence black hole mass. But it is possible that Mg II does not have a well defined intrinsic BLR size-luminosity relation for individual quasars, even though a global one for the general population may still exist. The dramatic changes in broad H$α$/H$β$ emission in the observationally-rare changing-look quasars are fully consistent with photoionization responses to extreme continuum variability, and the LOC model provides natural explanations for the persistence of broad Mg II in changing-look quasars defined on H$α$/H$β$, and the rare population of broad Mg II emitters in the spectra of massive inactive galaxies.

preprint2016arXiv

On the origin of the HLX-1 outbursts

HLX-1, currently the best intermediate-mass black hole candidate, has undergone seven violent outbursts, each with a peak X-ray luminosity of $L_{\mathrm{peak},\mathrm{X}}\sim 10^{42}\ \rm{erg\ s^{-1}}$. Interestingly, the properties of the HLX-1 outbursts evolve with time. In this work, we aim to constrain the physical parameters of the central engine of the HLX-1 outbursts in the framework of the black hole accretion. We find that the physical properties of the HLX-1 outbursts are consistent with being driven by the radiation pressure instability. This scenario can explain the evolution of the recurrent timescales of the HLX-1 outbursts as a function of the durations.

preprint2015arXiv

Evolution in the Black Hole - Galaxy Scaling Relations and the Duty Cycle of Nuclear Activity in Star-Forming Galaxies

We measure the location and evolutionary vectors of 69 Herschel-detected broad-line active galactic nuclei (BLAGNs) in the M_BH-M_* plane. BLAGNs are selected from the COSMOS and CDF-S fields, and span the redshift range 0.2< z<2.1. Black-hole masses are calculated using archival spectroscopy and single-epoch virial mass estimators, and galaxy total stellar masses are calculated by fitting the spectral energy distribution (subtracting the BLAGN component). The mass-growth rates of both the black hole and galaxy are calculated using Chandra/XMM-Newton X-ray and Herschel far-infrared data, reliable measures of the BLAGN accretion and galaxy star formation rates, respectively. We use Monte Carlo simulations to account for biases in our sample, due to both selection limits and the steep slope of the massive end of the galaxy stellar-mass distribution. We find our sample is consistent with no evolution in the M_BH-M_* relation from z~2 to z~0. BLAGNs and their host galaxies which lie off the black hole mass galaxy total stellar mass relation tend to have evolutionary vectors anti-correlated with their mass ratios: that is, galaxies with over-massive (under-massive) black holes tend to have a low (high) ratio of the specific accretion rate to the specific star formation rate. We also use the measured growth rates to estimate the preferred AGN duty cycle for our galaxies to evolve most consistently onto the local M_BH-M_Bul relation. Under reasonable assumptions of exponentially declining star formation histories, the data suggest a non-evolving (no more than a factor of a few) BLAGN duty cycle among star-forming galaxies of 10% (1sigma range of 1-42% at z<1 and 2-34% at z>1).

preprint2015arXiv

The Accretion Wind Model of the Fermi Bubbles (II): Radiation

In a previous work, we have shown that the formation of the 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 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 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 CR electrons 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.

preprint2015arXiv

The Biases of Optical Line-Ratio Selection for Active Galactic Nuclei, and the Intrinsic Relationship between Black Hole Accretion and Galaxy Star Formation

We use 317,000 emission-line galaxies from the Sloan Digital Sky Survey to investigate line-ratio selection of active galactic nuclei (AGNs). In particular, we demonstrate that "star formation dilution" by HII regions causes a significant bias against AGN selection in low-mass, blue, star-forming, disk-dominated galaxies. This bias is responsible for the observed preference of AGNs among high-mass, green, moderately star-forming, bulge-dominated hosts. We account for the bias and simulate the intrinsic population of emission-line AGNs using a physically-motivated Eddington ratio distribution, intrinsic AGN narrow line region line ratios, a luminosity-dependent Lbol/L[OIII] bolometric correction, and the observed Mbh-sigma relation. These simulations indicate that, in massive (log(M*/Msun) > 10) galaxies, AGN accretion is correlated with specific star formation rate but is otherwise uniform with stellar mass. There is some hint of lower black hole occupation in low-mass (log(M*/Msun) < 10) hosts, although our modeling is limited by uncertainties in measuring and interpreting the velocity dispersions of low-mass galaxies. The presence of star formation dilution means that AGNs contribute little to the observed strong optical emission lines (e.g., [OIII] and Ha) in low-mass and star-forming hosts. However the AGN population recovered by our modeling indicates that feedback by typical (low- to moderate-accretion) low-redshift AGNs has nearly uniform efficiency at all stellar masses, star formation rates, and morphologies. Taken together, our characterization of the observational bias and resultant AGN occupation function suggest that AGNs are unlikely to be the dominant source of star formation quenching in galaxies, but instead are fueled by the same gas which drives star formation activity.

preprint2015arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: Ensemble Spectroscopic Variability of Quasar Broad Emission Lines

We explore the variability of quasars in the MgII and Hbeta broad emission lines and UV/optical continuum emission using the Sloan Digital Sky Survey Reverberation Mapping project (SDSS-RM). This is the largest spectroscopic study of quasar variability to date: our study includes 29 spectroscopic epochs from SDSS-RM over $6$ months, containing 357 quasars with MgII and 41 quasars with Hbeta . On longer timescales, the study is also supplemented with two-epoch data from SDSS-I/II. The SDSS-I/II data include an additional $2854$ quasars with MgII and 572 quasars with Hbeta. The MgII emission line is significantly variable ($Δf/f$ 10% on 100-day timescales), a necessary prerequisite for its use for reverberation mapping studies. The data also confirm that continuum variability increases with timescale and decreases with luminosity, and the continuum light curves are consistent with a damped random-walk model on rest-frame timescales of $\gtrsim 5$ days. We compare the emission-line and continuum variability to investigate the structure of the broad-line region. Broad-line variability shows a shallower increase with timescale compared to the continuum emission, demonstrating that the broad-line transfer function is not a $δ$-function. Hbeta is more variable than MgII (roughly by a factor of $1.5$), suggesting different excitation mechanisms, optical depths and/or geometrical configuration for each emission line. The ensemble spectroscopic variability measurements enabled by the SDSS-RM project have important consequences for future studies of reverberation mapping and black hole mass estimation of $1<z<2$ quasars.

preprint2015arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: First Broad-line Hbeta and MgII Lags at z>~0.3 from six-Month Spectroscopy

Reverberation mapping (RM) measurements of broad-line region (BLR) lags in z>0.3 quasars are important for directly measuring black hole masses in these distant objects, but so far there have been limited attempts and success given the practical difficulties of RM in this regime. Here we report preliminary results of 15 BLR lag measurements from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project, a dedicated RM program with multi-object spectroscopy designed for RM over a wide redshift range. The lags are based on the 2014 spectroscopic light curves alone (32 epochs over 6 months) and focus on the Hbeta and MgII broad lines in the 100 lowest-redshift (z<0.8) quasars included in SDSS-RM; they represent a small subset of the lags that SDSS-RM (including 849 quasars to z~4.5) is expected to deliver. The reported preliminary lag measurements are for intermediate-luminosity quasars at 0.3<~z<0.8, including 9 Hbeta lags and 6 MgII lags, for the first time extending RM results to this redshift-luminosity regime and providing direct quasar black hole mass estimates over ~ half of cosmic time. The MgII lags also increase the number of known MgII lags by several-fold, and start to explore the utility of MgII for RM at high redshift. The location of these new lags at higher redshifts on the observed BLR size-luminosity relationship is statistically consistent with previous Hbeta results at z<0.3. However, an independent constraint on the relationship slope at z>0.3 is not yet possible due to the limitations in our current sample. Our results demonstrate the general feasibility and potential of multi-object RM for z>0.3 quasars.

preprint2015arXiv

The Sloan Digital Sky Survey Reverberation Mapping Project: Post-Starburst Signatures in Quasar Host Galaxies at z < 1

Quasar host galaxies are key for understanding the relation between galaxies and the supermassive black holes (SMBHs) at their centers. We present a study of 191 broad-line quasars and their host galaxies at z < 1, using high signal-to-noise ratio (SNR) spectra produced by the Sloan Digital Sky Survey Reverberation Mapping project. Clear detection of stellar absorption lines allows a reliable decomposition of the observed spectra into nuclear and host components, using spectral models of quasar and stellar radiations as well as emission lines from the interstellar medium. We estimate age, mass (M*), and velocity dispersion (sigma*) of the host stars, the star formation rate (SFR), quasar luminosity, and SMBH mass (Mbh), for each object. The quasars are preferentially hosted by massive galaxies with M* ~ 10^{11} Msun characterized by stellar ages around a billion years, which coincides with the transition phase of normal galaxies from the blue cloud to the red sequence. The host galaxies have relatively low SFRs and fall below the main sequence of star-forming galaxies at similar redshifts. These facts suggest that the hosts have experienced an episode of major star formation sometime in the past billion years, which was subsequently quenched or suppressed. The derived Mbh - sigma* and Mbh - M* relations agree with our past measurements and are consistent with no evolution from the local Universe. The present analysis demonstrates that reliable measurements of stellar properties of quasar host galaxies are possible with high-SNR fiber spectra, which will be acquired in large numbers with future powerful instruments such as the Subaru Prime Focus Spectrograph.

preprint2014arXiv

Fermi Bubbles Inflated by Winds Launched from the Hot Accretion Flow in Sgr A*

A pair of giant gamma-ray bubbles have been revealed by the {\it Fermi} LAT. In this paper we investigate their formation mechanism. Observations have indicated that the activity of the supermassive black hole located at the Galactic center, Sgr A*, was much stronger than the present time. Specifically, one possibility is that while Sgr A* was also in the hot accretion regime, the accretion rate should be $10^3-10^4$ times higher during the past $\sim 10^7$ yr. On the other hand, recent MHD numerical simulations of hot accretion flows have unambiguously shown the existence of strong winds and obtained their properties. Based on these knowledge, by performing three-dimensional hydrodynamical simulations, we show in this paper that the Fermi bubbles could be inflated by winds launched from the ``past' hot accretion flow in Sgr A*. In our model, the active phase of Sgr A* is required to last for about 10 million years and it was quenched no more than 0.2 million years ago. The Central Molecular Zone (CMZ) is included and it collimates the wind orientation towards the Galactic poles. Viscosity suppresses the Rayleigh-Taylor and Kelvin-Helmholtz instabilities and results in the smoothness of the bubble edge. The main observational features of the bubbles can be well explained. Specifically, the {\it ROSAT} X-ray features are interpreted by the shocked interstellar medium and the interaction region between winds and CMZ gas. The thermal pressure and temperature obtained in our model are in good consistency with the recent {\it Suzaku} observations.

preprint2014arXiv

Origin of the Correlations Between Supermassive Black Holes and Their Host Galaxies

Observations have shown that supermassive black holes in nearby elliptical galaxies correlate tightly with the stellar velocity dispersion (the $\MBH - σ$ relation) and the stellar mass (the $\MBH - \Mhost$ relation) of their host spheroids. However, the origin of these correlations remains ambiguous. In a previous paper by Zhu et al., we proposed a model which links the M-$σ$ relation to the the dynamical state of the system and the $\MBH - \Mhost$ relation to the self-regulation of galaxy growth. To test this model, we compile a sample of observed galaxies with different properties and examine the dependence of the above correlations on these parameters. We find that galaxies that satisfy the the $\MBH - σ$ correlation appear to have reached virial equilibrium, as indicated by the ratio between kinetic energy and gravitational potential, 2K/U $\sim$ 1. Furthermore, the ratio of black hole accretion rate to star formation rate remains nearly constant, BHAR /SFR $\sim$ $10^{-3}$, in active galaxies over a wide range of mass in the redshift range z=0 - 3. These results confirm our theoretical model that the observed correlations have different origins: the $\MBH - σ$ relation may result from galaxy relaxation, while the $\MBH$ - $\Mhost$ relation may be due to self-regulated black hole accretion and star formation in galaxies.