Showing 1–50 of 480 results for author: Chen, Z

ISMGCC: Finding Gas Structures in Molecular Interstellar Medium Using Gaussian Decomposition and Graph Theory

Authors: Haoran Feng, Zhiwei Chen, Zhibo Jiang, James S. Urquhart

Abstract: Molecular line emissions are commonly used to trace the distribution and properties of molecular Interstellar Medium (ISM). However, the emissions are heavily blended on the Galactic disk toward the inner Galaxy because of the relatively large line widths and the velocity overlaps of spiral arms. Structure identification methods based on voxel connectivity in PPV data cubes often produce unrealist… ▽ More Molecular line emissions are commonly used to trace the distribution and properties of molecular Interstellar Medium (ISM). However, the emissions are heavily blended on the Galactic disk toward the inner Galaxy because of the relatively large line widths and the velocity overlaps of spiral arms. Structure identification methods based on voxel connectivity in PPV data cubes often produce unrealistically large structures, which is the ``over-linking'' problem. Therefore, identifying molecular cloud structures in these directions is not trivial. We propose a new method based on Gaussian decomposition and graph theory to solve the over-linking problem, named ISMGCC (InterStellar Medium Gaussian Component Clustering). Using the MWISP ${}^{13}\mathrm{CO}~(1-0)$ data in the range of $13.5^{\circ} \leq l \leq 14.5^{\circ}, |b| \leq 0.5^{\circ}$, and $-100\leq V_{\mathrm{LSR}} \leq +200~\mathrm{km~s^{-1}}$, our method identified three hundred molecular gas structures with at least 16 pixels. These structures contain $92\%$ of the total flux in the raw data cube. Meanwhile, these structures show single-peaked line profiles on more than $93\%$ of their pixels. This indicates that the ISMGCC method could distinguish gas structures in crowded regions with minimal flux loss. △ Less

Submitted 2 September, 2024; originally announced September 2024.

Comments: 22 pages, 13 figures, submitted to RAA

On the Cosmic Variance of the Merger Rate Density of Binary Neutron Stars

Authors: Zhiwei Chen, Youjun Lu, Jie Wang, Zhen Jiang, Qingbo Chu, Xianghao Ma

Abstract: The cosmic variance on the star formation history may lead to bias to the merger rate density estimation of binary neutron star (BNS) mergers by the compact binary population synthesis. In this paper, we take the advantage of the large boxsize of the Millennium Simulation combined with the semi-analytic galaxy formation model GABE, and the parameterized population binary star evolution (BSE) model… ▽ More The cosmic variance on the star formation history may lead to bias to the merger rate density estimation of binary neutron star (BNS) mergers by the compact binary population synthesis. In this paper, we take the advantage of the large boxsize of the Millennium Simulation combined with the semi-analytic galaxy formation model GABE, and the parameterized population binary star evolution (BSE) model to examine how much effect will the cosmic variance introduce on the estimation of merger rate density of BNS mergers. We find that for sub-box size of $100\rm Mpc$ and $200\rm Mpc$, the variance of merger rate density $σ_{\rm R}/\rm R$ at different redshift is about $23\%-35\%$ and $13\%-20\%$ respectively. On one hand, as for the variance of the detection rate on BNS mergers with current LIGO-Virgo-KAGRA (LVK) detector network, this value is very small $\lesssim 10\%$, which indicates ignoring the cosmic variance is reasonable for estimating the merger rate density from current LVK observation. On the other hand, with next-generation gravitational wave detectors, it is possible to localize BNS mergers within sub-boxes possessing length of $\rm 40 Mpc$ for source redshift $z_{s}<0.2$. In such a small box, the cosmic variance of the merger rate density is significant, i.e., the value of $σ_{\rm R}/\rm R$ is about $\sim 55\%$. This hints that estimating the merger rate density of BNS in different sky areas may provide useful information on the cosmic variance. △ Less

Submitted 29 August, 2024; originally announced August 2024.

Comments: 7 pages, 5 figures, Accepted for Publication in ApJ

Minute-Cadence Observations of the LAMOST Fields with the TMTS: IV -- Catalog of Cataclysmic Variables from the First 3-yr Survey

Authors: Qichun Liu, Jie Lin, Xiaofeng Wang, Zhibin Dai, Yongkang Sun, Gaobo Xi, Jun Mo, Jialian Liu, Shengyu Yan, Alexei V. Filippenko, Thomas G. Brink, Yi Yang, Kishore C. Patra, Yongzhi Cai, Zhihao Chen, Liyang Chen, Fangzhou Guo, Xiaojun Jiang, Gaici Li, Wenxiong Li, Weili Lin, Cheng Miao, Xiaoran Ma, Haowei Peng, Qiqi Xia , et al. (2 additional authors not shown)

Abstract: The Tsinghua University--Ma Huateng Telescopes for Survey (TMTS) started to monitor the LAMOST plates in 2020, leading to the discovery of numerous short-period eclipsing binaries, peculiar pulsators, flare stars, and other variable objects. Here, we present the uninterrupted light curves for a sample of 64 cataclysmic variables (CVs) observed/discovered using the TMTS during its first three-year… ▽ More The Tsinghua University--Ma Huateng Telescopes for Survey (TMTS) started to monitor the LAMOST plates in 2020, leading to the discovery of numerous short-period eclipsing binaries, peculiar pulsators, flare stars, and other variable objects. Here, we present the uninterrupted light curves for a sample of 64 cataclysmic variables (CVs) observed/discovered using the TMTS during its first three-year observations, and we introduce new CVs and new light-variation periods (from known CVs) revealed through the TMTS observations. Thanks to the high-cadence observations of TMTS, diverse light variations, including superhumps, quasi-periodic oscillations, large-amplitude orbital modulations, and rotational modulations, are able to be detected in our CV samples, providing key observational clues for understanding the fast-developing physical processes in various CVs. All of these short-timescale light-curve features help further classify the subtypes of CV systems. We highlight the light-curve features observed in our CV sample and discuss further implications of minute-cadence light curves for CV identifications and classifications. Moreover, we examine the H$α$ emission lines in the spectra from our nonmagnetic CV samples (i.e., dwarf novae and nova-like subclasses) and find that the distribution of H$α$ emission strength shows significant differences between the sources with orbital periods above and below the period gap, which agrees with the trend seen from the SDSS nonmagnetic CV sample. △ Less

Submitted 21 August, 2024; originally announced August 2024.

Comments: 27 pages, 12 figures in main text, accepted for the publication in Universe

The Transition from Galaxy-wide Gas Inflow to Outflow in Quasar Host Galaxies

Authors: Zhicheng He, Zhifu Chen, Guilin Liu, Tinggui Wang, Luis C. Ho, Junxian Wang, Weihao Bian, Zheng Cai, Guobin Mou, Qiusheng Gu, Zhiwen Wang

Abstract: Galactic-wide outflows driven by active galactic nuclei (AGNs) is a routinely invoked feedback mechanism in galaxy evolution models. Hitherto, the interplay among the interstellar gas on galactic scales, the propagation of AGN outflows and the fundamental AGN parameters during evolution remains elusive. Powerful nuclear outflows are found to favorably exist at early AGN stages usually associated w… ▽ More Galactic-wide outflows driven by active galactic nuclei (AGNs) is a routinely invoked feedback mechanism in galaxy evolution models. Hitherto, the interplay among the interstellar gas on galactic scales, the propagation of AGN outflows and the fundamental AGN parameters during evolution remains elusive. Powerful nuclear outflows are found to favorably exist at early AGN stages usually associated with high accretion rates and weak narrow emission lines. In a sample of quasars emitting Mg II narrow absorption lines (NALs) from the Sloan Digital Sky Survey, we discover an unprecedented phenomenon where galaxy-scale inflow-dominated transforming into outflow-dominated gas accompanied by an increasing strength of the narrow [O III] line, at a confidence level of 6.7σ. The fact that nuclear outflows diminish while galaxy-wide outflows intensifies as AGNs evolve implies that early-stage outflows interact with interstellar medium on galactic scales and trigger the gradual transformation into galaxy-wide outflows, providing observational links to the hypothetical multi-stage propagation of AGN outflows that globally regulates galaxy evolution. △ Less

Submitted 8 August, 2024; originally announced August 2024.

Comments: Accepted by SCIENCE CHINA Physics, Mechanics & Astronomy; 15 pages, 4 figures and 6 appendix figures

Revisiting the fundamental metallicity relation with observation and simulation

Authors: Chengyu Ma, Kai Wang, Enci Wang, Yingjie Peng, Haochen Jiang, Haoran Yu, Cheng Jia, Zeyu Chen, Haixin Li, Xu Kong

Abstract: The gas-phase metallicity of galaxies is regulated by multiple astrophysical processes, which makes it a crucial diagnostic of galaxy formation and evolution. Beyond the fundamental mass-metallicity relation, a debate about the secondary galaxy property to predict the metallicity of galaxies arises. Motivated by this, we systematically examine the relationship between gas-phase metallicity and oth… ▽ More The gas-phase metallicity of galaxies is regulated by multiple astrophysical processes, which makes it a crucial diagnostic of galaxy formation and evolution. Beyond the fundamental mass-metallicity relation, a debate about the secondary galaxy property to predict the metallicity of galaxies arises. Motivated by this, we systematically examine the relationship between gas-phase metallicity and other galaxy properties, i.e. star formation rate (SFR) and galaxy size, in addition to stellar mass in both observation and simulation. We utilize the data from the MaNGA (Mapping Nearby Galaxies at Apache Point Observatory) survey and the TNG50 simulations. We find that the combination of $M_*/R_{\rm e}^β$ with $β\sim 0.6-1$ is in much stronger correlation to the metallicity than stellar mass alone, regardless of whether the SFR is included or not, in both observation and simulation. This indicates that galaxy size plays a more important role in determining gas-phase metallicity of galaxies than SFR. In addition, the TNG simulation predicts that the SFR, although being a subdominant role, becomes increasingly important in high-$z$ universe. Finally, we speculate that SFR modulates metallicity on the temporal dimension, synchronized with time-varying gas inflows, and galaxy size regulates metallicity on the spatial dimension by affecting the gravitational potential and the mass loading factor. △ Less

Submitted 31 July, 2024; originally announced July 2024.

Comments: 11 pages, 6 figures, ApJL accepted

MeerKLASS L-band deep-field intensity maps: entering the HI dominated regime

Authors: MeerKLASS Collaboration, Matilde Barberi-Squarotti, José L. Bernal, Philip Bull, Stefano Camera, Isabella P. Carucci, Zhaoting Chen, Steven Cunnington, Brandon N. Engelbrecht, José Fonseca, Keith Grainge, Melis O. Irfan, Yichao Li, Aishrila Mazumder, Sourabh Paul, Alkistis Pourtsidou, Mario G. Santos, Marta Spinelli, Jingying Wang, Amadeus Witzemann, Laura Wolz

Abstract: We present results from MeerKAT single-dish HI intensity maps, the final observations to be performed in L-band in the MeerKAT Large Area Synoptic Survey (MeerKLASS) campaign. The observations represent the deepest single-dish HI intensity maps to date, produced from 41 repeated scans over $236\,{\rm deg}^2$, providing 62 hours of observational data for each of the 64 dishes before flagging. By in… ▽ More We present results from MeerKAT single-dish HI intensity maps, the final observations to be performed in L-band in the MeerKAT Large Area Synoptic Survey (MeerKLASS) campaign. The observations represent the deepest single-dish HI intensity maps to date, produced from 41 repeated scans over $236\,{\rm deg}^2$, providing 62 hours of observational data for each of the 64 dishes before flagging. By introducing an iterative self-calibration process, the estimated thermal noise of the reconstructed maps is limited to ${\sim}\,1.21\,$mK ($1.2\,\times$ the theoretical noise level). This thermal noise will be sub-dominant relative to the HI fluctuations on large scales ($k\,{\lesssim}\,0.15\,h\,\text{Mpc}^{-1}$), which demands upgrades to power spectrum analysis techniques, particularly for covariance estimation. In this work, we present the improved MeerKLASS analysis pipeline, validating it on both a suite of mock simulations and a small sample of overlapping spectroscopic galaxies from the Galaxy And Mass Assembly (GAMA) survey. Despite only overlapping with ${\sim}\,25\%$ of the MeerKLASS deep field, and a conservative approach to covariance estimation, we still obtain a ${>}\,4\,σ$ detection of the cross-power spectrum between the intensity maps and the 2269 galaxies at the narrow redshift range $0.39\,{<}\,z\,{<}\,0.46$. We briefly discuss the HI auto-power spectrum from this data, the detection of which will be the focus of follow-up work. For the first time with MeerKAT single-dish intensity maps, we also present evidence of HI emission from stacking the maps onto the positions of the GAMA galaxies. △ Less

Submitted 31 July, 2024; originally announced July 2024.

Comments: 29 pages, 27 figures. See Fig.14 and 15 for auto-HI and HI-galaxy cross power spectra. Fig 19 and 20 for stacking detections. Submitted to MNRAS

Reducing instrumental errors in Parkes Pulsar Timing Array data

Authors: Axl F. Rogers, Willem van Straten, Sergei Gulyaev, Aditya Parthasarathy, George Hobbs, Zu-Cheng Chen, Yi Feng, Boris Goncharov, Agastya Kapur, Xiaojin Liu, Daniel Reardon, Christopher J. Russell, Andrew Zic

Abstract: This paper demonstrates the impact of state-of-the-art instrumental calibration techniques on the precision of arrival times obtained from 9.6 years of observations of millisecond pulsars using the Murriyang 64-m CSIRO Parkes Radio Telescope. Our study focuses on 21-cm observations of 25 high-priority pulsars that are regularly observed as part of the Parkes Pulsar Timing Array (PPTA) project, inc… ▽ More This paper demonstrates the impact of state-of-the-art instrumental calibration techniques on the precision of arrival times obtained from 9.6 years of observations of millisecond pulsars using the Murriyang 64-m CSIRO Parkes Radio Telescope. Our study focuses on 21-cm observations of 25 high-priority pulsars that are regularly observed as part of the Parkes Pulsar Timing Array (PPTA) project, including those predicted to be the most susceptible to calibration errors. We employ Measurement Equation Template Matching (METM) for instrumental calibration and Matrix Template Matching (MTM) for arrival time estimation, resulting in significantly improved timing residuals with up to a sixfold reduction in white noise compared to arrival times estimated using Scalar Template Matching and conventional calibration based on the Ideal Feed Assumption. The median relative reduction in white noise is 33 percent, and the maximum absolute reduction is 4.5 microseconds. For PSR J0437-4715, METM and MTM reduce the best-fit power-law amplitude (2.7 sigma) and spectral index (1.7 sigma) of the red noise in the arrival time residuals, which can can be tentatively interpreted as mitigation of 1/f noise due to otherwise unmodeled steps in polarimetric response. These findings demonstrate the potential to directly enhance the sensitivity of pulsar timing array experiments through more accurate methods of instrumental calibration and arrival time estimation. △ Less

Submitted 29 July, 2024; originally announced July 2024.

Comments: 30 pages, 9 figures, accepted for publication in The Astrophysical Journal

Deep rest-UV JWST/NIRSpec spectroscopy of early galaxies: the demographics of CIV and N-emitters in the reionization era

Authors: Michael W. Topping, Daniel P. Stark, Peter Senchyna, Zuyi Chen, Adi Zitrin, Ryan Endsley, Stéphane Charlot, Lukas J. Furtak, Michael V. Maseda, Adele Plat, Renske Smit, Ramesh Mainali, Jacopo Chevallard, Stephen Molyneux, Jane R. Rigby

Abstract: JWST has recently discovered a subset of reionization era galaxies with ionized gas that is metal poor in oxygen and carbon but heavily-enriched in nitrogen. This abundance pattern is almost never seen in lower redshift galaxies but is commonly observed in globular cluster stars. We have recently demonstrated that this peculiar abundance pattern appears in a compact ($\simeq 20$ pc) metal-poor gal… ▽ More JWST has recently discovered a subset of reionization era galaxies with ionized gas that is metal poor in oxygen and carbon but heavily-enriched in nitrogen. This abundance pattern is almost never seen in lower redshift galaxies but is commonly observed in globular cluster stars. We have recently demonstrated that this peculiar abundance pattern appears in a compact ($\simeq 20$ pc) metal-poor galaxy undergoing a strong burst of star formation. This galaxy was originally selected based on strong CIV emission, indicating a hard radiation field rarely seen locally. In this paper, we present JWST/NIRSpec observations of another reionization-era galaxy known to power strong CIV emission, the $z=7.04$ gravitationally-lensed galaxy A1703-zd6. The emission line spectrum reveals this is a metal poor galaxy ($12+\log(\rm O/H) = 7.47\pm0.19$) dominated by a young stellar population ($1.6^{+0.5}_{-0.4}$ Myr) that powers a very hard ionizing spectrum (CIV EW = 19.4 $\unicode{x212B}$, He II EW = 2.2 $\unicode{x212B}$). The ISM is highly-enriched in nitrogen ($\log(\rm N/O)=-0.6$) with very high electron densities ($8-19\times10^4$ cm$^{-3}$) and extreme ionization conditions rarely seen at lower redshift. We also find intense CIV emission (EW$\gtrsim20$ $\unicode{x212B}$) in two new $z\gtrsim 6$ metal poor galaxies. To put these results in context, we search for UV line emission in a sample of 737 $z\gtrsim 4$ galaxies with NIRSpec spectra, establishing that 40(30)% of systems with [OIII]+H$β$ EW $>2000\unicode{x212B}$ have NIV] (CIV) detections with EW$>5$ $\unicode{x212B}$ ($>10$ $\unicode{x212B}$). These results suggest high N/O ratios and hard ionizing sources appear in a brief phase following a burst of star formation in compact high density stellar complexes. △ Less

Submitted 26 July, 2024; originally announced July 2024.

Comments: 29 pages, 18 figures; submitted to ApJ

Observational Evidence of Particles Acceleration by Relativistic Magnetic Reconnection in Gamma-ray Bursts

Authors: Cheng-Feng Peng, Rui-Jing Lu, Wen-Qiang Liang, Zhe-geng Chen

Abstract: Gamma-ray bursts (GRBs) as the most energetic explosions in the modern universe have been studied over half a century, but the physics of the particle acceleration and radiation responsible for their observed spectral behaviors are still not well understood. Based on the comprehensive analysis of the pulse properties in both bright GRB~160625B and GRB~160509A, for the first time, we identify evide… ▽ More Gamma-ray bursts (GRBs) as the most energetic explosions in the modern universe have been studied over half a century, but the physics of the particle acceleration and radiation responsible for their observed spectral behaviors are still not well understood. Based on the comprehensive analysis of the pulse properties in both bright GRB~160625B and GRB~160509A, for the first time, we identify evidences of particle acceleration by relativistic magnetic reconnection from the evolutionary behavior of the two spectral breaks ($E_{\rm p}$ and $E_{\rm cut}$). Meanwhile, the adiabatic cooling process of the emitting particles in the magnetic reconnection regions produces a relation between the spectral index and the flux. We also discuss the physics behind spectral energy correlations. Finally, we argue that the identification of an anticorrelation between $E_{\rm cut}$ and $L_{\rm iso}$ may opens a new avenue for diagnostics of the physics of the particle acceleration and radiation in a variety of astrophysical sources. △ Less

Submitted 19 July, 2024; originally announced July 2024.

Comments: 17 pages, 8 figures. Accepted for publication in ApJ

Revealing Gas Inflows towards the Galactic Central Molecular Zone

Authors: Yang Su, Shiyu Zhang, Yan Sun, Ji Yang, Qing-Zeng Yan, Shaobo Zhang, Zhiwei Chen, Xuepeng Chen, Xin Zhou, Lixia Yuan

Abstract: We study the gas inflows towards the Galactic Central Molecular Zone (CMZ) based on the gas morphological and kinematic features from the MWISP in the region of l=1.2 deg--19.0 deg and |b|<3.0 deg. We find that the near dust lane extends to l~15 deg, in which the end of the large-scale gas structure intersects with the 3 kpc-ring at a distance of ~5 kpc. Intriguingly, many filamentary MCs, togethe… ▽ More We study the gas inflows towards the Galactic Central Molecular Zone (CMZ) based on the gas morphological and kinematic features from the MWISP in the region of l=1.2 deg--19.0 deg and |b|<3.0 deg. We find that the near dust lane extends to l~15 deg, in which the end of the large-scale gas structure intersects with the 3 kpc-ring at a distance of ~5 kpc. Intriguingly, many filamentary MCs, together with the bow-like/ballistic-like clouds and continuous CO features with notable velocity gradient, are finely outlined along the long structure. These MCs also have relatively large velocity dispersions, indicating the shocked gas generated by local continuous accretion and thus the enhanced turbulence along the entire gas structure. We suggest that the ~3.1--3.6 kpc long CO structure originates from the accretion molecular gas driven by the Galactic bar. The gas near the bar end at the 3 kpc-ring becomes an important reservoir for the large-scale accreting flows inwards to the CMZ through the bar channel. The inclination angle of the bar is estimated to be 20--26 deg, while the pattern speed of the bar is 30--35 km/s. The total mass of the whole gas lane is about (0.9-1.7)x10^7 Msun according to the calculated X_CO=(0.6-1.4)x10^20 cm^-2 (Kkm/s)^-1 from the large-scale CO data and the complementary HI data. The mean gas inflow rate is about 0.8-1.4 Msun/yr, which seems to be comparable to the outflow's rate of the Galactic nuclear winds after applying the updated lower X-factor value above. △ Less

Submitted 4 August, 2024; v1 submitted 15 July, 2024; originally announced July 2024.

Comments: 16 pages, 3 figures

Journal ref: 2024,ApJL,971,L6

The neutron star mass, distance, and inclination from precision timing of the brilliant millisecond pulsar J0437$-$4715

Authors: Daniel J. Reardon, Matthew Bailes, Ryan M. Shannon, Chris Flynn, Jacob Askew, N. D. Ramesh Bhat, Zu-Cheng Chen, Małgorzata Curyło, Yi Feng, George B. Hobbs, Agastya Kapur, Matthew Kerr, Xiaojin Liu, Richard N. Manchester, Rami Mandow, Saurav Mishra, Christopher J. Russell, Mohsen Shamohammadi, Lei Zhang, Andrew Zic

Abstract: The observation of neutron stars enables the otherwise impossible study of fundamental physical processes. The timing of binary radio pulsars is particularly powerful, as it enables precise characterization of their (three-dimensional) positions and orbits. PSR~J0437$-$4715 is an important millisecond pulsar for timing array experiments and is also a primary target for the Neutron Star Interior Co… ▽ More The observation of neutron stars enables the otherwise impossible study of fundamental physical processes. The timing of binary radio pulsars is particularly powerful, as it enables precise characterization of their (three-dimensional) positions and orbits. PSR~J0437$-$4715 is an important millisecond pulsar for timing array experiments and is also a primary target for the Neutron Star Interior Composition Explorer (NICER). The main aim of the NICER mission is to constrain the neutron star equation of state by inferring the compactness ($M_p/R$) of the star. Direct measurements of the mass $M_p$ from pulsar timing therefore substantially improve constraints on the radius $R$ and the equation of state. Here we use observations spanning 26 years from Murriyang, the 64-m Parkes radio telescope, to improve the timing model for this pulsar. Among the new precise measurements are the pulsar mass $M_p=1.418\pm 0.044$ $M_{\odot}$, distance $D=156.96 \pm 0.11$ pc, and orbital inclination angle $i=137.506 \pm 0.016^\circ$, which can be used to inform the X-ray pulse profile models inferred from NICER observations. We demonstrate that these results are consistent between multiple data sets from the Parkes Pulsar Timing Array (PPTA), each modeled with different noise assumptions. Using the longest available PPTA data set, we measure an apparent second derivative of the pulsar spin frequency and discuss how this can be explained either by kinematic effects due to the proper motion and radial velocity of the pulsar or excess low-frequency noise such as a gravitational-wave background. △ Less

Submitted 30 July, 2024; v1 submitted 9 July, 2024; originally announced July 2024.

Comments: 13 pages, 3 figures, accepted for publication in Astrophysical Journal Letters

Submillimeter and Mid-Infrared Variability of Young Stellar Objects in the M17SWex Intermediate-Mass Star-Forming Region

Authors: Geumsook Park, Doug Johnstone, Carlos Contreras Pena, Jeong-Eun Lee, Sheng-Yuan Liu, Gregory Herczeg, Steve Mairs, Zhiwei Chen, Jennifer Hatchell, Kee-Tae Kim, Mi-Ryang Kim, Keping Qiu, Yao-Te Wang, Xu Zhang, The JCMT Transient Team

Abstract: We present a comprehensive analysis of young stellar object (YSO) variability within the M17 Southwest Extension (M17 SWex), using 3.5 years of monitoring data from the JCMT Transient Survey at sub-millimeter (sub-mm) and 9 years from the NEOWISE mission at mid-infrared (mid-IR). Our study encompasses observations of 147 bright sub-mm peaks identified within our deep JCMT co-added map as well as 1… ▽ More We present a comprehensive analysis of young stellar object (YSO) variability within the M17 Southwest Extension (M17 SWex), using 3.5 years of monitoring data from the JCMT Transient Survey at sub-millimeter (sub-mm) and 9 years from the NEOWISE mission at mid-infrared (mid-IR). Our study encompasses observations of 147 bright sub-mm peaks identified within our deep JCMT co-added map as well as 156 YSOs in NEOWISE W1 and 179 in W2 that were previously identified in Spitzer surveys. We find three robust sub-mm variables: two are candidate YSOs and one is a likely extragalactic source. At mid-IR wavelengths, our analysis reveals secular and stochastic variability in 47 YSOs, with the highest fraction of secular variability occurring at the earliest evolutionary stage. This is similar to what has previously been observed for low-mass YSO variability within the Gould Belt. However, we observe less overall variability in M17SWex at both the sub-mm and mid-IR. We suspect that this lower fraction is due to the greater distance to M17 SWex. Our findings showcase the utility of multi-wavelength observations to better capture the complex variability phenomena inherent to star formation processes and demonstrate the importance of years-long monitoring of a diverse selection of star-forming environments. △ Less

Submitted 3 July, 2024; originally announced July 2024.

Comments: Accepted for Publication in The Astronomical Journal

Probing the connection between IceCube neutrinos and MOJAVE AGN

Authors: R. Abbasi, M. Ackermann, J. Adams, S. K. Agarwalla, J. A. Aguilar, M. Ahlers, J. M. Alameddine, N. M. Amin, K. Andeen, C. Argüelles, Y. Ashida, S. Athanasiadou, L. Ausborm, S. N. Axani, X. Bai, A. Balagopal V., M. Baricevic, S. W. Barwick, S. Bash, V. Basu, R. Bay, J. J. Beatty, J. Becker Tjus, J. Beise, C. Bellenghi , et al. (399 additional authors not shown)

Abstract: Active Galactic Nuclei (AGN) are prime candidate sources of the high-energy, astrophysical neutrinos detected by IceCube. This is demonstrated by the real-time multi-messenger detection of the blazar TXS 0506+056 and the recent evidence of neutrino emission from NGC 1068 from a separate time-averaged study. However, the production mechanism of the astrophysical neutrinos in AGN is not well establi… ▽ More Active Galactic Nuclei (AGN) are prime candidate sources of the high-energy, astrophysical neutrinos detected by IceCube. This is demonstrated by the real-time multi-messenger detection of the blazar TXS 0506+056 and the recent evidence of neutrino emission from NGC 1068 from a separate time-averaged study. However, the production mechanism of the astrophysical neutrinos in AGN is not well established which can be resolved via correlation studies with photon observations. For neutrinos produced due to photohadronic interactions in AGN, in addition to a correlation of neutrinos with high-energy photons, there would also be a correlation of neutrinos with photons emitted at radio wavelengths. In this work, we perform an in-depth stacking study of the correlation between 15 GHz radio observations of AGN reported in the MOJAVE XV catalog, and ten years of neutrino data from IceCube. We also use a time-dependent approach which improves the statistical power of the stacking analysis. No significant correlation was found for both analyses and upper limits are reported. When compared to the IceCube diffuse flux, at 100 TeV and for a spectral index of 2.5, the upper limits derived are $\sim3\%$ and $\sim9\%$ for the time-averaged and time-dependent case, respectively. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Comments: 14 Pages 7 Figures

Reionization Parameter Inference from 3D Minkowski Functionals of the 21 cm Signals

Authors: Kangning Diao, Zhaoting Chen, Xuelei Chen, Yi Mao

Abstract: The Minkowski Functionals (MFs), a set of topological summary statistics, have emerged as a powerful tool for extracting non-Gaussian information. We investigate the prospect of constraining the reionization parameters using the MFs of the 21 cm brightness temperature field from the epoch of reionization (EoR). Realistic effects, including thermal noise, synthesized beam, and foreground avoidance,… ▽ More The Minkowski Functionals (MFs), a set of topological summary statistics, have emerged as a powerful tool for extracting non-Gaussian information. We investigate the prospect of constraining the reionization parameters using the MFs of the 21 cm brightness temperature field from the epoch of reionization (EoR). Realistic effects, including thermal noise, synthesized beam, and foreground avoidance, are applied to the mock observations from the radio interferometric array experiments such as the Hydrogen Epoch of Reionization Array (HERA) and the Square Kilometre Array (SKA). We demonstrate that the MFs of the 21 cm signal measured with SKA-Low can be used to distinguish different reionization models, whereas the MF measurement with a HERA-like array cannot be made accurately enough. We further forecast the accuracies with which the MF measurements can place constraints on reionization parameters, using the standard MCMC analysis for parameter inference based on forward modeling. We find that for SKA-Low observation, MFs provide unbiased estimations of the reionization parameters with accuracies comparable to the power spectrum (PS) analysis. Furthermore, joint constraints using both MFs and PS can improve the constraint accuracies by up to $30\%$ compared to those with the PS alone. Nevertheless, the constraint accuracies can be degraded if the EoR window is shrunk with strong foreground avoidance. Our analysis demonstrates the promise of MFs as a set of summary statistics that extract complementary information from the 21 cm EoR field to the two-point statistics, which suggests a strong motivation for incorporating the MFs into the data analysis of future 21 cm observations. △ Less

Submitted 28 June, 2024; originally announced June 2024.

Comments: 20 pages, 10 figures, submitted to ApJ, comments welcome

Design and Implementation of a Scalable Correlator Based on ROACH2+GPU Cluster for Tianlai 96-Dual-Polarization Antenna Array

Authors: Zhao Wang, Ji-Xia Li, Ke Zhang, 1 Feng-Quan Wu, Hai-Jun Tian, Chen-Hui Niu, Ju-Yong Zhang, Zhi-Ping Chen, Dong-Jin Yu, Xue-Lei Chen

Abstract: The digital correlator is one of the most crucial data processing components of a radio telescope array. With the scale of radio interferometeric array growing, many efforts have been devoted to developing a cost-effective and scalable correlator in the field of radio astronomy. In this paper, a 192-input digital correlator with six CASPER ROACH2 boards and seven GPU servers has been deployed as t… ▽ More The digital correlator is one of the most crucial data processing components of a radio telescope array. With the scale of radio interferometeric array growing, many efforts have been devoted to developing a cost-effective and scalable correlator in the field of radio astronomy. In this paper, a 192-input digital correlator with six CASPER ROACH2 boards and seven GPU servers has been deployed as the digital signal processing system for Tianlai cylinder pathfinder located in Hongliuxia observatory. The correlator consists of 192 input signals (96 dual-polarization), 125-MHz bandwidth, and full-Stokes output. The correlator inherits the advantages of the CASPER system, for example, low cost, high performance, modular scalability, and a heterogeneous computing architecture. With a rapidly deployable ROACH2 digital sampling system, a commercially expandable 10 Gigabit switching network system, and a flexible upgradable GPU computing system, the correlator forms a low-cost and easily-upgradable system, poised to support scalable large-scale interferometeric array in the future. △ Less

Submitted 27 June, 2024; originally announced June 2024.

Comments: 12 pages, 13 figures, accepted for publication in Frontiers in Astronomy and Space Sciences, section Astronomical Instrumentation

The Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER). VI. The High-Mass Stellar Initial Mass Function of M33

Authors: Tobin M. Wainer, Benjamin F. Williams, L. Clifton Johnson, Daniel R. Weisz, Julianne J. Dalcanton, Anil C. Seth, Andrew Dolphin, Meredith J. Durbin, Eric F. Bell, Zhuo Chen, Puragra Guhathakurta, Eric W. Koch, Christina W. Lindberg, Erik Rosolowsky, Karin M. Sandstrom, Evan D. Skillman, Adam Smercina, Estephani E. TorresVillanueva

Abstract: We measure the high-mass stellar initial mass function (IMF) from resolved stars in M33 young stellar clusters. Leveraging \textit{Hubble Space Telescope's} high resolving power, we fully model the IMF probabilistically. We first model the optical CMD of each cluster to constrain its power-law slope $Γ$, marginalized over other cluster parameters in the fit (e.g., cluster age, mass, and radius). W… ▽ More We measure the high-mass stellar initial mass function (IMF) from resolved stars in M33 young stellar clusters. Leveraging \textit{Hubble Space Telescope's} high resolving power, we fully model the IMF probabilistically. We first model the optical CMD of each cluster to constrain its power-law slope $Γ$, marginalized over other cluster parameters in the fit (e.g., cluster age, mass, and radius). We then probabilistically model the distribution of MF slopes for a highly strict cluster sample of 9 clusters more massive than log(Mass/M$_{\odot}$)=3.6; above this mass, all clusters have well-populated main sequences of massive stars and should have accurate recovery of their MF slopes, based on extensive tests with artificial clusters. We find the ensemble IMF is best described by a mean high-mass slope of $\overlineΓ = 1.49\pm0.18$, with an intrinsic scatter of $σ^{2}_Γ = 0.02^{+0.16}_{0.00}$, consistent with a universal IMF. We find no dependence of the IMF on environmental impacts such as the local star formation rate or galactocentric radius within M33, which serves as a proxy for metallicity. This $\overlineΓ$ measurement is consistent with similar measurements in M31, despite M33 having a much higher star formation rate intensity. While this measurement is formally consistent with the canonical Kroupa ($Γ= 1.30$) IMF, as well as the Salpeter ($Γ= 1.35)$) value, it is the second Local Group cluster sample to show evidence for a somewhat steeper high-mass IMF slope. We explore the impacts a steeper IMF slope has on a number of astronomical sub-fields. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: Accepted for publication in ApJ. 9 Figures, 1 Table

Constraints on Ultra Heavy Dark Matter Properties from Dwarf Spheroidal Galaxies with LHAASO Observations

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: In this work we try to search for signals generated by ultra-heavy dark matter at the Large High Altitude Air Shower Observatory (LHAASO) data. We look for possible gamma-ray by dark matter annihilation or decay from 16 dwarf spheroidal galaxies in the field of view of LHAASO. Dwarf spheroidal galaxies are among the most promising targets for indirect detection of dark matter which have low fluxes… ▽ More In this work we try to search for signals generated by ultra-heavy dark matter at the Large High Altitude Air Shower Observatory (LHAASO) data. We look for possible gamma-ray by dark matter annihilation or decay from 16 dwarf spheroidal galaxies in the field of view of LHAASO. Dwarf spheroidal galaxies are among the most promising targets for indirect detection of dark matter which have low fluxes of astrophysical $γ$-ray background while large amount of dark matter. By analyzing more than 700 days observational data at LHAASO, no significant dark matter signal from 1 TeV to 1 EeV is detected. Accordingly we derive the most stringent constraints on the ultra-heavy dark matter annihilation cross-section up to EeV. The constraints on the lifetime of dark matter in decay mode are also derived. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: 17 pages, 12 figures, accepted by PRL

IceCube Search for Neutrino Emission from X-ray Bright Seyfert Galaxies

Authors: R. Abbasi, M. Ackermann, J. Adams, S. K. Agarwalla, J. A. Aguilar, M. Ahlers, J. M. Alameddine, N. M. Amin, K. Andeen, C. Argüelles, Y. Ashida, S. Athanasiadou, L. Ausborm, S. N. Axani, X. Bai, A. Balagopal V., M. Baricevic, S. W. Barwick, S. Bash, V. Basu, R. Bay, J. J. Beatty, J. Becker Tjus, J. Beise, C. Bellenghi , et al. (400 additional authors not shown)

Abstract: The recent IceCube detection of TeV neutrino emission from the nearby active galaxy NGC 1068 suggests that active galactic nuclei (AGN) could make a sizable contribution to the diffuse flux of astrophysical neutrinos. The absence of TeV $γ$-rays from NGC 1068 indicates neutrino production in the vicinity of the supermassive black hole, where the high radiation density leads to $γ$-ray attenuation.… ▽ More The recent IceCube detection of TeV neutrino emission from the nearby active galaxy NGC 1068 suggests that active galactic nuclei (AGN) could make a sizable contribution to the diffuse flux of astrophysical neutrinos. The absence of TeV $γ$-rays from NGC 1068 indicates neutrino production in the vicinity of the supermassive black hole, where the high radiation density leads to $γ$-ray attenuation. Therefore, any potential neutrino emission from similar sources is not expected to correlate with high-energy $γ$-rays. Disk-corona models predict neutrino emission from Seyfert galaxies to correlate with keV X-rays, as they are tracers of coronal activity. Using through-going track events from the Northern Sky recorded by IceCube between 2011 and 2021, we report results from a search for individual and aggregated neutrino signals from 27 additional Seyfert galaxies that are contained in the BAT AGN Spectroscopic Survey (BASS). Besides the generic single power-law, we evaluate the spectra predicted by the disk-corona model. Assuming all sources to be intrinsically similar to NGC 1068, our findings constrain the collective neutrino emission from X-ray bright Seyfert galaxies in the Northern Hemisphere, but, at the same time, show excesses of neutrinos that could be associated with the objects NGC 4151 and CGCG 420-015. These excesses result in a 2.7$σ$ significance with respect to background expectations. △ Less

Submitted 11 June, 2024; originally announced June 2024.

Comments: 17 pages, 9 figures

Search for neutrino emission from hard X-ray AGN with IceCube

Authors: R. Abbasi, M. Ackermann, J. Adams, S. K. Agarwalla, J. A. Aguilar, M. Ahlers, J. M. Alameddine, N. M. Amin, K. Andeen, C. Argüelles, Y. Ashida, S. Athanasiadou, L. Ausborm, S. N. Axani, X. Bai, A. Balagopal V., M. Baricevic, S. W. Barwick, S. Bash, V. Basu, R. Bay, J. J. Beatty, J. Becker Tjus, J. Beise, C. Bellenghi , et al. (401 additional authors not shown)

Abstract: Active Galactic Nuclei (AGN) are promising candidate sources of high-energy astrophysical neutrinos since they provide environments rich in matter and photon targets where cosmic ray interactions may lead to the production of gamma rays and neutrinos. We searched for high-energy neutrino emission from AGN using the $\textit{Swift}$-BAT Spectroscopic Survey (BASS) catalog of hard X-ray sources and… ▽ More Active Galactic Nuclei (AGN) are promising candidate sources of high-energy astrophysical neutrinos since they provide environments rich in matter and photon targets where cosmic ray interactions may lead to the production of gamma rays and neutrinos. We searched for high-energy neutrino emission from AGN using the $\textit{Swift}$-BAT Spectroscopic Survey (BASS) catalog of hard X-ray sources and 12 years of IceCube muon track data. First, upon performing a stacked search, no significant emission was found. Second, we searched for neutrinos from a list of 43 candidate sources and found an excess from the direction of two sources, Seyfert galaxies NGC 1068 and NGC 4151. We observed NGC 1068 at flux $φ_{ν_μ+\barν_μ}$ = $4.02_{-1.52}^{+1.58} \times 10^{-11}$ TeV$^{-1}$ cm$^{-2}$ s$^{-1}$ normalized at 1 TeV, with power-law spectral index, $γ$ = 3.10$^{+0.26}_{-0.22}$, consistent with previous IceCube results. The observation of a neutrino excess from the direction of NGC 4151 is at a post-trial significance of 2.9$σ$. If interpreted as an astrophysical signal, the excess observed from NGC 4151 corresponds to a flux $φ_{ν_μ+\barν_μ}$ = $1.51_{-0.81}^{+0.99} \times 10^{-11}$ TeV$^{-1}$ cm$^{-2}$ s$^{-1}$ normalized at 1 TeV and $γ$ = 2.83$^{+0.35}_{-0.28}$. △ Less

Submitted 12 June, 2024; v1 submitted 10 June, 2024; originally announced June 2024.

M17 MIR: A Massive Star is Forming via Episodic Mass Accretion

Authors: Wei Zhou, Zhiwei Chen, Zhibo Jiang, Haoran Feng, Yu Jiang

Abstract: We analyzed the Atacama Large Millimeter/submillimeter Array (ALMA) band 6 data for the outbursting massive protostar M17~MIR. The ALMA CO $J=2-1$ data reveal a collimated and bipolar north-south outflow from M17~MIR. The blue-shifted outflow exhibits four CO knots (N1 to N4) along the outflow axis, while the red-shifted outflow appears as a single knot (S1). The extremely high velocity (EHV) emis… ▽ More We analyzed the Atacama Large Millimeter/submillimeter Array (ALMA) band 6 data for the outbursting massive protostar M17~MIR. The ALMA CO $J=2-1$ data reveal a collimated and bipolar north-south outflow from M17~MIR. The blue-shifted outflow exhibits four CO knots (N1 to N4) along the outflow axis, while the red-shifted outflow appears as a single knot (S1). The extremely high velocity (EHV) emissions of N1 and S1 are jet-like and contain sub-knots along the outflow axis. Assuming the nearest EHV sub-knots trace the ejecta from the accretion outbursts in the past decades, a tangential ejection velocity of $\sim421\,\mathrm{km\,s^{-1}}$ is derived for M17~MIR. Assuming the same velocity, the dynamical times of the multiple ejecta, traced by the four blue-shifted CO knots, range from 20 to 364 years. The four blue-shifted CO knots imply four clustered accretion outbursts with a duration of tens of years in the past few hundred years. The intervals between the four clustered accretion outbursts are also about tens of years. These properties of the four clustered accretion outbursts are in line with the disk gravitational instability and fragmentation model. The episodic accretion history of M17~MIR traced by episodic outflow suggests that a massive star can form from a lower-mass protostar via frequent episodic accretion events triggered by disk gravitational instability and fragmentation. The first detection of the knotty outflow from an outbursting massive protostar suggests that mass ejections accompanied with accretion events could serve as an effective diagnostic tool for the episodic accretion histories of massive protostars. △ Less

Submitted 17 June, 2024; v1 submitted 7 June, 2024; originally announced June 2024.

Comments: Accepted for publication in ApJL; typos corrected

On the capability of high redshift kSZ measurement with galaxy surveys

Authors: Ziyang Chen, Pengjie Zhang

Abstract: The kSZ effect has been detected at z<1 using various techniques and data sets. The ongoing and upcoming spectroscopic galaxy surveys such as DESI and PFS will push the detection beyond z = 1, and therefore map the baryon distribution at high redshifts. Such detection can be achieved by both the kSZ stacking and tomography methods. While the two methods are theoretically equivalent, they differ si… ▽ More The kSZ effect has been detected at z<1 using various techniques and data sets. The ongoing and upcoming spectroscopic galaxy surveys such as DESI and PFS will push the detection beyond z = 1, and therefore map the baryon distribution at high redshifts. Such detection can be achieved by both the kSZ stacking and tomography methods. While the two methods are theoretically equivalent, they differ significantly in the probed physics and scales, and required data sets. Taking the combination of PFS and ACT as an example, we build mocks of kSZ and galaxies, quantify the kSZ detection S/N, and compare between the two methods. We segment the PFS galaxies into three redshift bins: 0.6 < z < 1.0, 1.0 < z < 1.6, and 1.6 < z < 2.4. For tomography method, our analysis reveals that the two higher redshift bins exhibit higher S/N, with values of 32 and 28, respectively, compared to the first redshift bin (S/N = 8). This is attributed to not only the increasing of electron density with redshifts, but also the larger survey volume and the reduced non-linearity, facilitating velocity reconstruction at higher redshifts. Therefore, the capability of the PFS survey to measure high redshift kSZ effect stands as a substantial advantage over other spectroscopic surveys at lower redshift. The S/N of kSZ stacking largely depends on the number of clusters/groups available from another photometric survey. But in general, its S/N is lower than that of kSZ tomography. Incorporating next-generation CMB surveys like CMB-S4, characterized by significantly reduced instrument noise and improved angular resolution, is expected to enhance tomographic detection by a factor of ten and stacking detection by five. This future high S/N detection holds the promise of not only providing precise constraints on the overall baryon abundance but also initiating a new insight into baryon distribution. △ Less

Submitted 7 June, 2024; originally announced June 2024.

Comments: 23 pages, 6 figures

Imputation of Missing Photometric Data and Photometric Redshift Estimation for CSST

Authors: Zhijian Luo, Zhirui Tang, Zhu Chen, Liping Fu, Wei Du, Shaohua Zhang, Yan Gong, Chenggang Shu, Junhao Lu, Yicheng Li, Xian-Min Meng, Xingchen Zhou, Zuhui Fan

Abstract: Accurate photometric redshift (photo-$z$) estimation requires support from multi-band observational data. However, in the actual process of astronomical observations and data processing, some sources may have missing observational data in certain bands for various reasons. This could greatly affect the accuracy and reliability of photo-$z$ estimation for these sources, and even render some estimat… ▽ More Accurate photometric redshift (photo-$z$) estimation requires support from multi-band observational data. However, in the actual process of astronomical observations and data processing, some sources may have missing observational data in certain bands for various reasons. This could greatly affect the accuracy and reliability of photo-$z$ estimation for these sources, and even render some estimation methods unusable. The same situation may exist for the upcoming Chinese Space Station Telescope (CSST). In this study, we employ a deep learning method called Generative Adversarial Imputation Networks (GAIN) to impute the missing photometric data in CSST, aiming to reduce the impact of data missing on photo-$z$ estimation and improve estimation accuracy. Our results demonstrate that using the GAIN technique can effectively fill in the missing photometric data in CSST. Particularly, when the data missing rate is below 30\%, the imputation of photometric data exhibits high accuracy, with higher accuracy in the $g$, $r$, $i$, $z$, and $y$ bands compared to the $NUV$ and $u$ bands. After filling in the missing values, the quality of photo-$z$ estimation obtained by the widely used Easy and Accurate Zphot from Yale (EAZY) software is notably enhanced. Evaluation metrics for assessing the quality of photo-$z$ estimation, including the catastrophic outlier fraction ($f_{out}$), the normalized median absolute deviation ($\rm {σ_{NMAD}}$), and the bias of photometric redshift ($bias$), all show some degree of improvement. Our research will help maximize the utilization of observational data and provide a new method for handling sample missing values for applications that require complete photometry data to produce results. △ Less

Submitted 3 June, 2024; originally announced June 2024.

Data quality control system and long-term performance monitor of the LHAASO-KM2A

Authors: Zhen Cao, F. Aharonian, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, W. Bian, A. V. Bukevich, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, H. X. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. Chen , et al. (263 additional authors not shown)

Abstract: The KM2A is the largest sub-array of the Large High Altitude Air Shower Observatory (LHAASO). It consists of 5216 electromagnetic particle detectors (EDs) and 1188 muon detectors (MDs). The data recorded by the EDs and MDs are used to reconstruct primary information of cosmic ray and gamma-ray showers. This information is used for physical analysis in gamma-ray astronomy and cosmic ray physics. To… ▽ More The KM2A is the largest sub-array of the Large High Altitude Air Shower Observatory (LHAASO). It consists of 5216 electromagnetic particle detectors (EDs) and 1188 muon detectors (MDs). The data recorded by the EDs and MDs are used to reconstruct primary information of cosmic ray and gamma-ray showers. This information is used for physical analysis in gamma-ray astronomy and cosmic ray physics. To ensure the reliability of the LHAASO-KM2A data, a three-level quality control system has been established. It is used to monitor the status of detector units, stability of reconstructed parameters and the performance of the array based on observations of the Crab Nebula and Moon shadow. This paper will introduce the control system and its application on the LHAASO-KM2A data collected from August 2021 to July 2023. During this period, the pointing and angular resolution of the array were stable. From the observations of the Moon shadow and Crab Nebula, the results achieved using the two methods are consistent with each other. According to the observation of the Crab Nebula at energies from 25 TeV to 100 TeV, the time averaged pointing errors are estimated to be $-0.003^{\circ} \pm 0.005^{\circ}$ and $0.001^{\circ} \pm 0.006^{\circ}$ in the R.A. and Dec directions, respectively. △ Less

Submitted 13 June, 2024; v1 submitted 20 May, 2024; originally announced May 2024.

Comments: 15 pages, 9 figures

Constraining the nonstandard propagating gravitational waves in the cosmological background with GWTC-3

Authors: Zu-Cheng Chen, Lang Liu

Abstract: The detection of gravitational waves (GWs) has opened a new window to test the fundamental nature of gravity. We present constraints on the nonstandard propagation of GWs using the spectral siren method applied to binary black hole (BBH) mergers from the third Gravitational-Wave Transient Catalog (GWTC-3). The spectral siren method exploits the redshift distribution of BBHs to probe the cosmic exp… ▽ More The detection of gravitational waves (GWs) has opened a new window to test the fundamental nature of gravity. We present constraints on the nonstandard propagation of GWs using the spectral siren method applied to binary black hole (BBH) mergers from the third Gravitational-Wave Transient Catalog (GWTC-3). The spectral siren method exploits the redshift distribution of BBHs to probe the cosmic expansion history and break degeneracies between cosmology and modified gravity effects. We focus on the friction term $ν$ in the nonstandard GW propagation equation, which characterizes the running of the Planck mass. Assuming the standard $Λ$CDM cosmology, we find $ν= 0.5^{+3.5}_{-2.6}$ (median and $90\%$ credible interval), improving upon previous constraints from the bright siren event GW170817 by an order of magnitude. This improvement is due to the higher redshifts of BBHs in GWTC-3, reaching up to $z \sim 1$. Our result suggests that the propagation of GWs is consistent with the predictions of general relativity, placing limits on modified gravity theories that predict a time-varying Planck mass. As the sensitivity of GW detectors improves, the spectral siren method will provide a powerful tool for testing gravity on cosmological scales and probing the physics of the early Universe. △ Less

Submitted 20 May, 2024; v1 submitted 16 May, 2024; originally announced May 2024.

Comments: v1, 7 pages, 1 figure; v2, references added

Discovery of Very-high-energy Gamma-ray Emissions from the Low Luminosity AGN NGC 4278 by LHAASO

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: The first source catalog of Large High Altitude Air Shower Observatory reported the detection of a very-high-energy gamma ray source, 1LHAASO J1219+2915. In this paper a further detailed study of the spectral and temporal behavior of this point-like source have been carried. The best-fit position of the TeV source ($\rm{RA}=185.05^{\circ}\pm0.04^{\circ}$, $\rm{Dec}=29.25^{\circ}\pm0.03^{\circ}$) i… ▽ More The first source catalog of Large High Altitude Air Shower Observatory reported the detection of a very-high-energy gamma ray source, 1LHAASO J1219+2915. In this paper a further detailed study of the spectral and temporal behavior of this point-like source have been carried. The best-fit position of the TeV source ($\rm{RA}=185.05^{\circ}\pm0.04^{\circ}$, $\rm{Dec}=29.25^{\circ}\pm0.03^{\circ}$) is compatible with NGC 4278 within $\sim0.03$ degree. Variation analysis shows an indication of the variability at a few months level in the TeV band, which is consistent with low frequency observations. Based on these observations, we report the detection of TeV $γ$-ray emissions from this low-luminosity AGN NGC 4278. The observations by LHAASO-WCDA during active period has a significance level of 8.8\,$σ$ with best-fit photon spectral index $\varGamma=2.56\pm0.14$ and a flux $f_{1-10\,\rm{TeV}}=(7.0\pm1.1_{\rm{sta}}\pm0.35_{\rm{syst}})\times10^{-13}\,\rm{photons\,cm^{-2}\,s^{-1}}$, or approximately $5\%$ of the Crab Nebula. The discovery of VHE from NGC 4278 indicates that the compact, weak radio jet can efficiently accelerate particles and emit TeV photons. △ Less

Submitted 13 May, 2024; originally announced May 2024.

Comments: 11 pages, 5 figures

Power spectrum multipoles and clustering wedges during the Epoch of Reionization

Authors: Zhaoting Chen, Alkistis Pourtsidou

Abstract: We study the viability of using power spectrum clustering wedges as summary statistics of 21cm surveys during the Epoch of Reionization (EoR). For observations along a large lightcone $z\sim 7-9$, the power spectrum is subject to large anisotropic effects due to the evolution along the light-of-sight. Information on the physics of reionization can be extracted from the anisotropy using the power s… ▽ More We study the viability of using power spectrum clustering wedges as summary statistics of 21cm surveys during the Epoch of Reionization (EoR). For observations along a large lightcone $z\sim 7-9$, the power spectrum is subject to large anisotropic effects due to the evolution along the light-of-sight. Information on the physics of reionization can be extracted from the anisotropy using the power spectrum multipoles. Signals of the power spectrum monopole are highly correlated at scales smaller than the typical ionization bubble, which can be disentangled by including higher-order multipoles. By simulating observations of the low frequency part of the SKA Observatory, we find that the sampling of the cylindrical wavenumber $k$-space is highly non-uniform due to the baseline distribution. Measurements in clustering wedges can be used for isolating parts of $k$-space with relatively uniform sampling, allowing for more precise parameter inference. Using Fisher Matrix forecasts, we find that the reionization model can be inferred with per-cent level precision with $\sim 120$hrs of integration time using SKA-Low. Comparing to using only the power spectrum monopole above the foreground wedge, model inference using multipole power spectra in clustering wedges yields a factor of $\sim 3$ improvement on parameter constraints. △ Less

Submitted 8 May, 2024; originally announced May 2024.

Comments: 14 pages, 11 figures. Submitted to MNRAS

Search for joint multimessenger signals from potential Galactic PeVatrons with HAWC and IceCube

Authors: R. Alfaro, C. Alvarez, J. C. Arteaga-Velázquez, D. Avila Rojas, H. A. Ayala Solares, R. Babu, E. Belmont-Moreno, K. S. Caballero-Mora, T. Capistrán, A. Carramiñana, S. Casanova, U. Cotti, J. Cotzomi, S. Coutiño de León, E. De la Fuente, D. Depaoli, N. Di Lalla, R. Diaz Hernandez, J. C. Díaz-Vélez, K. Engel, T. Ergin, K. L. Fan, K. Fang, N. Fraija, S. Fraija , et al. (469 additional authors not shown)

Abstract: Galactic PeVatrons are sources that can accelerate cosmic rays to PeV energies. The high-energy cosmic rays are expected to interact with the surrounding ambient material or radiation, resulting in the production of gamma rays and neutrinos. To optimize for the detection of such associated production of gamma rays and neutrinos for a given source morphology and spectrum, a multi-messenger analysis… ▽ More Galactic PeVatrons are sources that can accelerate cosmic rays to PeV energies. The high-energy cosmic rays are expected to interact with the surrounding ambient material or radiation, resulting in the production of gamma rays and neutrinos. To optimize for the detection of such associated production of gamma rays and neutrinos for a given source morphology and spectrum, a multi-messenger analysis that combines gamma rays and neutrinos is required. In this study, we use the Multi-Mission Maximum Likelihood framework (3ML) with IceCube Maximum Likelihood Analysis software (i3mla) and HAWC Accelerated Likelihood (HAL) to search for a correlation between 22 known gamma-ray sources from the third HAWC gamma-ray catalog and 14 years of IceCube track-like data. No significant neutrino emission from the direction of the HAWC sources was found. We report the best-fit gamma-ray model and 90% CL neutrino flux limit from the 22 sources. From the neutrino flux limit, we conclude that the gamma-ray emission from five of the sources can not be produced purely from hadronic interactions. We report the limit for the fraction of gamma rays produced by hadronic interactions for these five sources. △ Less

Submitted 6 May, 2024; originally announced May 2024.

Acceptance Tests of more than 10 000 Photomultiplier Tubes for the multi-PMT Digital Optical Modules of the IceCube Upgrade

Authors: R. Abbasi, M. Ackermann, J. Adams, S. K. Agarwalla, J. A. Aguilar, M. Ahlers, J. M. Alameddine, N. M. Amin, K. Andeen, C. Argüelles, Y. Ashida, S. Athanasiadou, L. Ausborm, S. N. Axani, X. Bai, A. Balagopal V., M. Baricevic, S. W. Barwick, S. Bash, V. Basu, R. Bay, J. J. Beatty, J. Becker Tjus, J. Beise, C. Bellenghi , et al. (399 additional authors not shown)

Abstract: More than 10,000 photomultiplier tubes (PMTs) with a diameter of 80 mm will be installed in multi-PMT Digital Optical Modules (mDOMs) of the IceCube Upgrade. These have been tested and pre-calibrated at two sites. A throughput of more than 1000 PMTs per week with both sites was achieved with a modular design of the testing facilities and highly automated testing procedures. The testing facilities… ▽ More More than 10,000 photomultiplier tubes (PMTs) with a diameter of 80 mm will be installed in multi-PMT Digital Optical Modules (mDOMs) of the IceCube Upgrade. These have been tested and pre-calibrated at two sites. A throughput of more than 1000 PMTs per week with both sites was achieved with a modular design of the testing facilities and highly automated testing procedures. The testing facilities can easily be adapted to other PMTs, such that they can, e.g., be re-used for testing the PMTs for IceCube-Gen2. Single photoelectron response, high voltage dependence, time resolution, prepulse, late pulse, afterpulse probabilities, and dark rates were measured for each PMT. We describe the design of the testing facilities, the testing procedures, and the results of the acceptance tests. △ Less

Submitted 20 June, 2024; v1 submitted 30 April, 2024; originally announced April 2024.

Comments: 24 pages, 19 figures, 2 tables, submitted to JINST

Revisiting the Constraint on Equation of State of Neutron Star based on the Binary Neutron Star Mergers

Authors: Yun-Peng Li, Zhi-Lin Chen, Da-Bin Lin, En-Wei Liang

Abstract: The merger of neutron star (NS)-NS binary can form different production of the compact remnant, among which the supramassive NS (SMNS) could create an internal plateau and the followed steep decay marks the collapse of the SMNS. The proportion of SMNS and the corresponding collapse-time are often used to constrain the NS equation of state (EoS). This paper revisits this topic by considering the ef… ▽ More The merger of neutron star (NS)-NS binary can form different production of the compact remnant, among which the supramassive NS (SMNS) could create an internal plateau and the followed steep decay marks the collapse of the SMNS. The proportion of SMNS and the corresponding collapse-time are often used to constrain the NS equation of state (EoS). This paper revisits this topic by considering the effect of an accretion disk on the compact remnant, which is not considered in previous works. Compared with previous works, the collapse-time distribution (peaks $\sim$100 s) of the SMNSs formed from NS-NS merger is almost unaffected by the initial surface magnetic ($B_{{\rm s},i}$) of NS, but the total energy output of the magnetic dipole radiation from the SMNSs depends on $B_{{\rm s},i}$ significantly. Coupling the constraints from the SMNS fraction, we exclude some EoSs and obtain three candidate EoSs, i.e., DD2, ENG, and MPA1. By comparing the distributions of the collapse-time and the luminosity of the internal plateau (in the short gamma-ray bursts) for those from observations with those obtained based on the three candidate EoSs, it is shown that only the EoS of ENG is favored. Our sample based on the ENG EOS and a mass distribution motivated by Galactic systems suggests that approximately $99\%$ of NS-NS mergers collapse to form a black hole within $10^7$s. This includes scenarios forming a BH promptly ($36.5\%$), a SMNS ($60.7\%$), or a stable NS that transitions into a BH or a SMNS following accretion ($2.1\%$). It also indicates that the remnants for GW170817 and GW190425, and the second object of GW190814 are more likely to be BHs. △ Less

Submitted 30 April, 2024; originally announced April 2024.

Deep Learning for Cosmological Parameter Inference from Dark Matter Halo Density Field

Authors: Zhiwei Min, Xu Xiao, Jiacheng Ding, Liang Xiao, Jie Jiang, Donglin Wu, Qiufan Lin, Yin Li, Yang Wang, Shuai Liu, Zhixin Chen, Xiangru Li, Jinqu Zhang, Le Zhang, Xiao-Dong Li

Abstract: We propose a lightweight deep convolutional neural network (lCNN) to estimate cosmological parameters from simulated three-dimensional DM halo distributions and associated statistics. The training dataset comprises 2000 realizations of a cubic box with a side length of 1000 $h^{-1}{\rm Mpc}$, and interpolated over a cubic grid of $300^3$ voxels, with each simulation produced using $512^3$ DM parti… ▽ More We propose a lightweight deep convolutional neural network (lCNN) to estimate cosmological parameters from simulated three-dimensional DM halo distributions and associated statistics. The training dataset comprises 2000 realizations of a cubic box with a side length of 1000 $h^{-1}{\rm Mpc}$, and interpolated over a cubic grid of $300^3$ voxels, with each simulation produced using $512^3$ DM particles and $512^3$ neutrinos . Under the flat $Λ$CDM model, simulations vary standard six cosmological parameters including $Ω_m$, $Ω_b$, $h$, $n_s$, $σ_8$, $w$, along with the neutrino mass sum, $M_ν$. We find that: 1) within the framework of lCNN, extracting large-scale structure information is more efficient from the halo density field compared to relying on the statistical quantities including the power spectrum, the two-point correlation function, and the coefficients from wavelet scattering transform; 2) combining the halo density field with its Fourier transformed counterpart enhances predictions, while augmenting the training dataset with measured statistics further improves performance; 3) achieving high accuracy in inferring $Ω_m$, $h$, $n_s$, and $σ_8$ by the neural network model, while being inefficient in predicting $Ω_b$,$M_ν$ and $w$; 4) compared to the simple random forest network trained with three statistical quantities, lCNN yields unbiased estimations with reduced statistical errors: approximately 33.3\% for $Ω_m$, 20.0\% for $h$, 8.3\% for $n_s$, and 40.0\% for $σ_8$. Our study emphasizes this lCNN-based novel approach in extracting large-scale structure information and estimating cosmological parameters. △ Less

Submitted 1 September, 2024; v1 submitted 15 April, 2024; originally announced April 2024.

Comments: v2: matches version accepted for publication in PRD,17 pages,12 figures

Detecting a Gravitational-Wave Background from Inflation with Null Energy Condition Violation: Prospects for Taiji

Authors: Zu-Cheng Chen, Lang Liu

Abstract: The null energy condition (NEC) is a fundamental principle in general relativity, and its violation could leave discernible signatures in gravitational waves (GWs). A violation of the NEC during the primordial era would imprint a blue-tilted spectrum on the stochastic gravitational wave background (SGWB) at nanohertz frequencies, potentially accounting for the recently detected signal by pulsar ti… ▽ More The null energy condition (NEC) is a fundamental principle in general relativity, and its violation could leave discernible signatures in gravitational waves (GWs). A violation of the NEC during the primordial era would imprint a blue-tilted spectrum on the stochastic gravitational wave background (SGWB) at nanohertz frequencies, potentially accounting for the recently detected signal by pulsar timing arrays. Remarkably, models of NEC violation during inflation also predict a nearly scale-invariant GW spectrum in the millihertz frequency range, which could be detectable by upcoming space-based GW detectors such as Taiji. The observation of this distinctive spectrum would provide compelling evidence for new physics beyond the standard cosmological paradigm. In this study, we explore Taiji's ability to detect an SGWB arising from NEC violation during inflation, considering various foregrounds and noise sources, including an extragalactic foreground from binary black hole mergers throughout the universe, a galactic foreground from white dwarf binaries, and the intrinsic noise of the Taiji detector. Employing comprehensive Bayesian parameter estimation techniques to analyze simulated Taiji data, we demonstrate a remarkable precision improvement of three orders of magnitude compared to the NANOGrav 15-year data set for measuring the tensor power spectrum amplitude, $P_{T,2}$, during the second inflationary stage. This substantial enhancement in measurement capabilities underscores Taiji's potential as a powerful probe for investigating the NEC violation in the early Universe. △ Less

Submitted 25 August, 2024; v1 submitted 12 April, 2024; originally announced April 2024.

Comments: 16 pages, 5 figures

Constraints on Inflation with Null Energy Condition Violation from Advanced LIGO and Advanced Virgo's First Three Observing Runs

Authors: Zu-Cheng Chen, Lang Liu

Abstract: The null energy condition (NEC) is a cornerstone of general relativity, and its violation could leave observable imprints in the cosmic gravitational wave spectrum. Theoretical models suggest that NEC violations during inflation can amplify the primordial tensor power spectrum, leading to distinct features in the stochastic gravitational wave background (SGWB). In this work, we search for these NE… ▽ More The null energy condition (NEC) is a cornerstone of general relativity, and its violation could leave observable imprints in the cosmic gravitational wave spectrum. Theoretical models suggest that NEC violations during inflation can amplify the primordial tensor power spectrum, leading to distinct features in the stochastic gravitational wave background (SGWB). In this work, we search for these NEC-violating signatures in the SGWB using data from Advanced LIGO and Advanced Virgo's first three observing runs. Our analysis reveals no statistically significant evidence of such signals, allowing us to place stringent upper limits on the tensor power spectrum amplitude, $P_{T,2}$, during the second inflationary stage. Specifically, we find that $P_{T,2} \lesssim 0.15$ at a $95\%$ confidence level. Notably, this upper limit is consistent with constraints derived from pulsar timing array observations, reinforcing the hypothesis that NEC violations during inflation could explain the signal detected by pulsar timing arrays. Our findings contribute to a deeper understanding of the early Universe and highlight the potential of current and future gravitational wave experiments in probing the physics of inflation and NEC violations. △ Less

Submitted 19 May, 2024; v1 submitted 10 April, 2024; originally announced April 2024.

Comments: 11 pages, 1 figures, title changed as required by the Editor, accepted by JCAP;

Probing the shape of the brown dwarf desert around main-sequence A-F-G-type stars using post-common-envelope WD$-$BD binaries

Authors: Zhangliang Chen, Yizhi Chen, Chen Chen, Hongwei Ge, Bo Ma

Abstract: Brown dwarfs (BDs) possessing masses within the range $40-60 M_{\rm Jup}$ are rare around solar-type main-sequence (MS) stars, which gives rise to the brown dwarf desert (BDD). One caveat associated with previous studies of BDD is the relatively limited sample size of MS$-$BD binaries with accurately determined BD masses. We aim to produce a large sample of brown dwarf companions with precisely de… ▽ More Brown dwarfs (BDs) possessing masses within the range $40-60 M_{\rm Jup}$ are rare around solar-type main-sequence (MS) stars, which gives rise to the brown dwarf desert (BDD). One caveat associated with previous studies of BDD is the relatively limited sample size of MS$-$BD binaries with accurately determined BD masses. We aim to produce a large sample of brown dwarf companions with precisely determined mass around main-sequence A-F-G type stars using observations of post common-envelope white dwarf (WD)$-$BD binaries. We employ the rapid binary evolution code COMPAS to deduce the properties of MS$-$BD binary progenitors from post common-envelope WD$-$BD binaries. This method supplements the directly observed MS$-$BD binary sample, enriching the data available for analyzing BDD around main-sequence A-F-G type stars. Our study opens a new window for studying the shape of BDD around A-F-G type main-sequence stars in the short period regime. We find tentative evidence that the `driest' part of BDD around A-F-G type stars may extend into an orbital period of several hundred days, albeit with a small sample size. More post common-envelope WD$-$BD binaries detected in the future will advance our understanding of the BDD around A-F-G type stars. △ Less

Submitted 10 April, 2024; originally announced April 2024.

Comments: 7 pages, 4 figures. Accepted for publication in A&A

GW230529_181500: A Potential Primordial Binary Black Hole Merger in the Mass Gap

Authors: Qing-Guo Huang, Chen Yuan, Zu-Cheng Chen, Lang Liu

Abstract: During the fourth observing run of the LIGO-Virgo-KAGRA detector network, the LIGO Livingston observatory detected a coalescing compact binary, GW230529_181500, with component masses of $2.5-4.5\, M_\odot$ and $1.2-2.0\, M_\odot$ at the $90\%$ credible level. The gravitational-wave data alone is insufficient to determine whether the components are neutron stars or black holes. In this paper, we pr… ▽ More During the fourth observing run of the LIGO-Virgo-KAGRA detector network, the LIGO Livingston observatory detected a coalescing compact binary, GW230529_181500, with component masses of $2.5-4.5\, M_\odot$ and $1.2-2.0\, M_\odot$ at the $90\%$ credible level. The gravitational-wave data alone is insufficient to determine whether the components are neutron stars or black holes. In this paper, we propose that GW230529_181500 originated from the merger of two primordial black holes (PBHs). We estimate a merger rate of $5.0^{+47.0}_{-4.9} \mathrm{Gpc}^{-3}\,\mathrm{yr}^{-1}$ for compact binary coalescences with properties similar to GW230529_181500. Assuming the source is a PBH-PBH merger, GW230529-like events lead to approximately $1.7^{+36.2}_{-1.5} \times 10^{-3}$ of the dark matter in the form of PBHs. The required abundance of PBHs to explain this event is consistent with existing upper limits derived from microlensing, cosmic microwave background observations and the null detection of gravitational wave background by LIGO-Virgo-KAGRA. △ Less

Submitted 1 August, 2024; v1 submitted 8 April, 2024; originally announced April 2024.

Comments: 15 pages, 3 figures; Accepted by JCAP

LHAASO-KM2A detector simulation using Geant4

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (254 additional authors not shown)

Abstract: KM2A is one of the main sub-arrays of LHAASO, working on gamma ray astronomy and cosmic ray physics at energies above 10 TeV. Detector simulation is the important foundation for estimating detector performance and data analysis. It is a big challenge to simulate the KM2A detector in the framework of Geant4 due to the need to track numerous photons from a large number of detector units (>6000) with… ▽ More KM2A is one of the main sub-arrays of LHAASO, working on gamma ray astronomy and cosmic ray physics at energies above 10 TeV. Detector simulation is the important foundation for estimating detector performance and data analysis. It is a big challenge to simulate the KM2A detector in the framework of Geant4 due to the need to track numerous photons from a large number of detector units (>6000) with large altitude difference (30 m) and huge coverage (1.3 km^2). In this paper, the design of the KM2A simulation code G4KM2A based on Geant4 is introduced. The process of G4KM2A is optimized mainly in memory consumption to avoid memory overffow. Some simpliffcations are used to signiffcantly speed up the execution of G4KM2A. The running time is reduced by at least 30 times compared to full detector simulation. The particle distributions and the core/angle resolution comparison between simulation and experimental data of the full KM2A array are also presented, which show good agreement. △ Less

Submitted 7 April, 2024; originally announced April 2024.

The polarization-angle flip in GRB prompt emission

Authors: Kangfa Cheng, Xiaohong Zhao, Jirong Mao, Zhifu Chen

Abstract: Aims. We aim to provide an explanation for the PA rotation in GRBs and find the physical conditions that lead to the rotation by 90 degrees in the toroidal magnetic-field (MF) model. Moreover, we present some observable polarization properties in the MF model that can be tested in the future. Results. We find that the PA rotation in the toroidal MF is primarily related to three critical factors: t… ▽ More Aims. We aim to provide an explanation for the PA rotation in GRBs and find the physical conditions that lead to the rotation by 90 degrees in the toroidal magnetic-field (MF) model. Moreover, we present some observable polarization properties in the MF model that can be tested in the future. Results. We find that the PA rotation in the toroidal MF is primarily related to three critical factors: the viewing angle, the jet opening angle, and the jet Lorentz factor. Additionally, the PA can experience flips of 90 degrees twice. The conditions for the flips are $q \gtrsim 0.5$ (except for $q\simeq 1$) and $y_j =(Γθ_j)^2 \gtrsim 4$. However, the two flips in the PA might not be concurrently observable due to the constraint of flux. Taking these conditions into account and assuming a random orientation between the jet axis and the line of sight (LOS), we obtain a theoretical upper limit (without any constraints) for the observed rate of GRBs in the X-ray or $γ$-ray band displaying the flips in PA as $R_{ch} \lesssim 80\%$. We further constrain the observed rate as $R_{ch} \sim 16\%$ according to the maximal post-flip polarized flux level, where the observed rate of single and double flips each account for $\sim 8\%$. It should be noted that the observed rates are different in various wavebands. The observed rate of the second PA flip in the optical bands should be higher than that in the X-ray or $γ$-ray band since the flux in the optical band declines much slower than that in the X-ray or $γ$-ray band. Moreover, when the LOS is close to the jet edge ($q\to 1$), it is the easiest case in which to observe the 90-degree PA flip due to the relatively high post-flip polarized flux level. △ Less

Submitted 6 May, 2024; v1 submitted 24 March, 2024; originally announced March 2024.

Comments: 17 pages, 13 figures, Accepted by Astronomy & Astrophysics

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

Strong Coupling of Hydrodynamics and Reactions in Nuclear Statistical Equilibrium for Modeling Convection in Massive Stars

Authors: Michael Zingale, Zhi Chen, Eric T. Johnson, Max P. Katz, Alexander Smith Clark

Abstract: We build on the simplified spectral deferred corrections (SDC) coupling of hydrodynamics and reactions to handle the case of nuclear statistical equilibrium (NSE) and electron/positron captures/decays in the cores of massive stars. Our approach blends a traditional reaction network on the grid with a tabulated NSE state from a very large, $\mathcal{O}(100)$ nuclei, network. We demonstrate how to a… ▽ More We build on the simplified spectral deferred corrections (SDC) coupling of hydrodynamics and reactions to handle the case of nuclear statistical equilibrium (NSE) and electron/positron captures/decays in the cores of massive stars. Our approach blends a traditional reaction network on the grid with a tabulated NSE state from a very large, $\mathcal{O}(100)$ nuclei, network. We demonstrate how to achieve second-order accuracy in the simplified-SDC framework when coupling NSE to hydrodynamics, with the ability to evolve the star on the hydrodynamics timestep. We discuss the application of this method to convection in massive stars leading up to core-collapse. We also show how to initialize the initial convective state from a 1D model in a self-consistent fashion. All of these developments are done in the publicly available Castro simulation code and the entire simulation methodology is fully GPU accelerated. △ Less

Submitted 21 March, 2024; originally announced March 2024.

Comments: submitted to ApJ

Measurements of All-Particle Energy Spectrum and Mean Logarithmic Mass of Cosmic Rays from 0.3 to 30 PeV with LHAASO-KM2A

Authors: The LHAASO Collaboration, Zhen Cao, F. Aharonian, Q. An, A. Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen , et al. (256 additional authors not shown)

Abstract: We present the measurements of all-particle energy spectrum and mean logarithmic mass of cosmic rays in the energy range of 0.3-30 PeV using data collected from LHAASO-KM2A between September 2021 and December 2022, which is based on a nearly composition-independent energy reconstruction method, achieving unprecedented accuracy. Our analysis reveals the position of the knee at… ▽ More We present the measurements of all-particle energy spectrum and mean logarithmic mass of cosmic rays in the energy range of 0.3-30 PeV using data collected from LHAASO-KM2A between September 2021 and December 2022, which is based on a nearly composition-independent energy reconstruction method, achieving unprecedented accuracy. Our analysis reveals the position of the knee at $3.67 \pm 0.05 \pm 0.15$ PeV. Below the knee, the spectral index is found to be -$2.7413 \pm 0.0004 \pm 0.0050$, while above the knee, it is -$3.128 \pm 0.005 \pm 0.027$, with the sharpness of the transition measured with a statistical error of 2%. The mean logarithmic mass of cosmic rays is almost heavier than helium in the whole measured energy range. It decreases from 1.7 at 0.3 PeV to 1.3 at 3 PeV, representing a 24% decline following a power law with an index of -$0.1200 \pm 0.0003 \pm 0.0341$. This is equivalent to an increase in abundance of light components. Above the knee, the mean logarithmic mass exhibits a power law trend towards heavier components, which is reversal to the behavior observed in the all-particle energy spectrum. Additionally, the knee position and the change in power-law index are approximately the same. These findings suggest that the knee observed in the all-particle spectrum corresponds to the knee of the light component, rather than the medium-heavy components. △ Less

Submitted 26 March, 2024; v1 submitted 15 March, 2024; originally announced March 2024.

Comments: 8 pages, 3 figures

Journal ref: Physical Review Letters 132, 131002 (2024)

Relative velocities between $^{13}$CO structures within $^{12}$CO Molecular clouds

Authors: Lixia Yuan, Ji Yang, Xuepeng Chen, Yang Su, Shaobo Zhang, Xin Zhou, Zhiwei Chen, Qing-Zeng Yan, Min Fang, Fujun Du, Yan Sun, Hongchi Wang, Ye Xu

Abstract: Velocity fields of molecular clouds (MCs) can provide crucial information on the merger and split between clouds, as well as their internal kinematics and maintenance, energy injection and redistribution, even star formation within clouds. Using the CO spectral lines data from the Milky Way Imaging Scroll Painting (MWISP) survey, we measure the relative velocities along the line of sight ($Δ$V… ▽ More Velocity fields of molecular clouds (MCs) can provide crucial information on the merger and split between clouds, as well as their internal kinematics and maintenance, energy injection and redistribution, even star formation within clouds. Using the CO spectral lines data from the Milky Way Imaging Scroll Painting (MWISP) survey, we measure the relative velocities along the line of sight ($Δ$V$_{\rm LOS}$) between $^{13}$CO structures within $^{12}$CO MCs. Emphasizing MCs with double and triple $^{13}$CO structures, we find that approximately 70$\%$ of $Δ$V$_{\rm LOS}$ values are less than $\sim$ 1 km s$^{-1}$, and roughly 10$\%$ of values exceed 2 km s$^{-1}$, with a maximum of $\sim$ 5 km s$^{-1}$. Additionally, we compare $Δ$V$_{\rm LOS}$ with the internal velocity dispersion of $^{13}$CO structures ($σ_{\rm ^{13}CO,in}$) and find that about 40$\%$ of samples in either double or triple regime display distinct velocity discontinuities, i.e. the relative velocities between $^{13}$CO structures are larger than the internal linewidths of $^{13}$CO structures. Among these 40$\%$ samples in the triple regime, 33$\%$ exhibit signatures of combinations through the two-body motion, whereas the remaining 7$\%$ show features of configurations through the multiple-body motion. The $Δ$V$_{\rm LOS}$ distributions for MCs with double and triple $^{13}$CO structures are similar, as well as their $Δ$V$_{\rm LOS}$/$σ_{\rm ^{13}CO,in}$ distributions. This suggests that relative motions of $^{13}$CO structures within MCs are random and independent of cloud complexities and scales. △ Less

Submitted 11 March, 2024; originally announced March 2024.

Comments: 16 pages, 11 figures, accepted for publication in AJ

Observation of Seven Astrophysical Tau Neutrino Candidates with IceCube

Authors: IceCube Collaboration, R. Abbasi, M. Ackermann, J. Adams, S. K. Agarwalla, J. A. Aguilar, M. Ahlers, J. M. Alameddine, N. M. Amin, K. Andeen, G. Anton, C. Argüelles, Y. Ashida, S. Athanasiadou, S. N. Axani, X. Bai, A. Balagopal V., M. Baricevic, S. W. Barwick, V. Basu, R. Bay, J. J. Beatty, J. Becker Tjus, J. Beise, C. Bellenghi , et al. (380 additional authors not shown)

Abstract: We report on a measurement of astrophysical tau neutrinos with 9.7 years of IceCube data. Using convolutional neural networks trained on images derived from simulated events, seven candidate $ν_τ$ events were found with visible energies ranging from roughly 20 TeV to 1 PeV and a median expected parent $ν_τ$ energy of about 200 TeV. Considering backgrounds from astrophysical and atmospheric neutrin… ▽ More We report on a measurement of astrophysical tau neutrinos with 9.7 years of IceCube data. Using convolutional neural networks trained on images derived from simulated events, seven candidate $ν_τ$ events were found with visible energies ranging from roughly 20 TeV to 1 PeV and a median expected parent $ν_τ$ energy of about 200 TeV. Considering backgrounds from astrophysical and atmospheric neutrinos, and muons from $π^\pm/K^\pm$ decays in atmospheric air showers, we obtain a total estimated background of about 0.5 events, dominated by non-$ν_τ$ astrophysical neutrinos. Thus, we rule out the absence of astrophysical $ν_τ$ at the $5σ$ level. The measured astrophysical $ν_τ$ flux is consistent with expectations based on previously published IceCube astrophysical neutrino flux measurements and neutrino oscillations. △ Less

Submitted 26 March, 2024; v1 submitted 4 March, 2024; originally announced March 2024.

Comments: Accepted for publication in Physical Review Letters. This version includes full author list metadata

Journal ref: Phys.Rev.Lett. 132 (2024) 15, 151001

Improved modeling of in-ice particle showers for IceCube event reconstruction

Authors: R. Abbasi, M. Ackermann, J. Adams, S. K. Agarwalla, J. A. Aguilar, M. Ahlers, J. M. Alameddine, N. M. Amin, K. Andeen, G. Anton, C. Argüelles, Y. Ashida, S. Athanasiadou, L. Ausborm, S. N. Axani, X. Bai, A. Balagopal V., M. Baricevic, S. W. Barwick, S. Bash, V. Basu, R. Bay, J. J. Beatty, J. Becker Tjus, J. Beise , et al. (394 additional authors not shown)

Abstract: The IceCube Neutrino Observatory relies on an array of photomultiplier tubes to detect Cherenkov light produced by charged particles in the South Pole ice. IceCube data analyses depend on an in-depth characterization of the glacial ice, and on novel approaches in event reconstruction that utilize fast approximations of photoelectron yields. Here, a more accurate model is derived for event reconstr… ▽ More The IceCube Neutrino Observatory relies on an array of photomultiplier tubes to detect Cherenkov light produced by charged particles in the South Pole ice. IceCube data analyses depend on an in-depth characterization of the glacial ice, and on novel approaches in event reconstruction that utilize fast approximations of photoelectron yields. Here, a more accurate model is derived for event reconstruction that better captures our current knowledge of ice optical properties. When evaluated on a Monte Carlo simulation set, the median angular resolution for in-ice particle showers improves by over a factor of three compared to a reconstruction based on a simplified model of the ice. The most substantial improvement is obtained when including effects of birefringence due to the polycrystalline structure of the ice. When evaluated on data classified as particle showers in the high-energy starting events sample, a significantly improved description of the events is observed. △ Less

Submitted 22 April, 2024; v1 submitted 4 March, 2024; originally announced March 2024.

Comments: 28 pages, 18 figures, 1 table, submitted to JINST, updated to account for comments received

Journal ref: 2024 JINST 19 P06026

The Multilayer Nature of Molecular Gas toward the Cygnus Region

Authors: Shiyu Zhang, Yang Su, Xuepeng Chen, Min Fang, Qingzeng Yan, Shaobo Zhang, Yan Sun, Xiaolong Wang, Haoran Feng, Yuehui Ma, Miaomiao Zhang, Zi Zhuang, Xin Zhou, Zhiwei Chen, Ji Yang

Abstract: We study the physical properties and 3D distribution of molecular clouds (MCs) toward the Cygnus region using the MWISP CO survey and Gaia DR3 data. Based on Gaussian decomposition and clustering for $\rm ^{13}CO$ lines, over 70% of the fluxes are recovered. With the identification result of $\rm ^{13}CO$ structures, two models are designed to measure the distances of the molecular gas in velocity… ▽ More We study the physical properties and 3D distribution of molecular clouds (MCs) toward the Cygnus region using the MWISP CO survey and Gaia DR3 data. Based on Gaussian decomposition and clustering for $\rm ^{13}CO$ lines, over 70% of the fluxes are recovered. With the identification result of $\rm ^{13}CO$ structures, two models are designed to measure the distances of the molecular gas in velocity crowding regions. The distances of more than 200 large $\rm ^{13}CO$ structures are obtained toward the 150 square degree region. Additionally, tens of the identified MC structures coincide well with masers and/or intense mid-IR emission. We find multiple gas layers toward the region: (1) the extensive gas structures composing the Cygnus Rift from 700 pc to 1 kpc across the whole region; (2) the $\sim$ 1.3 kpc gas layer mainly in the Cygnus X South region; and (3) the 1.5 kpc dense filament at the Cygnus X North region and many cometary clouds shaped by Cygnus OB2. We also note that the spatial distribution of YSO candidates is generally consistent with the molecular gas structures. The total molecular mass of the Cygnus region is estimated to be $\sim 2.7\times10^{6}~M_{\odot}$ assuming an X-factor ratio $X_{\rm CO} = 2 \times 10^{20} \rm cm^{-2} (K\cdot km\cdot s^{-1})^{-1}$. The foreground Cygnus Rift contributes $\sim$25% of the molecular mass in the whole region. Our work presents a new 3D view of the MCs' distribution toward the Cygnus X region, as well as the exact molecular gas mass distribution in the foreground Cygnus Rift. △ Less

Submitted 23 April, 2024; v1 submitted 29 February, 2024; originally announced March 2024.

Comments: 51 pages, 26 figures, 4 tables, to match the AJ version (2024 AJ 167 220Z). The data can be found at doi: 10.57760/sciencedb.16716

Characterization of the Astrophysical Diffuse Neutrino Flux using Starting Track Events in IceCube

Authors: R. Abbasi, M. Ackermann, J. Adams, S. K. Agarwalla, J. A. Aguilar, M. Ahlers, J. M. Alameddine, N. M. Amin, K. Andeen, G. Anton, C. Argüelles, Y. Ashida, S. Athanasiadou, L. Ausborm, S. N. Axani, X. Bai, A. Balagopal V., M. Baricevic, S. W. Barwick, S. Bash, V. Basu, R. Bay, J. J. Beatty, J. Becker Tjus, J. Beise , et al. (394 additional authors not shown)

Abstract: A measurement of the diffuse astrophysical neutrino spectrum is presented using IceCube data collected from 2011-2022 (10.3 years). We developed novel detection techniques to search for events with a contained vertex and exiting track induced by muon neutrinos undergoing a charged-current interaction. Searching for these starting track events allows us to not only more effectively reject atmospher… ▽ More A measurement of the diffuse astrophysical neutrino spectrum is presented using IceCube data collected from 2011-2022 (10.3 years). We developed novel detection techniques to search for events with a contained vertex and exiting track induced by muon neutrinos undergoing a charged-current interaction. Searching for these starting track events allows us to not only more effectively reject atmospheric muons but also atmospheric neutrino backgrounds in the southern sky, opening a new window to the sub-100 TeV astrophysical neutrino sky. The event selection is constructed using a dynamic starting track veto and machine learning algorithms. We use this data to measure the astrophysical diffuse flux as a single power law flux (SPL) with a best-fit spectral index of $γ= 2.58 ^{+0.10}_{-0.09}$ and per-flavor normalization of $φ^{\mathrm{Astro}}_{\mathrm{per-flavor}} = 1.68 ^{+0.19}_{-0.22} \times 10^{-18} \times \mathrm{GeV}^{-1} \mathrm{cm}^{-2} \mathrm{s}^{-1} \mathrm{sr}^{-1}$ (at 100 TeV). The sensitive energy range for this dataset is 3 - 550 TeV under the SPL assumption. This data was also used to measure the flux under a broken power law, however we did not find any evidence of a low energy cutoff. △ Less

Submitted 27 February, 2024; originally announced February 2024.

Comments: 27 pages, 28 figures

Journal ref: Phys. Rev. D 110, 022001 (2024)

New method for estimating molecular cloud distances based on Gaia, 2MASS, and the TRILEGAL galaxy model

Authors: Juan Mei, Zhiwei Chen, Zhibo Jiang, Sheng Zheng, Haoran Feng

Abstract: We propose a new method for estimating the distances of molecular clouds traced by CO line emission. Stars from 2MASS and Gaia EDR3 are selected as on-cloud stars when they are projected on a cloud. The background on-cloud stars have redder colors on average than the foreground stars. Instead of searching for stars projected away from the cloud, we employed the TRILEGA galaxy model to mimic the st… ▽ More We propose a new method for estimating the distances of molecular clouds traced by CO line emission. Stars from 2MASS and Gaia EDR3 are selected as on-cloud stars when they are projected on a cloud. The background on-cloud stars have redder colors on average than the foreground stars. Instead of searching for stars projected away from the cloud, we employed the TRILEGA galaxy model to mimic the stellar population without cloud extinction along the sightline toward the cloud. Our method does not require an exact boundary of a cloud. The boundaries are highly variable and depend on the sensitivity of the molecular line data. For each cloud, we compared the distributions of on-cloud stars to the TRILEGAL stellar populations in the diagram of $J-K_s$ color versus distance. The intrinsic $J-K_s$ colors of main-sequence and evolved stars from TRILEGAL were considered separately, and they were used as the baseline for subtracting the observed $J-K_s$ colors. The baseline-corrected $J-K_s$ color was deployed with the Bayesian analysis and Markov chain Monte Carlo sampling to determine the distance at which the $J-K_s$ color jump is largest. This method was successfully applied to measure the distances of 27 molecular clouds, which were selected from previously published cloud samples. By replacing TRILEGAL with the GALAXIA galaxy model, we were able to measure the distances for 21 of the 27 clouds. The distances of the 21 clouds based on the GALAXIA model agree well with those based on the TRILEGAL model. The distances of the 27 clouds estimated by this method are consistent with previous estimates. We will apply this new method to a larger region of the gaseous galactic plane, in particular, for the inner galactic region, where a region free of CO emission is hard to separate from the crowded field of clouds. △ Less

Submitted 5 May, 2024; v1 submitted 27 February, 2024; originally announced February 2024.

Comments: 11 pages, 11 figures, accepted by A&A

Journal ref: A&A 685, A39 (2024)

Can we distinguish the adiabatic fluctuations and isocurvature fluctuations with pulsar timing arrays?

Authors: Zu-Cheng Chen, Lang Liu

Abstract: Understanding the nature of primordial fluctuations is critical to our comprehension of the Universe's early stages. While these fluctuations are known to be nearly scale-invariant, quasi-adiabatic, and nearly Gaussian on large scales, their behavior at smaller scales remains less well-defined and may offer insights into new physics. Recent observations by the NANOGrav, PPTA, EPTA, and CPTA collab… ▽ More Understanding the nature of primordial fluctuations is critical to our comprehension of the Universe's early stages. While these fluctuations are known to be nearly scale-invariant, quasi-adiabatic, and nearly Gaussian on large scales, their behavior at smaller scales remains less well-defined and may offer insights into new physics. Recent observations by the NANOGrav, PPTA, EPTA, and CPTA collaborations suggest the presence of a stochastic gravitational wave background, which, while consistent with the contribution from supermassive black hole binaries, also opens the possibility of probing new physics. This paper explores whether this signal could stem from primordial isocurvature and adiabatic fluctuations. We adopt parameterized spectra for both types of fluctuations to fit the observations from the latest NANOGrav data. Furthermore, we employ Bayesian analysis to assess the distinguishability of these models in light of current PTA sensitivities. Our findings indicate that with the capabilities, PTAs cannot conclusively differentiate between isocurvature and adiabatic fluctuations. △ Less

Submitted 26 February, 2024; originally announced February 2024.

Comments: 10 pages

On the initial spin period distribution of neutron stars

Authors: Shen-Shi Du, Xiao-Jin Liu, Zu-Cheng Chen, Zhi-Qiang You, Xing-Jiang Zhu, Zong-Hong Zhu

Abstract: We derive the initial spin period distribution of neutron stars by studying the population of young pulsars associated with supernova remnants. Our hierarchical Bayesian approach accounts for the measurement uncertainties of individual observations and selection effects. Without correcting for selection effects, as done in previous studies, we find that pulsar initial spin periods follow a Weibull… ▽ More We derive the initial spin period distribution of neutron stars by studying the population of young pulsars associated with supernova remnants. Our hierarchical Bayesian approach accounts for the measurement uncertainties of individual observations and selection effects. Without correcting for selection effects, as done in previous studies, we find that pulsar initial spin periods follow a Weibull distribution, peaking at 40~ms, which is favoured against the log-normal distribution with a Bayes factor of 200. The known selection effects in radio pulsar surveys, including pulse broadening and a period-dependent beaming fraction, have been quantitatively investigated. We show that, based on measurements of pulsar luminosity and spin period from the ATNF Pulsar Catalogue, the impact of pulse broadening on the inference of pulsar period distribution is likely to be insignificant. Correcting for the beaming selection effect, a Weibull distribution remains to be the preferred model, while its peak slightly shifts to longer periods at 50~ms. Our method will prove useful in constraining the birth properties of neutron stars in the Square Kilometre Array era. △ Less

Submitted 29 April, 2024; v1 submitted 20 February, 2024; originally announced February 2024.

Comments: Accepted for publication in The Astrophysical Journal

Journal ref: APJ, 2024, 968:105

Ly$α$ emission in galaxies at $z\simeq5-6$: new insight from JWST into the statistical distributions of Ly$α$ properties at the end of reionization

Authors: Mengtao Tang, Daniel P. Stark, Richard S. Ellis, Fengwu Sun, Michael Topping, Brant Robertson, Sandro Tacchella, Santiago Arribas, William M. Baker, Rachana Bhatawdekar, Kristan Boyett, Andrew J. Bunker, Stéphane Charlot, Zuyi Chen, Jacopo Chevallard, Gareth C. Jones, Nimisha Kumari, Jianwei Lyu, Roberto Maiolino, Michael V. Maseda, Aayush Saxena, Lily Whitler, Christina C. Williams, Chris Willott, Joris Witstok

Abstract: JWST has recently sparked a new era of Ly$α$ spectroscopy, delivering the first measurements of the Ly$α$ escape fraction and velocity profile in typical galaxies at $z\simeq6-10$. These observations offer new prospects for insight into the earliest stages of reionization. But to realize this potential, we need robust models of Ly$α$ properties in galaxies at $z\simeq5-6$ when the IGM is mostly io… ▽ More JWST has recently sparked a new era of Ly$α$ spectroscopy, delivering the first measurements of the Ly$α$ escape fraction and velocity profile in typical galaxies at $z\simeq6-10$. These observations offer new prospects for insight into the earliest stages of reionization. But to realize this potential, we need robust models of Ly$α$ properties in galaxies at $z\simeq5-6$ when the IGM is mostly ionized. Here we use new JWST observations from the JADES and FRESCO surveys combined with VLT/MUSE and Keck/DEIMOS data to characterize statistical distributions of Ly$α$ velocity offsets, escape fractions, and EWs in $z\simeq5-6$ galaxies. We find that galaxies with large Ly$α$ escape fractions (> 0.2) are common at $z\simeq5-6$, comprising 30 per cent of Lyman break selected samples. Comparing to literature studies, our census suggests that Ly$α$ becomes more prevalent in the galaxy population toward higher redshift from $z\sim3$ to $z\sim6$, although we find that this evolution slows considerably between $z\sim5$ and $z\sim6$, consistent with modest attenuation from residual HI in the mostly ionized IGM at $z\simeq5-6$. We find significant evolution in Ly$α$ velocity profiles between $z\simeq2-3$ and $z\simeq5-6$, likely reflecting the influence of resonant scattering from residual intergalactic HI on the escape of Ly$α$ emission near line center. This effect will make it challenging to use Ly$α$ peak offsets as a probe of Lyman continuum leakage at $z\simeq5-6$. We use our $z\simeq5-6$ Ly$α$ distributions to make predictions for typical Ly$α$ properties at $z\gtrsim8$ and discuss implications of a recently-discovered Ly$α$ emitter at $z\simeq8.5$ with a small peak velocity offset (156 km s$^{-1}$). △ Less

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

Comments: 31 pages, 12 figures, accepted by MNRAS

Bridging the gap between luminous red novae and common envelope evolution: the role of recombination energy and radiation force

Authors: Zhuo Chen, Natalia Ivanova

Abstract: Luminous red novae (LRNe) and their connection to common envelope evolution (CEE) remain elusive in astrophysics. Here, we present a radiation hydrodynamic model capable of simulating the light curves of material ejected during a CEE. For the first time, the radiation hydrodynamic model incorporates complete recombination physics for hydrogen and helium. The radiation hydrodynamic equations are so… ▽ More Luminous red novae (LRNe) and their connection to common envelope evolution (CEE) remain elusive in astrophysics. Here, we present a radiation hydrodynamic model capable of simulating the light curves of material ejected during a CEE. For the first time, the radiation hydrodynamic model incorporates complete recombination physics for hydrogen and helium. The radiation hydrodynamic equations are solved with Guangqi. With time-independent ejecta simulations, we show that the peaks in the light curves are attributed to radiation-dominated ejecta, while the extended plateaus are produced by matter-dominated ejecta. To showcase our model's capability, we fit the light curve of AT2019zhd. The central mass object of $6M_{\odot}$ is assumed based on observations and scaling relations. Our model demonstrates that the ejecta mass of AT2019zhd falls within the range of $0.04M_{\odot}$ to $0.1M_{\odot}$. Additionally, we demonstrate that recombination energy and radiation force acceleration significantly impact the light curves, whereas dust formation has a limited effect during the peak and plateau phases. △ Less

Submitted 17 February, 2024; v1 submitted 8 February, 2024; originally announced February 2024.

Comments: accepted by the ApJL

The Stellar "Snake" -- II: The Mass Function

Authors: Xiang-Ming Yang, Sarah A. Bird, Jiadong Li, Hai-Jun Tian, Dan Qiu, Jia-Peng Li, Cheng-Yuan Li, Gao-Chao Liu, Peng Zhang, Ju-Yong Zhang, Zhi-Ping Chen

Abstract: We present a comprehensive investigation on the mass function (MF) of a snake-like stellar structure in the solar neighbourhood, building on our previous discovery. To ensure the reliability of the data, we reselect the member stars of the Stellar ``Snake'' in the latest {\it Gaia} Data Release 3 using the same approach as the initial series of articles. We also precisely measure the physical para… ▽ More We present a comprehensive investigation on the mass function (MF) of a snake-like stellar structure in the solar neighbourhood, building on our previous discovery. To ensure the reliability of the data, we reselect the member stars of the Stellar ``Snake'' in the latest {\it Gaia} Data Release 3 using the same approach as the initial series of articles. We also precisely measure the physical parameters of the clusters within the Stellar Snake. In light of the high completeness of the member stars in the cluster regions, we develop a simulated model color-magnitude diagram-based inference method to derive the mass function, binary fraction, and mass-ratio distribution of the clusters in the Stellar Snake. Notably, despite their similar ages and metallicity, we discover systematic variations in the MFs along the elongation direction of the Snake in the mass range of 0.5 to 2.0 M$_\odot$. The ``head'' of the Snake conforms to a canonical initial mass function with a power-law slope of $α\sim-2.3$. Extending towards the ``tail,'' the MF becomes more top-light, indicating a deficiency of massive stars within these clusters. This result provides evidence for the delayed formation of massive stars in the clusters. Such clues give support to the hypothesis that the Stellar Snake constitutes as a hierarchically primordial structure. △ Less

Submitted 6 February, 2024; originally announced February 2024.

Comments: 19 pages, 15 figures, Accepted for publication on MNRAS

Confronting primordial black holes with LIGO-Virgo-KAGRA and the Einstein Telescope

Authors: Zu-Cheng Chen, Alex Hall

Abstract: The detection of gravitational waves (GWs) from binary black hole (BBH) coalescences by the LIGO-Virgo-KAGRA (LVK) Collaboration has raised fundamental questions about the genesis of these events. In this chapter, we explore the possibility that PBHs, proposed candidates for dark matter, may serve as the progenitors of the BBHs observed by LVK. Employing a Bayesian analysis, we constrain the PBH m… ▽ More The detection of gravitational waves (GWs) from binary black hole (BBH) coalescences by the LIGO-Virgo-KAGRA (LVK) Collaboration has raised fundamental questions about the genesis of these events. In this chapter, we explore the possibility that PBHs, proposed candidates for dark matter, may serve as the progenitors of the BBHs observed by LVK. Employing a Bayesian analysis, we constrain the PBH model using the LVK third GW Transient Catalog (GWTC-3), revealing that stellar-mass PBHs cannot dominate cold dark matter. Considering a mixed population of astrophysical black holes (ABHs) and PBHs, we determine that approximately $1/4$ of the detectable events in the GWTC-3 can be attributed to PBH binaries. We also forecast detectable event rate distributions for PBH and ABH binaries by the third-generation ground-based GW detectors, such as the Einstein Telescope, offering a potential avenue to distinguish PBHs from ABHs based on their distinct redshift evolutions. △ Less

Submitted 6 February, 2024; originally announced February 2024.

Comments: to appear [with updates] in the book "Primordial Black Holes", ed. Chris Byrnes, Gabriele Franciolini, Tomohiro Harada, Paolo Pani, Misao Sasaki; Springer (2024); comments are welcome