Showing 1–50 of 59 results for author: Yi, T

Optical intraday variability analysis for the BL Lacertae object 1ES 1426+42.8

Authors: X. Chang, D. R. Xiong, T. F. Yi, C. X. Liu, G. Bhatta, J. R. Xu, Y. L. Gong

Abstract: The observation data of blazar 1ES 1426 + 42.8 were obtained using the 1.02 m optical telescope of Yunnan Observatories during $2021$ to $2023$. Intraday variability (IDV) is detected on seven nights. We use the turbulent model to investigate the mechanism of IDV in 1ES 1426 + 42.8. The fitting light curves match the actual IDV curves well. Using this model, we obtain the parameters such as the si… ▽ More The observation data of blazar 1ES 1426 + 42.8 were obtained using the 1.02 m optical telescope of Yunnan Observatories during $2021$ to $2023$. Intraday variability (IDV) is detected on seven nights. We use the turbulent model to investigate the mechanism of IDV in 1ES 1426 + 42.8. The fitting light curves match the actual IDV curves well. Using this model, we obtain the parameters such as the size of turbulent cells and the width of pulses in the jet. A possible short-lived Quasi-periodic oscillation (QPO) of $58.55 \pm 8.09$ minutes was detected on April 26, 2022 whose light curve exhibits eight cycles at $>3σ$ global significance and confirmed by several different techniques. Through a more detailed analysis of the light curve of this night, we find that the period is shortened from 54.23 minutes ($4σ$) to 29.71 minutes ($3σ$). The possible QPO and period shortening phenomenon are best explained by the processes of magnetic reconnections. △ Less

Submitted 10 September, 2024; originally announced September 2024.

Two Dynamically Discovered Compact Object Candidate Binary Systems from LAMOST Low-resolution Survey

Authors: Senyu Qi, Wei-Min Gu, Zhi-Xiang Zhang, Tuan Yi, Jin-Zhong Liu, Ling-Lin Zheng

Abstract: We report two binary systems, LAMOST J035540+381550 (hereafter J035540) and LAMOST J035916+400732 (hereafter J035916), identified through the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) low-resolution survey (LRS). Each of these two systems contains an M-type star orbiting with a invisible compact object candidate. Follow-up spectroscopic observations of Palomar 200-inch tel… ▽ More We report two binary systems, LAMOST J035540+381550 (hereafter J035540) and LAMOST J035916+400732 (hereafter J035916), identified through the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) low-resolution survey (LRS). Each of these two systems contains an M-type star orbiting with a invisible compact object candidate. Follow-up spectroscopic observations of Palomar 200-inch telescope (P200) enhance radial velocity measurements. We use radial velocities from LAMOST and P200, as well as light curves from Zwicky Transient Facility (ZTF) to constrain orbital parameters. The masses of the visible M-type stars are estimated by fitting the MIST isochrones and SEDs. The mass functions for the unseen companions are: $0.22\pm 0.01 M_{\odot}$ for J035540 and $0.16\pm 0.01 M_{\odot}$ for J035916. With the orbital and stellar parameters derived above and assuming the orbital inclination is 90 degree (edge-on), we find that the minimum masses of the invisible companions exceeds that of the visible stars. The single-lined feature and the dynamical evidence suggest the presence of compact objects. J035540's ZTF light curve, modeled with PHOEBE, yields a compact object mass of $0.70^{+0.12}_{-0.05} M_{\odot}$. For J035916, ellipsoidal modulation analysis constrains the light curve amplitude, yielding a compact object mass range of $0.57-0.90 M_{\odot}$. The mass estimates indicate that both are likely white dwarfs. These findings underscore the efficiency of optical time-domain surveys and dynamical methods in identifying faint, massive white dwarfs, along with other compact objects in binaries. △ Less

Submitted 25 June, 2024; originally announced June 2024.

Comments: 12 pages, 9 figures, accepted for publication in MNRAS

A Sample of Compact Object Candidates in Single-lined Spectroscopic Binaries from LAMOST Medium Resolution Survey

Authors: Hao-Bin Liu, Wei-Min Gu, Zhi-Xiang Zhang, Tuan Yi, Jin-Zhong Liu, Mouyuan Sun

Abstract: The stellar spectra from LAMOST Medium Resolution Survey can be used to search for compact objects in binaries. The LAMOST DR10 catalog includes > 980, 000 targets with multiple medium resolution spectra. We select the targets with large or rapid radial velocity variation, and obtained an input-sample of 1822 sources. We use light curves and spectra to identify and exclude eclipsing binaries and d… ▽ More The stellar spectra from LAMOST Medium Resolution Survey can be used to search for compact objects in binaries. The LAMOST DR10 catalog includes > 980, 000 targets with multiple medium resolution spectra. We select the targets with large or rapid radial velocity variation, and obtained an input-sample of 1822 sources. We use light curves and spectra to identify and exclude eclipsing binaries and double-lined spectroscopic binaries in the input-sample. We finally derive a catalog of 89 candidates with well-folded radial velocity, which are all single-lined spectroscopic binaries, indicating an unseen companion residing in each system. The mass function of each system can be well constrained based on the radial velocity curve. In our sample, 26 sources have mass function higher than 0.1 $M_{\odot}$, among which 18 sources have ellipsoidal type light curves. In our opinion, compact objects are likely existent in all these 26 binaries, which are worth follow-up identification. △ Less

Submitted 16 May, 2024; originally announced May 2024.

Comments: 20 pages, 10 figures, accepted for publication in ApJ

The Binary Broadening Function

Authors: Tuan Yi

Abstract: We propose an extended formalism for the spectral broadening function (BF) based on the multiplication rule of block matrices. The formalism, which we named the binary broadening function (BBF), can produce decomposed BFs for individual components of a binary star system by using two spectral templates. The decomposed BFs can be used to derive precise rotational profiles and radial velocities for… ▽ More We propose an extended formalism for the spectral broadening function (BF) based on the multiplication rule of block matrices. The formalism, which we named the binary broadening function (BBF), can produce decomposed BFs for individual components of a binary star system by using two spectral templates. The decomposed BFs can be used to derive precise rotational profiles and radial velocities for individual components. We test the BBF on simulated spectra and actual observational spectra to show that the method is feasible on spectroscopic binaries, even when the spectral lines of two stellar components are heavily blended. To demonstrate the capability of the method, we conduct a simulation of `sketching' (imaging) a transiting circumbinary exoplanet using the BBF. We also discuss issues of implementation such as the variation of BBF with biased templates, the pros and cons of BBF, and cases when the method is not applicable. △ Less

Submitted 4 March, 2024; v1 submitted 30 January, 2024; originally announced January 2024.

Comments: 14 pages, 10 figures, 1 table. Accepted for publication in The Astrophysical Journal

The hidden companion in J1527: a 0.69 solar-mass white dwarf?

Authors: Zhi-Xiang Zhang, Hao-Bin Liu, Tuan Yi, Mouyuan Sun, Wei-Min Gu

Abstract: Finding nearby neutron stars can probe the supernova and metal-enrichment histories near our Solar system. Recently, Lin et al. 2023 reported an exciting neutron star candidate, 2MASS J15274848+3536572 (hereafter J1527), with a small Gaia distance of 118 parsecs. They claim that J1527 harbors an unseen neutron star candidate with an unusually low mass of $0.98\pm0.03\,M_{\odot}$. In this work, we… ▽ More Finding nearby neutron stars can probe the supernova and metal-enrichment histories near our Solar system. Recently, Lin et al. 2023 reported an exciting neutron star candidate, 2MASS J15274848+3536572 (hereafter J1527), with a small Gaia distance of 118 parsecs. They claim that J1527 harbors an unseen neutron star candidate with an unusually low mass of $0.98\pm0.03\,M_{\odot}$. In this work, we use the Canada-France-Hawaii Telescope high-resolution spectrum to measure J1527's orbital inclination independently. Our spectral fitting suggests an orbital inclination of $63\pm2$ degrees. Instead, by fitting a complex hybrid variability model consisting of the ellipsoidal-variation component and the star-spot modulation to the observed light curve, Lin et al. 2023 obtains an orbital inclination of $45.2_{-0.20}^{+0.13}$ degrees. We speculate that the orbital inclination obtained by the light-curve fitting is underestimated, since J1527's light curves are obviously not pure ellipsoidal variations. According to our new inclination ($i\sim 63$ degrees), the mass of the unseen compact object is reduced to $0.69\pm0.02\,M_\odot$, which is as massive as a typical white dwarf. △ Less

Submitted 16 January, 2024; originally announced January 2024.

Comments: 7 pages, 1 figure, 1 table, accepted for publication in The Astrophysical Journal Letters

Gaia22dkvLb: A Microlensing Planet Potentially Accessible to Radial-Velocity Characterization

Authors: Zexuan Wu, Subo Dong, Tuan Yi, Zhuokai Liu, Kareem El-Badry, Andrew Gould, L. Wyrzykowski, K. A. Rybicki, Etienne Bachelet, Grant W. Christie, L. de Almeida, L. A. G. Monard, J. McCormick, Tim Natusch, P. Zielinski, Huiling Chen, Yang Huang, Chang Liu, A. Merand, Przemek Mroz, Jinyi Shangguan, Andrzej Udalski, J. Woillez, Huawei Zhang, Franz-Josef Hambsch , et al. (28 additional authors not shown)

Abstract: We report discovering an exoplanet from following up a microlensing event alerted by Gaia. The event Gaia22dkv is toward a disk source rather than the traditional bulge microlensing fields. Our primary analysis yields a Jovian planet with M_p = 0.59^{+0.15}_{-0.05} M_J at a projected orbital separation r_perp = 1.4^{+0.8}_{-0.3} AU, and the host is a ~1.1 M_sun turnoff star at ~1.3 kpc. At r'~14,… ▽ More We report discovering an exoplanet from following up a microlensing event alerted by Gaia. The event Gaia22dkv is toward a disk source rather than the traditional bulge microlensing fields. Our primary analysis yields a Jovian planet with M_p = 0.59^{+0.15}_{-0.05} M_J at a projected orbital separation r_perp = 1.4^{+0.8}_{-0.3} AU, and the host is a ~1.1 M_sun turnoff star at ~1.3 kpc. At r'~14, the host is far brighter than any previously discovered microlensing planet host, opening up the opportunity of testing the microlensing model with radial velocity (RV) observations. RV data can be used to measure the planet's orbital period and eccentricity, and they also enable searching for inner planets of the microlensing cold Jupiter, as expected from the ''inner-outer correlation'' inferred from Kepler and RV discoveries. Furthermore, we show that Gaia astrometric microlensing will not only allow precise measurements of its angular Einstein radius theta_E, but also directly measure the microlens parallax vector and unambiguously break a geometric light-curve degeneracy, leading to definitive characterization of the lens system. △ Less

Submitted 30 May, 2024; v1 submitted 7 September, 2023; originally announced September 2023.

Comments: Accepted by AJ

Wide binaries with white dwarf or neutron star companions discovered from Gaia DR3 and LAMOST

Authors: Xinlin Zhao, Huijun Mu, Song Wang, Xue Li, Junhui Liu, Bowen Huang, Weimin Gu, Junfeng Wang, Tuan Yi, Zhixiang Zhang, Haibo Yuan, Zhongrui Bai, Hailong Yuan, Haotong Zhang, Jifeng Liu

Abstract: Gaia DR3 mission has identified and provided about 440,000 binary systems with orbital solutions, offering a valuable resource for searching binaries including a compact component. By combining the Gaia DR3 data with radial velocities (RVs) from the LAMOST spectroscopic survey, we identify three wide binaries possibly containing a compact object. For two of these sources with a main-sequence compa… ▽ More Gaia DR3 mission has identified and provided about 440,000 binary systems with orbital solutions, offering a valuable resource for searching binaries including a compact component. By combining the Gaia DR3 data with radial velocities (RVs) from the LAMOST spectroscopic survey, we identify three wide binaries possibly containing a compact object. For two of these sources with a main-sequence companion, no obvious excess is observed in the blue/red band of the Gaia DR3 XP spectra, and the LAMOST medium-resolution spectra exhibit clear single-lined features. The absence of an additional component from spectral disentangling analysis further suggests the presence of compact objects within these systems. On the other hand, the visible star of the third source is a stripped giant star. In contrast to most binaries including stripped stars, no emission line is detected in the optical spectra. The unseen star could potentially be a massive white dwarf or neutron star, but the possibility of an F-type dwarf star scenario cannot be ruled out. An examination of about ten binaries containing white dwarfs or neutron stars using both kinematic and chemical methods suggest most of these systems are located in the thin disk of the Milky Way. △ Less

Submitted 20 October, 2023; v1 submitted 6 August, 2023; originally announced August 2023.

Comments: 20 pages, 13 figures, accepted by APJ

Comprehensive study of the blazars from Fermi-LAT LCR: The log-normal flux distribution and linear RMS-Flux relation

Authors: Na Wang, Ting-Feng Yi, Liang Wang, Li-Sheng Mao, Zhi-Yuan Pu, Gong-Ming Ning, Wei-Tian Huang, He Lu, Shun Zhang, Yu-Tong Chen, Liang Dong

Abstract: Fermi-LAT LCR provide continuous and regularly-sampled gamma-ray light curves, spanning about 14 years, for a large sample of blazars. The log-normal flux distribution and linear RMS-Flux relation of the light curves for a few of Fermi blazar have been examined in previous studies. However, the probability that blazars exhibit log-normal flux distribution and linear RMS-Flux relation in their gamm… ▽ More Fermi-LAT LCR provide continuous and regularly-sampled gamma-ray light curves, spanning about 14 years, for a large sample of blazars. The log-normal flux distribution and linear RMS-Flux relation of the light curves for a few of Fermi blazar have been examined in previous studies. However, the probability that blazars exhibit log-normal flux distribution and linear RMS-Flux relation in their gamma-ray light curves has not been systematically explored. In this study, we comprehensively research on the distribution of gamma-ray flux and the statistical characteristics on a large sample of 1414 variable blazars from the Fermi-LAT LCR catalog, including 572 FSRQs, 477 BL Lacs, and 365 BCUs, and statistically compare their flux distributions with normal and log-normal distributions. The results indicate that the probability of not reject log-normal is 42.05% for the large sample, and there is still 2.05% probability of not reject normality, based on the joint of Kolmogorov-Smirnov, Shapiro-Wilk and Normality tests. We further find that the probability that BL Lacs conforms to the log-normal distribution is higher than that of FSRQs. Besides, after removing sources with less than 200 data points from this large sample, a sample of 549 blazars, which is still a large sample comparing to the previous studies, was obtained. Basing on dividing the light curves into segments every 20 points (or 40 points, or one year), we fitted the linear RMS-Flux relation of this three different sets, and found that the Pearson correlation coefficients are all close to 1 of the most blazars. This result indicates a strong linear correlation between the RMS and the flux of this 549 blazars. The log-normal distribution and linear RMS-Flux relation indicate that the variability of gamma-ray flux for most blazars is non-linear and multiplicative process. △ Less

Submitted 19 July, 2023; originally announced July 2023.

Comments: 13pages, 5figures, Accepted for publication in RAA

Two Transient Quasi-periodic Oscillations in $γ$-Ray Emission from the Blazar S4 0954+658

Authors: Yunlu Gong, Shiting Tian, Liancheng Zhou, Tingfeng Yi, Jun Fang

Abstract: In this work, we report periodicity search analyses in the gamma-ray light curve of the blazar S4 0954+658 monitoring undertaken by the Fermi Large Area Telescope (LAT). Four analytical methods and a tool are adopted to detect any periodic flux modulation and corresponding significance level, revealing that (i) a 66 d quasi-periodic oscillation (QPO) with the significance level of $> 5σ$ spanning… ▽ More In this work, we report periodicity search analyses in the gamma-ray light curve of the blazar S4 0954+658 monitoring undertaken by the Fermi Large Area Telescope (LAT). Four analytical methods and a tool are adopted to detect any periodic flux modulation and corresponding significance level, revealing that (i) a 66 d quasi-periodic oscillation (QPO) with the significance level of $> 5σ$ spanning over 600 d from 2015 to 2016 (MJD 57145--57745), resulting in continuous observation of nine cycles, which is one of the highest cycles discerned in blazar gamma-ray light curve; (ii) a possible QPO of 210 d at a moderate significance of $\sim3.5σ$ lasted for over 880 d from 2020 to 2022 (MJD 59035--59915), which lasted for four cycles. In addition, we discuss several physical models to explain the origin of the two transient QPOs and conclude that a geometrical scenario involving a plasma blob moving helically inside the jet can explain the time scale of the QPO. △ Less

Submitted 6 April, 2023; originally announced April 2023.

Comments: 10 pages, 11 figures; accepted for publication in ApJ

Searching for Compact Object Candidates from LAMOST Time-Domain Survey of Four K2 Plates

Authors: Senyu Qi, Wei-Min Gu, Tuan Yi, Zhi-Xiang Zhang, Song Wang, Jifeng Liu

Abstract: The time-domain (TD) surveys of the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) yield high-cadence radial velocities, paving a new avenue to study binary systems including compact objects. In this work, we explore LAMOST TD spectroscopic data of four K2 plates and present a sample of six single-lined spectroscopic binaries that may contain compact objects. We conduct analyse… ▽ More The time-domain (TD) surveys of the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) yield high-cadence radial velocities, paving a new avenue to study binary systems including compact objects. In this work, we explore LAMOST TD spectroscopic data of four K2 plates and present a sample of six single-lined spectroscopic binaries that may contain compact objects. We conduct analyses using phase-resolved radial velocity measurements of the visible star, to characterize each source and to infer the properties of invisible companion. By fitting the radial velocity curves for the six targets, we obtain accurate orbital periods, ranging from $\sim$ (0.6-6) days, and radial velocity semi-amplitudes, ranging from $\sim$ (50-130) km s$^{-1}$. We calculate the mass function of the unseen companions to be between 0.08 and 0.17 $M_{\odot}$. Based on the mass function and the estimated stellar parameters of the visible star, we determine the minimum mass of the hidden star. Three targets, J034813, J063350, and J064850, show ellipsoidal variability in the light curves from K2, ZTF, and TESS surveys. Therefore, we can put constraints on the mass of the invisible star using the ellipsoidal variability. We identify no X-ray counterparts for these targets except for J085120, of which the X-ray emission can be ascribed to stellar activity. We note that the nature of these six candidates is worth further characterization utilizing multi-wavelength follow-up observations. △ Less

Submitted 10 March, 2023; originally announced March 2023.

Comments: 13 pages, 6 figures, accepted for publication in The Astronomical Journal

Multicolour Optical Variability Monitoring of Blazars with High Time Resolution

Authors: X. Chang, T. F. Yi, D. R. Xiong, C. X. Liu, X. Yang, H. Z. Li, Y. L. Gong, W. W. Na, Y. Li, Z. H. Chen, J. P. Chen, L. S. Mao

Abstract: We carried out a high time-resolution, multicolour optical observing campaign for eight $γ$-ray detected blazars during 2010-2020. We analyze flux variations, correlations between magnitudes and colours on different timescales. Intraday variability (IDV) is detected in all eight sources of our sample. A bluer-when-brighter (BWB) chromatic trend is dominant on intraday timescales. On the short time… ▽ More We carried out a high time-resolution, multicolour optical observing campaign for eight $γ$-ray detected blazars during 2010-2020. We analyze flux variations, correlations between magnitudes and colours on different timescales. Intraday variability (IDV) is detected in all eight sources of our sample. A bluer-when-brighter (BWB) chromatic trend is dominant on intraday timescales. On the short timescales, the BWB trend only shows up in ON 231, 3C 279, BL Lacertae and 1E 1458.8+2249. There is a BWB trend in 3C 279 on the long timescale. We estimate the upper limits of black hole mass for three blazars (i.e. ON 321, 1ES 1426+42.8, PKS 1510-089) using variability timescales. On April 13, 2010 a potential quasi-periodic oscillation (QPO) with the period of $P=48.67\pm13.90$ minutes is found in 1ES 1426+42.8. The light curve on March 16, 2021 further shows the existence of the QPO phenomenon. The QPO in this target deserves further observation and confirmation. △ Less

Submitted 2 February, 2023; originally announced February 2023.

Comments: 18 pages, 8 figures, accepted in MNRAS

Effects of vertical advection on multimessenger signatures of black hole neutrino-dominated accretion flows in compact binary coalescences

Authors: Bing-Guang Chen, Tong Liu, Yan-Qing Qi, Bao-Quan Huang, Yun-Feng Wei, Tuan Yi, Wei-Min Gu, Li Xue

Abstract: In the coalescence events of binary neutron star (NS) or a black hole (BH) and an NS, a BH hyperaccretion disk might be eventually formed. At very high mass accretion rates, MeV neutrinos will be emitted from this disk, which is called a neutrino-dominated accretion flow (NDAF). Neutrino annihilation in the space out of the disk is energetic enough to launch ultrarelativistic jets to power gamma-r… ▽ More In the coalescence events of binary neutron star (NS) or a black hole (BH) and an NS, a BH hyperaccretion disk might be eventually formed. At very high mass accretion rates, MeV neutrinos will be emitted from this disk, which is called a neutrino-dominated accretion flow (NDAF). Neutrino annihilation in the space out of the disk is energetic enough to launch ultrarelativistic jets to power gamma-ray bursts. Moreover, vertical advection might exist in NDAFs, which can generate the magnetic buoyancy bubbles to release gamma-ray photons. In this paper, we visit the effects of the vertical advection in NDAFs on the disk structure and gamma-ray and neutrino luminosities for different accretion rates. Then we study the anisotropic emission of kilonovae and the following gravitational waves (GWs) driven by the gamma-ray photons and neutrinos from NDAFs. Comparing NDAFs without vertical advection, the neutrino luminosity and GW strains slightly decrease for the case with vertical advection, and the kilonovae will be brightened by the injected gamma-ray photons. The future joint multimessenger observations might distinguish whether the vertical advection exists in NDAFs or not after compact binary coalescences. △ Less

Submitted 16 November, 2022; originally announced November 2022.

Comments: 14 pages, 7 figures, accepted for publication in ApJ

The Nearest Neutron Star Candidate in a Binary Revealed by Optical Time-domain Surveys

Authors: Ling-Lin Zheng, Mouyuan Sun, Wei-Min Gu, Tuan Yi, Zhi-Xiang Zhang, Pei Wang, Junfeng Wang, Jianfeng Wu, Shan-Shan Weng, Song Wang, Sen-Yu Qi, Jia Zhang, Chun-Qian Li, Jian-Rong Shi, Yong Shao, Xiang-Dong Li, Jin-Bo Fu, Fan Yang, Zhongrui Bai, Yu Bai, Haotong Zhang, Jifeng Liu

Abstract: The near-Earth (within $\sim 100$ pc) supernova explosions in the past several million years can cause the global deposition of radioactive elements (e.g., $^{60}$Fe) on Earth. The remnants of such supernovae are too old to be easily identified. It is therefore of great interest to search for million-year-old near-Earth neutron stars or black holes, the products of supernovae. However, neutron sta… ▽ More The near-Earth (within $\sim 100$ pc) supernova explosions in the past several million years can cause the global deposition of radioactive elements (e.g., $^{60}$Fe) on Earth. The remnants of such supernovae are too old to be easily identified. It is therefore of great interest to search for million-year-old near-Earth neutron stars or black holes, the products of supernovae. However, neutron stars and black holes are challenging to find even in our Solar neighbourhood if they are not radio pulsars or X-ray/$γ$-ray emitters. Here we report the discovery of one of the nearest ($127.7 \pm 0.3$ pc) neutron star candidates in a detached single-lined spectroscopic binary LAMOST J235456.73+335625.9 (hereafter J2354). Utilizing the time-resolved ground-based spectroscopy and space photometry, we find that J2354 hosts an unseen compact object with $M_{\mathrm{inv}}$ being $1.4 \sim 1.6\ M_{\odot}$. The follow-up Swift ultraviolet (UV) and X-ray observations suggest that the UV and X-ray emission is produced by the visible star rather than the compact object. Hence, J2354 probably harbours a neutron star rather than a hot ultramassive white dwarf. Two-hour exceptionally sensitive radio follow-up observations with Five-hundred-meter Aperture Spherical radio Telescope fail to reveal any pulsating radio signals at the $6σ$ flux upper limit of $12.5\ μ\mathrm{Jy}$. Therefore, the neutron star candidate in J2354 can only be revealed via our time-resolved observations. Interestingly, the distance between J2354 and our Earth can be as close as $\sim 50$ pc around $2.5$ Myrs ago, as revealed by the Gaia kinematics. Our discovery demonstrates a promising way to unveil the hidden near-Earth neutron stars in binaries by exploring the optical time domain, thereby facilitating understanding of the metal-enrichment history in our Solar neighbourhood. △ Less

Submitted 27 November, 2023; v1 submitted 7 October, 2022; originally announced October 2022.

Comments: 26 pages, 13 figures

Journal ref: Sci. China-Phys. Mech. Astron. 66, 129512 (2023)

Compact Object Candidates with K/M-dwarf Companions from LAMOST Low-resolution Survey

Authors: Hui-Jun Mu, Wei-Min Gu, Tuan Yi, Ling-Lin Zheng, Hao Sou, Zhong-Rui Bai, Hao-Tong Zhang, Ya-Juan Lei, Cheng-Ming Li

Abstract: Searching for compact objects (black holes, neutron stars, or white dwarfs) in the Milky Way is essential for understanding the stellar evolution history, the physics of compact objects, and the structure of our Galaxy. Compact objects in binaries with a luminous stellar companion are perfect targets for optical observations. Candidate compact objects can be achieved by monitoring the radial veloc… ▽ More Searching for compact objects (black holes, neutron stars, or white dwarfs) in the Milky Way is essential for understanding the stellar evolution history, the physics of compact objects, and the structure of our Galaxy. Compact objects in binaries with a luminous stellar companion are perfect targets for optical observations. Candidate compact objects can be achieved by monitoring the radial velocities of the companion star. However, most of the spectroscopic telescopes usually obtain stellar spectra at a relatively low efficiency, which makes a sky survey for millions of stars practically impossible. The efficiency of a large-scale spectroscopic survey, the Large Sky Area Multi-Object Fiber Spectroscopy Telescope (LAMOST), presents a specific opportunity to search for compact object candidates, i.e., simply from the spectroscopic observations. Late-type K/M stars are the most abundant populations in our Galaxy. Owing to the relatively large Keplerian velocities in the close binaries with a K/M-dwarf companion, a hidden compact object could be discovered and followed-up more easily. In this study, compact object candidates with K/M-dwarf companions are investigated with the LAMOST low-resolution stellar spectra. Based on the LAMOST Data Release 5, we obtained a sample of $56$ binaries, each containing a K/M-dwarf with a large radial velocity variation $ΔV_{\rm R} > 150~{\rm km~s}^{-1}$. Complemented with the photometric information from the Transiting Exoplanet Survey Satellite, we derived a sample of $35$ compact object candidates, among which, the orbital periods of $16$ sources were revealed by the light curves. Considering two sources as examples, we confirmed that a compact object existed in the two systems by fitting the radial velocity curve. This study demonstrates the principle and the power of searching for compact objects through LAMOST. △ Less

Submitted 28 September, 2022; originally announced September 2022.

Comments: 20 pages, 7 figures, 4 tables, Science China Physics, Mechanics & Astronomy, Vol.65 No.2:229711

One-off and Repeating Fast Radio Bursts: A Statistical Analysis

Authors: Hao-Yan Chen, Wei-Min Gu, Mouyuan Sun, Tuan Yi

Abstract: According to the number of detected bursts, fast radio bursts (FRBs) can be classified into two categories, i.e., one-off FRBs and repeating ones. We make a statistical comparison of these two categories based on the first FRB catalog of the Canadian Hydrogen Intensity Mapping Experiment Fast Radio Burst Project. Using the Anderson-Darling, Kolmogrov-Smirnov, and Energy statistic tests, we find si… ▽ More According to the number of detected bursts, fast radio bursts (FRBs) can be classified into two categories, i.e., one-off FRBs and repeating ones. We make a statistical comparison of these two categories based on the first FRB catalog of the Canadian Hydrogen Intensity Mapping Experiment Fast Radio Burst Project. Using the Anderson-Darling, Kolmogrov-Smirnov, and Energy statistic tests, we find significant statistical differences ($p$-value $<$ 0.001) of the burst properties between the one-off FRBs and the repeating ones. More specifically, after controlling for distance, we find that the peak luminosities of one-off FRBs are, on average, higher than the repeating ones; the pulse temporal widths of repeating FRBs are, on average, longer than the one-off ones. The differences indicate that these two categories could have distinct physical origins. Moreover, we discuss the sub-populations of FRBs and provide statistical evidence to support the existence of sub-populations in one-off FRBs and in repeating ones. △ Less

Submitted 28 September, 2022; originally announced September 2022.

Comments: 10 pages, 6 figures, accepted for publication in ApJ

A White Dwarf-Main Sequence Binary Unveiled by Time-Domain Observations from LAMOST and TESS

Authors: Ling-Lin Zheng, Wei-Min Gu, Mouyuan Sun, Zhixiang Zhang, Tuan Yi, Jianfeng Wu, Junfeng Wang, Jin-Bo Fu, Sen-Yu Qi, Fan Yang, Song Wang, Liang Wang, Zhongrui Bai, Haotong Zhang, Chun-Qian Li, Jian-Rong Shi, Weikai Zong, Yu Bai, Jifeng Liu

Abstract: We report a single-lined white dwarf-main sequence binary system, LAMOST J172900.17+652952.8, which is discovered by LAMOST's medium resolution time-domain surveys. The radial velocity semi-amplitude and orbital period of the optical visible star are measured by using the Palomar 200-inch telescope follow-up observations and the light curves from TESS. Thus the mass function of the invisible candi… ▽ More We report a single-lined white dwarf-main sequence binary system, LAMOST J172900.17+652952.8, which is discovered by LAMOST's medium resolution time-domain surveys. The radial velocity semi-amplitude and orbital period of the optical visible star are measured by using the Palomar 200-inch telescope follow-up observations and the light curves from TESS. Thus the mass function of the invisible candidate white dwarf is derived, $f(M_{\rm{2}}) = 0.120\,\pm\,0.003\,M_{\odot}$. The mass of the visible star is measured based on the spectral energy distribution fitting, $M_{\mathrm{1}}$ = $0.81^{+0.07}_{-0.06}\,M_{\odot}$. Hence, the mass of its invisible companion is $M_{\rm{2}}\,\gtrsim\,0.63\,M_{\odot}$. The companion ought to be a compact object rather than a main-sequence star owing to the mass ratio $q = M_{\rm{2}} / M_{\rm 1} \gtrsim 0.78$ and the single-lined spectra. The compact object is likely to be a white dwarf except for small inclination angle, $i\,\lesssim\,40^{\circ}$. By using the GALEX NUV flux, the effective temperature of the white dwarf candidate is constrained as $T_{\rm eff}^{\rm WD}\,\lesssim\,12000-13500$ K. It is difficult to detect white dwarfs which are outshone by their bright companions via single-epoch optical spectroscopic surveys. Therefore, the optical time-domain surveys can play an important role in unveiling invisible white dwarfs and other compact objects in binaries. △ Less

Submitted 28 September, 2022; originally announced September 2022.

Comments: 15 pages, 6 figures, ApJ, 936, 33

A long-period pre-ELM system discovered from LAMOST medium-resolution survey

Authors: Zhi-Xiang Zhang, Ling-Lin Zheng, Wei-Min Gu, Mouyuan Sun, Tuan Yi, Jian-Rong Shi, Song Wang, Zhong-Rui Bai, Hao-Tong Zhang, Wen-Yuan Cui, Junfeng Wang, Jianfeng Wu, Xiang-Dong Li, Yong Shao, Kai-Xing Lu, Yu Bai, Chunqian Li, Jin-Bo Fu, Jifeng Liu

Abstract: We present LAMOST~J041920.07+072545.4 (hereafter J0419), a close binary consisting of a bloated extremely low mass pre-white dwarf (pre-ELM WD) and a compact object with an orbital period of 0.607189~days. The large-amplitude ellipsoidal variations and the evident Balmer and He~I emission lines suggest a filled Roche lobe and ongoing mass transfer. No outburst events were detected in the 15 years… ▽ More We present LAMOST~J041920.07+072545.4 (hereafter J0419), a close binary consisting of a bloated extremely low mass pre-white dwarf (pre-ELM WD) and a compact object with an orbital period of 0.607189~days. The large-amplitude ellipsoidal variations and the evident Balmer and He~I emission lines suggest a filled Roche lobe and ongoing mass transfer. No outburst events were detected in the 15 years of monitoring of J0419, indicating a very low mass transfer rate. The temperature of the pre-ELM, $T_\mathrm{eff} = 5793_{-133}^{+124}\,\rm K$, is cooler than the known ELMs, but hotter than most CV donors. Combining the mean density within the Roche lobe and the radius constrained from our SED fitting, we obtain the mass of the pre-ELM, $M_1 = 0.176\pm 0.014\,M_\odot$. The joint fitting of light and radial velocity curves yields an inclination angle of $i = 66.5_{-1.7}^{+1.4}$ degrees, corresponding to the compact object mass of $M_2 = 1.09\pm 0.05\,M_\odot$. The very bloated pre-ELM has a smaller surface gravity ($\log g = 3.9\pm 0.01$, $R_1 = 0.78 \pm 0.02\,R_\odot$) than the known ELMs or pre-ELMs. The temperature and the luminosity ($L_\mathrm{bol} = 0.62_{-0.10}^{+0.11}\,L_\odot$) of J0419 are close to the main sequence, which makes the selection of such systems through the HR diagram inefficient. Based on the evolutionary model, the relatively long period and small $\log g$ indicate that J0419 could be close to the "bifurcation period" in the orbit evolution, which makes J0419 to be a unique source to connect ELM/pre-ELM WD systems, wide binaries and cataclysmic variables. △ Less

Submitted 28 September, 2022; originally announced September 2022.

Comments: 25 pages, 13 figures, ApJ, 933, 193

A dynamically discovered and characterized non-accreting neutron star -- M dwarf binary candidate

Authors: Tuan Yi, Wei-Min Gu, Zhi-Xiang Zhang, Ling-Lin Zheng, Mouyuan Sun, Junfeng Wang, Zhongrui Bai, Pei Wang, Jianfeng Wu, Yu Bai, Song Wang, Haotong Zhang, Yize Dong, Yong Shao, Xiang-Dong Li, Jia Zhang, Yang Huang, Fan Yang, Qingzheng Yu, Hui-Jun Mu, Jin-Bo Fu, Senyu Qi, Jing Guo, Xuan Fang, Chuanjie Zheng , et al. (4 additional authors not shown)

Abstract: Optical time-domain surveys can unveil and characterize exciting but less-explored non-accreting and/or non-beaming neutron stars (NS) in binaries. Here we report the discovery of such a NS candidate using the LAMOST spectroscopic survey. The candidate, designated LAMOST J112306.9+400736 (hereafter J1123), is in a single-lined spectroscopic binary containing an optically visible M star. The star's… ▽ More Optical time-domain surveys can unveil and characterize exciting but less-explored non-accreting and/or non-beaming neutron stars (NS) in binaries. Here we report the discovery of such a NS candidate using the LAMOST spectroscopic survey. The candidate, designated LAMOST J112306.9+400736 (hereafter J1123), is in a single-lined spectroscopic binary containing an optically visible M star. The star's large radial velocity variation and ellipsoidal variations indicate a relatively massive unseen companion. Utilizing follow-up spectroscopy from the Palomar 200-inch telescope and high-precision photometry from TESS, we measure a companion mass of $1.24_{-0.03}^{+0.03}~M_{\odot}$. Main-sequence stars with this mass are ruled out, leaving a NS or a massive white dwarf (WD). Although a massive WD cannot be ruled out, the lack of UV excess radiation from the companion supports the NS hypothesis. Deep radio observations with FAST yielded no detections of either pulsed or persistent emission. J1123 is not detected in numerous X-ray and gamma-ray surveys. These non-detections suggest that the NS candidate is not presently accreting and pulsing. Our work exemplifies the capability of discovering compact objects in non-accreting close binaries by synergizing the optical time-domain spectroscopy and high-cadence photometry. △ Less

Submitted 25 September, 2022; originally announced September 2022.

Comments: 53 pages, 15 figures, publication in Nature Astronomy

A 31.3 day Transient Quasiperiodic Oscillation in Gamma-ray Emission from Blazar S5 0716+714

Authors: Junping Chen, Tingfeng Yi, Yunlu Gong, Xing Yang, Zhihui Chen, Xin Chang, Lisheng Mao

Abstract: We systematically search for quasiperiodic oscillatory (QPO) signals on the month timescale among the 1525 sources given in the Fermi Large Area Telescope Light Curve Repository. We find a transient QPO of 31.3$\pm$1.8 days in the gamma-ray band light curve of the TeV blazar S5 0716+714, which has seven cycles (MJD 55918-56137) for the first time by weighted wavelet Z-transform and Lomb-Scargle pe… ▽ More We systematically search for quasiperiodic oscillatory (QPO) signals on the month timescale among the 1525 sources given in the Fermi Large Area Telescope Light Curve Repository. We find a transient QPO of 31.3$\pm$1.8 days in the gamma-ray band light curve of the TeV blazar S5 0716+714, which has seven cycles (MJD 55918-56137) for the first time by weighted wavelet Z-transform and Lomb-Scargle periodogram methods. Monte Carlo simulations based on the power spectral density and probability distribution function were used to evaluate the confidence level of the QPO, and the result is $\sim 4.1σ$. Seasonal autoregressive integrated moving average modeling of the light curve revealed it is a significant physical QPO. The physical models to explain the sporadic QPO of the month-timescale QPOs in blazar were discussed. Our studies indicate that the helical jet model and blob move helically in a curved jet model to properly explain this kind of transient QPO. △ Less

Submitted 17 October, 2022; v1 submitted 14 September, 2022; originally announced September 2022.

Comments: 13 pages, 5 figures, accepted for publication in ApJ

Searching for compact objects in binaries with Gaia DR3

Authors: Jin-Bo Fu, Wei-Min Gu, Zhi-Xiang Zhang, Tuan Yi, Sen-Yu Qi, Ling-Lin Zheng, Jifeng Liu

Abstract: We search for compact objects in binaries based on Gaia DR3. A sample of ten targets is derived under the conditions: radial velocity variable, low temperature ($T_{\rm eff} < 6000$ K), high mass function ($f(M_2) > 1 M_\odot$), and ellipsoidal-like light curves. Two targets have LAMOST spectroscopic observations, one of which is a double-lined spectroscopic binary. The observational data of seven… ▽ More We search for compact objects in binaries based on Gaia DR3. A sample of ten targets is derived under the conditions: radial velocity variable, low temperature ($T_{\rm eff} < 6000$ K), high mass function ($f(M_2) > 1 M_\odot$), and ellipsoidal-like light curves. Two targets have LAMOST spectroscopic observations, one of which is a double-lined spectroscopic binary. The observational data of seven targets are not self-consistent, since their photometric periods are even shorter than the theoretical minimum orbital periods calculated by the stellar parameters from Gaia DR3. After excluding these seven inconsistent targets and another target contaminated by a near-bright star, the remaining two targets may contain compact objects worth follow-up observations. This work may serve as an example to demonstrate the feasibility of searching for compact objects in the massive Gaia data. △ Less

Submitted 18 October, 2022; v1 submitted 12 July, 2022; originally announced July 2022.

Comments: 15 pages, 8 figures, accepted for publication in ApJ

Quasi-periodic behaviour in the $γ$-ray light curve of the blazar PKS 0405-385

Authors: Yunlu Gong, Liancheng Zhou, Min Yuan, Haiyun Zhang, Tingfeng Yi, Jun Fang

Abstract: We analyze the quasi-periodic oscillation (QPO) of the historical light curve of FSRQs PKS 0405-385 detected by the Fermi LAT from August 2008 to November 2021. To identify and determine the QPO signal of PKS 0405-385 in the $γ$-ray light curve, we use four time series analysis techniques based on frequency and time domains, i.e., the Lomb-Scargle periodogram (LSP), the weighted wavelet z-transfor… ▽ More We analyze the quasi-periodic oscillation (QPO) of the historical light curve of FSRQs PKS 0405-385 detected by the Fermi LAT from August 2008 to November 2021. To identify and determine the QPO signal of PKS 0405-385 in the $γ$-ray light curve, we use four time series analysis techniques based on frequency and time domains, i.e., the Lomb-Scargle periodogram (LSP), the weighted wavelet z-transform (WWZ), the REDFIT and the epoch folding. The results show that PKS 0405-385 has a quasi-periodic behavior of $\sim$2.8 yr with the significance of $\sim$4.3$σ$ in Fermi long-term monitoring. Remarkably, we also performed QPO analysis in the G-band light curve observed from October 2014 to October 2021 using LSP and WWZ technology, and the results ($\sim$4$σ$ of significance) are consistent with the periodic detection in $γ$-ray. This may imply that the optical emission is radiated by an electron population same as the \gr\ emission. In discussing the possible mechanism of quasi-periodic behavior, either the helical motion within a jet or the supermassive black hole binary system provides a viable explanation for the QPO of 2.8 yr, and the relevant parameters have been estimated. △ Less

Submitted 4 May, 2022; originally announced May 2022.

Comments: 10 pages, 8 figures, Accepted for publication in ApJ

Effects of Gravitational Wave Radiation of Eccentric Neutron Star-White Dwarf Binaries on the Periodic Activity of Fast Radio Burst Sources

Authors: Yi-Qing Lin, Hao-Yan Chen, Wei-Min Gu, Tuan Yi

Abstract: We revisit the eccentric neutron star (NS)-white dwarf (WD) binary model for the periodic activity of fast radio burst (FRB) sources, by including the effects of gravitational wave (GW) radiation. In this model, the WD fills its Roche lobe at the periastron and mass transfer occurs from the WD to the NS. The accreted materials can be fragmented and arrive at the NS episodically, resulting in multi… ▽ More We revisit the eccentric neutron star (NS)-white dwarf (WD) binary model for the periodic activity of fast radio burst (FRB) sources, by including the effects of gravitational wave (GW) radiation. In this model, the WD fills its Roche lobe at the periastron and mass transfer occurs from the WD to the NS. The accreted materials can be fragmented and arrive at the NS episodically, resulting in multiple bursts through curvature radiation. Consequently, the WD may be kicked away owing to the conservation of angular momentum. To initiate the next mass transfer, the WD has to refill its Roche lobe through GW radiation. In this scenario, whether the periodic activity can show up relies on three timescales, i.e., the orbital period $P_{\rm orb}$, the timescale $T_{\rm GW}$ for the Roche lobe to be refilled, and the time span $T_{\rm frag}$ for all the episodic events corresponding to each mass transfer process. Only when the two conditions $T_{\rm GW} \gtrsim P_{\rm orb}$ and $T_{\rm frag} < P_{\rm orb}$ are both satisfied, the periodic activity will manifest itself and the period should be equal to $P_{\rm orb}$. In this spirit, the periodic activity is more likely to show up for relatively long periods ($P_{\rm orb} \gtrsim$ several days). Thus, it is reasonable that FRBs 180916 and 121102, the only two sources having been claimed to manifest periodic activity, both correspond to relatively long periods. △ Less

Submitted 11 March, 2022; originally announced March 2022.

Comments: 7 pages, 2 figures, accepted for pubilcation in ApJ

GRB 211227A as a peculiar long gamma-ray burst from compact star merger

Authors: Hou-Jun Lü, Hao-Yu Yuan, Ting-Feng Yi, Xiang-Gao Wang, You-Dong Hu, Yong Yuan, Jared Rice, Jian-Guo Wang, Jia-Xin Cao, De-Feng Kong, Emilio Fernandez-Garc�a, Alberto J. Castro-Tirado, Ji-Shun Lian, Wen-Pei Gan, Shan-Qin Wang, Li-Ping Xin, M. D. Caballero-Garc�a, Yu-Feng Fan, En-Wei Liang

Abstract: Long-duration gamma-ray bursts (GRBs) associated with supernovae (SNe) are believed to originate from massive star core-collapse events, whereas short-duration GRBs that are related to compact star mergers are expected to be accompanied by kilonovae. GRB 211227A, which lasted about 84 s, had an initial short/hard spike followed by a series of soft gamma-ray extended emission at redshift $z=$0.228.… ▽ More Long-duration gamma-ray bursts (GRBs) associated with supernovae (SNe) are believed to originate from massive star core-collapse events, whereas short-duration GRBs that are related to compact star mergers are expected to be accompanied by kilonovae. GRB 211227A, which lasted about 84 s, had an initial short/hard spike followed by a series of soft gamma-ray extended emission at redshift $z=$0.228. We performed follow-up observations of the optical emission using BOOTES, LCOGT, and the Lijiang 2.4m telescope, but we detected no associated supernova signature, even down to very stringent limits at such a low redshift. We observed the host galaxy within a large error-circle and roughly estimate the physical offset of GRB 211227A as $20.47\pm14.47$ kpc from the galaxy center. These properties are similar to those of GRB 060614, and suggest that the progenitor of GRB 211227A is not favored to be associated with the death of massive stars. Hence, we propose that GRB 211227A originates from a compact star merger. Calculating pseudo-kilonova emission for this case by adopting the typical parameters, we find that any associated pseudo-kilonova is too faint to be detected. If this is the case, it explains naturally the characteristics of the prompt emission, the lack of SN and kilonova emission, and the large physical offset from the galaxy center. △ Less

Submitted 30 May, 2022; v1 submitted 17 January, 2022; originally announced January 2022.

Comments: 12 pages, 7 figures, and 1 Table. Matched with the published verison

Journal ref: ApJL, 931:L23, 2022

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

Authors: Wan-Min Zheng, Qiaoya Wu, Jianfeng Wu, Song Wang, Mouyuan Sun, Jing Guo, Junhui Liu, Tuan Yi, Zhi-Xiang Zhang, Wei-Min Gu, Junfeng Wang, Lijun Gou, Jifeng Liu, Paul J. Callanan, Luis C. Ho, Pen�lope Longa-Pe�a, Jerome A. Orosz, Mark T. Reynolds

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

Submitted 14 December, 2021; originally announced December 2021.

Comments: Accepted for publication in ApJ

Journal ref: 2022ApJ...925...83Z

LTD064402+245919: A Subgiant with a 1-3 M$_{\odot}$ Undetected Companion Identified from LAMOST-TD Data

Authors: Fan Yang, Bo Zhang, Richard J. Long, You-Jun Lu, Su-Su Shan, Xing Wei, Jian-Ning Fu, Xian-Fei Zhang, Zhi-Chao Zhao, Yu Bai, Tuan Yi, Ling-Lin Zheng, Ze-Ming Zhou, Ji-Feng Liu

Abstract: Single-line spectroscopic binaries recently contribute to the stellar-mass black hole discovery, independently of the X-ray transient method. We report the identification of a single-line binary system LTD064402+245919, with an orbital period of 14.50 days. The observed component is a subgiant with a mass of 2.77$\pm$0.68M$_{\odot}$, radius 15.5$\pm$2.5R$_{\odot}$, effective temperature… ▽ More Single-line spectroscopic binaries recently contribute to the stellar-mass black hole discovery, independently of the X-ray transient method. We report the identification of a single-line binary system LTD064402+245919, with an orbital period of 14.50 days. The observed component is a subgiant with a mass of 2.77$\pm$0.68M$_{\odot}$, radius 15.5$\pm$2.5R$_{\odot}$, effective temperature $T_{\rm eff}$ 4500$\pm$200K, and surface gravity log\emph{g} 2.5$\pm$0.25dex. The discovery makes use of the LAMOST time-domain (LAMOST-TD) and ZTF survey. Our general-purpose software pipeline applies the Lomb-Scargle periodogram to determine the orbital period and uses machine-learning to classify the variable type from the folded light curves. We apply a combined model to estimate the orbital parameters from both the light and radial velocity curves, taking constraints on the primary star mass, mass function, and detection limit of secondary luminosity into consideration. We obtain a radial velocity semi-amplitude of 44.6$\pm$1.5 km s$^{-1}$, mass ratio of 0.73$\pm$0.07, and an undetected component mass of 2.02$\pm$0.49M$_{\odot}$ when the type of the undetected component is not set. We conclude that the inclination is not well constrained, and that the secondary mass is larger than 1M$_{\odot}$ when the undetected component is modelled as a compact object. According to our investigations using an MCMC simulation, increasing the spectra SNR by a factor of 3 would enable the secondary light to be distinguished (if present). The algorithm and software in this work are able to serve as general-purpose tools for the identification of compact objects quiescent in X-rays. △ Less

Submitted 24 October, 2021; v1 submitted 15 October, 2021; originally announced October 2021.

Comments: Accepted by ApJ

Reconciling the 16.35-day period of FRB 20180916B with jet precession

Authors: Hao-Yan Chen, Wei-Min Gu, Mouyuan Sun, Tong Liu, Tuan Yi

Abstract: A repeating fast radio burst (FRB), FRB 20180916B (hereafter FRB 180916), was reported to have a 16.35-day period. This period might be related to a precession period. In this paper, we investigate two precession models to explain the periodic activity of FRB 180916. In both models, the radio emission of FRB 180916 is produced by a precessing jet. For the first disk-driven jet precession model, an… ▽ More A repeating fast radio burst (FRB), FRB 20180916B (hereafter FRB 180916), was reported to have a 16.35-day period. This period might be related to a precession period. In this paper, we investigate two precession models to explain the periodic activity of FRB 180916. In both models, the radio emission of FRB 180916 is produced by a precessing jet. For the first disk-driven jet precession model, an extremely low viscous parameter (i.e., the dimensionless viscosity parameter $α\lesssim 10^{-8}$) is required to explain the precession of FRB 180916, which implies its implausibility. For the second tidal force-driven jet precession model, we consider a compact binary consists of a neutron star/black hole and a white dwarf; the white dwarf fills its Roche lobe and mass transfer occurs. Due to the misalignment between the disk and orbital plane, the tidal force of the white dwarf can drive jet precession. We show that the relevant precession periods are several days to hundreds of days, depending on the specific accretion rates and component masses. The duration of FRB 180916 generation in the binary with extremely high accretion rate will be several thousand years. △ Less

Submitted 19 August, 2021; originally announced August 2021.

Comments: 12 pages, 4 figures, accepted for publication in ApJ

An Empirical Bayesian Approach to Limb-darkening in Modeling WASP-121b Transit Light Curves

Authors: Fan Yang, Richard J. Long, Ji-Feng Liu, Su-Su Shan, Rui Guo, Bo Zhang, Tuan Yi, Ling-Lin Zheng, Zhi-Chao Zhao

Abstract: We present a novel, iterative method using an empirical Bayesian approach for modeling the limb darkened WASP-121b transit from the TESS light curve. Our method is motivated by the need to improve $R_{p}/R_{\ast}$ estimates for exoplanet atmosphere modeling, and is particularly effective with the limb darkening (LD) quadratic law requiring no prior central value from stellar atmospheric models. Wi… ▽ More We present a novel, iterative method using an empirical Bayesian approach for modeling the limb darkened WASP-121b transit from the TESS light curve. Our method is motivated by the need to improve $R_{p}/R_{\ast}$ estimates for exoplanet atmosphere modeling, and is particularly effective with the limb darkening (LD) quadratic law requiring no prior central value from stellar atmospheric models. With the non-linear LD law, the method has all the advantages of not needing atmospheric models but does not converge. The iterative method gives a different $R_{p}/R_{\ast}$ for WASP-121b at a significance level of 1$σ$ when compared with existing non-iterative methods. To assess the origins and implications of this difference, we generate and analyze light curves with known values of the limb darkening coefficients (LDCs). We find that non-iterative modeling with LDC priors from stellar atmospheric models results in an inconsistent $R_{p}/R_{\ast}$ at 1.5$σ$ level when the known LDC values are as those previously found when modeling real data by the iterative method. In contrast, the LDC values from the iterative modeling yields the correct value of $R_{p}/R_{\ast}$ to within 0.25$σ$. For more general cases with different known inputs, Monte Carlo simulations show that the iterative method obtains unbiased LDCs and correct $R_{p}/R_{\ast}$ to within a significance level of 0.3$σ$. Biased LDC priors can cause biased LDC posteriors and lead to bias in the $R_{p}/R_{\ast}$ of up to 0.82$\%$, 2.5$σ$ for the quadratic law and 0.32$\%$, 1.0$σ$ for the non-linear law. Our improvement in $R_{p}/R_{\ast}$ estimation is important when analyzing exoplanet atmospheres. △ Less

Submitted 15 April, 2021; originally announced April 2021.

Comments: Accepted by AJ

Cosmological Evolution of Fermi-LAT Gamma-ray Blazars Using Novel Nonparametric Methods

Authors: Houdun Zeng, Vahé Petrosian, Tingfeng Yi

Abstract: Multi-wavelength analyses of spectra of active galactic nuclei (AGNs) provide useful information on the physical processes in the accretion disk and jets of black holes. This, however, is limited to bright sources and may not represent the population as a whole. Another approach is through the investigation of the cosmological evolution of the luminosity function (LF) that shows varied evolution (… ▽ More Multi-wavelength analyses of spectra of active galactic nuclei (AGNs) provide useful information on the physical processes in the accretion disk and jets of black holes. This, however, is limited to bright sources and may not represent the population as a whole. Another approach is through the investigation of the cosmological evolution of the luminosity function (LF) that shows varied evolution (luminosity and density) at different wavelengths. These differences and the correlations between luminosities can shed light on the Jet-accretion disk connection. Most such studies use forward fitting parametric methods that involve several functions and many parameters. We use non-parametric, non-binning methods developed by Efron and Petrosian, and Lynden-Bell, for obtaining unbiased description of the evolution of the LF, from data truncated by observational selection effects. We present analysis of the evolution of gamma-ray LF of blazars with main focus on flat spectrum radio quasars (FSRQs). This requires analysis of both gamma-ray and optical data, essential for redshift measurements, and a description of the joint LF. We use a new approach which divides the sample into two sub-samples, each with its own flux limit. We use the Fermi-LAT and GAIA observations, and present results on the gamma-ray LF and its evolution, and determine the intrinsic correlation between the gamma-ray and optical luminosities corrected for the well known false correlation induced by their similar redshift dependence and evolution of the two luminosities. We also present a direct estimation of the contribution of blazars to the spectrum of the extragalactic gamma-ray background. △ Less

Submitted 10 April, 2021; originally announced April 2021.

Comments: 46 pages, 11 figures, 2 tables, accepted for publication in ApJ

Highly Variable $γ$-Ray Emission of CTD 135 and Implications for its Compact Symmetric Structure

Authors: Ying-Ying Gan, Hai-Ming Zhang, Jin Zhang, Xing Yang, Ting-Feng Yi, Yun-Feng Liang, En-Wei Liang

Abstract: The $γ$-ray emission properties of CTD 135, a typical compact symmetric object (CSO), are investigated with $\sim$11-year Fermi/LAT observations. We show that it has bright and significantly variable GeV emission, with the $γ$-ray luminosity of $L_γ\sim10^{47}$ erg s$^{-1}$ and a variation index of TS$_{\rm var}=1002$. A quasi-periodic oscillation (QPO) with a periodicity of $\sim$460 days is dete… ▽ More The $γ$-ray emission properties of CTD 135, a typical compact symmetric object (CSO), are investigated with $\sim$11-year Fermi/LAT observations. We show that it has bright and significantly variable GeV emission, with the $γ$-ray luminosity of $L_γ\sim10^{47}$ erg s$^{-1}$ and a variation index of TS$_{\rm var}=1002$. A quasi-periodic oscillation (QPO) with a periodicity of $\sim$460 days is detected in the global 95\% false-alarm level. These $γ$-ray emission features are similar to that of blazars. Its broadband spectral energy distribution (SED) can be attributed to the radiations of the relativistic electrons accelerated in the core region and the extended region. The SED modeling shows that the $γ$-rays are from the core region, which has a Doppler boosting factor of $δ\sim10.8$ and relativistically moves with a small viewing angle, being similar to blazar jets. On the base of the analysis results, we propose that the episodic activity of the central engine in CTD 135 results in a blazar-like jet and the bubble-like lobes as the Fermi bubbles in the Galaxy. The strong $γ$-ray emission with obvious variability is from the jet radiations and the symmetric radio structure is attributed to the bubbles. The jet radiation power and disk luminosity in units of Eddington luminosity of CTD 135 follow the same relation as other young radio galaxies, indicating that its jet radiation may be also driven by the Eddington ratio. △ Less

Submitted 9 April, 2021; originally announced April 2021.

Comments: 14 pages, 7 figures, accepted for publication in RAA

Jet Properties of Compact Steep-Spectrum Sources and an Eddington-Ratio-Driven Unification Scheme of Jet Radiation in Active Galactic Nuclei

Authors: Jin Zhang, Hai-Ming Zhang, Ying-Ying Gan, Ting-Feng Yi, Jun-Feng Wang, En-Wei Liang

Abstract: Compact steep-spectrum sources (CSSs) likely represent a population of young radio-loud active galactic nuclei (AGNs) and have been identified as gamma-ray emitting sources. We present a comprehensive analysis of their gamma-ray emission observed with Fermi/LAT and establish their broadband spectral energy distributions (SEDs). We derive their jet properties by the SED fits with a two-zone leptoni… ▽ More Compact steep-spectrum sources (CSSs) likely represent a population of young radio-loud active galactic nuclei (AGNs) and have been identified as gamma-ray emitting sources. We present a comprehensive analysis of their gamma-ray emission observed with Fermi/LAT and establish their broadband spectral energy distributions (SEDs). We derive their jet properties by the SED fits with a two-zone leptonic model for radiations from the compact core and large-scale extended region, and explore the possible signature of a unification picture of jet radiation among subclasses of AGNs. We show that the observed gamma-rays of CSSs with significant variability are contributed by the radiation of their compact cores via the inverse Compton process of the torus photons. The derived power-law distribution index of the radiating electrons is p_1~1.5-1.8, magnetic field strength is B~0.15-0.6 G, and Doppler boosting factor is δ~2.8-8.9. Assuming that the jet is composed of electron-positron pairs, the compact cores of CSSs are magnetized and have a high radiation efficiency, similar to that of flat spectrum radio quasars. The six CSSs on average have higher Eddington ratio and black hole mass than those non-GeV-detected CSSs, and they follow the correlation between the jet power in units of Eddington luminosity (P_jet/L_Edd) and Eddington ratio (R_Edd) with other sub-classes of AGNs, P_jet/L_Edd~R_Edd^0.52, indicating that R_Edd would be a key physical driver for the unification scheme of AGN jet radiation. △ Less

Submitted 2 July, 2020; v1 submitted 23 May, 2020; originally announced May 2020.

Comments: 37 Pages, 11 Figures, 4 Tables, accepted for publication in ApJ

The gamma-ray and Optical Variability Analysis of the BL Lac Object 3FGL J0449.4-4350

Authors: Xing Yang, Tingfeng Yi, Yan Zhang, Huaizhen Li, Lisheng Mao, Haiming Zhang, Li Ma

Abstract: We have assembled the historical light curves of the BL Lac Object 3FGL J0449.4-4350 at optical and gamma-ray bands, the time spanning about 10 years, analyzed the periodic variability of the light curves by using four different methods (Lomb-Scargle periodogram, REDFIT38, Jurkevich and DACF). We detected a marginally possible quasi-periodic oscillation (QPO) of ~450 days. Assuming it originates f… ▽ More We have assembled the historical light curves of the BL Lac Object 3FGL J0449.4-4350 at optical and gamma-ray bands, the time spanning about 10 years, analyzed the periodic variability of the light curves by using four different methods (Lomb-Scargle periodogram, REDFIT38, Jurkevich and DACF). We detected a marginally possible quasi-periodic oscillation (QPO) of ~450 days. Assuming it originates from the helical motion jet in a supermassive binary black hole (SMBBH) system undergoing major merger, we estimate the primary black hole mass M~7.7*10^{9} M_sun. To explore the origin of the gamma-ray, we investigated the optical-gamma-ray correlations using discrete correlation function (DCF) method, and found that the correlation between the two bands is very significant. This strong correlation tends to imply lepton self-synchro-Compton (LSSC) model to produce the gamma-ray. △ Less

Submitted 25 March, 2020; originally announced March 2020.

Comments: 19 pages, 6 figures, 2 tables, Accepted for publication in PASP

Journal ref: PASP,132:044101(2020)

Diversity of gamma-ray and Radio Variabilities of Bright Blazars and Implications for gamma-ray Emission Location

Authors: Hai-Ming Zhang, Zhen-Jie Wang, Jin Zhang, Ting-Feng Yi, Liang Chen, Rui-Jing Lu, En-Wei Liang

Abstract: Violent multi-wavelength variabilities are observed in gamma-ray-selected blazars. We present an analysis of long-term light curves for eight bright blazars to explore the co-variation pattern in the gamma-ray and radio bands. We extract their gamma-ray light curves and spectra with data observed by the Fermi/LAT since 2008. We find diverse co-variation patterns between the gamma-ray and radio (at… ▽ More Violent multi-wavelength variabilities are observed in gamma-ray-selected blazars. We present an analysis of long-term light curves for eight bright blazars to explore the co-variation pattern in the gamma-ray and radio bands. We extract their gamma-ray light curves and spectra with data observed by the Fermi/LAT since 2008. We find diverse co-variation patterns between the gamma-ray and radio (at 43 GHz) fluxes in these sources. The gamma-ray and radio fluxes of 3C 454.3 and PKS 1633+382 are correlated without any time-lag, suggesting that they are from the same radiation region. Similar correlation is also observed in 3C 273 and PKS 1222+216, but the radio flux is lag behind the gamma-ray flux approximately ~160 days and ~290 days, respectively. This likely suggests that their gamma-ray emission regions are located at the upstream of their radio cores at 43 GHz. The gamma-ray and radio fluxes of the other four blazars are not correlated, implying that the gamma-ray and radio emission may be from different regions in their jets. The gamma-ray light curves of the eight blazars can be decomposed into some long timescale variability components and fast spike flares. We propose that they may be attributed to the central engine activity and the magnetic reconnection process or turbulence in the local emission region, respectively. △ Less

Submitted 24 March, 2020; originally announced March 2020.

Comments: 34 pages, 11 figures, 2 tables, accepted for publication in PASJ

Confirmed short periodic variability of subparsec supermassive binary black hole candidate Mrk 231

Authors: Andjelka B. Kovacevic, Tignfeng Yi, Xinyu Dai, Xing Yang, Iva Cvorovic-Hajdinjak, Luka C. Popovic

Abstract: Here we confirm the short periodic variability of a subparsec supermassive binary black hole (SMBBH) candidate Mrk 231 in the extended optical photometric data set collected by the Catalina Real-Time Transient Survey (CRTS) and All-Sky Automated Survey for Supernovae (ASAS-SN). Using the Lomb-Scargle periodogram and 2DHybrid method, we detected the significant periodicity of ~ 1.1 yr beyond a damp… ▽ More Here we confirm the short periodic variability of a subparsec supermassive binary black hole (SMBBH) candidate Mrk 231 in the extended optical photometric data set collected by the Catalina Real-Time Transient Survey (CRTS) and All-Sky Automated Survey for Supernovae (ASAS-SN). Using the Lomb-Scargle periodogram and 2DHybrid method, we detected the significant periodicity of ~ 1.1 yr beyond a damped random walk model in the CRTS+ASAS-SN optical data set. Mrk 231 has been previously proposed as a SMBBH candidate with a highly unequal mass ratio (q~ 0.03), very tight mutual separation of ~590 AU, and an orbital period of ~1.2 yr. Hence, our result further supports, even though not prove, the intriguing hypothesis that SMBBHs with low mass ratios may be more common than close-equal mass SMBBHs. This result, however, was obtained from the contribution of CRTS data with limited sampling cadence and photometric accuracy, and further monitoring of Mrk 231 is crucial to confirm the periodicity. △ Less

Submitted 13 March, 2020; originally announced March 2020.

Comments: Accepted for publication in MNRAS

A Neutron Star-White Dwarf Binary Model for Periodically Active Fast Radio Burst Sources

Authors: Wei-Min Gu, Tuan Yi, Tong Liu

Abstract: We propose a compact binary model with an eccentric orbit to explain periodically active fast radio burst (FRB) sources, where the system consists of a neutron star (NS) with strong dipolar magnetic fields and a magnetic white dwarf (WD). In our model, the WD fills its Roche lobe at periastron, and mass transfer occurs from the WD to the NS around this point. The accreted material may be fragmente… ▽ More We propose a compact binary model with an eccentric orbit to explain periodically active fast radio burst (FRB) sources, where the system consists of a neutron star (NS) with strong dipolar magnetic fields and a magnetic white dwarf (WD). In our model, the WD fills its Roche lobe at periastron, and mass transfer occurs from the WD to the NS around this point. The accreted material may be fragmented into a number of parts, which arrive at the NS at different times. The fragmented magnetized material may trigger magnetic reconnection near the NS surface. The electrons can be accelerated to an ultra-relativistic speed, and therefore the curvature radiation of the electrons can account for the burst activity. In this scenario, the duty cycle of burst activity is related to the orbital period of the binary. We show that such a model may work for duty cycles roughly from ten minutes to two days. For the recently reported 16.35-day periodicity of FRB 180916.J0158+65, our model does not naturally explain such a long duty cycle, since an extremely high eccentricity ($e>0.95$) is required. △ Less

Submitted 30 June, 2020; v1 submitted 24 February, 2020; originally announced February 2020.

Comments: 5 pages, 2 figures, accepted for publication in MNRAS

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

Authors: Tuan Yi, Mouyuan Sun, Wei-Min Gu

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

Submitted 2 October, 2019; originally announced October 2019.

Comments: 13 pages, 4 figures, accepted for publication in ApJ

Searching for Black Hole Candidates by LAMOST and ASAS-SN

Authors: Ling-Lin Zheng, Wei-Min Gu, Tuan Yi, Jin-Bo Fu, Hui-Jun Mu, Fan Yang, Song Wang, Zhong-Rui Bai, Hao Sou, Yu Bai, Yi-Ze Dong, Hao-Tong Zhang, Ya-Juan Lei, Junfeng Wang, Jianfeng Wu, Jifeng Liu

Abstract: Most dynamically confirmed stellar-mass black holes and the candidates were originally selected from X-ray outbursts. In the present work, we search for black hole candidates in the LAMOST survey by using the spectra along with photometry from the ASAS-SN survey, where the orbital period of the binary may be revealed by the periodic light curve, such as the ellipsoidal modulation type. Our sample… ▽ More Most dynamically confirmed stellar-mass black holes and the candidates were originally selected from X-ray outbursts. In the present work, we search for black hole candidates in the LAMOST survey by using the spectra along with photometry from the ASAS-SN survey, where the orbital period of the binary may be revealed by the periodic light curve, such as the ellipsoidal modulation type. Our sample consists of 9 binaries, where each source contains a giant star with large radial velocity variation ($ΔV_{\rm R} > 70~{\rm km~s^{-1}}$) and periods known from light curves. We focus on the 9 sources with long periods ($T_{\rm ph} > 5$ days) and evaluate the mass $M_2$ of the optically invisible companion. Since the observed $ΔV_{\rm R}$ from only a few repeating spectroscopic observations is a lower limit of the real amplitude, the real mass $M_2$ can be significantly higher than the current evaluation. It is likely an efficient method to place constraints on $M_2$ by combining $ΔV_{\rm R}$ from LAMOST and $T_{\rm ph}$ from ASAS-SN, particularly by the ongoing LAMOST Medium Resolution Survey. △ Less

Submitted 13 September, 2019; originally announced September 2019.

Comments: 12 pages, 3 figures, accepted for publication in Astronomical Journal

Thick-disc model to explain the spectral state transition in NGC 247

Authors: Jing Guo, Mouyuan Sun, Wei-Min Gu, Tuan Yi

Abstract: We propose the thick-disc model of Gu et al. 2016 to interpret the transition between soft ultraluminous state (SUL) and supersoft ultraluminous (SSUL) state in NGC 247. As accretion rate increases, the inner disc will puff up and act as shield to block the innermost X-ray emission regions and absorb both soft and hard X-ray photons. The absorbed X-ray emission will be re-radiated as a much softer… ▽ More We propose the thick-disc model of Gu et al. 2016 to interpret the transition between soft ultraluminous state (SUL) and supersoft ultraluminous (SSUL) state in NGC 247. As accretion rate increases, the inner disc will puff up and act as shield to block the innermost X-ray emission regions and absorb both soft and hard X-ray photons. The absorbed X-ray emission will be re-radiated as a much softer blackbody X-ray spectrum. Hence NGC 247 shows flux dips in the hard X-ray band and transits from the SUL state to the SSUL state. The $\sim 200$s transition timescale can be explained by the viscous timescale. According to our model, the inner disc in the super-soft state is thicker and has smaller viscous timescale than in the soft state. X-ray flux variability, which is assumed to be driven by accretion rate fluctuations, might be viscous time-scale invariant. Therefore, in the SSUL state, NGC 247 is more variable. The bolometric luminosity is saturated in the thick disc; the observed radius-temperature relation can therefore be naturally explained. △ Less

Submitted 25 February, 2019; originally announced February 2019.

Comments: 4 pages, 3 figures, Accepted for publication in MNRAS

A possible feedback mechanism of outflows from a black hole hyperaccretion disk in the center of jet-driven iPTF14hls

Authors: Tong Liu, Cui-Ying Song, Tuan Yi, Wei-Min Gu, Xiao-Feng Wang

Abstract: iPTF14hls is an unusually bright, long-lived II-P supernova (SN), whose light curve has at least five peaks. We propose that the outflows from the black hole hyperaccretion systems in the center of the collapsars should continuously inject into the envelope. For a jet-driven core-collapsar model, the outflow feedback results in prolonging the accretion timescale and fluctuating accretion rates in… ▽ More iPTF14hls is an unusually bright, long-lived II-P supernova (SN), whose light curve has at least five peaks. We propose that the outflows from the black hole hyperaccretion systems in the center of the collapsars should continuously inject into the envelope. For a jet-driven core-collapsar model, the outflow feedback results in prolonging the accretion timescale and fluctuating accretion rates in our analytic solutions. Thus, the long period of luminous, varying SN iPTF14hls might originate from the choked jets, which are regulated by the feedback of the strong disk outflows in a massive core-collapsar. One can expect that jet-driven iPTF14hls may last no more than approximately 3,000 days, and the luminosity may quickly decrease in the later stages. Moreover, the double-peak light curves in some SNe might be explained by the outflow feedback mechanism. △ Less

Submitted 9 February, 2019; originally announced February 2019.

Comments: 15 pages, 2 figures, accepted for publication in Journal of High Energy Astrophysics

A Method to Search for Black Hole Candidates with Giant Companions by LAMOST

Authors: Wei-Min Gu, Hui-Jun Mu, Jin-Bo Fu, Ling-Lin Zheng, Tuan Yi, Zhong-Rui Bai, Song Wang, Hao-Tong Zhang, Ya-Juan Lei, Yu Bai, Jianfeng Wu, Junfeng Wang, Jifeng Liu

Abstract: We propose a method to search for stellar-mass black hole (BH) candidates with giant companions from spectroscopic observations. Based on the stellar spectra of LAMOST Data Release 6, we obtain a sample of seven giants in binaries with large radial velocity variation $ΔV_R > 80~{\rm km~s^{-1}}$. With the effective temperature, surface gravity, and metallicity provided by LAMOST, and the parallax g… ▽ More We propose a method to search for stellar-mass black hole (BH) candidates with giant companions from spectroscopic observations. Based on the stellar spectra of LAMOST Data Release 6, we obtain a sample of seven giants in binaries with large radial velocity variation $ΔV_R > 80~{\rm km~s^{-1}}$. With the effective temperature, surface gravity, and metallicity provided by LAMOST, and the parallax given by {\it Gaia}, we can estimate the mass and radius of the giant, and therefore evaluate the possible mass of the optically invisible star in the binary. We show that the sources in our sample are potential BH candidates, and are worthy of dynamical measurement by further spectroscopic observations. Our method may be particularly valid for the selection of BH candidates in binaries with unknown orbital periods. △ Less

Submitted 7 February, 2019; originally announced February 2019.

Comments: 13 pages, 2 figures, accepted for publication in ApJ Letters

Compact binary merger and kilonova: outflows from remnant disc

Authors: Tuan Yi, Wei-Min Gu, Tong Liu, Rajiv Kumar, Hui-Jun Mu, Cui-Ying Song

Abstract: Outflows launched from a remnant disc of compact binary merger may have essential contribution to the kilonova emission. Numerical calculations are conducted in this work to study the structure of accretion flows and outflows. By the incorporation of limited-energy advection in the hyper-accretion discs, outflows occur naturally from accretion flows due to imbalance between the viscous heating and… ▽ More Outflows launched from a remnant disc of compact binary merger may have essential contribution to the kilonova emission. Numerical calculations are conducted in this work to study the structure of accretion flows and outflows. By the incorporation of limited-energy advection in the hyper-accretion discs, outflows occur naturally from accretion flows due to imbalance between the viscous heating and the sum of the advective and radiative cooling. Following this spirit, we revisit the properties of the merger outflow ejecta. Our results show that around $10^{-3} \sim 10^{-1} M_\odot$ of the disc mass can be launched as powerful outflows. The amount of unbound mass varies with the disc mass and the viscosity. The outflow-contributed peak luminosity is around $10^{40} \sim 10^{41}$ erg/s. Such a scenario can account for the observed kilonovae associated with short gamma-ray bursts, including the recent event AT2017gfo (GW170817). △ Less

Submitted 29 January, 2018; originally announced January 2018.

Comments: 8 pages, 6 figures, accepted for publication in MNRAS

The Radio and $γ$-ray Variability Analysis of S5 0716+714

Authors: H. Z. Li, Y. G. Jiang, T. F. Yi, D. F. Guo, X. Chen, H. M. Zhang, Q. G. Gao, F. W. Lu, J. Y. Ren

Abstract: We have studied the variability of \object{S5 0716+714} at radio 15 GHz and $γ$-ray band using three different methods. A possible periodicity of $P_{15 GHz}=266.0\pm11.5$ and $P_γ=344.0\pm16.4$ days are obtained for radio 15 GHz and $γ$-ray light curves, respectively. The variability may be related to the intrinsically emission mechanism. The difference between the variability timescales of radio… ▽ More We have studied the variability of \object{S5 0716+714} at radio 15 GHz and $γ$-ray band using three different methods. A possible periodicity of $P_{15 GHz}=266.0\pm11.5$ and $P_γ=344.0\pm16.4$ days are obtained for radio 15 GHz and $γ$-ray light curves, respectively. The variability may be related to the intrinsically emission mechanism. The difference between the variability timescales of radio 15 GHz and $γ$-ray may be due to that the emission of radio 15 GHz is produced via the synchrotron process, while the $γ$-ray is produced by both the SSC and EC processes. △ Less

Submitted 20 January, 2018; originally announced January 2018.

Comments: 8 pages, 5 figures, Accepted for publication in ApSS

Flux and Spectral Variation Characteristics of 3C 454.3 at the GeV Band

Authors: Hai-Ming Zhang, Jin Zhang, Rui-Jing Lu, Ting-Feng Yi, Xiao-Li Huang, En-Wei Liang

Abstract: We analyze the long-term lightcurve of 3C 454.3 observed with Fermi/LAT and investigate its relation to the flux in the radio, optical, and X-ray bands. By fitting the 1-day binned GeV lightcurve with multiple Gaussian function (MGF), we propose that the typical variability timescale in the GeV band is 1--10 days. The GeV flux variation is accompanied by the spectral variation characterized as flu… ▽ More We analyze the long-term lightcurve of 3C 454.3 observed with Fermi/LAT and investigate its relation to the flux in the radio, optical, and X-ray bands. By fitting the 1-day binned GeV lightcurve with multiple Gaussian function (MGF), we propose that the typical variability timescale in the GeV band is 1--10 days. The GeV flux variation is accompanied by the spectral variation characterized as flux-tracking, i.e., "harder when brighter". The GeV flux is correlated with the optical and X-ray fluxes, and a weak correlation between gamma-ray flux and radio flux is also observed. The gamma-ray flux is not correlated with the optical linear polarization degree for the global lightcurves, but they show a correlation for the lightcurves before MJD 56000. The power density spectrum of the global lightcurve shows an obvious turnover at ~7.7 days, which may indicate a typical variability timescale of 3C 454.3 in the gamma-ray band. This is also consistent with the derived timescales by fitting the global lightcurve with MGF. The spectral evolution and an increase of the optical linear polarization degree along with the increase of the gamma-ray flux may indicate that the radiation particles are accelerated and the magnetic field is ordered by the shock processes during the outbursts. In addition, the nature of 3C 454.3 may be consistent with the self-organized criticality system, similar to Sagittarius A, and thus the outbursts are from plasmoid ejections driven by magnetic reconnection. This may further support the idea that the jet radiation regions are magnetized. △ Less

Submitted 21 December, 2017; originally announced December 2017.

Comments: 29 pages, 9 figures, accepted for publication in RAA

Multiple Wavelength Variability and Quasi-Periodic Oscillation of PMN J0948+0022

Authors: Jin Zhang, Hai-Ming Zhang, Yong-Kai Zhu, Ting-Feng Yi, Su Yao, Rui-Jing Lu, En-Wei Liang

Abstract: We present a comprehensive analysis of multiple wavelength observational data of the first GeV-selected narrow-line Seyfert 1 galaxy PMN J0948+0022. We derive its lightcurves in the γ-ray and X-ray bands from the data observed with Fermi/LAT and Swift/XRT, and make the optical and radio lightcurves by collecting the data from the literature. These lightcurves show significant flux variations. With… ▽ More We present a comprehensive analysis of multiple wavelength observational data of the first GeV-selected narrow-line Seyfert 1 galaxy PMN J0948+0022. We derive its lightcurves in the γ-ray and X-ray bands from the data observed with Fermi/LAT and Swift/XRT, and make the optical and radio lightcurves by collecting the data from the literature. These lightcurves show significant flux variations. With the LAT data we show that this source is analogue to typical flat spectrum radio quasars in the $L_γ-Γ_γ$ plane, where $L_γ$ and $Γ_γ$ are the luminosity and spectral index in the LAT energy band. The gamma-ray flux is correlated with the V-band flux with a lag of ~44 days, and a moderate quasi-periodic oscillation (QPO) with a periodicity of ~490 days observed in the LAT lightcurve. A similar QPO signature is also found in the V-band lightcurve. The gamma-ray flux is not correlated with the radio flux in 15 GHz, and no similar QPO signature is found in a confidence level of 95%. Possible mechanisms of the QPO are discussed. We propose that gravitational wave observations in the future may clarify the current plausible models for the QPO. △ Less

Submitted 7 September, 2017; originally announced September 2017.

Comments: 18 pages, 7 figures, accepted for publication in ApJ

The correlations of multi-wave band luminosity and BLR luminosity in Fermi 2LAC blazars

Authors: Zerui Wang, Rui Xue, Zhaohua Xie, Leiming Du, Tingfeng Yi, Yunbing Xu, Wenguang Liu

Abstract: We collect a sample of 78 Fermi detected blazars with broad line region (BLR) data and the quasi-simultaneous multi-wave band data. By analyzing the sample, we find: (1) For whole blazar sample, there exist significant correlations between radio, $γ$-ray luminosity and BLR luminosity. The slope of the best-fit linear regression equation is close to the theoretical value 1. Our results provided ano… ▽ More We collect a sample of 78 Fermi detected blazars with broad line region (BLR) data and the quasi-simultaneous multi-wave band data. By analyzing the sample, we find: (1) For whole blazar sample, there exist significant correlations between radio, $γ$-ray luminosity and BLR luminosity. The slope of the best-fit linear regression equation is close to the theoretical value 1. Our results provided another support of jet-disk symbiosis. (2) For FSRQ sample, we find a significant correlation between radio luminosity and optical luminosity, and weaker correlations between X-ray luminosity and radio, optical luminosity. It indicates that the X-ray band radiation mechanism is slightly different from radio and optical bands. The X-ray emission is diluted by inverse Compton scattering (IC). There is a significant correlation between BLR luminosity and $γ$-ray luminosity, but no correlations between X-ray luminosity and BLR, $γ$-ray luminosity, which suggests that the seed soft photons are different for IC process. From these results we suggest that the X-ray band emission is diluted by synchrotron self-Compton (SSC) process, $γ$-ray band emission is produced by external Compton (EC) process. △ Less

Submitted 28 August, 2017; originally announced August 2017.

Comments: 10 pages, 1 figure, accepted for publication in Ap&SS

On estimators of the jet bolometric luminosity of Fermi 2LAC blazars

Authors: Zerui Wang, Rui Xue, Leiming Du, Zhaohua Xie, Dingrong Xiong, Tingfeng Yi, Yunbing Xu, Wenguang Liu

Abstract: Bolometric luminosity is a basic physical parameter that is widely used in the study of blazars. Due to the lack of simultaneous full wavelength data, several estimators of the bolometric luminosity are being used in practice. In this paper, we study and evaluate the reliability and significance of six estimators, the $5GHz$ luminosity, the 1keV luminosity, the $γ$-ray luminosity, the $5GHz$ lumin… ▽ More Bolometric luminosity is a basic physical parameter that is widely used in the study of blazars. Due to the lack of simultaneous full wavelength data, several estimators of the bolometric luminosity are being used in practice. In this paper, we study and evaluate the reliability and significance of six estimators, the $5GHz$ luminosity, the 1keV luminosity, the $γ$-ray luminosity, the $5GHz$ luminosity+the 1keV luminosity, the $5GHz$ luminosity+the $γ$-ray luminosity and the 1keV luminosity+the $γ$-ray luminosity, by analyzing the linear correlations between the integrated bolometric luminosity and them. Our main results are as follows. (i) All the six estimators are reliable in the sense that they are all significant correlated with the bolometric luminosity. (ii) Ranking from the higher significance of the reliability to lower one the six estimators are the the $5GHz$ luminosity+the $γ$-ray luminosity, the 1keV luminosity+the $γ$-ray luminosity, $γ$-ray luminosity, the $5GHz$ luminosity+the 1keV luminosity, the $5GHz$ luminosity and the 1keV luminosity. (iii) We suggest that the bolometric luminosity can be well estimated by the $γ$-ray luminosity using the best linear equation that given in this paper for Fermi FSRQs. (iv) According to the linear regressions obtained in the analysis, we provide calibration for each estimator. △ Less

Submitted 20 July, 2017; originally announced July 2017.

Comments: 10 pages, 4 figures, accepted for publication in Ap&SS

Chemical Enhancements in Shock-accelerated Particles: Ab-initio Simulations

Authors: Damiano Caprioli, Dennis T. Yi, Anatoly Spitkovsky

Abstract: We study the thermalization, injection, and acceleration of ions with different mass/charge ratios, $A/Z$, in non-relativistic collisionless shocks via hybrid (kinetic ions-fluid electrons) simulations. In general, ions thermalize to a post-shock temperature proportional to $A$. When diffusive shock acceleration is efficient, ions develop a non-thermal tail whose extent scales with $Z$ and whose n… ▽ More We study the thermalization, injection, and acceleration of ions with different mass/charge ratios, $A/Z$, in non-relativistic collisionless shocks via hybrid (kinetic ions-fluid electrons) simulations. In general, ions thermalize to a post-shock temperature proportional to $A$. When diffusive shock acceleration is efficient, ions develop a non-thermal tail whose extent scales with $Z$ and whose normalization is enhanced as $(A/Z)^2$, so that incompletely-ionized heavy ions are preferentially accelerated. We discuss how these findings can explain observed heavy-ion enhancements in Galactic cosmic rays. △ Less

Submitted 26 April, 2017; originally announced April 2017.

Comments: 6 pages, 4 figures, submitted to PRL

Journal ref: Phys. Rev. Lett. 119, 171101 (2017)

Multi-color optical monitoring of the quasar 3C 273 from 2005 to 2016

Authors: Dingrong Xiong, Jinming Bai, Haojing Zhang, Junhui Fan, Minfeng Gu, Tingfeng Yi, Xiong Zhang

Abstract: We have monitored the quasar 3C 273 in optical $V$, $R$ and $I$ bands from 2005 to 2016. Intraday variability (IDV) is detected on seven nights. The variability amplitudes for most of nights are less than 10\% and four nights more than 20\%. When considering the nights with time spans $>4$ hours, the value of duty cycle (DC) is 14.17 per cent. Over the twelve years, the overall magnitude and color… ▽ More We have monitored the quasar 3C 273 in optical $V$, $R$ and $I$ bands from 2005 to 2016. Intraday variability (IDV) is detected on seven nights. The variability amplitudes for most of nights are less than 10\% and four nights more than 20\%. When considering the nights with time spans $>4$ hours, the value of duty cycle (DC) is 14.17 per cent. Over the twelve years, the overall magnitude and color index variabilities are $\bigtriangleup I=0^{\rm m}.67$, $\bigtriangleup R=0^{\rm m}.72$, $\bigtriangleup V=0^{\rm m}.68$, and $\bigtriangleup (V-R)=0^{\rm m}.25$ respectively. The largest clear IDV has an amplitude of 42% over just 5.8 minutes and the weakest detected IDV is 5.4% over 175 minutes. The BWB (bluer when brighter) chromatic trend is dominant for 3C 273 and appears at different flux levels on intraday timescales. The BWB trend exists for short-term timescales and intermediate-term timescales but different timescales have different correlations. There is no BWB trend for our whole time-series data sets. A significant anti-correlation between BWB trend and length of timescales is found. Combining with $V$-band data from previous works, we find a possible quasi-periodicity of $P=3918\pm1112$ days. The possible explanations for the observed variability, BWB chromatic trend and periodicity are discussed. △ Less

Submitted 6 March, 2017; v1 submitted 5 March, 2017; originally announced March 2017.

Comments: 63 pages, 11 figures, 6 tables. Accepted for publication in ApJS

Evaluating Optical Classification for {\em Fermi} Blazar Candidates with a Statistical method using Broadband Spectral Indices

Authors: Ting-Feng Yi, Jin Zhang, Rui-Jing Lu, Rui Huang, En-Wei Liang

Abstract: We aim to test if a blazar candidate of uncertain-type (BCU) in the third Fermi active galactic nuclei catalog (3LAC) can be potentially classified as a BL Lac object or a flat spectrum radio quasar (FSRQ) by performing a statistical analysis of its broadband spectral properties. We find that 34% of the radio-selected BCUs (583 BCUs) are BL Lac-like and 20% of them are FSRQ-like, which maybe withi… ▽ More We aim to test if a blazar candidate of uncertain-type (BCU) in the third Fermi active galactic nuclei catalog (3LAC) can be potentially classified as a BL Lac object or a flat spectrum radio quasar (FSRQ) by performing a statistical analysis of its broadband spectral properties. We find that 34% of the radio-selected BCUs (583 BCUs) are BL Lac-like and 20% of them are FSRQ-like, which maybe within 90% level of confidence. Similarly, 77.3% of the X-ray selected BCUs (176 BCUs) are evaluated as BL Lac-like and 6.8% of them may be FSRQ-like sources. And 88.7% of the BL Lac-like BCUs that have synchrotron peak frequencies available are high synchrotron peaked BL Lacs in the X-ray selected BCUs. The percentages are accordingly 62% and 7.3% in the sample of 124 optical-selected BCUs. The high ratio of source numbers of the BL Lac-like to the FSRQ-like BCUs in the X-ray and optically selected BCU samples is due to the selection effect. Examining the consistency between our evaluation and spectroscopic identification case by case with a sample of 78 radio-selected BCUs, it is found that the statistical analysis and its resulting classifications agree with the results of the optical follow-up spectroscopic observations. Our observation campaign for high-|ρ_{s}| BCUs , i.e., |ρ_{s}|>0.8, selected with our method is ongoing. △ Less

Submitted 27 February, 2017; originally announced February 2017.

Comments: 38 pages, 6 figures, 3 tables, accepted for publication in ApJ

Vertical Advection Effects on Hyper-accretion Disks and Potential Link between Gamma-ray Bursts and Kilonovae

Authors: Tuan Yi, Wei-Min Gu, Feng Yuan, Tong Liu, Hui-Jun Mu

Abstract: Recent simulations on super-Eddington accretion flows have shown that, apart from the diffusion process, the vertical advection based on magnetic buoyancy can be a more efficient process to release the trapped photons in the optically thick disk. As a consequence, the radiative luminosity from the accretion disk can be far beyond the Eddington value. Following this spirit, we revisit the structure… ▽ More Recent simulations on super-Eddington accretion flows have shown that, apart from the diffusion process, the vertical advection based on magnetic buoyancy can be a more efficient process to release the trapped photons in the optically thick disk. As a consequence, the radiative luminosity from the accretion disk can be far beyond the Eddington value. Following this spirit, we revisit the structure and radiation of hyper-accretion disks with mass accretion rates in the range $10^{-3}\sim 10~M_{\sun}~{\rm s}^{-1}$. Our results show that, due to the strong cooling through the vertical advection, the disk temperature becomes lower than that in the classic model without the vertical advection process, and therefore the neutrino luminosity from the disk is lower. On the other hand, the gamma-ray photons released through the vertical advection can be extremely super-Eddington. We argue that the large amount of escaped gamma-ray photons may have more significant contribution to the primordial fireball than the neutrino annihilation, and may hint a link between gamma-ray bursts and kilonovae in the black hole hyper-accretion scenario. △ Less

Submitted 25 January, 2017; originally announced January 2017.

Comments: 13 pages, 7 figures, accepted for publication in ApJ

Central Engine of Late-Time X-ray Flares with Internal Origin

Authors: Hui-Jun Mu, Wei-Min Gu, Shu-Jin Hou, Tong Liu, Da-Bin Lin, Tuan Yi, En-Wei Liang, Ju-Fu Lu

Abstract: This work focuses on a sample of seven extremely late-time X-ray flares with peak time $t_{\rm p} > 10^4 {\rm s}$, among which two flares can be confirmed as the late-time activity of central engine. The main purpose is to investigate the mechanism of such late-time flares based on the internal origin assumption. In the hyper-accreting black hole (BH) scenario, we study the possibility of two well… ▽ More This work focuses on a sample of seven extremely late-time X-ray flares with peak time $t_{\rm p} > 10^4 {\rm s}$, among which two flares can be confirmed as the late-time activity of central engine. The main purpose is to investigate the mechanism of such late-time flares based on the internal origin assumption. In the hyper-accreting black hole (BH) scenario, we study the possibility of two well-known mechanisms as the central engine to power such X-ray flares, i.e., the neutrino-antineutrino annihilation and the Blandford-Znajek (BZ) process. Our results show that the annihilation luminosity is far below the observational data. Thus, the annihilation mechanism cannot account for such late-time flares. For the BZ process, if the role of outflows is taken into consideration, the inflow mass rate near the horizon will be quite low such that the magnetic field will probably be too weak to power the observed X-ray flares. We therefore argue that, for the late-time flares with internal origin, the central engine is unlikely to be associated with BHs. On the contrary, a fast rotating neutron star with strong bipolar magnetic fields may be responsible for such flares. △ Less

Submitted 27 September, 2016; originally announced September 2016.

Comments: 16 pages, 6 figures, accepted for publication in ApJ