Showing 1–50 of 1,178 results for author: LI, Y

Flow-Based Generative Emulation of Grids of Stellar Evolutionary Models

Authors: Marc Hon, Yaguang Li, Joel Ong

Abstract: We present a flow-based generative approach to emulate grids of stellar evolutionary models. By interpreting the input parameters and output properties of these models as multi-dimensional probability distributions, we train conditional normalizing flows to learn and predict the complex relationships between grid inputs and outputs in the form of conditional joint distributions. Leveraging the exp… ▽ More We present a flow-based generative approach to emulate grids of stellar evolutionary models. By interpreting the input parameters and output properties of these models as multi-dimensional probability distributions, we train conditional normalizing flows to learn and predict the complex relationships between grid inputs and outputs in the form of conditional joint distributions. Leveraging the expressive power and versatility of these flows, we showcase their ability to emulate a variety of evolutionary tracks and isochrones across a continuous range of input parameters. In addition, we describe a simple Bayesian approach for estimating stellar parameters using these flows and demonstrate its application to asteroseismic datasets of red giants observed by the Kepler mission. By applying this approach to red giants in open clusters NGC 6791 and NGC 6819, we illustrate how large age uncertainties can arise when fitting only to global asteroseismic and spectroscopic parameters without prior information on initial helium abundances and mixing length parameter values. We also conduct inference using the flow at a large scale by determining revised estimates of masses and radii for 15,388 field red giants. These estimates show improved agreement with results from existing grid-based modelling, reveal distinct population-level features in the red clump, and suggest that the masses of Kepler red giants previously determined using the corrected asteroseismic scaling relations have been overestimated by 5-10%. △ Less

Submitted 12 July, 2024; originally announced July 2024.

Comments: 27 pages, 18 figures. Accepted for publication in ApJ. Code, animation, and interactive visualizations are available at https://github.com/mtyhon/modelflows/. Table 4 is also available as ancillary file attached to this submission

An impulsive geomagnetic effect from an early-impulsive flare

Authors: Hugh S. Hudson, Edward. W. Cliver, Lyndsay Fletcher, Declan A. Diver, Peter T. Gallagher, Ying Li, Christopher M. J. Osborne, Craig Stark, Yang Su

Abstract: The geomagnetic "solar flare effect" (SFE) results from excess ionization in the Earth's ionosphere, famously first detected at the time of the Carrington flare in 1859. This indirect detection of a flare constituted one of the first cases of "multimessenger astronomy," whereby solar ionizing radiation stimulates ionospheric currents. Well-observed SFEs have few-minute time scales and perturbation… ▽ More The geomagnetic "solar flare effect" (SFE) results from excess ionization in the Earth's ionosphere, famously first detected at the time of the Carrington flare in 1859. This indirect detection of a flare constituted one of the first cases of "multimessenger astronomy," whereby solar ionizing radiation stimulates ionospheric currents. Well-observed SFEs have few-minute time scales and perturbations of >10 nT, with the greatest events reaching above 100 nT. In previously reported cases the SFE time profiles tend to resemble those of solar soft X-ray emission, which ionizes the D-region; there is also a less-well-studied contribution from Lyman-alpha. We report here a specific case, from flare SOL2024-03-10 (M7.4), in which an impulsive SFE deviated from this pattern. This flare contained an "early impulsive" component of exceptionally hard radiation, extending up to gamma-ray energies above 1 MeV, distinctly before the bulk of the flare soft X-ray emission. We can characterize the spectral distribution of this early-impulsive component in detail, thanks to the modern extensive wavelength coverage. A more typical gradual SFE occurred during the flare's main phase. We suggest that events of this type warrant exploration of the solar physics in the "impulse response" limit of very short time scales. △ Less

Submitted 12 July, 2024; originally announced July 2024.

Comments: MNRAS to be published

X-ray spectral and timing evolution during the 2018 outburst of MAXI J1820+070

Authors: YaXing Li, Zhen Yan, ChenXu Gao, Wenfei Yu

Abstract: We made use high-cadence observations from the $Insight$-HXMT and $NICER$ to scrutinize the spectral and timing evolution during the 2018 outburst of the black hole X-ray binary (BHXRB) MAXI J1820+070. It's hardness-intensity diagram (HID) displays a ''q''-like track including all the spectral states, along a unique loop in the hard state. The tracks observed in the HID is anticipated in the evolu… ▽ More We made use high-cadence observations from the $Insight$-HXMT and $NICER$ to scrutinize the spectral and timing evolution during the 2018 outburst of the black hole X-ray binary (BHXRB) MAXI J1820+070. It's hardness-intensity diagram (HID) displays a ''q''-like track including all the spectral states, along a unique loop in the hard state. The tracks observed in the HID is anticipated in the evolution of the components responsible for Compton and reflection emission. This is substantiated by the relationship between the X-ray luminosity $L_\mathrm{X}$ and photon index $Γ$, as well as the relationship between X-ray luminosity $L_\mathrm{X}$ and the ratio of Compton to disk luminosities $L_\mathrm{C}/L_\mathrm{D}$. Both of these relationships exhibit a pattern reminiscent of HID. During the hard state, the hardness (also $Γ$) is determined by either reflection component ($R_{f}>1$ ) or Compton component ($R_{f}<1$) depending on the value reflection fraction $R_{f}$. So the distinctive evolution of $R_{f}$ leads to the unique loop in the HID (also in the $L_\mathrm{X}$--$Γ$ plane) of hard state. Additionally, we found a negative correlation between frequency of the type-C quasi-periodic oscillation (QPO) ($ν_{\mathrm{C,QPO}}$) and the optical depth of the Compton emission ($τ$), and a positive correlation between $ν_{\mathrm{C,QPO}}$ and $Γ$. These correlations strongly suggest a coupling between the QPO properties and the underlying process responsible for Comptonization. △ Less

Submitted 11 July, 2024; originally announced July 2024.

Comments: 14 pages, 10 figures, submitted to MNRAS

Spectroastrometry and Reverberation Mapping (SARM) of Active Galactic Nuclei. I. The H$β$ Broad-line Region Structure and Black Hole Mass of Five Quasars

Authors: Yan-Rong Li, Chen Hu, Zhu-Heng Yao, Yong-Jie Chen, Hua-Rui Bai, Sen Yang, Pu Du, Feng-Na Fang, Yi-Xin Fu, Jun-Rong Liu, Yue-Chang Peng, Yu-Yang Songsheng, Yi-Lin Wang, Ming Xiao, Shuo Zhai, Hartmut Winkler, Jin-Ming Bai, Luis C. Ho, Romain G. Petrov, Jesus Aceituno, Jian-Min Wang

Abstract: We conduct a reverberation mapping (RM) campaign to spectroscopically monitor a sample of selected bright active galactic nuclei with large anticipated broad-line region (BLR) sizes adequate for spectroastrometric observations by the GRAVITY instrument on the Very Large Telescope Interferometer. We report the first results for five objects, IC 4329A, Mrk 335, Mrk 509, Mrk 1239, and PDS 456, among… ▽ More We conduct a reverberation mapping (RM) campaign to spectroscopically monitor a sample of selected bright active galactic nuclei with large anticipated broad-line region (BLR) sizes adequate for spectroastrometric observations by the GRAVITY instrument on the Very Large Telescope Interferometer. We report the first results for five objects, IC 4329A, Mrk 335, Mrk 509, Mrk 1239, and PDS 456, among which Mrk 1239 and PDS 456 are for the first time spectroscopically monitored. We obtain multi-year monitoring data and perform multi-component spectral decomposition to extract the broad H$β$ profiles. We detect significant time lags between the H$β$ and continuum variations, generally obeying the previously established BLR size-luminosity relation. Velocity-resolved H$β$ time lags illustrate diverse, possibly evolving BLR kinematics. We further measure the H$β$ line widths from mean and rms spectra and the resulting virial products show good consistency among different seasons. Adopting a unity virial factor and the full width at half maximum of the broad H$β$ line from the mean spectrum as the measure of velocity, the obtained black hole mass averaged over seasons is $\log M_\bullet/M_\odot=8.02_{-0.14}^{+0.09}$, $6.92_{-0.12}^{+0.12}$, $8.01_{-0.25}^{+0.16}$, $7.44_{-0.14}^{+0.13}$, and $8.59_{-0.11}^{+0.07}$ for the five objects, respectively. The black hole mass estimations using other line width measures are also reported (up to the virial factors). For objects with previous RM campaigns, our mass estimates are in agreement with earlier results. In a companion paper, we will employ BLR dynamical modeling to directly infer the black hole mass and thereby determine the virial factors. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Comments: 32 pages, 6 tables, 20 figures. To appear in ApJ

Using Galaxy Evolution as Source of Physics-Based Ground Truth for Generative Models

Authors: Yun Qi Li, Tuan Do, Evan Jones, Bernie Boscoe, Kevin Alfaro, Zooey Nguyen

Abstract: Generative models producing images have enormous potential to advance discoveries across scientific fields and require metrics capable of quantifying the high dimensional output. We propose that astrophysics data, such as galaxy images, can test generative models with additional physics-motivated ground truths in addition to human judgment. For example, galaxies in the Universe form and change ove… ▽ More Generative models producing images have enormous potential to advance discoveries across scientific fields and require metrics capable of quantifying the high dimensional output. We propose that astrophysics data, such as galaxy images, can test generative models with additional physics-motivated ground truths in addition to human judgment. For example, galaxies in the Universe form and change over billions of years, following physical laws and relationships that are both easy to characterize and difficult to encode in generative models. We build a conditional denoising diffusion probabilistic model (DDPM) and a conditional variational autoencoder (CVAE) and test their ability to generate realistic galaxies conditioned on their redshifts (galaxy ages). This is one of the first studies to probe these generative models using physically motivated metrics. We find that both models produce comparable realistic galaxies based on human evaluation, but our physics-based metrics are better able to discern the strengths and weaknesses of the generative models. Overall, the DDPM model performs better than the CVAE on the majority of the physics-based metrics. Ultimately, if we can show that generative models can learn the physics of galaxy evolution, they have the potential to unlock new astrophysical discoveries. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: 20 pages, 14 figures, 1 Table, code: https://github.com/astrodatalab/li2024_public, training data: https://zenodo.org/records/11117528

Unraveling the role of merger histories in the population of Insitu stars: linking IllustrisTNG cosmological simulation to H3 survey

Authors: Razieh Emami, Lars Hernquist, Randall Smith, James F. Steiner, Grant Tremblay, Douglas Finkbeiner, Mark Vogelsberger, Josh Grindlay, Federico Marinacci, Kung-Yi Su, Cecilia Garraffo, Yuan-Sen Ting, Phillip A. Cargile, Rebecca L. Davies, Chloë E. Benton, Yijia Li, Letizia Bugiani, Amir H. Khoram, Sownak Bose

Abstract: We undertake a comprehensive investigation into the distribution of insitu stars within Milky Way-like galaxies, leveraging TNG50 simulations and comparing their predictions with data from the H3 survey. Our analysis reveals that 28% of galaxies demonstrate reasonable agreement with H3, while only 12% exhibit excellent alignment in their profiles, regardless of the specific spatial cut employed to… ▽ More We undertake a comprehensive investigation into the distribution of insitu stars within Milky Way-like galaxies, leveraging TNG50 simulations and comparing their predictions with data from the H3 survey. Our analysis reveals that 28% of galaxies demonstrate reasonable agreement with H3, while only 12% exhibit excellent alignment in their profiles, regardless of the specific spatial cut employed to define insitu stars. To uncover the underlying factors contributing to deviations between TNG50 and H3 distributions, we scrutinize correlation coefficients among internal drivers(e.g., virial radius, star formation rate [SFR]) and merger-related parameters (such as the effective mass-ratio, mean distance, average redshift, total number of mergers, average spin-ratio and maximum spin alignment between merging galaxies). Notably, we identify significant correlations between deviations from observational data and key parameters such as the median slope of virial radius, mean SFR values, and the rate of SFR change across different redshift scans. Furthermore, positive correlations emerge between deviations from observational data and parameters related to galaxy mergers. We validate these correlations using the Random Forest Regression method. Our findings underscore the invaluable insights provided by the H3 survey in unraveling the cosmic history of galaxies akin to the Milky Way, thereby advancing our understanding of galactic evolution and shedding light on the formation and evolution of Milky Way-like galaxies in cosmological simulations. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: 20 pages, 13 figures

Search for Classical Cepheids in Galactic Open Clusters and Calibration of the Period Wesenheit Metallicity Relation in the Gaia Bands

Authors: Huajian Wang, Ye Xu, Zehao Lin, Chaojie Hao, Dejian Liu, Yingjie Li

Abstract: It is beneficial to calibrate the period Wesenheit metallicity relation (PWZR) of Delta Cephei stars (DCEPs), i.e., classical Cepheids, using accurate parallaxes of associated open clusters (OCs) from Gaia data release 3 (DR3). To this aim, we obtain a total of 43 OC-DCEPs (including 33 fundamental mode, 9 first overtone mode, and 1 multimode DCEPs.) and calibrate the PWZR as… ▽ More It is beneficial to calibrate the period Wesenheit metallicity relation (PWZR) of Delta Cephei stars (DCEPs), i.e., classical Cepheids, using accurate parallaxes of associated open clusters (OCs) from Gaia data release 3 (DR3). To this aim, we obtain a total of 43 OC-DCEPs (including 33 fundamental mode, 9 first overtone mode, and 1 multimode DCEPs.) and calibrate the PWZR as $W_G=(-3.356 \,\pm\, 0.033) \,(\log{P-1})+(-5.947 \,\pm\, 0.025)+(-0.285 \,\pm\, 0.064)[\textrm{Fe/H}]$. The concurrently obtained residual parallax offset in OC, $zp = -4\pm5\,μ\textrm{as}$, demonstrate the adequacy of the parallax corrections within the magnitude range of OC member stars. By comparing the field DCEPs' DR3 parallaxes with their photometric parallaxes derived by our PWZR, we estimated the residual parallax offset in field DCEPs as $zp = -15\pm3\,μ\textrm{as}$. Using our PWZR, we estimate the distance modulus of the Large Magellanic Cloud to be $18.482 \,\pm\, 0.040$ mag, which aligns well with the most accurate published value obtained through geometric methods. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Journal ref: 2024 AJ 168 34

Spatial distribution of C4H and c-C3H2 in cold molecular cores

Authors: Yijia Liu, Junzhi Wang, Shu Liu, Ningyu Tang, Yan Gong, Yuqiang Li, Juan LI, Rui Luo, Yani Xu

Abstract: C$_4$H and $c$-C$_3$H$_2$, as unsaturated hydrocarbon molecules, are important for forming large organic molecules in the interstellar medium. We present mapping observations of C$_4$H ($N$=9$-8$) lines, $c$-C$_3$H$_2$ ($J_{Ka,Kb}$=2$_{1,2}$-1$_{0,1}$) %at 85338.894 MHz and H$^{13}$CO$^+$ ($J$=1$-0$) %at 86754.2884 MHz toward 19 nearby cold molecular cores in the Milky Way with the IRAM 30m telesc… ▽ More C$_4$H and $c$-C$_3$H$_2$, as unsaturated hydrocarbon molecules, are important for forming large organic molecules in the interstellar medium. We present mapping observations of C$_4$H ($N$=9$-8$) lines, $c$-C$_3$H$_2$ ($J_{Ka,Kb}$=2$_{1,2}$-1$_{0,1}$) %at 85338.894 MHz and H$^{13}$CO$^+$ ($J$=1$-0$) %at 86754.2884 MHz toward 19 nearby cold molecular cores in the Milky Way with the IRAM 30m telescope. C$_4$H 9--8 was detected in 13 sources, while $c$-C$_3$H$_2$ was detected in 18 sources. The widely existing C$_4$H and $c$-C$_3$H$_2$ molecules in cold cores provide material to form large organic molecules. Different spatial distributions between C$_4$H 9--8 and $c$-C$_3$H$_2$ 2--1 were found. The relative abundances of these three molecules were obtained under the assumption of local thermodynamic equilibrium conditions with a fixed excitation temperature. The abundance ratio of C$_4$H to $c$-C$_3$H$_2$ ranged from 0.34 $\pm$ 0.09 in G032.93+02 to 4.65 $\pm$ 0.50 in G008.67+22. A weak correlation between C$_4$H/H$^{13}$CO$^+$ and $c$-C$_3$H$_2$/H$^{13}$CO$^+$ abundance ratios was found, with a correlation coefficient of 0.46, which indicates that there is no tight astrochemical connection between C$_4$H and $c$-C$_3$H$_2$ molecules. △ Less

Submitted 28 June, 2024; originally announced June 2024.

Comments: 17 pages, 2 figures

A lensed FRB candidate in the first CHIME/FRB Catalogue and its potential implications

Authors: Chenming Chang, Songbo Zhang, Di Xiao, Zhenfan Tang, Ye Li, Junjie Wei, Xuefeng Wu

Abstract: Fast radio bursts (FRBs) are immensely energetic radio pulses with durations of milliseconds. Given their high all-sky rate, the probability of an FRB being lensed by an intervening massive object is non-negligible. In this study, we search for possible lensing candidates within the first Canadian Hydrogen Intensity Mapping Experiment FRB catalogue using an autocorrelation algorithm and verificati… ▽ More Fast radio bursts (FRBs) are immensely energetic radio pulses with durations of milliseconds. Given their high all-sky rate, the probability of an FRB being lensed by an intervening massive object is non-negligible. In this study, we search for possible lensing candidates within the first Canadian Hydrogen Intensity Mapping Experiment FRB catalogue using an autocorrelation algorithm and verification through signal simulations. We identify FRB 20190308C as a lensed candidate with a significance of 3.4 sigma. Furthermore, we constrain the mass of the lensing object using the Chang-Refsdal lens model, based on the flux ratio and time delay between the substructures of FRB 20190308C. Future long-term and high-precision observations are expected to reveal more lensed FRBs. △ Less

Submitted 28 June, 2024; originally announced June 2024.

Comments: 5 pages, 4 figures, submitted

Detection, sky localization and early warning for binary neutron star mergers by detectors located in China of different configurations in third generation detector network

Authors: Yufeng Li, Ik Siong Heng, Man Leong Chan, Xilong Fan, Lijun Gou

Abstract: This work shows the results of an evaluation of the impact that a detector located in China, with a noise budget comparable to that of a proposed high-frequency detector with a 20 km arm length, an Einstein Telescope (ET) or a Cosmic Explorer (CE), could have on the network of ET-CE in terms of detection rate, localization, and providing early warning alert for simulated binary neutron star (BNS)s… ▽ More This work shows the results of an evaluation of the impact that a detector located in China, with a noise budget comparable to that of a proposed high-frequency detector with a 20 km arm length, an Einstein Telescope (ET) or a Cosmic Explorer (CE), could have on the network of ET-CE in terms of detection rate, localization, and providing early warning alert for simulated binary neutron star (BNS)s. The results indicate that a three-detector network including a Chinese detector could identify at least 4.4% more BNS mergers than an ET-CE network alone. The localization uncertainty could be reduced by a factor of more than 5 on average compared to the ET-CE network. With a three-detector network involving a Chinese detector, up to 89% of BNS mergers could be located within 10 square degrees of the sky 10 minutes prior to the merger. The assessment suggests that the potential for early warning signals is highest when the Chinese detector is similar to ET, whereas the sources are detected with the highest signal-to-noise ratio and localized to the smallest regions when the detector is more akin to CE. Interestingly, the C20N network (comprising ET+CE+C20) can achieve comparable localization performance as the ET network while outperforming the ETCN network (featuring the ET+CE+ an ET-like detector in China) in terms of detection capabilities, especially at large distances, indicating that adding a 20 km kilohertz detector in China to ET-CE network would make significant contributions at least as adding an ET-like detector in China to multi-messenger astronomy for almost all BNS observations. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: accepted by PRD

Timing and Scintillation Studies of Pulsars in Globular Cluster M3 (NGC 5272) with FAST

Authors: Baoda Li, Li-yun Zhang, Jumei Yao, Dejiang Yin, Ralph P. Eatough, Minghui Li, Yifeng Li, Yujie Lian, Yu Pan, Yinfeng Dai, Yaowei Li, Xingnan Zhang, Tianhao Su, Yuxiao Wu, Tong Liu, Kuo Liu, Lin Wang, Lei Qian, Zhichen Pan

Abstract: We present the phase-connected timing solutions of all the five pulsars in globular cluster (GC) M3 (NGC 5272), namely PSRs M3A to F (PSRs J1342+2822A to F), with the exception of PSR M3C, from FAST archival data. In these timing solutions, those of PSRs M3E, and F are obtained for the first time. We find that PSRs M3E and F have low mass companions, and are in circular orbits with periods of 7.1… ▽ More We present the phase-connected timing solutions of all the five pulsars in globular cluster (GC) M3 (NGC 5272), namely PSRs M3A to F (PSRs J1342+2822A to F), with the exception of PSR M3C, from FAST archival data. In these timing solutions, those of PSRs M3E, and F are obtained for the first time. We find that PSRs M3E and F have low mass companions, and are in circular orbits with periods of 7.1 and 3.0 days, respectively. For PSR M3C, we have not detected it in all the 41 observations. We found no X-ray counterparts for these pulsars in archival Chandra images in the band of 0.2-20 keV. We noticed that the pulsars in M3 seem to be native. From the Auto-Correlation Function (ACF) analysis of the M3A's and M3B's dynamic spectra, the scintillation timescale ranges from $7.0\pm0.3$ min to $60.0\pm0.6$ min, and the scintillation bandwidth ranges from $4.6\pm0.2$ MHz to $57.1\pm1.1$ MHz. The measured scintillation bandwidths from the dynamic spectra indicate strong scintillation, and the scattering medium is anisotropic. From the secondary spectra, we captured a scintillation arc only for PSR M3B with a curvature of $649\pm23 {\rm m}^{-1} {\rm mHz}^{-2}$. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: 14 pages, 4 figures, accepted for publication in The Astrophysical Journal

Searching anomalies using nonlinear dimensionality reduction techniques

Authors: X. Yang, G. Hobbs, S. -B. Zhang, A. Zic, Lawrence Toomey, Y. Li, J. -S. Wang, S. Dai, X. -F. Wu

Abstract: We have searched for anomalous events using 2,520 hours of archival observations from Murriyang, CSIRO's Parkes radio telescope. These observations were originally undertaken to search for pulsars. We used a machine-learning algorithm based on ResNet and Uniform Manifold Approximation and Projection (UMAP) in order to identify parts of the data stream that potentially contain anomalous signals. Ma… ▽ More We have searched for anomalous events using 2,520 hours of archival observations from Murriyang, CSIRO's Parkes radio telescope. These observations were originally undertaken to search for pulsars. We used a machine-learning algorithm based on ResNet and Uniform Manifold Approximation and Projection (UMAP) in order to identify parts of the data stream that potentially contain anomalous signals. Many of these anomalous events are radio frequency interference, which were subsequently filtered using multibeam information. We detected 202 anomalous events and provide their positions and event times. We discuss the possibility that one of the events comes from radio emission from a white dwarf star.The other events are currently of unknown type. △ Less

Submitted 1 July, 2024; v1 submitted 26 June, 2024; originally announced June 2024.

Comments: 9 pages, 7 figures

CCAT: Detector Noise Limited Performance of the RFSoC-based Readout Electronics for mm/sub-mm/far-IR KIDs

Authors: Adrian K. Sinclair, James Burgoyne, Anthony I. Huber, Colin Murphy, Steve K. Choi, Cody J. Duell, Zachary B. Huber, Yaqiong Li, Scott C. Chapman, Michael D. Niemack, Thomas Nikola, Eve M. Vavagiakis, Samantha Walker, Jordan D. Wheeler, Jason Austermann, Lawrence Lin, Ruixuan Xie, Bugao Zou, Philip D. Mauskopf

Abstract: The Fred Young Submillimeter Telescope (FYST), on Cerro Chajnantor in the Atacama desert of Chile, will conduct wide-field and small deep-field surveys of the sky with more than 100,000 detectors on the Prime-Cam instrument. Kinetic inductance detectors (KIDs) were chosen as the primary sensor technology for their high density focal plane packing. Additionally, they benefit from low cost, ease of… ▽ More The Fred Young Submillimeter Telescope (FYST), on Cerro Chajnantor in the Atacama desert of Chile, will conduct wide-field and small deep-field surveys of the sky with more than 100,000 detectors on the Prime-Cam instrument. Kinetic inductance detectors (KIDs) were chosen as the primary sensor technology for their high density focal plane packing. Additionally, they benefit from low cost, ease of fabrication, and simplified cryogenic readout, which are all beneficial for successful deployment at scale. The cryogenic multiplexing complexity is pulled out of the cryostat and is instead pushed into the digital signal processing of the room temperature electronics. Using the Xilinx Radio Frequency System on a Chip (RFSoC), a highly multiplexed KID readout was developed for the first light Prime-Cam and commissioning Mod-Cam instruments. We report on the performance of the RFSoC-based readout with multiple detector arrays in various cryogenic setups. Specifically we demonstrate detector noise limited performance of the RFSoC-based readout under the expected optical loading conditions. △ Less

Submitted 21 June, 2024; originally announced June 2024.

Comments: draft submitted to SPIE

Multi-messenger modeling of the Monogem pulsar halo

Authors: Youyou Li, Oscar Macias, Shinichiro Ando, Jacco Vink

Abstract: The High-Altitude Water Cherenkov Telescope (HAWC) has detected TeV halos associated with two nearby pulsars/pulsar wind nebulae (PWN) -- Geminga and B0656+14. These TeV halos extend up to tens of pc from the central accelerators, indicating that the diffusion of ultrarelativistic electrons and positrons in the interstellar medium has been suppressed by two orders of magnitude. Although Geminga an… ▽ More The High-Altitude Water Cherenkov Telescope (HAWC) has detected TeV halos associated with two nearby pulsars/pulsar wind nebulae (PWN) -- Geminga and B0656+14. These TeV halos extend up to tens of pc from the central accelerators, indicating that the diffusion of ultrarelativistic electrons and positrons in the interstellar medium has been suppressed by two orders of magnitude. Although Geminga and B0656+14 are at similar distances and in the same field of view, they have distinct histories. Notably, B0656+14 probably still resides within its parent supernova remnant, the Monogem Ring, which can be observed in X-rays. In this work, we perform high-resolution simulations of the propagation and emission of relativistic lepton pairs around B0656+14 using a two-zone diffusion model using the GALPROP numerical code. We compared the predicted inverse-Compton spectrum to the observations made by HAWC and Fermi-LAT and found physically plausible model parameters that resulted in a good fit to the data. Additionally, we estimated the contribution of this TeV-halo to the positron flux observed on Earth and found it to be smaller than 10\% of the measured flux. We conclude that future observations of the TeV halo and its synchrotron emission counterpart in radio and X-ray frequencies will be crucial to distinguish between various possible models. △ Less

Submitted 19 June, 2024; originally announced June 2024.

Concurrent Accretion and Migration of Giant Planets in their Natal Disks with Consistent Accretion Torque

Authors: Ya-Ping Li, Yi-Xian Chen, Douglas N. C. Lin

Abstract: Migration commonly occurs during the epoch of planet formation. For emerging gas giant planets, it proceeds concurrently with their growth through the accretion of gas from their natal protoplanetary disks. Similar migration process should also be applied to the stellar-mass black holes embedded in active galactic nucleus disks. In this work, we perform high resolution 3D and 2D numerical hydrodyn… ▽ More Migration commonly occurs during the epoch of planet formation. For emerging gas giant planets, it proceeds concurrently with their growth through the accretion of gas from their natal protoplanetary disks. Similar migration process should also be applied to the stellar-mass black holes embedded in active galactic nucleus disks. In this work, we perform high resolution 3D and 2D numerical hydrodynamical simulations to study the migration dynamics for accreting embedded objects over the disk viscous timescales in a self-consistent manner. We find that an accreting planet embedded in a predominantly viscous disk has a tendency to migrate outward, in contrast to the inward orbital decay of non-accreting planets. 3D and 2D simulations find the consistent outward migration results for the accreting planets. Under this circumstance, the accreting planet's outward migration is mainly due to the asymmetric spiral arms feeding from the global disk into the Hill radius. This is analogous to the unsaturated corotation torque although the imbalance is due to material accretion within the libration timescale rather than diffusion onto the inner disk. In a disk with a relatively small viscosity, the accreting planets clear deep gaps near their orbits. The tendency of inward migration is recovered, albeit with suppressed rates. By performing a parameter survey with a range of disks' viscosity, we find that the transition from outward to inward migration occurs with the effective viscous efficiency factor $α\sim 0.003$ for Jupiter-mass planets. △ Less

Submitted 18 June, 2024; originally announced June 2024.

Comments: 17 pages, 9 figures, Accepted by ApJ

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

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

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

Submitted 17 June, 2024; originally announced June 2024.

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

Association between a Failed Prominence Eruption and the Drainage of Mass from Another Prominence

Authors: Jianchao Xue, Li Feng, Hui Li, Ping Zhang, Jun Chen, Guanglu Shi, Kaifan Ji, Ye Qiu, Chuan Li, Lei Lu, Beili Ying, Ying Li, Yu Huang, Youping Li, Jingwei Li, Jie Zhao, Dechao Song, Shuting Li, Zhengyuan Tian, Yingna Su, Qingmin Zhang, Yunyi Ge, Jiahui Shan, Qiao Li, Gen Li , et al. (9 additional authors not shown)

Abstract: Sympathetic eruptions of solar prominences have been studied for decades, however, it is usually difficult to identify their causal links. Here we present two failed prominence eruptions on 26 October 2022 and explore their connections. Using stereoscopic observations, the south prominence (PRO-S) erupts with untwisting motions, flare ribbons occur underneath, and new connections are formed during… ▽ More Sympathetic eruptions of solar prominences have been studied for decades, however, it is usually difficult to identify their causal links. Here we present two failed prominence eruptions on 26 October 2022 and explore their connections. Using stereoscopic observations, the south prominence (PRO-S) erupts with untwisting motions, flare ribbons occur underneath, and new connections are formed during the eruption. The north prominence (PRO-N) rises up along with PRO-S, and its upper part disappears due to catastrophic mass draining along an elongated structure after PRO-S failed eruption. We suggest that the eruption of PRO-S initiates due to a kink instability, further rises up, and fails to erupt due to reconnection with surrounding fields. The elongated structure connecting PRO-N overlies PRO-S, which causes the rising up of PRO-N along with PRO-S and mass drainage after PRO-S eruption. This study suggests that a prominence may end its life through mass drainage forced by an eruption underneath. △ Less

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

Comments: 15 pages, 7 figures, has been accepted by Solar Physics

Dense Outflowing Molecular Gas in Massive Star-forming Regions

Authors: Yani Xu, Junzhi Wang, Shu Liu, Juan Li, Yuqiang LI, Rui Luo, Chao Ou, Siqi Zheng, Yijia Liu

Abstract: Dense outflowing gas, traced by transitions of molecules with large dipole moment, is important for understanding mass loss and feedback of massive star formation. HCN 3-2 and HCO$^+$ 3-2 are good tracers of dense outflowing molecular gas, which are closely related to active star formation. In this study, we present on-the-fly (OTF) mapping observations of HCN 3-2 and HCO$^+$ 3-2 toward a sample o… ▽ More Dense outflowing gas, traced by transitions of molecules with large dipole moment, is important for understanding mass loss and feedback of massive star formation. HCN 3-2 and HCO$^+$ 3-2 are good tracers of dense outflowing molecular gas, which are closely related to active star formation. In this study, we present on-the-fly (OTF) mapping observations of HCN 3-2 and HCO$^+$ 3-2 toward a sample of 33 massive star-forming regions using the 10-m Submillimeter Telescope (SMT). With the spatial distribution of line wings of HCO$^+$ 3-2 and HCN 3-2, outflows are detected in 25 sources, resulting in a detection rate of 76$\%$. The optically thin H$^{13}$CN and H$^{13}$CO$^+$ 3-2 lines are used to identify line wings as outflows and estimate core mass. The mass $M_{out}$, momentum $P_{out}$, kinetic energy $E_{K}$, force $F_{out}$ and mass loss rate $\dot M_{out}$ of outflow and core mass, are obtained for each source. A sublinear tight correlation is found between the mass of dense molecular outflow and core mass, with an index of $\sim$ 0.8 and a correlation coefficient of 0.88. △ Less

Submitted 13 June, 2024; originally announced June 2024.

Comments: 22 pages, 5 figures, 4 tables, accepted in AJ

Report number: ad47c4

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

AGN Feedback in Quiescent Galaxies at Cosmic Noon Traced by Ionized Gas Emission

Authors: Letizia Bugiani, Sirio Belli, Minjung Park, Rebecca L. Davies, J. Trevor Mendel, Benjamin D. Johnson, Amir H. Khoram, Chloë Benton, Andrea Cimatti, Charlie Conroy, Razieh Emami, Joel Leja, Yijia Li, Gabriel Maheson, Elijah P. Mathews, Rohan P. Naidu, Erica J. Nelson, Sandro Tacchella, Bryan A. Terrazas, Rainer Weinberger

Abstract: We analyze ionized gas emission lines in deep rest-frame optical spectra of 16 quiescent galaxies at redshift $1.7<z<3.5$ observed with JWST/NIRSpec by the Blue Jay survey. Robust detection of emission lines in $75\%$ of the sample indicates the presence of ongoing ionizing sources in this passive population. The H$α$ line luminosities confirm that the population is quiescent, with star formation… ▽ More We analyze ionized gas emission lines in deep rest-frame optical spectra of 16 quiescent galaxies at redshift $1.7<z<3.5$ observed with JWST/NIRSpec by the Blue Jay survey. Robust detection of emission lines in $75\%$ of the sample indicates the presence of ongoing ionizing sources in this passive population. The H$α$ line luminosities confirm that the population is quiescent, with star formation rates that are at least ten times lower than the main sequence of star formation. The quiescent sample is clearly separate from the star-forming population in line diagnostic diagrams, and occupies a region usually populated by active galactic nuclei (AGN). Analysis of the observed line ratios, equivalent widths, and velocity dispersions leads us to conclude that in most cases the gas is ionized by AGN activity, despite the lack of X-ray detections. A subset of the sample also hosts ionized and/or neutral outflows. Our results show, for the first time using a representative sample, that low luminosity AGN are extremely common among quiescent galaxies at high redshift. These low luminosity AGN may play a key role in quenching star formation and in maintaining massive galaxies quiescent from Cosmic Noon to $z\sim0$. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: 20 pages, 18 figures. Submitted to ApJ

CCAT: Comparisons of 280 GHz TiN and Al Kinetic Inductance Detector Arrays

Authors: Cody J. Duell, Jason Austermann, James Beall, James R. Burgoyne, Scott C. Chapman, Steve K. Choi, Rodrigo G. Freundt, Jiansong Gao, Christopher Groppi, Anthony I. Huber, Zachary B. Huber, Johannes Hubmayr, Ben Keller, Yaqiong Li, Lawrence T. Lin, Justin Matthewson, Philip Mauskopf, Alicia Middleton, Colin C. Murphy, Michael D. Niemack, Thomas Nikola, Adrian K. Sinclair, Ema Smith, Jeff van Lanen, Anna Vaskuri , et al. (5 additional authors not shown)

Abstract: The CCAT Collaboration's six-meter Fred Young Submillimeter Telescope is scheduled to begin observing in the Chilean Atacama in 2025, targeting a variety of science goals throughout cosmic history. Prime-Cam is a 1.8-meter diameter cryostat that will host up to seven independent instrument modules designed for simultaneous spectroscopic and broadband, polarimetric surveys at millimeter to submilli… ▽ More The CCAT Collaboration's six-meter Fred Young Submillimeter Telescope is scheduled to begin observing in the Chilean Atacama in 2025, targeting a variety of science goals throughout cosmic history. Prime-Cam is a 1.8-meter diameter cryostat that will host up to seven independent instrument modules designed for simultaneous spectroscopic and broadband, polarimetric surveys at millimeter to submillimeter wavelengths. The first of these instrument modules, the 280 GHz module, will include ${\sim}$10,000 kinetic inductance detectors (KIDs) across three arrays. While the first array was fabricated out of tri-layer TiN/Ti/TiN, the other two arrays were fabricated out of a single layer of Al. This combination of materials within the same instrument provides a unique opportunity to directly compare the performance and noise properties of two different detector materials that are seeing increasing use within the field. We present preliminary comparisons here based on lab testing, along with a discussion of the potential impacts on operation when observing and translating raw data to science-grade maps. △ Less

Submitted 10 June, 2024; originally announced June 2024.

Comments: 6 pages, 3 figures, conference proceedings submitted to the Journal of Low Temperature Physics

Prevalence of non-standard collapsing of strong Langmuir turbulence in solar corona plasmas

Authors: Yaokun Li, Haomin Sun, Hao Ning, Sulan Ni, Xiangliang Kong, Jiansen He, Yao Chen

Abstract: We present a fully-kinetic simulation of the full life cycle of strong Langmuir turbulence (SLT) excited by electron beams that are accelerated under the solar corona conditions. We find that (1) most packets ($\sim$80%) are affected by their neighbors during their collapse, as a result, their spatial scale variations present non-standard evolutionary features, i.e., deviating away from what was p… ▽ More We present a fully-kinetic simulation of the full life cycle of strong Langmuir turbulence (SLT) excited by electron beams that are accelerated under the solar corona conditions. We find that (1) most packets ($\sim$80%) are affected by their neighbors during their collapse, as a result, their spatial scale variations present non-standard evolutionary features, i.e., deviating away from what was predicted by the Zakharov model; (2) the collapsing cavity is too shallow to trap the wave packet due to the growth of the Coulomb force, as a result a majority ($\sim$70%) of the packet energy runs away and a secondary localization may occur. The study indicates that the non-standard Langmuir collapse may play an important role in coronal plasmas interacting with an intense electron beam, that may be eventually confirmed by humanity's first mission to fly through the corona. △ Less

Submitted 8 June, 2024; originally announced June 2024.

The Heavier the Faster: A Sub-population of Heavy, Rapidly Spinning and Quickly Evolving Binary Black Holes

Authors: Wei-Hua Guo, Yin-Jie Li, Yuan-Zhu Wang, Yong Shao, Shichao Wu, Tao Zhu, Yi-Zhong Fan

Abstract: The spins of binary black holes (BBHs) measured from gravitational waves carry notable information of the formation pathways. Here we propose a quantity "dimensionless net spin" ($χ_{\rm N}$), which is related to the sum of angular momentum of component black holes in the system, to provide a novel perspective to study the origin(s) of BBHs. By performing hierarchical Bayesian inference on… ▽ More The spins of binary black holes (BBHs) measured from gravitational waves carry notable information of the formation pathways. Here we propose a quantity "dimensionless net spin" ($χ_{\rm N}$), which is related to the sum of angular momentum of component black holes in the system, to provide a novel perspective to study the origin(s) of BBHs. By performing hierarchical Bayesian inference on $χ_{\rm N}$, we find strong evidence that the marginal distribution of this quantity can be better fitted by two Gaussian components rather than one: there is a narrow peak at $χ_{\rm N} \sim 0.15$ and another extended peak at $χ_{\rm N} \sim 0.47$. We also find that the rapidly spinning systems likely dominate the high-mass end of the population and they evolve with redshift much quicker. These findings bring new challenges to the field binary scenario, and suggest that dynamical process should plays a key role in forming high total mass BBHs. △ Less

Submitted 5 June, 2024; originally announced June 2024.

Comments: Submitted on 25 April 2024

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.

FAST Discovery of Eight Isolated Millisecond Pulsars in NGC 6517

Authors: Dejiang Yin, Li-yun Zhang, Lei Qian, Ralph P. Eatough, Baoda Li, Duncan R. Lorimer, Yinfeng Dai, Yaowei Li, Xingnan Zhang, Minghui Li, Tianhao Su, Yuxiao Wu, Yu Pan, Yujie Lian, Tong Liu, Zhen Yan, Zhichen Pan

Abstract: We present the discovery of 8 isolated millisecond pulsars in Globular Cluster (GC) NGC 6517 using the Five-Hundred-meter Aperture Spherical radio Telescope (FAST). The spin periods of those pulsars (namely PSR J1801-0857K to R, or, NGC 6517K to R) are all shorter than 10 ms. With these discoveries, NGC 6517 is currently the GC with the most known pulsars in the FAST sky. The largest difference in… ▽ More We present the discovery of 8 isolated millisecond pulsars in Globular Cluster (GC) NGC 6517 using the Five-Hundred-meter Aperture Spherical radio Telescope (FAST). The spin periods of those pulsars (namely PSR J1801-0857K to R, or, NGC 6517K to R) are all shorter than 10 ms. With these discoveries, NGC 6517 is currently the GC with the most known pulsars in the FAST sky. The largest difference in dispersion measure of the pulsars in NGC 6517 is 11.2 cm$^{-3}$ pc, the second among all GCs. The fraction of isolated pulsars in this GC (16 of 17, 94$\%$) is consistent with previous studies indicating an overabundance of isolated pulsars in the densest GCs, especially in those undergoing cluster core collapse. Considering the FAST GC pulsar discoveries, we modeled the GC pulsar population using the empirical Bayesian method described by Turk and Lorimer with the recent counts. Using this approach, we find that the expected number of potential pulsars in GCs seems to be correlated with the central escape velocity, hence, the GCs Liller 1, NGC 6441, M54 (NGC 6715), and $ω$-Cen (NGC 5139) are expected to host the largest numbers of pulsars. △ Less

Submitted 28 May, 2024; originally announced May 2024.

Comments: 21 pages, 2 figures, accepted for publication in The Astrophysical Journal Letters

Five parameters are all you need (in $Λ$CDM)

Authors: Paulo Montero-Camacho, Yin Li, Miles Cranmer

Abstract: The standard cosmological model, with its six independent parameters, has proven remarkably successful in describing the evolution of the Universe. One of these parameters, the optical depth to reionization $τ_{\mathrm{reio}}$, represents the scatterings that Cosmic Microwave Background (CMB) photons will experience after decoupling from the primordial plasma as the intergalactic medium transition… ▽ More The standard cosmological model, with its six independent parameters, has proven remarkably successful in describing the evolution of the Universe. One of these parameters, the optical depth to reionization $τ_{\mathrm{reio}}$, represents the scatterings that Cosmic Microwave Background (CMB) photons will experience after decoupling from the primordial plasma as the intergalactic medium transitions from neutral to ionized. $τ_{\mathrm{reio}}$ depends on the neutral hydrogen fraction $x_{\mathrm{HI}}$, which, in turn, should theoretically depend on cosmology. We present a novel method to establish the missing link between cosmology and reionization timeline. We discover the timeline has a universal shape well described by the Gompertz mortality law, applicable to any cosmology within our simulated data. This enables us to map cosmology to reionization using symbolic regression and to treat $τ_{\mathrm{reio}}$ as a derived parameter. Reanalyzing CMB data with our universal $x_{\mathrm{HI}}$ tightens the constraint on $τ_{\mathrm {reio}}$ by more than one order of magnitude to $\approx 1\%$ and reduced the error on the amplitude of the primordial fluctuations by a factor of 2.5 compared to Planck's PR3 constraint. While our results rely on the astrophysical assumptions in our simulations, the methodology itself is independent of these assumptions; after all, $τ_{\mathrm{reio}}$ is fundamentally a function of cosmology. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: 9 figures. Comments welcomed

A Hamiltonian, post-Born, three-dimensional, on-the-fly ray tracing algorithm for gravitational lensing

Authors: Alan Junzhe Zhou, Yin Li, Scott Dodelson, Rachel Mandelbaum, Yucheng Zhang, Xiangchong Li, Giulio Fabbian

Abstract: The analyses of the next generation cosmological surveys demand an accurate, efficient, and differentiable method for simulating the universe and its observables across cosmological volumes. We present Hamiltonian ray tracing (HRT) -- the first post-Born (accounting for lens-lens coupling and without relying on the Born approximation), three-dimensional (without assuming the thin-lens approximatio… ▽ More The analyses of the next generation cosmological surveys demand an accurate, efficient, and differentiable method for simulating the universe and its observables across cosmological volumes. We present Hamiltonian ray tracing (HRT) -- the first post-Born (accounting for lens-lens coupling and without relying on the Born approximation), three-dimensional (without assuming the thin-lens approximation), and on-the-fly (applicable to any structure formation simulations) ray tracing algorithm based on the Hamiltonian formalism. HRT performs symplectic integration of the photon geodesics in a weak gravitational field, and can integrate tightly with any gravity solver, enabling co-evolution of matter particles and light rays with minimal additional computations. We implement HRT in the particle-mesh library \texttt{pmwd}, leveraging hardware accelerators such as GPUs and automatic differentiation capabilities based on \texttt{JAX}. When tested on a point-mass lens, HRT achieves sub-percent accuracy in deflection angles above the resolution limit across both weak and moderately strong lensing regimes. We also test HRT in cosmological simulations on the convergence maps and their power spectra. △ Less

Submitted 21 May, 2024; originally announced May 2024.

Comments: 25 pages, 8 figures

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

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

Cue: A Fast and Flexible Photoionization Emulator for Modeling Nebular Emission Powered By Almost Any Ionizing Source

Authors: Yijia Li, Joel Leja, Benjamin D. Johnson, Sandro Tacchella, Rebecca Davies, Sirio Belli, Minjung Park, Razieh Emami

Abstract: The complex physics governing nebular emission in galaxies, particularly in the early universe, often defy simple low-dimensional models. This has proven to be a significant barrier in understanding the (often diverse) ionizing sources powering this emission. We present Cue, a highly flexible tool for interpreting nebular emission across a wide range of abundances and ionizing conditions of galaxi… ▽ More The complex physics governing nebular emission in galaxies, particularly in the early universe, often defy simple low-dimensional models. This has proven to be a significant barrier in understanding the (often diverse) ionizing sources powering this emission. We present Cue, a highly flexible tool for interpreting nebular emission across a wide range of abundances and ionizing conditions of galaxies at different redshifts. Unlike typical nebular models used to interpret extragalactic nebular emission, our model does not require a specific ionizing spectrum as a source, instead approximating the ionizing spectrum with a 4-part piece-wise power-law. We train a neural net emulator based on the CLOUDY photoionization modeling code and make self-consistent nebular continuum and line emission predictions. Along with the flexible ionizing spectra, we allow freedom in [O/H], [N/O], [C/O], gas density, and total ionizing photon budget. This flexibility allows us to either marginalize over or directly measure the incident ionizing radiation, thereby directly interrogating the source of the ionizing photons in distant galaxies via their nebular emission. Our emulator demonstrates a high accuracy, with $\sim$1% uncertainty in predicting the nebular continuum and $\sim$5% uncertainty in the emission lines. Mock tests suggest Cue is well-calibrated and produces useful constraints on the ionizing spectra when $S/N (\mathrm{H}_α) \gtrsim 10$, and furthermore capable of distinguishing between the ionizing spectra predicted by single and binary stellar models. The compute efficiency of neural networks facilitates future applications of Cue for rapid modeling of the nebular emission in large samples and Monte Carlo sampling techniques. △ Less

Submitted 7 May, 2024; originally announced May 2024.

Comments: submitted to ApJ

CURLING - I. The Influence of Point-like Image Approximation on the Outcomes of Cluster Strong Lens Modeling

Authors: Yushan Xie, Huanyuan Shan, Nan Li, Ran Li, Eric Jullo, Chen Su, Xiaoyue Cao, Jean-Paul Kneib, Ana Acebron, Mengfan He, Ji Yao, Chunxiang Wang, Jiadong Li, Yin Li

Abstract: Cluster-scale strong lensing is a powerful tool for exploring the properties of dark matter and constraining cosmological models. However, due to the complex parameter space, pixelized strong lens modeling in galaxy clusters is computationally expensive, leading to the point-source approximation of strongly lensed extended images, potentially introducing systematic biases. Herein, as the first pap… ▽ More Cluster-scale strong lensing is a powerful tool for exploring the properties of dark matter and constraining cosmological models. However, due to the complex parameter space, pixelized strong lens modeling in galaxy clusters is computationally expensive, leading to the point-source approximation of strongly lensed extended images, potentially introducing systematic biases. Herein, as the first paper of the ClUsteR strong Lens modelIng for the Next-Generation observations (CURLING) program, we use lensing ray-tracing simulations to quantify the biases and uncertainties arising from the point-like image approximation for JWST-like observations. Our results indicate that the approximation works well for reconstructing the total cluster mass distribution, but can bias the magnification measurements near critical curves and the constraints on the cosmological parameters, the total matter density of the Universe $Ω_{\rm m}$, and dark energy equation of state parameter $w$. To mitigate the biases, we propose incorporating the extended surface brightness distribution of lensed sources into the modeling. This approach reduces the bias in magnification from 46.2 per cent to 0.09 per cent for $μ\sim 1000$. Furthermore, the median values of cosmological parameters align more closely with the fiducial model. In addition to the improved accuracy, we also demonstrate that the constraining power can be substantially enhanced. In conclusion, it is necessary to model cluster-scale strong lenses with pixelized multiple images, especially for estimating the intrinsic luminosity of highly magnified sources and accurate cosmography in the era of high-precision observations. △ Less

Submitted 5 May, 2024; originally announced May 2024.

Comments: 12 pages, 8 figures

Spectral and Imaging Observations of a C2.3 White-Light Flare from the Advanced Space-Based Solar Observatory (ASO-S) and the Chinese H$α$ Solar Explorer (CHASE)

Authors: Qiao Li, Ying Li, Yang Su, Dechao Song, Hui Li, Li Feng, Yu Huang, Youping Li, Jingwei Li, Jie Zhao, Lei Lu, Beili Ying, Jianchao Xue, Ping Zhang, Jun Tian, Xiaofeng Liu, Gen Li, Zhichen Jing, Shuting Li, Guanglu Shi, Zhengyuan Tian, Wei Chen, Yingna Su, Qingmin Zhang, Dong Li , et al. (5 additional authors not shown)

Abstract: Solar white-light flares are characterized by an enhancement in the optical continuum, which are usually large flares (say X- and M-class flares). Here we report a small C2.3 white-light flare (SOL2022-12-20T04:10) observed by the \emph{Advanced Space-based Solar Observatory} and the \emph{Chinese H$α$ Solar Explorer}. This flare exhibits an increase of $\approx$6.4\% in the photospheric Fe \texts… ▽ More Solar white-light flares are characterized by an enhancement in the optical continuum, which are usually large flares (say X- and M-class flares). Here we report a small C2.3 white-light flare (SOL2022-12-20T04:10) observed by the \emph{Advanced Space-based Solar Observatory} and the \emph{Chinese H$α$ Solar Explorer}. This flare exhibits an increase of $\approx$6.4\% in the photospheric Fe \textsc{i} line at 6569.2\,Å and {$\approx$3.2\%} in the nearby continuum. The continuum at 3600\,Å also shows an enhancement of $\approx$4.7\%. The white-light brightening kernels are mainly located at the flare ribbons and co-spatial with nonthermal hard X-ray sources, which implies that the enhanced white-light emissions are related to nonthermal electron-beam heating. At the brightening kernels, the Fe \textsc{i} line displays an absorption profile that has a good Gaussian shape, with a redshift up to $\approx$1.7 km s$^{-1}$, while the H$α$ line shows an emission profile though having a central reversal. The H$α$ line profile also shows a red or blue asymmetry caused by plasma flows with a velocity of several to tens of km s$^{-1}$. It is interesting to find that the H$α$ asymmetry is opposite at the conjugate footpoints. It is also found that the CHASE continuum increase seems to be related to the change of photospheric magnetic field. Our study provides comprehensive characteristics of a small white-light flare that help understand the energy release process of white-light flares. △ Less

Submitted 2 May, 2024; originally announced May 2024.

Comments: 23 pages, 6 figures, accepted by Solar Physics

JWST/NIRCam Detection of the Fomalhaut C Debris Disk in Scattered Light

Authors: Kellen Lawson, Joshua E. Schlieder, Jarron M. Leisenring, Ell Bogat, Charles A. Beichman, Geoffrey Bryden, András Gáspár, Tyler D. Groff, Michael W. McElwain, Michael R. Meyer, Thomas Barclay, Per Calissendorff, Matthew De Furio, Yiting Li, Marcia J. Rieke, Marie Ygouf, Thomas P. Greene, Julien H. Girard, Mario Gennaro, Jens Kammerer, Armin Rest, Thomas L. Roellig, Ben Sunnquist

Abstract: Observations of debris disks offer important insights into the formation and evolution of planetary systems. Though M dwarfs make up approximately 80% of nearby stars, very few M-dwarf debris disks have been studied in detail -- making it unclear how or if the information gleaned from studying debris disks around more massive stars extends to the more abundant M dwarf systems. We report the first… ▽ More Observations of debris disks offer important insights into the formation and evolution of planetary systems. Though M dwarfs make up approximately 80% of nearby stars, very few M-dwarf debris disks have been studied in detail -- making it unclear how or if the information gleaned from studying debris disks around more massive stars extends to the more abundant M dwarf systems. We report the first scattered-light detection of the debris disk around the M4 star Fomalhaut C using JWST's Near Infrared Camera (NIRCam; 3.6$~μ$m and 4.4$~μ$m). This result adds to the prior sample of only four M-dwarf debris disks with detections in scattered light, and marks the latest spectral type and oldest star among them. The size and orientation of the disk in these data are generally consistent with the prior ALMA sub-mm detection. Though no companions are identified, these data provide strong constraints on their presence -- with sensitivity sufficient to recover sub-Saturn mass objects in the vicinity of the disk. This result illustrates the unique capability of JWST for uncovering elusive M-dwarf debris disks in scattered light, and lays the groundwork for deeper studies of such objects in the 2--5$~μ$m regime. △ Less

Submitted 1 May, 2024; originally announced May 2024.

Comments: 16 pages, 5 figures

Polarization Perspectives on Hercules X-1: Further Constraining the Geometry

Authors: Qingchang Zhao, Hancheng Li, Lian Tao, Hua Feng, Shuangnan Zhang, Roland Walter, Mingyu Ge, Hao Tong, Long Ji, Liang Zhang, Jinlu Qu, Yue Huang, Xiang Ma, Shu Zhang, Qianqing Yin, Hongxing Yin, Ruican Ma, Shujie Zhao, Panping Li, Zixu Yang, Hexin Liu, Wei Yu, Yiming Huang, Zexi Li, Yajun Li , et al. (2 additional authors not shown)

Abstract: We conduct a comprehensive analysis of the accreting X-ray pulsar, Hercules X-1, utilizing data from IXPE and NuSTAR. IXPE performed five observations of Her X-1, consisting of three in the Main-on state and two in the Short-on state. Our time-resolved analysis uncovers the linear correlations between the flux and polarization degree as well as the pulse fraction and polarization degree. Geometry… ▽ More We conduct a comprehensive analysis of the accreting X-ray pulsar, Hercules X-1, utilizing data from IXPE and NuSTAR. IXPE performed five observations of Her X-1, consisting of three in the Main-on state and two in the Short-on state. Our time-resolved analysis uncovers the linear correlations between the flux and polarization degree as well as the pulse fraction and polarization degree. Geometry parameters are rigorously constrained by fitting the phase-resolved modulations of Cyclotron Resonance Scattering Feature and polarization angle with a simple dipole model and Rotating Vector Model respectively, yielding roughly consistent results. The changes of $χ_{\rm p}$ (the position angle of the pulsar's spin axis on the plane of the sky) between different Main-on observations suggest the possible forced precession of the neutron star crust. Furthermore, a linear association between the energy of Cyclotron Resonance Scattering Feature and polarization angle implies the prevalence of a dominant dipole magnetic field, and their phase-resolved modulations likely arise from viewing angle effects. △ Less

Submitted 1 May, 2024; originally announced May 2024.

Comments: Accepted for MNRAS

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.

Potential Paradigm Shift in Hazard Risk Management: AI-Based Weather Forecast for Tropical Cyclone Hazards

Authors: Kairui Feng, Dazhi Xi, Wei Ma, Cao Wang, Yuanlong Li, Xuanhong Chen

Abstract: The advents of Artificial Intelligence (AI)-driven models marks a paradigm shift in risk management strategies for meteorological hazards. This study specifically employs tropical cyclones (TCs) as a focal example. We engineer a perturbation-based method to produce ensemble forecasts using the advanced Pangu AI weather model. Unlike traditional approaches that often generate fewer than 20 scenario… ▽ More The advents of Artificial Intelligence (AI)-driven models marks a paradigm shift in risk management strategies for meteorological hazards. This study specifically employs tropical cyclones (TCs) as a focal example. We engineer a perturbation-based method to produce ensemble forecasts using the advanced Pangu AI weather model. Unlike traditional approaches that often generate fewer than 20 scenarios from Weather Research and Forecasting (WRF) simulations for one event, our method facilitates the rapid nature of AI-driven model to create thousands of scenarios. We offer open-source access to our model and evaluate its effectiveness through retrospective case studies of significant TC events: Hurricane Irma (2017), Typhoon Mangkhut (2018), and TC Debbie (2017), affecting regions across North America, East Asia, and Australia. Our findings indicate that the AI-generated ensemble forecasts align closely with the European Centre for Medium-Range Weather Forecasts (ECMWF) ensemble predictions up to seven days prior to landfall. This approach could substantially enhance the effectiveness of weather forecast-driven risk analysis and management, providing unprecedented operational speed, user-friendliness, and global applicability. △ Less

Submitted 29 April, 2024; originally announced April 2024.

Widespread rapid quenching at cosmic noon revealed by JWST deep spectroscopy

Authors: Minjung Park, Sirio Belli, Charlie Conroy, Benjamin D. Johnson, Rebecca L. Davies, Joel Leja, Sandro Tacchella, J. Trevor Mendel, Chloë Benton, Letizia Bugiani, Razieh Emami, Amirhossein Khoram, Yijia Li, Gabriel Maheson, Elijah P. Mathews, Rohan P. Naidu, Erica J. Nelson, Bryan A. Terrazas, Rainer Weinberger

Abstract: Massive quiescent galaxies in the young universe are expected to be quenched rapidly, but it is unclear whether they all experience starbursts before quenching and what physical mechanism drives rapid quenching. We study 16 massive quiescent galaxies ($\log(M_\star/M_\odot) > 10$) at $z\sim2$ selected from a representative sample of the Blue Jay survey. We reconstruct their star formation historie… ▽ More Massive quiescent galaxies in the young universe are expected to be quenched rapidly, but it is unclear whether they all experience starbursts before quenching and what physical mechanism drives rapid quenching. We study 16 massive quiescent galaxies ($\log(M_\star/M_\odot) > 10$) at $z\sim2$ selected from a representative sample of the Blue Jay survey. We reconstruct their star formation histories by fitting spectral energy distribution models to the JWST/NIRSpec $R\sim1000$ spectra. We find that massive quiescent galaxies can be split into three categories with roughly equal numbers of galaxies according to their SFHs: 1) Relatively old galaxies quenched at early epochs; 2) Galaxies that are rapidly and recently quenched after a flat or bursty formation history (depending on the assumed prior); 3) Galaxies that are rapidly and recently quenched after a major starburst. Most recently quenched galaxies show neutral gas outflows, probed by blueshifted $\rm Na\,I\,D$ absorption, and ionized gas emission, with line ratios consistent with active galactic nucleus (AGN) diagnostics. This suggests that AGN activity drives multi-phase gas outflows, leading to rapid quenching. By tracing back the SFHs of the entire sample, we predict the number density of massive quiescent galaxies at $z=4-6$: $n=3.0\pm1.4\times10^{-5}\,\rm Mpc^{-3}$. The two oldest massive quiescent galaxies in our sample appear to have extremely early formation and quenching ($z\gtrsim6$), possibly descendants of early post-starbursts at $z>3$. These galaxies still show neutral gas reservoirs and low-level star formation, consistent with weak H$α$ emission, perhaps because the ejective AGN feedback that caused rapid quenching has weakened over time. △ Less

Submitted 27 April, 2024; originally announced April 2024.

Comments: Submitted to ApJ; 21 pages, 9 figures, 3 tables (excluding appendix)

Radio Frequency Interference from Radio Navigation Satellite Systems: simulations and comparison to MeerKAT single-dish data

Authors: Brandon Engelbrecht, Mario G. Santos, José Fonseca, Yichao Li, Jingying Wang, Melis O. Irfan, Stuart E. Harper, Keith Grainge, Philip Bull, Isabella P. Carucci, Steven Cunnington, Alkistis Pourtsidou, Marta Spinelli, Laura Wolz

Abstract: Radio Frequency Interference (RFI) is emitted from various sources, terrestrial or orbital, and create a nuisance for ground-based 21cm experiments. In particular, single-dish 21cm intensity mapping experiments will be highly susceptible to contamination from these sources due to its wide primary beam and sensitivity. This work aims to simulate the contamination effects emitted from orbital source… ▽ More Radio Frequency Interference (RFI) is emitted from various sources, terrestrial or orbital, and create a nuisance for ground-based 21cm experiments. In particular, single-dish 21cm intensity mapping experiments will be highly susceptible to contamination from these sources due to its wide primary beam and sensitivity. This work aims to simulate the contamination effects emitted from orbital sources in the Radio Navigational Satellite System within the 1100-1350 MHz frequency. This simulation can be split into two parts: (I) satellite positioning, emission power, and beam response on the telescope and (II) fitting of the satellite signal to data in order to improve the original model. We use previously observed single dish MeerKAT L-band data which needs to be specially calibrated to include data contaminated by satellite-based RFI. We find that due to non-linearity effects, it becomes non-trivial to fit the satellite power. However, when masking regions where this non-linearity is problematic, we can recreate the satellite contamination with high accuracy around its peak frequencies. The simulation can predict satellite movements and signal for past and future observations, which can help in RFI avoidance and testing novel cleaning methods. The predicted signal from simulations sits below the noise in the target cosmology window for the L-band (970 - 1015 MHz) making it difficult to confirm any out-of-band emission from satellites. However, a power spectrum analysis shows that such signal can still contaminate the 21cm power spectrum at these frequencies. In our simulations, this contamination overwhelms the auto-power spectrum but still allows for a clean detection of the signal in cross-correlations with mild foreground cleaning. Whether such contamination does exist one will require further characterization of the satellite signals far away from their peak frequencies. △ Less

Submitted 27 April, 2024; originally announced April 2024.

Comments: 20 pages, 27 figures, 4 tables, Comments welcome

Sensitivity-Improved Polarization Maps at 40 GHz with CLASS and WMAP data

Authors: Rui Shi, John W. Appel, Charles L. Bennett, Ricardo Bustos, David T. Chuss, Sumit Dahal, Jullianna Denes Couto, Joseph R. Eimer, Thomas Essinger-Hileman, Kathleen Harrington, Jeffrey Iuliano, Yunyang Li, Tobias A. Marriage, Matthew A. Petroff, Karwan Rostem, Zeya Song, Deniz A. N. Valle, Duncan J. Watts, Janet L. Weiland, Edward J. Wollack, Zhilei Xu

Abstract: Improved polarization measurements at frequencies below 70 GHz with degree-level angular resolution are crucial for advancing our understanding of the Galactic synchrotron radiation and the potential polarized anomalous microwave emission and ultimately benefiting the detection of primordial $B$ modes. In this study, we present sensitivity-improved 40 GHz polarization maps obtained by combining th… ▽ More Improved polarization measurements at frequencies below 70 GHz with degree-level angular resolution are crucial for advancing our understanding of the Galactic synchrotron radiation and the potential polarized anomalous microwave emission and ultimately benefiting the detection of primordial $B$ modes. In this study, we present sensitivity-improved 40 GHz polarization maps obtained by combining the CLASS 40 GHz and WMAP $Q$-band data through a weighted average in the harmonic domain. The decision to include WMAP $Q$-band data stems from similarities in the bandpasses. Leveraging the accurate large-scale measurements from WMAP $Q$ band and the high-sensitivity information from CLASS 40 GHz band at intermediate scales, the noise level at $\ell\in[30, 100]$ is reduced by a factor of $2-3$ in the map space. A pixel domain analysis of the polarized synchrotron spectral index ($β_s$) using WMAP $K$ band and the combined maps (mean and 16/84th percentile across the $β_s$ map: $-3.08_{-0.20}^{+0.20}$) reveals a stronger preference for spatial variation (PTE for a uniform $β_s$ hypothesis smaller than 0.001) than the results obtained using WMAP $K$ and $Ka$ bands ($-3.08_{-0.14}^{+0.14}$). The cross-power spectra of the combined maps follow the same trend as other low-frequency data, and validation through simulations indicates negligible bias introduced by the combination method (sub-percent level in the power spectra). The products of this work are publicly available on $\mathtt{LAMBDA}$. △ Less

Submitted 26 April, 2024; originally announced April 2024.

Comments: 19 pages, 16 figures, 1 table

Calibrating non-parametric morphological indicators from {\it JWST} images for galaxies over $0.5<z<3$

Authors: Jian Ren, F. S. Liu, Nan Li, Qifan Cui, Pinsong Zhao, Yubin Li, Qi Song, Hassen M. Yesuf, Xian Zhong Zheng

Abstract: The measurements of morphological indicators of galaxies are often influenced by a series of observational effects. In this study, we utilize a sample of over 800 TNG50 simulated galaxies with log($M_*$/M$_\odot$)$>9$ at $0.5<z<3$ to investigate the differences in non-parametric morphological indicators ($C$, $S$, $Gini$, $M_{\rm 20}$, $A_{\rm O}$, and $D_{\rm O}$) derived from noise-free and high… ▽ More The measurements of morphological indicators of galaxies are often influenced by a series of observational effects. In this study, we utilize a sample of over 800 TNG50 simulated galaxies with log($M_*$/M$_\odot$)$>9$ at $0.5<z<3$ to investigate the differences in non-parametric morphological indicators ($C$, $S$, $Gini$, $M_{\rm 20}$, $A_{\rm O}$, and $D_{\rm O}$) derived from noise-free and high-resolution TNG50 images and mock images simulated to have the same observational conditions as {\it JWST}/NIRCam. We quantify the relationship between intrinsic and observed values of the morphological indicators and accordingly apply this calibration to over 4600 galaxies in the same stellar mass and redshift ranges observed in {\it JWST} CEERS and JADES surveys. We find a significant evolution of morphological indicators with rest-frame wavelength ($λ_{\rm rf}$) at $λ_{\rm rf}<1$\,$μ$m, while essentially no obvious variations occur at $λ_{\rm rf}>1$\,$μ$m. The morphological indicators of star-forming galaxies (SFGs) and quiescent galaxies (QGs) are significantly different. The morphologies of QGs exhibit a higher sensitivity to rest-frame wavelength than SFGs. After analyzing the evolution of morphological indicators in the rest-frame V-band (0.5-0.7\,$μ$m) and rest-frame J-band (1.1-1.4\,$μ$m), we find that the morphologies of QGs evolve substantially with both redshift and stellar mass. For SFGs, the $C$, $Gini$ and $M_{\rm 20}$ show a rapid evolution with stellar mass at log($M_*$/M$_\odot$)$\geq10.5$, while the $A_{\rm O}$, $D_{\rm O}$ and $A$ evolve with both redshift and stellar mass. Our comparison shows that TNG50 simulations effectively reproduce the morphological indicators we measured from {\it JWST} observations when the impact of dust attenuation is considered. △ Less

Submitted 25 April, 2024; originally announced April 2024.

Comments: 21 pages, 14 figures, 1 table. Accepted for publication in ApJ

Soft X-ray prompt emission from a high-redshift gamma-ray burst EP240315a

Authors: Y. Liu, H. Sun, D. Xu, D. S. Svinkin, J. Delaunay, N. R. Tanvir, H. Gao, C. Zhang, Y. Chen, X. -F. Wu, B. Zhang, W. Yuan, J. An, G. Bruni, D. D. Frederiks, G. Ghirlanda, J. -W. Hu, A. Li, C. -K. Li, J. -D. Li, D. B. Malesani, L. Piro, G. Raman, R. Ricci, E. Troja , et al. (170 additional authors not shown)

Abstract: Long gamma-ray bursts (GRBs) are believed to originate from core collapse of massive stars. High-redshift GRBs can probe the star formation and reionization history of the early universe, but their detection remains rare. Here we report the detection of a GRB triggered in the 0.5--4 keV band by the Wide-field X-ray Telescope (WXT) on board the Einstein Probe (EP) mission, designated as EP240315a,… ▽ More Long gamma-ray bursts (GRBs) are believed to originate from core collapse of massive stars. High-redshift GRBs can probe the star formation and reionization history of the early universe, but their detection remains rare. Here we report the detection of a GRB triggered in the 0.5--4 keV band by the Wide-field X-ray Telescope (WXT) on board the Einstein Probe (EP) mission, designated as EP240315a, whose bright peak was also detected by the Swift Burst Alert Telescope and Konus-Wind through off-line analyses. At a redshift of $z=4.859$, EP240315a showed a much longer and more complicated light curve in the soft X-ray band than in gamma-rays. Benefiting from a large field-of-view ($\sim$3600 deg$^2$) and a high sensitivity, EP-WXT captured the earlier engine activation and extended late engine activity through a continuous detection. With a peak X-ray flux at the faint end of previously known high-$z$ GRBs, the detection of EP240315a demonstrates the great potential for EP to study the early universe via GRBs. △ Less

Submitted 25 April, 2024; originally announced April 2024.

Comments: 41 pages, 8 figures, 7 tables

New Timing Results of MSPs from NICER Observations

Authors: Shijie Zheng, Dawei Han, Heng Xu, Kejia Lee, Jianping Yuan, Haoxi Wang, Mingyu Ge, Liang Zhang, Yongye Li, Yitao Yin, Xiang Ma, Yong Chen, Shuangnan Zhang

Abstract: Millisecond pulsars (MSPs) are known for their long-term stability. Using six years of observations from the Neutron Star Interior Composition Explorer (NICER), we have conducted an in-depth analysis of the X-ray timing results for six MSPs: PSRs B1937+21, B1821$-$24, J0437$-$4715, J0030+0451, J0218+4232, and J2124$-$3358. The timing stability parameter $σ_z$ has been calculated, revealing remarka… ▽ More Millisecond pulsars (MSPs) are known for their long-term stability. Using six years of observations from the Neutron Star Interior Composition Explorer (NICER), we have conducted an in-depth analysis of the X-ray timing results for six MSPs: PSRs B1937+21, B1821$-$24, J0437$-$4715, J0030+0451, J0218+4232, and J2124$-$3358. The timing stability parameter $σ_z$ has been calculated, revealing remarkable timing precision on the order of $10^{-14}$ for PSRs B1937+21 and J0437$-$4715, and $10^{-13}$ for PSRs B1821$-$24, J0218+4232, and J0030+0451 over a timescale of 1000 days. These findings underscore the feasibility of autonomous in-orbit timekeeping using X-ray observations of MSPs. In addition, the consistency of long-term spin-down noise in the X-ray and radio bands has been investigated by comparison with IPTA radio data. △ Less

Submitted 24 April, 2024; originally announced April 2024.

Millicharged Condensates on Earth

Authors: Asher Berlin, Roni Harnik, Ying-Ying Li, Bin Xu

Abstract: We demonstrate that long-ranged terrestrial electric fields can be used to exclude or discover ultralight bosonic particles with extremely small charge, beyond that probed by astrophysics. Bound condensates of scalar millicharged particles can be rapidly produced near electrostatic generators or in the atmosphere. If such particles directly couple to the photon, they quickly short out such electri… ▽ More We demonstrate that long-ranged terrestrial electric fields can be used to exclude or discover ultralight bosonic particles with extremely small charge, beyond that probed by astrophysics. Bound condensates of scalar millicharged particles can be rapidly produced near electrostatic generators or in the atmosphere. If such particles directly couple to the photon, they quickly short out such electrical activity. Instead, for interactions mediated by a kinetically-mixed dark photon, the effects of this condensate are suppressed depending on the size of the kinetic mixing, but may still be directly detected with precision electromagnetic sensors. Analogous condensates can also develop in other theories involving new long-ranged forces, such as those coupled to baryon and lepton number. △ Less

Submitted 24 April, 2024; originally announced April 2024.

Comments: 8 pages, 2 figures

Report number: FERMILAB-PUB-24-0002-SQMS-T, USTC-ICTS/PCFT-24-08

Cross-Temporal Spectrogram Autoencoder (CTSAE): Unsupervised Dimensionality Reduction for Clustering Gravitational Wave Glitches

Authors: Yi Li, Yunan Wu, Aggelos K. Katsaggelos

Abstract: The advancement of The Laser Interferometer Gravitational-Wave Observatory (LIGO) has significantly enhanced the feasibility and reliability of gravitational wave detection. However, LIGO's high sensitivity makes it susceptible to transient noises known as glitches, which necessitate effective differentiation from real gravitational wave signals. Traditional approaches predominantly employ fully s… ▽ More The advancement of The Laser Interferometer Gravitational-Wave Observatory (LIGO) has significantly enhanced the feasibility and reliability of gravitational wave detection. However, LIGO's high sensitivity makes it susceptible to transient noises known as glitches, which necessitate effective differentiation from real gravitational wave signals. Traditional approaches predominantly employ fully supervised or semi-supervised algorithms for the task of glitch classification and clustering. In the future task of identifying and classifying glitches across main and auxiliary channels, it is impractical to build a dataset with manually labeled ground-truth. In addition, the patterns of glitches can vary with time, generating new glitches without manual labels. In response to this challenge, we introduce the Cross-Temporal Spectrogram Autoencoder (CTSAE), a pioneering unsupervised method for the dimensionality reduction and clustering of gravitational wave glitches. CTSAE integrates a novel four-branch autoencoder with a hybrid of Convolutional Neural Networks (CNN) and Vision Transformers (ViT). To further extract features across multi-branches, we introduce a novel multi-branch fusion method using the CLS (Class) token. Our model, trained and evaluated on the GravitySpy O3 dataset on the main channel, demonstrates superior performance in clustering tasks when compared to state-of-the-art semi-supervised learning methods. To the best of our knowledge, CTSAE represents the first unsupervised approach tailored specifically for clustering LIGO data, marking a significant step forward in the field of gravitational wave research. The code of this paper is available at https://github.com/Zod-L/CTSAE △ Less

Submitted 23 April, 2024; originally announced April 2024.

VLBI with SKA: Possible Arrays and Astrometric Science

Authors: Yingjie Li, Ye Xu, Jingjing Li, Shuaibo Bian, Zehao Lin, Chaojie Hao, Dejian Liu

Abstract: The next generation of very long baseline interferometry (VLBI) is stepping into the era of microarcsecond ($μ$as) astronomy, and pushing astronomy, especially astrometry, to new heights. VLBI with the Square Kilometre Array (SKA), SKA-VLBI, will increase current sensitivity by an order of magnitude, and reach astrometric precision routinely below 10 $μ$as, even challenging 1 $μ$as. This advanceme… ▽ More The next generation of very long baseline interferometry (VLBI) is stepping into the era of microarcsecond ($μ$as) astronomy, and pushing astronomy, especially astrometry, to new heights. VLBI with the Square Kilometre Array (SKA), SKA-VLBI, will increase current sensitivity by an order of magnitude, and reach astrometric precision routinely below 10 $μ$as, even challenging 1 $μ$as. This advancement allows precise parallax and proper motion measurements of various celestial objects. Such improvements can be used to study objects (including isolated objects, and binary or multiple systems) in different stellar stages (such as star formation, main-sequence stars, asymptotic giant branch stars, pulsars, black holes, white dwarfs, etc.), unveil the structure and evolution of complex systems (such as the Milky Way), benchmark the international celestial reference frame, and reveal cosmic expansion. Furthermore, the theory of general relativity can also be tested with SKA-VLBI using precise measurements of light deflection under the gravitational fields of different solar system objects and the perihelion precession of solar system objects. △ Less

Submitted 22 April, 2024; originally announced April 2024.

Comments: 41 pages, 12 figures, 4 tables. Accepted to RAA (Review)

Reconstructing Intrinsic Stellar Noise with Stellar Atmospheric Parameters and Chromospheric Activity

Authors: Jinghua Zhang, Maosheng Xiang, Jie Yu, Jian Ge, Ji-Wei Xie, Hui Zhang, Yaguang Li, You Wu, Chun-Qian Li, Shaolan Bi, Hong-Liang Yan, Jian-Rong Shi

Abstract: Accurately characterizing intrinsic stellar photometric noise induced by stellar astrophysics, such as stellar activity, granulation, and oscillations, is of crucial importance for detecting transiting exoplanets. In this study, we investigate the relation between the intrinsic stellar photometric noise, as quantified by the Kepler rrmsCDPP measurement, and the level of stellar chromospheric activ… ▽ More Accurately characterizing intrinsic stellar photometric noise induced by stellar astrophysics, such as stellar activity, granulation, and oscillations, is of crucial importance for detecting transiting exoplanets. In this study, we investigate the relation between the intrinsic stellar photometric noise, as quantified by the Kepler rrmsCDPP measurement, and the level of stellar chromospheric activity, as indicated by the S-index of Ca II HK lines derived from the LAMOST spectra. Our results reveal a clear positive correlation between S-index and rrmsCDPP, and the correlation becomes more significant at higher activity levels and on longer timescales. We have therefore built an empirical relation between rrmsCDPP and S-index as well as Teff, logg, [Fe/H], and apparent magnitude with the XGBoost regression algorithm, using the LAMOST-Kepler common star sample as the training set. This method achieves a precision of ~20 ppm for inferring the intrinsic noise from the S-index and other stellar labels on a 6-hour integration duration. We have applied this empirical relation to the full LAMOST DR7 spectra database, and obtained the intrinsic noise predictions for 1,358,275 stars. The resultant catalog is publicly available and expected to be valuable for optimizing target selection for future exoplanet-hunting space missions, such as the Earth 2.0 mission. △ Less

Submitted 21 April, 2024; originally announced April 2024.

Comments: Accepted for publication in ApJS

Deep Learning and LLM-based Methods Applied to Stellar Lightcurve Classification

Authors: Yu-Yang Li, Yu Bai, Cunshi Wang, Mengwei Qu, Ziteng Lu, Roberto Soria, Jifeng Liu

Abstract: Light curves serve as a valuable source of information on stellar formation and evolution. With the rapid advancement of machine learning techniques, it can be effectively processed to extract astronomical patterns and information. In this study, we present a comprehensive evaluation of deep-learning and large language model (LLM) based models for the automatic classification of variable star ligh… ▽ More Light curves serve as a valuable source of information on stellar formation and evolution. With the rapid advancement of machine learning techniques, it can be effectively processed to extract astronomical patterns and information. In this study, we present a comprehensive evaluation of deep-learning and large language model (LLM) based models for the automatic classification of variable star light curves, based on large datasets from the Kepler and K2 missions. Special emphasis is placed on Cepheids, RR Lyrae, and eclipsing binaries, examining the influence of observational cadence and phase distribution on classification precision. Employing AutoDL optimization, we achieve striking performance with the 1D-Convolution+BiLSTM architecture and the Swin Transformer, hitting accuracies of 94\% and 99\% correspondingly, with the latter demonstrating a notable 83\% accuracy in discerning the elusive Type II Cepheids-comprising merely 0.02\% of the total dataset.We unveil StarWhisper LightCurve (LC), an innovative Series comprising three LLM-based models: LLM, multimodal large language model (MLLM), and Large Audio Language Model (LALM). Each model is fine-tuned with strategic prompt engineering and customized training methods to explore the emergent abilities of these models for astronomical data. Remarkably, StarWhisper LC Series exhibit high accuracies around 90\%, significantly reducing the need for explicit feature engineering, thereby paving the way for streamlined parallel data processing and the progression of multifaceted multimodal models in astronomical applications. The study furnishes two detailed catalogs illustrating the impacts of phase and sampling intervals on deep learning classification accuracy, showing that a substantial decrease of up to 14\% in observation duration and 21\% in sampling points can be realized without compromising accuracy by more than 10\%. △ Less

Submitted 16 April, 2024; originally announced April 2024.

Comments: 35 pages, 20 figures

Dynamical Model of Rotation and Orbital Coupling for Deimos

Authors: Kai Huang, Lijun Zhang, Yongzhang Yang, Mao Ye, Yuqiang Li

Abstract: This paper introduces a novel dynamical model, building upon the existing dynamical model for Deimos in the current numerical ephemerides, which only encompasses the simple libration effects of Deimos. The study comprehensively incorporates the rotational dynamics of Deimos influenced by the torque exerted by the major celestial bodies (Mars, the Sun) in the solar system within the inertial space.… ▽ More This paper introduces a novel dynamical model, building upon the existing dynamical model for Deimos in the current numerical ephemerides, which only encompasses the simple libration effects of Deimos. The study comprehensively incorporates the rotational dynamics of Deimos influenced by the torque exerted by the major celestial bodies (Mars, the Sun) in the solar system within the inertial space. Consequently, a full dynamical model is formulated to account for the complete coupling between the rotation and orbit of Deimos. Simultaneously, employing precision orbit determination methods used for artificial satellites, we develop an adjustment model for fitting data to the complete model. The 12-order Adams--Bashforth--Moulton (ABM) integration algorithm is employed to synchronously integrate the 12 state variables of the full model to obtain the orbit of Deimos. Numerical simulation results indicate that the full dynamical model and adjustment model are stable and reliable. Compared to the simple model, the polar axis of Deimos in the inertial space exhibits a more complex oscillation in the full model. This work further advances the current dynamical model for Deimos and establishes the foundational model for the generation of a new set of precise numerical ephemerides for Deimos. △ Less

Submitted 15 April, 2024; originally announced April 2024.

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

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

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

Submitted 15 April, 2024; originally announced April 2024.

Comments: 26 pages, 16 figures

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 15 April, 2024; originally announced April 2024.

Comments: 10 pages,9 figures