Source author record

Jian-Min Wang

Jian-Min Wang 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

62works
11topics
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

62 published item(s)

preprint2024arXiv

AGN STORM 2: V. Anomalous Behavior of the CIV Light Curve in Mrk 817

An intensive reverberation mapping campaign on the Seyfert 1 galaxy Mrk817 using the Cosmic Origins Spectrograph (COS) on the Hubble Space Telescope (HST) revealed significant variations in the response of the broad UV emission lines to fluctuations in the continuum emission. The response of the prominent UV emission lines changes over a $\sim$60-day duration, resulting in distinctly different time lags in the various segments of the light curve over the 14 months observing campaign. One-dimensional echo-mapping models fit these variations if a slowly varying background is included for each emission line. These variations are more evident in the CIV light curve, which is the line least affected by intrinsic absorption in Mrk817 and least blended with neighboring emission lines. We identify five temporal windows with distinct emission line response, and measure their corresponding time delays, which range from 2 to 13 days. These temporal windows are plausibly linked to changes in the UV and X-ray obscuration occurring during these same intervals. The shortest time lags occur during periods with diminishing obscuration, whereas the longest lags occur during periods with rising obscuration. We propose that the obscuring outflow shields the ultraviolet broad lines from the ionizing continuum. The resulting change in the spectral energy distribution of the ionizing continuum, as seen by clouds at a range of distances from the nucleus, is responsible for the changes in the line response.

preprint2022arXiv

Accretion Disk Size Measurements of Active Galactic Nuclei Monitored by the Zwicky Transient Facility

We compile a sample of 92 active galactic nuclei (AGNs) at z<0.75 with $gri$ photometric light curves from the archival data of the Zwicky Transient Facility and measure the accretion disk sizes via continuum reverberation mapping. We employ Monte Carlo simulation tests to assess the influences of data sampling and broad emission lines and select out the sample with adequately high sampling cadences (3 days apart in average) and minimum contaminations of broad emission lines. The inter-band time delays of individual AGNs are calculated using the interpolated cross-correlation function and then these delays are fitted with a generalized accretion disk model, in which inter-band time delays are a power function of wavelength, black hole mass, and luminosity. A Markov-chain Monte Carlo method is adopted to determine the best parameter values. Overall the inter-band time delays can be fitted with the $τ\ \propto λ^{4/3}$ relation as predicted from a steady-state, optically thick, geometrically thin accretion disk, however, the yielded disk size is systematically larger than expected, although the ratio of the measured to theoretical disk sizes depend on using the emissivity -- or responsivity -- weighted disk radius. These results are broadly consistent with previous studies, all together raising a puzzle about the "standard" accretion disk model.

preprint2022arXiv

Cosmology Intertwined: A Review of the Particle Physics, Astrophysics, and Cosmology Associated with the Cosmological Tensions and Anomalies

In this paper we will list a few important goals that need to be addressed in the next decade, also taking into account the current discordances between the different cosmological probes, such as the disagreement in the value of the Hubble constant $H_0$, the $σ_8$--$S_8$ tension, and other less statistically significant anomalies. While these discordances can still be in part the result of systematic errors, their persistence after several years of accurate analysis strongly hints at cracks in the standard cosmological scenario and the necessity for new physics or generalisations beyond the standard model. In this paper, we focus on the $5.0\,σ$ tension between the {\it Planck} CMB estimate of the Hubble constant $H_0$ and the SH0ES collaboration measurements. After showing the $H_0$ evaluations made from different teams using different methods and geometric calibrations, we list a few interesting new physics models that could alleviate this tension and discuss how the next decade's experiments will be crucial. Moreover, we focus on the tension of the {\it Planck} CMB data with weak lensing measurements and redshift surveys, about the value of the matter energy density $Ω_m$, and the amplitude or rate of the growth of structure ($σ_8,fσ_8$). We list a few interesting models proposed for alleviating this tension, and we discuss the importance of trying to fit a full array of data with a single model and not just one parameter at a time. Additionally, we present a wide range of other less discussed anomalies at a statistical significance level lower than the $H_0$--$S_8$ tensions which may also constitute hints towards new physics, and we discuss possible generic theoretical approaches that can collectively explain the non-standard nature of these signals.[Abridged]

preprint2022arXiv

Detection of eccentric close-binary supermassive black holes with incomplete interferometric data

Recent studies have proposed that GRAVITY+ instrument is able to trace the circular orbit of the subparsec close-binary supermassive black holes (CB-SMBHs) by measuring the photocentre variation of the hot dust emission. However, the CB-SMBHs orbit may become highly eccentric throughout the evolution of these objects, and the orbital period may be far longer than the observational time baseline. We investigate the problem of detecting the CB-SMBH with hot dust emission and high eccentricity (eCBSMBH, e=0.5) when the observed time baselines of their astrometric data and radial velocities are considerably shorter than the orbital period. The parameter space of the Keplerian model of the eCBSMBH is large for exploratory purposes. We therefore applied the Bayesian method to fit orbital elements of the eCBSMBH to combined radial velocity and astrometric data covering a small fraction of the orbital period. We estimate that a number of potential eCBSMBH systems within reach of GRAVITY + will be similar to the number of the planned circular targets. We show that using observational time baselines that cover ~ 10% of the orbit increases the possibility of determining the period, eccentricity, and total mass of an eCBSMBH. When the observational time baseline becomes too short (~ 5%), the quality of the retrieved eCBSMBH parameters degrades. We also illustrate how interferometry may be used to estimate the photo-centre at the eCBSMBH emission line, which could be relevant for GRAVITY+ successors. Even if the astrometric signal for eCBSMBH systems is reduced by a factor of sqrt{1-e^{2}} compared to circular ones, we find that the hot dust emission of eCBSMBHs can be traced by GRAVITY+ at the elementary level.

preprint2022arXiv

Monitoring AGNs with H$β$ Asymmetry. III. Long-term Reverberation Mapping Results of 15 Palomar-Green Quasars

In this third paper of the series reporting on the reverberation mapping (RM) campaign of active galactic nuclei with asymmetric H$β$ emission-line profiles, we present results for 15 Palomar-Green (PG) quasars using spectra obtained between the end of 2016 to May 2021. This campaign combines long time spans with relatively high cadence. For 8 objects, both the time lags obtained from the entire light curves and the measurements from individual observing seasons are provided. Reverberation mapping of 9 of our targets has been attempted for the first time, while the results for 6 others can be compared with previous campaigns. We measure the H$β$ time lags over periods of years and estimate their black hole masses. The long duration of the campaign enables us to investigate their broad line region (BLR) geometry and kinematics for different years by using velocity-resolved lags, which demonstrate signatures of diverse BLR geometry and kinematics. The BLR geometry and kinematics of individual objects are discussed. In this sample, the BLR kinematics of Keplerian/virialized motion and inflow is more common than outflow.

preprint2022arXiv

Spectroastrometry and Reverberation Mapping: the Mass and Geometric Distance of the Supermassive Black Hole in the Quasar 3C 273

The quasar 3C 273 has been observed with infrared spectroastrometry (SA) on broad Pa$α$ line and optical reverberation mapping (RM) on broad H$β$ line. SA delivers information about the angular size and structure of the Pa$α$ broad-line region (BLR), while RM delivers information about the physical size and structure of the H$β$ BLR. Based on the fact that the two BLRs share the mass of the supermassive black hole (SMBH) and viewing inclination, a combination of SA and velocity-resolved RM (SARM) thereby allows us to simultaneously determine the SMBH mass and geometric distance through dynamically modeling the two BLRs. We construct a suite of dynamical models with different geometric configurations and apply a Bayesian approach to obtain the parameter inferences. Overall the obtained masses and distances are insensitive to specific BLR configurations but more or less depend on parameterizations of the vertical distributions. The most probable model, chosen in light of the Bayes factor, yields an angular-size distance of $\log\,(D_{\rm A}/{\rm Mpc}) = 2.83_{-0.28}^{+0.32}$ and SMBH mass of $\log\,(M_\bullet/M_\odot)=9.06_{-0.27}^{+0.21}$, which agrees with the relationships between SMBH masses and bulge properties. The BLRs have an inclination of $5_{-1}^{+1}$ degrees, consistent with that of the large-scale jet in 3C 273. Our approach reinforces the capability of SARM analysis to measure SMBH mass and distance of AGNs even though SA and RM observations are undertaken with different emission lines and/or in different periods.

preprint2022arXiv

Systematically smaller single-epoch quasar black hole masses using a radius-luminosity relationship corrected for spectral bias

Determining black hole masses and accretion rates with better accuracy and precision is crucial for understanding quasars as a population. These are fundamental physical properties that underpin models of active galactic nuclei. A primary technique to measure the black hole mass employs the reverberation mapping of low-redshift quasars, which is then extended via the radius-luminosity relationship for the broad-line region to estimate masses based on single-epoch spectra. An updated radius-luminosity relationship incorporates the flux ratio of optical Fe ii to H$β$ ($\equiv \mathcal{R}_{\rm Fe}$) to correct for a bias in which more highly accreting systems have smaller line-emitting regions than previously realized. In this current work, we demonstrate and quantify the effect of using this Fe-corrected radius-luminosity relationship on mass estimation by employing archival data sets possessing rest-frame optical spectra over a wide range of redshifts. We find that failure to use a Fe-corrected radius predictor results in overestimated single-epoch black hole masses for the most highly accreting quasars. Their accretion rate measures ($L_{\rm Bol}/ L_{\rm Edd}$ and $\dot{\mathscr{M}}$), are similarly underestimated. The strongest Fe-emitting quasars belong to two classes: high-z quasars with rest-frame optical spectra, which given their extremely high luminosities, require high accretion rates, and their low-z analogs, which given their low black holes masses, must have high accretion rates to meet survey flux limits. These classes have mass corrections downward of about a factor of two, on average. These results strengthen the association of the dominant Eigenvector 1 parameter $\mathcal{R}_{\rm Fe}$ with the accretion process.

preprint2020arXiv

A parallax distance to 3C 273 through spectroastrometry and reverberation mapping

Distance measurements for extragalactic objects are a fundamental problem in astronomy and cosmology. In the era of precision cosmology, we urgently need better measurements of cosmological distances to observationally test the increasing $H_{0}$ tension of the Hubble constant measured from different tools. Using spectroastrometry, GRAVITY at The Very Large Telescope Interferometer successfully revealed the structure, kinematics and angular sizes of the broad-line region (BLR) of 3C 273 with an unprecedentedly high spatial resolution. Fortunately, reverberation mapping (RM) of active galactic nuclei (AGNs) reliably provides linear sizes of their BLRs. Here we report a joint analysis of spectroastrometry and RM observations to measure AGN distances. We apply this analysis to 3C 273 observed by both GRAVITY and an RM campaign,and find an angular distance of $551.5_{-78.7}^{+97.3}\, {\rm Mpc}$ and $H_{0}=71.5_{-10.6}^{+11.9}\,{\rm km\,s^{-1}\,Mpc^{-1}}$. Advantages of the analysis are 1) its pure geometrical measurements and 2) it simultaneously yields mass of the central black hole in the BLR. Moreover, we can conveniently repeat measurements of selected AGNs to efficiently reduce the statistical and systematic errors. Future observations of a reasonably sized sample ($\sim 30$ AGNs) will provide distances of the AGNs and hence a new way of measuring $H_{0}$ with a high precision $\left(\lesssim 3\%\right)$ to test the $H_{0}$ tension.

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

Constraints on individual supermassive binary black holes using observations of PSR J1909-3744

We perform a search for gravitational waves (GWs) from several supermassive binary black hole (SMBBH) candidates (NGC 5548, Mrk 231, OJ 287, PG 1302-102, NGC 4151, Ark 120 and 3C 66B) in long-term timing observations of the pulsar PSR J1909$-$3744 obtained using the Parkes radio telescope. No statistically significant signals were found. We constrain the chirp masses of those SMBBH candidates and find the chirp mass of NGC 5548 and 3C 66B to be less than $2.4 \times 10^9\,\rm M_{\odot}$ and $2.5 \times 10^9\,\rm M_{\odot}$ (with 95% confidence), respectively. Our upper limits remain a factor of 3 to 370 above the likely chirp masses for these candidates as estimated from other approaches. The observations processed here provide upper limits on the GW strain amplitude that improve upon the results from the first Parkes Pulsar Timing Array data release by a factor of 2 to 7. We investigate how information about the orbital parameters can help improve the search sensitivity for individual SMBBH systems. Finally, we show that these limits are insensitive to uncertainties in the Solar System ephemeris model.

preprint2020arXiv

On the relation between hard X-ray photon index versus accretion rate for super-Eddington accreting quasars

We investigate whether the hard X-ray photon index ($Γ$) versus accretion rate correlation for super-Eddington accreting quasars is different from that for sub-Eddington accreting quasars. We construct a sample of 113 bright quasars from the Sloan Digital Sky Survey Data Release 14 quasar catalog, including 38 quasars as the super-Eddington subsample and 75 quasars as the sub-Eddington subsample. We derive black-hole masses using a simple-epoch virial mass formula based on the ${\rm Hβ}$ lines, and we use the standard thin disk model to derive the dimensionless accretion rates ($\dot{\mathscr{M}}$) for our sample. The X-ray data for these quasars are collected from the Chandra and XMM-Newton archives. We fit the hard X-ray spectra using a single power-law model to obtain $Γ$ values. We find a statistically significant ($R_{\rm S}=0.43$, $p=7.75\times{10}^{-3}$) correlation between $Γ$ and $\dot{\mathscr{M}}$ for the super-Eddington subsample. The $Γ$-$\dot{\mathscr{M}}$ correlation for the sub-Eddington subsample is also significant, but weaker ($R_{\rm S}=0.30$, $p=9.98\times{10}^{-3}$). Linear regression analysis shows that ${\rm Γ}=(0.34\pm0.11){\rm log}{\dot{\mathscr{M}}}+(1.71\pm0.17)$ and ${\rm Γ}=(0.09\pm0.04){\rm log}{\dot{\mathscr{M}}}+(1.93\pm0.04)$ for the super- and sub-Eddington subsamples, respectively. The $Γ$-$\dot{\mathscr{M}}$ correlations of the two subsamples are different, suggesting different disk-corona connections in these two types of systems. We propose one qualitative explanation of the steeper $Γ$-$\dot{\mathscr{M}}$ correlation in the super-Eddington regime that involves larger seed photon fluxes received by the compact coronae from the thick disks in super-Eddington accreting quasars.

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

Untangling Optical Emissions of the Jet and Accretion Disk in the Flat-Spectrum Radio Quasar 3C 273 with Reverberation Mapping Data

3C 273 is an intensively monitored flat-spectrum radio quasar with both a beamed jet and blue bump together with broad emission lines. The coexistence of the comparably prominent jet and accretion disk leads to complicated variability properties. Recent reverberation mapping monitoring for 3C 273 revealed that the optical continuum shows a distinct long-term trend that does not have a corresponding echo in the Hbeta fluxes. We compile multi-wavelength monitoring data from the Swift archive and other ground-based programs and clearly find two components of emissions at optical wavelength. One component stems from the accretion disk itself and the other component can be ascribed to the jet contribution, which also naturally accounts for the non-echoed trend in reverberation mapping data. We develop an approach to decouple the optical emissions from the jet and accretion disk in 3C 273 with the aid of multi-wavelength monitoring data. By assuming the disk emission has a negligible polarization in consideration of the low inclination of the jet, the results show that the jet contributes a fraction of ~10% at the minimum and up to ~40% at the maximum to the total optical emissions. This is the first time to provide a physical interpretation to the "detrending" manipulation conventionally adopted in reverberation mapping analysis. Our work also illustrates the importance of appropriately analyzing variability properties in cases of coexisting jets and accretion disks.

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

Kinematic signatures of reverberation mapping of close binaries of supermassive black holes in active galactic nuclei III. The case of elliptical orbits

An intense search for kinematic signatures of reverberation mapping of close binary SMBH (SMBBH) is currently ongoing. Elliptical SMBBH systems (i.e. both orbits and disc-like broad-line regions (BLR) are elliptic) have not been assessed in two-dimensional transfer functions (2DTF) studies. We aim to numerically infer such a 2DTF because the geometry of the unresolved region is imprinted on their optical emission. Through this, we determine their specific kinematical signatures. We simulated the geometry and kinematics of SMBBH whose components are on elliptical orbits. Each SMBH had a disc-like elliptical BLR. The SMBHs were active and orbited each other tightly at a subparsec distance. Here we calculate for the first time 2DTF, as defined in the velocity-time delay plane, for several elliptical configurations of SMBBH orbits and their BLRs. We find that these very complex configurations are clearly resolved in maps. These results are distinct from those obtained from circular and disc-wind geometry. We calculate the expected line variability for all SMBBH configurations. We show that the line shapes are influenced by the orbital phase of the SMBBH. Some line profiles resemble observed profiles, but they can also be much deformed to look like those from the disc-wind model. First, our results imply that using our 2DTF, we can detect and quantify kinematic signatures of elliptical SMBBH. Second, the calculated expected line profiles share some intriguing similarities with observed profiles, but also with some profiles that are synthesised in disc-wind models. To overcome the non-uniqueness of the spectral line shapes as markers of SMBBH, they must be accompanied with 2DTF.

preprint2019arXiv

Kinematic signatures of reverberation mapping of close binaries of supermassive black holes in active galactic nuclei. II. Atlas of two-dimensional transfer functions

Most large galaxies harbor supermassive black holes (SMBHs) in their centers, and galaxies merge. Consequently, binary SMBHs should be common in galactic nuclei. However, close binaries of SMBH (CB-SMBHs) with sub-parsec separation cannot be imaged directly using current facilities. Some indirect signatures, such as periodic signals in light curves and double peaks in emission-line profile, have been used to find CB-SMBH candidates, but ambiguities still exist and no definitive conclusions can be made. We have recently proposed a new method focusing on kinematic signatures that can be derived from reverberation mapping of CB-SMBHs, one that offers a promising avenue to address this important problem. In this paper, we calculated models for a wide range of parameters, but BLRs of two BHs are close but still not merged. The purpose of this supplementary paper is to provide an atlas of two-dimensional transfer functions of CB-SMBHs with a wide range of orbital and geometrical parameters to aid more efficient identification of CB-SMBH candidates in reverberation mapping data.

preprint2019arXiv

The Radius-Luminosity Relationship Depends on Optical Spectra in Active Galactic Nuclei

The radius-luminosity (R-L) relationship of active galactic nuclei (AGNs) established by the reverberation mapping (RM) observations has been widely used as a single-epoch black hole mass estimator in the research of large AGN samples. However, the recent RM campaigns discovered that the AGNs with high accretion rates show shorter time lags by factors of a few comparing with the predictions from the R-L relationship. The explanation of the shortened time lags has not been finalized yet. We collect 8 different single-epoch spectral properties to investigate how the shortening of the time lags correlate with those properties and to understand what is the origin of the shortened lags. We find that the flux ratio between Fe II and H$β$ emission lines shows the most prominent correlation, thus confirm that accretion rate is the main driver for the shortened lags. In addition, we establish a new scaling relation including the relative strength of Fe II emission. This new scaling relation can provide less biased estimates of the black hole mass and accretion rate from the single-epoch spectra of AGNs.

preprint2016arXiv

A Non-parametric Approach to Constrain the Transfer Function in Reverberation Mapping

Broad emission lines of active galactic nuclei stem from a spatially extended region (broad-line region, BLR) that is composed of discrete clouds and photoionized by the central ionizing continuum. The temporal behaviors of these emission lines are blurred echoes of the continuum variations (i.e., reverberation mapping, RM) and directly reflect the structures and kinematic information of BLRs through the so-called transfer function (also known as the velocity-delay map). Based on the previous works of Rybicki & Press (1992) and Zu et al. (2011), we develop an extended, non-parametric approach to determine the transfer function for RM data, in which the transfer function is expressed as a sum of a family of relatively displaced Gaussian response functions. Therefore, arbitrary shapes of transfer functions associated with complicated BLR geometry can be seamlessly included, enabling us to relax the presumption of a specified transfer function frequently adopted in previous studies and to let it be determined by observation data. We formulate our approach in a previously well-established framework that incorporates the statistical modeling of the continuum variations as a damped random walk process and takes into account the long-term secular variations which are irrelevant to RM signals. The Application to RM data shows the fidelity of our approach.

preprint2016arXiv

A note on periodicity of long-term variations of optical continuum in active galactic nuclei

Graham et al. found a sample of active galactic nuclei (AGNs) and quasars from the Catalina Real-time Transient Survey (CRTS) that have long-term periodic variations in optical continuum, the nature of the periodicity remains uncertain. We investigate the periodic variability characteristics of the sample by testing the relations of the observed variability periods with AGN optical luminosity, black hole mass and accretion rates, and find no significant correlations. We also test the observed periods in several different aspects related to accretion disks surrounding single black holes, such as the Keplerian rotational periods of 5100~Å photon-emission regions and self-gravity dominated regions and the precessing period of warped disks. These tests shed new lights on understanding AGN variability in general. Under the assumption that the periodic behavior is associated with SMBHB systems in particular, we compare the separations ($\mathscr{D}_{\bullet}$) against characteristic radii of broad-line regions ($R_{\rm BLR}$) of the binaries and find $\mathscr{D}_{\bullet}$$\approx$ 0.05 $R_{\rm BLR}$. This interestingly implies that these binaries have only circumbinary BLRs.

preprint2016arXiv

Improving the Flux Calibration in Reverberation Mapping by Spectral Fitting: Application to the Seyfert Galaxy MCG--6-30-15

We propose a method for the flux calibration of reverberation mapping spectra based on accurate measurement of [O III] $λ5007$ emission by spectral fitting. The method can achieve better accuracy than the traditional method of van Groningen & Wanders (1992), allowing reverberation mapping measurements for object with variability amplitudes as low as $\sim$ 5%. As a demonstration, we reanalyze the data of the Seyfert 1 galaxy MCG--6-30-15 taken from the 2008 campaign of the Lick AGN Monitoring Project, which previously failed to obtain a time lag for this weakly variable object owing to a relatively large flux calibration uncertainty. We detect a statistically significant rest-frame time lag of $6.38_{-2.69}^{+3.07}$ days between the H$β$ and $V$-band light curves. Combining this lag with FWHM(H$β$) = $1933\pm81$ $\rm km~s^{-1}$ and a virial coefficient of $f$ = 0.7, we derive a virial black hole mass of $3.26_{-1.40}^{+1.59}\times10^6$ $M_{\odot}$, which agrees well with previous estimates by other methods.

preprint2016arXiv

Reverberation Mapping of the Broad-line Region in NGC 5548: Evidence for Radiation Pressure?

NGC 5548 is the best-observed reverberation-mapped active galactic nucleus with long-term, intensive monitoring. Here we report results from a new observational campaign between January and July, 2015. We measure the centroid time lag of the broad H$β$ emission line with respect to the 5100 Åcontinuum and obtain $τ_{\rm cent} = 7.20^{+1.33}_{-0.35}$ days in the rest frame. This yields a black hole mass of $M_{\bullet}=8.71^{+3.21}_{-2.61} $x$ 10^{7}M_{\odot}$ using a broad H$β$ line dispersion of $3124\pm302$ km s$^{-1}$ and a virial factor of $f_{_{\rm BLR}}=6.3\pm1.5$ for the broad-line region (BLR), consistent with the mass measurements from previous H$β$ campaigns. The high-quality data allow us to construct a velocity-binned delay map for the broad H$β$ line, which shows a symmetric response pattern around the line center, a plausible kinematic signature of virialized motion of the BLR. Combining all the available measurements of H$β$ time lags and the associated mean 5100 Å luminosities over 18 campaigns between 1989 and 2015, we find that the H$β$ BLR size varies with the mean optical luminosity, but, interestingly, with a possible delay of $2.35_{-1.25}^{+3.47}$ yrs. This delay coincides with the typical BLR dynamical timescale of NGC 5548, indicating that the BLR undergoes dynamical changes, possibly driven by radiation pressure.

preprint2016arXiv

Spectroscopic Indication for a Centi-parsec Supermassive Black Hole Binary in the Galactic Center of NGC 5548

As a natural consequence of cosmological hierarchical structure formation, sub-parsec supermassive black hole binaries (SMBHBs) should be common in galaxies but thus far have eluded spectroscopic identification. Based on four decades of optical spectroscopic monitoring, we report that the nucleus of NGC 5548, a nearby Seyfert galaxy long suspected to have experienced a major merger about one billion years ago, exhibits long-term variability with a period of 14 years in the optical continuum and broad Hbeta emission line. Remarkably, the double-peaked profile of Hbeta shows systematic velocity changes with a similar period. These pieces of observations plausibly indicate that a SMBHB resides in the center of NGC 5548. The complex, secular variations in the line profiles can be explained by orbital motion of a binary with equal mass and a semi-major axis of 22 light-days (corresponding to 18 milli-parsec). At a distance of 75 Mpc, NGC 5548 is one of the nearest sub-parsec SMBHB candidates that offers an ideal laboratory for gravitational wave detection.

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.

preprint2016arXiv

Test of the Formation mechanism of the Broad Line Region in Active Galactic Nuclei

The origin of the Broad Line Region (BLR) in active galaxies remains unknown. It seems to be related to the underlying accretion disk but an efficient mechanism is required to rise the material from the disk surface without giving too strong signatures of the outflow in the case of the low ionization lines. We discuss in detail two proposed mechanisms: (i) radiation pressure acting on dust in the disk atmosphere creating a failed wind (ii) the gravitational instability of the underlying disk. We compare the predicted location of the inner radius of the BLR in those two scenarios with the observed position obtained from the reverberation studies of several active galaxies. The failed dusty outflow model well represents the observational data while the predictions of the self-gravitational instability are not consistent with observations. The issue remains why actually we do not see any imprints of the underlying disk instability in the BLR properties.

preprint2016arXiv

The Detailed Science Case for the Maunakea Spectroscopic Explorer: the Composition and Dynamics of the Faint Universe

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

preprint2016arXiv

The Fundamental Plane of the Broad-line Region in Active Galactic Nuclei

Broad emission lines in active galactic nuclei (AGNs) mainly arise from gas photoionized by continuum radiation from an accretion disk around a central black hole. The shape of the broad-line profile, described by ${\cal D}_{_{\rm Hβ}}={\rm FWHM}/σ_{_{\rm Hβ}}$, the ratio of full width at half maximum to the dispersion of broad H$β$, reflects the dynamics of the broad-line region (BLR) and correlates with the dimensionless accretion rate ($\dot{\mathscr{M}}$) or Eddington ratio ($L_{\rm bol}/L_{\rm Edd}$). At the same time, $\dot{\mathscr{M}}$ and $L_{\rm bol}/L_{\rm Edd}$ correlate with ${\cal R}_{\rm Fe}$, the ratio of optical Fe II to H$β$ line flux emission. Assembling all AGNs with reverberation mapping measurements of broad H$β$, both from the literature and from new observations reported here, we find a strong bivariate correlation of the form $\log(\dot{\mathscr{M}},L_{\rm bol}/L_{\rm Edd})=α+β{\cal D}_{_{\rm Hβ}}+γ{\cal R}_{\rm Fe},$ where $α=(2.47,0.31)$, $β=-(1.59,0.82)$ and $γ=(1.34,0.80)$. We refer to this as the fundamental plane of the BLR. We apply the plane to a sample of $z < 0.8$ quasars to demonstrate the prevalence of super-Eddington accreting AGNs are quite common at low redshifts.

preprint2015arXiv

A bias-free quantum random number generation using photon arrival time selectively

We present a high-quality, bias-free quantum random number generator (QRNG) using photon arrival time selectively in accordance with the number of photon detection events within a sampling time interval in attenuated light. It is well showed in both theoretical analysis and experiments verification that this random number production method eliminates both bias and correlation perfectly without more post processing and the random number can clearly pass the standard randomness tests. We fulfill theoretical analysis and experimental verification of the method whose rate can reach up to 45Mbps.

preprint2015arXiv

Alignments Of Black Holes With Their Warped Accretion Disks And Episodic Lifetimes Of Active Galactic Nuclei

Warped accretion disks have attracted intensive attention because of their critical role on shaping the spin of supermassive massive black holes (SMBHs) through the Bardeen-Petterson effect, a general relativistic effect that leads to final alignments or anti-alignments between black holes and warped accretion disks. We study such alignment processes by explicitly taking into account the finite sizes of accretion disks and the episodic lifetimes of AGNs that delineate the duration of gas fueling onto accretion disks. We employ an approximate global model to simulate the evolution of accretion disks, allowing to determine the gravitomagnetic torque that drives the alignments in a quite simple way. We then track down the evolutionary paths for mass and spin of black holes both in a single activity episode and over a series of episodes. Given with randomly and isotropically oriented gas fueling over episodes, we calculate the spin evolution with different episodic lifetimes and find that it is quite sensitive to the lifetimes. We therefore propose that spin distribution of SMBHs can place constraints on the episodic lifetimes of AGNs and vice versa. Applications of our results on the observed spin distributions of SMBHs and the observed episodic lifetimes of AGNs are discussed, although both the measurements at present are yet ambiguous to draw a firm conclusion. Our prescription can be easily incorporated into semi-analytic models for black hole growth and spin evolution.

preprint2015arXiv

Supermassive Black Holes with High Accretion Rates in Active Galactic Nuclei. III. Detection of Fe II Reverberation in Nine Narrow-Line Seyfert 1 Galaxies

This is the third in a series of papers reporting on a large reverberation-mapping campaign aimed to study the properties of active galactic nuclei (AGNs) with high accretion rates. We present new results on the variability of the optical Fe II emission lines in 10 AGNs observed by the Yunnan Observatory 2.4m telescope during 2012--2013. We detect statistically significant time lags, relative to the AGN continuum, in nine of the sources. This accurate measurement is achieved by using a sophisticated spectral fitting scheme that allows for apparent flux variations of the host galaxy, and several narrow lines, due to the changing observing conditions. Six of the newly detected lags are indistinguishable from the Hbeta lags measured in the same sources. Two are significantly longer and one is slightly shorter. Combining with Fe II lags reported in previous studies, we find a Fe II radius--luminosity relationship similar to the one for Hbeta, although our sample by itself shows no clear correlation. The results support the idea that Fe II emission lines originate in photoionized gas which, for the majority of the newly reported objects, is indistinguishable from the Hbeta-emitting gas. We also present a tentative correlation between the lag and intensity of Fe II and Hbeta and comment on its possible origin.

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

A Bayesian Approach to Estimate the Size and Structure of the Broad-line Region in Active Galactic Nuclei Using Reverberation Mapping Data

This is the first paper in a series devoted to systematic study of the size and structure of the broad-line region (BLR) in active galactic nuclei (AGNs) using reverberation mapping (RM) data. We employ a recently developed Bayesian approach that statistically describes the variabibility as a damped random walk process and delineates the BLR structure using a flexible disk geometry that can account for a variety of shapes, including disks, rings, shells, and spheres. We allow for the possibility that the line emission may respond non-linearly to the continuum, and we detrend the light curves when there is clear evidence for secular variation. We use a Markov Chain Monte Carlo implementation based on Bayesian statistics to recover the parameters and uncertainties for the BLR model. The corresponding transfer function is obtained self-consistently. We tentatively constrain the virial factor used to estimate black hole masses; more accurate determinations will have to await velocity-resolved RM data. Application of our method to RM data with Hbeta monitoring for about 40 objects shows that the assumed BLR geometry can reproduce quite well the observed emission-line fluxes from the continuum light curves. We find that the Hbeta BLR sizes obtained from our method are on average ~20% larger than those derived from the traditional cross-correlation method. Nevertheless, we still find a tight BLR size-luminosity relation with a slope of alpha=0.55\pm0.03 and an intrinsic scatter of ~0.18 dex. In particular, we demonstrate that our approach yields appropriate BLR sizes for some objects (such as Mrk 142 and PG 2130+099) where traditional methods previously encountered difficulties.

preprint2014arXiv

A Bayesian Method for the Intercalibration of Spectra in Reverberation Mapping

Flux calibration of spectra in reverberation mapping (RM) is most often performed by assuming the flux constancy of some specified narrow emission lines, which stem from an extended region that is sometimes partially spatially resolved, in contrast to the point-like broad-line region and the central continuum source. The inhomogeneous aperture geometries used among different observation sets in a joint monitoring campaign introduce systematic deviations to the fluxes of broad lines and central continuum, and intercalibration over these data sets is required. As an improvement to the previous empirical correction performed by comparing the (nearly) contemporaneous observation points, we describe a feasible Bayesian method that obviates the need for (nearly) contemporaneous observations, naturally incorporates physical models of flux variations, and fully takes into account the measurement errors. In particular, it fits all the data sets simultaneously regardless of samplings and makes use of all of the information in the data sets. A Markov Chain Monte Carlo implementation is employed to recover the parameters and uncertainties for intercalibration. Application to the RM data sets of NGC 5548 with joint monitoring shows the high fidelity of our method.

preprint2014arXiv

A New Approach to Constrain Black Hole Spins in Active Galaxies Using Optical Reverberation Mapping

A tight relation between the size of the broad-line region (BLR) and optical luminosity has been established in about 50 active galactic nuclei studied through reverberation mapping of the broad Hbeta emission line. The R_blr-L relation arises from simple photoionization considerations. Using a general relativistic model of an optically thick, geometrically thin accretion disk, we show that the ionizing luminosity jointly depends on black hole mass, accretion rate, and spin. The non-monotonic relation between the ionizing and optical luminosity gives rise to a complicated relation between the BLR size and the optical luminosity. We show that the reverberation lag of Hbeta to the varying continuum depends very sensitively to black hole spin. For retrograde spins, the disk is so cold that there is a deficit of ionizing photons in the BLR, resulting in shrinkage of the hydrogen ionization front with increasing optical luminosity, and hence shortened Hbeta lags. This effect is specially striking for luminous quasars undergoing retrograde accretion, manifesting in strong deviations from the canonical R_blr-L relation. This could lead to a method to estimate black hole spins of quasars and to study their cosmic evolution. At the same time, the small scatter of the observed R_blr-L relation for the current sample of reverberation-mapped active galaxies implies that the majority of these sources have rapidly spinning black holes.

preprint2014arXiv

Possible hard X-ray shortages in bursts from KS 1731-260 and 4U 1705-44

Aims: A hard X-ray shortage, implying the cooling of the corona, was observed during bursts of IGR J17473-272, 4U 1636-536, Aql X-1, and GS 1826-238. Apart from these four sources, we investigate here an atoll sample, in which the number of bursts for each source is larger than 5, to explore the possible additional hard X-ray shortage during {\it Rossi X-ray timing explorer (RXTE)} era. Methods: According to the source catalog that shows type-I bursts, we analyzed all the available pointing observations of these sources carried out by the {\it RXTE} proportional counter array (PCA). We grouped and combined the bursts according to their outburst states and searched for the possible hard X-ray shortage while bursting. Results: We found that the island states of KS 1731-260 and 4U 1705-44 show a hard X-ray shortage at significant levels of 4.5 and 4.7 $σ$ and a systematic time lag of $0.9 \pm 2.1$ s and $2.5 \pm 2.0$ s with respect to the soft X-rays, respectively. While in their banana branches and other sources, we did not find any consistent shortage.

preprint2014arXiv

Supermassive black holes with high accretion rates in active galactic nuclei: II. the most luminous standard candles in the Universe

This is the second in a series of papers reporting on a large reverberation mapping (RM) campaign to measure black hole (BH) mass in high accretion rate active galactic nuclei (AGNs). The goal is to identify super-Eddington accreting massive black holes (SEAMBHs) and to use their unique properties to construct a new method for measuring cosmological distances. Based on theoretical models, the saturated bolometric luminosity of such sources is proportional to the BH mass which can be used to obtain their distance. Here we report on five new RM measurements and show that in four of the cases we can measure the BH mass and three of these sources are SEAMBHs. Together with the three sources from our earlier work, we now have six new sources of this type. We use a novel method based on a minimal radiation efficiency to identify nine additional SEAMBHs from earlier RM-based mass measurements. We use a Bayesian analysis to determine the parameters of the new distance expression, and the method uncertainties, from the observed properties of the objects in the sample. The ratio of the newly measured distances to the standard cosmological ones has a mean scatter of 0.14 dex, indicating that SEAMBHs can be use as cosmological distance probes. With their high luminosity, long period of activity and large numbers at high redshifts, SEAMBHs have a potential to extend the cosmic distance ladder beyond the range now explored by type Ia supernovae.

preprint2013arXiv

Evolution of Warped Accretion Disks in Active Galactic Nuclei. I. Roles of Feeding at the Outer Boundaries

We investigate the alignment processes of spinning black holes and their surrounding warped accretion disks in a frame of two different types of feeding at the outer boundaries. We consider (1) fixed flows in which gas is continually fed with a preferred angular momentum, and (2) free flows in which there is no gas supply and the disks diffuse freely at their outer edges. As expected, we find that for the cases of fixed flows the black hole disk systems always end up aligning on timescales of several 1e6 yr, irrespective of the initial inclinations. If the initial inclination angles are larger than pi/2, the black hole accretion transits from retrograde to prograde fashion, and the accreted mass onto the black holes during these two phases is comparable. On the other hand, for the cases of free flows, both alignments and anti-alignments can occur, depending on the initial inclinations and the ratios of the angular momentum of the disks to that of the black holes. In such cases, the disks will be consumed within timescales of 1e6 yr by black holes accreting at the Eddington limit. We propose that there is a close connection between the black hole spin and the lifetime for which the feeding persists, which determines the observable episodic lifetimes of active galactic nuclei. We conclude that careful inclusion of the disk feeding at the outer boundaries is crucial for modeling the evolution of the black hole spin.

preprint2013arXiv

Outflows from active galactic nuclei: The BLR-NLR metallicity correlation

The metallicity of active galactic nuclei (AGNs), which can be measured by emission line ratios in their broad and narrow line regions (BLRs and NLRs), provides invaluable information about the physical connection between the different components of AGNs. From the archival databases of the International Ultraviolet Explorer, the Hubble Space Telescope and the Sloan Digital Sky Survey, we have assembled the largest sample available of AGNs which have adequate spectra in both the optical and ultraviolet bands to measure the narrow line ratio [N II]/Hα and also, in the same objects, the broad-line N V/C IV ratio. These permit the measurement of the metallicities in the NLRs and BLRs in the same objects. We find that neither the BLR nor the NLR metallicity correlate with black hole masses or Eddington ratios, but there is a strong correlation between NLR and BLR metallicities. This metallicity correlation implies that outflows from BLRs carry metal-rich gas to NLRs at characteristic radial distances of ~ 1.0 kiloparsec. This chemical connection provides evidence for a kinetic feedback of the outflows to their hosts. Metals transported into the NLR enhance the cooling of the ISM in this region, leading to local star formation after the AGNs turn to narrow line LINERs. This post-AGN star formation is predicted to be observable as an excess continuum emission from the host galaxies in the near infrared and ultraviolet, which needs to be further explored.

preprint2013arXiv

The hard X-ray behavior of Aql X-1 during type-I bursts

We report the discovery of an anti-correlation between the soft and the hard X-ray lightcurves of the X-ray binary Aql X-1 when bursting. This behavior may indicate that the corona is cooled by the soft X-ray shower fed by the type-I X-ray bursts, and that this process happens within a few seconds. Stacking the Aql X-1 lightcurves of type-I bursts, we find a shortage in the 40--50 keV band, delayed by 4.5$\pm$1.4 s with respect to the soft X-rays. The photospheric radius expansion (PRE) bursts are different in that neither a shortage nor an excess shows up in the hard X-ray lightcurve.

preprint2013arXiv

Type-I X-ray bursts reveal a fast co-evolving behavior of the corona in an X-ray binary

The coronae in X-ray binaries (XRBs) still remain poorly understood, although they have been believed for a long time to play a key role in modeling the characteristic outbursts of XRBs. Type-I X-ray bursts, the thermonuclear flashes happening on the surface of a neutron star (NS), can be used as a probe to the innermost region of a NS XRB, where the corona is believed to be located very close to the NS. We report the discovery of a tiny life cycle of the corona that is promptly co-evolved with the type-I bursts superimposed on the outburst of the NS XRB IGR J17473$-$2721. This finding may serve as the first evidence of directly seeing the rapid disappearance and formation of a corona in an XRB with a cooling/heating timescale of less than a second, which can strongly constrain the accretion models in XRBs at work.

preprint2013arXiv

X-ray bursts as a probe of the corona: the case of XRB 4U 1636-536

To investigate the possible cooling of the corona by soft X-rays bursts, we have studied 114 bursts embedded in the known X-ray evolution of 4U 1636-536. We have grouped these bursts according to the ratio of the flux in the 1.5--12 keV band with respect to that in the 15--50 keV band, as monitored by RXTE/ASM and Swift/BAT, respectively. We have detected a shortage at hard X-rays while bursting. This provides hints for a corona cooling process driven by soft X-rays fed by the bursts that occurred on the surface of neutron star. The flux shortage at 30--50 keV has a time lag of 2.4$\pm$1.5 seconds with respect to that at 2--10 keV, which is comparable to that of 0.7$\pm$0.5 seconds reported in bursts of IGR 17473-2721. We comment on the possible origin of these phenomena and on the implications for the models on the location of the corona.

preprint2012arXiv

Cosmological Evolution of Supermassive Black Holes. II. Evidence for Downsizing of Spin Evolution

The spin is an important but poorly constrained parameter for describing supermassive black holes (SMBHs). Using the continuity equation of SMBH number density, we explicitly obtain the mass-dependent cosmological evolution of the radiative efficiency for accretion, which serves as a proxy for SMBH spin. Our calculations make use of the SMBH mass function of active and inactive galaxies (derived in the first paper of this series), the bolometric luminosity function of active galactic nuclei (AGNs), corrected for the contribution from Compton-thick sources, and the observed Eddington ratio distribution. We find that the radiative efficiency generally increases with increasing black hole mass at high redshifts (z>~1), roughly as η\propto M_bh^0.5, while the trend reverses at lower redshifts, such that the highest efficiencies are attained by the lowest mass black holes. Black holes with M_bh>~10^8.5M_sun maintain radiative efficiencies as high as η~0.3-0.4 at high redshifts, near the maximum for rapidly spinning systems, but their efficiencies drop dramatically (by an order of magnitude) by z~0. The pattern for lower mass holes is somewhat more complicated but qualitatively similar. Assuming that the standard accretion disk model applies, we suggest that the accretion history of SMBHs and their accompanying spins evolve in two distinct regimes: an early phase of prolonged accretion, plausibly driven by major mergers, during which the black hole spins up, then switching to a period of random, episodic accretion, governed by minor mergers and internal secular processes, during which the hole spins down. The transition epoch depends on mass, mirroring other evidence for "cosmic downsizing" in the AGN population; it occurs at z~2 for high-mass black holes, and somewhat later, at z~1, for lower-mass systems.

preprint2012arXiv

Double-peaked narrow emission-line galaxies from Sloan Digital Sky Survey. I. Sample and basic properties

Recently, much attention has been given to double-peaked narrow emission-line galaxies, some of which are suggested to be related with merging galaxies. We make a systematic search to build the largest sample of these sources from Data Release 7 of Sloan Digital Sky Survey (SDSS DR7). With reasonable criteria of fluxes, full-width-half-maximum of emission lines and separations of the peaks, we select 3,030 double-peaked narrow emission-lines galaxies. In light of the existence of broad Balmer lines and the locations of the two components of double-peaked narrow emission lines distinguished by the Kauffmann et al. (2003) criteria in the Baldwin-Phillips-Terlevich (BPT) diagram, we find that there are 81 type I AGN, 837 double-type II AGN (2-type II), 708 galaxies with double star forming components (2-SF), 400 with mixed star forming and type II AGN components (type II + SF) and 1,004 unknown-type objects. As a by-product, a sample of galaxies (12,582) with asymmetric or top-flat profiles of emission lines is established. After inspecting the SDSS images of the two samples visually, we find 54 galaxies with dual cores. The present samples can be used to study the dynamics of merging galaxies, the triggering mechanism of black hole activity, the hierarchical growth of galaxies and dynamics of narrow line regions driven by outflows and rotating disk.

preprint2012arXiv

INTEGRAL and Swift observations of the Be X-ray binary 4U 1036-56 (RX J1037.5-5647) and its possible relation with gamma-ray transients

We present timing, spectral, and long-term temporal analysis of the high mass X-ray binary (HMXB) 4U 1036-56 using INTEGRAL and Swift observations. We show that it is a weak hard X-ray source spending a major fraction of the time in quiescence, and only occasionally characterized by X-ray outbursts. The outburst activity we report here lasts several days, with a dynamic range spanned by the luminosity in quiescence and in outburst as high as ~30. We report the detection of pulse period at 854.75+/-4.39 s during an outburst, which is consistent with previous measurements. Finally, we analyze the possibility of 4U 1036--56's association with the unidentified transient gamma-ray sources AGL J1037--5708 & GRO J1036--55, as prompted by its positional correlation.

preprint2012arXiv

Intervening BLR Clouds' Effects on Optical/UV Spectrum

Recent x-ray observations of Mrk 766 suggest that broad emission line region clouds cross our line of sight and produce variable x-ray absorption. Here we investigate what optical/ultraviolet spectroscopic features would be produced by such "Intervening BLR Clouds" (IBC) crossing our line of sight to the accretion disk, the source of the optical/UV continuum. Although the emission spectrum produced by intervening clouds is identical to the standard BLR model, they may produce absorption features on the optical or UV continuum. Single clouds will have little effect on the optical/UV spectrum because BLR clouds are likely to be much smaller than the accretion disk. This is unlike the X-ray case, where the radiation source is considerably smaller. However, an ensemble of intervening clouds will produce spectroscopic features in the FUV including a strong depression between the Lyman limit and Lyα. The amount of the depression will indicate the line-of-sight covering factor of clouds, an unknown quantity that is important for the ionization of the intergalactic medium and the energy budget of AGN. Comparison with observations suggests that the SED of Mrk 766 may be affected by intervening BLR clouds and IBC may exist in most of AGNs.

preprint2012arXiv

The CIV Baldwin effect in QSOs from Seventh Data Release of the Sloan Digital Sky Survey

Using the properties of SDSS DR7 QSOs catalog from Shen et al., the Baldwin effect, its slope evolution, the underlying drive for a large sample of 35019 QSOs with reliable spectral analysis are investigated. We find that the Baldwin effect exists in this large QSOs sample, which is almost the same in 11 different redshift bins, up to $z\sim 5$. The slope is -0.238 by the BCES (\civ\ EW depends on the continuum), -0.787 by the BCES bisector. For 11 redshift-bins, there is an increasing of the Baldwin effect slope from $z\sim1.5$ to $z\sim2.0$. From $z\sim2.0$ to $z\sim5.0$, the slope change is not clear considering their uncertainties or larger redshift bins. There is a strong correlation between the rest-frame \civ\ EW and \civ-based \mbh while the relation between the \civ\ EW and \mgii-based \mbh is very weak. With the correction of \civ-based \mbh from the \civ\ blueshift relative to \mgii, we suggest that this strong correlation is due to the bias of the \civ-based \mbh, with respect to that from the \mgii\ line. Considering the \mgii-based \mbh, a medium strong correlation is found between the \civ\ EW and the Eddington ratio, which implies that the Eddington ratio seems to be a better underlying physical parameter than the central black hole mass.

preprint2012arXiv

Two-Component Structure of the Hbeta Broad-Line Region in Quasars. I. Evidence from Spectral Principal Component Analysis

We report on a spectral principal component analysis (SPCA) of a sample of 816 quasars, selected to have small Fe II velocity shifts with spectral coverage in the rest wavelength range 3500--5500 Å. The sample is explicitly designed to mitigate spurious effects on SPCA induced by Fe II velocity shifts. We improve the algorithm of SPCA in the literature and introduce a new quantity, \emph{the fractional-contribution spectrum}, that effectively identifies the emission features encoded in each eigenspectrum. The first eigenspectrum clearly records the power-law continuum and very broad Balmer emission lines. Narrow emission lines dominate the second eigenspectrum. The third eigenspectrum represents the Fe II emission and a component of the Balmer lines with kinematically similar intermediate velocity widths. Correlations between the weights of the eigenspectra and parametric measurements of line strength and continuum slope confirm the above interpretation for the eigenspectra. Monte Carlo simulations demonstrate the validity of our method to recognize cross talk in SPCA and firmly rule out a single-component model for broad Hbeta. We also present the results of SPCA for four other samples that contain quasars in bins of larger Fe II velocity shift; similar eigenspectra are obtained. We propose that the Hbeta-emitting region has two kinematically distinct components: one with very large velocities whose strength correlates with the continuum shape, and another with more modest, intermediate velocities that is closely coupled to the gas that gives rise to Fe II emission.

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.

preprint2011arXiv

A Systematic Study on Energy Dependence of Quasi-Periodic Oscillation Frequency in GRS 1915+105

Systematically studying all the RXTE/PCA observations for GRS 1915+105 before November 2010, we have discovered three additional patterns in the relation between Quasi-Periodic Oscillation (QPO) frequency and photon energy, extending earlier outcomes reported by Qu et al. (2010). We have confirmed that as QPO frequency increases, the relation evolves from the negative correlation to positive one. The newly discovered patterns provide new constraints on the QPO models.

preprint2011arXiv

Cosmological evolution of supermassive black holes. I. mass function at 0<z<2

We present the mass function of supermassive black holes (SMBHs) over the redshift range z=0-2, using the latest deep luminosity and mass functions of field galaxies to constrain the masses of their spheroids, which we relate to SMBH mass through the empirical correlation between SMBH and spheroid mass (the M_bh-M_sph relation). In ddition to luminosity fading of the stellar content of the spheroids, we carefully consider the variation of the bulge-to-total luminosity ratio of the galaxy populations and the M_bh/M_sph ratio, which, according to numerous recent studies, evolves rapidly with redshift. The SMBH mass functions derived from the galaxy luminosity and mass functions show very good agreement, both in shape and in normalization. The resultant SMBH mass function and integrated mass density for the local epoch (z~0) match well those derived independently by other studies. Consistent with other evidence for cosmic downsizing, the upper end of the mass function remains roughly constant since $z\approx2$, while the space density of lower mass black holes undergoes strong evolution. We carefully assess the impact of various sources of uncertainties on our calculations. A companion paper uses the mass function derived in this work to determine the radiative efficiency of black hole accretion and constraints that can be imposed on the cosmological evolution of black hole spin.

preprint2011arXiv

The 2008 outburst of IGR J17473--2721: evidence for a disk corona?

The 2008 outburst of the atoll source IGR J17473--2721 was observed by INTEGRAL, RXTE and Swift. Tens of type-I X-ray bursts were found in this outburst. Joint observations provide sufficient data to look into the behavior of IGR J17473--2721 at the rising part of the outburst. We find that the joint energy spectra can be well fitted with a model composed of a blackbody and a cutoff power-law, with a cutoff energy decreasing from $ \sim$ 150 keV to $\sim$ 40 keV as the source leaves the quiescent state toward the low hard state. This fits into a scenario in which the corona is cooled by the soft X-rays along the outburst evolution, as observed in several other atoll sources. By using the flux measured in the 1.5--30 keV band of the type-I bursts during the outburst, we find that the linear relationship between the burst duration and the flux still holds for those bursts that occur at the decaying part of the low hard state, but with a different slope than the overall one that was estimated with the bursts happening in the whole extent of, and for the rest of the low hard state. The significance of such a dichotomy in the type-I X-ray bursts is $\sim$ 3 $σ$ under an F-test. Similar results are hinted at as well with the broader energy-band that was adopted recently. This dichotomy may be understood in a scenario where part of the accreting material forms a corona on the way of falling onto the surface of the neutron star during the decaying part of the low hard state. Based on the accretion rates of the preceding LHS, estimated from type-I X-ray bursts and from persistent emission, at least for IGR J17473-2721, most of the accretion material may fall on the neutron star (NS) surface in the LHS. Considering the burst behavior in the context of the outburst indicates a corona formed on top of the disk rather than on the NS surface.

preprint2010arXiv

Evolution of gaseous disk viscosity driven by supernova explosion. II. Structure and emissions from star-forming galaxies at high redshift

(Abridged) High redshift galaxies are undergoing intensive evolution of dynamical structure and morphologies. We incorporate the feedback into the dynamical equations through mass dropout and angular momentum transportation driven by the SNexp-excited turbulent viscosity. We numerically solve the equations and show that there can be intensive evolution of structure of the gaseous disk. Secular evolution of the disk shows interesting characteristics that are 1) high viscosity excited by SNexp can efficiently transport the gas from 10kpc to $\sim 1$kpc forming a stellar disk whereas a stellar ring forms for the case with low viscosity; 2) starbursts trigger SMBH activity with a lag $\sim 10^8$yr depending on star formation rates, prompting the joint evolution of SMBHs and bulges; 3) the velocity dispersion is as high as $\sim 100~\kms$ in the gaseous disk. In order to compare the present models with the observed dynamical structure and images, we use the incident continuum from the simple stellar synthesis (GALAXEV) and CLOUDY to calculate emission line ratios of H$α$, H$β$, $\OIII$ and $\NII$, and H$α$ brightness of gas photoionized by young massive stars formed on the disks. The models can produce the main features of emission from star forming galaxies and the observed relation between turbulent velocity and the H$α$ brightness. We successfully apply the present model to BX 389 and BX 482 observed in SINS high$-z$ sample, which are bulge and disk-dominated, respectively. High viscosity excited by SNexp is able to efficiently transport the gas into a bulge to maintain high star formation rates, or, to form a stellar ring close enough to the bulge so that it immigrates into the bulge of its host galaxy. This leads to a fast growing bulge. Implications and future work of the present models have been extensively discussed for galaxy formation.

preprint2010arXiv

INTEGRAL and Swift/XRT observations of the source PKS 0208-512

The active galaxy PKS 0208-512, detected at lower energies by COMPTEL, has been claimed to be a MeV blazar from EGRET. We report on the most recent INTEGRAL observations of the blazar PKS 0208-512, which are supplemented by Swift ToO observations. The high energy X-ray and gamma-ray emission of PKS 0208-512 during August - December 2008 has been studied using 682 ks of INTEGRAL guest observer time and ~ 56 ks of Swift/XRT observations. These data were collected during the decay of a gamma-ray flare observed by Fermi/LAT. At X-ray energies (0.2 - 10 keV) PKS 0208-512 is significantly detected by Swift/XRT, showing a power-law spectrum with a photon index of ~ 1.64. Its X-ray luminosity varied by roughly 30% during one month. At hard X-/soft gamma-ray energies PKS 0208-512 shows a marginally significant (~ 3.2 sigma) emission in the 0.5-1 MeV band when combining all INTEGRAL/SPI data. Non-detections at energies below and above this band by INTEGRAL/SPI may indicate intrinsic excess emission. If this possible excess is produced by the blazar, one possible explanation could be that its jet consists of an abundant electron-positron plasma, which may lead to the emission of an annihilation radiation feature. Assuming this scenario, we estimate physical parameters of the jet of PKS 0208-512.

preprint2010arXiv

Type-I bursts within outbursts of IGR J17473-2721

Two outbursts were observed by RXTE in the history of the atoll source IGR J17473-2721. During the most recent outburst in 2008, the source showed a complete series of spectral states/transitions. The neutron star system was prolific in type-I X-ray bursts, and we investigate them in the context of complete outbursts evolution. A total exposure of ~ 309 ks was collected by RXTE during the two outbursts of IGR J17473-2721. We carried out a systematic search for type-I bursts in this data set. For each burst found, we investigated the burst profile, the peak flux, and their dependence on the accretion rate along the evolution of the outbursts. Eighteen type-I X-ray bursts were found from IGR J17473-2721: two from the outburst in 2005 and the other 16 from the recent outburst in 2008. Among them, 3 bursts show photospheric radius expansion (PRE). The distance to the source is estimated as 6.4 kpc with a 15% uncertainty based on the three bursts that show PRE. In the recent outburst, there are 6 bursts showing up in the low/hard state prior to the state transition to a high/soft state, 3 bursts at the end phase of the high/soft state, and 7 in the following low/hard state. The blackbody radius of these bursts presents a variety of interesting features. We find that at the end of the recent outburst, the profile of the blackbody radius is anti-correlated with the blackbody temperature and the burst flux. The durations of the type-I burst are found to correlate with the Eddington ratio and to have two parallel evolution groups. Along the decreasing Eddington ratio, the burst duration decreases and ends in each group the PRE bursts occurred. This provides new clues to the type-I bursts in the context of outbursts for atoll XRBs.

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

Implications of infalling Fe II - emitting clouds in active galactic nuclei: anisotropic properties

We investigate consequences of the discovery that Fe II emission in quasars, one of the spectroscopic signatures of "Eigenvector 1", may originate in infalling clouds. Eigenvector 1 correlates with the Eddington ratio L/L_Edd so that Fe II/Hbeta increases as L/L_Edd increases. We show that the "force multiplier", the ratio of gas opacity to electron scattering opacity, is ~ 10^3 - 10^4 in Fe II-emitting gas. Such gas would be accelerated away from the central object if the radiation force is able to act on the entire cloud. As had previously been deduced, infall requires that the clouds have large column densities so that a substantial amount of shielded gas is present. The critical column density required for infall to occur depends on L/L_Edd, establishing a link between Eigenvector 1 and the Fe II/Hbeta ratio. We see predominantly the shielded face of the infalling clouds rather than the symmetric distribution of emitters that has been assumed. The Fe II spectrum emitted by the shielded face is in good agreement with observations thus solving several long-standing mysteries in quasar emission lines.

preprint2009arXiv

INTEGRAL and Swift/XRT observations of IGR J19405-3016

IGR J19405-3016 is reported in the 3rd IBIS catalog as being one of its lowest significance sources (~4.6 sigma under an exposure of about 371 ks). This leads to a caveat in multi-wavelength study although the source was identified in optical as a Seyfert 1. The currently increased INTEGRAL data stimulate us to investigate the reality of this source. We analyze all available observations carried out by INTEGRAL and Swift on IGR J19405-3016. We find that IGR J19405-3016 has a detection significance of ~ 9.4 sigma in the 20-60 keV band during the observational period between March 2003 and March 2008. Thus confirms a real source detection reported previously. The source position and error location are therefore updated. The source is found to be constant over years at the hard X-rays. Over the three XRT observations, the source flux varied by up to 39% from the average, and the spectrum is generally soft. The combined XRT/ISGRI spectrum is well fitted with a simple power law model (photon index 2.11+-0.03). Such a photon index is well consistent with the mean value 1.98 (dispersion 0.27) as obtained from Swift/BAT AGN samples at 14-195 keV. The spectral slope of IGR J19405-3016 is larger than the average spectral slope found by Molina et al. (2009). A similar discrepancy is found with the results of Beckmann et al. (2009) regarding Seyfert 1 AGNs. A possible explanation of this simple spectral description may be that the low level of the column density allows for the `true' spectrum to appear at soft X-rays as well.

preprint2009arXiv

Outbursts from IGR J17473-2721

We have investigated the outbursts of IGR J17473-2721. We analyzed all available observations carried out by RXTE on IGR J17473-2721 during its later outburst and as well all the available SWIFT/BAT data. The flux of the latter outburst rose in ~ one month and then kept roughly constant for the following ~ two months. During this time period, the source was in a low/hard state. The source moved to a high/soft state within the following three days, accompanied by the occurrence of an additional outburst at soft X-rays and the end of the preceding outburst in hard X-rays. During the decay of this soft outburst, the source went back to a low/hard state within 6 days, with a luminosity 4 times lower than the first transition. This shows a full cycle of the hysteresis in transition between the hard and the soft states. The fact that the flux remained roughly constant for ~ two months at times prior to the spectral transition to a high/soft state might be regarded as the result of balancing the evaporation of the inner disk and the inward accretion flow, in a model in which the state transition is determined by the mass flow rate. Such a balance might be broken via an additional mass flow accreting onto the inner disk, which lightens the extra soft outburst and causes the state transition. However, the possibility of an origin of the emission from the jet during this time period cannot be excluded. The spectral analysis suggests an inclined XRB system for IGR J17473-2721. Such a long-lived preceding low/hard state makes IGR J17473-2721 resemble the behavior of outbursts seen in black hole X-ray binaries like GX 339-4.

preprint2008arXiv

High energy properties of PKS 1830-211

We report on an analysis of X- and $γ$-ray observations of PKS 1830-211, based on the long-term campaigns carried out by \emph{INTEGRAL} and COMPTEL. The \emph{INTEGRAL} data currently available present a $33σ$ significance detection in the 20-100 keV band, while the COMPTEL 6-years data provide a $5.2σ$ significance detection in the 1-3 MeV energy band. At hard X-rays, \emph{INTEGRAL} and supplementary \emph{SWIFT} observations show flux variability on timescales of months. At $γ$-rays, the source shows persistent emission over years. The hard X-ray spectrum is well represented by a power-law model, with $Γ\sim 1.3$ in the 20-250 keV band. This photon index is well consistent with the previous report of $Γ\sim 1.3$ obtained at $E > 3.5$ keV from the best fit of \emph{XMM-Newton} data with a broken power law model. The joint \emph{XMM-Newton}/\emph{INTEGRAL} spectrum presented here is then fit with a broken power-law model and the parameters are refined compared to the previous. The results show the photon index changes from $\sim 1.0$ to $\sim 1.3$ at a break energy $\sim 4$ keV. At MeV energies, the spectrum softens to $Γ\sim 2.2$. These results, together with the EGRET measurement at $E \ge 100$ MeV, constitute a broad-band spectrum containing the peak of the power output at MeV energies, similar to most high-luminosity $γ$-ray blazars. The measured spectral characterstics are then discussed in the framework of the gravitational lens effects.

preprint2008arXiv

Integral and Rxte/Asm Observations on Igr J17098-3628

To probe further the possible nature of the unidentified source IGR J17098-3628, we have carried out a detailed analysis of its long-term time variability as monitored by RXTE/ASM, and of its hard X-ray properties as observed by INTEGRAL. INTEGRAL has monitored this sky region over years and significantly detected IGR J17098-3628 only when the source was in this dubbed active state. In particular, at $\ge$ 20 keV, IBIS/ISGRI caught an outburst in March 2005, lasting for $\sim$5 days with detection significance of 73$σ$ (20-40 keV) and with the emission at $< $200 keV. The ASM observations reveal that the soft X-ray lightcurve shows a similar outburst to that detected by INTEGRAL, however the peak of the soft X-ray lightcurve either lags, or is preceded by, the hard X-ray ($>$20 keV) outburst by $\sim$2 days. This resembles the behavior of X-ray novae like XN 1124-683, hence it further suggests a LMXB nature for IGR J17098-3628. While the quality of the ASM data prevents us from drawing any definite conclusions, these discoveries are important clues that, coupled with future observations, will help to resolve the as yet unknown nature of IGR J17098-3628.

preprint2008arXiv

Stochastic Electron Acceleration in Shell-Type Supernova Remnants

We study the stochastic electron acceleration by fast mode waves in the turbulent downstream of weakly magnetized collisionless astrophysical shocks. The acceleration is most efficient in a dissipative layer, and the model characteristics are determined by the shock speed, density, magnetic field, and turbulence decay length. The model explains observations of shell-type supernova remnants RX J1713.7-3946 and J0852.0-4622 and can be tested by observations in hard X-rays with the HXMT and NuSTAR or gamma-rays with the GLAST.

preprint2007arXiv

The Unified Model of Active Galactic Nuclei: II. Evolutionary Connection

(Abbreviated) We assemble a sample composed of 243 nearby Seyfert galaxies with redshifts $z\le 0.05$ to test the unification scheme. The sample contains 94 BLS1s, 44 NLS1s, 36 X-ray absorbed HBLR S2s, 42 X-ray absorbed non-HBLR S2s and 27 X-ray unabsorbed Seyfert 2s (unabsorbed non-HBLR S2s and HBLR S2s). We find that: 1) NLS1s have less massive black hole masses than BLS1s; 2) HBLRS2s have the same mass distribution of the black holes with BLS1s; 3) the absorbed non-HBLR S2s have less massive black holes than HBLR S2s; 4) unabsorbed non-HBLR S2s have the most massive black holes. We thus have a queue of black hole masses from small to large: narrow to broad line Seyfert galaxies, providing new evidence for the evolutionary sequence of Seyfert galaxies. We find that the opening angles of the torus in NLS1s and absorbed non-HBLR S2s are significantly smaller than that in BLS1s and HBLR S2s. The growth of the black holes and increases of the opening angles of the tori determine the various appearances of Seyfert galaxies. We also find that the unabsorbed Seyfert 2 galaxies could be caused by low gas-to-dust ratios in the present sample. This indicates that the star formation histories could be different in the unabsorbed from in absorbed Seyfert 2 galaxies, showing evidence for suppressed star formation by black hole activities. We outline a new unification scheme based on the orientation hypothesis: Seyfert galaxies can be unified by including growth of black holes, Eddington ratios, changing opening angles of tori and gas-to-dust ratios in the tori. Seyfert galaxies are tending to finally evolve to unabsorbed non-HBLR Seyfert 2 galaxies, in which the black holes are accreting with low accretion rates and both the broad line region and dusty torus disappear.