2025-12-30 | | Total: 73
We introduce Galaxy Zoo Evo, a labeled dataset for building and evaluating foundation models on images of galaxies. GZ Evo includes 104M crowdsourced labels for 823k images from four telescopes. Each image is labeled with a series of fine-grained questions and answers (e.g. "featured galaxy, two spiral arms, tightly wound, merging with another galaxy"). These detailed labels are useful for pretraining or finetuning. We also include four smaller sets of labels (167k galaxies in total) for downstream tasks of specific interest to astronomers, including finding strong lenses and describing galaxies from the new space telescope Euclid. We hope GZ Evo will serve as a real-world benchmark for computer vision topics such as domain adaption (from terrestrial to astronomical, or between telescopes) or learning under uncertainty from crowdsourced labels. We also hope it will support a new generation of foundation models for astronomy; such models will be critical to future astronomers seeking to better understand our universe.
Over the past two decades, the \textit{Swift} and \textit{Fermi} missions have identified a rare class of ``double-trigger'' gamma-ray bursts (GRBs) that produce two independent trigger events. These events are characterized by a sufficiently long quiescent period during which the on-board trigger system can reset, resulting in the subsequent emission being recorded as a second independent event. Consistent sky localization confirms that both trigger events originated from the same astrophysical source. Here, we present a systematic classification and characteristics study of three such cases: GRB 091024A, GRB 110709B, and GRB 220627A. We investigate each trigger episode emission independently using standard classification diagnostics, including duration ($T_{90}$), hardness ratio, minimum variability timescale (MVT), spectral lag ($τ_{\rm lag}$), peak energy ($E_{\rm p}$), and energetics. We compare these properties with those of typical long GRBs (LGRBs) and with single-trigger LGRBs that exhibit extended quiescent periods. Our analysis reveals that all sub-bursts from the three double-trigger events consistently lie within the LGRB classification region. These results indicate that double-trigger GRBs are not a physically distinct subclass, but rather products of LGRB central engines that undergo extended dormancy and subsequent reactivation.
This work analyzes four Sun-like double-lined spectroscopic binary (SB2) systems by combining visual and spectroscopic observational data with Al-Wardat's atmospheric modeling method to accurately determine their fundamental parameters. For each system, we determine stellar masses, orbital parallaxes, effective temperatures, spectral types, semimajor axes, and eccentricities with high precision, resolving discrepancies between astrometric and spectroscopic measurements. Moreover, we assess the potential for stable planetary orbits in these systems. We also calculate habitable zones around these binaries based on the orbital evolution of planetary orbits. These systems may represent promising targets for future extrasolar planet searches around Sun-like stars due to their robust physical and orbital parameters that can be used to determine planetary habitability and stability.
As its periastron passage occurred during the third quarter of 2020, system 24 Aqr is of particular significance. New visual solutions for the latest speckle interferometry observations collected by the Lowell Discovery Telescope (LTD) with its new QWSSI speckle camera are presented here. A variety of techniques were used to analyze the system, including ORBITX code for orbital solution, Al-Wardat's method for analyzing multiple stellar systems, and Edwards' method for analyzing visual and spectroscopic binaries. We derive precise masses and the complete set of its fundamental parameters for the three components, and we introduce a new orbital solution, and a new dynamical parallax, which is very close to the measured value given by Hipparcos 2007 and from that of Gaia DR2. In the next section, we discuss the possibility of a coplanar orbit. In conclusion, we demonstrate that we need a 65-m telescope to resolve the inner binary visually, although an array of telescopes could be used instead.
Dark matter admixed neutron stars (DANSs) serve as a specific astrophysical laboratory for probing the features of dark matter (DM) and have emerged as a promising candidate for interpreting recent astrophysical observations (e.g., by NICER and LIGO/Virgo). Accurately constraining the internal DM content of DANSs is therefore of critical importance. In this work, we construct the equations of state (EoS) for DANS matter by employing twelve nuclear matter (NM) models within the covariant density functional (CDF) theory and a self-interacting fermionic model for DM. Using these EoSs as input, we solve the two-fluid Tolman-Oppenheimer-Volkov (TOV) equations to systematically investigate the influence of DM on the global properties of neutron stars (NSs). By incorporating recent observational constraints on NS properties, the maximum DM mass fraction $f_χ^{\mathrm{max}}$ in DANSs is determined for each NM EoS model. Our analysis reveals a strong linear correlation (Pearson coefficient $r=0.98$) between $f_χ^{\mathrm{max}}$ and the maximum mass of a pure NS, $M_{\rm{NS}}^{\mathrm{max}}$, described by $f_χ^{\mathrm{max}} = 0.22 M_{\mathrm{NS}}^{\mathrm{max}} - 0.44$. Leveraging this correlation and the observed NS maximum mass distribution, $P(M_{\text{NS}}^{\max} \mid \text{EM})$, we derive the probability distribution function (PDF) for the maximum DM mass, $P(M_χ^{\max} \mid \text{EM})$, in DANSs. We find that at the 68\% confidence level, $M_χ^{\mathrm{max}}=0.150^{+0.070}_{-0.051}\ M_{\odot}$. This quantitative constraint on the DM mass provides a critical prior for interpreting potential observational signatures of DANSs, such as anomalous tidal deformabilities and distinctive gravitational-wave signals.
Surviving rapid inward orbital migration is a crucial aspect of formation models for the Jupiter's Galilean moons. The primary aim of this study is to investigate the orbital migration of the Galilean moons by incorporating self-consistent solid dynamics in circumjovian disk models. We perform two-fluid simulations using the FARGO3D code on a 2D polar grid. The simulations model a satellite with the mass of a proto-moon, Europa, or Ganymede interacting with a circumjovian disk. The dust component, coupled to the gas via a drag force, is characterized by the dust-to-gas mass ratio ($ε$) and the Stokes number ($T_s$). The effect of solids fundamentally alter the satellites' evolution. We identify a vast parameter space where migration is slowed, halted, robustly reversed -leading to outward migration-, or significantly accelerated inward. The migration rate is dependent on satellite mass, providing a natural source of differential migration. Solid dynamics provides a robust and self-consistent mechanism that fundamentally alters the migration of the Galilean moons, potentially addressing the long-standing migration catastrophe. This mechanism critically affects the survival of satellites and could offer a viable physical process to explain the establishment of resonances through differential migration. These findings establish that solid torques are a critical, non-negligible factor in shaping the final architecture of satellite systems.
In past work, we described the use of a reconfigurable intelligent surface (RIS) mounted on the rim of an axisymmetric prime focus-fed reflector to create nulls in the close-in sidelobes. In this paper, we show that similar performance is possible in an offset Gregorian reflector system using a RIS on the rim of the subreflector. Applications include radio astronomy, where offset Gregorian reflectors are common and observations are subject to deleterious levels of interference from satellites entering through sidelobes. We show that an efficient RIS replacing the outer one-third of the subreflector surface, employing passive elements with 1-bit phase-only control, can create a null in the peak of the second sidelobe in the quiescent pattern. This is achieved using a simple unconstrained optimization algorithm to set the states of the RIS elements. The algorithm yields a deep null with just 0.2~dB reduction in main lobe directivity, despite lacking any constraints on main lobe pattern. Compared to our previous approach of mounting the RIS on the rim of the main reflector, the subreflector-based approach demonstrated in this paper requires a much smaller RIS and can implemented in existing systems by replacing the subreflector.
Context. Processing radio interferometric data often requires storing forward-predicted model data. In direction-dependent calibration, these data may have a volume an order of magnitude larger than the original data. Existing lossy compression techniques work well for observed, noisy data, but cause issues in calibration when applied to forward-predicted model data. Aims. To reduce the volume of forward-predicted model data, we present a lossless compression method called Simulated Signal Compression (Sisco) for noiseless data that integrates seamlessly with existing workflows. We show that Sisco can be combined with baseline-dependent averaging for further size reduction. Methods. Sisco decomposes complex floating-point visibility values and uses polynomial extrapolation in time and frequency to predict values, groups bytes for efficient encoding, and compresses residuals using the Deflate algorithm. We evaluate Sisco on diverse LOFAR, MeerKAT, and MWA datasets with various extrapolation functions. Implemented as an open-source Casacore storage manager, it can directly be used by any observatory that makes use of this format. Results. We find that a combination of linear and quadratic prediction yields optimal compression, reducing noiseless forward-predicted model data to 24% of its original volume on average. Compression varies by dataset, ranging from 13% for smooth data to 38% for less predictable data. For pure noise data, compression achieves just a size of 84% due to the unpredictability of such data. With the current implementation, the achieved compression throughput is with 534 MB/s mostly dominated by I/O on our testing platform, but occupies the processor during compression or decompression. Finally, we discuss the extension to a lossy algorithm.
The DArk Matter Particle Explorer (DAMPE) is a space-borne high-energy particle detector that surveys the $γ$-ray sky above$\sim 2~\rm GeV$ with a peak acceptance of $\sim 0.2~\rm m^2\,sr$. With the 102 months of data collected by DAMPE, we show that the Fermi bubbles are detected at a significance of $\sim 26σ$ and identify a GeV excess in the direction of Galactic center at $\sim 7 σ$ confidence. Both spectra and morphology are consistent with those observed by Fermi-LAT and the GeV excess component can be interpreted by the dark matter annihilation with a mass of $\sim 50$ GeV and a velocity-averaged cross section of $\sim 10^{-26}~{\rm cm^{3}~s^{-1}}$ for the $χχ\rightarrow b\bar{b}$ channel. Our results thus provide the first independent detection of these two intriguing diffuse gamma-ray sources besides Fermi-LAT.
High-mass protostars are deeply embedded in dust inside their natal cores and are not easily detectable. However, maser emission at centimeter wavelengths, owing to its high brightness, enables us to study gas kinematics in protostars' circumstellar regions. We aim to understand the origin of the ring-like structures outlined by the 6.7 GHz methanol maser emission in six high-mass young stellar objects by performing a sensitive search of the associated radio-continuum emission and derive its properties. We used the Karl G. Jansky Very Large Array in the A configuration at C and K bands in order to image radio-continuum as well as 6.7 GHz methanol and 22 GHz water maser emission. We present the first images of the thermal jets towards four targets in our sample, G23.389+00.185, G23.657-00.127, G28.817+00.365, and G30.400-00.296. In a further target, G23.207-00.377, the complex K band continuum emission makes it unclear whether the detected peaks trace jet knots from a single young protostar or mark multiple compact young protostars. The remaining source G31.047+00.356 shows radio continuum emission associated with an evolved H II region.
Using Gaia DR3 data, binary star catalogs have been created containing information on a total of more than 2.6 million pairs. This increases by more than an order of magnitude the ensemble of binary stars with known characteristics, which previously numbered about 140 thousand pairs. To perform statistical analysis of the complete ensemble of binary stars, including both previously known and newly discovered pairs, cross-identification by coordinates was carried out between the most complete pre-Gaia publication compilative binary star catalog ILB and data from binary star catalogs based on Gaia DR3 results. An analysis of the results of this identification was performed, showing the dependence of its characteristics both on the data from the source catalogs and on coordinates. It is shown that in dense stellar fields, particularly in the Galactic disk, an increase in the fraction of false positive identifications can be expected. At the same time, for systems with large proper motion, there is a high probability of a false negative outcome. Possible modifications to the identification method are proposed to reduce the role of the described systematic errors and increase the reliability of its results.
We present comprehensive photometric and spectroscopic observations of Supernova (SN) 2022eyw, a luminous member of the Type Iax SN subclass. SN 2022eyw reached a peak absolute magnitude of $M_g = -17.80\pm0.15$ mag and exhibited a rise time of $\sim$15 days, placing it among the brighter Iax events. The bolometric light curve indicates a synthesized $^{56}$Ni mass of $0.120\pm0.003~\text{M}_{\odot}$, with an estimated ejecta mass of $0.79\pm0.09~\text{M}_{\odot}$ and kinetic energy of $0.19\times10^{51}$ erg. The spectral evolution from -8 to +110 days past maximum reveals features characteristic of bright Type Iax Supernovae, including a transition from Fe III to Fe II dominance, moderate expansion velocities, and a lack of strong C III absorption. TARDIS spectral modelling of the early-phase spectra indicates a well-mixed ejecta dominated by Fe-group elements. In addition, traces of unburnt carbon are detected, pointing to incomplete burning as expected in pure deflagration models. Late-time spectral evolution shows a blend of permitted and forbidden lines. Comparison with deflagration models suggests that SN 2022eyw originated from a partial deflagration of a Chandrasekhar-mass white dwarf, with explosion properties intermediate between the N3-def and N5-def models. These observations support pure deflagration of a CO white dwarf as a viable explosion mechanism for its luminous members.
This paper constructs and analyzes a three channel dissipative framework for Warm Higgs Inflation, wherein the total dissipation coefficient, $Υ(h,T)$, is decomposed into low temperature, high temperature, and threshold activated contributions. A genetic algorithm is employed for the global numerical solution and statistical inference of the background field dynamics. To overcome the single channel dominance degeneracy in high dimensional parameter scans, two classes of structural priors are introduced into the objective function: a \texttt{mixing} prior to suppress extreme channel fractions and an \texttt{entropy} prior to favor multi channel coexistence. Furthermore, the adoption of a layered warmness criterion (e.g., $Q \equiv Υ/3H$) decouples model selection from cosmological observables, thereby enhancing analytical transparency. The complete workflow is demonstrated on a $14$ dimensional phenomenological model. An ablation study of the priors (\texttt{noprior} vs. \texttt{mixing} vs. \texttt{mixing+entropy}) yields $18871$ viable parameter points, revealing that the priors significantly enhance the discovery probability of non-trivial multi channel solutions within a parameter space naturally biased towards single channel dominance. After imposing the observational constraint $r < 0.036$, the number of retained solutions for each scenario is $14485$, $1889$, and $1971$, respectively. A typical best fit solution exhibits a "channel relay" dynamical feature during its evolution and a genuinely mixed state at the pivot scale (e.g., $f_{{\rm HT},*}\simeq 0.399480$, $f_{{\rm th},*}\simeq 0.600517$), implying that the microscopic origin of dissipation need not be unique within a single inflationary history.
In this work, we report evidence suggesting the potential future detection of a month-scale quasi-periodic oscillation (QPO) in the gamma-ray light curve of OP 313. We analysed almost 16.8 years of Fermi-LAT gamma-ray data and applied the Bayesian block method to the monthly-binned light curve. We identified four high-flux states and investigated the possibility of a QPO in the fourth high-flux state (MJD 59482-60832). Using the Weighted Wavelet Z-transform (WWZ) and Lomb-Scargle Periodogram (LSP) methods, we find tentative evidence for a month-scale QPO; however, its detection significance is limited by the small number of observed cycles. With a sufficiently long data set, the QPO may be detected with higher significance in the future. We further explored possible physical origins of this potential QPO and examined several models. We found that a curved-jet model can explain the observed behaviour.
Fast radio bursts (FRBs) are millisecond-duration radio transients whose phys- ical origin remains uncertain. Magnetar-based models, motivated by observed properties such as polarization and large rotation measures, suggest that FRB emission may be modulated by the magnetar spin period. We present an efficient method to search for periodic signals in repeating FRBs by combining phase folding and Markov Chain Monte Carlo (MCMC) parameter estimation. Our method accelerates period searches. We test the method using observational data from repeater FRB 20201124A, and show that it can recover reported candidate periods.
Filamentary infrared dark clouds (IRDCs) are believed to represent the initial conditions for massive star and cluster formation. We investigate the IRDC G035.39-00.33 using SiO, H13CO+, CH3OH, and CS emission observed with ALMA at 3.5\arcsec\ resolution (0.05 pc). The SiO emission traces shock activity within the cloud, providing insights into current star formation and cloud formation mechanisms. We identify several regions with broad SiO emission clearly associated with outflows, pinpointing the locations of ongoing star formation across the cloud. The ALMA images also reveal a series of spatially extended SiO emission spots with narrow line profiles, aligned along an arc-like path that is also seen in CS and CH3OH emission. While the broad SiO emission is mainly associated with the main cloud filament, as seen in visual extinction, the narrow SiO arch is located at the edge of the cloud, far from the identified sites of star formation activity. The presence of these arc-like morphologies suggests that large-scale shocks may have compressed the gas in the surroundings of the G035.39-00.33 cloud, shaping its filamentary structure. By inspecting the large-scale radio continuum emission around G035.39-00.33, we find that this IRDC is part of a larger star-forming complex where the densest and coolest material appears at the interacting regions between a Supernova Remnant (SNR) and an expanding HII region. In particular, we hypothesize that this IRDC may be spatially coincident with the ionized expanding gas associated with the previously identified SNR G35.6-0.4. We suggest that collisions between giant molecular clouds and expanding gas flows from interacting SNRs and HII regions may be responsible for the observed arc-like structures. Such shock compressions could play an important role in the formation of IRDCs and in the potential triggering of star formation.
Chinese Space Station Telescope (CSST), which will begin its scientific operations around 2027, is going to survey the sky area of the median-to-high Galactic latitude and median-to-high ecliptic latitude. The high astrometric precision of the CSST Survey Camera for faint objects enables the detection of a number of giant planets and brown dwarfs around M-dwarfs and brown dwarfs via differential astrometry in its optical survey. In this paper, we predict the number of giant planets and brown dwarfs around stars and brown dwarfs detectable with CSST astrometry. We generate synthetic samples of CSST stellar and substellar sources, and carry out companion injection-recovery simulations in the samples using different occurrence rates for FGK-dwarfs, M-dwarfs, and brown dwarfs. We calculate companion yields based on CSST astrometric precision. Our analysis reveals that over its 10-year mission, the CSST Survey Camera could barely discover giant planets and low-mass BDs around FGK-dwarfs, but is projected to detect 20 - 170 giant planets and low-mass brown dwarfs around M-dwarfs within 300 pc, and 300 - 570 brown dwarf binaries within 600 pc. Therefore, CSST astrometry is likely to significantly increase the current sample of substellar companions around M-dwarfs and brown dwarfs. This sample will deepen our understanding of planet formation and evolution around low-mass stars and brown dwarfs.
Gamma-Ray Burst (GRB) prompt and afterglow emission, as well as a kilonova (KN), are the expected electromagnetic (EM) counterparts of Binary Neutron Star (BNS) and Neutron Star -- Black Hole (NSBH) mergers. We aim to infer the KN ejecta parameters and the progenitor properties by modeling merger-driven GRBs with a claim of KN, good data and robust redshift measurement. We model the afterglow and KN, and perform a Bayesian analysis, within the Nuclear physics and Multi-Messenger Astrophysics (NMMA) framework. The KN emission is modeled with the radiative transfer code POSSIS and for afterglow we use the afterglowpy library. In contrast to previous approaches, our methodology simultaneously models both afterglow and KN. We find that all GRBs in our sample have a KN, but we were unable to confirm or exclude its presence in GRB 150101B. A BNS progenitor is favored for GRB 160821B, GRB 170817A/AT2017gfo, GRB 211211A, and GRB 230307A. For GRB 150101B and GRB 191019A, we obtain a slight preference for NSBH scenario, while a BNS is also viable. For KN emission, we find that the median wind mass $\langle M_{\rm wind}\rangle=0.027^{+0.046}_{-0.019}$ $M_{\odot}$ is larger than the dynamical $\langle M_{\rm dyn}\rangle = 0.012^{+0.007}_{-0.006}$ $M_{\odot}$. We find that $M_{\rm wind}$ and the beaming corrected kinetic energy of the jet can be attributed as $log(M_{\rm wind})=-20.23+0.38\,log(E_{0,J})$. We confirm the results of numerical simulation that $\tildeΛ$ increases with decrease in $\mathcal{M}_{\rm \,Chirp}$. Our work shows that EM modeling can be effective for probing the progenitors, and for the first time presents the progenitor properties of a sizable sample of merger-driven GRBs.
Cepheids are fundamental distance indicators, playing a crucial role not only in the cosmic distance ladder but also in mapping the structure, kinematics, and extinction properties of the Milky Way. Using high-precision photometry and parallaxes from $Gaia$, we identify a significant anti-correlation between the $G$-band extinction coefficient and reddening for Galactic Cepheids, quantified as $R_G = 1.921 \pm 0.060 - (0.107 \pm 0.022)\,E(G_{\mathrm{BP}} - G_{\mathrm{RP}})$. We propose that this anti-correlation is partly driven by the pronounced non-linear effects inherent to the broad $Gaia$ bands, while the remaining part arise from the $R_V$ variations caused by diverse interstellar medium. Adopting a fixed $R_G$ would not only lead to an overestimation of the metallicity dependence of Cepheid luminosities, but also systematically underestimate the distances to highly reddened Cepheids. Moreover, the strong reddening dependence of $R_G$ makes Wesenheit function based on it unsuitable for highly reddened Cepheids, since the definition of Wesenheit magnitudes requires a fixed extinction coefficient. In contrast, infrared-based distances, being less affected by non-linear effects and insensitive to $R_V$, provide the most reliable Cepheid distances at present. This work emphasizes the importance of accurately determining $R_V$ for Galactic Cepheids and accounting for non-linear effects in distance measurements, particularly in the optical bands.
Traditional direct detection experiments lack the sensitivity to probe the sub-GeV dark matter (DM), primarily due to the low energy of the expected nuclear recoils. In this work, we investigate cosmic-ray (CR) upscattering as a mechanism to accelerate DM particles to detectable velocities in underground experiments. By analyzing four models of DM-nucleon interactions -- namely scalar, vector, pseudoscalar, and axial-vector mediators -- we derive constraints on the coupling parameters using data from the LZ, XENON, and Borexino experiments, covering mediator mass from $10^{-6}$ to $1$ GeV. As the mediator mass varies, the shift in dominance between momentum transfer and mediator mass leads to a turnover in the constraints around $10^{-2}$--$10^{-3}~\mathrm{GeV}$. Our results extend the reach of direct detection into the sub-GeV window and clarify the critical role of momentum dependence in light-mediator scenarios.
We present optical-ultraviolet photometry and optical spectra for the type II supernova (SN) 2022acko. The spectroscopic observations span phases from $\sim$ 1.5 to $\sim$ 60 days after the explosion, while the light curve was observed up to $\sim$ 300 days. The V-band peak is $-15.5 \pm 0.3$ mag, suggesting that SN 2022acko is a low-luminosity SN II (LLSN). The overall observed properties of SN 2022acko are consistent with those produced by a lower mass progenitor ($\rm M_{ZAMS} \sim $9-10M$_{\odot}$). The spectra at $t=1.5$d and $t=2.5$d exhibit a broad emission feature peaking near 4600 Å(the ``ledge'' feature), which we interpret as blueshifted He II 4686 Ålines arising from the ionized ejecta. Moreover, a possible flash-ionized (FI) emission line of H$α$ (FWHM $\sim 1100\ \rm km \ s^{-1}$) was superposed on the broad emission component of H$α$ P-Cgyni profile in the $t=1.5$d spectrum. Assuming an ejecta velocity of $\rm 12000\ km\ s^{-1}$, the rapid disappearance of this narrow H$α$ emission line within two days suggests highly confined CSM within $\sim \rm 2\times10^{14}\, cm$. Assuming a spherically symmetric CSM, the mass loss rate within this radius is estimated to be $\rm \sim 5 \times 10^{-4} M_{\odot} \ year^{-1}$ based on our hybrid light curve model. The early ``ledge'' feature observed in SN 2022acko have also been observed in other SNe II, suggesting that early-phase circumstellar interaction (CSI) is more common than previously thought.
Tidal disruption events (TDEs) occur when a star is gravitationally disrupted by the tidal field of a supermassive black hole during a close encounter. Radio emission has recently been detected in TDEs and is commonly attributed to synchrotron radiation from both wind and jetted outflows. However, several TDEs exhibit bright radio flares at late times, which cannot be easily explained if the wind is launched promptly after the stellar disruption. In this study, we model the radio light curves of TDEs with delayed radio flares using three scenarios: an instantaneous wind, a delayed wind, and a delayed relativistic jet. We show that the instantaneous wind model struggles to reproduce delayed radio flare events, indicating the necessity of an additional delayed outflow component. In contrast, the delayed wind model provides a consistent explanation for the observed radio phenomenology, successfully reproducing events both with and without delayed radio flares. For some delayed radio flare events (e.g., ASASSN-15oi and AT 2019dsg), both the delayed wind and delayed jet models can reproduce the observed radio light curves. The delayed jet model produces x-ray and optical emission that is detectable at typical TDE distances, in contrast to wind-driven scenarios. This highlights how multiwavelength observations offer an effective means of distinguishing among possible outflow mechanisms.
Microquasars are radio-emitting X-ray binaries accompanied by relativistic jets. They are established sources of 100~TeV gamma rays and are considered promising candidates for cosmic-ray acceleration. Motivated by recent detections of $\sim 100~$TeV photons from Cygnus~X-1 and $\sim~$PeV photons from Cygnus~X-3 by the Large High Altitude Air Shower Observatory (LHAASO), we employ the Astrophysical Multimessenger Emission Simulator (AMES) to model their multimessenger emission considering compact outflow regions as cosmic-ray accelerators, spanning from radio to ultra-high-energy gamma rays. Our results show that the observed $>$TeV gamma rays can originate from either $pγ$ or $pp$ interactions, depending on the location and physical conditions of the emission region, while also reproducing the lower-energy spectra. The different configurations yield unique, observationally testable predictions. In the $0.1-10$~TeV energy range, where current observations provide only upper limits, they predict either a deep dip, a mild suppression, or a power-law spectrum. Additionally, models involving AU-scale blob regions predict strong variability, while those invoking more extended and static external zones show more stable behavior. We also provide a possible qualitative explanation for the distinct modulation patterns across different energy bands, which relies primarily on changes in the Doppler factor and external $γγ$ absorption. Finally, our neutrino predictions, which properly account for muon and pion cooling effects, reveal a significantly suppressed flux, indicating that detecting these sources may be more challenging than previously anticipated.
The MAMMOTH-LyC survey is a cycle 30 Hubble Space Telescope (HST) medium program obtaining 18-orbit-deep WFC3/UVIS F225W imaging in two massive galaxy protocluster fields at $z\sim2.2$. We introduce this survey by reporting the discovery of J1244-LyC1, a strong Lyman continuum (LyC) leaker at $z = 2.39$, exhibiting clear merger signatures. J1244-LyC1 has a highly significant ($10σ$) LyC detection, corresponding to an absolute escape fraction of $f_{\mathrm{esc}} \! =\!36\%\pm4\%$ ($1σ$). The LyC emission is spatially resolved into multiple peaks that coincide with the system's disturbed morphology, confirming genuine multi-site LyC leakage. With a stellar mass of $10^{10.2}{M_\odot}$, J1244-LyC1 is both the first confirmed high-redshift LyC-leaking merger and the most massive LyC emitter known to date. We interpret J1244-LyC1 as a merger-driven starburst system in which tidal interactions have disrupted the interstellar medium, creating multiple low-column-density pathways that facilitate LyC escape. This discovery provides the first direct evidence of spatially resolved LyC escape in a merging system, offering new insight into the potential role of major mergers in driving the cosmic reionization.
We present H$α$ luminosity function (LF) measurements at redshifts $z\sim1.3$ and $z\sim2.0$ using JWST NIRISS grism data from the GLASS-JWST survey. Based on emission lines spectroscopically identified in the F115W, F150W and F200W filters, we select 99 H$α$ emitters. Through detailed effective volume and completeness analysis for each source, we construct the H$α$ LF in two redshift bins. Thanks to the sensitivity of NIRISS WFSS and gravitational lensing magnification, our sample reaches intrinsic H$α$ luminosities $\sim$10 times deeper than previous grism surveys, down to $L_{\rm Hα}\sim10^{40.5}~\rm erg~s^{-1}$ at $z\sim1.3$ and $L_{\rm Hα}\sim10^{40.9}~\rm erg~s^{-1}$ at $z\sim2.0$ with completeness larger than 0.8, corresponding to star formation rates of 0.4 and 1.0 $M_{\odot}~\rm yr^{-1}$, respectively. We robustly constrain the faint-end slope of the H$α$ luminosity function to be $-1.50^{+0.14}_{-0.08}$ at $z\sim1.3$ and $-1.60^{+0.17}_{-0.09}$ at $z\sim2.0$ after considering the cosmic variance of $\sim 20\%$, consistent with previous estimations. The emission-line samples presented here will enable further detailed studies of galaxy properties including metallicities. We find a negligible contribution from bright active galactic nuclei in our sample. We estimate integrated cosmic star formation rate densities of $0.097^{+0.015}_{-0.016}~M_{\odot}~\rm yr^{-1}~Mpc^{-3}$ at $z\sim1.3$ and $0.129^{+0.025}_{-0.030}~M_{\odot}~\rm yr^{-1}~Mpc^{-3}$ at $z\sim2.0$. The methodology presented here can be readily applicable to other JWST slitless spectroscopic datasets and future wide-field slitless surveys, including those from Euclid, Roman, and the Chinese Space Station Telescope.