Most interesting arXiv papers
Wednesday 11-03-2026
http://arxiv.org/abs/2603.09905 (pdf) (57)
The One-hundred-deg^2 DECam Imaging in Narrowbands (ODIN) survey is conducting the widest-field deep narrow-band imaging of the equatorial and southern skies. ODIN uses three custom-built narrow-band (NB) filters that sample Lya-emitting galaxies (LAEs) within thin cosmic slices centered at z=2.4, 3.1, and 4.5. In this work, we utilize extensive DESI spectroscopy of ODIN-selected galaxies in the COSMOS and XMM-LSS fields to validate our LAE selection. 2-4 hr exposures with DESI yielded redshift confirmation of 3,075 ODIN LAE candidates with NB magnitudes brighter than 26 mag. Restricting to objects that yield high-confidence redshifts, the confirmation rates are (93, 96, 92)% at z=(2.4, 3.1, 4.5). The primary contaminants consist of active galactic nuclei at the expected Lya redshift range and lower redshifts (C IV, C III]), with the remainder being star-forming galaxies ([O II] and [O III]). We find minimal contamination from [O II] emitters in our sample (< 1%), implying that our REW>20 A narrow-band excess photometry requirement is sufficient to remove them.
http://arxiv.org/abs/2603.09828 (pdf) (51)
We present a trade study of possible ultra-deep surveys with the Nancy Grace Roman Space Telescope, optimizing the depth-area-filter parameter space for high-redshift galaxy science. Using a mock galaxy catalog derived from a 2 sq. degree lightcone created using the Santa Cruz semi-analytic model and populated with over 7.6 million galaxies at 0<z<10 with M_UV < -15, with realistic clustering and synthetic photometry, we evaluate sixteen 500-hour survey configurations spanning 0.28-2 sq. degrees and four filter combinations. We demonstrate that even a single Roman pointing dramatically reduces cosmic variance compared to HST-like observations, more faithfully recovering the true UV luminosity function. For each survey configuration, we explore photometric redshift recovery, sample contamination, and measurements of the rest-UV luminosity function and non-ionizing UV luminosity density across four redshift bins at z 6-9. We find that inclusion of the R062 filter is essential for studies at z 5-6, reducing sample contamination from nearly 100% to negligible levels and recovering the bright end of the luminosity function. The F184 filter improves galaxy recovery at z>9 and is critical for stellar contamination removal at all redshifts. Based on these results, we recommend that a Roman ultra-deep survey cover at least two Roman pointings (0.56 sq. degrees) with all six filters (R062, Z087, Y106, J129, H158, F184), reducing uncertainties on the rest-UV luminosity density by factors of 2-4 relative to the deepest existing JWST programs. Building off of the Deep Tier of the High Latitude Time Domain Survey to add depth and filter coverage to existing (or planned) observations is an excellent option.
http://arxiv.org/abs/2603.08830 (pdf) (46)
Luminous quasars at z > 6 are key probes of early supermassive black hole (SMBH) growth, massive galaxy evolution, and intergalactic medium properties during cosmic reionization. However, their discovery is very challenging due to their scarcity and overwhelming contamination, as foreground ultracool dwarfs (UCDs) outnumber z>6 quasars by 2-4 orders of magnitude. In this work, we leverage the extensive coverage of DESI Legacy Survey DR10 to conduct a self-supervised search for quasars at z > 6, directly analyzing multiband optical images and minimizing the biases of traditional catalog-driven color-color selection criteria. By applying a contrastive learning (CL) method followed by spectral energy distribution (SED) fitting prioritization, we identified 1139 high-priority quasar candidates, for which we expect a competitive ∼1:1 quasar-to-UCD ratio based on literature samples. We spectroscopically confirm 16 new quasars at z = 5.94-6.45, achieving a 45
http://arxiv.org/abs/2603.09805 (pdf) (41)
Galaxies may suffer some starburst and quenched periods in their history due to e.g. galaxy mergers and feedback. However, semi-numerical simulations of the Epoch of Reionization (EoR) typically do not accurately model the effects of the star formation history (SFH) of galaxies. Keeping the same total ionizing photon budget from galaxies, we investigate how the ionization and heating of the Intergalactic Medium (IGM), as well as the associated 21-cm signal during the EoR, depends on the variations in the modeling of the SFH of galaxies. We adopt the Jiutian-300 N-body dark matter simulation and the semi-analytic model L-Galaxies 2020 to model galaxy formation. Using the galaxy catalog from L-Galaxies 2020 as input, we post-process the Jiutian-300 density field with the one-dimensional radiative transfer code Grizzly to model the reionization process and the 21-cm signal. We find that the ionized regions produced by galaxies with a SFH derived from L-Galaxies 2020 are slightly larger and warmer than the ones obtained with a constant SFR. For a fixed stellar mass, galaxies produce smaller ionized regions with increasing stellar mass weighted stellar age τ_age. This results in a different topology and timing of the IGM ionization and heating obtained from Grizzly. The SFH of galaxies is highly dependent on τ_age and redshift. Different models of the galactic SFH affect the gas heating and ionizing processes during the EoR, and as a consequence also the 21-cm global signal and power spectrum.
http://arxiv.org/abs/2603.09866 (pdf) (40)
JWST can pierce through dusty molecular clouds to study the early stages of star formation, where young star clusters are actively driving stellar feedback and still emerging from their natal cloud. We present a first look of the JWST/NIRSpec multiplex spectroscopy observations acquired by the Feedback in Emerging extrAgalactic Star clusTers (FEAST) program for the nearby spiral galaxy NGC628. We showcase JWST&aposs ability to resolve the spectral properties of emerging young star clusters (eYSCs) and their immediate interstellar medium (ISM) by focusing on a bright star-forming complex (0.5×0.5 kpc^2) in the northern spiral arm as a science proof-of-concept. The eYSC spectra are rich in ionized gas (from HII regions), as well as warm H_2 and polycyclic aromatic hydrocarbon (PAH) emission from photodissociation regions (PDRs), consistent with young star formation. Paα equivalent widths and H/He ionizing photon fluxes both indicate the presence of hot, young massive stars (O8.5V-O8V), consistent with photometry SED estimates. The ionized gas is highly correlated with H_2 and PAH emission, suggesting that the PDR morphology evolves as clusters emerge from their natal cloud. We find a photoionization-dominated regime from independent line diagnostics, with little contribution from Supernovae-driven shocks, highlighting the importance of pre-Supernovae feedback when massive stars are present. This pilot study showcases how JWST&aposs multiplex spectroscopy mode can disentangle the mechanisms present in the youngest stages of star formation for the first time outside the Local Group.
http://arxiv.org/abs/2603.09389 (pdf) (33)
Accurate and automated galaxy redshift determination is essential for maximizing the scientific return of spectroscopic surveys. In this paper, we propose a data-driven method to address this challenge. The method first learns a rest-frame representation of galaxy spectra using Non-negative Matrix Factorization (NMF). The method then reconstructs new spectra using this representation at different trial redshifts, and identifies the correct redshift by selecting the one that minimizes the reconstruction error. We apply our method to galaxy spectra from the Multi Unit Spectroscopic Explorer (MUSE), covering redshifts from 0 to 6.7. Our method achieves an overall success rate of 93.7%. We further demonstrate two applications: (i) the separation between true and false sources, and (ii) the detection of blended sources from one-dimensional spectra. Our results demonstrate that NMF-based representations provide a powerful and physically motivated framework for redshift estimation in current and future large spectroscopic surveys.
http://arxiv.org/abs/2603.09207 (pdf) (31)
Giant low surface brightness galaxies (gLSBs) are galaxies with extremely extended, faint, optical disks over 50 kpc in radius and have high total masses which can reach 10^12 solar masses. The existence of such galaxies is problematic for current models of galaxy formation, since the major mergers responsible for the large total mass would likely have destroyed the extended optical disk. Examining the gas content of these galaxies is an important step in determining their formation mechanism, whether it be through slow gas accretion or the large disk (re)forming after a major merger. We present neutral atomic hydrogen (HI) observations of 19 gLSBs identified with the Hyper Suprime-Cam Subaru Strategic Program survey. Although most have high HI masses, they are generally lower than expected based on their large optical sizes, and we do identify some gLSBs with unusually low gas content. The HI spectra of these galaxies show evidence for a rotational disk, though these disks are more asymmetric than other galaxies with comparable mass. Four galaxies with similar surface brightness profiles to the gLSBs have also been selected from the Numerical Investigation of a Hundred Astrophysical Objects (NIHAO) simulation for comparison. There is evidence for significant galaxy mergers in the past for three of these NIHAO galaxies and these three galaxies show similar asymmetry in their HI spectra. Together, these results could indicate the large optical disk of gLSBs are the result of a recent merger.
http://arxiv.org/abs/2603.09917 (pdf) (31)
Recent studies have uncovered a rare population of super-massive (M* > 1e11 Msun) star-forming spiral galaxies, Super Spiral Galaxies (SSGs), whose existence challenges classical quenching scenarios. We investigate the resolved star-forming and molecular-gas properties of UGC 8179 (z=0.052, log(M*/Msun)=11.62) and assess whether its local star formation (SF) follows the same physical processes as typical Star-Forming Main Sequence (SFMS) spirals. We combined the first NOEMA CO(1-0) interferometric observations of an SSG with pixel-by-pixel SED fitting, based on archival UV-to-mid-IR imaging. Our 3"x3" pixel maps provide resolved measurements of M*, SFR and molecular gas surface densities across its extended disc. UGC 8179 hosts a massive rotating molecular gas reservoir of M_H2 = 1.02 1e10 Msun, yielding a standard molecular gas fraction, with typical depletion time ∼1 Gyr in the observed region, despite its extreme mass. We derived lower limits of log(fmol) > -1.61 ±0.06 and log(tdep) > -8.82 ±0.13 at the scale of the galaxy. The large spatial extent of UGC 8179 enables us to probe low surface-density regimes hardly accessible in nearby disks (Σ* < 1e7 Msun / kpc2 ; Σ_SFR < 1e-3.5 Msun/yr/kpc2). All three resolved scaling relations (rSFMS, rKS and rMGMS) are well defined. The rKS slope (0.87 ±0.09) is broadly consistent with unity, indicating standard local SF processes. The rSFMS shows a shallower global slope (0.80 ±0.02) due to a central suppression in sSFR ( -0.5 dex). This break suggests the influence of a bulge, driving a transition to a more dynamically regulated SF regime in the inner disc. UGC 8179 provides evidence that SSGs can sustain standard local SF processes while exhibiting central dynamical regulation at high stellar surface densities.
http://arxiv.org/abs/2603.09739 (pdf) (30)
Protoclusters represent sites of accelerated galaxy formation and extreme astrophysical activity characteristic of dense environments. Identifying massive protoclusters and mapping their spatial structures are therefore crucial first steps in understanding how the large-scale environment influences galaxy evolution. We combine wide-field Lyα imaging from the ODIN survey with extensive DESI and ancillary spectroscopy across the extended COSMOS and XMM-LSS fields (≈14 deg^2) to search for massive protoclusters. We confirm six systems at z≈ 2.4 and z≈ 3.1, reconstruct their three-dimensional structures, estimate descendant halo masses, and, for one structure at z≈ 3.12, demonstrate that overlapping narrowband filters (NB497 and N501) provide accurate redshift tomography for emission-line galaxies. One protocluster at z≈ 2.45 overlaps with one of the LATIS tomographic fields, enabling direct comparison between galaxy and H i overdensities traced by Lyα forest absorption. Another at z≈ 3.12 hosts a massive quiescent galaxy (M_∗≈ 1.2 × 10^11M_⊙), indicating early quenching in a dense environment. By comparing Lyα emission properties across environments, we find that protocluster galaxies exhibit higher median line fluxes and a deficit of faint emitters relative to the field. The effect is strongest when both 2D and 3D density information are combined, indicating that galaxies in the densest protocluster cores are most affected by environmental processes. This effect is stronger at z≈3.1 than at z≈2.4, suggesting possible redshift evolution.
http://arxiv.org/abs/2603.09308 (pdf) (30)
Dark matter comprises 85% of the entire mass of the Universe, but the fundamental nature of its constituent particles remains elusive. In this thesis, I test for two competitive dark matter models: the conventional heavy particle paradigm, and dark matter being ultralight bosons of mass ∼ 10^-22eV (ψDM). More specifically, I test for the faint-end turnover induced by ψDM models, exploiting the strong lensing power by massive galaxy clusters to probe intrinsically fainter magnitudes. A key challenge for such an analysis would be contamination by low-z galaxies sharing similar observed SEDs as high-z galaxies. As I will demonstrate, such a contamination issue is generally severe and may wash out the faint-end turnover signatures. I also show that ∼ 50
http://arxiv.org/abs/2603.08940 (pdf) (29)
Understanding the formation and evolution of the Milky Way (MW) requires detailed knowledge of its satellite galaxies. In this study, we focus on the Sextans dwarf spheroidal (dSph) galaxy, a faint, dark matter (DM)-dominated satellite, to investigate the role of tidal and baryonic effects in shaping its observed properties. Using tailored N-body simulations, we explore possible orbits of Sextans in different MW models to reconstruct its progenitor&aposs properties. Our simulations demonstrate the stars in Sextans are only mildly affected by galactic tides and the stellar kinematics provide robust constraints on its dynamical mass within the half-light radius, while the tidal mass loss of its DM component depends primarily on MW mass. The recovered infall mass of Sextans ranges from 1.22 to 3.14×10^9 M_⊙ for MW masses from 0.8 to 2×10^12 M_⊙. If the DM density remained as cuspy as NFW profile, the infall mass would be smaller by a factor of 2. Although with large ranges, the possible infall masses of Sextans recovered by our simulations are consistent with the stellar mass-halo mass relation in TNG50 and abundance matching results. We find some cases for the cuspy DM density profile where the infall mass is smaller than 10^9 M_⊙, possibly indicating that star formation in Sextans is more efficient than in other satellites. The recovered DM halo structural parameters from our simulations provide valuable constraints for future studies on the DM content and formation history of Sextans.
http://arxiv.org/abs/2603.08996 (pdf) (28)
Measurements of the electron density of populations of extragalactic HII regions in nearby galaxies remain limited, despite the relevance of this quantity for characterizing the porosity of the interstellar medium and the escape of the ionizing radiation. We initiated a project aimed at analyzing the root-mean-square electron density ne_rms, the in-situ density ne and the volume filling factor (phi) of extragalactic HII regions, investigating the dependence of these attributes on nebular and host galaxy properties. We present an image-segmentation methodology for constructing homogeneous HII region catalogues, and apply it to two pilot galaxies: NGC 2403 and NGC 628. We derive ne_rms from their Halpha luminosities and equivalent radii (R_eq), and obtain ne and phi for spectroscopic subsamples. While ne is below 300 cm^-3, ne_rms is typically one to two orders of magnitude lower, implying that phi is in the range 10^-4 to 10^-1. The two galaxies exhibit a similar size-density relation, which breaks for R_eq > 50 pc, show at most a weak dependence of ne_rms on galactocentric radius for NGC 2403, and no clear dependence of ne or phi on these parameters. Combining these results with published data, ne_rms presents tentative scaling relations with the median HII region size, the fraction of large regions in the parent galaxy, and the star formation rate surface density. These trends, if confirmed, would provide new constraints for massive cluster formation models and important clues for interpreting dependencies observed at high redshift, underscoring the necessity of consistently extending this analysis to larger samples.
http://arxiv.org/abs/2603.09976 (pdf) (28)
We present measurements of black hole (BH) properties of 861 Type-1 active galactic nuclei (AGNs) in the North Ecliptic Pole (NEP)-Wide field. These AGNs are detected in both optical and mid-infrared (MIR) surveys and are identified as Type-1 AGNs in optical spectroscopic surveys. By performing spectral energy distribution (SED) and line fitting, we obtained their MIR continuum luminosities (L_MIR) as well as full width at half maximum (FWHM) values for the C4, Mg2, Hβ, and Hα lines. Using these measurements, we derived bolometric luminosities (10^43.20–10^47.27 erg s^-1), BH masses (10^7.29–10^9.67 M_⊙), and Eddington ratios (10^-2.74–10^-0.08) for ∼450 objects over a wide redshift range (z=0.09–4.71). The use of L_MIR and FWHM values effectively alleviates the effects of dust extinction, enabling reliable estimates of BH properties even for dust-obscured AGNs. Moreover, we find that 34
http://arxiv.org/abs/2603.09879 (pdf) (28)
The capability of reconstructing the past star formation history of dwarf elliptical galaxies out of the Local Volume relies on modelling bright stellar populations currently evolving through the red giant branch (RGB) and the asymptotic giant branch (AGB) phases. Recent studies proposed the use of the relative fractions of RGB and AGB stars populating specific boxes in the observational colour-magnitude plane to infer the epoch within which 90 We aim at understanding the physical process of stellar evolution that are constrained by the relationship between the relative fraction of AGB and RGB stars of dwarf galaxies and the T90 epoch. We use updated stellar models that include the description of dust formation in the wind, to undertake a population synthesis approach, to allow monitoring the variation of the NAGB/NRGB ratio with time. The effects of some specific ingredients, such as the mass loss experienced by low-mass stars during the RGB phase, and the details of the time variation of the star formation rate, are extensively explored and tested against data. The mass lost by low-mass stars during the RGB evolution proves the most relevant ingredients affecting the time variation of NAGB/NRGB: at metallicities 1/10 solar, a mass loss 0.25Msun is required to reproduce the observations. This analysis allows to derive a relationship between NAGB/NRGB and T90, with a 1 Gyr uncertainty on T90.
http://arxiv.org/abs/2603.08804 (pdf) (27)
Low-mass stars from the first epoch of star formation may still persist in the Milky Way and its satellite dwarf galaxies today; however, their detection is confounded by surface pollution from interstellar accretion and internal mixing, which obscure their primordial composition and blur their distinction from second-generation stars. Asteroseismology offers a probe of the internal structure and evolutionary state of stars, and hence may aid in the search for primordial stars. In this second paper of the series, we present the first non-radial adiabatic pulsation analysis of low-mass, metal-free stellar models. We use a 0.85 M_⊙ red giant as a case study and compare its seismic signatures with those of higher-metallicity models. At the same central hydrogen fractions, Pop III main-sequence models display systematically higher r_02_02/Δν ratio and lower Δν than metal-enriched analogues, a direct consequence of their larger internal sound speeds and steeper core-envelope stratification. To interpret the structural dependence during giant evolution, we introduce a composite asteroseismic diagnostic, ψ/ΔΠ_1, which traces how metallicity influences the balance between acoustic and buoyancy cavities through its imprint on opacity, core contraction, and mean molecular-weight gradients. Pop III models occupy a distinct locus in the ψ-ΔΠ_1 plane due to their radiative interiors with lower mean densities and delayed development of core mean molecular weight gradients. We find that asteroseismology is a powerful diagnostic for identifying relic Pop III stars despite potentially polluted surfaces, providing a clear pathway for future searches of the Galaxy&aposs oldest surviving stars with upcoming surveys.
http://arxiv.org/abs/2603.09339 (pdf) (26)
A large fraction of the extragalactic sky is obscured by foreground dust and stars along the plane of the Milky Way, leaving a major gap ( 20%) in whole-sky maps of large-scale structures – an incompleteness that is even more severe for peculiar velocity samples. This has long limited an unambiguous interpretation of observed cosmic flows and their connection to the underlying mass-density field. We present a new hybrid reconstruction methodology which combines 65,518 galaxy peculiar velocity distances from the CF4++ catalogue (Courtois2025) with 8283 new galaxy redshifts observed near the southern Galactic plane (|b| <= 10 degrees) Zone of Avoidance. A major advance is the inclusion of 2176 high-sensitivity, interferometric HI redshifts obtained with the SARAO MeerKAT telescope which for the first time provide coverage of the innermost 3degrees-wide strip of the southern ZOA and to unprecedented depth. This hybrid redshifts peculiar velocities approach yields a substantially revised view of the inferred overdensities in and around the ZOA. In particular, the Vela supercluster emerges as a dominant mass concentration, rivaling the Shapley concentration and exceeding the mass associated with Laniakea and the Great Attractor region. With a total mass of 33.8 10^16 Msol, a characteristic radius of 70 hmpc, and a double core morphology at a distance of 189 hmpc, Vela dominates the mass budget and gravitational influence of the southern Zone of Avoidance. These results provide the most complete and dynamically consistent picture to date of the southern Zone of Avoidance and demonstrate the transformative potential of hybrid reconstruction techniques tailored for the next generation of large-scale surveys.
http://arxiv.org/abs/2603.08787 (pdf) (21)
We present new JWST/NIRCam 4-5 μm (F410M, F430M) and JWST/MIRI 18-25 μm (F1800W, F2100W, F2550W) imaging detections of the nearby (3.6 pc) cold (275 K) gas giant exoplanet ε Ind Ab. The F2550W detection of ε Ind Ab constitutes the longest wavelength image of an exoplanet acquired to date. Combining three decades of radial velocity monitoring, Gaia-Hipparcos absolute astrometry, and relative astrometry from direct imaging (including the new NIRCam astrometry), we conduct a comprehensive re-analysis of ε Ind Ab&aposs orbit and obtain a dynamical mass M_Ab = 6.5^+0.7_-0.6 M_Jup. Using ε Ind Ab&aposs NIRCam and MIRI photometry, we assemble the first 4-25 μm spectral energy distribution (SED) of a cold gas giant outside the Solar System. The NIRCam photometry supports a metal-enriched atmosphere for ε Ind Ab based on analysis with atmospheric model grids, consistent with predictions from the giant planet mass-metallicity relation. While the current data do not provide definitive evidence for or against the presence of water ice clouds, we tentatively find that the H_2O vapor absorption-dominated F2550W photometry is systematically brighter (>1σ, but <2σ) than predictions from cloud-free/rainout chemistry models and better explained by a cloudy model. We calculate a bolometric luminosity of log L_bol/L_⊙ = -7.23 ± 0.03 dex by directly integrating ε Ind Ab&aposs SED. Combining this with the planet&aposs dynamical mass and age (3.5 ± 1.0 Gyr), we demonstrate excellent agreement with evolutionary model predictions in a new regime of low luminosities, low masses, and old ages. Our results establish ε Ind Ab as a benchmark system for planetary evolution studies and set the stage for the detailed atmospheric characterization of this frigid extrasolar world.
http://arxiv.org/abs/2603.08976 (pdf) (20)
Low-metallicity massive stars are assumed to be progenitors of certain supernovae, gamma-ray bursts, and gravitational wave emitting mergers. These exotic phenomena contribute to their host galaxies through strong ionizing radiation and mechanical feedback. Here we investigate a certain type of very metal-poor (0.02 Zsun) hot massive single stars that rotate fast and evolve chemically homogeneously. Combining state-of-the-art theories of stellar evolution and stellar atmospheres modelling we predict synthetic spectra of these stars corresponding to different masses and evolutionary phases. The predicted spectra in early evolutionary phases is classified mainly as an early-O type giant or supergiant while in later evolutionary phases most of our model spectra are assigned to the WO-type spectral class. The Hubble Space Telescope&aposs (HST) Ultraviolet Legacy Library of Young Stars as Essential Standards (ULLYSES) program will enable us to compare our predicted spectra with observations of stars of similar nature (e.g., metal-poor stars in Sextant A).
http://arxiv.org/abs/2603.09969 (pdf) (20)
In the late Universe, and on cosmological scales, dark matter is conventionally assumed to be collisionless, as a consequence of the strong existing bounds on dark matter interactions at the Cosmic Microwave Background last-scattering surface. Challenging this lore, here we show that dark matter interactions can be naturally weak at early times, but then grow to observationally relevant strengths at very late times, even significantly after reionization. This is realized if dark matter recouples to a dark radiation species in the range of redshifts probed by the current generation of galaxy surveys. We systematically study, for the first time, the phenomenology of this dark matter recoupling scenario. A combination of Cosmic Microwave Background and Baryon Acoustic Oscillation data show that the interaction needs to be weak at present, if the entirety of dark matter couples to dark radiation. From a complementary perspective, a 4% fraction of dark matter could still be strongly interacting today. Implications for a microscopic model realizing the recoupling dynamics are discussed.
http://arxiv.org/abs/2603.09762 (pdf) (20)
Large time-domain sky surveys generate extensive multi-year catalogs of light curves in which scientifically valuable transients, such as supernovae (SNe), are vastly outnumbered by artifacts and routine star variability. While supervised machine learning models can efficiently filter known classes, they struggle with extreme class imbalance and may overlook rare or novel events. Conversely, unsupervised anomaly detection provides broad discovery potential but lacks targeted sensitivity. We present a hybrid strategy that integrates a supervised SN probability score (SN-score) into the PineForest active anomaly detection framework to enhance SN discovery rate in the 23rd data release of the Zwicky Transient Facility. We train a binary classifier using light-curve features of spectroscopically confirmed SNe from the ZTF Bright Transient Survey, achieving a ROC-AUC approximately 0.98. Incorporating the SN-score as an additional feature, together with a small set of labeled priors, significantly accelerates the discovery of SN-like transients across ten extragalactic ZTF fields. This method increases discovery efficiency without compromising the ability to identify diverse astrophysical anomalies. Application of the combined methodology resulted in the discovery of seven previously unreported SN candidates, one AGN candidate, one unusual Galactic variable star SNAD283, as well as two host galaxies exhibiting multiple supernova events. These results demonstrate its value for scalable and expert-guided transient search in current and future surveys, including the Vera C. Rubin Observatory Legacy Survey of Space and Time.
http://arxiv.org/abs/2603.09670 (pdf) (20)
We present a study of recent star formation in the dwarf irregular galaxy DDO 43 using GALEX FUV and WISE NIR imaging. We identify regions of elevated FUV flux, indicating unobscured star-forming activity across much of the galaxy. To further characterize the stellar content, we compare the FUV fluxes to archival WISE W1 and W2 infrared data across 56 regions of interest. A general correlation is found between the FUV and infrared fluxes, suggesting spatially coherent star formation throughout the galaxy. A few regions, however, show elevated infrared fluxes but little or no UV emission, potentially indicating localized, dust-obscured star formation.
http://arxiv.org/abs/2309.02298 (pdf) (19)
We use kinematic data of proper motions from Gaia of forty-two globular and open clusters from Large Magellanic Cloud (LMC) to explore the possibility of them having extragalactic origins. We find the difference between the proper motions of cluster stars and a surrounding patch of young LMC stars in each case. We find five globular clusters towards the north-east showing a high difference (> 0.11 mas/yr, or > 25 km/s). We also examine the statistical significance of this difference taking into account both measurement errors of cluster and surrounding stars as well as inherent dispersion of stellar motions in the local galactic environment. The five globular clusters (NGC 2005, NGC 2210, NGC 1978, Hodge 3 and Hodge 11) have mean proper motions that lie outside the 85% confidence interval of the mean of surrounding young stars, with a clear outlier (NGC 1978 outside 99.96% confidence) whose difference cannot be accounted for by statistical noise. A young cluster (NGC 2100) also fitting the criteria is ruled out owing to contrary evidence from literature. This indicates a possible interaction with a dwarf galaxy resulting in the accretion/disruption in path of the five globular clusters, or possibly one or more past merger(s) of smaller galaxy/galaxies with LMC from its north-eastern region. This direction also coincides with the location of Tarantula Nebula, suggesting the possibility of the interaction event or merger having triggered its star formation activity.
http://arxiv.org/abs/2603.08917 (pdf) (18)
We present a broadband spectral analysis of the black hole X-ray binary GX 339-4 during its 2021 outburst, covering both hard and soft spectral states. Using simultaneous observations from NuSTAR, NICER, and Insight-HXMT, we investigate the evolution of the accretion disk with a focus on the disk normalization derived from the diskbb component, which serves as a proxy for the apparent inner disk radius. In the standard single Comptonization model, the disk normalization in the hard state is more than an order of magnitude lower than in the soft state (∼0.3×10^3 vs. ∼3.0×10^3). This result contradicts the widely accepted view that the disk radius is smaller in the soft state than in the hard state. By incorporating an additional warm Comptonization component, the disk normalization in the hard state increases to values (≳ 10^4) exceeding those in the soft state (∼ 10^3), yielding results consistent with a physically truncated, cooler accretion disk. The results of this work support the presence of a dual-corona geometry in the hard state, comprising both a hot, optically thin corona and a warm, optically thick corona, while the soft state spectrum is well described by a single hot Comptonization component alone. Our findings emphasize the importance of including a warm corona in hard-state spectra, as it leads to a more physically consistent picture of the accretion geometry across spectral states.
http://arxiv.org/abs/2603.09489 (pdf) (17)
We present the first on-sky demonstration of dual-field interferometry at the CHARA Array and the first direct resolution of the inner Ba–Bb subsystem in the bright hierarchical triple α Piscium. Using H-band fringe tracking on component A with MIRC-X to stabilize K-band science fringes on component B with MYSTIC, we detected a companion at a projected separation of 7 mas, confirming a long-suspected but previously unresolved short-period subsystem within the B component. The nearly equal H/K-band flux ratio indicates that Ba and Bb are near-twin F-type stars, consistent with the two narrow-lined components seen in optical spectra of B. By combining CHARA interferometry with archival VLTI/GRAVITY astrometry and radial velocities from archival and new spectroscopy (NARVAL and ARCES), we derive a well-constrained orbit with a period of P = 25 d, eccentricity e ≃ 0.6, and inclination i ≃ 65^∘, yielding precise dynamical masses of 1.668±0.033 M_⊙ and 1.646±0.029 M_⊙. No additional companion is detected down to ΔH ≈ 5 at separations of 0.2–2 AU. We also obtained dual-field differential astrometry of the wide A–B pair with a precision of 0.234 mas at a separation of 1.85”, with an error budget dominated by internal delay-line actuators, fringe-tracking performance and chromatic dispersion. While the long-period outer orbit is not refined by these measurements, their agreement with the published astrometric orbit provides an on-sky validation of the CHARA dual-field mode. These results establish α Psc as a well-characterized hierarchical system suitable for future benchmark studies and demonstrate CHARA&aposs new capability for off-axis interferometry and sub-mas astrometry on arcsecond-scale binaries.
http://arxiv.org/abs/2603.09457 (pdf) (16)
The parameter S_8, a key probe of cosmic structure growth, exhibits a persistent ∼3σ tension between high-redshift measurements from cosmic microwave background (CMB) anisotropies and low-redshift weak gravitational lensing observations. This discrepancy may indicate either unaccounted systematic effects or new physics beyond the standard ΛCDM cosmology. In this work, we constrain S_8 using the high purity CatNorth 1.5 million quasar candidate catalog and the Planck DR4 CMB lensing data across the broad redshift ranges through auto-correlation and cross-correlation analyses. To address the spatial incompleteness, we develop a machine-learning-based selection function that effectively suppresses the systematics-induced power spectrum excess on large scales. Our robust low-redshift measurements at z<1.5 yield S_8 = 0.844^+0.058_-0.056, consistent with the Planck 2018 CMB anisotropies constraints of S_8=0.834±0.016 but lower than the 0.879^+0.055_-0.055 reported by a previous work using the Quaia quasar candidate catalog. However, for high-redshift faint quasars at z>1.5, we find a lower value of S_8=0.724^+0.058_-0.054, likely due to the sample incompleteness and/or the foreground contamination. Further tests on the volume-limited samples exhibit a consistent trend: S_8 = 0.835^+0.053_-0.049 for z < 2, 0.824^+0.061_-0.062 for 0.4 < z < 1.5, and a lower value of 0.789^+0.062_-0.062 for the higher redshift range of 1.5 < z < 2.5. While future data may refine these results, our current measurements based on a large sample of quasar candidates show less evidence of the S_8 tension.
http://arxiv.org/abs/2603.09540 (pdf) (16)
We aim to develop a model-driven deep learning approach to age determination, by training neural networks on stellar evolutionary grids. Contrary to the usual data-driven deep learning approach of using prior age estimates as training data, our method has the potential for a wider and less biased range of application. The low computational cost of deep learning methods compared to bayesian isochrone-fitting allows for a broad analysis of large spectroscopic catalogues. We train multilayer perceptrons on different stellar evolutionary grids to map [M/H], MG, (GBP - GRP) to stellar age τ. We combine Gaia photometry and parallaxes, metallicities and α elements from spectroscopic surveys and extinction maps, which are passed through the neural networks to estimate stellar ages. We apply our method to the LAMOST DR10, GALAH DR3 DR4 and APOGEE DR17 spectroscopic surveys, for which we estimate the ages using the BaSTI tracks, along with other stellar evolutionary models. We leverage this novel technique to study, for the first time, differences in age estimates from several evolutionary grids applied on very large datasets. In addition, we date 13 open clusters and one globular cluster and find a median absolute deviation with literature ages of 0.20 Gyr. Along with the stellar ages catalogues from our estimates, we release NEST (Neural Estimator of Stellar Times), a python package to estimate stellar age based on this work, as well as a web interface. We show that, when using the same evolutionary grid, our method retrieves the same ages as a bayesian approach like SPInS, for only a fraction of the computational cost, with a 60,000 speedup factor for a typical star. This model-driven deep learning technique thus opens up the way for broad galactic archeology studies on the largest datasets available today and in the near future with upcoming surveys such as 4MOST.
http://arxiv.org/abs/2603.08808 (pdf) (16)
Axions and axion-like particles (ALPs) arise naturally in many extensions of the Standard Model and are among the well-motivated candidates for dark matter. In the presence of magnetic fields of galaxy clusters, Cosmic Microwave Background (CMB) photons can convert to ALPs, with the efficiency of the process governed by the cluster electron density and magnetic field profiles, the photon–ALP coupling strength (g_aγ), as well as the frequency (ν) of the photon at the redshift of the cluster. The CMB blackbody spectrum suggests that this resonant conversion also takes place at radio wavelengths, following the spectral behaviour of the ALP distortion signal. This opens a new window to search for ALPs using cosmic variance cancellation (CVC), with multi-frequency tracers of the same phenomenon in CMB photon–ALP resonant conversion. The constraints on the ALP signal ratios from different combinations of microwave and radio bands of the Simons Observatory (SO) and the Square Kilometre Array (SKA) can be significantly improved using CVC compared to the case of using auto-only spectra from the two experiments. With the large number of galaxy clusters that will be observed by SO and SKA, we will be able to obtain much more information using CVC, especially for low-mass ALPs with stronger signals. Using the auto-only spectra from galaxy clusters up to redshift z = 1 for inference of the normalized ratio parameter, we obtain a standard deviation of 5.9 × 10^-2 for an ALP mass m_a = 10^-14 eV, which improves to 1 × 10^-2 using CVC. This method provides a universal probe of the ALP distortion signal using its spectral dependence and can also invalidate false detections of the ALP signal based on its frequency behaviour in different bands.
http://arxiv.org/abs/2603.08984 (pdf) (15)
We present astrometric measurements for 13 cold brown dwarfs in the solar neighborhood (d < 20pc). By combining archival Spitzer data with our own Hubble Space Telescope (HST) observations, we achieve parallax uncertainties typically around 10%. Using Spitzer and HST photometry we compare our targets with other known late T and Y dwarfs in the Solar neighborhood, confirming that there is large intrinsic scatter in the near- and mid-infrared absolute magnitudes and colors of this population, further highlighting the diversity observed spectroscopically by several James Webb Space Telescope (JWST) programs. This scatter makes photometric distance estimates highly unreliable and, therefore, makes astrometric parallax measurements fundamental for a meaningful characterization of even the nearest cold brown dwarfs.
http://arxiv.org/abs/2603.09562 (pdf) (15)
Tidal Disruption Events are exotic astrophysical phenomena where matter from a star or the interstellar medium is captured by a supermassive black hole. The process liberates enormous energy, within a few months to a year timescale, enough to detect dormant black holes in near as well as the farthest galaxies. We revisit the long-term spectral variabilities associated with the jetted Tidal Disruption Event by exploring the archival X-ray data obtained with MAXI/GSC, Swift/XRT, and XMM-Newton observatories. Our analysis reveals that the spectral indices decrease non-monotonically as evolves with time. We also find that the soft (0.3-1.5 keV) and hard (1.5-10 keV) X-ray photon counts are highly correlated with a maximum correlation coefficient of 0.95 and peak at zero lag. Moreover, the soft and hard band variabilities obtained from XMM-Newton observations are highly correlated with a Pearson cross-correlation coefficient of 0.96. This indicates that the soft and hard X-ray photons are emitted from the same site, which is most likely a Compton cloud, i.e., the corona. Assuming the hard X-ray photons originate from the corona, we find that the coronal parameter undergoes rapid expansion during the early phases when accompanied by a relativistic jet launching and subsequently evolves toward a state of saturation with minor fluctuations in the latter stages. The temporal variation of the coronal size is consistent with a simple theoretical conjecture. We also discuss the application of our analytical outcomes to other jetted and non-jetted tidal disruption events.
http://arxiv.org/abs/2603.09228 (pdf) (14)
We present a multi-wavelength investigation of radio sources in the globular cluster M22 (NGC6656) using VLA, HST, and Chandra observations. Among the eight identified counterparts, we highlight VLA22 as the most promising stellar-mass black hole (BH) candidate. Its radio and X-ray luminosities follow the established L_R-L_X correlation for quiescent black hole low-mass X-ray binaries (BH-LMXBs), while its moderately steep radio spectrum and X-ray spectral hardening further support this classification. Analysis of two potential optical counterparts-a bright main-sequence star and a faint subgiant/red giant-suggests a binary system with a relatively long orbital period. The discovery of VLA22 consistent with recent retention models that stellar-mass BH can be retained within globular clusters over Hubble timescales. Additionally, VLA19 exhibits a characteristically inverted radio spectrum (α= 0.79 ± 0.39, S_ν∝ ν^α) indicative of a compact jet, while VLA40 also aligns with the BH L_R-L_X track, though both require further observations to definitively confirm their nature.
http://arxiv.org/abs/2603.09066 (pdf) (14)
Photospheric emission from relativistic gamma-ray burst (GRB) jets is a promising mechanism for producing the Band-like spectra observed in the prompt phase, yet the connections between jet structure, dissipation location, and polarization signatures remain unclear. We investigate time-dependent photospheric radiation transfer in structured relativistic jets by coupling two-dimensional axisymmetric special relativistic hydrodynamic (SRHD) simulations with Monte Carlo photon propagation. Photon escape and subphotospheric dissipation are characterized using the residual line-of-sight optical depth tau_out evaluated along each photon trajectory, allowing a direction-dependent treatment of photon decoupling in structured jets. The radiative transfer includes Klein-Nishina Compton scattering and polarization evolution using the Mueller matrix formalism. We perform a systematic parameter study exploring the effects of viewing angle, electron-positron pair loading (Z_pm), and the optical-depth window of subphotospheric dissipation. The model produces time-resolved spectra, peak-energy evolution E_pk(t), Band parameters, polarization degree Pi(E,t), and last-scattering statistics. We find that jet angular structure and the geometry of the line-of-sight optical depth strongly regulate spectral evolution and polarization signatures. The dissipation depth and pair loading jointly control the stability of E_pk, the formation of high-energy spectral tails, and the energy dependence of polarization. These results provide quantitative predictions for GRB prompt-emission spectra and polarization that can be tested with current and upcoming high-energy polarimeters.
http://arxiv.org/abs/2603.09867 (pdf) (13)
Quantifying the timescales of star cluster emergence from their natal clouds remains one of the main challenges in our understanding of the star formation process. These timescales are fundamental measurements of the star formation cycle within galaxies, yet they are difficult to constrain due to the complex interplay between stellar feedback and star formation across a vast range of physical scales. Here we present HST and JWST observations of thousands of young star clusters in four nearby galaxies (M51, M83, NGC 628, and NGC 4449). A substantial fraction of these clusters are still embedded within their natal gas and remain invisible at optical wavelengths. We constrain their emergence process by measuring the timescales required to disperse the surrounding material. We find a strong correlation between dispersal timescale and cluster stellar mass, with massive clusters emerging more rapidly than their lower mass counterparts. This is a critical constraint on simulations of star formation and stellar feedback, which struggle to fully reproduce the formation and emergence of star clusters. Our results emphasize the central role of massive clusters in driving the escape of ionizing radiation into the galactic medium. Finally, they impose important limitations to the time available for planet formation in massive cluster environments where disks get exposed to UV irradiation and further gas infall is shut off.
http://arxiv.org/abs/2603.09162 (pdf) (13)
Radio interferometry enables high-resolution imaging of astronomical radio sources by synthesizing a large effective aperture from an array of antennas and solving a deconvolution problem to reconstruct the image. Deep learning has emerged as a promising solution to the imaging problem, reducing computational costs and enabling super-resolution. However, existing DL-based methods often fall short of the requirements for real-world deployment due to limitations in handling high dynamic range, large field of view, and mismatches between training and test conditions. In this work, we build upon and extend the POLISH framework, a recent DL model for radio interferometric imaging. We introduce key improvements to enable robust reconstruction and super-resolution under real-world conditions: (1) a patch-wise training and stitching strategy for scaling to wide-field imaging and (2) a nonlinear arcsinh-based intensity transformation to manage high dynamic range. We conduct comprehensive evaluations using the T-RECS simulation suite with realistic sky models and point spead functions (PSF), and demonstrate that our approach significantly improves reconstruction quality and robustness. We test the model on realistic simulated strong gravitational lenses and show that lens systems with Einstein radii near the PSF scale can be recovered after deconvolution with our POLISH model, potentially yielding 10× more galaxy-galaxy lensing systems from the Deep Synoptic Array (DSA) survey than with image-plane CLEAN. Our results highlight the potential of DL models as practical, scalable tools for next-generation radio astronomy.
http://arxiv.org/abs/2603.09962 (pdf) (12)
We present a comprehensive polarimetric study of 26 nearby molecular clouds in four far-infrared bands (53 μm to 214 μm) using 52 archival SOFIA/HAWC+ datasets. Far-infrared dust polarization observations probe the plane-of-sky magnetic field. To investigate scale-dependent trends, we group the molecular clouds by distance and analyze the data at common angular (25”) and common physical (0.052 pc and 0.32 pc) resolutions. The two shorter wavelengths are more impacted by smoothing, exhibiting a larger decrease in percent polarization. We analyze the polarization spectrum – the polarization fraction as a function of wavelength – and find that it depends more strongly on column density than dust temperature. We find a "falling" spectrum at the 0.052 pc resolution, but find a "flat" spectrum at the 0.32 pc resolution, suggesting that resolution plays an important role in the observed polarization spectra. We propose that warm dust grain emission in small-scale structures (≲ 0.1 pc) traces different magnetic field geometries only resolved in our close regime data. There is no preferred magnetic field orientation across our data, which suggests that the magnetic field in our ∼ parsec scale regions is decoupled from the large-scale field that is primarily parallel to the Galactic plane. The relationship between percent polarization and column density varies between clouds, but the correlation between percent polarization and angular dispersion is consistent across regions. This compendium of dust polarization maps highlights the value of observing at multiple far-infrared wavelengths and will enable additional population-level studies of magnetic fields and dust across star-forming environments.
http://arxiv.org/abs/2603.09880 (pdf) (12)
The advent of high-precision Gaia parallaxes for Milky Way Cepheids enables per cent-level calibration of the local distance ladder and H_0. We revisit the Milky Way Cepheid calibration from Gaia EDR3 parallaxes using a fully forward-modelled Bayesian framework that simultaneously infers the period–luminosity relation, the Gaia parallax zero-point offset, and individual stellar distances while explicitly incorporating the disk geometry of the Galaxy through the distance prior and the selection functions specified in two distinct HST SH0ES campaigns. We derive an analytic treatment of the detection probability that accounts for magnitude, parallax, period, and extinction cuts and reduces the selection treatment to a tractable integral over distance and sky position. Posterior predictive checks show that this generative model matches well the observed distributions of parallaxes, magnitudes, and periods. Modelling Galactic structure and survey truncation self-consistently in a Bayesian framework yields period–luminosity parameters that agree with the SH0ES maximum-likelihood values at the <0.5 σ level, a consequence, we show, of the small intrinsic scatter of the Cepheid period–luminosity relation. Adopting, as recently advocated, a uniform-in-volume prior without simultaneously accounting for selection leads to a ∼ 0.05 mag bias in the period–luminosity zero-point and posterior predictive distributions incompatible with the observed data; this shift is mostly driven by the omission of the selection model. A consistent Bayesian treatment of Galactic structure and selection effects reinforces the local distance-ladder determination of H_0, and hence the Hubble tension with early-Universe inferences.
http://arxiv.org/abs/2603.09271 (pdf) (12)
Using multi-wavelength observations from the Solar Dynamics Observatory (SDO), we investigated the six-month decay process of the solar active region NOAA AR 12738 from April to October 2019. We systematically analyzed the region&aposs evolution by examining extreme ultraviolet (EUV) intensity variations, quantifying magnetic flux diffusion, and investigating thermodynamic changes via Differential Emission Measure (DEM) analysis. This study presents the first long-term tracking of a peripheral dimming region (dark moat), revealing its continuous areal decrease over time. DEM results reveal cooling plasma signatures and thermal restructuring, with the dimming region exhibiting a distinct temperature deficit in range 10^5.5 – 10^5.9 K. Potential field extrapolation identifies two dominant magnetic configurations: low-lying loops with cool plasma (<10^5.5 K), and high-arching structures connecting to the AR core, contributing to localized emission reduction. We found that the dimming is dominated by high-lying loops extending from the AR core, which are heated to temperatures above the main response of the 171 Å passband (> 10^5.8 K), consequently lacking plasma at the typical 10^5.8 K formation temperature. The thermal deficit, not just the absence of material, is the key driver of the reduced emission. Our results demonstrate that long-duration dimming provides a valuable diagnostic for understanding active region decay, thermal evolution, and coronal magnetic restructuring.
http://arxiv.org/abs/2603.08884 (pdf) (12)
Gamma-ray bursts (GRBs) are extremely powerful stellar explosions that have been observed to huge distances with redshifts exceeding 9. Although GRBs are not standard candles, one may standardize them by calibrating certain correlations that link an intrinsic parameter to an observed one. Two such correlations that have been discovered are the Amati relation and the Yonetoku relation. In this paper, we compile a large sample of 241 Swift long GRBs for the purpose of examining whether the Amati and Yonetoku relations are immune to redshift evolution. Our methodology encompasses two approaches: the first involves binning the data by redshift and fitting the two relations for each bin, then checking whether the fitting parameters evolve with redshift; the second approach involves using a redshift cutoff to divide the data into a low-redshift group and a high-redshift group, then checking whether the fitting parameters for the two relations are consistent with one another. Our results indicate that the Amati and Yonetoku relations are robust in the sense that they do not show any systematic or significant redshift evolution. Moreover, our results indicate that the high redshift bins show better fits compared to the low redshift bins, which indicates that the Amati and Yonetoku relations are more reliable for high redshift and hence are promising cosmological probes.
http://arxiv.org/abs/2603.09794 (pdf) (12)
If a gravitational wave event is lensed by a cluster or galaxy in our line-of-sight, it is expected that its host galaxy would also be lensed. Therefore, connecting lensed gravitational wave events even without direct optical counterpart could be feasible by identifying matching lenses in electromagnetic data and surveys. Seminal work has demonstrated the potential of this approach in LVK, Euclid, HST, JWST, and CSST mock data, motivating the need for a dedicated software package to perform such analyses in practice. Here, we present the alpha-version of silmarel, the first software package designed to bridge these cosmic signals and enable us analysis of real LVK gravitational-wave binaries together with telescope observations from instruments like Euclid or Hubble Space Telescope, and the future of multimessenger binary black hole lensing.
http://arxiv.org/abs/2603.09239 (pdf) (11)
We initiate the systematic search for planets in the 2023 data of the Korea Microlensing Telescope Network (KMTNet), focusing on those planets found by the KMTNet AnomalyFinder with low preliminary estimates of the mass-ratio, q<2× 10^-4. The 2023 season is the first for which the photometry of all events was re-reduced prior to the AnomalyFinder search, potentially increasing its sensitivity to planets. We find three strong low-q planet candidates, KMT-2023-BLG-0164 (q∼ 1.3× 10^-4), KMT-2023-BLG-1286 (q∼ 1.9× 10^-4), and KMT-2023-BLG-1746 (q∼ 8× 10^-5). KMT-2023-BLG-0164 is notable in that the source is projected on a very bright (I=16.0) foreground star, which is either the planet&aposs host or (more likely) a companion to the host. We obtain a spectrum, finding that its mass and distance are M∼ 1.0 M_⊙ and D∼ 1.5 kpc, the latter being the distance of the lens (D_L) regardless of whether the spectroscopic target is the host or its companion. We also analyze two other candidates, KMT-2023-BLG-0614 and KMT-2023-BLG-1593, which are unlikely to enter the statistical sample due to their ambiguous interpretations as possible non-planetary events.
http://arxiv.org/abs/2603.09810 (pdf) (11)
We present the first systematic study of neutrino emissions from massive stars, continuously tracking the late evolutionary stages through the early core-collapse supernova phase. Using progenitor and supernova models, we analyze the neutrino luminosities and spectra for progenitors with initial masses of 10–40 M_⊙. Our systematic analysis reveals that the compactness parameter (ξ_2.5) and carbon-oxygen core mass (M_CO) exhibit strong correlations with neutrino emission. In the pre-supernova phase, the time-integrated number of neutrinos correlates with ξ_2.5 when integrated over the final day and with M_CO for longer durations. For the early supernova phase (<200 ms post-bounce), the neutrino properties are relatively insensitive to the specific stellar evolution code used, allowing for a reliable extraction of physical correlations. We confirm that the neutrino emission features, including the electron neutrino burst properties and accretion-powered luminosity of other species, reflect the progenitor&aposs compactness. An evaluation of the observational feasibility for a nearby progenitor using a False Alarm Rate approach suggests that these correlations can persist even under practical detection conditions. Such a joint analysis of both phases provides complementary constraints on the internal structure. All calculated time-series data will be made publicly available.
http://arxiv.org/abs/2603.09305 (pdf) (11)
Longitudinal density fluctuations observed in stellar streams can result from gravitational interactions with massive perturbers in the Milky Way, such as dark matter subhalos. Analysing these density variations provides a powerful probe of properties (motion, mass, size, etc.) of the perturbing objects. However, caution is needed because density variations may arise naturally from internal dynamics of streams, namely epicycles. In this work, we focus on the Indus stellar stream, a remnant of an ancient dwarf satellite of the Galaxy. An Indus stream spanning ∼ 90^∘ is revealed in the southern Galactic sky using a comprehensive matched-filter analysis utilizing data from the Gaia mission. A spatial density model is fitted to the filtered map to quantitatively characterize the morphology, which demonstrates episodic density peaks and gaps in the stream. Through N-body simulations, we show that there are strong epicyclic motions of stars happening during tidal disruptions. The present-day longitudinal densities from simulations are comparable to the measurement from data, with similar numbers and locations of peaks and gaps, suggesting that the observed density should mainly be caused by epicycles. We also find that a cuspy dark matter halo for the Indus dwarf is likely to produce milder stellar epicyclic peaks compared to a cored halo which results in steeper peaks. This arises from different instantaneous mass loss due to distinct central mass distributions of halos, where a cored halo usually leads to severer tidal stripping. The observed density exhibits moderate peak sharpness, implying that Indus may have originally possessed a cuspy halo.
http://arxiv.org/abs/2603.08780 (pdf) (10)
With JWST, we are directly imaging cold ( 200-300K), solar-age giant exoplanets for the first time. At these temperatures many molecular features appear and water-ice clouds may condense and affect the emission spectrum; early photometric measurements of cold giant planets are already showing some tension with the predictions of cloud-free, solar-metallicity atmosphere models. Here we present new JWST/MIRI coronagraphic observations of the cold giant exoplanet Eps Ind Ab at 11.3um. Together with archival data, we use these new observations to study the atmosphere of this cold exoplanet, and we also re-fit its orbit, finding an updated mass of 7.6±0.7 Mj and an eccentricity of 0.24^+0.11_-0.08. The planet is significantly brighter (by 0.88±0.08 mag) at 11.3um than at 10.6um, indicating the presence of ammonia. However, this ammonia feature is shallower than expected. This could indicate a low-metallicity or nitrogen-depleted atmosphere, but our preferred explanation is the presence of thick water-ice clouds that suppress the ammonia feature and the near-IR emission of Eps Ind Ab. Photometry of the small but growing sample of cold, giant exoplanets demonstrates that they are consistently fainter than expected between 3 to 5um, consistent with the water-ice cloud hypothesis. 10.6um and 11.3um photometry of this cold exoplanet sample would be valuable to determine whether the suppressed ammonia feature is universal, and to frame a new open question about the underlying physical cause.
http://arxiv.org/abs/2603.08985 (pdf) (9)
Investigating the apparent anomalies in lithium (Li) surface abundance observed in the Sun and young stellar globular clusters holds significant promise for advancing our understanding of the mechanisms influencing Li depletion. This study delves into the intricate interplay between rotational mixing and rotational hydrostatic effects in pre-main-sequence (PMS) and main-sequence (MS) solar-type stars by employing grids of rotating models. We implement a novel approach in which both the magnetic field strength (B) and the mixing-length parameter (alpha-MLT) vary dynamically with stellar parameters. This avoids fixed values and aims to reduce free parameters while capturing key physical variability. Our models reproduce the observed Li abundance of Sun-like stars (A(Li) = 1.12 dex) consistent with the present-day solar value (1.1 +- 0.1 dex) and yield qualitatively consistent rotational spin-down trends across PMS and MS phases. However, at the solar age (4.57 Gyr), the same models over-predict the equatorial rotation rate (v = 4.72 kms-1 vs. 2.0 kms-1) and the mean surface magnetic field (B = 36.9 G vs. 1 G). These discrepancies reflect the omission of additional angular momentum loss mechanisms and possible oversimplifications in magnetic saturation physics. While the adaptive alpha-MLT converges to the solar-calibrated value ([1.76, 1.78]) at the present age, its variability during earlier phases significantly influences Li depletion. We compare predictions with observational data from 64 open clusters. The results demonstrate that incorporating time-dependent B and alpha-MLT improves Li predictions and captures rotational evolution trends, but cannot yet reproduce the present-day solar rotation and magnetic flux without additional physics. We discuss these limitations and outline future work for a more complete model of solar-type stars.
http://arxiv.org/abs/2603.09289 (pdf) (9)
As plans for the construction of third-generation gravitational wave (GW) detectors advance, research into strongly lensed GWs has become increasingly critical. It is anticipated that hundreds of multi-image lensed GWs will be detected annually. We present a comprehensive suite of lensed GW mock catalog derived from a composite lens mass model incorporating dark matter halos, galaxies, and subhalos. We analyze three source populations with four detector network configurations considering the earth rotation. Our simulations encompass not only conventional doublets and quadruplets but also subhalo-lensed events, highly magnified systems, and complete three or five image systems with a detectable central image, a feature distinct from optical lensing. For the joint ET+CE network, we forecast an annual detection rate of approximately 400 doublets and 36 quadruplets. Notably, this population includes roughly 107 events lensed by subhalos and 20 complete systems with detectable central images. Furthermore, we analyze high-magnification events (μ> 3), predicting approximately 360 such cases. Under a more relaxed selection criterion that requires only at least one lensed signal to exceed the detection threshold, we estimate a total of approximately 617 lensed events. We also investigate the impact of variations in lens mass models and stellar evolution models on event rates, as well as the distributions of SNR pairs and time delays. These results establish a more physically grounded statistical prior for the future identification and authentication of lensed GW signals. The Gravitational Waves-Lensing Mock Catalog (GW-LMC) have been made publicly available.
http://arxiv.org/abs/2603.08788 (pdf) (8)
The recent pericentric passage of the Large Magellanic Cloud (LMC) through the Milky Way (MW) has dislodged its centre of mass, inducing a state of dynamical disequilibrium, the reflex motion, in the kinematics of outer stellar halo stars. Using data out to distances of 160 kpc from the combined H3+SEGUE+MagE outer halo survey, we constrain the mass of the MW and LMC, as well as the resulting reflex motion and the velocity anisotropy of the stellar halo. Using a suite of 32,000 rigid MW–LMC simulations, each with a MW stellar halo evolved to the present day in the combined MW–LMC potential, we perform Simulation Based Inference by training a neural posterior estimator on the means and dispersions of the radial and tangential velocities of stars from the combined H3+SEGUE+MagE outer halo sample. Relative to halo stars at 100 kpc, we find the magnitude of the reflex velocity to be v_travel = 39.4^+7.6_-7.2 km s^-1. Simultaneously, we determine the enclosed MW mass to be M_MW(< 50 kpc) = 3.63 ± 0.16 × 10^11 M_⊙ and the enclosed LMC mass to be M_LMC(< 50 kpc) = 9.74^+2.07_-1.81× 10^10 M_⊙. Our results suggest that the total LMC mass must be at least ∼20
http://arxiv.org/abs/2603.09394 (pdf) (8)
Stochastic particle acceleration in turbulent plasmas plays a key role in shaping high-energy emission from relativistic outflows, such as those in Active Galactic Nuclei (AGN) and microquasars. While traditional Fermi-II models provide a foundational framework, they often oversimplify the complex nature of realistic magnetohydrodynamic (MHD) turbulence, especially in high-amplitude (δB/B_0 ∼ 1) and relativistic regimes. Recent plasma simulations for these conditions have revealed highly non-linear energization effects, such as sudden, large momentum jumps, that remain largely unexplored in astrophysical applications. We present a novel Monte Carlo framework STRIPE that models particle acceleration as a continuous-time random walk (CTRW), capturing both intermittent energy gains and radiative losses. The stochastic evolution of particle momenta is driven by jumps with random magnitudes determined by a distribution of magnetic-field-line velocity gradients, with synchrotron and inverse Compton cooling incorporated self-consistently. Using STRIPE, we explore particle acceleration under physical conditions characteristic of TeV-PeV γ-ray emitting microquasars recently identified by Large High Altitude Air Shower Observatory (LHAASO). We find that relativistic, high-amplitude turbulence naturally produces particle spectra with steep low-energy cutoffs, and hard extended power-law high-energy tails reaching tens of PeV. These features differ markedly from standard quasi-linear theory and are well suited to explaining the unexpectedly hard TeV-PeV spectra of LHAASO-detected microquasars. These results highlight turbulent acceleration in the relativistic regime as a promising mechanism for particle energization in microquasar systems, as well as potentially other extreme astrophysical environments.
http://arxiv.org/abs/2603.09179 (pdf) (8)
High-energy gamma-ray (>GeV) emission of gamma-ray bursts (GRBs) is very important in probing the jet evolution and particle acceleration of GRBs. The observations of high-energy photons are limited except for a few very bright GRBs, hindering precise measurements of the spectral and temporal evolutions of GRBs. Here we report the detection of high-energy gamma-ray emission up to 100 GeV with Fermi-LAT using a stacking analysis of a collection of 330 GRBs. High significance detection of the emission has been found, and the precise light curves and energy spectra can be measured. The light curves and time-resolved spectra of the sub-sample of 220 LAT individually detected GRBs can be well explained by the standard afterglow emission from a population of GRBs with both synchrotron and synchrotron self-Compton mechanisms, assuming a distribution of initial Lorentz factors. However, the emission of the relatively weak sample of the 110 LAT individually undetected GRBs cannot be well reproduced in the same framework, indicating the existence of possible energy injection effect in the GeV band for the first time. The observations hence provide new insights in understanding the high-energy emission of GRBs.
http://arxiv.org/abs/2603.09178 (pdf) (8)
The instrumental profile (IP) of a telescope is of great significance for spectroscopic analyses, especially for wavelength calibration and stellar parameter measurements. The Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) employs arc lamps for wavelength calibration. These lamps produce sharp emission lines with known wavelengths, and the observed arc lamp spectra can well characterize the IP. However, IPs are influenced by multiple factors, making them difficult to model accurately with traditional methods. Neural networks, which can automatically capture complex patterns and nonlinear features in data, provide a promising approach for high-precision IP measurement. We therefore construct a multi-layer perceptron (MLP) based on The Payne neural network to derive IPs for LAMOST. After training, the model can retrieve the IP for any fiber, at any wavelength, and at any time. We then apply the derived IP to stellar radial velocity (RV) measurements and analyze the impact of different IP center localization methods on the results. Finally, the dispersion of the measured RVs is reduced by approximately 3 km/s. This improvement will facilitate the search for long-period binary stars via RV variations.
http://arxiv.org/abs/2603.09876 (pdf) (7)
Polaris is a highly unusual Cepheid with observed properties that are difficult to reconcile with stellar evolutionary models. Since the initial detection of Polaris&apos magnetic field in 2020, we have conducted a magnetic monitoring campaign with the ESPaDOnS spectropolarimeter at the Canada-France-Hawaii Telescope. We compute Stokes V least-squares deconvolution profiles and measure the associated mean longitudinal magnetic field strengths ⟨ B_z⟩. The surface magnetic field has remained remarkably stable over five years of observations, with ⟨ B_z⟩ varying between approximately -3 G and +0.6 G. From the periodic modulation of ⟨ B_z⟩ we infer a stellar rotation period of P_rot=100.29±0.19 days. This is the first direct measurement of P_rot for a classical Cepheid. Previous interferometric radius measurements and P_rot imply an equatorial rotation velocity of v_eq=23.3±0.2 km s^-1. We set a conservative upper bound on the projected equatorial rotational velocity of v_eqsin i_⋆ < 13.5 km s^-1 and constrain the stellar inclination angle to be i_⋆<37^∘. Using the previously determined orbital solution, we find a high likelihood of a strong spin-orbit misalignment. We determine the lower bound on the obliquity angle between the stellar rotation and orbital axes to be β>18.7^∘ at 99% confidence. We discuss the challenges in interpreting the origin and properties of the surface magnetic field in the context of Polaris&apos uncertain evolutionary history and the merger hypothesis.
http://arxiv.org/abs/2603.09674 (pdf) (7)
Fast X-ray transients (FXTs) have been detected for over a decade, yet their origins are still enigmatic. The observed association between FXTs and broad-lined Type Ic supernovae (SNe Ic-BL) suggests that some may share the same progenitor with Long Gamma-Ray Bursts. In this work, we numerically simulate the long-term evolution of a relativistic jet propagating from inside the progenitor star up to the photon diffusion radius of the cocoon. Then we post-process the hydrodynamic results and calculate the cocoon cooling emission for various viewing angles from the jet axis. We find that, for viewing angles θ_v=10^∘–20^∘, the off-axis cocoon emission can produce FXTs with luminosity L_X≃ 10^47-48 erg s^-1 and duration t_X≃ 10-100s. The observed spectra are quasi-thermal with the peak energy E_peak≃0.8 keV. These properties naturally explain FXTs&apos observational features, including their high luminosity, soft spectra, and lack of gamma-ray counterparts. The Rayleigh-Jeans tail of the FXT spectra extends to the UV, producing an early UV flash simultaneously. As the cocoon expands and cools, the emission peak shifts to UV and optical bands, resulting in a bright optical plateau lasting for ∼1 day with color temperature T_UV/opt≃ (1-3)×10^4K, before the emergence of supernova emission. Although our model underpredicts the UV/optical luminosity at ∼1 day, it still provides useful diagnostics for identifying the origins of FXTs.
http://arxiv.org/abs/2603.09559 (pdf) (7)
Turbulence plays a critical role in the atmosphere, oceans, engineering, and astrophysics. The dissipation (heating) induced by turbulent flows is particularly important for the thermodynamics and chemistry of interstellar clouds, yet its structure and statistics remain poorly understood. Using high-resolution turbulence simulations with controlled explicit viscosity, we study the kinetic energy dissipation rate, ε_kin, across subsonic and supersonic regimes. We find that dissipation lags large-scale kinetic energy injection events by 1.64±0.21 and 0.48±0.07 turbulent turnover times in subsonic and supersonic turbulence, respectively. Correlations show ε_kin∝|∇×𝐯|^2 (vorticity squared) in the subsonic regime, where density fluctuations are negligible, while in the supersonic regime dissipation is primarily correlated with density, ε_kin^3/2. A spectral analysis demonstrates that achieving numerical convergence of ε_kin across all scales is challenging, especially in the subsonic case, even at 2048^3 resolution. Nonetheless, subsonic dissipation is clearly localised on small vorticity-dominated scales, while supersonic dissipation spans many scales, combining elongated, thin shocks with small-scale vorticity. Finally, we determine the fractal dimension of ε_kin. In the subsonic regime, intense dissipation is predominantly organised in flattened vortex filaments embedded in thin shearing layers on small scales, becoming more volume-filling at larger scales. In the supersonic regime, ε_kin exhibits a fractal dimension between 1 and 2 across nearly all scales, likely reflecting shock surfaces and their intersections forming filaments.
http://arxiv.org/abs/2603.09210 (pdf) (7)
The determination of unseen companion masses (M_1) is essential for identifying compact objects in binary systems, yet obtaining reliable orbital inclinations remains one of the most difficult challenges. In this study, we focus on ten targets selected from a sample of 89 compact object candidates characterized by large mass functions, rapid rotation, and high-quality Large Sky Area Multi-object Fiber Spectroscopic Telescope (LAMOST) spectra. We measure their projected rotational velocities (v sin i) from the LAMOST medium-resolution spectra and, combined with stellar radii, derive orbital inclinations and the corresponding companion masses. Our results show that five sources exhibit mass ratios M_1 / M_2 > 2/3, with no detectable spectral signatures of the unseen companions, providing strong evidence for their compact nature. Two particularly notable cases, J0341 and J0359, host companions with inferred masses of 1.39^+0.09_-0.10 M_⊙ and 1.34^+0.08_-0.09 M_⊙, respectively. These masses suggest that the invisible objects are either neutron stars or massive white dwarfs with masses close to the Chandrasekhar limit. If they are white dwarfs, these two targets are highly likely to be Type Ia supernova progenitors. This study highlights the potential of v sin i measurements as a systematic approach to unveiling compact objects in binaries.
http://arxiv.org/abs/2603.07714 (pdf) (7)
We investigate elastic scattering by a compact, horizonless body in curved spacetime, considering a massless scalar wave incident on a static, spherically symmetric, uniform-density star of radius R and mass M with a Schwarzschild exterior. We introduce an exact Debye-series decomposition of the scattering matrix, in the spirit of Debye expansions in Mie scattering. This decomposition separates direct surface reflection from contributions involving transmission into the interior and subsequent propagation, and admits a natural trajectory interpretation. We then determine the associated Regge-Debye pole spectrum in the complex angular-momentum plane. For neutron-star-like tenuities (R>3M), the spectrum exhibits two pole families: a surface-wave branch associated with the surface matching condition and a broad-resonance branch associated with the interior regularity condition. For ultracompact objects (R<3M), the surface-wave branch persists, while the interior-resonance sector splits into broad- and narrow-resonance branches. We next reconstruct the scattering amplitude from the Debye partial-wave contributions and find excellent agreement with direct partial-wave calculations. Finally, we develop complex angular-momentum representations order by order in the Debye series, making explicit how the pole families and non-pole sectors contribute to each Debye term. In the neutron-star-like regime, we find a genuine competition between Regge-Debye pole sums and branch-cut contributions, and show that, at high frequency, the rainbow-like enhancement already arises from the first interior-transmission contribution and is dominated by the interior-resonance Regge-Debye poles. By contrast, in the ultracompact regime, the Debye amplitudes are overwhelmingly pole dominated.
http://arxiv.org/abs/2603.09873 (pdf) (6)
Recent models and simulations of cluster formation within molecular clumps consider multi-scale, hierarchical accretion, which leads to clump mass growth over time. This mode of mass accumulation could have implications regarding the evolution of observable properties such as mass and radius, bringing into question the interpretation of commonly cited thresholds for high-mass star formation. In this paper, we use the conveyor belt model of cluster formation to create synthetic cores/clumps and derive physical and observational properties. We show that while this model successfully predicts many observed trends, modifications are required to match properties of high-mass prestellar clumps. When the model clumps are observationally classified as intermediate- or high-mass star-forming, the threshold delineating these two groups agrees with those found in the literature; however, results show that high-mass clumps at early evolutionary stages can be misclassified using standard surface density thresholds. Our logistic regression analysis reveals the quantity of material to ever enter a star-forming region is the most important factor in differentiating intermediate- and high-mass star-forming regions. This implies observations characterising the environment surrounding star-forming regions are crucial, especially at early evolutionary stages.
http://arxiv.org/abs/2603.09839 (pdf) (6)
Aims. We investigated the transport of solar energetic particles (SEPs) during the relativistic widespread event of 28 October 2021, quantifying the role of parallel and perpendicular diffusion and constraining the spatial extent of the injection region. Methods. We employed inverse modeling of particle focused transport and 2D numerical simulations including cross-field diffusion. Multi-spacecraft observations from STEREO-A, Solar Orbiter, and near-Earth spacecraft are used to reproduce particle intensity profiles and anisotropies across a wide range of electron and proton energies. Simulated flux profiles are compared across different heliolongitudes to derive consistent transport parameters. Results. The analysis yields parallel mean free paths within or slightly above the Palmer consensus range, and perpendicular mean free paths that correspond to ∼ 1–3
http://arxiv.org/abs/2603.09330 (pdf) (6)
The QCD Equation of State with N_f=3 massless quark flavours is determined non-perturbatively over a broad range of temperatures, extending from the electroweak scale down to 3 GeV, and smoothly connecting to the low-temperature regime. The comparison with perturbative predictions shows that, even at temperatures approaching the electroweak scale, the Equation of State can be accurately described only by adding terms beyond the known perturbative series, including non-perturbative contributions. The strategy that allows this investigation in the previously unexplored high-temperature regime combines shifted boundary conditions with a determination of the lines of constant physics based on the running of a non-perturbatively defined renormalized coupling. This methodology is general and can be applied to QCD with four or five massive quark flavours.
http://arxiv.org/abs/2603.08878 (pdf) (6)
The STIS echelle gratings can be used with a variety of different central wavelength settings. "Secondary" wavelength settings, designed to cover select absorption or emission lines, have not been calibrated as precisely as their primary mode counterparts. In particular, secondary echelle mode traces (and subsequent extraction regions) have been previously defined using straight line fits to each spectral order. In this work, we define a new general method for defining echelle traces that utilizes Gaussian process regression and accounts for the detailed curvature of each order across the detector. Across a variety of echelle grating and central wavelength settings, we find that this method can improve flux throughput by ∼4