Source author record

Youjun Lu

Youjun Lu 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

44works
9topics
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

44 published item(s)

preprint2024arXiv

On using the counting method to constrain the anisotropy of kilonova radiation

A large number of binary neutron star (BNS) mergers are expected to be detected by gravitational wave (GW) detectors and the electromagnetic (EM) counterparts (e.g., kilonovae) of a fraction of these mergers may be detected in multi-bands by large area survey telescopes. For a given number of BNS mergers detected by their GW signals, the expected numbers of their EM counterparts that can be detected by a survey with given selection criteria depend on the kilonova properties, including the anisotropy. In this paper, we investigate whether the anisotropy of kilonova radiation and the kilonova model can be constrained statistically by the counting method, i.e., using the numbers of BNS mergers detected via GW and multi-band EM signals. Adopting simple models for the BNS mergers, afterglows, and a simple two (blue and red)-component model for kilonovae, we generate mock samples for GW detected BNS mergers, their associated kilonovae and afterglows detected in multi-bands. By assuming some criteria for searching the EM counterparts, we simulate the observations of these EM counterparts and obtain the EM observed samples in different bands. With the numbers of BNS mergers detected by GW detectors and EM survey telescopes in different bands, we show that the anisotropy of kilonova radiation and the kilonova model can be well constrained by using the Bayesian analysis. Our results suggest that the anisotropy of kilonova radiation may be demographically and globally constrained by simply using the detection numbers of BNS mergers by GW detectors and EM survey telescopes in multi-bands.

preprint2022arXiv

On Dark Gravitational Wave Standard Sirens as Cosmological Inference and Forecasting the Constraint on Hubble Constant using Binary Black Holes Detected by Deci-hertz Observatory

Gravitational wave (GW) signals from compact binary coalescences can be used as standard sirens to constrain cosmological parameters if their redshift can be measured independently. However, mergers of stellar binary black holes (BBHs) may not have electromagnetic counterparts and thus have no direct redshift measurements. These dark sirens may be still used to statistically constrain cosmological parameters by combining their GW measured luminosity distances and localization with deep redshift surveys of galaxies around it. We investigate this dark siren method in detail by using mock BBH and galaxy samples. We find that the Hubble constant can be constrained well with an accuracy $\lesssim1\%$ with a few tens or more of BBH mergers at redshift up to $1$ if GW observations can provide accurate estimates of their luminosity distance (with relative error of $\lesssim0.01$) and localization ($\lesssim0.1~\rm{deg}^2$), though the constraint may be significantly biased if the luminosity distance and localization errors are larger. We also introduce a simple method to correct this bias and find it is valid when the luminosity distance and localization errors are modestly large. We further generate mock BBH samples, according to current constraints on BBH merger rate and the distributions of BBH properties, and find that the Deci-hertz Observatory (DO) in a half year observation period may detect about one hundred BBHs with signal-to-noise ratio $\varrho\gtrsim30$, relative luminosity distance error $\lesssim0.02$, and localization error $\lesssim0.01\rm{deg}^2$. By applying the dark standard siren method, we find that the Hubble constant can be constrained to the $\sim0.1-1\%$ level using these DO BBHs, an accuracy comparable to the constraints obtained by using electromagnetic observations in the near future, thus it may provide insight into the Hubble tension.

preprint2022arXiv

On Detecting Stellar Binary Black Holes via the LISA-Taiji Network

The detection of gravitational waves (GWs) by ground-based laser interferometer GW observatories (LIGO/Virgo) reveals a population of stellar binary black holes (sBBHs) with (total) masses up to $\sim 150M_\odot$, which are potential sources for space-based GW detectors, such as LISA and Taiji. In this paper, we investigate in details on the possibility of detecting sBBHs by the LISA-Taiji network in future. We adopt the sBBH merger rate density constrained by LIGO/VIRGO observations to randomly generate mock sBBHs samples. Assuming an observation period of $4$ years, we find that the LISA-Taiji network may detect several tens (or at least several) sBBHs with signal-to-noise ratio (SNR) $>8$ (or $>15$), a factor $2-3$ times larger than that by only using LISA or Taiji observations. Among these sBBHs, no more than a few that can merge during the $4$-year observation period. If extending the observation period to $10$ years, then the LISA-Taiji network may detect about one hundred (or twenty) sBBHs with SNR $>8$ (or $>15$), among them about twenty (or at least several) can merge within the observation period. Our results suggest that the LISA-Taiji network may be able to detect at least a handful to twenty or more sBBHs even if assuming a conservative SNR threshold ($15$) for ``detection'', which enables multi-band GW observations by space and ground-based GW detectors. We also further estimate the uncertainties in the parameter estimations of the sBBH systems ``detected'' by the LISA-Taiji network. We find that the relative errors in the luminosity distance measurements and sky localization are mostly in the range of $0.05-0.2$ and $1-100°^2$, respectively, for these sBBHs.

preprint2022arXiv

Probing the nature of dark matter via gravitational waves lensed by small dark matter halos

Dark matter (DM) occupies the majority of matter content in the universe and is probably cold (CDM). However, modifications to the standard CDM model may be required by the small-scale observations, and DM may be self-interacting (SIDM) or warm (WDM). Here we show that the diffractive lensing of gravitational waves (GWs) from binary black hole mergers by small halos ($\sim10^3-10^6M_\odot$; mini-halos) may serve as a clean probe to the nature of DM, free from the contamination of baryonic processes in the DM studies based on dwarf/satellite galaxies. The expected lensed GW signals and event rates resulting from CDM, WDM, and SIDM models are significantly different from each other, because of the differences in halo density profiles and abundances predicted by these models. We estimate the detection rates of such lensed GW events for a number of current and future GW detectors, such as the Laser Interferometer Gravitational Observatories (LIGO), the Einstein Telescope (ET), the Cosmic Explorer (CE), Gravitational-wave Lunar Observatory for Cosmology (GLOC), the Deci-Hertz Interferometer Gravitational Wave Observatory (DECIGO), and the Big Bang Observer (BBO). We find that GLOC may detect one such events per year assuming the CDM model, DECIGO (BBO) may detect more than several (hundreds of) such events per year, by assuming the CDM, WDM (with mass $>30$\,keV) or SIDM model, suggesting that the DM nature may be strongly constrained by DECIGO and BBO via the detection of diffractive lensed GW events by mini-halos. Other GW detectors are unlikely to detect a significant number of such events within a limited observational time period. However, if the inner slope of the mini-halo density profile is sufficiently steeper than the Navarro-Frenk-White (NFW) profile, e.g., the pseudo-Jaffe profile, one may be able to detect one to more than hundred such GW events by ET and CE.

preprint2020arXiv

A Tidal Disruption Event Candidate Discovered in the Active Galactic Nucleus SDSS J022700.77-042020.6

We report the discovery of a Tidal Disruption Event (TDE) candidate occurring in the Active Galactic Nucleus SDSS J022700.77-042020.6. A sudden increase in flux of J0227-0420 during the second half of 2009 is clearly shown in the long-term optical, UV, and NIR light curves. A plateau phase, following an initial decline, is seen in the NUV and optical u, g, r, i light curves. Moreover, we find possible evidence that the plateau phase in the NUV band may lag behind the optical ones by approximately 70-80 days with also a much shorter duration, i.e. $\sim$7-15 days against $\sim$40-50 days. The long-term NUV/optical (after the plateau phase), NIR and MIR light curves can be well characterized with a form of $L(t)\propto t^{-β}$, consistent with the expectation of a TDE. The plateaus can be explained if the stellar streams collide with the pre-existing AGN disk at different radii. Though the overall fallback rate decreases, the material in the outer disk gradually drifts inward and increases the local accretion rate at the inner region, producing the optical and UV plateaus. The possible lag between the optical and NUV plateaus can then be attributed to viscosity delay. The index $β$ of the NIR $J, H, K_s$ bands ($\sim1.4-3.3$) is steeper than that of the UV/optical ($\sim0.7-1.3$) and MIR bands ($\sim0.9-1.8$), which may suggest that a certain fraction of the dust in the inner region of the dusty torus may be sublimated during the TDE phase. Our results indicate that, due to collisions between stellar debris and pre-existing disk, the light curves of TDEs occurring in AGN may show distinctive features, which may shed new light on the accretion process.

preprint2020arXiv

Broad-line region configuration of the supermassive binary black hole candidate PG1302-102 in the relativistic Doppler boosting scenario

PG1302-102 is thought to be a supermassive binary black hole (BBH) system according to the periodical variations of its optical and UV photometry, which may be interpreted as being due to the relativistic Doppler boosting of the emission mainly from the disk around the secondary black hole (BH) modulated by its orbital motion. In this paper, we investigate several broad emission lines of PG1302-102 using archived UV spectra obtained by IUE, GALEX, and Hubble, to reveal the broad-line region (BLR) emission properties of this BBH system under the Doppler boosting scenario. We find that the broad lines Ly$α$, NV, CIV, and CIII] all show Gaussian profiles, and none of these lines exhibits obvious periodical variation. Adopting a simple model for the BLR, we perform Markov chain Monte Carlo fittings to these broad lines, and find that the BLR must be viewed at an orientation angle of $\sim33^{\circ}$, close to face-on. If the Doppler boosting interpretation is correct, then the BLR is misaligned with the BBH orbital plane by an angle of $\sim51^\circ$, which suggests that the Doppler boosted continuum variation has little effect on the broad-line emission and thus does not lead to periodical line variation. We further discuss the possible implications for such a BLR configuration with respect to the BBH orbital plane.

preprint2020arXiv

Convergence and Efficiency of Different Methods to Compute the Diffraction Integral for Gravitational Lensing of Gravitational Waves

Wave optics may need to be considered when studying the lensed waveforms of gravitational waves (GWs). However, the computation of the diffraction integral (amplification factor) in wave optics is challenging and time-consuming. It is vital to develop an accurate and efficient method to calculate the amplification factor for detecting lensed GW systems. In this paper, we investigate the convergence of the diffraction integral for gravitational lensing of GWs and analyze the accuracy and efficiency of a number of numerical methods that can be used to calculate this integral, including the integral mean method, asymptotic expansion method, Levin's method, zero points integral method, etc. We further introduce a new method by combining the zero points integral and the asymptotic expansion methods to calculate the diffraction integral, which provides an efficient and accurate way to calculate the lensed waveform of GWs.

preprint2020arXiv

Dynamical evolution of cosmic supermassive binary black holes and their gravitational wave radiation

We investigate the evolution of supermassive binary black holes (BBHs) in galaxies with realistic property distributions and the gravitational-wave (GW) radiation from the cosmic population of these BBHs. We incorporate a comprehensive treatment of the dynamical interactions of the BBHs with their environments by including the effects of galaxy triaxial shapes and inner stellar distributions, and generate a large number of BBH evolution tracks. By combining these BBH evolution tracks, galaxy mass functions, galaxy merger rates, and supermassive black hole-host galaxy relations into our model, we obtain the statistical distributions of surviving BBHs, BBH coalescence rates, the strength of their GW radiation, and the stochastic GW background (GWB) contributed by the cosmic BBH population. About ~1%-3% (or ~10%) of supermassive BHs at nearby galactic centers are expected to be binaries with mass ratio >1/3 (or >1/100). The characteristic strain amplitude of the GWB at frequency 1/yr is estimated to be ~$2.0^{+1.4}_{-0.8}\times 10^{-16}$, and the upper bound of its results obtained with the different BH-host galaxy relations can be up to $5.4\times 10^{-16}$, which await testing by future experiments (e.g., the Square Kilometer Array, FAST, Next-Generation Very Large Array). The turnover frequency of the GWB spectrum is at ~0.25nHz. The uncertainties on the above estimates and prospects for detecting individual sources are also discussed. The application of the cosmic BBH population to the Laser Interferometer Space Antenna (LISA) band provides a lower limit to the detection rate of BBHs by LISA, ~0.9/yr.

preprint2020arXiv

Non-interacting black hole binaries with Gaia and LAMOST

Until recently, black holes (BHs) could be discovered only through accretion from other stars in X-ray binaries, or in merging double compact objects. Improvements in astrometric and spectroscopic measurements have made it possible to detect BHs also in non-interacting BH binaries (nBHB) through a precise analysis of the companion's motion. In this study, using an updated version of the Startrack binary-star population modelling code and a detailed model of the Milky Way (MW) galaxy we calculate the expected number of detections for Gaia and LAMOST surveys. We develop a formalism to convolve the binary population synthesis output with a realistic stellar density distribution, star-formation history (SFH), and chemical evolution for the MW, which produces a probability distribution function of the predicted compact-binary population over the MW. This avoids the additional statistical uncertainty which is introduced by methods which Monte Carlo sample from binary population synthesis output to produce one potential specific realisation of the MW compact-binary distribution, and our method is also comparatively fast to such Monte Carlo realisations. Specifically, we predict $\sim41$-$340$ nBHBs to be observed by Gaia, although the numbers may drop to $\sim10$-$70$ if the recent ($\lesssim100\;$ Myr) star formation is low ($\sim1\;M_\odot$/yr ). For LAMOST we predict $\lesssim14$ detectable nBHBs, which is lower partially because its field-of-view covers just $\sim6\%$ of the Galaxy.

preprint2020arXiv

On the Power to Constrain the Accretion History of Massive Black Holes via Spin Measurements by Upcoming X-Ray Telescopes

The spin distribution of massive black holes (MBHs) contains rich information on their assembly history. However, only limited information can be extracted from currently available spin measurements of MBHs owing to the small sample size and large measurement uncertainties. Upcoming X-ray telescopes with improved spectral resolution and larger effective area are expected to provide new insights into the growth history of MBHs. Here we investigate, at a proof of concept level, how stringent constraints can be placed on the accretion history of MBHs by the spin measurements from future X-ray missions. We assume a toy model consisting of a two-phase accretion history composed of an initial coherent phase with a constant disk orientation, followed by a chaotic phase with random disk orientations in each accretion episode. By utilizing mock spin data generated from such models and performing Bayesian Markov Chain Monte Carlo simulations, we find that most accretion models of MBHs can be reconstructed provided that $\gtrsim100$ MBH spins are measured with an accuracy of $\lesssim0.1$. We also quantify the precision of the reconstructed parameters by adopting various combinations of sample sizes and spin accuracies, and find that the sample size is more crucial to model reconstruction once the spin accuracy reaches $\sim 0.1$. To some extent, a better spin accuracy will compensate for a small sample size and vice versa. Future X-ray missions such as the Advanced Telescope for High Energy Astrophysics and the enhanced X-ray Timing and Polarimetry mission, may provide spin measurements of $\gtrsim100$ MBHs with an uncertainty of $\sim0.04-0.1$ and will thus put strong constraints on the MBH growth history.

preprint2020arXiv

Phase-dependent study of near-infrared disk emission lines in LB-1

The mass, origin and evolutionary stage of the binary system LB-1 has been the subject of intense debate, following the claim that it hosts an $\sim$70$M_{\odot}$ black hole, in stark contrast with the expectations for stellar remnants in the Milky Way. We conducted a high-resolution, phase-resolved spectroscopic study of the near-infrared Paschen lines in this system, using the 3.5-m telescope at Calar Alto Observatory. We find that Pa$β$ and Pa$γ$ (after proper subtraction of the stellar absorption component) are well fitted with a standard double-peaked model, typical of disk emission. We measured the velocity shifts of the red and blue peaks at 28 orbital phases: the line center has an orbital motion in perfect antiphase with the stellar motion, and the radial velocity amplitude ranges from 8 to 13 km/s for different choices of lines and profile modelling. We interpret this curve as proof that the disk is tracing the orbital motion of the primary, ruling out the circumbinary disk and the hierarchical triple scenarios. The phase-averaged peak-to-peak half-separation (proxy for the projected rotational velocity of the outer disk) is $\sim$70 km s$^{-1}$, larger than the stellar orbital velocity and also inconsistent with a circumbinary disk. From those results, we infer a primary mass 4--8 times higher than the secondary mass. Moreover, we show that the ratio of the blue and red peaks (V/R intensity ratio) has a sinusoidal behaviour in phase with the secondary star, which can be interpreted as the effect of external irradiation by the secondary star on the outer disk. Finally, we briefly discuss our findings in the context of alternative scenarios recently proposed for LB-1. Definitive tests between alternative solutions will require further astrometric data from $Gaia$.

preprint2020arXiv

Testing Ghasemi-Nodehi-Bambi metric parameters with quasi-periodic oscillations

QPOs are seen as peak features in the X-ray power spectral density of stellar mass black holes and neutron stars, of which frequencies can be measured with high precision. These QPO frequencies are believed to be related to fundamental frequencies of test particles, which are mainly determined by the background metric. We consider the metric introduced in Ghasemi-Nodehi and Bambi (Eur Phys J C 76:290, 2016). The fundamental frequencies in this metric are sensitive to some of the metric parameters but insensitive to other parameters, which means, the differences in fundamental frequencies in this metric and those in the Kerr ones can be significant for small changes of some but not all parameters around the Kerr value. By comparing with the QPO observations of GRO J1655-40, we find that only one parameter of the Ghasemi-Nodehi-Bambi metric can be strongly constrained, but other parameters cannot. We also use nested algorithm to investigate whether better constraints on the Ghasemi-Nodehi-Bambi metric parameters can be obtained from QPO observations of multiple objects by simulations. We find that four parameters can be strongly constrained while other parameters cannot . Our results suggest that QPOs may be important tools for testing the Kerr metric.

preprint2020arXiv

The Formation of a 70 Msun Black Hole at High Metallicity

A 70Msun BH was discovered in Milky Way disk in a long period and almost circular detached binary system (LB-1) with a high metallicity 8Msun B star companion. Current consensus on the formation of BHs from high metallicity stars limits the black hole mass to be below 20Msun. Using simple evolutionary model, we show that the formation of a 70Msun BH in high metallicity environment is possible if stellar wind mass loss rates are reduced by factor of 5. As observations indicate, a fraction of massive stars have surface magnetic fields which may quench the wind mass-loss, independently of stellar mass and metallicity. We also computed detailed stellar evolution models and we confirm such a scenario. A non-rotating 85Msun model at Z=0.014 with decreased winds ends up as a 71Msun star prior core-collapse with a 32Msun helium core and a 28Msun CO core. Such star avoids pair-instability pulsation supernova mass loss and may form a 70Msun BH in the direct collapse. Stars that can form such BHs expand to significant size with radius of R>600Rsun, exceeding the size of LB-1 orbit. Therefore, we can explain the formation of BHs upto 70Msun at high metallicity and this result is independent from LB-1. However, if LB-1 hosts a massive BH we are unable to explain how such a binary star system could have formed without invoking some exotic scenarios.

preprint2019arXiv

On using inspiralling supermassive binary black holes in the PTA frequency band as standard sirens to constrain dark energy

Supermassive binary black holes (SMBBHs) in galactic centers may radiate gravitational wave (GW) in the nano-Hertz frequency band, which are expected to be detected by pulsar timing arrays (PTAs) in the near future. GW signals from individual SMBBHs at cosmic distances, if detected by PTAs, are potentially powerful standard sirens that can be used to independently measure distances and thus put constraints on cosmological parameters. In this paper, we investigate the constraint that may be obtained on the equation of state ($w$) of dark energy by using those SMBBHs, expected to be detected by the PTAs in the Square Kilometre Array (SKA) era. By considering both the currently available SMBBH candidates and mock SMBBHs in the universe resulting from a simple galaxy major merger model, we find that $\sim 200$ to $3000$ SMBBHs with chirp mass $>10^9M_\odot$ are expected to be detected with signal-to-noise ratio $>10$ by SKA-PTA with conservative and optimistic settings and they can be used to put a constraint on $w$ to an uncertainty of $Δw\sim 0.02-0.1$. If further information on the mass and mass ratio of those SMBBHs can be provided by electromagnetic observations (e.g., chirp mass uncertainty $\lesssim 50\%$), the constraint may be further improved to $\lesssim 0.01$ level, as many more SMBBHs will be detected by SKA-PTA with relatively better distance measurements and can be used as the standard sirens.

preprint2019arXiv

Probing black hole accretion tracks, scaling relations and radiative efficiencies from stacked X-ray active galactic nuclei

The masses of supermassive black holes at the centres of local galaxies appear to be tightly correlated with the mass and velocity dispersions of their galactic hosts. However, the local Mbh-Mstar relation inferred from dynamically measured inactive black holes is up to an order-of-magnitude higher than some estimates from active black holes, and recent work suggests that this discrepancy arises from selection bias on the sample of dynamical black hole mass measurements. In this work we combine X-ray measurements of the mean black hole accretion luminosity as a function of stellar mass and redshift with empirical models of galaxy stellar mass growth, integrating over time to predict the evolving Mbh-Mstar relation. The implied relation is nearly independent of redshift, indicating that stellar and black hole masses grow, on average, at similar rates. Matching the de-biased local Mbh-Mstar relation requires a mean radiative efficiency ~0.15, in line with theoretical expectations for accretion onto spinning black holes. However, matching the "raw" observed relation for inactive black holes requires a mean radiative efficiency around 0.02, far below theoretical expectations. This result provides independent evidence for selection bias in dynamically estimated black hole masses, a conclusion that is robust to uncertainties in bolometric corrections, obscured active black hole fractions, and kinetic accretion efficiency. For our fiducial assumptions, they favour moderate-to-rapid spins of typical supermassive black holes, to achieve a mean radiative efficiency ~0.12-0.20. Our approach has similarities to the classic Soltan analysis, but by using galaxy-based data instead of integrated quantities we are able to focus on regimes where observational uncertainties are minimized.

preprint2019arXiv

Testing the relativistic Doppler boost hypothesis for supermassive binary black holes candidates via broad emission line profiles

Optical periodicity QSOs found by transient surveys are suggested to be sub-parsec supermassive binary black holes (BBHs). An intriguing interpretation for the periodicity of some of those QSOs is that the continuum is radiated from the accretion disk associated with the BBH secondary component and modulated by the periodical rotation of the secondary via Doppler-boost effect. Close to edge-on orbital orientation can lead to more significant Doppler-boost effect and thus are preferred for these systems, which is distinct from those normal type-1 QSOs with more or less face-on orientations. Therefore, the profiles of broad lines emitted from these Doppler-modulated systems may be significantly different from other systems that are not Doppler-modulated. We investigate the properties of the broad emission lines of optical-periodicity QSOs, including both a sample of QSOs that can be interpreted by the Doppler-modulated effects and a sample that cannot. We find that there is no obvious difference in the profiles and other properties of various (stacked) broad emission lines of these two samples, though a simple broad line region model would suggest significant differences. Our finding raises a challenge to the Doppler boost hypothesis for some of those BBHs candidates with optical periodicity.

preprint2016arXiv

A Baryonic Effect on the Merger Timescale of Galaxy Clusters

Accurate estimation of the merger timescale of galaxy clusters is important to understand the cluster merger process and further the formation and evolution of the large-scale structure of the universe. In this paper, we explore a baryonic effect on the merger timescale of galaxy clusters by using hydrodynamical simulations. We find that the baryons play an important role in accelerating the merger process. The merger timescale decreases with increasing the gas fraction of galaxy clusters. For example, the merger timescale is shortened by a factor of up to 3 for merging clusters with gas fractions 0.15, compared with the timescale obtained with zero gas fractions. The baryonic effect is significant for a wide range of merger parameters and especially more significant for nearly head-on mergers and high merging velocities. The baryonic effect on the merger timescale of galaxy clusters is expected to have impacts on the structure formation in the universe, such as the cluster mass function and massive substructures in galaxy clusters, and a bias of "no-gas" may exist in the results obtained from the dark matter-only cosmological simulations.

preprint2016arXiv

Dependence of the broad Fe K$α$ line on the physical parameters of AGN

In this paper, the dependence of the broad Fe K$α$ line on the physical parameters of AGN, such as the black hole mass $M_{BH}$, accretion rate (equivalently represented by Eddington ratio $λ_{Edd}$), and optical classification, is investigated by applying the X-ray spectra stacking method to a large sample of AGN which have well measured optical parameters. A broad line feature is detected ($>3σ$) in the stacked spectra of the high $λ_{Edd}$ sub-sample ($\logλ_{Edd}>-0.9$). The profile of the broad line can be well fitted with relativistic broad line model, with the line energy consistent with highly ionized Fe K$α$ line (i.e. Fe xxvi). A model consisting of multiple narrow lines cannot be ruled out, however. We found hints that the Fe K line becomes broader as the $λ_{Edd}$ increases. No broad line feature is shown in the sub-sample of broad-line Seyfert 1 (BLS1) galaxies and in the full sample, while a broad line might be present, though at low significance, in the sub-sample of narrow-line Seyfert 1 (NLS1) galaxies. We find no strong dependence of the broad line on black hole masses. Our results indicate that the detection/properties of the broad Fe K$α$ line may strongly depend on $λ_{Edd}$, which can be explained if the ionization state and/or truncation radius of the accretion disc changes with $λ_{Edd}$. The non-detection of the broad line in the BLS1 sub-sample can be explained if the the average EW of the relativistic Fe K$α$ line is weak or/and the fraction of sources with relativistic Fe K$α$ line is small in BLS1 galaxies.

preprint2016arXiv

Prospects for constraining the spin of the massive black hole at the Galactic center via the relativistic motion of a surrounding star

In this paper, we investigate the quality of constraining the spin of the massive black hole (MBH) at the Galactic center (GC) by using full general relativistic simulations of the motion of a surrounding star. We obtain the dependence mapping of the spin-induced signals on any spin direction of the MBH for given example stars, which indicates the feasibility to test whether the spin direction is the same as the normal of the young stellar disk located at the GC, and, further to provide insights into the assembly history of the MBH. We demonstrate the quality of constraining the MBH spin that may be achieved, given any set of the astrometric and the redshift precisions of observational facilities. We find that in the ranges of the astrometric and the velocity precisions with 1--30$μ$as and 0.1--10 km/s, an improvement in astrometric precision would be more effective at improving the quality of constraining the spin than an improvement in velocity precision. We obtain the parameter space of the semimajor axis and the eccentricity for the orbit of the target star that a high-precision constraint on the GC MBH spin can be obtained via the motion of the star. Our results show that the spin of the GC MBH can be constrained with 1-sigma error <~0.1 or even >~0.02 by monitoring the orbital motion of a star, if existing as expected, with semimajor axis <~300AU and eccentricity >~0.95 over a period shorter than a decade through future facilities.

preprint2015arXiv

A probable Milli-Parsec Supermassive Binary Black Hole in the Nearest Quasar Mrk 231

Supermassive binary black holes (BBHs) are unavoidable products of galaxy mergers and are expected to exist in the cores of many quasars. Great effort has been made during the past several decades to search for BBHs among quasars; however, observational evidence for BBHs remains elusive and ambiguous, which is difficult to reconcile with theoretical expectations. In this paper, we show that the distinct optical-to-UV spectrum of Mrk 231 can be well interpreted as emission from accretion flows onto a BBH, with a semimajor axis of ~590AU and an orbital period of ~1.2 year. The flat optical and UV continua are mainly emitted from a circumbinary disk and a mini-disk around the secondary black hole (BH), respectively; and the observed sharp drop off and flux deficit at wavelength lambda ~ 4000-2500 Angstrom is due to a gap (or hole) opened by the secondary BH migrating within the circumbinary disk. If confirmed by future observations, this BBH will provide a unique laboratory to study the interplay between BBHs and accretion flows onto them. Our result also demonstrates a new method to find sub-parsec scale BBHs by searching for deficits in the optical-to-UV continuum among the spectra of quasars.

preprint2015arXiv

On Testing the Kerr Metric of the Massive Black Hole in the Galactic Center via Stellar Orbital Motion: Full General Relativistic Treatment

The S-stars in the Galactic center (GC) are anticipated to provide unique dynamical constraint on the spin of the GC massive black hole (MBH). In this paper, we develop a fast full general relativistic method to simultaneously constrain the MBH mass, spin, and spin direction by considering both the motion of a star and the propagation of photons from the star to a distant observer. Assuming some example stars, we demonstrate that the spin-induced effects on the projected trajectory and redshift curve of a star depend on both the value and the direction of the spin. The maximum effects over a full orbit can differ by a factor upto more than one order of magnitude for cases with significantly different spin directions. Adopting the Markov Chain Monte Carlo fitting technique, we illustrate that the spin of the GC MBH is likely to be well constrained by using the motion of S0-2/S2 over a period of ~45yr if it is close to one and the astrometric and spectroscopic precisions (sigma_p,sigma_Z) can be as high as (10muas, 1km/s). In the mean time, the distance from the sun to the GC and the MBH mass can also be constrained to an unprecedented accuracy (0.01%-0.1%). If there exists a star with semimajor axis significantly smaller than that of S0-2/S2 and eccentricity larger than that of S0-2/S2, the MBH spin can be constrained with high accuracy over a period of <~10yr for (sigma_p,sigma_Z) ~ (10muas,1km/s), even if the spin is only moderately large (>~0.2).

preprint2015arXiv

Quiescent luminous red galaxies (LRGs) as cosmic chronometers: on the significance of mass and environmental dependence

Massive luminous red galaxies (LRGs) are believed to be evolving passively and can be used as cosmic chronometers to estimate the Hubble constant. However, different LRGs may be located in different environments. The environmental effects may limit the use of the LRGs as cosmic chronometers. We aim to investigate the environmental and mass dependence of the formation of 'quiescent' LRGs selected from the Sloan Digital Sky Survey Date Release 8 and to pave the way for using LRGs as cosmic chronometers. Using the population synthesis software STARLIGHT, we derive the stellar populations in each LRG through the full spectrum fitting and obtain the mean age distribution and the mean star formation history (SFH) of those LRGs. We find that there is no apparent dependence of the mean age and the SFH of quiescent LRGs on their environment, while the ages of those quiescent LRGs depend weakly on their mass. We compare the SFHs of the SDSS LRGs with those obtained from a semi-analytical galaxy formation model, and find that they are roughly consistent with each other if we consider the errors in the STARLIGHT-derived ages. We find that a small fraction of later star formation in LRGs leads to a systematical overestimation (~28 %) of the Hubble constant by the differential age method, and the systematical errors in the STARLIGHT-derived ages may lead to an underestimation (~ 16 %) of the Hubble constant. However, these errors can be corrected by a detailed study of the mean SFH of those LRGs and by calibrating the STARLIGHT-derived ages to those obtained independently by other methods. The environmental effects do not play significant role in the age estimates of quiescent LRGs, and the quiescent LRGs as a population can be used securely as cosmic chronometers.

preprint2015arXiv

Simulating the galaxy cluster "El Gordo" and identifying the merger configuration

The observational features of the massive galaxy cluster "El Gordo" (ACT-CL J0102-4915), such as the X-ray emission, the Sunyaev-Zel'dovich (SZ) effect, and the surface mass density distribution, indicate that they are caused by an exceptional ongoing high-speed collision of two galaxy clusters, similar to the well-known Bullet Cluster. We perform a series of hydrodynamical simulations to investigate the merging scenario and identify the initial conditions for the collision in ACT-CL J0102-4915. By surveying the parameter space of the various physical quantities that describe the two colliding clusters, including their total mass (M), mass ratio (ξ), gas fractions (f_b), initial relative velocity (V), and impact parameter (P), we find out an off-axis merger with P~800h_{70}^{-1}kpc, V~2500km/s, M~3x10^{15}Msun, and ξ=3.6 that can lead to most of the main observational features of ACT-CL J0102-4915. Those features include the morphology of the X-ray emission with a remarkable wake-like substructure trailing after the secondary cluster, the X-ray luminosity and the temperature distributions, and also the SZ temperature decrement. The initial relative velocity required for the merger is extremely high and rare compared to that inferred from currently available Lambda cold dark matter (LCDM) cosmological simulations, which raises a potential challenge to the LCDM model, in addition to the case of the Bullet Cluster.

preprint2014arXiv

Cosmological parameter fittings with the BICEP2 data

Combining the latest Planck, Wilkinson Microwave Anisotropy Probe (WMAP), and baryon acoustic oscillation (BAO) data, we exploit the recent cosmic microwave background (CMB) B-mode power spectra data released by the BICEP2 collaboration to constrain the cosmological parameters of the $Λ$CDM model, esp. the primordial power spectra parameters of the scalar and the tensor modes, $n_s,α_s, r, n_t$. We obtained constraints on the parameters for a lensed $Λ$CDM model using the Markov Chain Monte Carlo (MCMC) technique, the marginalized $68\%$ bounds are $r = 0.1043\pm_{0.0914}^{0.0307}$, $n_s = 0.9617\pm_{0.0061}^{ 0.0061}$, $α_s = -0.0175\pm_{0.0097}^{0.0105}$, $n_t = 0.5198\pm_{0.4579}^{0.4515}$. We found that a blue tilt for $n_t$ is favored slightly, but it is still well consistent with flat or even red tilt. Our $r$ value is slightly smaller than the one obtained by the BICEP team, as we allow $n_t$ as a free parameter without imposing the single-field slow roll inflation consistency relation. If we impose this relation, $r=0.2130\pm_{0.0609}^{0.0446}$. For most other parameters, the best fit values and measurement errors are not altered much by the introduction of the BICEP2 data.

preprint2014arXiv

Cyclic universe due to phantom and quintessence

We explore a cyclic universe due to phantom and quintessence fields. We find that, in every cycle of the evolution of the universe, the phantom dominates the cosmic early history and quintessence dominates the cosmic far future. In this model of universe, there are infinite cycles of expansion and contraction. Different from the inflationary universe, the corresponding cosmic space-time is geodesically complete and quantum stable. But similar to the Cyclic Model, the flatness problem, the horizon problem and the large scale structure of the universe can be explained in this cyclic universe.

preprint2014arXiv

Microlensing of Sub-parsec Massive Binary Black Holes in Lensed QSOs: Light Curves and Size-Wavelength Relation

Sub-parsec binary massive black holes (BBHs) are long anticipated to exist in many QSOs but remain observationally elusive. In this paper, we propose a novel method to probe sub-parsec BBHs through microlensing of lensed QSOs. If a QSO hosts a sub-parsec BBH in its center, it is expected that the BBH is surrounded by a circum-binary disk, each component of the BBH is surrounded by a small accretion disk, and a gap is opened by the secondary component in between the circum-binary disk and the two small disks. Assuming such a BBH structure, we generate mock microlensing light curves for some QSO systems that host BBHs with typical physical parameters. We show that microlensing light curves of a BBH QSO system at the infrared-optical-UV bands can be significantly different from those of corresponding QSO system with a single massive black hole (MBH), mainly because of the existence of the gap and the rotation of the BBH (and its associated small disks) around the center of mass. We estimate the half-light radii of the emission region at different wavelengths from mock light curves and find that the obtained half-light radius vs. wavelength relations of BBH QSO systems can be much flatter than those of single MBH QSO systems at a wavelength range determined by the BBH parameters, such as the total mass, mass ratio, separation, accretion rates, etc. The difference is primarily due to the existence of the gap. Such unique features on the light curves and half-light radius-wavelength relations of BBH QSO systems can be used to select and probe sub-parsec BBHs in a large number of lensed QSOs to be discovered by current and future surveys, including the Panoramic Survey Telescope and Rapid Response System (Pan-STARRS), the Large Synoptic Survey telescope (LSST) and Euclid.

preprint2014arXiv

Offsets between the X-ray and the Sunyaev-Zel'dovich-effect peaks in merging galaxy clusters and their cosmological implications

Observations reveal that the peaks of the X-ray map and the Sunyaev-Zel'dovich (SZ) effect map of some galaxy clusters are offset from each other. In this paper, we perform a set of hydrodynamical simulations of mergers of two galaxy clusters to investigate the spatial offset between the maxima of the X-ray and the SZ surface brightness of the merging clusters. We find that significantly large SZ-X-ray offsets (>100kpc) can be produced during the major mergers of galaxy clusters. The significantly large offsets are mainly caused by a `jump effect' occurred between the primary and secondary pericentric passages of the two merging clusters, during which the X-ray peak may jump to the densest gas region located near the center of the small cluster, but the SZ peak remains near the center of the large one. Our simulations show that merging systems with higher masses and larger initial relative velocities may result in larger offset sizes and longer offset time durations; and only nearly head-on mergers are likely to produce significantly large offsets. We further investigate the statistical distribution of the SZ-X-ray offset sizes and find that (1) the number distribution of the offset sizes is bimodal with one peak located at low offsets ~0 and the other at large offsets ~350-450kpc/h, but the objects with intermediate offsets are scarce; and (2) the probabilities of the clusters in the mass range higher than 2x10^{14}Msun/h that have offsets larger than 20, 50, 200, 300, and 500kpc/h are 34.0%, 11.1%, 8.0%, 6.5%, and 2.0% respectively at z=0.7. The probability is sensitive to the underlying pairwise velocity distribution and the merger rate of clusters. Future observations on the offsets for a large number of clusters may put strong constraints on the cosmic velocity fields on the cluster scale and the cluster merger rate. (Abridged)

preprint2014arXiv

On the Existence of Pulsars in the Vicinity of the Massive Black Hole in the Galactic Center

Pulsars, if existing and detectable in the immediate vicinity of the massive black hole (MBH) in the Galactic center (GC), may be used as a superb tool to probe both the environment and the metric of the central MBH. The recent discovery of a magnetized pulsar in the GC suggests that many more pulsars should exist near the MBH. In this paper, we estimate the number and the orbital distribution of pulsars in the vicinity of the MBH in the GC by assuming that the pulsar progenitors, similar to the GC S-stars, were captured to orbits tightly bound to the MBH through the tidal breakup of stellar binaries. We use the current observations on both the GC S-stars and the hypervelocity stars to calibrate the injection rate(s) of and the dynamical model(s) for the stellar binaries. By including the relaxation processes, supernova kicks, and gravitational wave radiation in our simulations, we estimate that ~97-190 (9-14) pulsars may presently orbit the central MBH with semimajor axes <=4000AU (<=1000AU), which is compatible with the current observational constraints on the number of the GC pulsars. The semimajor axis and the pericenter distance of the pulsar closest to the central MBH are probably in the range of ~120-460AU and ~2-230AU, respectively. Future telescopes, such as the SKA, may be able to detect a significant number of pulsars with semimajor axis smaller than a few thousand AU in the GC. Long-term monitoring of these pulsars would be helpful in constraining both the environment and the metric of the central MBH. Our preferred model also results in about ten hyperfast pulsars with velocity >~1500km/s moving away from the Milky Way.

preprint2014arXiv

Probing baryonic processes and gastrophysics in the formation of the Milky Way dwarf satellites: I. metallicity distribution properties

In this paper, we study the chemical properties of the stars in the dwarf satellites around the MW-like host galaxies, and explore the possible effects of several baryonic processes, including supernova (SN) feedback, the reionization of the universe and H$_2$ cooling, on them and how current and future observations may put some constraints on these processes. We use a semi-analytical model to generate MW-like galaxies, for which a fiducial model can reproduce the luminosity function and the stellar metallicity--stellar mass correlation of the MW dwarfs. Using the simulated MW-like galaxies, we focus on investigating three metallicity properties of their dwarfs: the stellar metallicity--stellar mass correlation of the dwarf population, and the metal-poor and metal-rich tails of the stellar metallicity distribution in individual dwarfs. We find that (1) the slope of the stellar metallicity--stellar mass correlation is sensitive to the SN feedback strength and the reionization epoch; (2) the extension of the metal-rich tails is mainly sensitive to the SN feedback strength; (3) the extension of the metal-poor tails is mainly sensitive to the reionization epoch; (4) none of the three chemical properties are sensitive to the H$_2$ cooling process; and (5) comparison of our model results with the current observational slope of the stellar metallicity--stellar mass relation suggests that the local universe is reionized earlier than the cosmic average and local sources may have a significant contribution to the reionization in the local region, and an intermediate to strong SN feedback strength is preferred. Future observations of metal-rich and metal-poor tails of stellar metallicity distributions will put further constraints on the SN feedback and the reionization processes.

preprint2014arXiv

Relativistic Fe K$α$ Line Revealed in the Composite X-ray Spectrum of Narrow Line Seyfert 1 Galaxies --- do their black holes have averagely low or intermediate spins?

While a broad profile of the Fe K$α$ emission line is frequently found in the X-ray spectra of typical Seyfert galaxies, the situation is unclear in the case of Narrow Line Seyfert 1 galaxies (NLS1s)---an extreme subset which are generally thought to harbor less massive black holes with higher accretion rates. In this paper, the ensemble property of the Fe K$α$ line in NLS1s is investigated by stacking the X-ray spectra of a large sample of 51 NLS1s observed with {\it XMM-Newton}. The composite X-ray spectrum reveals a prominent, broad emission feature over 4--7 keV, characteristic of the broad Fe K$α$ line. In addition, there is an indication for a possible superimposing narrow (unresolved) line, either emission or absorption, corresponding to Fe XXVI or Fe XXV, respectively. The profile of the broad emission feature can well be fitted with relativistic broad-line models, with the line energy consistent either with 6.4 keV (i.e., neutral Fe) or with 6.67 keV (i.e., highly ionized Fe), in the case of the narrow line being emission and absorption, respectively. Interestingly, there are tentative indications for low or intermediate values of the average spins of the black holes ($a<0.84$), as inferred from the profile of the composite broad line. If the observed feature is indeed a broad line rather than resulting from partial covering absorption, our results suggest that a relativistic Fe line may in fact be common in NLS1s; and there are tentative indications that black holes in NLS1s may not spin very fast in general.

preprint2013arXiv

On the Full Spectrum Fitting of Luminous Red Galaxies by Using ULySS and STARLIGHT

In this paper, we select a sample of quiescent luminous red galaxies (LRGs) from Sloan Digital Sky Survey Data Release 7 (SDSS DR7) with signal-to-noise ratio (S/N) substantially high to study the consistency of the full spectrum fitting method by using different packages, mainly, ULySS and STARLIGHT. The spectrum of each galaxy in the sample is fitted by the full spectrum fitting packages ULySS and STARLIGHT, respectively. We find: (1) for spectra with higher S/Ns, the ages of stellar populations obtained from ULySS are slightly older than that from STARLIGHT, and metallicities derived from ULySS are slightly richer than that from STARLIGHT. In general, both packages can give roughly consistent fitting results. (2) for low S/N spectra, it is possible that the fitting by ULySS can be trapped some local minimum parametric regions during execution and thus may give unreliable results, while STARLIGHT can still give reliable results. Based on the fitting results of LRGs, we further analyse their star formation history (SFH) and the relation between their age and velocity dispersion, and find that those highly agrees with conclusions from others' previous work.

preprint2013arXiv

Radiative Efficiency of Disk Accretion in Individual SDSS QSOs

We estimate the radiative efficiency (epsilon) of individual type 1 SDSS QSOs by using their bolometric luminosities (L_{bol}) and accretion rates (\dot{M}_{\bullet,acc}), which may be related to the assembly histories and spins of the central massive black holes (MBHs). We estimate L_{bol} by using the empirical spectral energy distributions of QSOs and \dot{M}_{\bullet,acc} by fitting the observed optical luminosity(/-ies) with the thin accretion disk model, assuming the MBH masses given by the virial mass estimator(s) (M_{\bullet,vir}). We find an apparent correlation between epsilon and M_{\bullet,vir}, which is strong at redshift z \la 1.8, weak at z \ga 2, and consistent with that found by Davis & Laor (2011) for 80 PG QSOs at z \leq 0.5. To investigate whether this correlation is intrinsic or not, we construct a mock sample of QSOs according to the true MBH mass and Eddington ratio distributions given in Kelly & Shen (2013). By comparing the results obtained from the mock sample with that from the SDSS sample, we demonstrate that the apparent epsilon-M_{\bullet,vir} correlation can be produced by and mainly due to the selection effects of the SDSS sample and the bias induced by the usage of M_{\bullet,vir} as the true MBH mass. The mean values of epsilon of those SDSS QSOs are consistent with being a constant \simeq 0.11-0.16 over the redshift range of 0.3 \la z \la 4. We conclude that the current SDSS QSO data is consistent with no strong intrinsic correlation between radiative efficiency and true MBH mass and no significant redshift evolution of radiative efficiencies.

preprint2013arXiv

S-stars in the Galactic center and hypervelocity stars in the Galactic halo: two faces of the tidal breakup of stellar binaries by the central massive black hole?

In this paper, we investigate the link between the hypervelocity stars (HVSs) discovered in the Galactic halo and the S-stars moving in the Galactic center (GC), under the hypothesis that they are both the products of the tidal breakup of the same population of stellar binaries by the central massive black hole (MBH). By adopting several hypothetical models for binaries to be injected into the vicinity of the MBH and doing numerical simulations, we realize the tidal breakup processes of the binaries and their follow-up evolution. We find that many statistical properties of the detected HVSs and S-stars can be reproduced under some binary injecting models, and their number ratio can be reproduced if the stellar initial mass function is top-heavy (e.g., with slope ~-1.6). The total number of the captured companions is ~50 that have masses in the range ~3-7Msun and semimajor axes <~4000 AU and survive to the present within their main-sequence lifetime. The innermost one is expected to have a semimajor axis ~300-1500 AU and a pericenter distance ~10-200 AU, with a significant probability of being closer to the MBH than S2. Future detection of such a closer star would offer an important test to general relativity. The majority of the surviving ejected companions of the S-stars are expected to be located at Galactocentric distances <~20 kpc, and have heliocentric radial velocities ~-500-1500 km/s and proper motions up to ~5-20 mas/yr. Future detection of these HVSs may provide evidence for the tidal-breakup formation mechanism of the S-stars.

preprint2012arXiv

The Age-Redshift Relation For Luminous Red Galaxies Obtained From the Full Spectrum Fitting and Its Cosmological Implications

The relative age of galaxies at different redshifts can be used to infer the Hubble parameter and put constraints on cosmological models. We select 23,883 quiescent luminous red galaxies (LRGs) from the SDSS DR7 and divide them into four sub-samples according to their velocity dispersions and each sub-sample is further divided into 12 redshift bins. The spectra of the LRGs in each redshift and velocity bin are co-added in order to obtain a combined spectrum with relatively high $S/N$. Adopting the GalexEV/SteLib model, we estimate the mean ages of the LRGs from these combined spectra by the full-spectrum fitting method. We check the reliability of the estimated age by using Monte-Carlo simulations and find that the estimates are robust and reliable. Assuming that the LRGs in each sub-sample and each redshift bin were on average formed at the same time, the Hubble parameter at the present time $H_0$ is estimated from the age--redshift relation obtained for each sub-sample, which is compatible with the $H_0$ value measured by other methods. We demonstrate that a systematic bias (up to $\sim 20%$) may be introduced to the $H_0$ estimation because of recent star formation in the LRGs due to the later major mergers at $z\la 0.4$, but this bias may be negligible for those sub-samples with large velocity dispersions. Using the age--redshift relations obtained from the sub-sample with the largest velocity dispersion or the two sub-samples with high velocity dispersions, we find $H_0= 65^{+7}_{-3}\kmsmpc$ or $H_0= 74^{+5}_{-4}\kmsmpc$ by assuming a spatially flat $Λ$CDM cosmology. With upcoming surveys, such as the Baryon Oscillation Spectroscopic Survey (BOSS), even larger samples of quiescent massive LRGs may be obtained, and thus the Hubble parameter can be measured with high accuracy through the age--redshift relation.

preprint2012arXiv

The cosmic evolution of massive black holes and galaxy spheroids: global constraints at redshift z<~1.2

We study the observational constraints on the cosmic evolution of the relationships between the massive black hole (MBH) mass (M_bh) and the stellar mass (M^*_sph; or velocity dispersion σ) of the host galaxy/spheroid. Assuming that the M_bh-M^*_sph (or M_bh-σ) relation evolves with redshift as \propto (1+z)^Γ, the MBH mass density can be obtained from either the observationally determined galaxy stellar mass functions or velocity dispersion distribution functions over redshift z~0-1.2 for any given Γ. The MBH mass density at different redshifts can also be inferred from the luminosity function of QSOs/AGNs provided known radiative efficiency ε. By matching the MBH density inferred from galaxies to that obtained from QSOs/AGNs, we find that Γ=0.64^{+0.27}_{-0.29} for the M_bh-M^*_sph relation and Γ=-0.21^{+0.28}_{-0.33} for the M_bh-σrelation, and ε=0.11^{+0.04}_{-0.03}. Our results suggest that the MBH mass growth precedes the bulge mass growth but the galaxy velocity dispersion does not increase with the mass growth of the bulge after the quench of nuclear activity, which is roughly consistent with the two-phase galaxy formation scenario proposed by Oser et al. (2012) in which a galaxy roughly double its masses after z=1 due to accretion and minor mergers while its velocity dispersion drops slightly.

preprint2011arXiv

Metal-poor hypervelocity star candidates from the Sloan Digital Sky Survey

Hypervelocity stars are believed to be ejected out from the Galactic center through dynamical interactions of (binary) stars with the central massive black hole(s). In this letter, we report 13 metal-poor F-type hypervelocity star candidates selected from 370,000 stars of the data release 7 of the Sloan Digital Sky Survey. With a detailed analysis of the kinematics of these stars, we find that seven of them were likely ejected from the Galactic center (GC) or the Galactic disk, four neither originated from the GC nor the Galactic disk, and the other two were possibly ejected from either the Galactic disk or other regions. Those candidates which unlikely originated from the GC or the Galactic disk, may be explained by other mechanisms, like the tidal disruption of the Milky Way's dwarf galaxies in the Galactic potential, or the gravitational interactions with a massive black hole at the center of M31 or M32.

preprint2011arXiv

On probing the properties of QSOs through their proximity effects on the intergalactic medium

The proximity effect (PE) of QSOs is believed to be useful in constraining the QSO lifetime. Observations on the PE so far, however, give apparently contradictory results -- some are consistent with a long QSO lifetime (>~ a few 10^7 yr), but others appear to be only consistent with a short QSO lifetime <~ 10^6 yr. In this paper, we show that this apparent contradiction may be solved by simultaneously taking into account both the effect due to the density enhancement in the QSO near zones and that due to the obscuration of the tori associated with the QSOs, using a large number of Monte-Carlo generated synthetic Lyman alpha forest spectra. We demonstrate that the QSO properties and environment can be constrained simultaneously by the transverse PE and the line of sight PE of bright type 1 QSOs together. The current available measurements on the PEs of type 1 QSOs suggest that (1) the density is significantly enhanced in the vicinity of the QSOs; (2) the QSO lifetime is consistent with being as large as a few 10^7 yr and a substantially shorter lifetime (e.g., <~10^6 yr) is not required; and (3) the half opening angle of the tori associated with QSOs is ~60 deg, consistent with some other independent estimates. Our simulations also show that the TPE of type 2 QSOs can be significantly different from that of type 1 QSOs, which may be useful to put further constraints on the QSO properties and the QSO environment.

preprint2011arXiv

The low frequency of dual AGNs versus the high merger rate of galaxies: A phenomenological model

Dual AGNs are natural byproducts of hierarchical mergers of galaxies in the LambdaCDM cosmogony. Recent observations have shown that only a small fraction (~ 0.1%-1%) of AGNs at redshift z<~ 0.3 are dual with kpc-scale separations, which is rather low compared to the high merger rate of galaxies. Here we construct a phenomenological model to estimate the number density of dual AGNs and its evolution according to the observationally-estimated major merger rates of galaxies and various scaling relations on the properties of galaxies and their central massive black holes. We show that our model reproduces the observed frequency and separation distribution of dual AGNs provided that significant nuclear activities are triggered only in gas-rich progenitor galaxies with central massive black holes and only when the nuclei of these galaxies are roughly within the half-light radii of their companion galaxies. Under these constraints, the observed low dual AGN frequency is consistent with the relatively high merger rate of galaxies and supports the hypothesis that major mergers lead to AGN/QSO activities. We also predict that the number of kpc-scale dual AGNs decreases with increasing redshift and only about 0.02%--0.06% of AGNs are dual AGNs with double-peaked narrow line features at redshifts of z 0.5-1.2. Future observations of high-redshift dual AGNs would provide a solid test for this prediction.

preprint2011arXiv

The spatial and velocity distributions of hypervelocity stars

Hypervelocity stars (HVSs) found in the Galactic halo are probably the dynamical products of interactions between (binary) stars and the massive black hole(s) (MBH) in the Galactic center (GC). It has been shown that the detected HVSs are spatially consistent with being located on two thin disks (Lu et al.), one of which has the same orientation as the clockwise-rotating stellar disk in the GC. Here we perform a large number of three-body experiments of the interactions between the MBH and binary stars bound to it, and find that the probability of ejecting HVSs is substantially enhanced by multiple encounters between the MBH and binary stars at distances substantially larger than the initial tidal breakup radii. Assuming that the HVS progenitors are originated from the two disks, the inclination distribution of the HVSs relative to the disk planes can be reproduced by either the mechanism of tidal breakup of binary stars or the mechanism of ejecting HVSs by a hypothetical binary black hole (BBH) in the GC. However, an isotropical origination of HVS progenitors is inconsistent with the observed inclination distribution. Assuming that the HVSs were ejected out by the tidal breakup mechanism, its velocity distribution can be reproduced if their progenitors diffuse onto low angular momentum orbits slowly and most of the progenitors were broken up at relatively large distances due to multiple encounters. Assuming that the HVSs were ejected out by a BBH within the allowed parameter space in the GC, our simulations produce relatively flatter velocity spectra compared to the observed ones; however, the BBH mechanism cannot be statistically ruled out, yet. Future surveys of HVSs and better statistics of their spatial and velocity distributions should enable to distinguish the ejection mechanisms of HVSs and shed new light on the dynamical environment of the MBH.(abridged)

preprint2009arXiv

On the spatial distribution and the origin of hypervelocity stars

Hypervelocity stars (HVSs) escaping away from the Galactic halo are dynamical products of interactions of stars with the massive black hole(s) (MBH) in the Galactic Center (GC). They are mainly B-type stars with their progenitors unknown. OB stars are also populated in the GC, with many being hosted in a clockwise-rotating young stellar (CWS) disk within half a parsec from the MBH and their formation remaining puzzles. In this paper, we demonstrate that HVSs can well memorize the injecting directions of their progenitors using both analytical arguments and numerical simulations, i.e., the ejecting direction of an HVS is almost anti-parallel to the injecting direction of its progenitor. Therefore, the spatial distribution of HVSs maps the spatial distribution of the parent population of their progenitors directly. We also find that almost all the discovered HVSs are spatially consistent with being located on two thin disk planes. The orientation of one plane is consistent with that of the (inner) CWS disk, which suggests that most of the HVSs originate from the CWS disk or a previously existed disk-like stellar structure with an orientation similar to it. The rest of HVSs may be correlated with the plane of the northern arm of the mini-spiral in the GC or the plane defined by the outer warped part of the CWS disk. Our results not only support the GC origin of HVSs but also imply that the central disk (or the disk structure with a similar orientation) should persist or be frequently rejuvenated over the past 200 Myr, which adds a new challenge to the stellar disk formation and provides insights to the longstanding problem of gas fueling into massive black holes.

preprint2008arXiv

Toward precise constraints on growth of massive black holes

Growth of massive black holes (MBHs) in galactic centers comes mainly from gas accretion during their QSO/AGN phases. In this paper we apply an extended Soltan argument, connecting the local MBH mass function with the time-integral of the QSO luminosity function, to the demography of MBHs and QSOs from recent optical and X-ray surveys, and obtain robust constraints on the luminosity evolution (or mass growth history) of individual QSOs (or MBHs). We find that the luminosity evolution probably involves two phases: an initial exponentially increasing phase set by the Eddington limit and a following phase in which the luminosity declines with time as a power law (with a slope of -1.2--1.3) set by a self-similar long-term evolution of disk accretion. Neither an evolution involving only the increasing phase with a single Eddington ratio nor an exponentially declining pattern in the second phase is likely. The period of a QSO radiating at a luminosity higher than 10% of its peak value is about (2-3)x10^8 yr, during which the MBH obtains ~80% of its mass. The mass-to-energy conversion efficiency is $0.16\pm0.04 ^{+0.05}_{-0}$, with the latter error accounting for the maximum uncertainty due to Compton-thick AGNs. The expected Eddington ratios in QSOs from the constrained luminosity evolution cluster around a single value close to 0.5-1 for high-luminosity QSOs and extend to a wide range of lower values for low-luminosity ones. The Eddington ratios for high luminosity QSOs appear to conflict with those estimated from observations (~0.25) by using some virial mass estimators for MBHs in QSOs unless the estimators systematically over-estimate MBH masses by a factor of 2-4. We also infer the fraction of optically obscured QSOs ~60-80%. Further applications of the luminosity evolution of individual QSOs are also discussed.

preprint2007arXiv

Hypervelocity binary stars: smoking gun of massive binary black holes

The hypervelocity stars recently found in the Galactic halo are expelled from the Galactic center through interactions between binary stars and the central massive black hole or between single stars and a hypothetical massive binary black hole. In this paper, we demonstrate that binary stars can be ejected out of the Galactic center with velocities up to 10^3 km/s, while preserving their integrity, through interactions with a massive binary black hole. Binary stars are unlikely to attain such high velocities via scattering by a single massive black hole or through any other mechanisms. Based on the above theoretical prediction, we propose a search for binary systems among the hypervelocity stars. Discovery of hypervelocity binary stars, even one, is a definitive evidence of the existence of a massive binary black hole in the Galactic center.

preprint2007arXiv

On the origin of kinematic distribution of the sub-parsec young stars in the Galactic center

Within a half-parsec from the Galactic center (GC), there is a population of coeval young stars which appear to reside in a coherent disk. Surrounding this dynamically-cool stellar system, there is a population of stars with a similar age and much larger eccentricities and inclinations relative to the disk. We propose a hypothesis for the origin of this dynamical dichotomy. Without specifying any specific mechanism, we consider the possibility that both stellar populations were formed within a disk some 6 Myr ago. But this orderly structure was dynamically perturbed outside-in by an intruding object with a mass ~10^4 Msun, which may be an intermediate-mass black hole (IMBH) or a dark stellar cluster hosting an IMBH. We suggest that the perturber migrated inward to ~0.15-0.3pc from the GC under the action of dynamical friction. Along the way, it captured many stars in the outer disk region into its mean-motion resonance, forced them to migrate with it, closely encountered with them, and induced the growth of their eccentricity and inclination. But stars in the inner regions of the disk retain their initial coplanar structure. We predict that some of the inclined and eccentric stars surrounding the disk may have similar Galactocentric semimajor axis. Future precision determination of their kinematic distribution of these stars will not only provide a test for this hypothesis but also evidences for the presence of an IMBH or a dark cluster at the immediate proximity of the massive black hole at the GC. (abridged)

preprint2002arXiv

The relativistic Iron K-alpha line from an accretion disc onto a static non-baryonic compact object

This paper continues the study of the properties of an accretion disc rotating around a non-baryonic (assumed super-massive) compact object. This kind of objects, generically known as boson stars, were earlier proposed as a possible alternative scenario to the existence of super-masive black holes in the center of every galaxy. A dilute boson star has also been proposed as a large part of the non-baryonic dark matter, flattening galactic rotational velocities curves. In this contribution, we compute the profile of the emission lines of Iron; its shape has been for long known as a useful diagnosis of the space-time geometry. We compare with the case of a Schwarzschild black hole, concluding that the differences are observationally distinguishable.