Researcher profile

Youjun Lu

Youjun Lu contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
16works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

16 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.