Showing 1–50 of 291 results for author: Lu, Y

A slightly oblate dark matter halo revealed by a retrograde precessing Galactic disk warp

Authors: Yang Huang, Qikang Feng, Tigran Khachaturyants, Huawei Zhang, Jifeng Liu, Juntai Shen, Timothy C. Beers, Youjun Lu, Song Wang, Haibo Yuan

Abstract: The shape of the dark matter (DM) halo is key to understanding the hierarchical formation of the Galaxy. Despite extensive efforts in recent decades, however, its shape remains a matter of debate, with suggestions ranging from strongly oblate to prolate. Here, we present a new constraint on its present shape by directly measuring the evolution of the Galactic disk warp with time, as traced by accu… ▽ More The shape of the dark matter (DM) halo is key to understanding the hierarchical formation of the Galaxy. Despite extensive efforts in recent decades, however, its shape remains a matter of debate, with suggestions ranging from strongly oblate to prolate. Here, we present a new constraint on its present shape by directly measuring the evolution of the Galactic disk warp with time, as traced by accurate distance estimates and precise age determinations for about 2,600 classical Cepheids. We show that the Galactic warp is mildly precessing in a retrograde direction at a rate of $ωおめが= -2.1 \pm 0.5 ({\rm statistical}) \pm 0.6 ({\rm systematic})$ km s$^{-1}$ kpc$^{-1}$ for the outer disk over the Galactocentric radius [$7.5, 25$] kpc, decreasing with radius. This constrains the shape of the DM halo to be slightly oblate with a flattening (minor axis to major axis ratio) in the range $0.84 \le q_Φふぁい \le 0.96$. Given the young nature of the disk warp traced by Cepheids (less than 200 Myr), our approach directly measures the shape of the present-day DM halo. This measurement, combined with other measurements from older tracers, could provide vital constraints on the evolution of the DM halo and the assembly history of the Galaxy. △ Less

Submitted 29 June, 2024; originally announced July 2024.

Comments: Published in Nature Astronomy on June 27th, 2024. Final published version here: https://www.nature.com/articles/s41550-024-02309-5

A high-velocity star recently ejected by an intermediate-mass black hole in M15

Authors: Yang Huang, Qingzheng Li, Jifeng Liu, Xiaobo Dong, Huawei Zhang, Youjun Lu, Cuihua Du

Abstract: The existence of intermediate-mass black holes (IMBHs) is crucial for understanding various astrophysical phenomena, yet their existence remains elusive, except for the LIGO-Virgo detection. We report the discovery of a high-velocity star J0731+3717, whose backward trajectory about 21 Myr ago intersects that of globular cluster M15 within the cluster tidal radius. Both its metallicity [Fe/H] and i… ▽ More The existence of intermediate-mass black holes (IMBHs) is crucial for understanding various astrophysical phenomena, yet their existence remains elusive, except for the LIGO-Virgo detection. We report the discovery of a high-velocity star J0731+3717, whose backward trajectory about 21 Myr ago intersects that of globular cluster M15 within the cluster tidal radius. Both its metallicity [Fe/H] and its alpha-to-iron abundance ratio [$αあるふぁ$/Fe] are consistent with those of M15. Furthermore, its location falls right on the fiducial sequence of the cluster M15 on the color-absolute magnitude diagram, suggesting similar ages. These support that J0731+3717 is originally associated with M15 at a confidence level of $5.4σしぐま$. We find that such a high-velocity star ($V_{\rm ej} = 548^{+6}_{-5}$ km s$^{-1}$) was most likely tidally ejected from as close as one astronomical unit to the center of M15, confirming an IMBH ($\ge 100 M_{\odot}$ with a credibility of 98%) as the exclusive nature of the central unseen mass proposed previously. △ Less

Submitted 2 June, 2024; originally announced June 2024.

Comments: 17 pages, 7 figures and 3 tables, submitted for review

Discovery of a long thermonuclear X-ray burst from the ultra-compact binary 4U 1850$-$087

Authors: Yongqi Lu, Zhaosheng Li, Wenhui Yu, Yuanyue Pan, Maurizio Falanga

Abstract: We report the detection of a long X-ray burst triggered on MJD 60171.65 from the ultra-compact binary 4U 1850$-$087 by the Monitor of All-sky X-ray Image and Neutron Star Interior Composition Explorer (NICER). We analyse the NICER data observed in between MJD 60095.19$-$60177.43, including one observation covered part of the long X-ray burst tail, i.e., $0.15-3.8$ hr after the trigger. The persist… ▽ More We report the detection of a long X-ray burst triggered on MJD 60171.65 from the ultra-compact binary 4U 1850$-$087 by the Monitor of All-sky X-ray Image and Neutron Star Interior Composition Explorer (NICER). We analyse the NICER data observed in between MJD 60095.19$-$60177.43, including one observation covered part of the long X-ray burst tail, i.e., $0.15-3.8$ hr after the trigger. The persistent spectra are quite similar and well described by a combination of multi-color disk blackbody, with the inner temperature of 0.5 keV, and a thermally comptonized continuum with the asymptotic power-law photon index of $Γがんま\sim2.2$, and electron temperature of $kT_{\rm e}\sim20-30$ keV. The persistent fluxes were around $3.8\times10^{-10}~{\rm erg~cm^{-2}~s^{-1}}$, corresponding to a local accretion rate of $1\%~\dot{m}_{\rm Edd}$. Part of time-resolved burst spectra show a clear deviation from the blackbody model, which can be improved by considering the enhanced persistent emission due to the Poynting-Robertson drag, or the reflected disk emission illuminated by the burst. From the burst flux during the cooling tail, we estimate the burst duration, $τたう\approx 0.78$ hr, the burst fluence, $E_\mathrm{b} \approx 4.1 \times 10^{41}$ ergs, and the ignition column depth, $y_{\rm ign}\approx 3.5\times10^{10}~{\rm g~cm^{-2}}$. We propose that the long X-ray burst is powered by unstable burning of pure helium in deep layer. Moreover, we identify significant 1 keV emission lines in the burst spectra, which may originate from the surrounding disk. △ Less

Submitted 22 April, 2024; originally announced April 2024.

Comments: Submitted to ApJ, revised version after the first referee report, 9 pages, 6 figures

On Variation of Light Curves and Broad Emission Lines for Periodic QSOs from co-rotating Supermassive binary black holes in elliptical orbits

Authors: Junqiang Ge, Youjun Lu, Changshuo Yan, Jifeng Liu

Abstract: Context. Periodic QSOs are considered as candidates of supermassive binary black hole (BBH) systems in galactic centers. Further confirmation of these candidates may require different lines of observational evidences. Aims. Assuming the Doopler boosting scenario, in this paper we investigate the (coherent) variations of both broad emission lines (BELs) and continuum light curves for active BBH sys… ▽ More Context. Periodic QSOs are considered as candidates of supermassive binary black hole (BBH) systems in galactic centers. Further confirmation of these candidates may require different lines of observational evidences. Aims. Assuming the Doopler boosting scenario, in this paper we investigate the (coherent) variations of both broad emission lines (BELs) and continuum light curves for active BBH systems surrounding by a circumbinary broad line region (cBLR) and focus on their dependence on the eccentric orbital configuration. Methods. We calculate the variation of continuum light and the Doppler enhanced/weakened photoionization of each BLR cloud according to the motion of BBHs in elliptical orbits, and finally obtain the coherent variation of the continuum and BELs. Results. We find that both the amplitude and variation pattern of the continuum light curves and the evolution of the BEL profiles sensitively depend on the eccentric orbital configuration of BBH systems. If only the secondary BH is active, both the variation amplitudes of continuum light curves and BELs increase with increasing BBH inclination angles and orbital eccentricities, but decrease with increasing BBH mass ratio. If both BHs are active, the asymmetry in the ionization of BLR clouds at different areas caused by the Doppler boosting effect of the secondary BH is weakened due to that of the primary BH at the opposite direction, which leads to systematically smaller variation amplitudes of both continuum light curves and BELs than the cases with only secondary BH activated. Conclusions. The coherent variations of the BEL profiles with the continuum light for those periodic QSOs provide an important way to confirm the existence of BBHs in their center. Future joint analysis of the light curves and multi-epoch observed BEL profiles for periodic QSOs may lead to the identification of a number of BBH systems. △ Less

Submitted 8 April, 2024; originally announced April 2024.

Comments: 13 pages, 9 figures, accepted for publication in A&A

Journal ref: A&A 687, A57 (2024)

Direct collapse supermassive black holes from relic particle decay

Authors: Yifan Lu, Zachary S. C. Picker, Alexander Kusenko

Abstract: We investigate the formation of high-redshift supermassive black holes (SMBHs) via the direct collapse of baryonic clouds, where the necessary formation of molecular hydrogen is supressed by a Lyman-Werner (LW) photon background from relic particle decay. We improve on existing studies by dynamically simulating the collapse, accounting for the adiabatic contraction of the DM halo, as well as the i… ▽ More We investigate the formation of high-redshift supermassive black holes (SMBHs) via the direct collapse of baryonic clouds, where the necessary formation of molecular hydrogen is supressed by a Lyman-Werner (LW) photon background from relic particle decay. We improve on existing studies by dynamically simulating the collapse, accounting for the adiabatic contraction of the DM halo, as well as the in-situ production of the LW photons within the cloud which reduce the impact of the cloud's shielding. We find a viable parameter space where the decay of either some or all of the dark matter could successfully lead to the direct collapse of a SMBH. △ Less

Submitted 5 April, 2024; originally announced April 2024.

Comments: 8 pages, 4 figures

Report number: IPMU24-0013

Impact of Starlink constellation on Early LSST: a Photometric Analysis of Satellite Trails with BRDF Model

Authors: Yao Lu

Abstract: We report a simulation and quantification of the impact of the Starlink constellation on LSST in terms of the trail surface brightness using a BRDF-based satellite photometric model. A total of 11,908 satellites from the Gen1 and Gen2A constellations are used to focus on the interference to the initial phase of LSST operation. The all-sky simulation shows that approximately 69.33% of the visible s… ▽ More We report a simulation and quantification of the impact of the Starlink constellation on LSST in terms of the trail surface brightness using a BRDF-based satellite photometric model. A total of 11,908 satellites from the Gen1 and Gen2A constellations are used to focus on the interference to the initial phase of LSST operation. The all-sky simulation shows that approximately 69.33% of the visible satellites over station have an apparent brightness greater than 7 mag with a v1.5 satellite model. The impact of satellite streaks exhibit a non-monotonic relationship to the solar altitude, with the worst moments occurring around $-15^{\circ}$ solar altitude. The assessment based on simulated schedules indicates that no trails can reach the saturation-level magnitude, but 71.61% trails show a surface brightness brighter than the best-case crosstalk correctable limits, and this percentage increases as the dodging weight increases. Therefore, avoiding satellites in the scheduler algorithm is an effective mitigation method, but both the number of streaks and their brightness should be taken into account simultaneously. △ Less

Submitted 17 March, 2024; originally announced March 2024.

Comments: 7 pages, 4 figures

CSST large-scale structure analysis pipeline: I. constructing reference mock galaxy redshift surveys

Authors: Yizhou Gu, Xiaohu Yang, Jiaxin Han, Yirong Wang, Qingyang Li, Zhenlin Tan, Wenkang Jiang, Yaru Wang, Jiaqi Wang, Antonios Katsianis, Xiaoju Xu, Haojie Xu, Wensheng Hong, Houjun Mo, Run Wen, Xianzhong Zheng, Feng Shi, Pengjie Zhang, Zhongxu Zhai, Chengze Liu, Wenting Wang, Ying Zu, Hong Guo, Youcai Zhang, Yi Lu , et al. (7 additional authors not shown)

Abstract: In this paper, we set out to construct a set of reference mock galaxy redshift surveys (MGRSs) for the future Chinese Space-station Survey Telescope (CSST) observation, where subsequent survey selection effects can be added and evaluated. This set of MGRSs is generated using the dark matter subhalos extracted from a high-resolution Jiutian $N$-body simulation of the standard $Λらむだ$CDM cosmogony with… ▽ More In this paper, we set out to construct a set of reference mock galaxy redshift surveys (MGRSs) for the future Chinese Space-station Survey Telescope (CSST) observation, where subsequent survey selection effects can be added and evaluated. This set of MGRSs is generated using the dark matter subhalos extracted from a high-resolution Jiutian $N$-body simulation of the standard $Λらむだ$CDM cosmogony with $Ωおめが_m=0.3111$, $Ωおめが_Λらむだ=0.6889$, and $σしぐま_8=0.8102$. The simulation has a boxsize of $1~h^{-1} {\rm Gpc}$, and consists of $6144^3$ particles with mass resolution $3.723 \times 10^{8} h^{-1} M_\odot $. In order to take into account the effect of redshift evolution, we first use all 128 snapshots in the Jiutian simulation to generate a light-cone halo/subhalo catalog. Next, galaxy luminosities are assigned to the main and subhalo populations using the subhalo abundance matching (SHAM) method with the DESI $z$-band luminosity functions at different redshifts. Multi-band photometries, as well as images, are then assigned to each mock galaxy using a 3-dimensional parameter space nearest neighbor sampling of the DESI LS observational galaxies and groups. Finally, the CSST and DESI LS survey geometry and magnitude limit cuts are applied to generate the required MGRSs. As we have checked, this set of MGRSs can generally reproduce the observed galaxy luminosity/mass functions within 0.1 dex for galaxies with $L > 10^8 L_\odot$ (or $M_* > 10^{8.5} M_\odot$) and within 1-$σしぐま$ level for galaxies with $L < 10^8L_\odot$ (or $M_* < 10^{8.5} M_\odot$). Together with the CSST slitless spectra and redshifts for our DESI LS seed galaxies that are under construction, we will set out to test various slitless observational selection effects in subsequent probes. △ Less

Submitted 15 March, 2024; originally announced March 2024.

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

New constraints on Triton's atmosphere from the 6 October 2022 stellar occultation

Authors: Ye Yuan, Chen Zhang, Fan Li, Jian Chen, Yanning Fu, Chunhai Bai, Xing Gao, Yong Wang, Tuhong Zhong, Yixing Gao, Liang Wang, Donghua Chen, Yixing Zhang, Yang Zhang, Wenpeng Xie, Shupi Zhang, Ding Liu, Jun Cao, Xiangdong Yin, Xiaojun Mo, Jing Liu, Xinru Han, Tong Liu, Yuqiang Chen, Zhendong Gao , et al. (25 additional authors not shown)

Abstract: The atmosphere of Triton was probed directly by observing a ground-based stellar occultation on 6 October 2022. This rare event yielded 23 positive light curves collected from 13 separate observation stations contributing to our campaign. The significance of this event lies in its potential to directly validate the modest pressure fluctuation on Triton, a phenomenon not definitively verified by pr… ▽ More The atmosphere of Triton was probed directly by observing a ground-based stellar occultation on 6 October 2022. This rare event yielded 23 positive light curves collected from 13 separate observation stations contributing to our campaign. The significance of this event lies in its potential to directly validate the modest pressure fluctuation on Triton, a phenomenon not definitively verified by previous observations, including only five stellar occultations, and the Voyager 2 radio occultation in 1989. Using an approach consistent with a comparable study, we precisely determined a surface pressure of $14.07_{-0.13}^{+0.21}~\mathrm{μみゅーbar}$ in 2022. This new pressure rules out any significant monotonic variation in pressure between 2017 and 2022 through direct observations, as it is in alignment with the 2017 value. Additionally, both the pressures in 2017 and 2022 align with the 1989 value. This provides further support for the conclusion drawn from the previous volatile transport model simulation, which is consistent with the observed alignment between the pressures in 1989 and 2017; that is to say, the pressure fluctuation is modest. Moreover, this conclusion suggests the existence of a northern polar cap extended down to at least $45^\circ$N$-60^\circ$N and the presence of nitrogen between $30^\circ$S and $0^\circ$. △ Less

Submitted 24 March, 2024; v1 submitted 14 March, 2024; originally announced March 2024.

Comments: Astronomy & Astrophysics, in press. 9 pages, 2 figures, 3 tables

Journal ref: A&A 684, L13 (2024)

BRDF-Based Photometric Modeling of LEO Constellation Satellite from Massive Observations

Authors: Yao Lu

Abstract: Modeling the brightness of satellites in large Low-Earth Orbit (LEO) constellations can not only assist the astronomical community in assessing the impact of reflected light from satellites, optimizing observing schedules and guiding data processing, but also motivate satellite operators to improve their satellite designs, thus facilitating cooperation and consensus among different stakeholders. T… ▽ More Modeling the brightness of satellites in large Low-Earth Orbit (LEO) constellations can not only assist the astronomical community in assessing the impact of reflected light from satellites, optimizing observing schedules and guiding data processing, but also motivate satellite operators to improve their satellite designs, thus facilitating cooperation and consensus among different stakeholders. This work presents a photometric model of the Starlink satellites based on the Bidirectional Reflectance Distribution Function (BRDF) using millions of photometric observations. To enhance model accuracy and computational efficiency, data filtering and reduction are employed, and chassis blocking on the solar array and the earthshine effect are taken into account. The assumptions of the model are also validated by showing that the satellite attitude is as expected, the solar array is nearly perpendicular to the chassis, and both the solar array pseudo-specular reflection and the chassis earthshine should be included in the model. Reflectance characteristics of the satellites and the apparent magnitude distributions over station are finally discussed based on the photometric predictions from the model. In addition to assessing the light pollution and guiding the development of response measures, accurate photometric models of satellites can also play an important role in areas such as space situational awareness. △ Less

Submitted 15 May, 2024; v1 submitted 9 March, 2024; originally announced March 2024.

Comments: 16 pages, 9 figures

Distribution Properties of the 6.7 GHz Methanol Masers and Their Surrounding Gases in the Milky Way

Authors: Tian Yang, Xi Chen, Yan-Kun Zhang, Xu-Jia Ouyang, Shi-Min Song, Jia-Liang Chen, Ying Lu

Abstract: An updated catalog consisting of 1092 6.7-GHz methanol maser sources was reported in this work. Additionally, the NH3 (1, 1), NH3 (2, 2), and NH3 (3, 3) transitions were observed towards 214 star forming regions using the Shanghai Tianma radio telescope (TMRT) in order to examine the differences in physical environments, such as excitation temperature and column density of molecular clouds associa… ▽ More An updated catalog consisting of 1092 6.7-GHz methanol maser sources was reported in this work. Additionally, the NH3 (1, 1), NH3 (2, 2), and NH3 (3, 3) transitions were observed towards 214 star forming regions using the Shanghai Tianma radio telescope (TMRT) in order to examine the differences in physical environments, such as excitation temperature and column density of molecular clouds associated with methanol masers on the Galactic scale. Statistical results reveal that the number of 6.7 GHz methanol masers in the Perseus arm is significantly lower than that in the other three main spiral arms. In addition, the Perseus arm also has the lowest gas column density among the main spiral arms traced by the NH3 observations. Both of these findings suggest that the Perseus arm has the lowest rate of high-mass star formation compared to the other three main spiral arms. We also observed a trend in which both the luminosity of 6.7 GHz methanol masers and the ammonia gas column density decreased as the galactocentric distances. This finding indicates that the density of material in the inner Milky Way is generally higher than that in the outer Milky Way. It further suggests that high-mass stars are more easily formed at the head of spiral arms. Furthermore, we found that the column density of ammonia gas is higher in the regions on the arms than that in the inter-arm regions, supporting that the former is more likely to be the birthplace of high-mass stars. △ Less

Submitted 29 February, 2024; originally announced February 2024.

Methods for the detection of stellar rotation periods in individual TESS sectors and results from the Prime mission

Authors: Isabel L. Colman, Ruth Angus, Trevor David, Jason Curtis, Soichiro Hattori, Yuxi Lucy Lu

Abstract: For ongoing studies of the role of rotation in stellar evolution, we require large catalogs of rotation periods for testing and refining gyrochronology. While there is a wealth of data from the Kepler and K2 missions, TESS presents both an opportunity and a challenge: despite its all-sky coverage, rotation periods remain hard to detect. We analyzed individual TESS sectors to detect short-period st… ▽ More For ongoing studies of the role of rotation in stellar evolution, we require large catalogs of rotation periods for testing and refining gyrochronology. While there is a wealth of data from the Kepler and K2 missions, TESS presents both an opportunity and a challenge: despite its all-sky coverage, rotation periods remain hard to detect. We analyzed individual TESS sectors to detect short-period stellar rotation, using only parameters measured from light curves for a robust and unbiased method of evaluating detections. We used random forest classifiers for vetting, trained on a large corpus of period measurements in KELT data from the Oelkers et al. (2018) catalog and using TESS full-frame image light curves generated by eleanor (Feinstein et al. 2019). Finally, using data from the first 26 sectors of TESS, we analyzed 432,704 2-minute cadence single-sector light curves for FGKM dwarfs. We detected 16,800 periods in individual sector light curves, covering 10,909 distinct targets, and we present a catalog of the median period for each target as measured by a Lomb-Scargle periodogram. △ Less

Submitted 22 February, 2024; originally announced February 2024.

Comments: 23 pages, 15 figures, accepted for publication in The Astronomical Journal

Constraining the Yukawa Gravity with Post Newtonian Approximation using S-star Orbits around the Supermassive Black Hole in our Galactic Center

Authors: Yuan Tan, Youjun Lu

Abstract: A number of modified gravity theories (e.g., $f(R)$-gravity) lead to a Yukawa-like metric in the weak field limit which can be described by two Yukawa parameters, i.e., the strength $κかっぱ$ and the length scale $λらむだ$. The S-stars, orbiting around the supermassive black hole in the Galactic Center, are unique probes to test these gravity theories in relatively strong gravitational field. The Newtonian Yu… ▽ More A number of modified gravity theories (e.g., $f(R)$-gravity) lead to a Yukawa-like metric in the weak field limit which can be described by two Yukawa parameters, i.e., the strength $κかっぱ$ and the length scale $λらむだ$. The S-stars, orbiting around the supermassive black hole in the Galactic Center, are unique probes to test these gravity theories in relatively strong gravitational field. The Newtonian Yukawa gravity potential or a simple approximation to the Yukawa metric was usually adopted in previous studies when using the orbital motion of S-stars to constrain such modified gravity theories, which may be not sufficiently accurate considering recent and future high resolution observations. In this paper, we first derive the Post-Newtonian (PN) Yukawa motion equation at the 2PN order, and then investigate the high order effects on the orbital motions by comparison with those from the Newtonian Yukawa gravity potential. We further obtain constraints on $κかっぱ$ by using the observations on the orbital motions of several S-stars (i.e., S2, S38, and S55). Our results show that the current observations of these stars are compatible with the General Relativity and $κかっぱ$ can be constrained to $|κかっぱ|<0.01$ with $95\%$ confidence if $λらむだ\in(100,250)$\,AUえーゆー. We also estimate the possible improvements (about an order of magnitude or more) to the constraints by future higher resolution observations and the inclusion of closer S-stars, such as S4716. △ Less

Submitted 31 January, 2024; originally announced February 2024.

Comments: 15 pages, 8 figures, accepted for publication in Phys. Rev. D

Dark matter search in dwarf irregular galaxies with ten years of data from the IceCube neutrino observatory

Authors: Yi-Fei Lü, Ben-Yang Zhu, Rong-Lan Li, Xue-Kang Guo, Tian-Ci Liu, Yong-Bo Huang, Yun-Feng Liang

Abstract: Dwarf irregular galaxies (dIrrs), as rotationally-supported systems, have more reliable J-factor measurements than dwarf spheroidal galaxies and have received attention as targets for dark matter detection in recent years. In this paper, we use 10 years of IceCube muon-track data and an unbinned maximum-likelihood-ratio method to search for neutrino signals beyond the background from the direction… ▽ More Dwarf irregular galaxies (dIrrs), as rotationally-supported systems, have more reliable J-factor measurements than dwarf spheroidal galaxies and have received attention as targets for dark matter detection in recent years. In this paper, we use 10 years of IceCube muon-track data and an unbinned maximum-likelihood-ratio method to search for neutrino signals beyond the background from the directions of 7 dIrrs, aiming to detect neutrinos produced by heavy annihilation dark matter. We do not detect any significant signal. Based on such null results, we calculate the upper limits on the velocity-averaged annihilation cross section for 1 TeV - 10 PeV dark matter. Our limits, although weaker than the strictest constraints in the literature in this mass range, are also a good complement to the existing results considering the more reliable J-factor measurements of dIrrs. △ Less

Submitted 12 January, 2024; originally announced January 2024.

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

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

Authors: Siqi Zhang, Furen Deng, Youjun Lu

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

Submitted 8 January, 2024; originally announced January 2024.

Comments: 12 pages, 7 figures, MNRAS accepted

Detectability of Strongly Gravitationally Lensed Tidal Disruption Events

Authors: Zhiwei Chen, Youjun Lu, Yunfeng Chen

Abstract: More than 100 tidal disruption events (TDEs) have been detected at multi-bands, which can be viewed as extreme laboratories to investigate the accretion physics and gravity in the immediate vicinity of massive black holes. Future transient surveys are expected to detect several tens of thousands of TDEs, among which a small fraction may be strongly gravitationally lensed by intervening galaxies. I… ▽ More More than 100 tidal disruption events (TDEs) have been detected at multi-bands, which can be viewed as extreme laboratories to investigate the accretion physics and gravity in the immediate vicinity of massive black holes. Future transient surveys are expected to detect several tens of thousands of TDEs, among which a small fraction may be strongly gravitationally lensed by intervening galaxies. In this paper, we statistically etsimate the detection rate of lensed TDEs, with dependence on the limiting magnitude of the transient all-sky surveys searching for them. We find that the requisite limiting magnitude for an all-sky transient survey to observe at least $1$ yr$^{-1}$ is $\gtrsim 21.3$, $21.2$, and $21.5$ mag in the $u$, $g$, and $z$ bands, respectively. If the limiting magnitude of the all-sky survey can reach $\sim 25-26$ mag in the $u$, $g$, and $z$ bands, the detection rate can be upto about several tens to hundreds per year. The discovery and identification of the first image of the lensed TDE can be taken as an early warning of the second and other subsequent images, which may enable detailed monitoring of the pre-peak photometry and spectroscopy evolution of the TDE. The additional early-stage information may help to constrain the dynamical and radiation processes involving in the TDEs. △ Less

Submitted 25 January, 2024; v1 submitted 1 January, 2024; originally announced January 2024.

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

NICER views moderate, strong, and extreme photospheric expansion bursts from the ultracompact X-ray binary 4U 1820$-$30

Authors: Wenhui Yu, Zhaosheng Li, Yongqi Lu, Yuanyue Pan, Xuejuan Yang, Yupeng Chen, Shu Zhang, Maurizio Falanga

Abstract: Type I X-ray bursts in the ultracompact X-ray binary 4U 1820$-$30 are powered by the unstable thermonuclear burning of hydrogen-deficient material. We report the detection of 15 type I X-ray bursts from 4U 1820$-$30 observed by NICER in between 2017 and 2023. All these bursts occurred in the low state for the persistent flux in the range of $2.5-8\times10^{-9}~{\rm erg~s^{-1}~cm^{-2}}$ in 0.1$-$25… ▽ More Type I X-ray bursts in the ultracompact X-ray binary 4U 1820$-$30 are powered by the unstable thermonuclear burning of hydrogen-deficient material. We report the detection of 15 type I X-ray bursts from 4U 1820$-$30 observed by NICER in between 2017 and 2023. All these bursts occurred in the low state for the persistent flux in the range of $2.5-8\times10^{-9}~{\rm erg~s^{-1}~cm^{-2}}$ in 0.1$-$250 keV. The burst spectra during the tail can be well explained by blackbody model. However, for the first $\sim$5 s after the burst onset, the time-resolved spectra showed strong deviations from the blackbody model. The significant improvement of the fit can be obtained by taking into account of the enhanced persistent emission due to the Poynting-Robterson drag, the extra emission modelled by another blackbody component or by the reflection from the surrounding accretion disk. The reflection model provides a self-consistent and physically motivated explanation. We find that the accretion disk density changed with 0.5 s delay as response to the burst radiation, which indicates the distortion of the accretion disk during X-ray bursts. From the time-resolved spectroscopy, all bursts showed the characteristic of photospheric radius expansion (PRE). We find one superexpansion burst with the extreme photospheric radius $r_{\rm ph}>10^3$ km and blackbody temperature of $\sim 0.2$ keV, thirteen strong PRE bursts for $r_{\rm ph}>10^2$ km, and one moderate PRE burst for $r_{\rm ph}\sim55$ km. △ Less

Submitted 27 December, 2023; originally announced December 2023.

Comments: 16 pages, 13 figures, 4 tables. Accepted for publication in Astronomy & Astrophysics

Feeding plankton to whales: high-redshift supermassive black holes from tiny black hole explosions

Authors: Yifan Lu, Zachary S. C. Picker, Alexander Kusenko

Abstract: Recent observations of the high-redshift universe have uncovered a significant number of active galactic nuclei, implying that supermassive black holes (SMBHs) would have to have been formed at much earlier times than expected. Direct collapse of metal-free gas clouds to SMBHs after recombination could help explain the early formation of SMBHs, but this scenario is stymied by the fragmentation of… ▽ More Recent observations of the high-redshift universe have uncovered a significant number of active galactic nuclei, implying that supermassive black holes (SMBHs) would have to have been formed at much earlier times than expected. Direct collapse of metal-free gas clouds to SMBHs after recombination could help explain the early formation of SMBHs, but this scenario is stymied by the fragmentation of the clouds due to efficient molecular hydrogen cooling. We show that a subdominant population of tiny, evaporating primordial black holes, with significant clustering in some gas clouds, can heat the gas sufficiently so that molecular hydrogen is not formed, and direct collapse to to black holes is possible even at high redshifts. △ Less

Submitted 11 June, 2024; v1 submitted 22 December, 2023; originally announced December 2023.

Comments: 15 pages, 4 figures

Report number: IPMU23-0052

Journal ref: Phys. Rev. D 109, 123016 (2024)

In this Day and Age: An Empirical Gyrochronology Relation for Partially and Fully Convective Single Field Stars

Authors: Yuxi Lu, Ruth Angus, Daniel Foreman-Mackey, Soichiro Hattori

Abstract: Gyrochronology, the field of age-dating stars using mainly their rotation periods and masses, is ideal for inferring the ages of individual main-sequence stars. However, due to the lack of physical understanding of the complex magnetic fields in stars, gyrochronology relies heavily on empirical calibrations that require consistent and reliable stellar age measurements across a wide range of period… ▽ More Gyrochronology, the field of age-dating stars using mainly their rotation periods and masses, is ideal for inferring the ages of individual main-sequence stars. However, due to the lack of physical understanding of the complex magnetic fields in stars, gyrochronology relies heavily on empirical calibrations that require consistent and reliable stellar age measurements across a wide range of periods and masses. In this paper, we obtain a sample of consistent ages using the gyro-kinematic age-dating method, a technique to calculate the kinematics ages of stars. Using a Gaussian Process model conditioned on ages from this sample (~ 1 - 14 Gyr) and known clusters (0.67 - 3.8 Gyr), we calibrate the first empirical gyrochronology relation that is capable of inferring ages for single, main-sequence stars between 0.67 Gyr to 14 Gyr. Cross-validating and testing results suggest our model can infer cluster and asteroseismic ages with an average uncertainty of just over 1 Gyr. With this model, we obtain gyrochronology ages for ~ 100,000 stars within 1.5 kpc of the Sun with period measurements from Kepler and ZTF, and 384 unique planet host stars. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: Submitted to AJ. Missing citations welcome

Sphaleron freeze-in baryogenesis with gravitational waves from the QCD transition

Authors: Fei Gao, Julia Harz, Chandan Hati, Yi Lu, Isabel M. Oldengott, Graham White

Abstract: A large primordial lepton asymmetry is capable of explaining the baryon asymmetry of the Universe (BAU) through suppression of the electroweak sphaleron rates (``sphaleron freeze-in") which can lead to a first-order cosmic QCD transition with an observable gravitational wave (GW) signal. With next-to-leading order dimensional reduction and the exact 1-loop fluctuation determinant, we accurately co… ▽ More A large primordial lepton asymmetry is capable of explaining the baryon asymmetry of the Universe (BAU) through suppression of the electroweak sphaleron rates (``sphaleron freeze-in") which can lead to a first-order cosmic QCD transition with an observable gravitational wave (GW) signal. With next-to-leading order dimensional reduction and the exact 1-loop fluctuation determinant, we accurately compute the lepton asymmetry needed to realize this paradigm, finding it to be an order of magnitude smaller than previous estimates. Further, we apply an improved QCD equation of state capable of describing the phase transition line together with the critical endpoint leading to better agreement with lattice and functional QCD results. Based on this, we identify the range of lepton flavor asymmetries inducing a first-order cosmic QCD transition. We then extract the parameters relevant to the prediction of GW signal from a first-order cosmic QCD transition. Our result showcases the possibility of probing the sphaleron freeze-in paradigm as an explanation of BAU by future gravitational wave experiments like $μみゅー$Ares. △ Less

Submitted 1 September, 2023; originally announced September 2023.

Comments: 7 pages, 3 figures

Report number: MITP-23-049, ULB-TH/23-08, IFIC/23-41

Structure Formation after Reheating: Supermassive Primordial Black Holes and Fermi Ball Dark Matter

Authors: Marcos M. Flores, Yifan Lu, Alexander Kusenko

Abstract: In the presence of (relatively) long-range forces, structures can form even during the radiation dominated era, leading to compact objects, such as Fermi balls or primordial black holes (PBHs), which can account for all or part of dark matter. We present a detailed analysis of a model in which fermions are produced from the inflaton decay developing some particle-antiparticle asymmetry. These ferm… ▽ More In the presence of (relatively) long-range forces, structures can form even during the radiation dominated era, leading to compact objects, such as Fermi balls or primordial black holes (PBHs), which can account for all or part of dark matter. We present a detailed analysis of a model in which fermions are produced from the inflaton decay developing some particle-antiparticle asymmetry. These fermions undergo clustering and structure formation driven by a Yukawa interaction. The same interaction provides a cooling channel for the dark halos via scalar radiation, leading to rapid collapse and the formation of a compact object. We discuss the criteria for the formation of either PBHs and Fermi balls. In the PBH formation regime, supermassive PBHs can seed the active galactic nuclei or quasars found at high redshift. Alternatively, Fermi balls can account for all of the cold dark matter, while evading microlensing constraints. △ Less

Submitted 7 December, 2023; v1 submitted 17 August, 2023; originally announced August 2023.

Comments: 9 pages, 6 figures

Report number: IPMU23-0027

Journal ref: Phys. Rev. D 108, 123511 (2023)

Detecting active latitudes of Sun-like stars using asteroseismic a-coefficients

Authors: Othman Benomar, Masao Takata, Michael Bazot, Takashi Sekii, Laurent Gizon, Yuting Lu

Abstract: We introduce a framework to measure the asphericity of Sun-like stars using $a_1$, $a_2$ and $a_4$ coefficients, and constrain their latitudes of magnetic activity. Systematic errors on the inferred coefficients are evaluated in function of key physical and seismic parameters (inclination of rotation axis, average rotation, height-to-noise ratio of peaks in power spectrum). The measured a-coeffici… ▽ More We introduce a framework to measure the asphericity of Sun-like stars using $a_1$, $a_2$ and $a_4$ coefficients, and constrain their latitudes of magnetic activity. Systematic errors on the inferred coefficients are evaluated in function of key physical and seismic parameters (inclination of rotation axis, average rotation, height-to-noise ratio of peaks in power spectrum). The measured a-coefficients account for rotational oblateness and the effect of surface magnetic activity. We use a simple model that assumes a single latitudinal band of activity. Using solar SOHO/VIRGO/SPM data, we demonstrate the capability of the method to detect the mean active latitude and its intensity changes between 1999-2002 (maximum of activity) and 2006-2009 (minimum of activity). We further apply the method to study the solar-analogue stars 16 Cyg A and B using Kepler observations. An equatorial band of activity, exhibiting intensity that could be comparable to that of the Sun, is detected in 16 Cyg A. However, 16 Cyg B exhibits a bi-modality in $a_4$ that is challenging to explain. We suggest that this could be a manifestation of the transition between a quiet and an active phase of activity. Validating or invalidating this hypothesis may require new observations. △ Less

Submitted 17 August, 2023; originally announced August 2023.

Comments: Accepted to A&A

Journal ref: A&A 680, A27 (2023)

Effect of Neutrinos on Angular Momentum of Dark Matter Halo

Authors: Yu Chen, Chang-Zhi Lu, Yu Lu, Tingting Zhang, Tong-Jie Zhang

Abstract: Massive neutrinos are expected to affect the large-scale structure formation, including the major component of solid substances, dark matter halos. How halos are influenced by neutrinos is vital and interesting, and angular momentum (AM) as a significant feature provides a statistical perspective for this issue. Exploring halos from TianNu N-body cosmological simulation with the co-evolving neutri… ▽ More Massive neutrinos are expected to affect the large-scale structure formation, including the major component of solid substances, dark matter halos. How halos are influenced by neutrinos is vital and interesting, and angular momentum (AM) as a significant feature provides a statistical perspective for this issue. Exploring halos from TianNu N-body cosmological simulation with the co-evolving neutrino particles, we obtain some concrete conclusions. First, by comparing the same halos with and without neutrinos, in contrast to the neutrino-free case, over 89.71\% of halos have smaller halo moduli, over 71.06\% have smaller particle-mass-reduced (PMR) AM moduli, and over 95.44\% change their orientations of less than $0.65^\circ$. Moreover, the relative variation of PMR modulus is more visible for low-mass halos. Second, to explore the PMR moduli of halos in dense or sparse areas, we divide the whole box into big cubes, and search for halos within a small spherical cell in a single cube. From the two-level divisions, we discover that in denser cubes, the variation of PMR moduli with massive neutrinos decreases more significantly. This distinction suggests that neutrinos exert heavier influence on halos' moduli in compact regions. With massive neutrinos, most halos (86.60\%) have lower masses than without neutrinos. △ Less

Submitted 18 July, 2023; v1 submitted 18 July, 2023; originally announced July 2023.

Comments: 12 pages(without ref.), 12 figures. The statistical research of the angular momentum moduli of dark halos in TianNu N-body simulation. Second update only corrects the author name on the arXiv page

Journal ref: Res. Astron. Astrophys. 23:085025 (2023)

Chemical Cartography of the Sagittarius Stream with Gaia

Authors: Emily C. Cunningham, Jason A. S. Hunt, Adrian M. Price-Whelan, Kathryn V. Johnston, Melissa K. Ness, Yuxi Lu, Ivanna Escala, Ioana A. Stelea

Abstract: The stellar stream connected to the Sagittarius (Sgr) dwarf galaxy is the most massive tidal stream that has been mapped in the Galaxy, and is the dominant contributor to the outer stellar halo of the Milky Way. We present metallicity maps of the Sgr stream, using 34,240 red giant branch stars with inferred metallicities from Gaia BP/RP spectra. This sample is larger than previous samples of Sgr s… ▽ More The stellar stream connected to the Sagittarius (Sgr) dwarf galaxy is the most massive tidal stream that has been mapped in the Galaxy, and is the dominant contributor to the outer stellar halo of the Milky Way. We present metallicity maps of the Sgr stream, using 34,240 red giant branch stars with inferred metallicities from Gaia BP/RP spectra. This sample is larger than previous samples of Sgr stream members with chemical abundances by an order of magnitude. We measure metallicity gradients with respect to Sgr stream coordinates $(Λらむだ, B)$, and highlight the gradient in metallicity with respect to stream latitude coordinate $B$, which has not been observed before. We find $\nabla \mathrm{[M/H]} = -2.48 \pm 0.08 \times 10^{-2}$ dex/deg above the stream track ($B>B_0$ where $B_0=1.5$ deg is the latitude of the Sgr remnant) and $\nabla \mathrm{[M/H]} =- 2.02 \pm 0.08 \times 10^{-2}$ dex/deg below the stream track ($B<B_0$). By painting metallicity gradients onto a tailored N-body simulation of the Sgr stream, we find that the observed metallicities in the stream are consistent with an initial radial metallicity gradient in the Sgr dwarf galaxy of $\sim -0.1$ to $-0.2$ dex/kpc, well within the range of observed metallicity gradients in Local Group dwarf galaxies. Our results provide novel observational constraints for the internal structure of the dwarf galaxy progenitor of the Sgr stream. Leveraging new large datasets in conjunction with tailored simulations, we can connect the present day properties of disrupted dwarfs in the Milky Way to their initial conditions. △ Less

Submitted 17 July, 2023; originally announced July 2023.

Comments: 20 pages, 12 figures. Submitted to ApJ; comments welcome!

On the detection of the electromagnetic counterparts from lensed gravitational wave events by binary neutron star mergers

Authors: Hao Ma, Youjun Lu, Xiao Guo, Siqi Zhang, Qingbo Chu

Abstract: Future ground-based gravitational wave (GW) detectors, i.e., Einstein telescope (ET) and Cosmic Explorer (CE), are expected to detect a significant number of lensed binary neutron star (BNS) mergers, which may provide a unique tool to probe cosmology. In this paper, we investigate the detectability of the optical/infrared electromagnetic (EM) counterparts (kilonovae/afterglows) from these lensed B… ▽ More Future ground-based gravitational wave (GW) detectors, i.e., Einstein telescope (ET) and Cosmic Explorer (CE), are expected to detect a significant number of lensed binary neutron star (BNS) mergers, which may provide a unique tool to probe cosmology. In this paper, we investigate the detectability of the optical/infrared electromagnetic (EM) counterparts (kilonovae/afterglows) from these lensed BNS mergers by future GW detectors and EM telescopes using simple kilonova, afterglow, and lens models. ET and CE are expected to detect $\sim5.32^{+26.1}_{-5.10}$ and $67.3^{+332}_{-64.7}$ lensed BNS mergers per year. We find that the EM counterparts associated with all these mergers will be detectable by an all sky-survey in the H-band with the limiting magnitude $m_{\textrm{lim}}\gtrsim27$, while the detectable fraction is $\lesssim0.4\%$ in the g-/z-band if with $m_{\textrm{lim}}\lesssim24$. Generally it is more efficient to search the lensed EM counterparts by adopting the infrared bands than the optical/UV bands with the same $m_{\textrm{lim}}$. Future telescopes like Vera C. Rubin Observatory, China Space Station Telescope, and Euclid can hardly detect the EM counterparts of even one lensed BNS merger. Roman Space Telescope (RST) and James Webb Space Telescope (JWST) have the capability to detect about a few or more such events per year. Moreover, the time delays and separations between the lensed image pairs are typically in the ranges from minutes to months and from $0.1$ to $1$\,arcsec, suggesting that both the GW and EM images of most lensed BNS mergers can be well resolved by not only CE/ET in the time domain but also RST/JWST spatially. △ Less

Submitted 5 July, 2023; originally announced July 2023.

Comments: 16 pages, 10 figures, MNRAS in press

Diffractive lensing of nano-Hertz gravitational waves emitted from supermassive binary black holes by intervening galaxies

Authors: Hao Ma, Youjun Lu, Zhiwei Chen, Yunfeng Chen

Abstract: Pulsar timing array (PTA) experiments are expected to detect nano-Hertz gravitational waves (GWs) emitted from individual inspiralling supermassive binary black holes (SMBBHs). The GW signals from a small fraction of these SMBBHs may be diffractively lensed by intervening galaxies. In this paper, we investigate the diffractive lensing effects on the continuous GW signals from the lensed SMBBHs and… ▽ More Pulsar timing array (PTA) experiments are expected to detect nano-Hertz gravitational waves (GWs) emitted from individual inspiralling supermassive binary black holes (SMBBHs). The GW signals from a small fraction of these SMBBHs may be diffractively lensed by intervening galaxies. In this paper, we investigate the diffractive lensing effects on the continuous GW signals from the lensed SMBBHs and estimate the detectable number of such signals by PTAs, such as the Chinese PTA (CPTA) and the Square Kilometer Array (SKA) PTA. We find that the amplitude of the lensed GW signals may be only amplified by a factor of $\sim 1.01-1.14$ ($16\%-84\%$ range) and the phase of the signals may shift somewhat due to the lensing, significantly different from those strongly lensed high frequency GW signals from compact binary mergers in the geometric optics. We estimate that $\sim 0.01\%$ of all detected nano-Hertz GW signals from individual SMBBHs by future PTA experiments are lensed by foreground galaxies (i.e., up to $\sim 106$ for CPTA and up to $\sim 289$ for SKA-PTA). However, the lensed nano-Hertz GW signals are difficult to be distinguished from those without lensing by the PTA observations only. We further discuss the possibility about the identification of the lensed nano-Hertz GW signals from SMBBHs via the electromagnetic detection of their host galaxies or active galactic nuclei. △ Less

Submitted 5 July, 2023; originally announced July 2023.

Comments: 12 pages, 8 figures, Accepted for publication in MNRAS

Pulsar timing array detections of supermassive binary black holes: implications from the detected common process signal and beyond

Authors: Yunfeng Chen, Qingjuan Yu, Youjun Lu

Abstract: Pulsar timing arrays (PTAs) are anticipated to detect the stochastic gravitational wave background (GWB) from supermassive binary black holes (BBHs) as well as the gravitational waves from individual BBHs. Recently, a common process signal was reported by several PTAs. In this paper, we investigate the constraints on the BBH population model(s) by current PTA observations and further study the det… ▽ More Pulsar timing arrays (PTAs) are anticipated to detect the stochastic gravitational wave background (GWB) from supermassive binary black holes (BBHs) as well as the gravitational waves from individual BBHs. Recently, a common process signal was reported by several PTAs. In this paper, we investigate the constraints on the BBH population model(s) by current PTA observations and further study the detections of both the GWB and individual BBHs by current/future PTAs. We find that the MBH--host galaxy scaling relation, an important ingredient of the BBH population model, is required to either evolve significantly with redshift or have a normalization $\sim0.86-1.1$ dex higher than the empirical ones, if the GWB is the same as the common process signal. For both cases, the estimated detection probability for individual BBHs is too small for a positive detection by current PTAs. By involving either the constrained scaling relations or those empirical ones into the BBH population models, we estimate that the GWB may be detected with a signal-to-noise ratio $\gtrsim3$ by the PTAs based on the Five hundred meter Aperture Spherical radio Telescope (CPTA) and the Square Kilometer Array (SKAPTA) after $\sim2-3$ (or $\sim6-11$) years' observation, if it is the same as (or an order of magnitude lower than) the common process signal. The detection time of individual BBHs by CPTA and SKAPTA is close to that of the GWB detection. We show that the BBH population model can be strongly constrained by the number and property distributions of BBHs to be detected by future PTAs. △ Less

Submitted 19 June, 2023; originally announced June 2023.

Comments: 25 pages, 9 figures, submitted to ApJ on April 10, 2023, and revised after reviewer's report

The correlations of stellar tidal disruption rates with properties of massive black holes and their host galaxies

Authors: Yunfeng Chen, Qingjuan Yu, Youjun Lu

Abstract: Stars can be either disrupted as tidal disruption events (TDEs) or swallowed as a whole by massive black holes (MBHs) at galactic centers when they approach sufficiently close to these MBHs. In this work, we investigate the correlations of such stellar consumption rates with both the MBH mass $M_{\rm bh}$ and the inner slope of the host galaxy mass density distribution $αあるふぁ$. We introduce a simplifi… ▽ More Stars can be either disrupted as tidal disruption events (TDEs) or swallowed as a whole by massive black holes (MBHs) at galactic centers when they approach sufficiently close to these MBHs. In this work, we investigate the correlations of such stellar consumption rates with both the MBH mass $M_{\rm bh}$ and the inner slope of the host galaxy mass density distribution $αあるふぁ$. We introduce a simplified analytical power-law model with a power-law stellar mass density distribution surrounding MBHs and separate the contributions of two-body relaxation and stellar orbital precession for the stellar orbital angular momentum evolution in nonspherical galaxy potentials. The stellar consumption rates derived from this simplified model can be well consistent with the numerical results obtained with a more realistic treatment of stellar distributions and dynamics around MBHs, providing an efficient way to estimate TDE rates. The origin of the correlations of stellar consumption rates with $M_{\rm bh}$ and $αあるふぁ$ are explained by the dependence of this analytical model on those MBH/host galaxy properties and by the separation of the stellar angular momentum evolution mechanisms. We propose that the strong positive correlation between the rates of stellar consumption due to two-body relaxation and $αあるふぁ$ provides one interpretation for the overrepresentation of TDEs found in some rare E+A/poststarburst galaxies. We find high TDE rates for giant stars, up to those for solar-type stars. The understanding of the origin of the correlations of the stellar consumption rates will be necessary for obtaining the demographics of MBHs and their host galaxies via TDEs. △ Less

Submitted 19 June, 2023; originally announced June 2023.

Comments: 17 pages, 6 figures, ApJ in press

An abrupt change in the stellar spin-down law at the fully convective boundary

Authors: Yuxi Lu, Victor See, Louis Amard, Ruth Angus, Sean P. Matt

Abstract: The importance of the existence of a radiative core in generating a solar-like magnetic dynamo is still unclear. Analytic models and magnetohydrodynamic simulations of stars suggest the thin layer between a star's radiative core and its convective zone can produce shearing that reproduces key characteristics of a solar-like dynamo. However, recent studies suggest fully and partially convective sta… ▽ More The importance of the existence of a radiative core in generating a solar-like magnetic dynamo is still unclear. Analytic models and magnetohydrodynamic simulations of stars suggest the thin layer between a star's radiative core and its convective zone can produce shearing that reproduces key characteristics of a solar-like dynamo. However, recent studies suggest fully and partially convective stars exhibit very similar period-activity relations, hinting that dynamos generated by stars with and without radiative cores hold similar properties. Here, using kinematic ages, we discover an abrupt change in the stellar spin-down law across the fully convective boundary. We found that fully convective stars exhibit a higher angular momentum loss rate, corresponding to a torque that is $\sim$ 2.25 times higher for a given angular velocity than partially convective stars around the fully convective boundary. This requires a dipole field strength that is larger by a factor of $\sim$2.5, a mass loss rate that is $\sim$4.2 times larger, or some combination of both of those factors. Since stellar-wind torques depend primarily on large-scale magnetic fields and mass loss rates, both of which derive from magnetic activity, the observed abrupt change in spin-down law suggests that the dynamos of partially and fully convective stars may be fundamentally different △ Less

Submitted 15 June, 2023; originally announced June 2023.

Comments: Under review at Nature Astronomy

Scalar-Induced Gravitational Waves from Ghost Inflation and Parity Violation

Authors: Sebastian Garcia-Saenz, Yizhou Lu, Zhiming Shuai

Abstract: We calculate the scalar-induced gravitational wave energy density in the theory of Ghost Inflation, assuming scale invariance and taking into account both the power spectrum- and trispectrum-induced contributions. For the latter we consider the leading cubic and quartic couplings of the comoving curvature perturbation in addition to two parity-violating quartic operators. In the parity-even case,… ▽ More We calculate the scalar-induced gravitational wave energy density in the theory of Ghost Inflation, assuming scale invariance and taking into account both the power spectrum- and trispectrum-induced contributions. For the latter we consider the leading cubic and quartic couplings of the comoving curvature perturbation in addition to two parity-violating quartic operators. In the parity-even case, we find the relative importance of the trispectrum-induced signal to be suppressed by the requirement of perturbativity, strengthening a no-go theorem recently put forth. The parity-odd signal, even though also bound to be small, is non-degenerate with the Gaussian contribution and may in principle be comparable to the parity-even non-Gaussian part, thus potentially serving as a probe of the Ghost Inflation scenario and of parity violating physics during inflation. △ Less

Submitted 16 November, 2023; v1 submitted 15 June, 2023; originally announced June 2023.

Comments: 36pages

Journal ref: Phys.Rev. D108 (2023) no.12, 123507

The Structure and Dynamics of Massive High-$z$ Cosmic-Web Filaments: Three Radial Zones in Filament Cross-Sections

Authors: Yue Samuel Lu, Nir Mandelker, S. Peng Oh, Avishai Dekel, Frank C. van den Bosch, Volker Springel, Daisuke Nagai, Freeke van de Voort

Abstract: We analyse the internal structure and dynamics of cosmic-web filaments that connect massive high-$z$ haloes. Our analysis is based on a high-resolution AREPO cosmological simulation zooming-in on a volume encompassing three ${\rm Mpc}$-scale filaments feeding three massive haloes of $\sim 10^{12}\,\text{M}_\odot$ at $z \sim 4$, embedded in a large-scale sheet. Each filament is surrounded by a cyli… ▽ More We analyse the internal structure and dynamics of cosmic-web filaments that connect massive high-$z$ haloes. Our analysis is based on a high-resolution AREPO cosmological simulation zooming-in on a volume encompassing three ${\rm Mpc}$-scale filaments feeding three massive haloes of $\sim 10^{12}\,\text{M}_\odot$ at $z \sim 4$, embedded in a large-scale sheet. Each filament is surrounded by a cylindrical accretion shock of radius $r_{\rm shock} \sim 50 \,{\rm kpc}$. The post-shock gas is in virial equilibrium with the potential well set by an isothermal dark-matter filament. The filament line-mass is $\sim 9\times 10^8\,\text{M}_\odot\,{\rm kpc}^{-1}$, the gas fraction within $r_{\rm shock}$ is the universal baryon fraction, and the virial temperature is $\sim 7\times 10^5 {\rm K}$. In the outer ''thermal'' (T) zone, $r \geq 0.65 \, r_{\rm shock}$, inward gravity and ram-pressure forces are over-balanced by outwards thermal pressure forces, decelerating the inflowing gas expanding the shock outward. In the intermediate ''vortex'' (V) zone, $0.25 \leq r/ r_{\rm shock} \leq 0.65$, the velocity field is dominated by a quadrupolar vortex structure due to offset inflow along the sheet through the post-shock gas. The outwards force is dominated by centrifugal forces associated with these vortices, with additional contributions from global rotation and thermal pressure. The shear and turbulent forces associated with the vortices act inward. The inner ''stream'' (S) zone, $r < 0.25 \, r_{\rm shock}$, is a dense isothermal core, $T\sim 3 \times 10^4 \, {\rm K}$ and $n_{\rm H}\sim 0.01 \,{\rm cm^{-3}}$, defining the cold streams that feed galaxies. The core is formed by an isobaric cooling flow and is associated with a decrease in outwards forces, though it exhibits both inflows and outflows. [abridged] △ Less

Submitted 20 November, 2023; v1 submitted 6 June, 2023; originally announced June 2023.

Comments: 32 pages, 22 figures, submitted to MNRAS

Searching for dark-matter induced neutrino signals in dwarf spheroidal galaxies using 10 years of IceCube public data

Authors: Xue-Kang Guo, Yi-Fei Lü, Yong-Bo Huang, Rong-Lan Li, Ben-Yang Zhu, Yun-Feng Liang

Abstract: This study searches for neutrino signals from 18 dwarf spheroidal galaxies (dSphs) using 10 years of publicly available muon-track data of the IceCube neutrino observatory. We apply an unbinned likelihood analysis on each of these dSphs to derive the significance the putative neutrino emission. To further enhance our sensitivity, we also stack all dSphs together to perform a joint analysis. Howeve… ▽ More This study searches for neutrino signals from 18 dwarf spheroidal galaxies (dSphs) using 10 years of publicly available muon-track data of the IceCube neutrino observatory. We apply an unbinned likelihood analysis on each of these dSphs to derive the significance the putative neutrino emission. To further enhance our sensitivity, we also stack all dSphs together to perform a joint analysis. However, no significant neutrino emission signal was detected in either the single-source or stacking analysis. Based on these null results, we derive constraints on the annihilation cross section of dark matter particles. Compared to the existing literature, our constraints via the channel $χかいχかい\rightarrowμみゅー^+μみゅー^-$ are comparable to the ones from the VERITAS observations of dSphs. △ Less

Submitted 14 July, 2023; v1 submitted 5 June, 2023; originally announced June 2023.

Multiband gravitational wave observations of stellar binary black holes at the low to middle and high frequencies

Authors: Yuetong Zhao, Youjun Lu, Changshuo Yan, Zhiwei Chen, Wei-Tou Ni

Abstract: The ground-based gravitational wave (GW) observatories discover a population of merging stellar binary black holes (BBHs), which are promising targets for multiband observations by the low-, middle-, and high-frequency GW detectors. In this paper, we investigate the multiband GW detections of BBHs and demonstrate the advantages of such observations in improving the localization and parameter estim… ▽ More The ground-based gravitational wave (GW) observatories discover a population of merging stellar binary black holes (BBHs), which are promising targets for multiband observations by the low-, middle-, and high-frequency GW detectors. In this paper, we investigate the multiband GW detections of BBHs and demonstrate the advantages of such observations in improving the localization and parameter estimates of the sources. We generate mock samples of BBHs by considering different formation models as well as the merger rate density constrained by the current observations (GWTC-3). We specifically consider the astrodynamical middle-frequency interferometer GW observatory (AMIGO) in the middle-frequency band and estimate that it may detect $21$-$91$ BBHs with signal-to-noise ratio $\varrho\geq8$ in a $4$-yr observation period. The multiband observations by the low-frequency detectors [Laser Interferometer Space Antenna (LISA) and Taiji] and AMIGO may detect $5$-$33$ BBHs with $\varrho_{\rm LT}\geq5$ and $\varrho_{\rm AMI}\geq5$, which can evolve to the high-frequency band within $4$ yr and can be detected by the Cosmic Explorer (CE) and Einstein Telescope (ET). The joint observations of LISA-Taiji-AMIGO-ET-CE may localize the majority of the detectable BBHs in sky areas of $7\times10^{-7}$ to $2\times10^{-3}$ deg$^2$, which is improved by a factor of $\sim120$, $\sim2.4\times10^{5}$, $\sim1.8\times10^{4}$, or $\sim1.2\times10^{4}$, comparing with those by only adopting CE-ET, AMIGO, LISA-Taiji, or LISA-Taiji-AMIGO. These joint observations can also lead to an improvement of the measurement precision of the chirp mass (symmetric mass ratio) by a factor of $\sim5.5\times10^{4}$ ($33$), $\sim16$ ($8$), $\sim120$ ($90$), or $\sim5$ ($5$), comparing with those by CE-ET, AMIGO, LISA-Taiji, or LISA-Taiji-AMIGO. △ Less

Submitted 5 June, 2023; originally announced June 2023.

Comments: 16 pages, 12 figures, MNRAS published

Journal ref: 2023, MNRAS, 522, 2951

On networks of space-based gravitational-wave detectors

Authors: Rong-Gen Cai, Zong-Kuan Guo, Bin Hu, Chang Liu, Youjun Lu, Wei-Tou Ni, Wen-Hong Ruan, Naoki Seto, Gang Wang, Yue-Liang Wu

Abstract: The space-based laser interferometers, LISA, Taiji and TianQin, are targeting to observe milliHz gravitational waves (GWs) in the 2030s. The joint observations from multiple space-based detectors yield significant advantages. In this work, we recap the studies and investigations for the joint space-based GW detector networks to highlight: 1) the high precision of sky localization for the massive b… ▽ More The space-based laser interferometers, LISA, Taiji and TianQin, are targeting to observe milliHz gravitational waves (GWs) in the 2030s. The joint observations from multiple space-based detectors yield significant advantages. In this work, we recap the studies and investigations for the joint space-based GW detector networks to highlight: 1) the high precision of sky localization for the massive binary black hole (BBH) coalescences and the GW sirens in the cosmological implication, 2) the effectiveness to test the parity violation in the stochastic GW background observations, 3) the efficiency of subtracting galactic foreground, 4) the improvement in stellar-mass BBH observations. We inspect alternative networks by trading off massive BBH observations and stochastic GW background observation. △ Less

Submitted 22 May, 2023; v1 submitted 8 May, 2023; originally announced May 2023.

Comments: 32 pages, 14 figures, reviews on recent studies about space-borne GW networks, comments and feedbacks are welcome

Inspiraling streams of enriched gas observed around a massive galaxy 11 billion years ago

Authors: Shiwu Zhang, Zheng Cai, Dandan Xu, Rhythm Shimakawa, Fabrizio Arrigoni Battaia, Jason Xavier Prochaska, Renyue Cen, Zheng Zheng, Yunjing Wu, Qiong Li, Liming Dou, Jianfeng Wu, Ann Zabludoff, Xiaohui Fan, Yanli Ai, Emmet Gabriel Golden-Marx, Miao Li, Youjun Lu, Xiangcheng Ma, Sen Wang, Ran Wang, Feng Yuan

Abstract: Stars form in galaxies, from gas that has been accreted from the intergalactic medium. Simulations have shown that recycling of gas-the reaccretion of gas that was previously ejected from a galaxy-could sustain star formation in the early Universe. We observe the gas surrounding a massive galaxy at redshift 2.3 and detect emission lines from neutral hydrogen, helium, and ionized carbon that extend… ▽ More Stars form in galaxies, from gas that has been accreted from the intergalactic medium. Simulations have shown that recycling of gas-the reaccretion of gas that was previously ejected from a galaxy-could sustain star formation in the early Universe. We observe the gas surrounding a massive galaxy at redshift 2.3 and detect emission lines from neutral hydrogen, helium, and ionized carbon that extend 100 kiloparsecs from the galaxy. The kinematics of this circumgalactic gas is consistent with an inspiraling stream. The carbon abundance indicates that the gas had already been enriched with elements heavier than helium, previously ejected from a galaxy. We interpret the results as evidence of gas recycling during high-redshift galaxy assembly. △ Less

Submitted 3 May, 2023; originally announced May 2023.

Comments: Published in Science, 5 May 2023 (accepted version), Main text 20 pages, four figures in the main text, and 13 figures and 4 tables in the supplementary materials;

The luminosity functions of kilonovae from binary neutron star mergers under different equation of states

Authors: Chunyang Zhao, Youjun Lu, Qingbo Chu, Wen Zhao

Abstract: Kilonovae produced by mergers of binary neutron stars (BNSs) are important transient events to be detected by time domain surveys with the alerts from the ground-based gravitational wave detectors. The observational properties of these kilonovae depend on the physical processes involved in the merging processes and the equation of state (EOS) of neutron stars (NSs). In this paper, we investigate t… ▽ More Kilonovae produced by mergers of binary neutron stars (BNSs) are important transient events to be detected by time domain surveys with the alerts from the ground-based gravitational wave detectors. The observational properties of these kilonovae depend on the physical processes involved in the merging processes and the equation of state (EOS) of neutron stars (NSs). In this paper, we investigate the dependence of kilonova luminosities on the parameters of BNS mergers, and estimate the distribution functions of kilonova peak luminosities (KLFs) at the u, g, r, i, y, and z bands as well as its dependence on the NS EOS, by adopting a comprehensive semi-analytical model for kilonovae (calibrated by the observations of GW170817), a population synthesis model for the cosmic BNSs, and the ejecta properties of BNS mergers predicted by numerical simulations. We find that the kilonova light curves depend on both the BNS properties and the NS EOS, and the KLFs at the considered bands are bimodal with the bright components mostly contributed by BNS mergers with total mass $\lesssim 3.2M_\odot$/$2.8M_\odot$ and fainter components mostly contributed by BNS mergers with total mass $\gtrsim 3.2M_\odot$/$2.8M_\odot$ by assuming a stiff/soft (DD2/SLy) EOS. The emission of the kilonovae in the KLF bright components is mostly due to the radiation from the wind ejecta by the remnant discs of BNS mergers, while the emission of the kilonovae in the KLF faint components is mostly due to the radiation from the dynamical ejecta by the BNS mergers. △ Less

Submitted 13 April, 2023; v1 submitted 12 April, 2023; originally announced April 2023.

Comments: 28 pages, 16 figures, to appear in MNRAS

Modeling meteorite craters by impacting melted tin on sand

Authors: H. Y. Huang, P. R. Tsai, C. Y. Lu, H. Hau, Y. L. Chen, Z. T. Ling, Y. R. Wu, Tzay-Ming Hong

Abstract: To simulate the heated exterior of a meteorite, we impact a granular bed with melted tin. The morphology of tin remnant and crater is found to be sensitive to the temperature and solidification of tin. By employing deep learning and convolutional neural network, we can quantify and map the complex impact patterns onto network systems based on feature maps and Grad-CAM results. This gives us unprec… ▽ More To simulate the heated exterior of a meteorite, we impact a granular bed with melted tin. The morphology of tin remnant and crater is found to be sensitive to the temperature and solidification of tin. By employing deep learning and convolutional neural network, we can quantify and map the complex impact patterns onto network systems based on feature maps and Grad-CAM results. This gives us unprecedented details on how the projectile deforms and interacts with the granules, which information can be used to trace the development of different remnant shapes. Furthermore, full dynamics of granular system is revealed by the use of Particle Image Velocimetry. Kinetic energy, temperature and diameter of the projectile are used to build phase diagrams for the morphology of both crater and tin remnant. In addition to successfully reproducing key features of simple and complex craters, we are able to detect a possible artifact when compiling crater data from field studies. The depth of craters from high-energy impacts in our work is found to be independent of their width. However, when mixing data from different energy, temperature and diameter of projectile, a bogus power-law relationship appears between them. Like other controlled laboratory researches, our conclusions have the potential to benefit the study of paint in industry and asteroid sampling missions on the surface of celestial bodies. △ Less

Submitted 31 March, 2023; originally announced March 2023.

Comments: 6 pages, 5 figures

X-ray morphology due to charge-exchange emissions used to study the global structure around Mars

Authors: G. Y. Liang, T. R. Sun, H. Y. Lu, X. L. Zhu, Y. Wu, S. B. Li, H. G. Wei, D. W. Yuan, W. Cui, X. W. Ma, G. Zhao

Abstract: Soft x-ray emissions induced by solar wind ions that collide with neutral material in the solar system have been detected around planets, and were proposed as a remote probe for the solar wind interaction with the Martian exosphere. A multi-fluid three-dimensional magneto-hydrodynamic model is adopted to derive the global distributions of solar wind particles. Spherically symmetric exospheric H, H… ▽ More Soft x-ray emissions induced by solar wind ions that collide with neutral material in the solar system have been detected around planets, and were proposed as a remote probe for the solar wind interaction with the Martian exosphere. A multi-fluid three-dimensional magneto-hydrodynamic model is adopted to derive the global distributions of solar wind particles. Spherically symmetric exospheric H, H$_2$, He, O, and CO$_2$ density profiles and a sophisticated hybrid model that includes charge-exchange and proton/neutral excitation processes are used to study the low triplet line ratio $G=\frac{i+f}{r}$ (0.77$\pm$0.58) of O VII and total x-ray luminosity around Mars. We further calculate the emission factor $αあるふぁ$-value with different neutrals over a wide ion abundance and velocity ranges. Our results are in good agreement with those of previous reports. The evolution of the charge stage of solar wind ions shows that sequential recombination due to charge-exchange can be negligible at the interaction region. This only appears below the altitude of 400~km. The anonymous low disk $G$ ratio can be easily explained by the collisional quenching effect at neutral densities higher than 10$^{11}$cm$^{-3}$. However, the quenching contribution is small in Mars' exosphere and only appears below 400~km. Charge-exchange with H$_2$ and N$_2$ is still the most likely reason for this low $G$-ratio. X-ray emissivity maps in collisions with different neutrals differ from each other. A clear bow shock in the collision with all the neutrals is in accordance with previous reports. The resulting total x-ray luminosity of 6.55~MW shows a better agreement with the XMM-Newton observation of 12.8$\pm$1.4~MW than that of previous predictions. △ Less

Submitted 19 December, 2022; originally announced December 2022.

Comments: 17 pages, 13 figures; Accpted by ApJ

Type I X-ray bursts' spectra and fuel composition from the atoll and transient source 4U 1730$-$22

Authors: Yongqi Lu, Zhaosheng Li, Yuanyue Pan, Wenhui Yu, Yupeng Chen, Long Ji, Mingyu Ge, Shu Zhang, Jinlu Qu, Liming Song, Maurizio Falanga

Abstract: NICER observed two outbursts from the neutron star low-mass X-ray binary 4U~1730$-$22 in 2021 and 2022, which showed a similar spectral evolution in the hardness-intensity diagram. Seventeen type I X-ray bursts were identified in both outbursts. The X-ray burst spectra showed clear deviations from the blackbody model, firstly $\sim10$ s after onset. Adding the enhanced persistent emission due to t… ▽ More NICER observed two outbursts from the neutron star low-mass X-ray binary 4U~1730$-$22 in 2021 and 2022, which showed a similar spectral evolution in the hardness-intensity diagram. Seventeen type I X-ray bursts were identified in both outbursts. The X-ray burst spectra showed clear deviations from the blackbody model, firstly $\sim10$ s after onset. Adding the enhanced persistent emission due to the Poynting-Robertson drag or the reflection from the accretion disk both significantly improved the fitting results. We found that 12 out of 17 X-ray bursts showed the photospheric radius expansion (PRE) characteristic. Considering the nine PRE bursts out of ten X-ray bursts observed by Insight-HXMT, 78\% of bursts from 4U~1730--22 exhibited PRE. According to the burst rise time, the duration, the local accretion rate, and the burst fuel composition estimated from recurrence time, we propose that these PRE bursts were powered by pure helium. From the touchdown flux of PRE bursts, we estimate the source distance of $ d=7.54\pm{0.46} (X=0)$ kpc for a canonical neutron star with $M_{\rm NS}=1.4M_\odot$ and $R_{\rm NS}=10~{\rm km}$. △ Less

Submitted 21 March, 2023; v1 submitted 14 December, 2022; originally announced December 2022.

Comments: 11 pages, 6 figures; Matches the version published in A&A

Journal ref: A&A 670, A87 (2023)

There is No Place Like Home -- Finding Birth Radii of Stars in the Milky Way

Authors: Yuxi Lu, Ivan Minchev, Tobias Buck, Sergey Khoperskov, Matthias Steinmetz, Noam Libeskind, Gabriele Cescutti, Ken C. Freeman, Bridget Ratcliffe

Abstract: Stars move away from their birthplaces over time via a process known as radial migration, which blurs chemo-kinematic relations used for reconstructing the Milky Way (MW) formation history. To understand the true time evolution of the MW, one needs to take into account the effects of this process. We show that stellar birth radii can be derived directly from the data with minimum prior assumptions… ▽ More Stars move away from their birthplaces over time via a process known as radial migration, which blurs chemo-kinematic relations used for reconstructing the Milky Way (MW) formation history. To understand the true time evolution of the MW, one needs to take into account the effects of this process. We show that stellar birth radii can be derived directly from the data with minimum prior assumptions on the Galactic enrichment history. This is done by first recovering the time evolution of the stellar birth metallicity gradient, $d\mathrm{[Fe/H]}(R, τたう)/dR$, through its inverse relation to the metallicity range as a function of age today, allowing us to place any star with age and metallicity measurements back to its birthplace, $R_b$. Applying our method to a large, high-precision data set of MW disk subgiant stars, we find a steepening of the birth metallicity gradient from 11 to 8 Gyr ago, which coincides with the time of the last massive merger, Gaia-Sausage-Enceladus (GSE). This transition appears to play a major role in shaping both the age-metallicity relation and the bimodality in the [$αあるふぁ$/Fe]-[Fe/H] plane. By dissecting the disk into mono-$R_b$ populations, clumps in the low-[$αあるふぁ$/Fe] sequence appear, which are not seen in the total sample and coincide in time with known star-formation bursts, possibly associated with the Sagittarius Dwarf Galaxy. We estimated that the Sun was born at $4.5\pm 0.4$~kpc from the Galactic center. Our $R_b$ estimates provide the missing piece needed to recover the Milky Way formation history. △ Less

Submitted 22 May, 2024; v1 submitted 8 December, 2022; originally announced December 2022.

Constraining the origin of stellar binary black hole mergers by detections of their lensed host galaxies and gravitational wave signals

Authors: Zhiwei Chen, Youjun Lu, Yuetong Zhao

Abstract: A significant number of stellar binary black hole (sBBH) mergers may be lensed and detected by the third generation gravitational wave (GW) detectors. Their lensed host galaxies may be detectable, which thus helps to accurately localize these sources and provide a new approach to study the origin of sBBHs. In this paper, we investigate the detectability of the lensed host galaxies for the lensed s… ▽ More A significant number of stellar binary black hole (sBBH) mergers may be lensed and detected by the third generation gravitational wave (GW) detectors. Their lensed host galaxies may be detectable, which thus helps to accurately localize these sources and provide a new approach to study the origin of sBBHs. In this paper, we investigate the detectability of the lensed host galaxies for the lensed sBBH mergers. We find that the detection fraction of the host galaxies to the lensed GW events can be significantly different for a survey with a given limiting magnitude if sBBHs are produced by different mechanisms, such as the evolution of massive binary stars, the dynamical interactions in dense star clusters, and that assisted by active galactic nuclei or massive black holes. Furthermore, we illustrate that the statistical spatial distribution of those lensed sBBHs in its hosts resulting from different sBBH formation channels can be different from each other. Therefore, with the third generation GW detectors and future large scale galaxy surveys, it is possible to independently constrain the sBBH origin via the detection fraction of those lensed events with identifiable lensing host signatures and/or even constrain the contribution fractions from different sBBH formation mechanisms △ Less

Submitted 18 October, 2022; originally announced October 2022.

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

Bridging the gap -- the disappearance of the intermediate period gap for fully convective stars, uncovered by new ZTF rotation periods

Authors: Yuxi Lu, Jason L. Curtis, Ruth Angus, Trevor J. David, Soichiro Hattori

Abstract: The intermediate period gap, discovered by Kepler, is an observed dearth of stellar rotation periods in the temperature-period diagram at $\sim$ 20 days for G dwarfs and up to $\sim$ 30 days for early-M dwarfs. However, because Kepler mainly targeted solar-like stars, there is a lack of measured periods for M dwarfs, especially those at the fully convective limit. Therefore it is unclear if the in… ▽ More The intermediate period gap, discovered by Kepler, is an observed dearth of stellar rotation periods in the temperature-period diagram at $\sim$ 20 days for G dwarfs and up to $\sim$ 30 days for early-M dwarfs. However, because Kepler mainly targeted solar-like stars, there is a lack of measured periods for M dwarfs, especially those at the fully convective limit. Therefore it is unclear if the intermediate period gap exists for mid- to late-M dwarfs. Here, we present a period catalog containing 40,553 rotation periods (9,535 periods $>$ 10 days), measured using the Zwicky Transient Facility (ZTF). To measure these periods, we developed a simple pipeline that improves directly on the ZTF archival light curves and reduces the photometric scatter by 26%, on average. This new catalog spans a range of stellar temperatures that connect samples from Kepler with MEarth, a ground-based time domain survey of bright M-dwarfs, and reveals that the intermediate period gap closes at the theoretically predicted location of the fully convective boundary ($G_{\rm BP} - G_{\rm RP} \sim 2.45$ mag). This result supports the hypothesis that the gap is caused by core-envelope interactions. Using gyro-kinematic ages, we also find a potential rapid spin-down of stars across this period gap. △ Less

Submitted 12 October, 2022; originally announced October 2022.

On Detecting Nearby Nano-Hertz Gravitational Wave Sources via Pulsar Timing Arrays

Authors: Xiao Guo, Youjun Lu, Qingjuan Yu

Abstract: Massive binary black holes (MBBHs) in nearby galactic centers, if any, may be nano-Hertz gravitational wave (GW) sources for pulsar timing arrays (PTAs) to detect. Normally the objective GWs for PTA experiments are approximated as plane waves because its sources are presumably located faraway. For nearby GW sources, however, this approximation may be inaccurate due to the curved GW wave front and… ▽ More Massive binary black holes (MBBHs) in nearby galactic centers, if any, may be nano-Hertz gravitational wave (GW) sources for pulsar timing arrays (PTAs) to detect. Normally the objective GWs for PTA experiments are approximated as plane waves because its sources are presumably located faraway. For nearby GW sources, however, this approximation may be inaccurate due to the curved GW wave front and the GW strength changes along the paths of PTA pulsar pulses. In this paper, we analyze the near-field effect in the PTA detection of nearby sources and find it is important if the source distance is less than a few tens Mpc, and ignoring this effect may lead to a significant signal-to-noise underestimation especially when the source distance is comparable to the pulsar distances. As examples, we assume a nano-Hertz MBBH source located at either the Galactic Center (GC) or the Large Magellanic Cloud (LMC) according to the observational constraints/hints on the MBBH parameter space, and estimate its detectability by current/future PTAs. We find that the GC MBBH may be detectable by the Square Kilometer Array (SKA) PTA. It is challenging for detecting the LMC MBBH; however, if a number ($N\gtrsim10$) of stable millisecond pulsars can be found in the LMC center, the MBBH may be detectable via a PTA formed by these pulsars. We further illustrate the near-field effects on the PTA detection of an isotropic GW background contributed mainly by nearby GW sources, and the resulting angular correlation is similar to the Hellings-Downs curve. △ Less

Submitted 13 October, 2022; v1 submitted 12 September, 2022; originally announced September 2022.

Comments: 24 pages, 11 figures, accepted for publication in The Astrophysical Journal

High spectral-resolution interferometry down to 1 micron with Asgard/BIFROST at VLTI: Science drivers and project overview

Authors: Stefan Kraus, Daniel Mortimer, Sorabh Chhabra, Yi Lu, Isabelle Codron, Tyler Gardner, Narsireddy Anugu, John Monnier, Jean-Baptiste Le Bouquin, Michael Ireland, Frantz Martinache, Denis Defrère, Marc-Antoine Martinod

Abstract: We present science cases and instrument design considerations for the BIFROST instrument that will open the short-wavelength (Y/J/H-band), high spectral dispersion (up to R=25,000) window for the VLT Interferometer. BIFROST will be part of the Asgard Suite of instruments and unlock powerful venues for studying accretion & mass-loss processes at the early/late stages of stellar evolution, for detec… ▽ More We present science cases and instrument design considerations for the BIFROST instrument that will open the short-wavelength (Y/J/H-band), high spectral dispersion (up to R=25,000) window for the VLT Interferometer. BIFROST will be part of the Asgard Suite of instruments and unlock powerful venues for studying accretion & mass-loss processes at the early/late stages of stellar evolution, for detecting accreting protoplanets around young stars, and for probing the spin-orbit alignment in directly-imaged planetary systems and multiple star systems. Our survey on GAIA binaries aims to provide masses and precision ages for a thousand stars, providing a legacy data set for improving stellar evolutionary models as well as for Galactic Archaeology. BIFROST will enable off-axis spectroscopy of exoplanets in the 0.025-1" separation range, enabling high-SNR, high spectral resolution follow-up of exoplanets detected with ELT and JWST. We give an update on the status of the project, outline our key technology choices, and discuss synergies with other instruments in the proposed Asgard Suite of instruments. △ Less

Submitted 9 August, 2022; originally announced August 2022.

Comments: 19 pages, 7 figures, SPIE 2022 "Astronomical Telescopes and Instrumentation" manuscript 12183-66

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

Authors: Xiao Guo, Youjun Lu

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

Submitted 1 July, 2022; originally announced July 2022.

Comments: 12 pages, 10 figures, accepted for publication in PRD

WS-Snapshot: An effective algorithm for wide-field and large-scale imaging

Authors: Yangfan Xie, Feng Wang, Hui Deng, Ying Mei, Ying-He Celeste Lu, Gabriella Hodosan, Vladislav Stolyarov, Oleg Smirnov, Xiaofeng Li, Tim Cornwell

Abstract: The Square Kilometre Array (SKA) is the largest radio interferometer under construction in the world. The high accuracy, wide-field and large size imaging significantly challenge the construction of the Science Data Processor (SDP) of SKA. We propose a hybrid imaging method based on improved W-Stacking and snapshots. The w range is reduced by fitting the snapshot $uv$ plane, thus effectively enhan… ▽ More The Square Kilometre Array (SKA) is the largest radio interferometer under construction in the world. The high accuracy, wide-field and large size imaging significantly challenge the construction of the Science Data Processor (SDP) of SKA. We propose a hybrid imaging method based on improved W-Stacking and snapshots. The w range is reduced by fitting the snapshot $uv$ plane, thus effectively enhancing the performance of the improved W-Stacking algorithm. We present a detailed implementation of WS-Snapshot. With full-scale SKA1-LOW simulations, we present the imaging performance and imaging quality results for different parameter cases. The results show that the WS-Snapshot method enables more efficient distributed processing and significantly reduces the computational time overhead within an acceptable accuracy range, which would be crucial for subsequent SKA science studies. △ Less

Submitted 1 July, 2022; originally announced July 2022.

Comments: 10 pages, 10 figures, 6 tables, accepted by MNRAS

Meta-analysis of photometric and asteroseismic measurements of stellar rotation periods: the Lomb-Scargle periodogram, autocorrelation function, wavelet and rotational splitting analysis for 92 Kepler asteroseismic targets

Authors: Yuting Lu, Othman Benomar, Shoya Kamiaka, Yasushi Suto

Abstract: We perform photometric (the Lomb-Scargle periodogram, autocorrelation, and wavelet) and asteroseismic analyses of 92 Kepler solar-like main-sequence stars to understand the reliability of the measured stellar rotation periods. We focus on the 70 stars without reported stellar companions, and classify them into four groups according to the quarter-to-quarter variance of the Lomb-Scargle period and… ▽ More We perform photometric (the Lomb-Scargle periodogram, autocorrelation, and wavelet) and asteroseismic analyses of 92 Kepler solar-like main-sequence stars to understand the reliability of the measured stellar rotation periods. We focus on the 70 stars without reported stellar companions, and classify them into four groups according to the quarter-to-quarter variance of the Lomb-Scargle period and the precision of the asteroseismic period. We present detailed individual comparison among photometric and asteroseismic constraints for these stars. We find that most of our targets exhibit significant quarter-to-quarter variances in the photometric periods, suggesting that the photometrically estimated period should be regarded as a simplified characterization of the true stellar rotation period, especially under the presence of the latitudinal differential rotation. On the other hand, there are a fraction of stars with a relatively small quarter-to-quarter variance in the photometric periods, most of which have consistent values for asteroseismically and photometrically estimated rotation periods. We also identify over ten stars whose photometric and asteroseismic periods significantly disagree, which would be potentially interesting targets for further individual investigations. △ Less

Submitted 25 October, 2022; v1 submitted 10 June, 2022; originally announced June 2022.

Discovery of a 584.65 Hzへるつ Burst Oscillation in the Low Mass X-Ray Binary 4U 1730--22

Authors: Zhaosheng Li, Wenhui Yu, Yongqi Lu, Yuanyue Pan, Maurizio Falanga

Abstract: Type-I X-ray burst oscillations are powered by thermonuclear released on the neutron star (NS) surface in low mass X-ray binaries (LMXBs), where the burst oscillation frequencies are close to the NS spin rates. In this work, we report the detection of oscillation at 584.65 Hzへるつ during the cooling tail of a type-I X-ray bursts observed from the accreting NS LMXB 4U~1730--22 in 2022 March 20, by the \… ▽ More Type-I X-ray burst oscillations are powered by thermonuclear released on the neutron star (NS) surface in low mass X-ray binaries (LMXBs), where the burst oscillation frequencies are close to the NS spin rates. In this work, we report the detection of oscillation at 584.65 Hzへるつ during the cooling tail of a type-I X-ray bursts observed from the accreting NS LMXB 4U~1730--22 in 2022 March 20, by the \textit{ Neutron star Interior Composition Explorer} (\textit{NICER}) telescope. The oscillation signal showed a strong Leahy power, $P_{\rm m}\sim54.04$, around 584.65 Hzへるつ, which has single trial and multiple trials confidence levels of $7.05σしぐま$ and $4.78σしぐま$, respectively. The folded pulse profile of the oscillation in the 0.2--10 keV band showed a sinusoidal shape with the fraction amplitude rms of $(12.5\pm1.8)\%$. We found the oscillation frequency showed insignificant upward drifting, i.e., less than 0.3 Hzへるつ, during the cooling tail, similar as the behavior appearing in accreting millisecond X-ray pulsars (AMXP), and indicate the source could be an AMXP spinning at 1.71 ms. △ Less

Submitted 7 September, 2022; v1 submitted 24 May, 2022; originally announced May 2022.

Comments: matches published version

Journal ref: ApJ, 935, 123, 2022

An implementation of Galactic white dwarf binary data analysis for MLDC-3.1

Authors: Yang Lu, En-Kun Li, Yi-Ming Hu, Jian-dong Zhang, Jianwei Mei

Abstract: The space-borne gravitational wave detectors will observe a large population of double white dwarf binaries in the Milky Way. However, the search for double white dwarfs in the gravitational wave data will be time-consuming due to the large number of templates involved and antenna response calculation. In this paper, we implement an iterative combinatorial algorithm to search for double white dwar… ▽ More The space-borne gravitational wave detectors will observe a large population of double white dwarf binaries in the Milky Way. However, the search for double white dwarfs in the gravitational wave data will be time-consuming due to the large number of templates involved and antenna response calculation. In this paper, we implement an iterative combinatorial algorithm to search for double white dwarfs in MLDC-3.1 data. To quickly determine the rough parameters of the target sources, the following algorithms are adopted in a coarse search process: (1) using the downsampling method to reduce the number of original data points; (2) using the undersampling method to speed up the generation of a single waveform template; (3) using the stochastic template bank method to quickly construct the waveform template bank while achieving high coverage of the parameter space; (4) Combining the FFT acceleration algorithm with the stochastic template bank to reduce the calculation time of a single template. A fine search process is applied to further determine the parameters of the signals based on the coarse search, for which we adopt the particle swarm optimization. Finally, we detected $\mathcal{O}(10^4)$ double white dwarf signals, validating the feasibility of our method. △ Less

Submitted 7 February, 2023; v1 submitted 4 May, 2022; originally announced May 2022.

Comments: 16 pages, 15 figures

Reliability and limitations of inferring birth radii in the Milky Way disk

Authors: Yuxi Lu, Tobias Buck, Ivan Minchev, Melissa K. Ness

Abstract: Recovering the birth radii of observed stars in the Milky Way is one of the ultimate goals of Galactic Archaeology. One method to infer the birth radius and the evolution of the ISM metallicity assumes a linear relation between the ISM metallicity with radius at any given look-back time. Here we test the reliability of this assumption by using 4 zoom-in cosmological hydrodynamic simulations from t… ▽ More Recovering the birth radii of observed stars in the Milky Way is one of the ultimate goals of Galactic Archaeology. One method to infer the birth radius and the evolution of the ISM metallicity assumes a linear relation between the ISM metallicity with radius at any given look-back time. Here we test the reliability of this assumption by using 4 zoom-in cosmological hydrodynamic simulations from the NIHAO-UHD project. We find that one can infer precise birth radii only when the stellar disk starts to form, which for our modeled galaxies happens ~ 10 Gyr ago, in agreement with recent estimates for the Milky Way. At later times the linear correlation between the ISM metallicity and radius increases, as stellar motions become more ordered and the azimuthal variations of the ISM metallicity start to drop. The formation of a central bar and perturbations from mergers can increase this uncertainty in the inner and outer disk, respectively. △ Less

Submitted 2 May, 2022; v1 submitted 30 April, 2022; originally announced May 2022.

Further Evidence of Modified Spin-down in Sun-like Stars: Pileups in the Temperature-Period Distribution

Authors: Trevor J. David, Ruth Angus, Jason L. Curtis, Jennifer L. van Saders, Isabel L. Colman, Gabriella Contardo, Yuxi Lu, Joel C. Zinn

Abstract: We combine stellar surface rotation periods determined from NASA's Kepler mission with spectroscopic temperatures to demonstrate the existence of pileups at the long-period and short-period edges of the temperature-period distribution for main-sequence stars with temperatures exceeding $\sim 5500$K. The long-period pileup is well-described by a curve of constant Rossby number, with a critical valu… ▽ More We combine stellar surface rotation periods determined from NASA's Kepler mission with spectroscopic temperatures to demonstrate the existence of pileups at the long-period and short-period edges of the temperature-period distribution for main-sequence stars with temperatures exceeding $\sim 5500$K. The long-period pileup is well-described by a curve of constant Rossby number, with a critical value of $\mathrm{Ro_{crit}} \lesssim 2$. The long-period pileup was predicted by van Saders et al. (2019) as a consequence of weakened magnetic braking, in which wind-driven angular momentum losses cease once stars reach a critical Rossby number. Stars in the long-period pileup are found to have a wide range of ages ($\sim 2-6$Gyr), meaning that, along the pileup, rotation period is strongly predictive of a star's surface temperature but weakly predictive of its age. The short-period pileup, which is also well-described by a curve of constant Rossby number, is not a prediction of the weakened magnetic braking hypothesis but may instead be related to a phase of slowed surface spin-down due to core-envelope coupling. The same mechanism was proposed by Curtis et al. (2020) to explain the overlapping rotation sequences of low-mass members of differently aged open clusters. The relative dearth of stars with intermediate rotation periods between the short- and long-period pileups is also well-described by a curve of constant Rossby number, which aligns with the period gap initially discovered by McQuillan et al. (2013a) in M-type stars. These observations provide further support for the hypothesis that the period gap is due to stellar astrophysics, rather than a non-uniform star-formation history in the Kepler field. △ Less

Submitted 10 May, 2022; v1 submitted 16 March, 2022; originally announced March 2022.

Comments: Accepted to ApJ. 29 pages, 21 figures. The data and code required to reproduce this work is available at http://github.com/trevordavid/rossby-ridge