Showing 1–50 of 125 results for author: Tang, S

Rich and diverse molecular gas environments of closely-separated dual quasars viewed by ALMA

Authors: Shenli Tang, John D. Silverman, Zhaoxuan Liu, Manda Banerji, Tomoko Suzuki, Seiji Fujimoto, Andy Goulding, Masatoshi Imanishi, Toshihiro Kawaguchi, Connor Bottrell, Tilman Hartwig, Knud Jahnke, Masafusa Onoue, Malte Schramm, Yoshihiro Ueda

Abstract: We present a study of the molecular gas in five closely-spaced ($R_{\perp}<20$ kpc) dual quasars ($L_{\rm bol}\gtrsim10^{44}~\mathrm{erg~s}^{-1}$) at redshifts $0.4<z<0.8$ with the Atacama Large Millimeter/submillimeter Array. The dual quasar phase represents a distinctive stage during the interaction between two galaxies for investigating quasar fueling and feedback effects on the gas reservoir.… ▽ More We present a study of the molecular gas in five closely-spaced ($R_{\perp}<20$ kpc) dual quasars ($L_{\rm bol}\gtrsim10^{44}~\mathrm{erg~s}^{-1}$) at redshifts $0.4<z<0.8$ with the Atacama Large Millimeter/submillimeter Array. The dual quasar phase represents a distinctive stage during the interaction between two galaxies for investigating quasar fueling and feedback effects on the gas reservoir. The dual quasars were selected from the Sloan Digital Sky Survey and Subaru/Hyper Suprime-Cam Subaru Strategic Program, with confirmatory spectroscopic validation. Based on the detection of the CO J=2--1 emission line with Band 4, we derived key properties including CO luminosities, line widths, and molecular gas masses for these systems. Among the ten quasars of the five pairs, eight have line detections exceeding $5σ$. The detected sources prominently harbor substantial molecular gas reservoirs, with molecular gas masses ($M_{\text{molgas}}$) between $10^{9.6-10.5}~\mathrm{M_{\odot}}$, and molecular gas-to-stellar mass ratios ($μ_{\text{molgas}}$) spanning $18-97\%$. The overall $μ_{\text{molgas}}$ of these dual quasars agrees with that of inactive star-forming main-sequence galaxies at comparable redshifts, indicating no clear evidence of quenching. However, intriguing features in each individual system show possible evidence of AGN feedback, matter transfer, and compaction processes. △ Less

Submitted 12 July, 2024; originally announced July 2024.

Measuring the Spot Variability of T Tauri Stars Using Near-IR Atomic Fe and Molecular OH Lines

Authors: Shih-Yun Tang, Christopher M. Johns-Krull, L. Prato, Asa G. Stahl

Abstract: As part of the Young Exoplanets Spectroscopic Survey (YESS), this study explores the spot variability of 13 T Tauri Stars (TTSs) in the near-infrared $H$ band, using spectra from the Immersion GRating INfrared Spectrometer (IGRINS). By analyzing effective temperature ($T_{\rm eff}$) sensitive lines of atomic FeI at ~1.56259 um and ~1.56362 um, and molecular OH at ~1.56310 um and ~1.56317 um, we de… ▽ More As part of the Young Exoplanets Spectroscopic Survey (YESS), this study explores the spot variability of 13 T Tauri Stars (TTSs) in the near-infrared $H$ band, using spectra from the Immersion GRating INfrared Spectrometer (IGRINS). By analyzing effective temperature ($T_{\rm eff}$) sensitive lines of atomic FeI at ~1.56259 um and ~1.56362 um, and molecular OH at ~1.56310 um and ~1.56317 um, we develop an empirical equivalent width ratio (EWR) relationship for $T_{\rm eff}$ in the range of 3400-5000 K. This relationship allows for precise relative $T_{\rm eff}$ estimates to within tens of Kelvin and demonstrates compatibility with solar metallicity target models. However, discrepancies between observational data and model predictions limit the extension of the $T_{\rm eff}$-EWR relationship to a broader parameter space. Our study reveals that both classical and weak-line TTSs can exhibit $T_{\rm eff}$ variations exceeding 150 K over a span of two years. The detection of a quarter-phase delay between the EWR and radial velocity phase curves in TTSs indicates spot-driven signals. A phase delay of 0.06 $\pm$ 0.13 for CI Tau, however, suggests additional dynamics, potentially caused by planetary interaction, inferred from a posited 1:1 commensurability between the rotation period and orbital period. Moreover, a positive correlation between $T_{\rm eff}$ variation amplitude and stellar inclination angle support the existence of high-latitude spots on TTSs, further enriching our understanding of stellar surface activity in young stars. △ Less

Submitted 18 June, 2024; originally announced June 2024.

Comments: 31 pages, 19 figures, and 5 tables. Accepted to ApJ

Multi-spectral sirens: Gravitational-wave cosmology with (multi-) sub-populations of binary black holes

Authors: Yin-Jie Li, Shao-Peng Tang, Yuan-Zhu Wang, Yi-Zhong Fan

Abstract: The cosmic expansion rate can be directly measured with gravitational waves (GWs) of the compact binary mergers, by jointly constraining the mass function of the population and the cosmological model via the so called spectral sirens. Such a method relies on the features in the mass functions, which may originate from some individual sub-populations, and hence become blurred/indistinct due to the… ▽ More The cosmic expansion rate can be directly measured with gravitational waves (GWs) of the compact binary mergers, by jointly constraining the mass function of the population and the cosmological model via the so called spectral sirens. Such a method relies on the features in the mass functions, which may originate from some individual sub-populations, and hence become blurred/indistinct due to the superposition of different sub-populations. In this work we propose a novel approach to constrain the cosmic expansion rate with sub-populations of GW events, named multi-spectral sirens. We illustrate the advantage of the multi-spectral sirens compared to the traditional spectral sirens by simulation with mock data. The application of this approach to the GWTC-3 data yields $H_0=73.25^{+29.87}_{-25.55}~{\rm Mpc}^{-1}~{\rm km}~{\rm s}^{-1}$ (median and symmetric 68.3\% credible level). The incorporation of the bright standard siren GW170817 with a uniform prior in $ [10,200] ~{\rm Mpc}^{-1}~{\rm km}~{\rm s}^{-1}$ gives $H_0=72.38^{+15.03}_{-9.13}~{\rm Mpc}^{-1}~{\rm km}~{\rm s}^{-1}$ (68.3\% confidence level), corresponding to an improvement of $\sim28\%$ with respect to the measurement from sole GW170817. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: 14 pages, 8 figures, 1 table; comments are welcome

Digging into the ultraviolet luminosity functions of galaxies at high redshifts: galaxies evolution, reionization, and cosmological parameters

Authors: Yi-Ying Wang, Lei Lei, Shao-Peng Tang, Guan-Wen Yuan, Yi-Zhong Fan

Abstract: Thanks to the successful performance of the James Webb Space Telescope, our understanding of the epoch of reionization of the Universe has been advanced. The ultraviolet luminosity functions (UV LFs) of galaxies span a wide range of redshift, not only revealing the connection between galaxies and dark matter (DM) halos but also providing the information during reionization. In this work, we develo… ▽ More Thanks to the successful performance of the James Webb Space Telescope, our understanding of the epoch of reionization of the Universe has been advanced. The ultraviolet luminosity functions (UV LFs) of galaxies span a wide range of redshift, not only revealing the connection between galaxies and dark matter (DM) halos but also providing the information during reionization. In this work, we develop a model connecting galaxy counts and apparent magnitude based on UV LFs, which incorporates redshift-dependent star formation efficiency (SFE) and corrections for dust attenuation. By synthesizing some observations across the redshift range $4\le z \le 10$ from various galaxy surveys, we discern the evolving SFE with increasing redshift and DM halo mass through model fitting. Subsequent analyses indicate that the Thomson scattering optical depth was $τ_{\rm e} = 0.052^{+0.003}_{-0.002}$ and the epoch of reionization started (ended) at $z=20.58^{+6.25}_{-6.75}$ ($z=5.38^{+0.65}_{-0.70}$) which is insensitive to the choice of the truncated magnitude of the UV LFs. Incorporating additional dataset and some reasonable constraints, the amplitude of matter perturbation is found to be $σ_8=0.79\pm0.05$, which is consistent with the standard $Λ$CDM model. Future galaxy surveys and the dynamical simulations of galaxy evolution will break the degeneracy between SFE and cosmological parameters, improving the accuracy and the precision of the UV LF model further. △ Less

Submitted 15 May, 2024; originally announced May 2024.

Comments: 15 pages, 9 figures and 3 tables

Exploring field-evolution and dynamical-capture coalescing binary black holes in GWTC-3

Authors: Yin-Jie Li, Shao-Peng Tang, Shi-Jie Gao, Dao-Cheng Wu, Yuan-Zhu Wang

Abstract: The continuously expanding sample of gravitational-wave observations is revealing the formation and evolutionary mechanism of merging compact binaries. Two primary channels, namely, isolated field binary evolution and dynamical capture, are widely accepted as potential producers of merging binary black holes (BBHs), which are distinguishable with the spin-orientation distributions of the BBHs. We… ▽ More The continuously expanding sample of gravitational-wave observations is revealing the formation and evolutionary mechanism of merging compact binaries. Two primary channels, namely, isolated field binary evolution and dynamical capture, are widely accepted as potential producers of merging binary black holes (BBHs), which are distinguishable with the spin-orientation distributions of the BBHs. We investigate the two formation channels in GWTC-3, with a dedicated semi-parametric population model, i.e., a mixture of two sub-populations with different spin-orientation distributions (one is nearly-aligned and the other is nearly-isotropic). It turns out that the two sub-populations have different mass and mass-ratio distributions. The nearly-aligned sub-population, which is consistent with the isolated field formation channels, has a less preference for symmetric systems, and likely dominate the 10-solar-mass peak in the primary-mass function. While the isotropic sub-population shows a stronger preference for symmetric systems, and mainly contribute to the 35-solar-mass peak in the primary-mass function, consistent with the dynamical channels. Moreover, our results show that the purely isotropic-spin and the single well-aligned (i.e., the width of $\cosθ$ distribution $σ_{\rm t}<0.5$) scenario are ruled out (by a Bayes factor of $\ln\mathcal{B}=5.2$ and $\ln\mathcal{B}=9.8$). △ Less

Submitted 15 April, 2024; originally announced April 2024.

Comments: 26 pages, 16 figures

Upper Limit of Sound Speed in Nuclear Matter: A Harmonious Interplay of Transport Calculation and Perturbative QCD Constraint

Authors: Shao-Peng Tang, Yong-Jia Huang, Ming-Zhe Han, Yi-Zhong Fan

Abstract: Very recently, it has been shown that there is an upper bound on the squared sound speed of nuclear matter from the transport, which reads $c_{\rm s}^2 \leq 0.781$ \citep{2024arXiv240214085H}. In this work, we demonstrate that this upper bound is corroborated by the reconstructed equation of state (EOS) for ultra-dense matter. The reconstruction integrates multi-messenger observation for neutron s… ▽ More Very recently, it has been shown that there is an upper bound on the squared sound speed of nuclear matter from the transport, which reads $c_{\rm s}^2 \leq 0.781$ \citep{2024arXiv240214085H}. In this work, we demonstrate that this upper bound is corroborated by the reconstructed equation of state (EOS) for ultra-dense matter. The reconstruction integrates multi-messenger observation for neutron stars (NSs), as well as the theoretical constraints, including that from chiral effective field theory and perturbative quantum chromodynamics (pQCD). In particular, the latest radius measurements for PSR J0437-4715 ($11.36^{+0.95}_{-0.63}$ km) and PSR J0030+0451 ($11.71^{+0.88}_{-0.83}$ km, in the ST+PDT model) by NICER have been adopted. The result shows in all cases, the $c_{\rm s}^2 \leq 0.781$ upper limit for EOS will naturally yield the properties of matter near center of the massive neutron star consistent with the causality-driven constraint from pQCD, where in practice, the constraint is applied at ten nuclear saturation density ($n_{\rm L}=10n_{\rm s}$). We also note that there is a strong correlation for the maximum square of sound speed $c_s^2$ with $n_{\rm L}$, and $c_{\rm s}^2 \leq 0.781$ is somehow violated when $n_{\rm L} = n_{\rm c,TOV}$. The result indicates that a higher density in implementing the pQCD constraint, even considering the uncertainties from statistics, is more natural. Moreover, the remarkable agreement between the outcomes derived from these two distinct and independent constraints (i.e., the transport calculation and pQCD boundary) lends strong support to their validity. △ Less

Submitted 15 April, 2024; originally announced April 2024.

Comments: 6+3 pages, 5+4 figures

Report number: RIKEN-iTHEMS-Report-24

Bulk properties of PSR J0030+0451 inferred with the compactness measurement of NICER

Authors: Chuan-Ning Luo, Shao-Peng Tang, Ming-Zhe Han, Jin-Liang Jiang, Wei-Hong Gao, Da-Ming Wei

Abstract: In 2019, Neutron star Interior Composition ExploreR (NICER) mission released its findings on the mass and radius of the isolated neutron star (INS) PSR J0030+0451, revealing a mass of approximately 1.4 solar masses ($M_{\odot}$) and a radius near 13 kilometers. However, the recent re-analysis by the NICER collaboration \citep{vinciguerra2024updated} suggests that the available data primarily yield… ▽ More In 2019, Neutron star Interior Composition ExploreR (NICER) mission released its findings on the mass and radius of the isolated neutron star (INS) PSR J0030+0451, revealing a mass of approximately 1.4 solar masses ($M_{\odot}$) and a radius near 13 kilometers. However, the recent re-analysis by the NICER collaboration \citep{vinciguerra2024updated} suggests that the available data primarily yields a precise inference of the compactness for this source while the resulting mass and radius are strongly model-dependent and diverse (the 68.3\% credible regions just overlap slightly for the ST+PDT and PDT-U models). By integrating this compactness data with the equation of state (EoS) refined by our latest investigations, we have deduced the mass and radius for PSR J0030+0451, delivering estimates of $M=1.48^{+0.09}_{-0.10}~M_\odot$ and $R=12.38_{-0.70}^{+0.51}~{\rm km}$ for the compactness found in ST+PDT model, alongside $M=1.47^{+0.14}_{-0.20}~M_\odot$ and $R=12.37_{-0.69}^{+0.50}~{\rm km}$ for the compactness in PDT-U model. These two groups of results are well consistent with each other and the direct X-ray data inference within the ST+PDT model seems to be favored. Additionally, we have calculated the tidal deformability, moment of inertia, and gravitational binding energy for this NS. Furthermore, employing these refined EoS models, we have updated mass-radius estimates for three INSs with established gravitational redshifts. △ Less

Submitted 2 April, 2024; v1 submitted 20 March, 2024; originally announced March 2024.

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

The Present-Day Mass Function of Star Clusters in the Solar Neighborhood

Authors: Xiaoying Pang, Siqi Liao, Jiadong Li, Zhiqiang Yan, Mingjie Jian, M. B. N. Kouwenhoven, Shih-Yun Tang, Yifan Wang

Abstract: This work analyses the present-day mass function (PDMF) of 93~star clusters utilizing Gaia DR3 data, with membership determined by the StarGo machine learning algorithm. The impact of unresolved binary systems on mass estimation is rigorously assessed, adopting three mass ratio profiles for correction. The PDMF is characterized by the power-law index, $α$, derived through a robust maximum likeliho… ▽ More This work analyses the present-day mass function (PDMF) of 93~star clusters utilizing Gaia DR3 data, with membership determined by the StarGo machine learning algorithm. The impact of unresolved binary systems on mass estimation is rigorously assessed, adopting three mass ratio profiles for correction. The PDMF is characterized by the power-law index, $α$, derived through a robust maximum likelihood method that avoids biases associated with data binning. The value of $α$ for stars between the completeness limited mass of Gaia with a mean 0.3 $M_\odot$ for our cluster samples and 2 $M_\odot$, exhibits stability for clusters younger than 200 Myr, decreasing for older clusters, particularly when considering stars within the half-mass radius. The PDMF of these star clusters is consistent with a dynamically evolved Kroupa IMF via the loss of low-mass stars. Cluster morphology shows a correlation with $α$, as $α$ values exhibit a decreasing trend from filamentary to tidal-tail clusters, mirroring the sequence of increasing cluster age. The dependence of $α$ on total cluster mass is weak, with a subtle increase for higher-mass clusters, especially outside the half-mass radius. We do not observe a correlation between $α$ and the mean metallicity of the clusters. Younger clusters have lower metallicity compared to their older counterparts, which indicates that the older clusters might have migrated to the solar neighbourhood from the inner disk. A comparison with numerical models incorporating a black hole population suggests the need for observations of distant, older, massive open clusters to determine whether or not they contain black holes. △ Less

Submitted 13 March, 2024; originally announced March 2024.

Comments: 19 pages, 10 figures, accepted to ApJ

The Bulk Properties of Isolated Neutron Stars Inferred from the Gravitational Redshift Measurements

Authors: Chuan-Ning Luo, Shao-Peng Tang, Jing-Liang Jiang, Wei-Hong Gao, Da-Ming Wei

Abstract: The measurements of the bulk properties of most isolated neutron stars (INSs) are challenging tasks. Tang et al. (2020) have developed a new method, based on the equation of state (EoS) of neutron star (NS) material constrained by the observational data, to infer the gravitational masses of a few INSs whose gravitational redshifts are available. However, in that work, the authors only considered t… ▽ More The measurements of the bulk properties of most isolated neutron stars (INSs) are challenging tasks. Tang et al. (2020) have developed a new method, based on the equation of state (EoS) of neutron star (NS) material constrained by the observational data, to infer the gravitational masses of a few INSs whose gravitational redshifts are available. However, in that work, the authors only considered the constraints on the EoS from nuclear experiments/theories and the gravitational wave data of GW170817; the possible phase transition has not been taken into account. In this work, we adopt three EoS models (including the one incorporates a first-order strong phase transition) that are constrained by the latest multimessenger NS data, including in particular the recent mass\textendash radius measurements of two NSs by Neutron Star Interior Composition Explorer, to update the estimation of the gravitational masses of RBS 1223, RX J0720.4-3125, and RX J1856.5-3754. In comparison to our previous approach, the new constraints are tighter, and the gravitational masses are larger by about $0.1M_\odot$. All the inferred gravitational masses are within the range of the NS masses measured in other ways. We have also calculated the radius, tidal-deformability, and moment of inertia of these sources. The inclusion of the first-order strong phase transition has little influence on modifying the results. △ Less

Submitted 8 March, 2024; v1 submitted 6 March, 2024; originally announced March 2024.

Comments: 8pages,4figures

Journal ref: ApJ 930 4 (2022)

Mass and radius of the most massive neutron star: The probe of the equation of state and perturbative QCD

Authors: Shao-Peng Tang, Ming-Zhe Han, Yong-Jia Huang, Yi-Zhong Fan, Da-Ming Wei

Abstract: Recently, an association of GW190425 and FRB 20190425A had been claimed and a highly magnetized neutron star (NS) remnant was speculated. Given the $\sim 2.5$-h delay of the occurrence of FRB 20190425A, a uniformly rotating supramassive magnetar is favored since the differential rotation would have been promptly terminated by the magnetic braking. The required maximum gravitational mass (… ▽ More Recently, an association of GW190425 and FRB 20190425A had been claimed and a highly magnetized neutron star (NS) remnant was speculated. Given the $\sim 2.5$-h delay of the occurrence of FRB 20190425A, a uniformly rotating supramassive magnetar is favored since the differential rotation would have been promptly terminated by the magnetic braking. The required maximum gravitational mass ($M_{\rm TOV}$) of the nonrotating NS is $\approx 2.77M_\odot$, which is strongly in tension with the relatively low $M_{\rm TOV}\approx 2.25M_\odot$ obtained in current equation of state (EOS) constraints incorporating perturbative quantum chromodynamics (pQCD) information. However, the current mass-radius and mass-tidal deformability measurements of NSs alone do not convincingly exclude the high $M_{\rm TOV}$ possibility. By performing EOS constraints with mock measurements, we find that with a $2\%$ determination for the radius of PSR J0740+6620-like NS it is possible to distinguish between the low and high $M_{\rm TOV}$ scenarios. We further explore the prospect to resolve the issue of the appropriate density to impose the pQCD constraints with future massive NS observations or determinations of $M_{\rm TOV}$ and/or $R_{\rm TOV}$. It turns out that measuring the radius of a PSR J0740+6620-like NS is insufficient to probe the EOSs around 5 nuclear saturation density, where the information from pQCD becomes relevant. The additional precise $M_{\rm TOV}$ measurements anyhow could provide insights into the EOS at such a density. Indeed, supposing the central engine of GRB 170817A is a black hole formed via the collapse of a supramassive NS, the resulting $M_{\rm TOV}\approx 2.2M_\odot$ considerably softens the EOS at the center of the most massive NS, which is in favor of imposing the pQCD constraint at density beyond the one achievable in the NSs. △ Less

Submitted 30 April, 2024; v1 submitted 22 November, 2023; originally announced November 2023.

Comments: 12 pages, 5 figures, PRD published

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

Measuring Mass and Radius of the Maximum-mass Nonrotating Neutron Star

Authors: Shao-Peng Tang, Bo Gao, Yin-Jie Li, Yi-Zhong Fan, Da-Ming Wei

Abstract: The mass ($M_{\rm TOV}$) and radius ($R_{\rm TOV}$) of the maximum-mass nonrotating neutron star (NS) play a crucial role in constraining the elusive equation of state of cold dense matter and in predicting the fate of remnants from binary neutron star (BNS) mergers. In this study, we introduce a novel method to deduce these parameters by examining the mergers of second-generation (2G) black holes… ▽ More The mass ($M_{\rm TOV}$) and radius ($R_{\rm TOV}$) of the maximum-mass nonrotating neutron star (NS) play a crucial role in constraining the elusive equation of state of cold dense matter and in predicting the fate of remnants from binary neutron star (BNS) mergers. In this study, we introduce a novel method to deduce these parameters by examining the mergers of second-generation (2G) black holes (BHs) with NSs. These 2G BHs are assumed to originate from supramassive neutron stars (SMNSs) formed in BNS mergers. Since the properties of the remnant BHs arising from the collapse of SMNSs follow a universal relation governed by $M_{\rm TOV}$ and $R_{\rm TOV}$, we anticipate that by analyzing a series ($\sim 100$ detections) of mass and spin measurements of the 2G BHs using the third-generation ground-based gravitational-wave detectors, $M_{\rm TOV}$ and $R_{\rm TOV}$ can be determined with a precision of $\sim 0.01M_\odot$ and $\sim 0.6$ km, respectively. △ Less

Submitted 22 December, 2023; v1 submitted 27 September, 2023; originally announced September 2023.

Comments: 11 pages, 5 figures, ApJ published

Journal ref: 2024, ApJ, 960, 67

Maximum gravitational mass $M_{\rm TOV}=2.25^{+0.08}_{-0.07}M_\odot$ inferred at about $3\%$ precision with multimessenger data of neutron stars

Authors: Yi-Zhong Fan, Ming-Zhe Han, Jin-Liang Jiang, Dong-Sheng Shao, Shao-Peng Tang

Abstract: The maximal gravitational mass of nonrotating neutron stars ($M_{\rm TOV}$) is one of the key parameters of compact objects and only loose bounds can be set based on the first principle. With reliable measurements of the masses and/or radii of the neutron stars, $M_{\rm TOV}$ can be robustly inferred from either the mass distribution of these objects or the reconstruction of the equation of state… ▽ More The maximal gravitational mass of nonrotating neutron stars ($M_{\rm TOV}$) is one of the key parameters of compact objects and only loose bounds can be set based on the first principle. With reliable measurements of the masses and/or radii of the neutron stars, $M_{\rm TOV}$ can be robustly inferred from either the mass distribution of these objects or the reconstruction of the equation of state (EoS) of the very dense matter. For the first time we take the advantages of both two approaches to have a precise inference of $M_{\rm TOV}=2.25^{+0.08}_{-0.07}~M_\odot$ (68.3\% credibility), with the updated neutron star mass measurement sample, the mass-tidal deformability data of GW170817, the mass-radius data of PSR J0030+0451 and PSR J0740+6620, as well as the theoretical information from the chiral effective theory ($χ$EFT) and perturbative quantum chromodynamics (pQCD) at low and very high energy densities, respectively. This narrow credible range is benefited from the suppression of the high $M_{\rm TOV}$ by the pQCD constraint and the exclusion of the low $M_{\rm TOV}$ by the mass function. Three different EoS reconstruction methods are adopted separately, and the resulting $M_{\rm TOV}$ and $R_{\rm TOV}$ are found to be almost identical, where $R_{\rm TOV}=11.90^{+0.63}_{-0.60}$ km is the radius of the most massive non-rotating NS. This precisely evaluated $M_{\rm TOV}$ suggests that the EoS of neutron star matter is just moderately stiff and the $\sim 2.5-3M_\odot$ compact objects detected by the second generation gravitational wave detectors are most likely the lightest black holes. △ Less

Submitted 29 February, 2024; v1 submitted 22 September, 2023; originally announced September 2023.

Comments: 17 pages, 9 figures. The latest mass measurements of the neutron stars up to April 2023 have been summarized in Table 1 of the Appendix. Published in PRD with moderate revision

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

IllustrisTNG in the HSC-SSP: image data release and the major role of mini mergers as drivers of asymmetry and star formation

Authors: Connor Bottrell, Hassen M. Yesuf, Gergö Popping, Kiyoaki Christopher Omori, Shenli Tang, Xuheng Ding, Annalisa Pillepich, Dylan Nelson, Lukas Eisert, Hua Gao, Andy D. Goulding, Boris S. Kalita, Wentao Luo, Jenny E. Greene, Jingjing Shi, John D. Silverman

Abstract: At fixed galaxy stellar mass, there is a clear observational connection between structural asymmetry and offset from the star forming main sequence, $Δ$SFMS. Herein, we use the TNG50 simulation to investigate the relative roles of major mergers (stellar mass ratios $μ\geq0.25$), minor ($0.1 \leq μ< 0.25$), and mini mergers ($0.01 \leq μ< 0.1$) in driving this connection amongst star forming galaxi… ▽ More At fixed galaxy stellar mass, there is a clear observational connection between structural asymmetry and offset from the star forming main sequence, $Δ$SFMS. Herein, we use the TNG50 simulation to investigate the relative roles of major mergers (stellar mass ratios $μ\geq0.25$), minor ($0.1 \leq μ< 0.25$), and mini mergers ($0.01 \leq μ< 0.1$) in driving this connection amongst star forming galaxies (SFGs). We use dust radiative transfer post-processing with SKIRT to make a large, public collection of synthetic Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) images of simulated TNG galaxies over $0.1\leq z \leq 0.7$ with $\log M_{\star} / \mathrm{M}_{\odot}\geq9$ ($\sim750$k images). Using their instantaneous SFRs, known merger histories/forecasts, and HSC-SSP asymmetries, we show (1) that TNG50 SFGs qualitatively reproduce the observed trend between $Δ$SFMS and asymmetry and (2) a strikingly similar trend emerges between $Δ$SFMS and the time-to-coalescence for mini mergers. Controlling for redshift, stellar mass, environment, and gas fraction, we show that individual mini merger events yield small enhancements in SFRs and asymmetries that are sustained on long timescales (at least $\sim3$ Gyr after coalescence, on average) -- in contrast to major/minor merger remnants which peak at much greater amplitudes but are consistent with controls only $\sim1$ Gyr after coalescence. Integrating the boosts in SFRs and asymmetries driven by $μ\geq0.01$ mergers since $z=0.7$ in TNG50 SFGs, we show that mini mergers are responsible for (i) $55$ per cent of all merger-driven star formation and (ii) $70$ per cent of merger-driven asymmetric structure. Due to their relative frequency and prolonged boost timescales, mini mergers dominate over their minor and major counterparts in driving star formation and asymmetry in SFGs. △ Less

Submitted 7 October, 2023; v1 submitted 28 August, 2023; originally announced August 2023.

Comments: 31 pages; 19 figures; MNRAS accepted; Image data available via the TNG website: www.tng-project.org/bottrell23

Binary Star Evolution in Different Environments: Filamentary, Fractal, Halo and Tidal-tail Clusters

Authors: Xiaoying Pang, Yifan Wang, Shih-Yun Tang, Yicheng Rui, Jing Bai, Chengyuan Li, Fabo Feng, M. B. N. Kouwenhoven, Wen-Ping Chen, Rwei-ju Chuang

Abstract: Using membership of 85 open clusters from previous studies (Pang et al. 2021a,b, 2022b; Li et al. 2021) based on Gaia DR3 data, we identify binary candidates in the color-magnitude diagram, for systems with mass ratio q > 0.4. The binary fraction is corrected for incompleteness at different distances due to the Gaia angular resolution limit. We find a decreasing binary fraction with increasing clu… ▽ More Using membership of 85 open clusters from previous studies (Pang et al. 2021a,b, 2022b; Li et al. 2021) based on Gaia DR3 data, we identify binary candidates in the color-magnitude diagram, for systems with mass ratio q > 0.4. The binary fraction is corrected for incompleteness at different distances due to the Gaia angular resolution limit. We find a decreasing binary fraction with increasing cluster age, with substantial scatter. For clusters with a total mass > 200$M_\odot$, the binary fraction is independent of cluster mass. The binary fraction depends strongly on stellar density. Among four types of cluster environments, the lowest-density filamentary and fractal stellar groups have the highest mean binary fraction: 23.6% and 23.2%, respectively. The mean binary fraction in tidal-tail clusters is 20.8%, and is lowest in the densest halo-type clusters: 14.8%. We find clear evidence of early disruptions of binary stars in the cluster sample. The radial binary fraction depends strongly on the cluster-centric distance across all four types of environments, with the smallest binary fraction within the half-mass radius $r_h$, and increasing towards a few $r_h$. Only hints of mass segregation is found in the target clusters. The observed amount of mass segregation is not significant to generate a global effect inside the target clusters. We evaluate the bias of unresolved binary systems (assuming a primary mass of 1$M_\odot$) in 1D tangential velocity, which is 0.1-1$\,\rm km\,s^{-1}$. Further studies are required to characterize the internal star cluster kinematics using Gaia proper motions. △ Less

Submitted 13 July, 2023; originally announced July 2023.

Black holes as the source of dark energy: a stringent test with high-redshift JWST AGNs

Authors: Lei Lei, Lei Zu, Guan-Wen Yuan, Zhao-Qiang Shen, Yi-Ying Wang, Yuan-Zhu Wang, Zhen-Bo Su, Wen-ke Ren, Shao-Peng Tang, Hao Zhou, Chi Zhang, Zhi-Ping Jin, Lei Feng, Yi-Zhong Fan, Da-Ming Wei

Abstract: Studies have proposed that there is evidence for cosmological coupling of black holes (BHs) with an index of $k\approx 3$; hence, BHs serve as the astrophysical source of dark energy. However, the data sample is limited for the redshifts of $\leq 2.5$. In recent years, the James Webb Space Telescope (JWST) has detected many high-redshift active galactic nuclei (AGNs) and quasars. Among the JWST NI… ▽ More Studies have proposed that there is evidence for cosmological coupling of black holes (BHs) with an index of $k\approx 3$; hence, BHs serve as the astrophysical source of dark energy. However, the data sample is limited for the redshifts of $\leq 2.5$. In recent years, the James Webb Space Telescope (JWST) has detected many high-redshift active galactic nuclei (AGNs) and quasars. Among the JWST NIRSpec-/NIRCam-resolved AGNs, three are determined to be in early-type host galaxies with a redshift of $z\sim 4.5--7$. However, their $M_{\star}$ and $M_{\rm BH}$ are in tension with the predicted cosmological coupling of black holes with $k = 3$ at a confidence level of $\sim 2σ$, which challenges the hypothesis that BHs serve as the origin of dark energy. Future work on high-redshift AGNs using the JWST will further assess such a hypothesis by identifying more early-type host galaxies in the higher mass range. △ Less

Submitted 17 January, 2024; v1 submitted 5 May, 2023; originally announced May 2023.

Comments: 9 pages, 3 figures, 2 tables; Comments are welcome!

Journal ref: Sci. China-Phys. Mech. Astron. 67, 229811 (2024)

Morphological asymmetries of quasar host galaxies with Subaru Hyper Suprime-Cam

Authors: Shenli Tang, John D. Silverman, Hassen M. Yesuf, Xuheng Ding, Junyao Li, Connor Bottrell, Andy Goulding, Kiyoaki Christopher Omori, Yoshiki Toba, Toshihiro Kawaguchi

Abstract: How does the host galaxy morphology influence a central quasar or vice versa? We address this question by measuring the asymmetries of 2424 SDSS quasar hosts at $0.2<z<0.8$ using broad-band ($grizy$) images from the Hyper Suprime-Cam Subaru Strategic Program. Control galaxies (without quasars) are selected by matching the redshifts and stellar masses of the quasar hosts. A two-step pipeline is run… ▽ More How does the host galaxy morphology influence a central quasar or vice versa? We address this question by measuring the asymmetries of 2424 SDSS quasar hosts at $0.2<z<0.8$ using broad-band ($grizy$) images from the Hyper Suprime-Cam Subaru Strategic Program. Control galaxies (without quasars) are selected by matching the redshifts and stellar masses of the quasar hosts. A two-step pipeline is run to decompose the PSF and \sersic\ components, and then measure asymmetry indices ($A_{\rm CAS}$, $A_{\rm outer}$, and $A_{\rm shape}$) of each quasar host and control galaxy. We find a mild correlation between host asymmetry and AGN bolometric luminosity ($L_{\rm bol}$) for the full sample (spearman correlation of 0.37) while a stronger trend is evident at the highest luminosities ($L_{\rm bol}>45$). This then manifests itself into quasar hosts being more asymmetric, on average, when they harbor a more massive and highly accreting black hole. The merger fraction also positively correlates with $L_{\rm bol}$ and reaches up to 35\% for the most luminous. Compared to control galaxies, quasar hosts are marginally more asymmetric (excess of 0.017 in median at 9.4$σ$ level) and the merger fractions are similar ($\sim 16.5\%$). We quantify the dependence of asymmetry on optical band which demonstrates that mergers are more likely to be identified with the bluer bands and the correlation between $L_{\rm bol}$ and asymmetry is also stronger in such bands. We stress that the band dependence, indicative of a changing stellar population, is an important factor in considering the influence of mergers on AGN activity. △ Less

Submitted 12 April, 2023; originally announced April 2023.

Comments: 27 pages, 28 figures

Journal ref: MNRAS, volume 521, number 4, pages 5272--5297 (2023)

Star-Crossed Lovers DI Tau A and B: Orbit Characterization and Physical Properties Determination

Authors: Shih-Yun Tang, Asa G. Stahl, L. Prato, G. H. Schaefer, Christopher M. Johns-Krull, Brian A. Skiff, Charles A. Beichman, Taichi Uyama

Abstract: The stellar companion to the weak-line T Tauri star DI Tau A was first discovered by the lunar occultation technique in 1989 and was subsequently confirmed by a speckle imaging observation in 1991. It has not been detected since, despite being targeted by five different studies that used a variety of methods and spanned more than 20 years. Here, we report the serendipitous rediscovery of DI Tau B… ▽ More The stellar companion to the weak-line T Tauri star DI Tau A was first discovered by the lunar occultation technique in 1989 and was subsequently confirmed by a speckle imaging observation in 1991. It has not been detected since, despite being targeted by five different studies that used a variety of methods and spanned more than 20 years. Here, we report the serendipitous rediscovery of DI Tau B during our Young Exoplanets Spectroscopic Survey (YESS). Using radial velocity data from YESS spanning 17 years, new adaptive optics observations from Keck II, and a variety of other data from the literature, we derive a preliminary orbital solution for the system that effectively explains the detection and (almost all of the) non-detection history of DI Tau B. We estimate the dynamical masses of both components, finding that the large mass difference (q $\sim$0.17) and long orbital period ($\gtrsim$35 years) make DI Tau system a noteworthy and valuable addition to studies of stellar evolution and pre-main-sequence models. With a long orbital period and a small flux ratio (f2/f1) between DI Tau A and B, additional measurements are needed for a better comparison between these observational results and pre-main-sequence models. Finally, we report an average surface magnetic field strength ($\bar B$) for DI Tau A, of $\sim$0.55 kG, which is unusually low in the context of young active stars. △ Less

Submitted 20 March, 2023; originally announced March 2023.

Comments: 21 pages, 10 figures. Accepted to ApJ

Resolving the stellar-collapse and hierarchical-merger origins of the coalescing black holes

Authors: Yin-Jie Li, Yuan-Zhu Wang, Shao-Peng Tang, Yi-Zhong Fan

Abstract: Spin and mass properties provide essential clues in distinguishing the origins of coalescing black holes (BHs). With a dedicated semi-parametric population model for the coalescing binary black holes (BBHs), we identify two distinct categories of BHs among the GWTC-3 events, which is {favored over the one population scenario by} a logarithmic Bayes factor ($\ln\mathcal{B}$) of 7.5. One category, w… ▽ More Spin and mass properties provide essential clues in distinguishing the origins of coalescing black holes (BHs). With a dedicated semi-parametric population model for the coalescing binary black holes (BBHs), we identify two distinct categories of BHs among the GWTC-3 events, which is {favored over the one population scenario by} a logarithmic Bayes factor ($\ln\mathcal{B}$) of 7.5. One category, with a mass ranging from $\sim 25M_\odot$ to $\sim 80M_\odot$, is distinguished by the high spin-magnitudes ($\sim0.75$) and consistent with the hierarchical merger origin. The other category, characterized by low spins, has a sharp mass cutoff at $\sim 40M_\odot$, which is natural for the stellar-collapse origin and in particular the pair-instability explosion of massive stars. We infer the local hierarchical merger rate density as $0.46^{+0.61}_{-0.24}~{\rm Gpc^{-3}yr^{-1}}$. Additionally, we find that a fraction of the BBHs has a cosine-spin-tilt-angle distribution concentrated preferentially around $1$, which is consistent with the expectation of the isolated field binaries and the dynamically formed BBHs in the accretion disks of Active Galactic Nuclei. △ Less

Submitted 13 April, 2024; v1 submitted 6 March, 2023; originally announced March 2023.

Comments: Main text: 6 pages and 3 figures; Supplementary: 26 pages, 34 figures, and 4 tables

A Multilevel Scheduling Framework for Distributed Time-domain Large-area Sky Survey Telescope Array

Authors: Yajie Zhang, Ce Yu, Chao Sun, Zhaohui Shang, Yi Hu, Huiyu Zhi, Jinmao Yang, Shanjiang Tang

Abstract: Telescope arrays are receiving increasing attention due to their promise of higher resource utilization, greater sky survey area, and higher frequency of full space-time monitoring than single telescopes. Compared with the ordinary coordinated operation of several telescopes, the new astronomical observation mode has an order of magnitude difference in the number of telescopes. It requires efficie… ▽ More Telescope arrays are receiving increasing attention due to their promise of higher resource utilization, greater sky survey area, and higher frequency of full space-time monitoring than single telescopes. Compared with the ordinary coordinated operation of several telescopes, the new astronomical observation mode has an order of magnitude difference in the number of telescopes. It requires efficient coordinated observation by large-domain telescopes distributed at different sites. Coherent modeling of various abstract environmental constraints is essential for responding to multiple complex science goals. Also, due to competing science priorities and field visibility, how the telescope arrays are scheduled for observations can significantly affect observation efficiency. This paper proposes a multilevel scheduling model oriented toward the problem of telescope-array scheduling for time-domain surveys. A flexible framework is developed with basic functionality encapsulated in software components implemented on hierarchical architectures. An optimization metric is proposed to self-consistently weight contributions from time-varying observation conditions to maintain uniform coverage and efficient time utilization from a global perspective. The performance of the scheduler is evaluated through simulated instances. The experimental results show that our scheduling framework performs correctly and provides acceptable solutions considering the percentage of time allocation efficiency and sky coverage uniformity in a feasible amount of time. Using a generic version of the telescope-array scheduling framework, we also demonstrate its scalability and its potential to be applied to other astronomical applications. △ Less

Submitted 3 February, 2023; v1 submitted 18 January, 2023; originally announced January 2023.

Comments: Accepted for publication in the Astronomical Journal

Hierarchical triple mergers: testing Hawking's area theorem with the inspiral signals

Authors: Shao-Peng Tang, Yi-Zhong Fan, Da-Ming Wei

Abstract: Hawking's area theorem is one of the fundamental laws of black holes (BHs), which has been tested at a confidence level of $\sim 95\%$ with gravitational wave (GW) observations by analysing the inspiral and ringdown portions of GW signals, independently. In this work, we propose to carry out the test in a new way with the hierarchical triple merger (i.e., two successive BH mergers occurred sequent… ▽ More Hawking's area theorem is one of the fundamental laws of black holes (BHs), which has been tested at a confidence level of $\sim 95\%$ with gravitational wave (GW) observations by analysing the inspiral and ringdown portions of GW signals, independently. In this work, we propose to carry out the test in a new way with the hierarchical triple merger (i.e., two successive BH mergers occurred sequentially within the observation window of GW detectors), for which the properties of the progenitor BHs and the remnant BH of the first coalescence can be inferred from the inspiral portions of the two mergers. As revealed in our simulations, the BH area law can be well confirmed for some plausible hierarchical triple merger events detectable in LIGO/Virgo/KAGRA/LIGO-India's O4/O5 runs. Our proposed method provides significant facilitation for testing the area law and complements previous methods. △ Less

Submitted 19 June, 2023; v1 submitted 8 September, 2022; originally announced September 2022.

Comments: 8 pages, 8 figures, 2 tables, published in MNRAS

Journal ref: 2023, MNRAS, 523, 4113

Dynamical Origin for the Collinder 132-Gulliver 21 Stream: A Mixture of three Co-Moving Populations with an Age Difference of 250 Myr

Authors: Xiaoying Pang, Yuqian Li, Shih-Yun Tang, Long Wang, Yanshu Wang, Zhaoyu Li, Danchen Wang, M. B. N. Kouwenhoven, Mario Pasquato

Abstract: We use Gaia DR3 data to study the Collinder 132-Gulliver 21 region via the machine learning algorithm StarGO, and find eight subgroups of stars (ASCC 32, Collinder 132 gp 1--6, Gulliver 21) located in close proximity. Three co-moving populations were identified among these eight subgroups: (i) a coeval 25 Myr-old moving group (Collinder 132); (ii) an intermediate-age (50--100 Myr) group; and (iii)… ▽ More We use Gaia DR3 data to study the Collinder 132-Gulliver 21 region via the machine learning algorithm StarGO, and find eight subgroups of stars (ASCC 32, Collinder 132 gp 1--6, Gulliver 21) located in close proximity. Three co-moving populations were identified among these eight subgroups: (i) a coeval 25 Myr-old moving group (Collinder 132); (ii) an intermediate-age (50--100 Myr) group; and (iii) the 275 Myr-old dissolving cluster Gulliver 21. These three populations form parallel diagonal stripe-shape over-densities in the U--V distribution, which differ from open clusters and stellar groups in the solar neighborhood. We name this kinematic structure the Collinder 132-Gulliver 21 stream, as it extends over 270 pc in the 3D space. The oldest population Gulliver21 is spatially surrounded by the Collinder 132 moving group and the intermediate-age group. Stars in the Collinder 132-Gulliver 21 stream have an age difference up to 250 Myr. Metallicity information shows a variation of 0.3 dex between the youngest and oldest populations. The formation of the Collinder132-Gulliver 21 stream involves both star formation and dynamical heating. The youngest population (Collinder 132 moving group) with homogeneous metallicity is probably formed through filamentary star formation. The intermediate-age and the oldest population were then scatted by the Galactic bar or spiral structure resonance to intercept Collinder 132's orbit. Without mutual interaction between each population, the three populations are flying by each other currently and will become distinct three groups again in approximately ~50Myr. △ Less

Submitted 30 August, 2022; originally announced August 2022.

Comments: Accepted to ApJL, 15 pages, 7 figures

Deciphering the extreme X-ray variability of the nuclear transient eRASSt J045650.3-203750: A likely repeating partial tidal disruption event

Authors: Zhu Liu, A. Malyali, M. Krumpe, D. Homan, A. J. Goodwin, I. Grotova, A. Kawka, A. Rau, A. Merloni, G. E. Anderson, J. C. A. Miller-Jones, A. G. Markowitz, S. Ciroi, F. Di Mille, M. Schramm, Shenli Tang, D. A. H. Buckley, M. Gromadzki, Chichuan Jin, J. Buchner

Abstract: (Abridged) In this paper, we present the results of an exceptional repeating X-ray nuclear transient, eRASSt J045650.3-203750 (hereafter J0456-20), uncovered by SRG/eROSITA in a quiescent galaxy at redshift of z~0.077. The main results are: 1) J0456-20 cycles through four distinctive phases: an X-ray rising phase leading into an X-ray plateau phase which lasts for ~2 months. This is terminated by… ▽ More (Abridged) In this paper, we present the results of an exceptional repeating X-ray nuclear transient, eRASSt J045650.3-203750 (hereafter J0456-20), uncovered by SRG/eROSITA in a quiescent galaxy at redshift of z~0.077. The main results are: 1) J0456-20 cycles through four distinctive phases: an X-ray rising phase leading into an X-ray plateau phase which lasts for ~2 months. This is terminated by a rapid X-ray flux drop phase during which the X-ray flux can drastically drop by more than a factor of 100 within 1 week followed by an X-ray faint state for about two months before it starts the X-ray rising phase again; 2) the X-ray spectra are generally soft in the rising phase with a photon index >3.0, and become harder as the X-ray flux increases. There is evidence of a multi-colour disk with inner region temperature of $T_\text{in}=70$ eV at the beginning of the X-ray rising phase. The high quality XMM-Newton data suggest that a warm and hot corona could be responsible for the X-ray emission, through inverse Comptonisation of soft disk seed photons, during the plateau phase and at the bright end of the rising phase; 3) J0456-20 shows only moderate UV variability and no significant optical variability; 4) radio emission is only detected (as yet) in the X-ray plateau phase, and shows a rapid decline on a time-scale of 2 weeks. We conclude that J0456-20 is likely a repeating nuclear transient with a tentative recurrence time of ~223 days. We discuss several possibilities to explain J0456-20's observational properties, and currently favour a repeating partial tidal disruption event (TDE) as the most likely scenario. The long-term X-ray evolution is explained as a transition between a thermal disk-dominated soft state and a steep power-law state, implying that the corona can be formed within a few months and destroyed within a few weeks. △ Less

Submitted 13 January, 2023; v1 submitted 26 August, 2022; originally announced August 2022.

Comments: 27 pages, 13 figures, published in A&A

Journal ref: A&A 669, A75 (2023)

Potential Subpopulations and Assembling Tendency of the Merging Black Holes

Authors: Yuan-Zhu Wang, Yin-Jie Li, Jorick S. Vink, Yi-Zhong Fan, Shao-Peng Tang, Ying Qin, Da-Ming Wei

Abstract: The origins of coalescing binary black holes (BBHs) detected by the advanced LIGO/Virgo are still under debate, and clues may be present in the joint mass-spin distribution of these merger events. Here we construct phenomenological models containing two sub-populations to investigate the BBH population detected in gravitational wave observations. We find that our models can explain the GWTC-3 data… ▽ More The origins of coalescing binary black holes (BBHs) detected by the advanced LIGO/Virgo are still under debate, and clues may be present in the joint mass-spin distribution of these merger events. Here we construct phenomenological models containing two sub-populations to investigate the BBH population detected in gravitational wave observations. We find that our models can explain the GWTC-3 data rather well, and several constraints to our model are required by the data: first, the maximum mass for the component with a stellar-origin, $m_{\rm max}$, is $39.1^{+2.4}_{-2.7}M_{\odot}$ at 90\% credibility; second, about $15\%$ of the mergers happen in dynamical environments, in which $7-16\%$ of events are hierarchical mergers, and these BHs have an average spin magnitude significantly larger than the first-generation mergers, with ${\rm d}μ_{\rm a} > 0.4 $ at $99\%$ credibility; third, the dynamical component BHs tend to pair with each other with larger total mass and higher mass ratio. An independent analysis focusing on spins is also carried out, and we find that the spin amplitude of component BHs can be divided into two groups according to a division mass $m_{\rm d} = 46.1^{+5.6}_{-5.1}M_{\odot}$. These constraints can be naturally explained by current formation channels, and our results suggest that some of the observed events were likely from AGN disks. △ Less

Submitted 12 January, 2023; v1 submitted 25 August, 2022; originally announced August 2022.

Comments: 21 pages, 8 figures. Published in The Astrophysical Journal Letters

The late afterglow of GW170817/GRB170817A: a large viewing angle and the shift of the Hubble constant to a value more consistent with the local measurements

Authors: Yi-Ying Wang, Shao-Peng Tang, Zhi-Ping Jin, Yi-Zhong Fan

Abstract: The multi-messenger data of neutron star merger events are promising for constraining the Hubble constant. So far, GW170817 is still the unique gravitational wave event with multi-wavelength electromagnetic counterparts. In particular, its radio and X-ray emission have been measured in the past $\sim 3-5$ years. In this work, we fit the long-lasting X-ray, optical, and radio afterglow light curves… ▽ More The multi-messenger data of neutron star merger events are promising for constraining the Hubble constant. So far, GW170817 is still the unique gravitational wave event with multi-wavelength electromagnetic counterparts. In particular, its radio and X-ray emission have been measured in the past $\sim 3-5$ years. In this work, we fit the long-lasting X-ray, optical, and radio afterglow light curves of GW170817/GRB 170817A, including the forward shock radiation from both the decelerating relativistic GRB outflow and the sub-relativistic kilonova outflow (though whether the second component contributes significantly is still uncertain), and find out a relatively large viewing angle ($\sim 0.5\, \rm rad$). Such a viewing angle has been taken as a prior in the gravitational wave data analysis, and the degeneracy between the viewing angle and the luminosity distance is broken. Finally, we have a Hubble constant $H_0=72.57^{+4.09}_{-4.17}\, \rm km\, s^{-1}\, Mpc^{-1}$, which is more consistent with that obtained by other local measurements. If rather similar values are inferred from multi-messenger data of future neutron star merger events, it will provide critical support to the existence of the Hubble tension. △ Less

Submitted 5 December, 2022; v1 submitted 18 August, 2022; originally announced August 2022.

Comments: 12 pages, 5 figures, and 2 tables. Accepted for publication in ApJ

Plausible presence of new state in neutron stars with masses above $0.98M_{\rm TOV}$

Authors: Ming-Zhe Han, Yong-Jia Huang, Shao-Peng Tang, Yi-Zhong Fan

Abstract: We investigate the neutron star (NS) equation of state (EOS) by incorporating multi-messenger data of GW170817, PSR J0030+0451, PSR J0740+6620, and state-of-the-art theoretical progresses, including the information from chiral effective field theory ($χ$EFT) and perturbative quantum chromodynamics (pQCD) calculation. Taking advantage of the various structures sampling by a single-layer feed-forwar… ▽ More We investigate the neutron star (NS) equation of state (EOS) by incorporating multi-messenger data of GW170817, PSR J0030+0451, PSR J0740+6620, and state-of-the-art theoretical progresses, including the information from chiral effective field theory ($χ$EFT) and perturbative quantum chromodynamics (pQCD) calculation. Taking advantage of the various structures sampling by a single-layer feed-forward neural network model embedded in the Bayesian nonparametric inference, the structure of NS matter's sound speed $c_{\rm s}$ is explored in a model-agnostic way. It is found that a peak structure is common in the $c_{\rm s}^2$ posterior, locating at $2.4-4.8ρ_{\rm sat}$ (nuclear saturation density) and $c_{\rm s}^2$ exceeds ${c^{2}}/{3}$ at 90\% credibility. The non-monotonic behavior suggests evidence of the state deviating from hadronic matter inside the very massive NSs. Assuming the new/exotic state is featured as it is softer than typical hadronic models or even with hyperons, we find that a sizable ($\geq 10^{-3}M_\odot$) exotic core, likely made of quark matter, is plausible for the NS with a gravitational mass above about $0.98M_{\rm TOV}$, where $M_{\rm TOV}$ represents the maximum gravitational mass of a non-rotating cold NS. The inferred $M_{\rm TOV} = 2.18^{+0.27}_{-0.13}M_\odot$ (90\% credibility) is well consistent with the value of $2.17^{+0.15}_{-0.12}M_\odot$ estimated independently with GW170817/GRB 170817A/AT2017gfo assuming a temporary supramassive NS remnant formed after the merger. PSR J0740+6620, the most massive NS detected so far, may host an exotic core with a probability of $\approx 0.36$. △ Less

Submitted 22 May, 2023; v1 submitted 27 July, 2022; originally announced July 2022.

Comments: 8+5 pages, 6+5 figures, published in Science Bulletin

Report number: RIKEN-iTHEMS-Report-23

Journal ref: Science Bulletin, Volume 68, Issue 9, 15 May 2023, Pages 913-919

Reduced Order Model for Chemical Kinetics: A case study with Primordial Chemical Network

Authors: Kwok Sun Tang, Matthew Turk

Abstract: Chemical kinetics plays an important role in governing the thermal evolution in reactive flows problems. The possible interactions between chemical species increase drastically with the number of species considered in the system. Various ways have been proposed before to simplify chemical networks with an aim to reduce the computational complexity of the chemical network. These techniques oftentim… ▽ More Chemical kinetics plays an important role in governing the thermal evolution in reactive flows problems. The possible interactions between chemical species increase drastically with the number of species considered in the system. Various ways have been proposed before to simplify chemical networks with an aim to reduce the computational complexity of the chemical network. These techniques oftentimes require domain-knowledge experts to handcraftedly identify important reaction pathways and possible simplifications. Here, we propose a combination of autoencoder and neural ordinary differential equation to model the temporal evolution of chemical kinetics in a reduced subspace. We demonstrated that our model has achieved a close-to 10-fold speed-up compared to commonly used astro-chemistry solver for a 9-species primordial network, while maintaining 1 percent accuracy across a wide-range of density and temperature. △ Less

Submitted 14 July, 2022; originally announced July 2022.

Comments: 10 pages, 8 figures, accepted to the ICML 2022 Machine Learning for Astrophysics workshop

First report of a solar energetic particle event observed by China's Tianwen-1 mission in transit to Mars

Authors: Shuai Fu, Zheyi Ding, Yongjie Zhang, Xiaoping Zhang, Cunhui Li, Gang Li, Shuwen Tang, Haiyan Zhang, Yi Xu, Yuming Wang, Jingnan Guo, Lingling Zhao, Yi Wang, Xiangyu Hu, Pengwei Luo, Zhiyu Sun, Yuhong Yu, Lianghai Xie

Abstract: Solar energetic particles (SEPs) associated with flares and/or coronal mass ejection (CME)-driven shocks can impose acute radiation hazards to space explorations. To measure energetic particles in near-Mars space, the Mars Energetic Particle Analyzer (MEPA) instrument onboard China's Tianwen-1 (TW-1) mission was designed. Here, we report the first MEPA measurements of the widespread SEP event occu… ▽ More Solar energetic particles (SEPs) associated with flares and/or coronal mass ejection (CME)-driven shocks can impose acute radiation hazards to space explorations. To measure energetic particles in near-Mars space, the Mars Energetic Particle Analyzer (MEPA) instrument onboard China's Tianwen-1 (TW-1) mission was designed. Here, we report the first MEPA measurements of the widespread SEP event occurring on 29 November 2020 when TW-1 was in transit to Mars. This event occurred when TW-1 and Earth were magnetically well connected, known as the Hohmann-Parker effect, thus offering a rare opportunity to understand the underlying particle acceleration and transport process. Measurements from TW-1 and near-Earth spacecraft show similar double-power-law spectra and a radial dependence of the SEP peak intensities. Moreover, the decay phases of the time-intensity profiles at different locations clearly show the reservoir effect. We conclude that the double-power-law spectrum is likely generated at the acceleration site, and that a small but finite cross-field diffusion is crucial to understand the formation of the SEP reservoir phenomenon. These results provide insight into particle acceleration and transport associated with CME-driven shocks, which may contribute to the improvement of relevant physical models. △ Less

Submitted 14 July, 2022; originally announced July 2022.

Comments: This manuscript has been accepted for publication in The Astrophysical Journal Letters

Galaxy Merger Reconstruction with Equivariant Graph Normalizing Flows

Authors: Kwok Sun Tang, Yuan-Sen Ting

Abstract: A key yet unresolved question in modern-day astronomy is how galaxies formed and evolved under the paradigm of the $Λ$CDM model. A critical limiting factor lies in the lack of robust tools to describe the merger history through a statistical model. In this work, we employ a generative graph network, E(n) Equivariant Graph Normalizing Flows Model. We demonstrate that, by treating the progenitors as… ▽ More A key yet unresolved question in modern-day astronomy is how galaxies formed and evolved under the paradigm of the $Λ$CDM model. A critical limiting factor lies in the lack of robust tools to describe the merger history through a statistical model. In this work, we employ a generative graph network, E(n) Equivariant Graph Normalizing Flows Model. We demonstrate that, by treating the progenitors as a graph, our model robustly recovers their distributions, including their masses, merging redshifts and pairwise distances at redshift z=2 conditioned on their z=0 properties. The generative nature of the model enables other downstream tasks, including likelihood-free inference, detecting anomalies and identifying subtle correlations of progenitor features. △ Less

Submitted 6 July, 2022; originally announced July 2022.

Comments: 6 pages, 3 figures, accepted to the ICML 2022 Machine Learning for Astrophysics workshop

Nonparametric Representation of Neutron Star Equation of State Using Variational Autoencoder

Authors: Ming-Zhe Han, Shao-Peng Tang, Yi-Zhong Fan

Abstract: We introduce a new nonparametric representation of the neutron star (NS) equation of state (EoS) by using the variational autoencoder (VAE). As a deep neural network, the VAE is frequently used for dimensionality reduction since it can compress input data to a low-dimensional latent space using the encoder component and then reconstruct the data using the decoder component. Once a VAE is trained,… ▽ More We introduce a new nonparametric representation of the neutron star (NS) equation of state (EoS) by using the variational autoencoder (VAE). As a deep neural network, the VAE is frequently used for dimensionality reduction since it can compress input data to a low-dimensional latent space using the encoder component and then reconstruct the data using the decoder component. Once a VAE is trained, one can take the decoder of the VAE as a generator. We employ 100,000 EoSs that are generated using the nonparametric representation method based on \citet{2021ApJ...919...11H} as the training set and try different settings of the neural network, then we get an EoS generator (trained VAE's decoder) with four parameters. We use the mass\textendash{}tidal-deformability data of binary neutron star (BNS) merger event GW170817, the mass\textendash{}radius data of PSR J0030+0451, PSR J0740+6620, PSR J0437-4715, and 4U 1702-429, and the nuclear constraints to perform the joint Bayesian inference. The overall results of the analysis that includes all the observations are $R_{1.4}=12.59^{+0.36}_{-0.42}\,\rm km$, $Λ_{1.4}=489^{+114}_{-110}$, and $M_{\rm max}=2.20^{+0.37}_{-0.19}\,\rm M_\odot$ ($90\%$ credible levels), where $R_{1.4}$/$Λ_{1.4}$ are the radius/tidal-deformability of a canonical $1.4\,\rm M_\odot$ NS, and $M_{\rm max}$ is the maximum mass of a non-rotating NS. The results indicate that the implementation of the VAE techniques can obtain the reasonable results, while accelerate calculation by a factor of $\sim$ 3\textendash10 or more, compared with the original method. △ Less

Submitted 14 June, 2023; v1 submitted 8 May, 2022; originally announced May 2022.

Comments: 10 pages, 4 figures, 1 table, published in ApJ

Journal ref: The Astrophysical Journal, Volume 950, Number 2, 77 (2023)

3D Morphology of Open Clusters in the Solar Neighborhood with Gaia EDR3 II: Hierarchical Star Formation Revealed by Spatial and Kinematic Substructures

Authors: Xiaoying Pang, Shih-Yun Tang, Yuqian Li, Zeqiu Yu, Long Wang, Jiayu Li, Yezhang Li, Yifan Wang, Yanshu Wang, Teng Zhang, Mario Pasquato, M. B. N. Kouwenhoven

Abstract: We identify members of 65 open clusters in the solar neighborhood using the machine-learning algorithm StarGO based on Gaia EDR3 data. After adding members of twenty clusters from previous studies (Pang et al. 2021a,b; Li et al. 2021) we obtain 85 clusters, and study their morphology and kinematics. We classify the substructures outside the tidal radius into four categories: filamentary (f1) and f… ▽ More We identify members of 65 open clusters in the solar neighborhood using the machine-learning algorithm StarGO based on Gaia EDR3 data. After adding members of twenty clusters from previous studies (Pang et al. 2021a,b; Li et al. 2021) we obtain 85 clusters, and study their morphology and kinematics. We classify the substructures outside the tidal radius into four categories: filamentary (f1) and fractal (f2) for clusters $<100$ Myr, and halo (h) and tidal-tail (t) for clusters $>100$ Myr. The kinematical substructures of f1-type clusters are elongated; these resemble the disrupted cluster Group X. Kinematic tails are distinct in t-type clusters, especially Pleiades. We identify 29 hierarchical groups in four young regions (Alessi 20, IC 348, LP 2373, LP 2442); ten among these are new. The hierarchical groups form filament networks. Two regions (Alessi 20, LP 2373) exhibit global "orthogonal" expansion (stellar motion perpendicular to the filament), which might cause complete dispersal. Infalling-like flows (stellar motion along the filament) are found in UBC 31 and related hierarchical groups in the IC 348 region. Stellar groups in the LP 2442 region (LP 2442 gp 1-5) are spatially well-mixed but kinematically coherent. A merging process might be ongoing in the LP 2442 subgroups. For younger systems ($\lesssim30$ Myr), the mean axis ratio, cluster mass and half-mass radius tend to increase with age values. These correlations between structural parameters may imply two dynamical processes occurring in the hierarchical formation scenario in young stellar groups: (1) filament dissolution and (2) sub-group mergers. △ Less

Submitted 12 April, 2022; originally announced April 2022.

Comments: Accepted to ApJ, 33 pages, 14 figures. OC 3D interactive visualization included on http://3doc-morphology.lowell.edu

Divergence in mass ratio distributions between low-mass and high-mass coalescing binary black holes

Authors: Yin-Jie Li, Yuan-Zhu Wang, Shao-Peng Tang, Qiang Yuan, Yi-Zhong Fan, Da-Ming Wei

Abstract: Coalescing binary black hole (BBH) systems are likely formed via several channels, and it is challenging to understand their formation / evolutionary processes. Some features in the mass function of the primary components ($m_1$), such as the distinct Gaussian-like peak located at $\sim 34M_\odot$, have been previously found. In this work, we investigate the possible dependence of the mass ratio (… ▽ More Coalescing binary black hole (BBH) systems are likely formed via several channels, and it is challenging to understand their formation / evolutionary processes. Some features in the mass function of the primary components ($m_1$), such as the distinct Gaussian-like peak located at $\sim 34M_\odot$, have been previously found. In this work, we investigate the possible dependence of the mass ratio ($q=m_2/m_1$) distribution on the primary mass. We find a Bayesian odds ratio of 18.1 in favor of divergence in the mass ratio distributions between the low- and high-mass ranges over an invariable mass ratio distribution. The BBHs with $m_1\gtrsim29M_{\odot}$ have a stronger preference to be symmetric compared to those with $m_1\lesssim29M_{\odot}$ at a 97.6\% credible level. Additionally, we find mild evidence that the BBHs with $m_1$ located in the Gaussian-like peak have a mass ratio distribution different from that of other BBHs. Our findings may be in favor of some formation channels, such as the chemically homogeneous evolution and the dynamical assembly in globular clusters/nuclear star clusters, which are more likely to provide symmetric BBHs in the high-mass range. △ Less

Submitted 1 July, 2022; v1 submitted 5 January, 2022; originally announced January 2022.

Comments: 11 pages, 6 figures, 2 tables, published In ApJL

Journal ref: The Astrophysical Journal Letters, 933:L14 (9pp), 2022 July 1

Total reflection of a flare-driven quasi-periodic EUV wave train at a coronal hole boundary

Authors: Xinping Zhou, Yuandeng Shen, Zehao Tang, Chengrui Zhou, Yadan Duan, Song Tang

Abstract: The reflection, refraction, and transmission of large-scale extreme ultraviolet (EUV) waves (collectively, secondary waves) have been observed during their interactions with coronal structures such as active regions (ARs) and coronal holes (CHs). However, the effect of the total reflection of EUV waves has not been reported in the literature. Here, we present the first unambiguous observational ev… ▽ More The reflection, refraction, and transmission of large-scale extreme ultraviolet (EUV) waves (collectively, secondary waves) have been observed during their interactions with coronal structures such as active regions (ARs) and coronal holes (CHs). However, the effect of the total reflection of EUV waves has not been reported in the literature. Here, we present the first unambiguous observational evidence of the total reflection of a quasi-periodic EUV wave train during its interaction with a polar CH. The event occurred in NOAA AR 12473, located close to the southeast limb of the solar disk, and was characterized by a jet-like CME. In this study, we focus in particular on the driving mechanism s of the quasi-periodic wave train and the total reflection effect at the CH boundary. We find that the periods of the incident and the reflected wave trains are both about 100 seconds. The excitation of the quasi-periodic wave train was possibly due to the intermittent energy release in the associated flare since its period is similar to that of the quasi-periodic pulsations in the associated flare. Our observational results showed that the reflection of the wave train at the boundary of the CH was a total reflection because the measured incidence and critical angles satisfy the theory of total reflection, i.e., the incidence angle is less than the critical angle. △ Less

Submitted 30 December, 2021; originally announced December 2021.

Comments: 8 pages, 5 figures, accepted for publication in A&A

Journal ref: A&A 659, A164 (2022)

Prospects of calibrating afterglow modeling of short GRBs with gravitational wave inclination angle measurements and resolving the Hubble tension with a GW-GRB association event

Authors: Yi-Ying Wang, Shao-Peng Tang, Xin-Yu Li, Zhi-Ping Jin, Yi-Zhong Fan

Abstract: In the numerical modeling of the GRB afterglow, some approximations have been made for simplicity, and different groups developed their codes. A robust test of these models/approaches is challenging because of the lack of directly measured physical parameters. Fortunately, the viewing angle inferred from the afterglow modeling is widely anticipated to be the same as the inclination angle of the bi… ▽ More In the numerical modeling of the GRB afterglow, some approximations have been made for simplicity, and different groups developed their codes. A robust test of these models/approaches is challenging because of the lack of directly measured physical parameters. Fortunately, the viewing angle inferred from the afterglow modeling is widely anticipated to be the same as the inclination angle of the binary neutron star (BNS) mergers that can be evaluated with the gravitational wave (GW) data. Therefore in the future, it is possible to calibrate the afterglow modeling with the GW inclination angle measurements. We take three methods, including both analytical estimations and direct simulations, to project the uncertainties of the inclination angle measurements. For some BNS mergers accompanied with electromagnetic counterparts detected in the O4/O5 runs of LIGO/Virgo/KAGRA/LIGO-India detectors, we show that the inclination angle can be determined within an uncertainty of $\leq 0.1$ rad, supposing that the Hubble constant is known with an accuracy of $\leq 3\%$ and the uncertainty of Hubble flow velocity is within $\sim 1 \%$. The off-axis GRB outflow will give rise to afterglow emission, and the most energetic ones may be detectable at the distance of $\sim$ 100-200 Mpc even for a viewing angle of $\geq 0.3$ rad. Such events can thus serve as a robust test of the afterglow modeling approach. We have also evaluated the prospect of resolving the so-called Hubble tension with a single GW/GRB association event. We find out that a $\sim 3\%$ precision Hubble constant is obtainable if the uncertainty of the viewing angle can be constrained to be within $\sim 0.1$ rad, which is expected to be the case for some nearby ($\leq 250$ Mpc) bright/on-axis GRBs with a well-behaved afterglow light curve displaying a clear achromatic break at early times. △ Less

Submitted 26 June, 2022; v1 submitted 3 November, 2021; originally announced November 2021.

Comments: 18 pages, 11 figures, and 1 table; PRD in press

Mass and effective spin distributions of binary black hole events: Evidence for a population of hierarchical black hole mergers in active galactic nuclei disks

Authors: Yuan-Zhu Wang, Yi-Zhong Fan, Shao-Peng Tang, Ying Qin, Da-Ming Wei

Abstract: The origins of the coalescing binary black holes (BBHs) detected by the advanced LIGO/Virgo are still in debate and clues may present in the mass and effective spin ($χ_{\rm eff}$) distributions of these merger events. Here we analyze the GWTC-1 and GWTC-2.1 BBH events, and find evidence for a primary mass dependent $χ_{\rm eff}$ distribution. Below (above) the mass of $\sim 50M_\odot$, a nearly z… ▽ More The origins of the coalescing binary black holes (BBHs) detected by the advanced LIGO/Virgo are still in debate and clues may present in the mass and effective spin ($χ_{\rm eff}$) distributions of these merger events. Here we analyze the GWTC-1 and GWTC-2.1 BBH events, and find evidence for a primary mass dependent $χ_{\rm eff}$ distribution. Below (above) the mass of $\sim 50M_\odot$, a nearly zero (relatively high) effective spin is typical. In the joint analysis of the primary mass and $χ_{\rm eff}$ distributions, we identify a sub-population of BBHs characterized by both a bump-like primary mass distribution extending to $\sim 100M_\odot$ and a relatively high $χ_{\rm eff}\sim 0.3$, consistent with the expectation for the $\geq 2$g merger events in the disks of active galactic nuclei (AGNs). Our results suggest that the merger rate in AGN disks is around $1.7~{\rm Gpc}^{-3}{\rm yr}^{-1}$. △ Less

Submitted 20 October, 2021; originally announced October 2021.

Comments: 7 pages, 3 figures

TESS Hunt for Young and Maturing Exoplanets (THYME) VI: an 11 Myr giant planet transiting a very low-mass star in Lower Centaurus Crux

Authors: Andrew W. Mann, Mackenna L. Wood, Stephen P. Schmidt, Madyson G. Barber, James E. Owen, Benjamin M. Tofflemire, Elisabeth R. Newton, Eric E. Mamajek, Jonathan L. Bush, Gregory N. Mace, Adam L. Kraus, Pa Chia Thao, Andrew Vanderburg, Joe Llama, Christopher M. Johns-Krull, L. Prato, Asa G. Stahl, Shih-Yun Tang, Matthew J. Fields, Karen A. Collins, Kevin I. Collins, Tianjun Gan, Eric L. N. Jensen, Jacob Kamler, Richard P. Schwarz , et al. (26 additional authors not shown)

Abstract: Mature super-Earths and sub-Neptunes are predicted to be $\simeq$Jovian radius when younger than 10 Myr. Thus, we expect to find 5-15$R_\oplus$ planets around young stars even if their older counterparts harbor none. We report the discovery and validation of TOI 1227 b, a $0.85\pm0.05R_J$ (9.5$R_\oplus$) planet transiting a very low-mass star ($0.170\pm0.015M_\odot$) every 27.4 days. TOI~1227's ki… ▽ More Mature super-Earths and sub-Neptunes are predicted to be $\simeq$Jovian radius when younger than 10 Myr. Thus, we expect to find 5-15$R_\oplus$ planets around young stars even if their older counterparts harbor none. We report the discovery and validation of TOI 1227 b, a $0.85\pm0.05R_J$ (9.5$R_\oplus$) planet transiting a very low-mass star ($0.170\pm0.015M_\odot$) every 27.4 days. TOI~1227's kinematics and strong lithium absorption confirm it is a member of a previously discovered sub-group in the Lower Centaurus Crux OB association, which we designate the Musca group. We derive an age of 11$\pm$2 Myr for Musca, based on lithium, rotation, and the color-magnitude diagram of Musca members. The TESS data and ground-based follow-up show a deep (2.5\%) transit. We use multiwavelength transit observations and radial velocities from the IGRINS spectrograph to validate the signal as planetary in nature, and we obtain an upper limit on the planet mass of $\simeq0.5 M_J$. Because such large planets are exceptionally rare around mature low-mass stars, we suggest that TOI 1227 b is still contracting and will eventually turn into one of the more common $<5R_\oplus$ planets. △ Less

Submitted 9 March, 2022; v1 submitted 18 October, 2021; originally announced October 2021.

Comments: Accepted to the Astronomical Journal. Minor updates during referee process and proofs

Journal ref: AJ 163 156 (2022)

Third Data Release of the Hyper Suprime-Cam Subaru Strategic Program

Authors: Hiroaki Aihara, Yusra AlSayyad, Makoto Ando, Robert Armstrong, James Bosch, Eiichi Egami, Hisanori Furusawa, Junko Furusawa, Sumiko Harasawa, Yuichi Harikane, Bau-Ching Hsieh, Hiroyuki Ikeda, Kei Ito, Ikuru Iwata, Tadayuki Kodama, Michitaro Koike, Mitsuru Kokubo, Yutaka Komiyama, Xiangchong Li, Yongming Liang, Yen-Ting Lin, Robert H. Lupton, Nate B Lust, Lauren A. MacArthur, Ken Mawatari , et al. (42 additional authors not shown)

Abstract: The paper presents the third data release of Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP), a wide-field multi-band imaging survey with the Subaru 8.2m telescope. HSC-SSP has three survey layers (Wide, Deep, and UltraDeep) with different area coverages and depths, designed to address a wide array of astrophysical questions. This third release from HSC-SSP includes data from 278 nights of ob… ▽ More The paper presents the third data release of Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP), a wide-field multi-band imaging survey with the Subaru 8.2m telescope. HSC-SSP has three survey layers (Wide, Deep, and UltraDeep) with different area coverages and depths, designed to address a wide array of astrophysical questions. This third release from HSC-SSP includes data from 278 nights of observing time and covers about 670 square degrees in all five broad-band filters at the full depth ($\sim26$~mag at $5σ$) in the Wide layer. If we include partially observed area, the release covers 1,470 square degrees. The Deep and UltraDeep layers have $\sim80\%$ of the originally planned integration times, and are considered done, as we have slightly changed the observing strategy in order to compensate for various time losses. There are a number of updates in the image processing pipeline. Of particular importance is the change in the sky subtraction algorithm; we subtract the sky on small scales before the detection and measurement stages, which has significantly reduced false detections. Thanks to this and other updates, the overall quality of the processed data has improved since the previous release. However, there are limitations in the data (for example, the pipeline is not optimized for crowded fields), and we encourage the user to check the quality assurance plots as well as a list of known issues before exploiting the data. The data release website is https://hsc-release.mtk.nao.ac.jp/. △ Less

Submitted 30 August, 2021; originally announced August 2021.

Comments: 25 pages, 19 figures, submitted to PASJ. Data available at https://hsc-release.mtk.nao.ac.jp/

Population properties of neutron stars in the coalescing compact binaries

Authors: Yin-Jie Li, Shao-Peng Tang, Yuan-Zhu Wang, Qiang Yuan, Yi-Zhong Fan, Da-Ming Wei

Abstract: We perform a hierarchical Bayesian inference to investigate the population properties of the coalescing compact binaries involving at least one neutron star (NS). With the current observation data, we can not rule out either of the Double Gaussian, Single Gaussian and Uniform NS mass distribution models, although the mass distribution of the Galactic NSs is slightly preferred by the gravitational… ▽ More We perform a hierarchical Bayesian inference to investigate the population properties of the coalescing compact binaries involving at least one neutron star (NS). With the current observation data, we can not rule out either of the Double Gaussian, Single Gaussian and Uniform NS mass distribution models, although the mass distribution of the Galactic NSs is slightly preferred by the gravitational wave (GW) observations. The mass distribution of black holes (BHs) in the neutron star-black hole (NSBH) population is found to be similar to that for the Galactic X-ray binaries. Additionally, the ratio of the merger rate densities between NSBHs and BNSs is estimated to be about 3 : 7. The spin properties of the binaries, though constrained relatively poor, play nontrivial role in reconstructing the mass distribution of NSs and BHs. We find that a perfectly aligned spin distribution can be ruled out, while a purely isotropic distribution of spin orientation is still allowed. △ Less

Submitted 28 October, 2021; v1 submitted 16 August, 2021; originally announced August 2021.

Evidence of Early-stage Tidal Structures of Open Clusters Revealed by Kinematics with Gaia EDR3

Authors: Yezhang Li, Xiaoying Pang, Shih-Yun Tang

Abstract: Blanco 1, a 100Myr open cluster in the solar neighborhood, is well known for its two 50pc-long tidal tails. Taking Blanco 1 as a reference, we find evidence of early-stage tidal disruption in two other open clusters of ~120Myr: the Pleiades and NGC 2516, via Gaia EDR3 data. These two clusters have a total mass of 2-6 times that of Blanco 1. Despite having a similar age as Blanco 1, the Pleiades an… ▽ More Blanco 1, a 100Myr open cluster in the solar neighborhood, is well known for its two 50pc-long tidal tails. Taking Blanco 1 as a reference, we find evidence of early-stage tidal disruption in two other open clusters of ~120Myr: the Pleiades and NGC 2516, via Gaia EDR3 data. These two clusters have a total mass of 2-6 times that of Blanco 1. Despite having a similar age as Blanco 1, the Pleiades and NGC 2516 have a larger fraction of their members bound: 86% of their mass is inside the tidal radius, versus 63% for Blanco 1. However, a correlation between Blanco 1's 50pc-long tidal tails and the "kinematic tails" in velocity space is also found for the Pleiades and NGC 2516. This evidence supports the idea that the modest elongation seen in the spatial distribution for the Pleiades and NGC 2516 is a result of early-stage tidal disruption. △ Less

Submitted 22 July, 2021; originally announced July 2021.

Comments: 4 pages, 1 figure

GRB 200716C: Evidence for a Short Burst Being Lensed

Authors: Yun Wang, Lu-Yao Jiang, Cheng-Kui Li, Jia Ren, Shao-Peng Tang, Zi-Min Zhou, Yun-Feng Liang, Yi-Zhong Fan

Abstract: A tiny fraction of observed gamma-ray bursts (GRBs) may be lensed. The time delays induced by the gravitational lensing are milliseconds to seconds if the point lenses are intermediate-mass black holes. The prompt emission of the lensed GRBs, in principle, should have repeated pulses with identical light curves and spectra but different fluxes and slightly offset positions. In this work, we search… ▽ More A tiny fraction of observed gamma-ray bursts (GRBs) may be lensed. The time delays induced by the gravitational lensing are milliseconds to seconds if the point lenses are intermediate-mass black holes. The prompt emission of the lensed GRBs, in principle, should have repeated pulses with identical light curves and spectra but different fluxes and slightly offset positions. In this work, we search for such candidates within the GRBs detected by Fermi/GBM, Swift/BAT, and HXMT/HE and report the identification of an attractive event GRB 200716C that consists of two pulses. Both the autocorrelation analysis and the Bayesian inference of the prompt emission light curve are in favor of the gravitational lensing scenario. Moreover, the spectral properties of the two pulses are rather similar and follow the so-called Amati relation of short GRBs rather than long duration bursts. The measured flux ratios between the two pulses are nearly constant in all channels, as expected from gravitational lensing. We therefore suggest that the long duration burst GRB 200716C was a short event being lensed. The redshifted mass of the lens was estimated to be $4.25^{+2.46}_{-1.36}$ $\times$ $10^5$ $M_{\odot}$ (90$\%$ credibility). If correct, this could point towards the existence of an intermediate-mass black hole along the line of sight of GRB 200716C. △ Less

Submitted 15 September, 2021; v1 submitted 22 July, 2021; originally announced July 2021.

Comments: 11 pages, 3 figures, 2 tables

Journal ref: ApJL 918 L34 (2021)

The JCMT Transient Survey: Four Year Summary of Monitoring the Submillimeter Variability of Protostars

Authors: Yong-Hee Lee, Doug Johnstone, Jeong-Eun Lee, Gregory Herczeg, Steve Mairs, Carlos Contreras-Peña, Jennifer Hatchell, Tim Naylor, Graham S. Bell, Tyler L. Bourke, Colton Broughton, Logan Francis, Aashish Gupta, Daniel Harsono, Sheng-Yuan Liu, Geumsook Park, Spencer Plovie, Gerald H. Moriarty-Schieven, Aleks Scholz, Tanvi Sharma, Paula Stella Teixeira, Yao-Te Wang, Yuri Aikawa, Geoffrey C. Bower, Huei-Ru Vivien Chen , et al. (27 additional authors not shown)

Abstract: We present the four-year survey results of monthly submillimeter monitoring of eight nearby ($< 500 $pc) star-forming regions by the JCMT Transient Survey. We apply the Lomb-Scargle Periodogram technique to search for and characterize variability on 295 submillimeter peaks brighter than 0.14 Jy beam$^{-1}$, including 22 disk sources (Class II), 83 protostars (Class 0/I), and 190 starless sources.… ▽ More We present the four-year survey results of monthly submillimeter monitoring of eight nearby ($< 500 $pc) star-forming regions by the JCMT Transient Survey. We apply the Lomb-Scargle Periodogram technique to search for and characterize variability on 295 submillimeter peaks brighter than 0.14 Jy beam$^{-1}$, including 22 disk sources (Class II), 83 protostars (Class 0/I), and 190 starless sources. We uncover 18 secular variables, all of them protostars. No single-epoch burst or drop events and no inherently stochastic sources are observed. We classify the secular variables by their timescales into three groups: Periodic, Curved, and Linear. For the Curved and Periodic cases, the detectable fractional amplitude, with respect to mean peak brightness, is $\sim4$ % for sources brighter than $\sim$ 0.5 Jy beam$^{-1}$. Limiting our sample to only these bright sources, the observed variable fraction is 37 % (16 out of 43). Considering source evolution, we find a similar fraction of bright variables for both Class 0 and Class I. Using an empirically motivated conversion from submillimeter variability to variation in mass accretion rate, six sources (7 % of our full sample) are predicted to have years-long accretion events during which the excess mass accreted reaches more than 40 % above the total quiescently accreted mass: two previously known eruptive Class I sources, V1647 Ori and EC 53 (V371 Ser), and four Class 0 sources, HOPS 356, HOPS 373, HOPS 383, and West 40. Considering the full protostellar ensemble, the importance of episodic accretion on few years timescale is negligible, only a few percent of the assembled mass. However, given that this accretion is dominated by events of order the observing time-window, it remains uncertain as to whether the importance of episodic events will continue to rise with decades-long monitoring. △ Less

Submitted 22 July, 2021; originally announced July 2021.

Comments: Accepted for publication in the Astrophysical Journal

Black Hole Gravitational Potential Enhanced Fallback Accretion onto the Nascent Lighter Compact Object: Tentative Evidence in the O3 Run Data of LIGO/Virgo

Authors: Shao-Peng Tang, Yin-Jie Li, Yuan-Zhu Wang, Yi-Zhong Fan, Da-Ming Wei

Abstract: In a binary system, the gravitational potential of the primary black hole may play an important role in enhancing the fallback accretion onto the lighter compact object newly formed in the second supernova explosion. As a result, the final masses of the binary compact objects would be correlated, as suggested recently by Safarzadeh \& Wysocki. In this work, we analyze the mass distribution of four… ▽ More In a binary system, the gravitational potential of the primary black hole may play an important role in enhancing the fallback accretion onto the lighter compact object newly formed in the second supernova explosion. As a result, the final masses of the binary compact objects would be correlated, as suggested recently by Safarzadeh \& Wysocki. In this work, we analyze the mass distribution of four gravitational-wave events, which are characterized by both a small mass ratio and a low mass ($\leq 5M_\odot$) of the light component, and find tentative evidence for a mass correlation among the objects. To evaluate the feasibility of testing such a hypothesis with upcoming observations, we carry out simulations with a mock population and perform Bayesian hierarchical inference for the mass distribution. We find that with dozens of low mass ratio events, whether there exists correlation in the component mass distributions or not can be robustly tested and the correlation, if it exists, can be well determined. △ Less

Submitted 16 November, 2021; v1 submitted 19 July, 2021; originally announced July 2021.

Comments: 8 pages, 3 figures, ApJ published

Journal ref: 2021, ApJ, 922, 3

Disruption of Hierarchical Clustering in the Vela OB2 Complex and the Cluster Pair Collinder 135 and UBC7 with Gaia EDR3: Evidence of Supernova Quenching

Authors: Xiaoying Pang, Zeqiu Yu, Shih-Yun Tang, Jongsuk Hong, Zhen Yuan, Mario Pasquato, M. B. N. Kouwenhoven

Abstract: We identify hierarchical structures in the Vela OB2 complex and the cluster pair Collinder 135 and UBC 7 with Gaia EDR3 using the neural network machine learning algorithm StarGO. Five second-level substructures are disentangled in Vela OB2, which are referred to as Huluwa 1 (Gamma Velorum), Huluwa 2, Huluwa 3, Huluwa 4 and Huluwa 5. For the first time, Collinder 135 and UBC 7 are simultaneously i… ▽ More We identify hierarchical structures in the Vela OB2 complex and the cluster pair Collinder 135 and UBC 7 with Gaia EDR3 using the neural network machine learning algorithm StarGO. Five second-level substructures are disentangled in Vela OB2, which are referred to as Huluwa 1 (Gamma Velorum), Huluwa 2, Huluwa 3, Huluwa 4 and Huluwa 5. For the first time, Collinder 135 and UBC 7 are simultaneously identified as constituent clusters of the pair with minimal manual intervention. We propose an alternative scenario in which Huluwa 1-5 have originated from sequential star formation. The older clusters Huluwa 1-3 with an age of 10-22 Myr, generated stellar feedback to cause turbulence that fostered the formation of the younger-generation Huluwa 4-5 (7-20 Myr). A supernova explosion located inside the Vela IRAS shell quenched star formation in Huluwa 4-5 and rapidly expelled the remaining gas from the clusters. This resulted in global mass stratification across the shell, which is confirmed by the regression discontinuity method. The stellar mass in the lower rim of the shell is $0.32\pm0.14$ $\rm M_\odot$ higher than in the upper rim. Local, cluster-scale mass segregation is observed in the lowest-mass cluster Huluwa 5. Huluwa 1-5 (in Vela OB2) are experiencing significant expansion, while the cluster pair suffers from moderate expansion. The velocity dispersions suggest that all five groups (including Huluwa 1A and Huluwa 1B) in Vela OB2 and the cluster pair are supervirial and are undergoing disruption, and also that Huluwa 1A and Huluwa 1B may be a coeval young cluster pair. N-body simulations predict that Huluwa 1-5 in Vela OB2 and the cluster pair will continue to expand in the future 100 Myr and eventually dissolve. △ Less

Submitted 5 September, 2021; v1 submitted 14 June, 2021; originally announced June 2021.

Comments: 28 pages, 16 figures. Accepted to ApJ

Constraints on the phase transition and nuclear symmetry parameters from PSR $\mathrm{J}0740+6620$ and multimessenger data of other neutron stars

Authors: Shao-Peng Tang, Jin-Liang Jiang, Ming-Zhe Han, Yi-Zhong Fan, Da-Ming Wei

Abstract: Recently, the radius of neutron star (NS) PSR J0740+6620 was measured by Neutron Star Interior Composition Explorer (NICER) and an updated measurement of neutron skin thickness of ${}^{208}$Pb ($R_{\rm skin}^{208}$) was reported by the PREX-II experiment. These new measurements can help us better understand the unknown equation of state (EOS) of dense matter. In this work, we adopt a hybrid parame… ▽ More Recently, the radius of neutron star (NS) PSR J0740+6620 was measured by Neutron Star Interior Composition Explorer (NICER) and an updated measurement of neutron skin thickness of ${}^{208}$Pb ($R_{\rm skin}^{208}$) was reported by the PREX-II experiment. These new measurements can help us better understand the unknown equation of state (EOS) of dense matter. In this work, we adopt a hybrid parameterization method, which incorporates the nuclear empirical parameterization and some widely used phenomenological parameterizations, to analyze the results of nuclear experiments and astrophysical observations. With the joint Bayesian analysis of GW170817, PSR J0030+0451, and PSR J0740+6620, the parameters that characterize the ultradense matter EOS are constrained. We find that the slope parameter $L$ is approximately constrained to $70_{-18}^{+21}$ MeV, which predicts $R_{\rm skin}^{208}=0.204^{+0.030}_{-0.026}\,{\rm fm}$ by using the universal relation between $R_{\rm skin}^{208}$ and $L$. The bulk properties of canonical $1.4\,M_\odot$ NS (e.g., $R_{1.4}$ and $Λ_{1.4}$) as well as the pressure ($P_{2ρ_{\rm sat}}$) at two times the nuclear saturation density are well constrained by the data; i.e., $R_{1.4}$, $Λ_{1.4}$, and $P_{2ρ_{\rm sat}}$ are approximately constrained to $12.3\pm0.7$ km, $330_{-100}^{+140}$, and $4.1_{-1.2}^{+1.5}\times10^{34}\,{\rm dyn\,cm^{-2}}$, respectively. Besides, we find that the Bayes evidences of the hybrid star and normal NS assumptions are comparable, which indicates that current observation data are compatible with quarkyonic matter existing in the core of massive star. Finally, in the case of normal NS assumption, we obtain a constraint for the maximum mass of nonrotating NS $M_{\rm TOV}=2.30^{+0.30}_{-0.18}$ $M_\odot$. All of the uncertainties reported above are for 68.3% credible levels. △ Less

Submitted 21 September, 2021; v1 submitted 8 June, 2021; originally announced June 2021.

Comments: 12 pages, 7 figures, published in PRD

Journal ref: Phys. Rev. D 104, 063032 (2021)

Optical Spectroscopy of Dual Quasar Candidates from the Subaru HSC-SSP program

Authors: Shenli Tang, John D. Silverman, Xuheng Ding, Junyao Li, Khee-Gan Lee, Michael A. Strauss, Andy Goulding, Malte Schramm, Lalitwadee Kawinwanichakij, J. Xavier Prochaska, Joseph F. Hennawi, Masatoshi Imanishi, Kazushi Iwasawa, Yoshiki Toba, Issha Kayo, Masamune Oguri, Yoshiki Matsuoka, Kohei Ichikawa, Tilman Hartwig, Nobunari Kashikawa, Toshihiro Kawaguchi, Kotaro Kohno, Yuichi Matsuda, Tohru Nagao, Yoshiaki Ono , et al. (8 additional authors not shown)

Abstract: We report on a spectroscopic program to search for dual quasars using Subaru Hyper Suprime-Cam (HSC) images of SDSS quasars which represent an important stage during galaxy mergers. Using Subaru/FOCAS and Gemini-N/GMOS, we identify three new physically associated quasar pairs having projected separations less than 20 kpc, out of 26 observed candidates. These include the discovery of the highest re… ▽ More We report on a spectroscopic program to search for dual quasars using Subaru Hyper Suprime-Cam (HSC) images of SDSS quasars which represent an important stage during galaxy mergers. Using Subaru/FOCAS and Gemini-N/GMOS, we identify three new physically associated quasar pairs having projected separations less than 20 kpc, out of 26 observed candidates. These include the discovery of the highest redshift ($z=3.1$) quasar pair with a separation $<$ 10 kpc. Based on the sample acquired to date, the success rate of identifying physically associated dual quasars is $19\%$ when excluding stars based on their HSC colors. Using the full sample of six spectroscopically confirmed dual quasars, we find that the black holes in these systems have black hole masses ($M_{BH} \sim 10^{8-9}M_{\odot}$) similar to single SDSS quasars as well as their bolometric luminosities and Eddington ratios. We measure the stellar mass of their host galaxies based on 2D image decomposition of the five-band ($grizy$) optical emission and assess the mass relation between supermassive black holes (SMBHs) and their hosts. Dual SMBHs appear to have elevated masses relative to their host galaxies. Thus mergers may not necessarily align such systems onto the local mass relation, as suggested by the Horizon-AGN simulation. This study suggests that dual luminous quasars are triggered prior to the final coalescence of the two SMBHs, resulting in early mass growth of the black holes relative to their host galaxies. △ Less

Submitted 21 May, 2021; originally announced May 2021.

Comments: 36 pages, 21 figures, 4 tables

Impacts of Water Latent Heat on the Thermal Structure of Ultra-Cool Objects: Brown Dwarfs and Free-Floating Planets

Authors: Shih-Yun Tang, Tyler D. Robinson, Mark S. Marley, Natasha E. Batalha, Roxana Lupu, L. Prato

Abstract: Brown dwarfs are essential targets for understanding planetary and sub-stellar atmospheres across a wide range of thermal and chemical conditions. As surveys continue to probe ever deeper, and as observing capabilities continue to improve, the number of known Y dwarfs -- the coldest class of sub-stellar objects, with effective temperatures below about 600 K -- is rapidly growing. Critically, this… ▽ More Brown dwarfs are essential targets for understanding planetary and sub-stellar atmospheres across a wide range of thermal and chemical conditions. As surveys continue to probe ever deeper, and as observing capabilities continue to improve, the number of known Y dwarfs -- the coldest class of sub-stellar objects, with effective temperatures below about 600 K -- is rapidly growing. Critically, this class of ultra-cool objects has atmospheric conditions that overlap with Solar System worlds and, as a result, tools and ideas developed from studying Earth, Jupiter, Saturn and other nearby worlds are well-suited for application to sub-stellar atmospheres. To that end, we developed a one-dimensional (vertical) atmospheric structure model for ultra-cool objects that includes moist adiabatic convection, as this is an important process for many Solar System planets. Application of this model across a range of effective temperatures (350, 300, 250, 200 K), metallicities ([M/H] of 0.0, 0.5, 0.7, 1.5), and gravities (log $g$ of 4.0, 4.5, 4.7, 5.0) demonstrates strong impacts of water latent heat release on simulated temperature-pressure profiles. At the highest metallicities, water vapor mixing ratios reach an Earth-like 3%, with associated major alterations to the thermal structure in the atmospheric regions where water condenses. Spectroscopic and photometric signatures of metallicity and moist convection should be readily detectable at near- and mid-infrared wavelengths, especially with James Webb Space Telescope observations, and can help indicate the formation history of an object. △ Less

Submitted 16 August, 2021; v1 submitted 14 May, 2021; originally announced May 2021.

Comments: 16 pages, 10 figures, ApJ accepted version, data public available on zenodo

Quasinormal-modes of the Kerr-Newman black hole: GW150914 and fundamental physics implications

Authors: Hai-Tian Wang, Shao-Peng Tang, Peng-Cheng Li, Yi-Zhong Fan

Abstract: We develop an analytical ringdown waveform model, including both the fundamental and the overtone quasinormal modes, for charged black holes and show that it is precise enough to analyze current gravitational wave data. Applying this waveform model to GW150914, the charge to mass ratio of the remnant black hole ($λ_f$) has been constrained with the sole ringdown gravitational wave data for the fir… ▽ More We develop an analytical ringdown waveform model, including both the fundamental and the overtone quasinormal modes, for charged black holes and show that it is precise enough to analyze current gravitational wave data. Applying this waveform model to GW150914, the charge to mass ratio of the remnant black hole ($λ_f$) has been constrained with the sole ringdown gravitational wave data for the first time and the $90\%$ upper limit is $λ_f\leq 0.38$. Correspondingly, the deviation parameter of the scalar-tensor vector gravity ($α_{\rm s}$) is limited to be $α_{\rm s}\leq 0.17$. Our approach can be directly applied to other binary black hole mergers with loud ringdown radiation. △ Less

Submitted 2 December, 2021; v1 submitted 15 April, 2021; originally announced April 2021.

Comments: 6 pages, 4 figures; Updated to match version published in PRD

Journal ref: Physical Review D 104, 104063 (2021)

Tight constraints on Einstein-dilation-Gauss-Bonnet gravity from GW190412 and GW190814

Authors: Hai-Tian Wang, Shao-Peng Tang, Peng-Cheng Li, Ming-Zhe Han, Yi-Zhong Fan

Abstract: Gravitational-wave (GW) data can be used to test general relativity in the highly nonlinear and strong field regime. Modified gravity theories such as Einstein-dilation-Gauss-Bonnet and dynamical Chern-Simons can be tested with the additional GW signals detected in the first half of the third observing run of Advanced LIGO/Virgo. Specifically, we analyze gravitational-wave data of GW190412 and GW1… ▽ More Gravitational-wave (GW) data can be used to test general relativity in the highly nonlinear and strong field regime. Modified gravity theories such as Einstein-dilation-Gauss-Bonnet and dynamical Chern-Simons can be tested with the additional GW signals detected in the first half of the third observing run of Advanced LIGO/Virgo. Specifically, we analyze gravitational-wave data of GW190412 and GW190814 to place constraints on the parameters of these two theories. Our results indicate that dynamical Chern-Simons gravity remains unconstrained. For Einstein-dilation-Gauss-Bonnet gravity, we find $\sqrt{α_{\rm EdGB}}\lesssim 0.40\,\rm km$ when considering GW190814 data, assuming it is a black hole binary. Such a constraint is improved by a factor of approximately $10$ in comparison to that set by the first Gravitational-Wave Transient Catalog events. △ Less

Submitted 8 July, 2021; v1 submitted 15 April, 2021; originally announced April 2021.

Comments: 9 pages, 4 figures, published in PRD

Journal ref: Phys. Rev. D 104, 024015 (2021)

A flexible Gaussian process reconstruction method and the mass function of the coalescing binary black hole systems

Authors: Yin-Jie Li, Yuan-Zhu Wang, Ming-Zhe Han, Shao-Peng Tang, Qiang Yuan, Yi-Zhong Fan, Da-Ming Wei

Abstract: We develop a new method based on Gaussian process to reconstruct the mass distribution of binary black holes (BBHs). Instead of prespecifying the formalisms of mass distribution, we introduce a more flexible and nonparametric model with which the distribution can be mainly determined by the observed data. We first test our method with simulated data, and find that it can well recover the injected… ▽ More We develop a new method based on Gaussian process to reconstruct the mass distribution of binary black holes (BBHs). Instead of prespecifying the formalisms of mass distribution, we introduce a more flexible and nonparametric model with which the distribution can be mainly determined by the observed data. We first test our method with simulated data, and find that it can well recover the injected distribution. Then we apply this method to analyze the data of BBHs' observations from LIGO-Virgo Gravitational-Wave Transient Catalog 2. By reconstructing the chirp mass distribution, we find that there is a peak or a platform locating at $20-30\,M_{\odot}$ rather than a single-power-law-like decrease from low mass to high mass. Moreover, one or two peaks in the chirp mass range of $\mathcal{M}<20\,M_{\odot}$ may be favored by the data. Assuming a mass-independent mass ratio distribution of $p(q)\propto q^{1.4}$, we further obtain a distribution of primary mass, and find that there is a feature locating in the range of $(30, 40)\,M_{\odot}$, which can be related to \textsc{Broken Power Law} and \textsc{Power Law + Peak} distributions described in The LIGO Scientific Collaboration et al. (2020). Besides, the merger rate of BBHs is estimated to $\mathcal{R} = 26.29^{+14.21}_{-8.96}~{\rm Gpc^{-3}~yr^{-1}}$ supposing there is no redshift evolution. △ Less

Submitted 8 June, 2021; v1 submitted 7 April, 2021; originally announced April 2021.

Black Hole Mass Function of Coalescing Binary Black Hole Systems: Is There a Pulsational Pair Instability Mass Cutoff?

Authors: Yuan-Zhu Wang, Shao-Peng Tang, Yun-Feng Liang, Ming-Zhe Han, Xiang Li, Zhi-Ping Jin, Yi-Zhong Fan, Da-Ming Wei

Abstract: We analyze the LIGO/Virgo GWTC-2 catalog to study the primary mass distribution of the merging black holes. We perform hierarchical Bayesian analysis, and examine whether the mass distribution has a sharp cutoff for primary black hole masses below $65 M_\odot$, as predicted in pulsational pair instability supernova model. We construct two empirical mass functions. One is a piece-wise function with… ▽ More We analyze the LIGO/Virgo GWTC-2 catalog to study the primary mass distribution of the merging black holes. We perform hierarchical Bayesian analysis, and examine whether the mass distribution has a sharp cutoff for primary black hole masses below $65 M_\odot$, as predicted in pulsational pair instability supernova model. We construct two empirical mass functions. One is a piece-wise function with two power-law segments jointed by a sudden drop. The other consists of a main truncated power-law component, a Gaussian component, and a third very massive component. Both models can reasonably fit the data and a sharp drop of the mass distribution is found at $\sim 50M_\odot$, suggesting that the majority of the observed black holes can be explained by the stellar evolution scenarios in which the pulsational pair-instability process takes place. On the other hand, the very massive sub-population, which accounts for at most several percents of the total, may be formed through hierarchical mergers or other processes. △ Less

Submitted 6 April, 2021; originally announced April 2021.

Comments: Accepted for publication in ApJ (This manuscript incorporating the GWTC-2 data is to replace arXiv:2009.03854 that was based on the GWTC-1 data.)

IGRINS RV: A Precision RV Pipeline for IGRINS Using Modified Forward Modeling in the Near-Infrared

Authors: Asa G. Stahl, Shih-Yun Tang, Christopher M. Johns-Krull, L. Prato, Joe Llama, Gregory N. Mace, Jae Joon Lee, Heeyoung Oh, Jessica Luna, Daniel T. Jaffe

Abstract: Application of the radial velocity (RV) technique in the near infrared is valuable because of the diminished impact of stellar activity at longer wavelengths, making it particularly advantageous for the study of late-type stars but also for solar-type objects. In this paper, we present the IGRINS RV open source python pipeline for computing infrared RV measurements from reduced spectra taken with… ▽ More Application of the radial velocity (RV) technique in the near infrared is valuable because of the diminished impact of stellar activity at longer wavelengths, making it particularly advantageous for the study of late-type stars but also for solar-type objects. In this paper, we present the IGRINS RV open source python pipeline for computing infrared RV measurements from reduced spectra taken with IGRINS, a R ~ 45,000 spectrograph with simultaneous coverage of the H band (1.49--1.80 $μ$m) and K band (1.96--2.46 $μ$m). Using a modified forward modeling technique, we construct high resolution telluric templates from A0 standard observations on a nightly basis to provide a source of common-path wavelength calibration while mitigating the need to mask or correct for telluric absorption. Telluric standard observations are also used to model the variations in instrumental resolution across the detector, including a yearlong period when the K band was defocused. Without any additional instrument hardware, such as a gas cell or laser frequency comb, we are able to achieve precisions of 26.8 $\rm m\,s^{-1}$ in the K band and 31.1 $\rm m\,s^{-1}$ in the H band for narrow-line hosts. These precisions are empirically determined by a monitoring campaign of two RV standard stars as well as the successful retrieval of planet-induced RV signals for both HD 189733 and $τ$ Boo A; furthermore, our results affirm the presence of the Rossiter-McLaughlin effect for HD 189733. The IGRINS RV pipeline extends another important science capability to IGRINS, with publicly available software designed for widespread use. △ Less

Submitted 7 April, 2021; v1 submitted 5 April, 2021; originally announced April 2021.

Comments: 35 pages, 14 figures, accepted to AJ