Showing 1–50 of 9,549 results for author: Li

Galaxy Mergers in the Epoch of Reionization I: A JWST Study of Pair Fractions, Merger Rates, and Stellar Mass Accretion Rates at $z = 4.5-11.5$

Authors: Qiao Duan, Christopher J. Conselice, Qiong Li, Duncan Austin, Thomas Harvey, Nathan J. Adams, Kenneth J. Duncan, James Trussler, Leonardo Ferreira, Lewi Westcott, Honor Harris, Rogier A. Windhorst, Benne W. Holwerda, Thomas J. Broadhurst, Dan Coe, Seth H. Cohen, Simon P. Driver, Brenda Frye, Norman A. Grogin, Nimish P. Hathi, Rolf A. Jansen, Anton M. Koekemoer, Madeline A. Marshall, Mario Nonino, Rafael Ortiz III , et al. (7 additional authors not shown)

Abstract: We present a full analysis of galaxy major merger pair fractions, merger rates, and mass accretion rates, thus uncovering the role of mergers in galaxy formation at the earliest previously unexplored epoch of $4.5<z<11.5$. We target galaxies with masses $\log_{10}(\mathrm{M}_*/\mathrm{M}_\odot) = 8.0 - 10.0$, utilizing data from eight JWST Cycle-1 fields (CEERS, JADES GOODS-S, NEP-TDF, NGDEEP, GLA… ▽ More We present a full analysis of galaxy major merger pair fractions, merger rates, and mass accretion rates, thus uncovering the role of mergers in galaxy formation at the earliest previously unexplored epoch of $4.5<z<11.5$. We target galaxies with masses $\log_{10}(\mathrm{M}_*/\mathrm{M}_\odot) = 8.0 - 10.0$, utilizing data from eight JWST Cycle-1 fields (CEERS, JADES GOODS-S, NEP-TDF, NGDEEP, GLASS, El-Gordo, SMACS-0723, MACS-0416), covering an unmasked area of 189.36 $\mathrm{arcmin}^2$. We develop a new probabilistic pair-counting methodology that integrates full photometric redshift posteriors and corrects for detection incompleteness to quantify close pairs with physical projected separations between 20 and 50 kpc. Our analysis reveals an increase in pair fractions up to $z = 8$, reaching $0.211 \pm 0.065$, followed by a statistically flat evolution to $z = 11.5$. We find that the galaxy merger rate increases from the local Universe up to $z = 6$ and then stabilizes at a value of $\sim 6$ Gyr$^{-1}$ up to $z = 11.5$. We fit both a power-law and a power-law + exponential model to our pair fraction and merger rate redshift evolution, finding that the latter model describes the trends more accurately, particularly at $z = 8.0 - 11.5$. In addition, we measure that the average galaxy increases its stellar mass due to mergers by a factor of $2.77 \pm 0.99$ from redshift $z = 10.5$ to $z = 5.0$. Lastly, we investigate the impact of mergers on galaxy stellar mass growth, revealing that mergers contribute $71 \pm 25\%$ as much to galaxy stellar mass increases as star formation from gas. This indicates that mergers drive about half of galaxy assembly at high redshift. △ Less

Submitted 12 July, 2024; originally announced July 2024.

Comments: 30 Pages, 15 Figures, Submitted to MNRAS

How coronal mass ejections are influenced by the morphology and toroidal flux of their source magnetic flux ropes?

Authors: J. H. Guo, L. Linan, S. Poedts, Y. Guo, B. Schmieder, A. Lani, Y. W. Ni, M. Brchnelova, B. Perri, T. Baratashvili, S. T. Li, P. F. Chen

Abstract: Coronal mass ejections (CMEs) stand as intense eruptions of magnetized plasma from the Sun, playing a pivotal role in driving significant changes of the heliospheric environment. Deducing the properties of CMEs from their progenitors in solar source regions is crucial for space weather forecasting. Deducing the properties of CMEs from their progenitors in solar source regions is crucial for space… ▽ More Coronal mass ejections (CMEs) stand as intense eruptions of magnetized plasma from the Sun, playing a pivotal role in driving significant changes of the heliospheric environment. Deducing the properties of CMEs from their progenitors in solar source regions is crucial for space weather forecasting. Deducing the properties of CMEs from their progenitors in solar source regions is crucial for space weather forecasting. The primary objective of this paper is to establish a connection between CMEs and their progenitors in solar source regions, enabling us to infer the magnetic structures of CMEs before their full development. To this end, we create a dataset comprising a magnetic flux rope series with varying projection shapes, sizes and toroidal fluxes, using the Regularized Biot-Savart Laws (RBSL). Thereafter, we simulate the propagation of these flux ropes from the solar surface to a distance of 25$R_{\odot}$ with our global coronal MHD model which is named COCONUT. Our parametric survey reveals significant impacts of source flux ropes on the consequent CMEs. We find that the projection shape can influence the magnetic structures of CMEs at 20$R_{\odot}$, albeit with minimal impacts on the propagation speed. However, these impacts diminish as source flux ropes become fat. In terms of toroidal flux, our simulation results demonstrate a pronounced correlation with the propagation speed of CMEs, as well as the successfulness in erupting. This work builds the bridge between the CMEs in the outer corona and their progenitors in solar source regions. Our parametric survey suggests that the projection shape, cross-section radius and toroidal flux of source flux ropes are crucial parameters in predicting magnetic structures and propagation speed of CMEs, providing valuable insights for space weather prediction. △ Less

Submitted 12 July, 2024; originally announced July 2024.

Comments: 11 pages, 10 figrues, accepted for publication by A&A

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

Detailed Mapping of the Galactic Disk Structure in the Solar Neighborhood through LAMOST K Dwarfs

Authors: Xi-Can Tang, Hao Tian, Jing Li, Bing-qiu Chen, Yi-Rong Chen, Chao Liu, Dan Qiu

Abstract: The Galactic disk is one of the main components of the Milky Way, which contributes most of the luminosity. Its structure is essential for understanding the formation and evolution of the Milky Way. Using 174,443 K-type dwarf stars observed by both LAMOST and Gaia DR3, we study the disk density profile in the local volume within 1,200 pc. In the azimuthal dimension, we find strong asymmetric signa… ▽ More The Galactic disk is one of the main components of the Milky Way, which contributes most of the luminosity. Its structure is essential for understanding the formation and evolution of the Milky Way. Using 174,443 K-type dwarf stars observed by both LAMOST and Gaia DR3, we study the disk density profile in the local volume within 1,200 pc. In the azimuthal dimension, we find strong asymmetric signal of the thin disk. The surface density and the scale height of the southern disk significantly change versus the azimuthal angle at the same galactocentric distance $R$. Meanwhile, in the vertical dimension, the scale height of the northern disk has quite different trend than that of the southern one. The scale height of the southern disk shows a decreasing trend with $φ\sim-2.5^\circ$, and change to an increasing one with $φ\sim5.0^°$. Meanwhile, the scale height of the northern disk has a consistently smaller increase. Finally, we divide the entire sample into three subsamples based on metallicity and all three subsamples show significant non-axisymmetric and north-south asymmetric signals in the Galactic disk. Furthermore, we find that the scale height of the metal-poor ([Fe/H] $<$ -0.4 dex) subsample in the northern disk is greater than that of the metal-rich ([Fe/H] $>$ -0.1 dex) subsample. However, in the southern disk, the scale height exhibits varying relationships across different metallicity slices. △ Less

Submitted 12 July, 2024; originally announced July 2024.

Comments: 15 pages, 24 figures, 6 tables; accepted for publication in MNRAS

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

The Potential Impact of Noise Correlation in Next-generation Gravitational Wave Detectors

Authors: Isaac C. F. Wong, Peter T. H. Pang, Milan Wils, Francesco Cireddu, Walter Del Pozzo, Tjonnie G. F. Li

Abstract: Building upon the statistical formulation for parameter estimation in the presence of correlated noise proposed by Cireddu et al., we present an initial study to incorporate the effects of correlated noise into the analyses of various detector designs' performance. We consider a two L-shaped detector configuration located in the European Union, and compare the expectation of parameter estimation b… ▽ More Building upon the statistical formulation for parameter estimation in the presence of correlated noise proposed by Cireddu et al., we present an initial study to incorporate the effects of correlated noise into the analyses of various detector designs' performance. We consider a two L-shaped detector configuration located in the European Union, and compare the expectation of parameter estimation between the non-colocated and a hypothetical colocated configurations. In our study, we posit the existence of low-frequency correlated noise within the $5\text{ Hz}$ to $10\text{ Hz}$ range for the colocated detector configuration, with a varying degree of correlation. In this specific detector setup, our observations indicate an enhancement in the precision of intrinsic parameter measurements as the degree of correlation increases. This trend suggests that higher degrees of noise correlation may beneficially influence the accuracy of parameter estimation. In particular, when the noise is highly correlated, the uncertainty on chirp mass decreases by up to $30\%$. The absence of an inter-European baseline does hinder the estimation of the extrinsic parameters. However, given a realistic global network with the additional detector located in the United States, the uncertainty of extrinsic parameters is significantly reduced. This reduction is further amplified as the degree of noise correlation increases. When noise correlation exceeds a certain level, the colocated configuration outperforms the non-colocated one, reducing the $90\%$ credible area of sky location by up to $10\%$. We conclude that noise correlation significantly impacts detector performance, potentially altering both quantitative and qualitative outcomes. Thus, we recommend including noise correlation in comprehensive assessments of third-generation gravitational wave detector designs. △ Less

Submitted 11 July, 2024; originally announced July 2024.

Comments: 10 pages, 2 figures

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

The white-light superflares from cool stars in GWAC triggers

Authors: Guang-Wei Li, Liang Wang, Hai-Long Yuan, Li-Ping Xin, Jing Wang, Chao Wu, Hua-Li Li, Hasitieer Haerken, Wei-Hua Wang, Hong-Bo Cai, Xu-Hui Han, Yang Xu, Lei Huang, Xiao-Meng Lu, Jian-Ying Bai, Xiang-Yu Wang, Zi-Gao Dai, En-Wei Liang, Jian-Yan Wei

Abstract: M-type stars are the ones that flare most frequently, but how big their maximum flare energy can reach is still unknown. We present 163 flares from 162 individual M2 through L1-type stars that triggered the GWAC, with flare energies ranging from $10^{32.2}$ to $10^{36.4}$ erg . The flare amplitudes range from $\triangle G = 0.84$ to $\sim 10$ mag. Flare energy increases with stellar surface temper… ▽ More M-type stars are the ones that flare most frequently, but how big their maximum flare energy can reach is still unknown. We present 163 flares from 162 individual M2 through L1-type stars that triggered the GWAC, with flare energies ranging from $10^{32.2}$ to $10^{36.4}$ erg . The flare amplitudes range from $\triangle G = 0.84$ to $\sim 10$ mag. Flare energy increases with stellar surface temperature ($T_{\rm eff}$) but both $\triangle G$ and equivalent duration $\log_{10}(ED)$ seem to be independent of $T_{\rm eff}$. Combining periods detected from light curves of TESS and K2, spectra from LAMOST, SDSS and the 2.16 m Telescope, and the Gaia DR3 data, we found that these GWAC flare stars are young. For the stars that have spectra, we found that these stars are in or very near to the saturation region, and $\log_{10}(L_{\rm Hα}/L_{\rm bol})$ is lower for M7-L1 stars than for M2-M6 stars. We also studied the relation between GWAC flare bolometric energy $E_{\rm bol}$ and stellar hemispherical area $S$, and found that $\log_{10}E_{\rm bol}$ (in erg) increases with increasing $S$ (in cm$^2$), and the maximum flare energy $\log_{10}E_{\rm bol, max} \geqslant \log_{10}S + 14.25$. For M7-L1 stars, there seem to be other factors limiting their maximum flare energies in addition to stellar hemispherical area. △ Less

Submitted 11 July, 2024; originally announced July 2024.

Comments: 18 pages, 11 figures, 4 tables

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

Deep learning solutions to telescope pointing and guiding

Authors: Jackson Zariski, Kaitlin Kratter, Sarah Logsdon, Chad Bender, Dan Li, Heidi Schweiker, Jayadev Rajagopal, Bill McBride, Emily Hunting

Abstract: The WIYN 3.5m Telescope at Kitt Peak National Observatory hosts a suite of optical and near infrared instruments, including an extreme precision, optical spectrograph, NEID, built for exoplanet radial velocity studies. In order to achieve sub ms precision, NEID has strict requirements on survey efficiency, stellar image positioning, and guiding performance, which have exceeded the native capabilit… ▽ More The WIYN 3.5m Telescope at Kitt Peak National Observatory hosts a suite of optical and near infrared instruments, including an extreme precision, optical spectrograph, NEID, built for exoplanet radial velocity studies. In order to achieve sub ms precision, NEID has strict requirements on survey efficiency, stellar image positioning, and guiding performance, which have exceeded the native capabilities of the telescope's original pointing and tracking system. In order to improve the operational efficiency of the telescope we have developed a novel telescope pointing system, built on a recurrent neural network, that does not rely on the usual pointing models (TPoint or other quasi physical bases). We discuss the development of this system, how the intrinsic properties of the pointing problem inform our network design, and show preliminary results from our best models. We also discuss plans for the generalization of this framework, so that it can be applied at other sites. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Rossby Wave Instability and Substructure Formation in 3D Non-Ideal MHD Wind-Launching Disks

Authors: Chun-Yen Hsu, Zhi-Yun Li, Yisheng Tu, Xiao Hu, Min-Kai Lin

Abstract: Rings and gaps are routinely observed in the dust continuum emission of protoplanetary discs (PPDs). How they form and evolve remains debated. Previous studies have demonstrated the possibility of spontaneous gas rings and gaps formation in wind-launching disks. Here, we show that such gas substructures are unstable to the Rossby Wave Instability (RWI) through numerical simulations. Specifically,… ▽ More Rings and gaps are routinely observed in the dust continuum emission of protoplanetary discs (PPDs). How they form and evolve remains debated. Previous studies have demonstrated the possibility of spontaneous gas rings and gaps formation in wind-launching disks. Here, we show that such gas substructures are unstable to the Rossby Wave Instability (RWI) through numerical simulations. Specifically, shorter wavelength azimuthal modes develop earlier, and longer wavelength ones dominate later, forming elongated (arc-like) anti-cyclonic vortices in the rings and (strongly magnetized) cyclonic vortices in the gaps that persist until the end of the simulation. Highly elongated vortices with aspect ratios of 10 or more are found to decay with time in our non-ideal MHD simulation, in contrast with the hydro case. This difference could be caused by magnetically induced motions, particularly strong meridional circulations with large values of the azimuthal component of the vorticity, which may be incompatible with the columnar structure preferred by vortices. The cyclonic and anti-cyclonic RWI vortices saturate at moderate levels, modifying but not destroying the rings and gaps in the radial gas distribution of the disk. In particular, they do not shut off the poloidal magnetic flux accumulation in low-density regions and the characteristic meridional flow patterns that are crucial to the ring and gap formation in wind-launching disks. Nevertheless, the RWI and their associated vortices open up the possibility of producing non-axisymmetric dust features observed in a small fraction of protoplanetary disks through non-ideal MHD, although detailed dust treatment is needed to explore this possibility. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Intensity-sensitive quality assessment of extended sources in astronomical images

Authors: X. Li, K. Adamek, W. Armour

Abstract: Radio astronomy studies the Universe by observing the radio emissions of celestial bodies. Different methods can be used to recover the sky brightness distribution (SBD), which describes the distribution of celestial sources from recorded data, with the output dependent on the method used. Image quality assessment (IQA) indexes can be used to compare the differences between restored SBDs produced… ▽ More Radio astronomy studies the Universe by observing the radio emissions of celestial bodies. Different methods can be used to recover the sky brightness distribution (SBD), which describes the distribution of celestial sources from recorded data, with the output dependent on the method used. Image quality assessment (IQA) indexes can be used to compare the differences between restored SBDs produced by different image reconstruction techniques to evaluate the effectiveness of different techniques. However, reconstructed images (for the same SBD) can appear to be very similar, especially when observed by the human visual system (HVS). Hence current structural similarity methods, inspired by the HVS, are not effective. In the past, we have proposed two methods to assess point source images, where low amounts of concentrated information are present in larger regions of noise-like data. But for images that include extended source(s), the increase in complexity of the structure makes the IQA methods for point sources over-sensitive since the important objects cannot be described by isolated point sources. Therefore, in this article we propose augmented Low-Information Similarity Index (augLISI), an improved version of LISI, to assess images including extended source(s). Experiments have been carried out to illustrate how this new IQA method can help with the development and study of astronomical imaging techniques. Note that although we focus on radio astronomical images herein, these IQA methods are also applicable to other astronomical images, and imaging techniques. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Bayesian Inference of Fine-Features of Nuclear Equation of State from Future Neutron Star Radius Measurements to 0.1km Accuracy

Authors: Bao-An Li, Xavier Grundler, Wen-Jie Xie, Nai-Bo Zhang

Abstract: To more precisely constrain the Equation of State (EOS) of supradense neutron-rich nuclear matter, future high-precision X-ray and gravitational wave observatories are proposed to measure the radii of neutron stars (NSs) with an accuracy better than about 0.1 km. However, it remains unclear what particular aspects (other than the stiffness generally spoken of in the literature) of the EOS and to w… ▽ More To more precisely constrain the Equation of State (EOS) of supradense neutron-rich nuclear matter, future high-precision X-ray and gravitational wave observatories are proposed to measure the radii of neutron stars (NSs) with an accuracy better than about 0.1 km. However, it remains unclear what particular aspects (other than the stiffness generally spoken of in the literature) of the EOS and to what precision they will be better constrained. In this work, within a Bayesian framework using a meta-model EOS for NSs, we infer the posterior probability distribution functions (PDFs) of incompressibility $K_{0}$ and skewness $J_{0}$ of symmetric nuclear matter (SNM) as well as the slope $L$, curvature $K_{\rm{sym}}$, and skewness $J_{\rm{sym}}$ characterizing the density dependence of nuclear symmetry energy $E_{\rm{sym}}(ρ)$, respectively, from mean values of NS radii consistent with existing observations and an expected accuracy $ΔR$ ranging from about 1.0 km to 0.1 km. We found that (1) the $ΔR$ has little effect on inferring the stiffness of SNM at suprasaturation densities, (2) smaller $ΔR$ reveals more accurately not only the PDFs but also pairwise correlations among parameters characterizing high-density $E_{\rm{sym}}(ρ)$, (3) a double-peak feature of the PDF($K_{\rm{sym}}$) corresponding to the strong $K_{\rm{sym}}-J_{\rm{sym}}$ and $K_{\rm{sym}}-L$ anti-correlations is revealed when $ΔR$ is less than about 0.2 km, (4) the high-precision radius measurement for canonical NSs is more useful than that for massive ones for constraining the EOS of nucleonic matter around (2-3) times the saturation density of SNM. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Comments: 21 pages with 14 figures

Revisiting the dead time effects of Insight-HXMT/ME on timing analysis

Authors: Youli Tuo, Xiaobo Li, Ying Tan, Baiyang Wu, Weichun Jiang, Liming Song, Jinlu Qu, Sudeep Gogate, Shuang-Nan Zhang, Andrea Santangelo

Abstract: Dead time is a common instrumental effect of X-ray detectors which would alter the behavior of timing properties of astronomical signals, such as distorting the shape of power density spectra (PDS), affecting the root-mean-square of potential quasi-periodic oscillation signals, etc. We revisit the effects of the dead time of Medium Energy X-ray telescope (ME) onboard Insight-HXMT, based on the sim… ▽ More Dead time is a common instrumental effect of X-ray detectors which would alter the behavior of timing properties of astronomical signals, such as distorting the shape of power density spectra (PDS), affecting the root-mean-square of potential quasi-periodic oscillation signals, etc. We revisit the effects of the dead time of Medium Energy X-ray telescope (ME) onboard Insight-HXMT, based on the simulation of electronic read-out mechanism that causes the dead time, and the real data. We investigate dead time effects on the pulse profile as well as the Quasi-Periodic Oscillation (QPO) signals. The dead time coefficient suggests a linear correlation with the observed count rate in each phase bin of the pulse profile according to the simulation of periodic signal as well as the real data observed on Swift J0243.6+6124. The Fourier-amplitude-difference (FAD) method could well recover the intrinsic shape of the observed PDS in the case that the PDS is from two identical detectors. We apply this technique on ME, by splitting the 9 FPGA modules into 2 groups. The results indicate that the FAD technique suits the case when two groups of detectors are not largely different; and the recovered PDS of Sco X-1 observed by ME slightly enhances the significance of the previously known QPO signal, meanwhile the root-mean-square of QPO is significantly improved. We provide the FAD correction tool implemented in HXMTDAS for users in the future to better analyze QPO signals. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Comments: 9 pages, 8 figures, accepted for publication in MNRAS main journal

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

Digging into the Interior of Hot Cores with ALMA (DIHCA). IV. Fragmentation in High-mass Star-Forming Clumps

Authors: Kosuke Ishihara, Patricio Sanhueza, Fumitaka Nakamura, Masao Saito, Huei-Ru V. Chen, Shanghuo Li, Fernando Olguin, Kotomi Taniguchi, Kaho Morii, Xing Lu, Qiuyi Luo, Takeshi Sakai, Qizhou Zhang

Abstract: Fragmentation contributes to the formation and evolution of stars. Observationally, high-mass stars are known to form multiple-star systems, preferentially in cluster environments. Theoretically, Jeans instability has been suggested to determine characteristic fragmentation scales, and thermal or turbulent motion in the parental gas clump mainly contributes to the instability. To search for such a… ▽ More Fragmentation contributes to the formation and evolution of stars. Observationally, high-mass stars are known to form multiple-star systems, preferentially in cluster environments. Theoretically, Jeans instability has been suggested to determine characteristic fragmentation scales, and thermal or turbulent motion in the parental gas clump mainly contributes to the instability. To search for such a characteristic fragmentation scale, we have analyzed ALMA 1.33 mm continuum observations toward 30 high-mass star-forming clumps taken by the Digging into the Interior of Hot Cores with ALMA (DIHCA) survey. We have identified 573 cores using the dendrogram algorithm and measured the separation of cores by using the Minimum Spanning Tree (MST) technique. The core separation corrected by projection effects has a distribution peaked around 5800 au. In order to remove biases produced by different distances and sensitivities, we further smooth the images to a common physical scale and perform completeness tests. Our careful analysis finds a characteristic fragmentation scale of $\sim$7000 au, comparable to the thermal Jeans length of the clumps. We conclude that thermal Jeans fragmentation plays a dominant role in determining the clump fragmentation in high-mass star-forming regions, without the need of invoking turbulent Jeans fragmentation. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: 30 pages, 18 figures, Accepted in ApJ

Inferring the spin distribution of binary black holes using deep learning

Authors: Li Tang, Xi-Long Fan

Abstract: The spin characteristics of black holes offer valuable insights into the evolutionary pathways of their progenitor stars, crucial for understanding the broader population properties of black holes. Traditional Hierarchical Bayesian inference techniques employed to discern these properties often entail substantial time investments, and consensus regarding the spin distribution of Binary Black Hole… ▽ More The spin characteristics of black holes offer valuable insights into the evolutionary pathways of their progenitor stars, crucial for understanding the broader population properties of black holes. Traditional Hierarchical Bayesian inference techniques employed to discern these properties often entail substantial time investments, and consensus regarding the spin distribution of Binary Black Hole (BBH) systems remains elusive. In this study, leveraging observations from GWTC-3, we adopt a machine learning approach to infer the spin distribution of black holes within BBH systems. Specifically, we develop a Deep Neural Network (DNN) and train it using data generated from a Beta distribution. Our training strategy, involving the segregation of data into 10 bins, not only expedites model training but also enhances the DNN's versatility and adaptability to accommodate the burgeoning volume of gravitational wave observations. Utilizing Monte Carlo-bootstrap (MC-bootstrap) to generate observation-simulated samples, we derive spin distribution parameters: $α=1.3^{+0.25}_{-0.18},β=1.70^{+0.24}_{-0.29}$ for the larger BH sample and $α=1.37^{+0.31}_{-0.20},β=1.63^{+0.30}_{-0.20}$ for the smaller BH sample. Within our constraints, the distributions of component spin magnitudes suggest the likelihood of both black holes in the BBH merger possessing non-zero spin. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: 10 pages, 8 figures, 1 table

On the Northward Shift of the Heliospheric Current Sheet at the End of Solar Cycle 24

Authors: Huichao Li, Xueshang Feng

Abstract: Since solar cycle 16, the { heliospheric} current sheet (HCS) has been found to be shifted southward during the late declining to minimum phase. However, this trend is broken at the end of solar cycle 24. In this paper, we analyze the shift of the HCS by using information obtained from coronal model and insitu data provide by the near-Earth OMNI database and the Parker Solar Probe (PSP). Coronal p… ▽ More Since solar cycle 16, the { heliospheric} current sheet (HCS) has been found to be shifted southward during the late declining to minimum phase. However, this trend is broken at the end of solar cycle 24. In this paper, we analyze the shift of the HCS by using information obtained from coronal model and insitu data provide by the near-Earth OMNI database and the Parker Solar Probe (PSP). Coronal potential field source surface (PFSS) modeling results show that the northward shift is established at the beginning of 2018 and remains stable for about two years. Interplanetary magnetic field data obtained from and within 1 au also support the northward shift, as the southern polarity T appears more frequently than the northern polarity A between 2018-2020. Both model results and insitu observation obtained by PSP imply that the HCS shift is established in the corona, and then propagates into the heliosphere. The quadrupole term still has a significant influence on the formation of the HCS shift. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: Accepted by MNRAS

GD-1 Stellar Stream and Cocoon in the DESI Early Data Release

Authors: Monica Valluri, Parker Fagrelius, Sergey. E. Koposov, Ting S. Li, Oleg Y. Gnedin, Eric F. Bell, Raymond G. Carlberg, Andrew P. Cooper, Jessia N. Aguilar, Carlos Allende Prieto, Vasily Belokurov, Leandro Beraldo e Silva, David Brooks, Amanda Byström, Todd Claybaugh, Kyle Dawson, Arjun Dey, Peter Doel, Jaime E. Forero-Romero, Enrique Gaztañaga, Satya Gontcho A Gontcho, Klaus Honscheid, T . Kisner, Anthony Kremin, A. Lambert , et al. (27 additional authors not shown)

Abstract: We present ~ 126 new spectroscopically identified members of the GD-1 tidal stream obtained with the 5000-fiber Dark Energy Spectroscopic Instrument (DESI). We confirm the existence of a ``cocoon'' which is broad (FWHM~2.932deg~460pc) and kinematically hot (velocity dispersion, sigma~5-8km/s) component that surrounds a narrower (FWHM~0.353deg~55pc) and colder (sigma~ 2.2-2.6km/s) thin stream compo… ▽ More We present ~ 126 new spectroscopically identified members of the GD-1 tidal stream obtained with the 5000-fiber Dark Energy Spectroscopic Instrument (DESI). We confirm the existence of a ``cocoon'' which is broad (FWHM~2.932deg~460pc) and kinematically hot (velocity dispersion, sigma~5-8km/s) component that surrounds a narrower (FWHM~0.353deg~55pc) and colder (sigma~ 2.2-2.6km/s) thin stream component (based on a median per star velocity precision of 2.7km/s). The cocoon extends over at least a ~ 20deg segment of the stream observed by DESI. The thin and cocoon components have similar mean values of [Fe/H]: -2.54+/- 0.04dex and -2.45+/-0.06dex suggestive of a common origin. The data are consistent with the following scenarios for the origin of the cocoon. The progenitor of the GD-1 stream was an accreted globular cluster (GC) and: (a) the cocoon was produced by pre-accretion tidal stripping of the GC while it was still inside its parent dwarf galaxy; (b) the cocoon is debris from the parent dwarf galaxy; (c) an initially thin GC tidal stream was heated by impacts from dark subhalos in the Milky Way; (d) an initially thin GC stream was heated by a massive Sagittarius dwarf galaxy; or a combination of some these. In the first two cases the velocity dispersion and mean metallicity are consistent with the parent dwarf galaxy having a halo mass of ~0^9\msun. Future DESI spectroscopy and detailed modeling may enable us to distinguish between these possible origins. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: Submitted to ApJ, 23 pages, 13 figures 4 tables

DESI Early Data Release Milky Way Survey Value-Added Catalogue

Authors: Sergey E. Koposov, C. Allende-Prieto, A. P. Cooper, T. S. Li, L. Beraldo e Silva, B. Kim, A. Carrillo, A. Dey, C. J. Manser, F. Nikakhtar, A. H. Riley, C. Rockosi, M. Valluri, J. Aguilar, S. Ahlen, S. Bailey, R. Blum, D. Brooks, T. Claybaugh, S. Cole, A. de la Macorra, B. Dey, J. E. Forero-Romero, E. Gaztañaga, J. Guy , et al. (18 additional authors not shown)

Abstract: We present the stellar value-added catalogue based on the Dark Energy Spectroscopic Instrument (DESI) Early Data Release. The catalogue contains radial velocity and stellar parameter measurements for $\simeq$ 400,000 unique stars observed during commissioning and survey validation by DESI. These observations were made under conditions similar to the Milky Way Survey (MWS) currently carried out by… ▽ More We present the stellar value-added catalogue based on the Dark Energy Spectroscopic Instrument (DESI) Early Data Release. The catalogue contains radial velocity and stellar parameter measurements for $\simeq$ 400,000 unique stars observed during commissioning and survey validation by DESI. These observations were made under conditions similar to the Milky Way Survey (MWS) currently carried out by DESI but also include multiple specially targeted fields, such as those containing well-studied dwarf galaxies and stellar streams. The majority of observed stars have $16<r<20$ with a median signal-to-noise ratio in the spectra of $\sim$ 20. In the paper, we describe the structure of the catalogue, give an overview of different target classes observed, as well as provide recipes for selecting clean stellar samples. We validate the catalogue using external high-resolution measurements and show that radial velocities, surface gravities, and iron abundances determined by DESI are accurate to 1 km/s, $0.3$ dex and $\sim$ 0.15 dex respectively. We also demonstrate possible uses of the catalogue for chemo-dynamical studies of the Milky Way stellar halo and Draco dwarf spheroidal. The value-added catalogue described in this paper is the very first DESI MWS catalogue. The next DESI data release, expected in less than a year, will add the data from the first year of DESI survey operations and will contain approximately 4 million stars, along with significant processing improvements. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: Submitted to MNRAS; Value added catalogue is available at https://data.desi.lbl.gov/public/edr/vac/edr/mws/fuji/

Fermi-LAT discovery of the GeV emission of the superluminous supernovae SN 2017egm

Authors: Shang Li, Yun-Feng Liang, Neng-Hui Liao, Lei Lei, Yi-Zhong Fan

Abstract: Superluminous supernovae (SLSNe) are a new class of transients with luminosities $\sim10 -100$ times larger than the usual core-collapse supernovae (SNe). Their origin is still unclear and one widely discussed scenario involves a millisecond magnetar central engine. The GeV-TeV emission of SLSNe has been predicted in the literature but has not been convincingly detected yet. Here we report the sea… ▽ More Superluminous supernovae (SLSNe) are a new class of transients with luminosities $\sim10 -100$ times larger than the usual core-collapse supernovae (SNe). Their origin is still unclear and one widely discussed scenario involves a millisecond magnetar central engine. The GeV-TeV emission of SLSNe has been predicted in the literature but has not been convincingly detected yet. Here we report the search for the $γ$-ray emission in the direction of SN 2017egm, one of the closest SLSNe detected so far, with the 15-year {\it Fermi}-LAT Pass 8 data. There is a transient $γ$-ray source appearing about 2 months after this event and lasting a few months. Both the peak time and the luminosity of the GeV emission are consistent with the magnetar model prediction, suggesting that such a GeV transient is the high-energy counterpart of SN 2017egm and the central engine of this SLSNe is a young magnetar. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Velocity-Resolved Ionization Mapping of Broad Line Region. I. Insights into Diverse Geometry and Kinematics

Authors: Sha-Sha Li, Hai-Cheng Feng, H. T. Liu, J. M. Bai, Xiang Ji, Cheng Cheng, Kai-Xing Lu, Jian-Guo Wang, Rui Li

Abstract: Broad emission lines of active galactic nuclei (AGNs) originate from the broad-line region (BLR), consisting of dense gas clouds in orbit around an accreting supermassive black hole. Understanding the geometry and kinematics of the region is crucial for gaining insights into the physics and evolution of AGNs. Conventional velocity-resolved reverberation mapping may face challenges in disentangling… ▽ More Broad emission lines of active galactic nuclei (AGNs) originate from the broad-line region (BLR), consisting of dense gas clouds in orbit around an accreting supermassive black hole. Understanding the geometry and kinematics of the region is crucial for gaining insights into the physics and evolution of AGNs. Conventional velocity-resolved reverberation mapping may face challenges in disentangling the degeneracy between intricate motion and geometry of this region. To address this challenge, new key constraints are required. Here, we report the discovery of an asymmetric BLR using a novel technique: velocity-resolved ionization mapping, which can map the distance of emitting gas clouds by measuring Hydrogen line ratios at different velocities. By analyzing spectroscopic monitoring data, we find that the Balmer decrement is anticorrelated with the continuum and correlated with the lags across broad emission line velocities. Some line ratio profiles deviate from the expectations for a symmetrically virialized BLR, suggesting that the red-shifted and blue-shifted gas clouds may not be equidistant from the supermassive black hole (SMBH). This asymmetric geometry might represent a formation imprint, provide new perspectives on the evolution of AGNs, and influence SMBH mass measurements. △ Less

Submitted 7 July, 2024; originally announced July 2024.

Comments: 20 pages, 10 figures, Accepted by ApJ

Tilted Disk Precession and Negative Superhumps in HS 2325+8205: A Multi-Window Analysis

Authors: Qi-Bin Sun, Sheng-Bang Qian, Li-Ying Zhu, Qin-Mei Li, Min-Yu Li, Ping Li

Abstract: Tilted disk precession exists in different objects. Negative superhumps (NSHs) in cataclysmic variable stars (CVs) are hypothesized to arise from the interaction between the reverse precession of a tilted disk and the streams from the secondary star. Utilizing TESS photometry, we present a comprehensive investigation into the tilted disk precession and NSHs in the dwarf nova (DN) HS 2325+8205, emp… ▽ More Tilted disk precession exists in different objects. Negative superhumps (NSHs) in cataclysmic variable stars (CVs) are hypothesized to arise from the interaction between the reverse precession of a tilted disk and the streams from the secondary star. Utilizing TESS photometry, we present a comprehensive investigation into the tilted disk precession and NSHs in the dwarf nova (DN) HS 2325+8205, employing eclipse minima, eclipse depths, NSH frequencies, and NSH amplitudes and the correlation between them as the windows. We report the discovery of NSHs in HS 2325+8205 with a period of 0.185671(17) d. The NSH frequency was found to vary with a period of 3.943(9) d, similar to the tilted disk precession period validated in novae-like star (NL, SDSS J0812) and intermediate polar (IP, TV Col). The eclipsing minima of O-C were similarly found to vary cyclically in period 4.135(5) d, showing faster rise than fall. Furthermore, NSH amplitude varies parabolically, linearly increasing with periodic variations, potentially linked to changes in disk radius, mass transfer rate, and apparent area of the hot spot. Additionally, for the first time in DNe, we observe bi-periodic variations in eclipse depth (P1= 4.131(4) d and P2= 2.065(2) d ~ Pprec/2), resembling those seen in IPs, suggesting that variations with P2 are not attributable to the accretion curtain. Moreover, NSH amplitude and eclipse depth decrease with increasing NSH frequency, while NSH amplitude correlates positively with eclipse depth.These complex variations observed across multiple observational windows provide substantial evidence for understanding of tilted disk precession and NSHs. △ Less

Submitted 5 July, 2024; originally announced July 2024.

Comments: 26 pages, 10 figures and 2 tables

MIRI MRS Observations of Beta Pictoris II. The Spectroscopic Case for a Recent Giant Collision

Authors: Christine H. Chen, Cicero X. Lu, Kadin Worthen, David R. Law, B. A. Sargent, Amaya Moro-Martin, G. C. Sloan, Carey M. Lisse, Dan M. Watson, Julien H. Girard, Yiwei Chai, Dean C. Hines, Jens Kammerer, Alexis Li, Marshall Perrin, Laurent Pueyo, Isabel Rebollido, Karl R. Stapelfeldt, Christopher Stark, Michael W. Werner

Abstract: Modeling observations of the archetypal debris disk around $β$ Pic, obtained in 2023 January with the MIRI MRS on board JWST, reveals significant differences compared with that obtained with the IRS on board Spitzer. The bright 5 - 15 $μ$m continuum excess modeled using a $\sim$600 K black body has disappeared. The previously prominent 18 and 23 $μ$m crystalline forsterite emission features, arisi… ▽ More Modeling observations of the archetypal debris disk around $β$ Pic, obtained in 2023 January with the MIRI MRS on board JWST, reveals significant differences compared with that obtained with the IRS on board Spitzer. The bright 5 - 15 $μ$m continuum excess modeled using a $\sim$600 K black body has disappeared. The previously prominent 18 and 23 $μ$m crystalline forsterite emission features, arising from cold dust ($\sim$100 K) in the Rayleigh limit, have disappeared and been replaced by very weak features arising from the hotter 500 K dust population. Finally, the shape of the 10 $μ$m silicate feature has changed, consistent with a shift in the temperature of the warm dust population from $\sim$300 K to $\sim$500 K and an increase in the crystalline fraction of the warm, silicate dust. Stellar radiation pressure may have blown both the hot and the cold crystalline dust particles observed in the Spitzer spectra out of the planetary system during the intervening 20 years between the Spitzer and JWST observations. These results indicate that the $β$ Pic system has a dynamic circumstellar environment, and that periods of enhanced collisions can create large clouds of dust that sweep through the planetary system. △ Less

Submitted 5 July, 2024; originally announced July 2024.

Comments: 15 pages, 8 figures, ApJ in press

Comparing metallicity correlations in nearby non-AGN and AGN-host galaxies

Authors: Song-lin Li, Zefeng Li, Emily Wisnioski, Mark R. Krumholz, Sebastián F. Sánchez

Abstract: The gas-phase metallicity distribution within galaxies records critical information about galactic evolution. In this work we investigate how active galactic nuclei (AGN) influence this distribution by measuring the two-point correlation functions of gas-phase metallicity in 95 non-AGN and 37 AGN-host galaxies from the Calar Alto Legacy Integral Field spectroscopy Area integral field spectrographi… ▽ More The gas-phase metallicity distribution within galaxies records critical information about galactic evolution. In this work we investigate how active galactic nuclei (AGN) influence this distribution by measuring the two-point correlation functions of gas-phase metallicity in 95 non-AGN and 37 AGN-host galaxies from the Calar Alto Legacy Integral Field spectroscopy Area integral field spectrographic survey. We measure metallicity using a novel Bayesian method that properly includes both stellar and AGN contributions to emission line fluxes and allows us to measure metallicities in both AGN-host and non-AGN galaxies in a single, consistent framework. We find that the two-point correlation functions of both AGN-host and non-AGN galaxies are well-fit by a simple injection-diffusion model, and that the correlation lengths $l_\mathrm{corr}$ we derive for the non-AGN galaxies are reasonably consistent with those obtained in earlier work. The AGN-host galaxies generally have smaller $l_\mathrm{corr}$ than non-AGN galaxies at fixed stellar mass, but similar $l_\mathrm{corr}$ at fixed star formation rate (SFR), suggesting that the primary effect of hosting an AGN in this sample is a reduction in SFR at fixed stellar mass, and that this in turn suppresses the correlation length. Our findings further indicate that, while both SFR and stellar mass are positively correlated with metallicity correlation length $l_\mathrm{corr}$, the former is more fundamental, implying that fluctuations in the metallicity distribution within galaxies are driven more by short-term responses to physical processes such as star formation that can change much faster than a Hubble time. △ Less

Submitted 5 July, 2024; originally announced July 2024.

Comments: 13 pages, 10 figures, submitted to MNRAS

Surviving in the Hot Neptune Desert: The Discovery of the Ultra-Hot Neptune TOI-3261b

Authors: Emma Nabbie, Chelsea X. Huang, Jennifer A. Burt, David J. Armstrong, Eric E. Mamajek, Vardan Adibekyan, Sérgio G. Sousa, Eric D. Lopez, Daniel P. Thorngren, Jorge Fernández, Gongjie Li, James S. Jenkins, Jose I. Vines, João Gomes da Silva, Robert A. Wittenmyer, Daniel Bayliss, César Briceño, Karen A. Collins, Xavier Dumusque, Keith D. Horne, Marcelo F. Keniger, Nicholas Law, Jorge Lillo-Box, Shang-Fei Liu, Andrew W. Mann , et al. (23 additional authors not shown)

Abstract: The recent discoveries of Neptune-sized ultra-short period planets (USPs) challenge existing planet formation theories. It is unclear whether these residents of the Hot Neptune Desert have similar origins to smaller, rocky USPs, or if this discrete population is evidence of a different formation pathway altogether. We report the discovery of TOI-3261b, an ultra-hot Neptune with an orbital period… ▽ More The recent discoveries of Neptune-sized ultra-short period planets (USPs) challenge existing planet formation theories. It is unclear whether these residents of the Hot Neptune Desert have similar origins to smaller, rocky USPs, or if this discrete population is evidence of a different formation pathway altogether. We report the discovery of TOI-3261b, an ultra-hot Neptune with an orbital period $P$ = 0.88 days. The host star is a $V = 13.2$ magnitude, slightly super-solar metallicity ([Fe/H] $\simeq$ 0.15), inactive K1.5 main sequence star at $d = 300$ pc. Using data from the Transiting Exoplanet Survey Satellite and the Las Cumbres Observatory Global Telescope, we find that TOI-3261b has a radius of $3.82_{-0.35}^{+0.42}$ $R_{\oplus}$. Moreover, radial velocities from ESPRESSO and HARPS reveal a mass of $30.3_{-2.4}^{+2.2}$ $M_{\oplus}$, more than twice the median mass of Neptune-sized planets on longer orbits. We investigate multiple mechanisms of mass loss that can reproduce the current-day properties of TOI-3261b, simulating the evolution of the planet via tidal stripping and photoevaporation. Thermal evolution models suggest that TOI-3261b should retain an envelope potentially enriched with volatiles constituting $\sim$5% of its total mass. This is the second highest envelope mass fraction among ultra-hot Neptunes discovered to date, making TOI-3261b an ideal candidate for atmospheric follow-up observations. △ Less

Submitted 4 July, 2024; originally announced July 2024.

Comments: 20 pages, 11 figures, accepted to AJ

Hadronuclear interactions in AGN jets as the origin of the diffuse high-energy neutrino background

Authors: Rui Xue, Ze-Rui Wang, Jagdish C. Joshi, Wei-Jian Li

Abstract: The origin of diffuse high-energy neutrinos from TeV to PeV energies detected by IceCube Observatory remains a mystery. In our previous work, we have shown that hadronuclear (p-p) interactions in AGN jets could be important and generate detectable very-high-energy emissions. Here, we further explore these interactions in the AGN jets based on their luminosity function. The diffuse neutrino flux an… ▽ More The origin of diffuse high-energy neutrinos from TeV to PeV energies detected by IceCube Observatory remains a mystery. In our previous work, we have shown that hadronuclear (p-p) interactions in AGN jets could be important and generate detectable very-high-energy emissions. Here, we further explore these interactions in the AGN jets based on their luminosity function. The diffuse neutrino flux and corresponding $γ$-ray flux have been calculated and compared with observational data. In our modeling, two beaming patterns are considered separately. To make sure that the corresponding $γ$-ray flux does not overshoot the diffuse $γ$-ray background, we find that if the neutrino production region in jet is opaque to $γ$ rays, p-p interactions in AGN jets with a small viewing angle (the blazar case) are able to interpret the PeV neutrino background. Similarly, AGN jets with a large viewing angle (the radio galaxy case) may interpret the TeV neutrino background. While, if the neutrino production region is transparent to $γ$ rays, only blazars have the potential to interpret the DNB around PeV band. Some caveats are also discussed. △ Less

Submitted 4 July, 2024; originally announced July 2024.

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

Quasi-periodic oscillations of flare loops and slipping motion of ribbon substructures during a C-class flare

Authors: Yining Zhang, Ting Li, Jing Ye

Abstract: Quasi-periodic oscillations in solar flaring emission have been observed over the past few decades. To date, the underpinning processes resulting in the quasi-periodic oscillations remain unknown. In this paper, we report a unique event that exhibits both the long-duration quasi-periodic intensity oscillations of flare loops and the quasi-periodic slipping motion of ribbon substructures during a C… ▽ More Quasi-periodic oscillations in solar flaring emission have been observed over the past few decades. To date, the underpinning processes resulting in the quasi-periodic oscillations remain unknown. In this paper, we report a unique event that exhibits both the long-duration quasi-periodic intensity oscillations of flare loops and the quasi-periodic slipping motion of ribbon substructures during a C9.1-class flare (SOL2015-03-15-T01:15), using the observations from Solar Dynamics Observatory and Interface Region Imaging Spectrograph. The high-temperature flare loops rooted in the straight part of ribbons display a "bright-dim" intensity oscillation, with a period of about 4.5 minutes. The oscillation starts just after the flare onset and lasts over 3 hours. Meanwhile, the substructures within the ribbon tip display the quasi-periodic slipping motion along the ribbon at 1400 Åimages which has a similar periodicity to the stationary intensity oscillation of the flare loops in the straight part of the flare ribbons. We suggest that the quasi-periodic pattern is probably related to the loop-top dynamics caused by the reconnection outflow impinging on the flare loops. △ Less

Submitted 4 July, 2024; originally announced July 2024.

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

A new subclass of gamma-ray burst originating from compact binary merger

Authors: Chen-Wei Wang, Wen-Jun Tan, Shao-Lin Xiong, Shu-Xu Yi, Rahim Moradi, Bing Li, Zhen Zhang, Yu Wang, Yan-Zhi Meng, Jia-Cong Liu, Yue Wang, Sheng-Lun Xie, Wang-Chen Xue, Zheng-Hang Yu, Peng Zhang, Wen-Long Zhang, Yan-Qiu Zhang, Chao Zheng

Abstract: Type I gamma-ray bursts (GRBs) are believed to originate from compact binary merger usually with duration less than 2 seconds for the main emission. However, recent observations of GRB 211211A and GRB 230307A indicate that some merger-origin GRBs could last much longer. Since they show strikingly similar properties (indicating a common mechanism) which are different from the classic "long"-short b… ▽ More Type I gamma-ray bursts (GRBs) are believed to originate from compact binary merger usually with duration less than 2 seconds for the main emission. However, recent observations of GRB 211211A and GRB 230307A indicate that some merger-origin GRBs could last much longer. Since they show strikingly similar properties (indicating a common mechanism) which are different from the classic "long"-short burst (e.g. GRB 060614), forming an interesting subclass of type I GRBs, we suggest to name them as type IL GRBs. By identifying the first peak of GRB 230307A as a quasi-thermal precursor, we find that the prompt emission of type IL GRB is composed of three episodes: (1) a precursor followed by a short quiescent (or weak emission) period, (2) a long-duration main emission, and (3) an extended emission. With this burst pattern, a good candidate, GRB 170228A, was found in the Fermi/GBM archive data, and subsequent temporal and spectral analyses indeed show that GRB 170228A falls in the same cluster with GRB 211211A and GRB 230307A in many diagnostic figures. Thus this burst pattern could be a good reference for rapidly identifying type IL GRB and conducting low-latency follow-up observation. We estimated the occurrence rate and discussed the physical origins and implications for the three emission episodes of type IL GRBs. Our analysis suggests the pre-merger precursor model, especially the super flare model, is more favored for type IL GRBs. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Turbulent Diffuse Molecular Media with Non-ideal Magnetohydrodynamics and Consistent Thermochemistry: Numerical Simulations and Dynamic Characteristics

Authors: Nannan Yue, Lile Wang, Thomas Bisbas, Donghui Quan, Di Li

Abstract: Turbulent diffuse molecular clouds can exhibit complicated morphologies caused by the interactions among radiation, chemistry, fluids, and fields. We performed full 3D simulations for turbulent diffuse molecular interstellar media, featuring time-dependent non-equilibrium thermochemistry co-evolved with magnetohydrodynamics (MHD). Simulation results exhibit the relative abundances of key chemical… ▽ More Turbulent diffuse molecular clouds can exhibit complicated morphologies caused by the interactions among radiation, chemistry, fluids, and fields. We performed full 3D simulations for turbulent diffuse molecular interstellar media, featuring time-dependent non-equilibrium thermochemistry co-evolved with magnetohydrodynamics (MHD). Simulation results exhibit the relative abundances of key chemical species (e.g., C, CO, OH) vary by more than one order of magnitude for the "premature" epoch of chemical evolution ($t\lesssim 2\times 10^5~{\rm yr}$). Various simulations are also conducted to study the impacts of physical parameters. Non-ideal MHD effects are essential in shaping the behavior of gases, and strong magnetic fields ($\sim 10~μ{\rm G}$) tend to inhibit vigorous compressions and thus reduce the fraction of warm gases ($T\gtrsim 10^2~{\rm K}$). Thermodynamical and chemical conditions of the gas are sensitive to modulation by dynamic conditions, especially the energy injection by turbulence. Chemical features, including ionization (cosmic ray and diffuse interstellar radiation), would not directly affect the turbulence power spectra. Nonetheless, their effects are prominent in the distribution profiles of temperatures and gas densities. Comprehensive observations are necessary and useful to eliminate the degeneracies of physical parameters and constrain the properties of diffuse molecular clouds with confidence. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: 21 pages, 13 figures; re-submitted the revised manuscript to the Astrophysical Journal

Unveiling Mass Transfer in Solar Flares: Insights from Elemental Abundance Evolutions Observed by Chang'E-2 Solar X-ray Monitor

Authors: Man-Hei Ng, Chi-Long Tang, Xiaoping Zhang, Kuan-Vai Tam, Peng-Fei Chen, Wudong Dong, Jing Li, Chi-Pui Tang

Abstract: Understanding how elemental abundances evolve during solar flares helps shed light on the mass and energy transfer between different solar atmospheric layers. However, prior studies have mostly concentrated on averaged abundances or specific flare phases, leaving a gap in exploring the comprehensive observations throughout the entire flare process. Consequently, investigations into this area are r… ▽ More Understanding how elemental abundances evolve during solar flares helps shed light on the mass and energy transfer between different solar atmospheric layers. However, prior studies have mostly concentrated on averaged abundances or specific flare phases, leaving a gap in exploring the comprehensive observations throughout the entire flare process. Consequently, investigations into this area are relatively scarce. Exploiting the Solar X-ray Monitor data obtained from the Chang'E-2 lunar orbiter, we present two comprehensive soft X-ray spectroscopic observations of flares in active regions, AR 11149 and 11158, demonstrating elemental abundance evolutions under different conditions. Our findings unveil the inverse first ionization potential (IFIP) effect during flares for Fe for the first time, and reaffirm its existence for Si. Additionally, we observed a rare depletion of elemental abundances, marking the second IFIP effect in flare decay phases. Our study offers a CSHKP model-based interpretation to elucidate the formation of both the FIP and IFIP effects in flare dynamics, with the inertia effect being incorporated into the ponderomotive force fractionation model. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: Accepted ApJ

Modeling the Nonlinear Power Spectrum in Low-redshift HI Intensity Mapping

Authors: Zhixing Li, Laura Wolz, Hong Guo, Steven Cunnington, Yi Mao

Abstract: We present a simulation-based framework to forecast the HI power spectrum on non-linear scales ($k\gtrsim 1\ {\rm Mpc^{-1}}$), as measured by interferometer arrays like MeerKAT in the low-redshift ($z\leq 1.0$) universe. Building on a galaxy-based HI mock catalog, we meticulously consider various factors, including the emission line profiles of HI discs and some observational settings, and explore… ▽ More We present a simulation-based framework to forecast the HI power spectrum on non-linear scales ($k\gtrsim 1\ {\rm Mpc^{-1}}$), as measured by interferometer arrays like MeerKAT in the low-redshift ($z\leq 1.0$) universe. Building on a galaxy-based HI mock catalog, we meticulously consider various factors, including the emission line profiles of HI discs and some observational settings, and explore their impacts on the HI power spectrum. While it is relatively insensitive to the profile shape of HI emission line at these scales, we identify a strong correlation with the profile width, that is, the Full Width at Half Maxima (FWHM, also known as $W_{\rm 50}$ in observations) in this work. By modeling the width function of $W_{50}$ as a function of $v_{\rm max}$, we assign each HI source a emission line profile and find that the resulting HI power spectrum is comparatively close to results from particles in the IllustrisTNG hydrodynamical simulation. After implementing $k$-space cuts matching the MeerKAT data, our prediction replicates the trend of the measurements obtained by MeerKAT at $z\approx 0.44$, though with a significantly lower amplitude. Utilizing a Monte Carlo Markov Chain sampling method, we constrain the parameter $A_{W_{\rm 50}}$ in the $W_{\rm 50}$ models and $Ω_{\rm HI}$ with the MeerKAT measurements and find that a strong degeneracy exists between these two parameters. △ Less

Submitted 2 July, 2024; originally announced July 2024.

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

Origin of the Chromospheric Umbral Waves in Sunspots

Authors: Xinsheng Zhang, Xiaoli Yan, Zhike Xue, Jincheng Wang, Zhe Xu, Qiaoling Li, Yang Peng, Liping Yang

Abstract: Oscillations are ubiquitous in sunspots and the associated higher atmospheres. However, it is still unclear whether these oscillations are driven by the external acoustic waves (p-modes) or generated by the internal magnetoconvection. To obtain clues about the driving source of umbral waves in sunspots, we analyzed the spiral wave patterns (SWPs) in two sunspots registered by IRIS MgII 2796 Å slit… ▽ More Oscillations are ubiquitous in sunspots and the associated higher atmospheres. However, it is still unclear whether these oscillations are driven by the external acoustic waves (p-modes) or generated by the internal magnetoconvection. To obtain clues about the driving source of umbral waves in sunspots, we analyzed the spiral wave patterns (SWPs) in two sunspots registered by IRIS MgII 2796 Å slit-jaw images. By tracking the motion of the SWPs, we find for the first time that two one-armed SWPs coexist in the umbra, and they can rotate either in the same or opposite directions. Furthermore, by analyzing the spatial distribution of the oscillation centers of the one-armed SWPs within the umbra (the oscillation center is defined as the location where the SWP first appears), we find that the chromospheric umbral waves repeatedly originate from the regions with high oscillation power and most of the umbral waves occur in the dark nuclei and strong magnetic field regions of the umbra. Our study results indicate that the chromospheric umbral waves are likely excited by the p-mode oscillations. △ Less

Submitted 3 July, 2024; v1 submitted 2 July, 2024; originally announced July 2024.

Understanding the Broad-line Region of Active Galactic Nuclei with Photoionization. I. the Moderate-Accretion Regime

Authors: Qiaoya Wu, Yue Shen, Hengxiao Guo, Scott F. Anderson, W. N. Brandt, Catherine J. Grier, Patrick B. Hall, Luis C. Ho, Yasaman Homayouni, Keith Horne, Jennifer I-Hsiu Li, Donald P. Schneider

Abstract: Over three decades of reverberation mapping (RM) studies on local broad-line active galactic nuclei (AGNs) have measured reliable black-hole (BH) masses for $> 100$ AGNs. These RM measurements reveal a significant correlation between the Balmer broad-line region size and the AGN optical luminosity (the $R-L$ relation). Recent RM studies for AGN samples with more diverse BH accretion parameters (e.… ▽ More Over three decades of reverberation mapping (RM) studies on local broad-line active galactic nuclei (AGNs) have measured reliable black-hole (BH) masses for $> 100$ AGNs. These RM measurements reveal a significant correlation between the Balmer broad-line region size and the AGN optical luminosity (the $R-L$ relation). Recent RM studies for AGN samples with more diverse BH accretion parameters (e.g., mass and Eddington ratio) reveal a substantial intrinsic dispersion around the average $R-L$ relation, suggesting variations in the overall spectral energy distribution shape as functions of accretion parameters. Here we perform a detailed photoionization investigation of expected broad-line properties as functions of accretion parameters, using the latest models for the AGN continuum implemented in {\tt qsosed}. We compare theoretical predictions with observations of a sample of 67 $z\lesssim0.5$ reverberation-mapped AGNs with both rest-frame optical and UV spectra in the moderate-accretion regime (Eddington ratio $λ_{\rm Edd}\equiv L/L_{\rm Edd}<0.5$). The UV/optical line strengths and their dependences on accretion parameters can be reasonably well reproduced by the locally-optimally-emitting cloud (LOC) photoionization models. We provide quantitative recipes that use optical/UV line flux ratios to infer the ionizing continuum, which is not directly observable. In addition, photoionization models with universal values of ionization parameter ($\log U_{\rm H}=-2$) and hydrogen density ($\log n({\rm H})=12$) can qualitatively reproduce the observed global $R-L$ relation for the current AGN sample. However, such models fail to reproduce the observed trend of decreasing BLR size with $L/L_{\rm Edd}$ at fixed optical luminosity, which may imply that the gas density increases with the accretion rate. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Comments: 23 pages, 15 figures. Submitted to ApJ. Comments welcome!

Galaxy clustering in modified gravity from full-physics simulations. I: two-point correlation functions

Authors: Michael Collier, Sownak Bose, Baojiu Li

Abstract: We present an in-depth investigation of galaxy clustering based on a new suite of realistic large-box galaxy-formation simulations in $f(R)$ gravity, with a subgrid physics model that has been recalibrated to reproduce various observed stellar and gas properties. We focus on the two-point correlation functions of the luminous red galaxies (LRGs) and emission line galaxies (ELGs), which are primary… ▽ More We present an in-depth investigation of galaxy clustering based on a new suite of realistic large-box galaxy-formation simulations in $f(R)$ gravity, with a subgrid physics model that has been recalibrated to reproduce various observed stellar and gas properties. We focus on the two-point correlation functions of the luminous red galaxies (LRGs) and emission line galaxies (ELGs), which are primary targets of ongoing and future galaxy surveys such as DESI. One surprising result is that, due to several nontrivial effects of modified gravity on matter clustering and the galaxy-halo connection, the clustering signal does not depend monotonically on the fifth-force strength. For LRGs this complicated behaviour poses a challenge to meaningfully constraining this model. For ELGs, in contrast, this can be straightforwardly explained by the time evolution of the fifth force, which means that weaker $f(R)$ models can display nearly the same -- up to $25\%$ -- deviations from $Λ$CDM as the strongest ones, albeit at lower redshifts. This implies that even very weak $f(R)$ models can be strongly constrained, unlike with most other observations. Our results show that galaxy formation acquires a significant environment dependence in $f(R)$ gravity which, if not properly accounted for, may lead to biased constraints on the model. This highlights the essential role of hydrodynamical simulations in future tests of gravity exploring precision galaxy-clustering data from the likes of DESI and Euclid. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Comments: 18 pages, 8 figures + 2 figures in appendix

Spectral evolution of RX J0440.9+4431 during the 2022-2023 giant outburst observed with Insight-HXMT

Authors: P. P. Li, Peter A. Becker, L. Tao

Abstract: In 2022-2023, the X-ray pulsar RX J0440.9+4431 underwent a Type II giant outburst, reaching a peak luminosity L_x ~ 4*10^{37} erg/s. In this work, we utilize Insight-HXMT data to analyze the spectral evolution of RX J0440.9+4431 during the giant outburst. By analysing the variation of the X-ray spectrum during the outburst using standard phenomenological models, we find that as the luminosity appr… ▽ More In 2022-2023, the X-ray pulsar RX J0440.9+4431 underwent a Type II giant outburst, reaching a peak luminosity L_x ~ 4*10^{37} erg/s. In this work, we utilize Insight-HXMT data to analyze the spectral evolution of RX J0440.9+4431 during the giant outburst. By analysing the variation of the X-ray spectrum during the outburst using standard phenomenological models, we find that as the luminosity approaches the critical luminosity, the spectrum became flatter, with the photon enhancement predominantly concentrated around ~ 2 keV and 20-40 keV. The same behavior has also been noted in Type II outbursts from other sources. While the phenomenological models provide good fits to the spectrum, this approach is sometimes difficult to translate into direct insight into the details of the fundamental accretion physics. Hence we have also analyzed spectra obtained during high and low phases of the outburst using a new, recently-developed physics-based theoretical model, which allows us to study the variations of physical parameters such as temperature, density, and magnetic field strength during the outburst. Application of the theoretical model reveals that the observed spectrum is dominated by Comptonized bremsstrahlung emission emitted from the column walls in both the high and low states. We show that the spectral flattening observed at high luminosities results from a decrease in the electron temperature, combined with a compactification of the emission zone, which reduces the efficiency of bulk Comptonization. We also demonstrate that when the source is at maximum luminosity, the spectrum tends to harden around the peak of the pulse profile, and we discuss possible theoretical explanations for this behavior. We argue that the totality of the behavior in this source can be explained if the accretion column is in a quasi-critical state when at the maximum luminosity observed during the outburst. △ Less

Submitted 21 June, 2024; originally announced July 2024.

Comments: 11 pages, 4 figures, accepted for publication in A&A

Instrument Overview of Taurus: A Balloon-borne CMB and Dust Polarization Experiment

Authors: Jared L. May, Alexandre E. Adler, Jason E. Austermann, Steven J. Benton, Rick Bihary, Malcolm Durkin, Shannon M. Duff, Jeffrey P. Filippini, Aurelien A. Fraisse, Thomas J. L. J. Gascard, Sho M. Gibbs, Suren Gourapura, Jon E. Gudmundsson, John W. Hartley, Johannes Hubmayr, William C. Jones, Steven Li, Johanna M. Nagy, Kate Okun, Ivan L. Padilla, L. Javier Romualdez, Simon Tartakovsky, Michael R. Vissers

Abstract: Taurus is a balloon-borne cosmic microwave background (CMB) experiment optimized to map the E-mode polarization and Galactic foregrounds at the largest angular scales ($\ell$ $\lt$ 30) and improve measurements of the optical depth to reionization ($τ$). This will pave the way for improved measurements of the sum of neutrino masses in combination with high-resolution CMB data while also testing the… ▽ More Taurus is a balloon-borne cosmic microwave background (CMB) experiment optimized to map the E-mode polarization and Galactic foregrounds at the largest angular scales ($\ell$ $\lt$ 30) and improve measurements of the optical depth to reionization ($τ$). This will pave the way for improved measurements of the sum of neutrino masses in combination with high-resolution CMB data while also testing the $ΛCDM$ model on large angular scales and providing high-frequency maps of polarized dust foregrounds to the CMB community. These measurements take advantage of the low-loading environment found in the stratosphere and are enabled by NASA's super-pressure balloon platform, which provides access to 70% of the sky with a launch from Wanaka, New Zealand. Here we describe a general overview of Taurus, with an emphasis on the instrument design. Taurus will employ more than 10,000 100 mK transition edge sensor bolometers distributed across two low-frequency (150, 220 GHz) and one high-frequency (280, 350 GHz) dichroic receivers. The liquid helium cryostat housing the detectors and optics is supported by a lightweight gondola. The payload is designed to meet the challenges in mass, power, and thermal control posed by the super-pressure platform. The instrument and scan strategy are optimized for rigorous control of instrumental systematics, enabling high-fidelity linear polarization measurements on the largest angular scales. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Preliminary results of sky brightness measurements in near-infrared at Lenghu, China

Authors: Jinji Li, Bin Ma, Zhongnan Dong, Haoran Zhang

Abstract: Low sky brightness is crucial for ground-based astronomical observations, because it limits the observational capability to detect fainter sources. Lenghu, located on the Tibetan Plateau in China, has been identified as an high-quality astronomical site in China, including dark sky in optical band. In this work, we will report the preliminary results of near-infrared sky brightness measurements at… ▽ More Low sky brightness is crucial for ground-based astronomical observations, because it limits the observational capability to detect fainter sources. Lenghu, located on the Tibetan Plateau in China, has been identified as an high-quality astronomical site in China, including dark sky in optical band. In this work, we will report the preliminary results of near-infrared sky brightness measurements at Lenghu. Utilizing a wide-field small telescope equipped with an InGaAs camera, we have been conducting long-term monitoring of near-infrared sky brightness in the J and H' bands, respectively, since January 2024. For each image, photometry and astrometry were performed, then sky background was calibrated by standard stars from the 2MASS catalog. This report includes preliminary results on the sky brightness at zenith in the J and H' bands, as well as their variations with solar elevation at Lenghu. Our initial results indicate that the near-infrared sky brightness at Lenghu is comparable to that of other world-class sites, and long-term monitoring will be continued. △ Less

Submitted 8 July, 2024; v1 submitted 1 July, 2024; originally announced July 2024.

Comments: 8 pages, 8 figures, presented at 2024 SPIE Astronomical Telescopes + Instrumentation conference

Direct observational evidence of multi-epoch massive star formation in G24.47+0.49

Authors: Anindya Saha, Anandmayee Tej, Hong-Li Liu, Tie Liu, Guido Garay, Paul F. Goldsmith, Chang Won Lee, Jinhua He, Mika Juvela, Leonardo Bronfman, Tapas Baug, Enrique Vazquez-Semadeni, Patricio Sanhueza, Shanghuo Li, James O. Chibueze, N. K. Bhadari, Lokesh K. Dewangan, Swagat Ranjan Das, Feng-Wei Xu, Namitha Issac, Jihye Hwang, L. Viktor Toth

Abstract: Using new continuum and molecular line data from the ALMA Three-millimeter Observations of Massive Star-forming Regions (ATOMS) survey and archival VLA, 4.86 GHz data, we present direct observational evidence of hierarchical triggering relating three epochs of massive star formation in a ring-like H II region, G24.47+0.49. We find from radio flux analysis that it is excited by a massive star(s) of… ▽ More Using new continuum and molecular line data from the ALMA Three-millimeter Observations of Massive Star-forming Regions (ATOMS) survey and archival VLA, 4.86 GHz data, we present direct observational evidence of hierarchical triggering relating three epochs of massive star formation in a ring-like H II region, G24.47+0.49. We find from radio flux analysis that it is excited by a massive star(s) of spectral type O8.5V-O8V from the first epoch of star formation. The swept-up ionized ring structure shows evidence of secondary collapse, and within this ring a burst of massive star formation is observed in different evolutionary phases, which constitutes the second epoch. ATOMS spectral line (e.g., HCO$^+$(1-0)) observations reveal an outer concentric molecular gas ring expanding at a velocity of $\sim$ 9 $\rm km\,s^{-1}$, constituting the direct and unambiguous detection of an expanding molecular ring. It harbors twelve dense molecular cores with surface mass density greater than 0.05 $\rm g\,cm^{-2}$, a threshold typical of massive star formation. Half of them are found to be subvirial, and thus in gravitational collapse, making them third epoch of potential massive star-forming sites. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Comments: 18 pages, 7 figures, accepted for publication in The Astrophysical Journal Letters

Evidence for particle acceleration approaching PeV energies in the W51 complex

Authors: LHAASO Collaboration

Abstract: The $γ$-ray emission from the W51 complex is widely acknowledged to be attributed to the interaction between the cosmic rays (CRs) accelerated by the shock of supernova remnant (SNR) W51C and the dense molecular clouds in the adjacent star-forming region, W51B. However, the maximum acceleration capability of W51C for CRs remains elusive. Based on observations conducted with the Large High Altitude… ▽ More The $γ$-ray emission from the W51 complex is widely acknowledged to be attributed to the interaction between the cosmic rays (CRs) accelerated by the shock of supernova remnant (SNR) W51C and the dense molecular clouds in the adjacent star-forming region, W51B. However, the maximum acceleration capability of W51C for CRs remains elusive. Based on observations conducted with the Large High Altitude Air Shower Observatory (LHAASO), we report a significant detection of $γ$ rays emanating from the W51 complex, with energies from 2 TeV to 200 TeV. The LHAASO measurements, for the first time, extend the $γ$-ray emission from the W51 complex beyond 100 TeV and reveal a significant spectrum bending at tens of TeV. By combining the ``$π^0$-decay bump" featured data from Fermi-LAT, the broadband $γ$-ray spectrum of the W51 region can be well-characterized by a simple pp-collision model. The observed spectral bending feature suggests an exponential cutoff at $\sim400$~TeV or a power-law break at $\sim200$~TeV in the CR proton spectrum, most likely providing the first evidence of SNRs serving as CR accelerators approaching the PeV regime. Additionally, two young star clusters within W51B could also be theoretically viable to produce the most energetic $γ$ rays observed by LHAASO. Our findings strongly support the presence of extreme CR accelerators within the W51 complex and provide new insights into the origin of Galactic CRs. △ Less

Submitted 30 June, 2024; originally announced July 2024.

Accurate Shear Recovery with Multi-Band Images of Hyper Suprime-Cam

Authors: Cong Liu, Jun Zhang, Hekun Li, Pedro Alonso Vaquero, Wenting Wang

Abstract: The existing large scale weak lensing surveys typically reserve the best seeing conditions for a certain optical band to minimize shape measurement errors and maximize the number of usable background galaxies. This is because most popular shear measurement methods contain explicit or implicit thresholds on the galaxy-to-PSF (point spread function) size ratio, below which their shape measurement er… ▽ More The existing large scale weak lensing surveys typically reserve the best seeing conditions for a certain optical band to minimize shape measurement errors and maximize the number of usable background galaxies. This is because most popular shear measurement methods contain explicit or implicit thresholds on the galaxy-to-PSF (point spread function) size ratio, below which their shape measurement errors increase abruptly. Using the DECaLS data, we have previously demonstrated that the Fourier\_Quad method performs very well on poorly resolved galaxy images in general. It is therefore a ready tool for shear measurement with multi-band images regardless of their seeing conditions. In this paper, we apply the Fourier\_Quad pipeline on the multi-band images from the third public data release of the Hyper Suprime-Cam Subaru Strategic Program. We show that the shear catalogs from the five optical bands (g/r/i/z/y) all pass the field-distortion test with very high accuracy. Using the LOWZ and CMASS galaxies as foreground lenses, we show that the errorbar in the galaxy-galaxy lensing measurement can be decreased by factors around 15\% by combining shear catalogs from different bands. This indicates that it is worthful to do multi-bands shear measurements for a better shear statistics. △ Less

Submitted 29 June, 2024; originally announced July 2024.

Comments: 14 pages, 12 figures

A non-accreting neutron star and $δ$ Scuti candidate revealed by asteroseismological modelling and radial velocity measurements

Authors: Ping Li

Abstract: SZ Lyn is an astrometric binary consisting of a $δ$ Scuti star and an unseen component. The two stars around each other in an almost circular orbit (e=0.186) with a period of about 3.32 years. We aim to explore the physical properties of the two components using the radial velocity method and asteroseismology. Our results show that the $δ$ Scuti variable is a post-main sequence star with a helium… ▽ More SZ Lyn is an astrometric binary consisting of a $δ$ Scuti star and an unseen component. The two stars around each other in an almost circular orbit (e=0.186) with a period of about 3.32 years. We aim to explore the physical properties of the two components using the radial velocity method and asteroseismology. Our results show that the $δ$ Scuti variable is a post-main sequence star with a helium core ($M_{\rm He}$= 0.1671$^{+0.0028}_{-0.0045}$ M$_{\odot}$ and $R_{\rm He}$= 0.1414$^{+0.0035}_{-0.0020}$ R$_{\odot}$) surrounded by a burning hydrogen shell; its fundamental parameters is accurately calculated to be $M$=1.82$^{+0.04}_{-0.02}$ M$_{\odot}$, $R$=2.91$^{+0.02}_{-0.02}$ R$_{\odot}$, $L$=17.07$^{+1.28}_{-2.61}$ L$_{\odot}$, and $X_c$=0.078$^{+0.021}_{-0.004}$; the mass of the invisible object is thus estimated to be $M_{\rm co}$=1.68 $^{+0.62}_{-0.60}$ M$_{\odot}$. Combined with the low-resolution LAMOST spectrum, we suggest that the unseen object is an neutron star candidate. Our work demonstrates the methodology of using the radial velocity method and asteroseismology to detect the compact object in such an astrometric binary. △ Less

Submitted 29 June, 2024; originally announced July 2024.

Comments: In preparation

FAST survey of H I and OH absorption towards extragalactic radio sources

Authors: Yogesh Chandola, D. J. Saikia, Yin-Zhe Ma, Zheng Zheng, Chao-Wei Tsai, Di Li, Denis Tramonte, Hengxing Pan

Abstract: Neutral atomic hydrogen and molecular gas in the host galaxies of radio active galactic nuclei (AGN) can be traced using H I 21-cm and OH-1667 MHz absorption lines to understand the fueling and feedback processes. We present the results of an H I and OH absorption survey with the Five-hundred-meter Aperture Spherical radio Telescope (FAST) towards 40 radio sources of low-intermediate radio luminos… ▽ More Neutral atomic hydrogen and molecular gas in the host galaxies of radio active galactic nuclei (AGN) can be traced using H I 21-cm and OH-1667 MHz absorption lines to understand the fueling and feedback processes. We present the results of an H I and OH absorption survey with the Five-hundred-meter Aperture Spherical radio Telescope (FAST) towards 40 radio sources of low-intermediate radio luminosity ($\sim$10$^{23}$-10$^{26}$ W Hz$^{-1}$ at 1.4 GHz), red mid-infrared color (W2[4.6 $μ$m]$-$W3[12 $μ$m] $>$ 2.5 mag) and redshift up to 0.35. From 13 sources with good data at H I observing frequencies, we report the detection of H I absorption towards 8 sources, 5 of which are new detections including 4 in the redshift range 0.25 to 0.35. Our detection rates are consistent with our previous results with dependence on the star-formation history of the host galaxy reflected in the mid-infrared \textit{WISE} W2$-$W3 colors and the compactness of the radio source. We find no significant dependence of detection rates on radio luminosity or redshift. We also find that H I column densities are anti-correlated with the low-frequency spectral indices ($α_{\rm 150 MHz}^{\rm 1.4 GHz}$, $S_ν\propto ν^{-α}$). We do not have any detection from 23 sources with good data at OH observing frequencies. However, by stacking the spectra we estimate the 3$σ$ upper limit of OH column density to be 2.27$\times$10$^{14}$$T_{\rm ex}$/10 K $\times$1/$f_{\rm c}$ cm$^{-2}$. By stacking the OH spectra for 7 associated H I absorbers, we get a 3$σ$ upper limit of 3.47$\times$10$^{14}$ $T_{\rm ex}$/10 K $\times$1/$f_{\rm c}$ cm$^{-2}$ on OH column density and 1.78$\times$10$^{-7}$ on [OH]/[H I] ratio. △ Less

Submitted 28 June, 2024; originally announced June 2024.

Comments: 20 pages, 8 figures, accepted for publication in The Astrophysical Journal (ApJ)

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

The impact of shear on the rotation of Galactic plane molecular clouds

Authors: Raffaele Rani, Jia-Lun Li, Toby J. T. Moore, David J. Eden, Andrew J. Rigby, Geumsook Park, Yueh-Ning Lee

Abstract: Stars form in the densest regions of molecular clouds, however, there is no universal understanding of the factors that regulate cloud dynamics and their influence on the gas-to-stars conversion. This study considers the impact of Galactic shear on the rotation of giant molecular clouds (GMCs) and its relation to the solenoidal modes of turbulence. We estimate the direction of rotation for a large… ▽ More Stars form in the densest regions of molecular clouds, however, there is no universal understanding of the factors that regulate cloud dynamics and their influence on the gas-to-stars conversion. This study considers the impact of Galactic shear on the rotation of giant molecular clouds (GMCs) and its relation to the solenoidal modes of turbulence. We estimate the direction of rotation for a large sample of clouds in the \ce{^{13}CO}/\ce{C^{18}O} (3-2) Heterodyne Inner Milky Way Plane Survey (CHIMPS) and their corresponding sources in a new segmentation of the \ce{^{12}CO}(3-2) High-Resolution Survey (COHRS). To quantify the strength of shear, we introduce a parameter that describes the shear's ability to disrupt growing density perturbations within the cloud. Although we find no correlation between the direction of cloud rotation, the shear parameter, and the magnitude of the velocity gradient, the solenoidal fraction of the turbulence in the CHIMPS sample is positively correlated with the shear parameter and behaves similarly when plotted over Galactocentric distance. GMCs may thus not be large or long-lived enough to be affected by shear to the point of showing rotational alignment. In theory, Galactic shear can facilitate the rise of solenoidal turbulence and thus contribute to suppressing star formation. These results also suggest that the rotation of clouds is not strictly related to the overall rotation of the disc, but is more likely to be the imprint of Kelvin-Helmholtz instabilities in the colliding flows that formed the clouds. △ Less

Submitted 27 June, 2024; originally announced June 2024.

Comments: Accepted, MNRAS

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

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

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

Submitted 27 June, 2024; originally announced June 2024.

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

A population of mid-infrared large-amplitude variable young stellar objects from unTimely

Authors: Jiaxun Li, Tinggui Wang

Abstract: Utilizing a decade-long unTimely dataset, supplemented by multi-band data from archives, we search for young stellar objects (YSOs) with variations larger than one magnitude in W1 band within a region of 110 square degrees in the Galactic plane, covered by VISTA Variables in the Via Lactea (VVV). A total of 641 candidate YSOs have been identified. We classified them into bursts, dips, faders, secu… ▽ More Utilizing a decade-long unTimely dataset, supplemented by multi-band data from archives, we search for young stellar objects (YSOs) with variations larger than one magnitude in W1 band within a region of 110 square degrees in the Galactic plane, covered by VISTA Variables in the Via Lactea (VVV). A total of 641 candidate YSOs have been identified. We classified them into bursts, dips, faders, seculars, and irregulars. Within the burst category, 18 sources were identified as FUor candidates and 1 as an EXor candidate. Irregulars are the most prevalent in the sample. In both bursts and faders, the redder sources tend to show a pattern of bluer when brighter, whereas the bluer sources display the opposite trend, possibly related to the accretion structure of YSOs at different stages. Finally, we obtained the recurrence time scale for FUor eruptions at various stages of YSO evolution. Our findings indicate that younger YSOs generally experience more frequent eruptions compared to older ones. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: 22 pages, 20 figures