Most interesting arXiv papers
Wednesday 28-01-2026
http://arxiv.org/abs/2601.18890 (pdf) (68)
We use the TNG50 galaxy formation simulation to generate mock JWST NIRCam and NIRSpec microshutter array (MSA) observations of Hα-emitting gas in M_*=10^8-10^11.5 M_⊙ star-forming galaxies at z=2-6. We measure morphological properties from the mock imaging through Sersic profile fitting, and gas rotational velocities (v) and velocity dispersions (σ) by fitting the mock spectra as thin, rotating discs. To test the efficacy of such simple parametric models in describing complex ionised gas kinematics, we compare the best-fit quantities to intrinsic simulation measurements. At z=3, we find that v and σ for aligned and resolved sources generally agree well with intrinsic measurements, within a factor of ∼2 and ∼1.5, respectively. The recovery of kinematics is robust for smooth, disc-like systems, but v and σ can be over- or underestimated by more than a factor of 2, respectively, for intrinsically elongated systems. The scatter in the recovery accuracy is larger at higher redshift, as TNG50 galaxies at z>3 deviate more strongly from the thin rotating disc assumption. Despite uncertain measurements for individual galaxies, we find that key population trends, such as the weak redshift evolution of σ and v/σ as well as the dependence of σ on the global star formation rate, are broadly recovered by our kinematic modelling. Our work provides the end-to-end framework needed to compare NIRSpec MSA observations to cosmological simulations and to quantify observational biases in measuring ionised gas kinematics, highlighting the need for the development of dedicated models for high-redshift galaxies.
http://arxiv.org/abs/2601.19725 (pdf) (46)
Aims. We investigate the impact of diffuse ionised gas (DIG) on the determination of emission line ratios and gas-phase metallicities in two local analogues of high-redshift star-forming galaxies: UM 462 and IIZw 40. Understanding how DIG affects these quantities is essential for interpreting unresolved observations of distant galaxies, where integrated spectra are often used to trace their chemical evolution. Methods. Using archival Very Large Telescope, Multi-Unit Spectroscopic Explorer (MUSE) data, we spatially resolved the warm ionised medium of both galaxies. We derived oxygen abundances through the direct method and several HII-based strong-line calibrators, and we used the Hα surface brightness (Σ(Hα)) to distinguish regions dominated by HII or DIG emission. Results. Oxygen abundances derived from the N2 and O3N2 indices show an inverse correlation with Σ(Hα), ionisation parameter, and EW(Hα), with DIG-dominated regions exhibiting higher 12+log(O/H) than the galaxy mean by ∼0.2 dex in UM 462 and ∼0.1 dex in IIZw 40. The metallicity differences between HII-dominated and DIG-dominated Σ(Hα) bins reach ∼0.4 dex and ∼0.3 dex in UM 462 and IIZw 40, respectively. The observed trends with Σ(Hα), metallicity, EW(Hα), and ionisation parameter indicate smoothly varying ionisation conditions rather than true abundance variations. These effects reflect different ionisation sources and levels, and can produce spurious metallicity gradients in galaxies with extended DIG structures, potentially mimicking signatures of metal-poor gas infall. In our sample, DIG ionisation is most likely dominated by photon leakage from H II regions, with additional contributions from feedback-driven shocks.
http://arxiv.org/abs/2601.18871 (pdf) (43)
Recent studies have revealed that bars can form as early as a few billion years after the Big Bang, already displaying similar characteristics of evolved bars in the Local Universe. Bars redistribute angular momentum across the galaxy, regulating star formation, AGN activity, and building new stellar structures such as nuclear discs. However, the effects of bar-driven evolution on young galaxies are not yet known, as no evidence of bar-built stellar structures has ever been found beyond z = 1, until now. In this work, we show evidence of a bar-built, star-forming nuclear disc, already present at redshift z = 1.5. This is the first evidence of a bar-built stellar structure at Cosmic Noon. We find that this nuclear disc is actively forming stars and has the same size as some nuclear discs in nearby galaxies. This evidence solidifies the now emerging picture in which bars are fundamental not only in the late evolution of galaxies, but also in their early evolutionary stages. It changes the current paradigm by urging a revision of our picture of galaxy evolution beyond redshift one, to include new considerations on the role played by bars as early as a few billion years after the Big Bang.
http://arxiv.org/abs/2601.19609 (pdf) (34)
We present a catalog of 971 FR-I radio galaxies (FR-Is) identified from the Very Large Array Faint Images of the Radio Sky at Twenty-Centimeters (FIRST) survey. The identifications were made using a hybrid method that combines deep learning with ridgeline flux density distribution analysis. Among these sources, 845 are new discoveries. The catalog comprises sources characterized by edge-darkened double jets, an absence of significant bent morphology, and angular sizes ranging from 23 to 159 arcseconds. Optical and/or infrared counterparts have been identified for 813 FR-Is. Among these, the host galaxies are predominantly (88.1
http://arxiv.org/abs/2601.19525 (pdf) (33)
We present a pilot study on the host galaxy environments of CHIME/FRBs by cross-matching baseband-localised events with the LOFAR Two-metre Sky Survey Data Release 2 (LoTSS DR2) at 144 MHz. Unlike traditional methods reliant on optical imaging, our radio-based selection allows for the identification of dust-obscured or optically faint star-forming galaxies. Of the 140 CHIME FRBs considered, 33 lie within the LoTSS DR2 footprint, and 16 show potential radio counterparts. Through multi-wavelength analysis, spectral energy distribution (SED) fitting, and redshift constraints from the Macquart relation, we identify two secure and one tentative host candidates, all consistent with active star formation. However, their Hα-derived star formation rates appear underestimated, likely due to significant dust attenuation, as suggested by infrared colours and compact optical morphologies. Our results highlight the value of low-frequency radio data in complementing optical host searches and demonstrate the feasibility of host identification even in the absence of optical confirmation. With forthcoming data from LoTSS DR3 and the full CHIME/FRB baseband release, this method offers a promising path toward statistically robust studies of FRB host galaxies and their environments.
http://arxiv.org/abs/2601.19842 (pdf) (30)
We present the first robust dynamical measurement of the supermassive black hole (SMBH) mass in the massive early-type galaxy NGC 4061 using high-spatial-resolution ALMA observations of the ^12CO(2-1) emission. By combining archival Cycle 6 data with new Cycle 7 observations, we achieve a synthesized beam of 0”.16 × 0”.13, comparable to the expected sphere of influence of the central black hole. The molecular gas forms a regularly rotating circumnuclear disk aligned with the prominent dust lane seen in HST imaging. We model the full three-dimensional ALMA data cube using the KinMS forward-modeling framework, exploring both data-driven and analytic prescriptions for the gas surface brightness distribution. Our Bayesian analysis yields a best-fitting SMBH mass of M_BH = (1.17^+0.08_-0.10 [ stat.] ± 0.43 [ syst.]) × 10^9 M_⊙ and an I-band stellar mass-to-light ratio of M/L_F814W = 3.46^+0.07_-0.06 [ stat.] ± 0.10 [ syst.] M_⊙/L_⊙. The inferred black hole mass is fully consistent across different modeling assumptions and remains insensitive to plausible radial variations in the M/L_F814W profile. Our results resolve the long-standing discrepancy between previous indirect mass estimates based on conflicting stellar velocity dispersion measurements and demonstrate that the exceptionally large dispersion reported in the literature is likely spurious. This study highlights the power of high-resolution ALMA molecular gas kinematics for precision SMBH mass measurements at the high-mass end of the local black hole mass function.
http://arxiv.org/abs/2601.19148 (pdf) (30)
In this paper of ERGO-ML (Extracting Reality from Galaxy Observables with Machine Learning), we develop a model that infers the merger/assembly histories of galaxies directly from optical images. We apply the self-supervised contrastive learning framework NNCLR (Nearest-Neighbor Contrastive Learning of visual Representations) on realistic HSC mock images (g,r,i - bands) produced from galaxies simulated within the TNG50 and TNG100 flagship runs of the IllustrisTNG project. The resulting representation is then used as conditional input for a cINN (conditional Invertible Neural Network) to gain posteriors for merger/assembly statistics, particularly the lookback time and stellar mass of the last major merger and the fraction of ex-situ stars. Through validation against the ground truth available for simulated galaxies, we assess the performance of our model, achieving good accuracy in inferring the stellar ex-situ fraction (≤± 10 per cent for 80 per cent of the test sample) and the mass of the last major merger (within ± 0.5 log for stellar masses >10^9.5 ). We successfully apply the TNG-trained model to simulated mocks from the EAGLE simulation, demonstrating that our model is applicable outside of the TNG domain. We use our simulation-based model to infer aspects of the history of observed galaxies, in particular for HSC images that are close to the domain of TNG ones. We recover the trend of increasing ex-situ stellar fraction with stellar mass and more spherical morphology, but we also identify a discrepancy between TNG and HSC: on average, observed galaxies generally exhibit lower ex-situ fractions. Despite challenges such as information loss (e.g. projection effects and surface brightness limits) and domain shifts (from simulations to observations), our results demonstrate the feasibility of extracting the merger past of galaxies from their optical images.
http://arxiv.org/abs/2601.18864 (pdf) (29)
JWST has revealed a population of compact galaxies in the early Universe with broad emission lines and strong Balmer breaks; among them the so-called ”little red dots” (LRDs). Their nature remains uncertain with hypotheses including exotic phenomena. We assemble a sample of LRD-like objects at z>3 and use self-consistent radiative-transfer calculations to show that a supermassive black hole accreting from a dense gas cocoon accurately reproduces the detailed spectra. We show that the cocoons must be non-spherical, with comparable amounts of inflowing and outflowing material. And we predict correlations between Balmer break strength, Balmer line-absorption and scattering line width, which we confirm in our observed sample. We reproduce all LRD-like properties without requiring star-like atmospheres and we determine the typical black hole in our sample to be of order a million solar masses, with ionized cocoon masses of tens of solar masses potentially supplied from a much larger cold-gas reservoir.
http://arxiv.org/abs/2601.19424 (pdf) (29)
The determination of distances using time-series photometry of Type Ia supernovae (SNe Ia) relies on a 5% empirical correction related to the properties of their host galaxies, e.g., global stellar mass. It is therefore crucial for unbiased cosmology inference that host galaxy properties be self-consistently determined across the full range of redshifts probed, which we undertake in this study for approximately 2000 SNe in the Union3 compilation (now Union3.1). We use aperture-matched, homogeneously-reduced, optical-infrared photometry from the DESI Legacy Imaging Surveys to derive global galaxy properties using the stellar population synthesis and SED-fitting code Prospector. We find that the host masses of z<0.10 SNe in Union3 were, on average, overestimated relative to the rest of the sample, while the opposite was true for z<0.15 SNe in Pantheon+. After correction, the two studies&apos average distance modulus estimated for low-redshift SNe, previously >0.03 mag discrepant, come into 0.01 mag agreement. We then update the UNITY SN analysis and find that the uncertainties on all standardization parameters shrink to 0.6-0.9x their previous sizes. For flat-ΛCDM, we find Ω_m=0.344^+0.026_-0.025, a -0.3σ shift from Union3. We then combine with measurements of Baryon Acoustic Oscillations (BAO) and the Cosmic Microwave Background (CMB) exactly as done by DESI DR2 and find w_0=-0.719±0.084, w_a=-0.95^+0.29_-0.26, corresponding to 3.4σ evidence against a cosmological constant (down from 3.8σ). We also update the DESI combined probe analysis using our correction to Pantheon+ and the recent DES-SN5YR Dovekie recalibration, finding 3.2σ (up from 2.8σ) and 3.4σ (down from 4.2σ) evidence against a cosmological constant in the w_0w_a plane, altogether marking a significantly improved consistency across SN analyses.
http://arxiv.org/abs/2601.19861 (pdf) (28)
Clusters of galaxies have been used to measure a subtle effect predicted by Einstein: gravitational redshift. This signal encodes pristine information about our Universe, since it is sensitive to the depth of the clusters&apos gravitational potential wells. In this work, we show how gravitational redshift can be used to test a fundamental physical principle: the weak equivalence principle. This principle stipulates that all matter falls in the same way in a gravitational potential, regardless of its nature. By comparing the amplitude of the gravitational redshift signal with the velocity dispersion in galaxy clusters, we build a novel test of this principle targeted to the unknown dark matter. Our test is sensitive to any additional interaction that would alter the way dark matter falls in gravitational potentials, hence leading to a violation of the equivalence principle. We show that currently available data can constrain the presence of a fifth force in clusters at the level of 7-14%, while the newest surveys will reach a precision of a few percents. This demonstrates the crucial role played by galaxy clusters in testing fundamental properties of dark matter.
http://arxiv.org/abs/2601.19594 (pdf) (27)
The JAGB method has been proposed as a distance indicator for Local Group galaxies. We investigate the populations of the J region in the (J-K,J) colour-magnitude plane of the Large and Small Magellanic Clouds (LMC and SMC), aiming to distinguish general features of the J luminosity function (JLF) from those sensitive to the stellar population of each galaxy. Using a population synthesis approach based on stellar evolution and dust formation modelling, we predict the distribution of stars within the J region and compare it with observations. Stars in this region are identified as recently formed C-stars that have not yet accumulated large amounts of carbon. Typically, 2–3 M_⊙ stars remain longer in the J region, while lower-mass stars evolve faster. The JLF of the LMC, peaked at the expected magnitudes for these stars, confirms this picture. In the SMC, the J flux distribution is shifted to higher magnitudes, indicating an older population with lower-mass progenitors.
http://arxiv.org/abs/2601.18885 (pdf) (26)
Galactic globular clusters (GCs) were born shortly after the Big Bang. For such old stellar systems the initial mass function (IMF) at the high mass regime can never be observed directly, because stars more massive than about 1 Mo have evolved since longtime. However, the hydrostatic to explosive alpha-element ratio (HEx ratio) offers a way to bypass the lack of observable high mass stars through the yields that massive stars released when exploding as supernovae, incorporated in the stars we presently observe in GCs. The HEx ratio measures the percentage of high mass stars over the total number of stars exploding as supernovae and it is an efficient probe of the ephemeral first phases of the GC evolution. We exploited a recently completed survey to assemble a dataset of very homogeneous abundances of alpha-elements in 27 GCs from [Fe/H] -2.4 to -0.3 dex. In agreement with previous results from APOGEE, we confirm that the HEx ratio is indistinguishable for GCs formed in situ and accreted in the Galaxy, and that this ratio decreases with increasing metallicity. However, we posit that this trend is better explained by a metallicity-dependent IMF deficient in the highest mass stars at high metallicity, as corroborated by the declining [O/Mg] ratio as a function of the [Mg/H] ratio. At odds with the previous analysis based on APOGEE data, we detect an anti-correlation of HEx ratio with both present day and initial GC masses. Finally, we hypothesise that in that analysis, the stars of the GC M 54 were probably confused with stars in the core of the Sagittarius dwarf galaxy, where the cluster is presently immersed.
http://arxiv.org/abs/2601.18866 (pdf) (25)
(Abridged) We present a systematic assessment of two major cloudy atmospheric model grids – SM08 (Saumon Marley 2008) and Sonora Diamondback – when applied to low-resolution near-infrared (0.8-2.5 μm) spectroscopy. Our analysis focuses on a uniform sample of 142 age-benchmark brown dwarfs and planetary-mass objects spanning late-M, L, and T spectral types, with independently determined ages from 10 Myr to 10 Gyr. We perform forward-model spectral fitting for all benchmarks&apos IRTF/SpeX spectra (R∼80-250) using both SM08 and Sonora Diamondback atmospheric models to infer effective temperatures, surface gravities, metallicities, radii, and cloud sedimentation efficiencies. The two model grids yield broadly consistent results. Among L4-L9 dwarfs, we identify a statistically significant, population-level age dependence of the cloud parameter f_sed, with young benchmarks (<300 Myr) exhibiting systematically lower f_sed values than older counterparts. This trend is absent across L0-T5 and T0-T5, demonstrating that cloud properties vary with age and surface gravity and offering explanations for the observed gravity-dependent photometric properties at the late-L end of the L/T transition. By comparing spectroscopically inferred parameters with predictions from evolution models, we quantify systematic errors in the fitted atmospheric parameters and establish empirical calibrations to anchor future studies using these atmospheric models. Stacked residuals of the sample reveal wavelength-dependent data-model mismatches associated with key atomic and molecular absorption bands, highlighting the need for improved opacities and rainout chemistry. Finally, we show that including an interstellar-medium-like extinction term significantly improves the spectral fits, confirming and broadening previous findings and suggesting missing opacity sources in current cloudy models.
http://arxiv.org/abs/2601.18986 (pdf) (24)
Using a combination of Hubble Space Telescope and James Webb Space Telescope imaging, a runaway supermassive black hole (RBH-1) was recently identified with an inferred velocity of 954^+110_-126 km s^-1, likely ejected from a compact star-forming galaxy (denoted as GX) at z ≈ 0.96. Assuming the runaway black hole was the outcome of the gravitational-wave-driven merger of two black holes, we use its measured runaway velocity together with gravitational-wave recoil predictions from numerical relativity and black hole perturbation theory to constrain the mass ratio and spin configuration of the progenitor SMBHs that overcame the final-parsec problem and merged ∼ 70 Myr ago. We find that the progenitor binary must have been precessing, with a mass ratio m_1/m_2≲ 6, and that the more massive SMBH must have possessed a high spin (dimensionless spin magnitude ∼ 0.75) in order to generate a recoil of this magnitude. This has important astrophysical implications as similar SMBH mergers can be an interesting source population for the upcoming LISA mission with signal-to-noise ratios ≳ 1000. Furthermore, the progenitor SMBH properties imply that GX was likely formed through a major, gas-rich (“wet”) merger between two galaxies of comparable mass, with a mass ratio ≲ 4.
http://arxiv.org/abs/2601.18954 (pdf) (24)
In the ΛCDM paradigm, stellar halos form through the accretion and disruption of satellite galaxies. We introduce new semi-analytic modeling within the SatGen framework to track the ex-situ stellar components of Milky Way–like galaxies across large ensembles of merger trees, enabling a statistical study of the stochastic nature of galaxy assembly. We find that accreted stellar halos are typically built by only a few progenitors and are highly sensitive to the fate of the most massive satellite, producing order-of-magnitude variations in accreted stellar halo mass even at fixed host halo mass. Different stellar components trace distinct phases of host halo growth: central and accreted stellar mass correlate most strongly with early assembly, while surviving satellites trace more recent accretion. Finally, using Random Forest Regression, we quantify how well observable galaxy properties can recover halo assembly histories, providing a framework for interpreting upcoming low-surface-brightness observations of stellar halos.
http://arxiv.org/abs/2601.18926 (pdf) (24)
We present multiwavelength observations and analysis of six luminous fast blue optical transients (LFBOTs) discovered in Zwicky Transient Facility (ZTF) survey data. We identified these LFBOTs from their fast light-curve evolution (t_1/2≤ 12d), blue colors at peak brightness (g-r≤-0.5mag), a visible host galaxy, high optical luminosity (M_g<-20), and an X-ray or radio detection. With the exception of AT2024aehp (ZTF24abygbss), these transients exhibit peaks in their 10GHz radio light curves at t_rest≈ 50-100 d, with peak radio luminosities ranging from 10^38-10^40 erg s^-1. Modeling the radio emission as synchrotron radiation indicates a fast (v=0.1-0.3c) shock in a dense (n_e≈10^3-10^4 cm^-3) medium. The X-ray emission varies by ≈2 orders of magnitude in luminosity (10^42-10^44 erg s^-1) at t_rest∼20d. Analysis of the host-galaxy photometry and spectroscopy for each transient shows that they are predominantly nonnuclear (a few kpc offset) with star-forming host galaxies of stellar masses 10^9-10^11 ,M_⊙. Unlike all other LFBOTs to date, AT2024aehp exhibited a luminous (M<-19mag) plateau in the optical light curve; spectra during this plateau phase showed a featureless blue continuum. The 6-15 GHz radio emission of AT2024aehp brightened by over an order of magnitude from t_rest≈70d to t_rest≈130d. The mostly consistent radio behavior between optically selected LFBOTs implies a similar circumburst medium, leading us to prefer a progenitor scenario in which mass is lost in a consistent way shortly prior to the terminal event, such as a massive star merging with a compact object.
http://arxiv.org/abs/2601.19854 (pdf) (22)
SNe Ia have been used to provide key constraints on the equation-of-state parameter of dark energy. They are generally standardized under the assumption that they belong to a single population, with luminosities standardized in a continuous (roughly linear) fashion using the observed light-curve timescale. We update the Union3+UNITY1.5 SN cosmology analysis in light of increasing evidence for at least two core populations of SNe Ia and apply this "UNITY1.8" model to the updated "Union3.1" compilation (Hoyt et al. 2026). In addition to finding evidence for two different light-curve-shape (x1) distributions, we also find that the color distributions are different, that the light-curve-shape/magnitude standardization relations are different, and that these populations have different distributions across host-galaxy stellar mass and redshift. Importantly, we find that the residual host-mass luminosity step found in prior SN Ia cosmology analyses is now consistent with zero for unreddened SNe. We report a significantly tightened constraint on the split in the red-color standardization between SNe in low- and high-mass galaxies. We find that the estimated uncertainties shrink on cosmological parameters when fitting the same SNe assuming two modes versus one mode. We confirm similar trends in simulated data when running both versions of UNITY on the same (two-mode) simulations. For a flat LambdaCDM cosmology, we find Om = 0.334+0.025-0.024 from SNe alone; for a flat w0-wa cosmology, we find w0 = -0.760+0.084-0.082 and wa = -0.79+0.28-0.30 when including SNe, BAO, and CMB. In the 2D w0-wa plane, adding SNe to BAO and compressed CMB increases the tension with flat LambdaCDM from 2.1 sigma to 2.6 sigma.
http://arxiv.org/abs/2601.19011 (pdf) (21)
We have observed six `dark&apos neutral hydrogen (HI) clouds discovered in the Virgo cluster by the Arecibo Galaxy Environment Survey (AGES) with the Karl G. Jansky Very Large Array (VLA), giving higher angular and velocity resolution than the original AGES observations. We detected compact HI emission in two of the sources, AGESVC1 231 and AGESVC1 274, allowing us to firmly associate them with faint (m_g > 18.5), blue (g-i < 0.1) optical counterparts with high M_HI/L_g ratios. In a further two sources, we detected low column-density extended HI emission, consistent with these being dispersing clouds from ram-pressure stripping or tidal interactions. The final two sources were not detected with the VLA, allowing us to set low column-density limits on the HI detected by AGES that are consistent with these clouds also being formed from HI that is dispersing into the intra-cluster medium. The four HI sources not associated with optical counterparts thus appear likely to be relatively short-lived objects. No evidence was found for either pressure-supported turbulent spheres or stable dark galaxies.
http://arxiv.org/abs/2601.18896 (pdf) (21)
The merger history of the Galaxy has been traced back firmly to redshift 2 (10 Billion years ago). While there have been claims of the existence of at least one more significant merger before this time, supporting evidence has been indirect and contentious. Here we show that the population of globular clusters around the Galaxy depicts three distinct age-metallicity sequences, one associated with the merger with Gaia-Enceladus 10 billion years ago, one to the progenitor of the Milky Way and a third intermediate sequence associated to at least one merger which we estimate took place merely 1.5 billion years after the Big Bang. This discovery has been possible thanks to exquisite Hubble Space Telescope data and sophisticated analysis that enables very precise relative age determination of globular clusters. The newly identified sequence reveals that this merger took place with an object of stellar mass similar to that of Gaia-Enceladus ( 5x10^8 M_⊙), and which deposited most of its mass in the inner 6 kpc of the Milky Way. The unambiguous identification of a third merger event in the inner Galaxy puts to rest earlier debates, and honoring previous work we name the progenitor system Low-energy-Kraken-Heracles, or LKH for short.
http://arxiv.org/abs/2601.19438 (pdf) (19)
Current and upcoming large-scale structure surveys are pushing toward increasingly wide angular coverage, where wide-angle effects (arising from the varying line of sight across the curved sky) become critical for accurate modeling of the three-dimensional galaxy power spectrum. At the same time, these survey&aposs broader redshift reach makes the effects of redshift evolution (beyond the effective-redshift approximation) non-negligible on large radial scales. Additional observational effects such as the survey window function and integral constraints also become significant on these large scales, necessitating a careful theoretical treatment to robustly constrain local primordial non-Gaussianities and relativistic effects. In this work, we present a consistent and accurate theoretical framework for modeling the commonly used power spectrum multipoles (PSM) on large scales using the discrete spherical Fourier-Bessel (dSFB) basis. This basis ensures numerical stability and allows an exact separation between angular and radial modes. Using the dSFB basis, we study the impact of wide-angle effects and redshift evolution on the PSM, and incorporate the effects of window function convolution and integral constraints. We validate our PSM modeling using lognormal mocks under radial integral constraints with realistic survey geometries, demonstrating the readiness of our framework for application to all-sky galaxy surveys.
http://arxiv.org/abs/2601.18862 (pdf) (19)
The redshifted 21-cm line from neutral hydrogen (Hi) enables volumetric intensity mapping of large-scale structure in the post-reionization Universe. In anticipation of &aposs wide redshift coverage and high signal-to-noise clustering measurements, we study the redshift-space 21-cm bispectrum and its spherical-harmonic multipoles as probes of anisotropic non-linear structure formation and departures from General Relativity. Using a tree-level perturbative description for the 21-cm brightness-temperature field in redshift space, and adopting the Hu–Sawicki f(R) model as a representative modified-gravity scenario, we forecast the detectability of configuration-dependent signatures with an –like survey. We derive the bispectrum-multipole covariance including sample variance and thermal noise and evaluate the expected signal-to-noise of deviations relative to ΛCDM. We find that the observable information is dominated by the lowest multipoles, while higher-order modes are strongly suppressed. This concentration in the lowest multipoles is well matched to sensitivity and to the quasi-linear modes that are expected to remain accessible in practice. The strongest modified-gravity sensitivity arises from squeezed and stretched triangle configurations on quasi-linear scales, where scale-dependent growth enhances the bispectrum relative to the total variance. Our results position 21-cm bispectrum multipoles as a practical, SKA-ready observable for testing gravity beyond ΛCDM in the post-reionization epoch.
http://arxiv.org/abs/2601.19866 (pdf) (19)
Hypervelocity stars (HVSs) are stars ejected from the Galactic Centre (GC) through tidal interactions with the central supermassive black hole. Formed in the immediate vicinity of Sgr A^∗, these stars are accelerated to velocities high enough to escape the GC and be observable in the Galactic halo. Using spectroscopy from the Dark Energy Spectroscopic Instrument (DESI) and astrometry from Gaia, we conducted a six-dimensional search for HVSs and identified a compelling candidate, hereafter DESI-312, whose bound trajectory can be confidently traced back to the GC. The star resides in the inner halo and exhibits supersolar metallicity ([Fe/H] = 0.27± 0.09), distinct from other known stellar populations with radial orbits. Its inferred GC ejection velocity of 698^+35_-27 is consistent with a Hills mechanism ejection, supporting an origin in the innermost regions of the Milky Way. We considered alternative origins for the star, including disk ejections from young clusters and globular clusters, but these scenarios fail to explain both its orbit and metallicity. Unlike previously identified A- and B-type HVSs, DESI-312 is a ∼ 1 M_⊙ star on the main sequence or early subgiant branch, thus enabling a detailed chemical analysis of its atmosphere and offering a rare window - unobscured by dust and crowding - into the composition of the central regions of the Galaxy.
http://arxiv.org/abs/2512.17539 (pdf) (19)
We present a spectroscopic analysis of XMM-Newton and NuSTAR observations of the &aposcomplex&apos NLS1 PG 1535+547 at redshift z=0.038. These observations span three epochs: 2002 and 2006 with XMM-Newton alone, covering the 0.3-10 keV energy range, and a coordinated XMM-Newton and NuSTAR observation in 2016, covering the 0.3-60 keV energy range. The X-ray spectra across all epochs exhibit both neutral and ionized absorption, along with reflection features from the accretion disc, including a prominent Compton hump in the broadband data. Notably, the spectral shape varies across epochs. Our analysis suggests this variability is attributed to changes in both line-of-sight absorption and the intrinsic emission from PG 1535+547. The source is obscured by multiple layers of partially and/or fully covering neutral and ionized absorbers, with neutral column densities ranging from undetectable levels in the least obscured phase to ∼0.3-5×10^23cm^-2 in the most obscured phase. A clear warm absorber is revealed during the least obscured phase. The continuum remains fairly consistent (Γ≈ 2.2±0.1) during the first two observations, followed by a substantial flux decrease (by a factor of ∼7 in the 2-10 keV band) in 2016 compared to 2006. The 2016 data indicates the source is in a reflection-dominated state during this epoch, with a reflection fraction of R>7 and an X-ray source located at a height ≤ 1.72r_g. Simultaneous fitting of the multi-epoch data suggests a rapidly rotating black hole with a spin parameter, a>0.99. These findings imply that strong light-bending effects may account for the observed continuum flux reduction.
http://arxiv.org/abs/2601.19522 (pdf) (18)
We present the first determination of the Galactic stellar mass function (MF) for low-mass stars (0.2-0.5 M_sun) at metallicities [Fe/H] < -1. A sample of 53,000 stars was selected as metal-poor on the basis of both their halo-like orbits and their spectroscopic [Fe/H] from Gaia DR3 BP/RP (XP) spectra. These metallicity estimates for low-mass stars were enabled by calibrating Gaia XP spectra with stellar parameters from SDSS-V. For -1.5 < [Fe/H] < -1, we find that the MF below 0.5 M_sun exhibits a "bottom-heavy" power-law slope of alpha -1.6. We tentatively find that at even lower metallicities, the MF becomes very bottom-light, with a near-flat power-law slope of alpha 0 that implies a severe deficit of low-mass stars. This metallicity-dependent variation is insensitive to the adopted stellar evolution model. These results show that the Galactic low-mass MF is not universal, with variations in the metal-poor regime. A further calibration of XP metallicities in the regime of M < 0.5 M_sun and [Fe/H] < -1.5 will be essential to verify these tentative low-metallicity trends.
http://arxiv.org/abs/2601.18965 (pdf) (16)
The MIRI Excess Around Degenerates (MEAD) Survey is a cycle 2 JWST program designed to image nearby white dwarfs with MIRI at 10 and 15 microns. This survey targeted 56 white dwarfs within 25 pc to search for mid-infrared excesses, flux deficits from collision-induced absorption, and resolved substellar companions. In this paper we present our analysis of WD 0644+025, an unusually massive white dwarf (0.95 Msun) and the MEAD target exhibiting the most significant mid-infrared excess. The observed JWST MIRI photometry shows a 7.3 sigma excess at 15 microns and a 3.6 sigma excess at 10 microns, which may be associated with either a planetary companion or a circumstellar dust disk. This excess corresponds to a companion mass of 6.8 Mjup (Teff=261 +/- 9 K) with orbital distance <11.8 au, although substantially lower masses are possible if we consider a closely orbiting insolated companion. No spatially resolved sources are detected within 200 au, with contrast curve analysis excluding planets more massive than 2 Mjup beyond 12 au. Metal pollution is confirmed in both archival Keck HIRES spectra from 1999 and new observations from 2025, with no evidence suggesting the accretion rate has substantially changed over the decades. We explore possible dust disk morphologies to describe the observed IR excess, and find that traditional debris disks struggle to fit our data. WD 0644+025 thus represents a compelling case study in the growing population of white dwarfs with cold infrared excesses, and highlights JWST&aposs ability to probe planetary system remnants inaccessible to prior infrared observatories.
http://arxiv.org/abs/2601.19480 (pdf) (16)
The origin of accretion disc winds remains disputed to date. High inclination, dipping, neutron star Low Mass X-Ray Binaries (LMXBs) provide an excellent testbed to study the launching mechanism of such winds due to being persistently accreting and showing a nearly ubiquitous presence of highly-ionised plasmas. We aim to establish or rule out the presence of a wind in the high inclination LMXB 4U 1624-49, for which a highly ionised plasma has been repeatedly observed in X-ray spectra by Chandra and XMM-Newton, and a thermal-radiative pressure wind is expected. We leverage the exquisite spectral resolution of XRISM to perform phase-resolved spectroscopy of the full binary orbit to characterise the highly ionised plasma at all phases except during absorption dips. An outflow is clearly detected via phase-resolved spectroscopy of the source with XRISM/Resolve. Based on analysis of the radial velocity curve we determine an average velocity of 200-320 km/s and a column density above 10^23 cm^-2. The line profiles are generally narrow, spanning from 50 to 100 km/s, depending on the orbital phase, pointing to a low velocity sheer or turbulence of the highly ionised outflow and a potential increase of turbulence as the absorption dip is approached, likely due to turbulent mixing. The line profiles, together with the derived launching radius and wind velocity are consistent with a wind being launched from the outskirts of the disc and without stratification, pointing to a thermal-radiative pressure origin.
http://arxiv.org/abs/2601.19150 (pdf) (16)
We numerically study the Q-ball formation triggered by a cosmological first-order phase transition within the Friedberg-Lee-Sirlin model. By performing lattice simulations, we track the nonequilibrium dynamics throughout the transition, providing a precise description of the Q-ball formation mechanism and the resulting mass spectrum. Collapsing false-vacuum regions first form thermal balls, which subsequently cool via dissipative interactions and stabilize into long-lived Q-balls with nonzero spin. We observe a large population of low-mass Q-balls, as well as rare, massive Q-balls that are several times larger than the analytical prediction. The final Q-ball population exhibits a broad mass spectrum spanning over two orders of magnitude, characterized by an exponential tail of number density at large masses. The simulations suggest that the Q-ball abundance is approximately 50
http://arxiv.org/abs/2601.19863 (pdf) (15)
Recent simulations have identified long-lived “prompt cusps” – compact remnants of early density peaks with inner profiles ρ∝ r^-3/2. They can survive hierarchical assembly and potentially enhance signals of dark matter annihilation. In this work, we incorporate prompt cusps into the semi-analytic substructure framework SASHIMI, enabling a fully hierarchical, environment-dependent calculation of the annihilation luminosity that consistently tracks subhalos, sub-subhalos, and tidal stripping. We assign prompt cusps to first-generation microhalos and propagate their survival through the merger history, including an explicit treatment of cusps associated with stripped substructure. We find that the substructure hierarchy converges rapidly once a few levels are included, and that prompt cusps can raise the total annihilation boost of Milky-Way–size hosts at z=0 to B∼ O(10) for fiducial cusp-occupation assumptions, compared to a subhalo-only baseline of B_sh∼few. Across a wide range of host masses and redshifts, prompt cusps increase the normalization of B(M_host,z) while largely preserving its mass and redshift trends. Compared to universal-average, peak-based estimates, our fiducial boosts are lower by about an order of magnitude, primarily reflecting a correspondingly smaller inferred cusp abundance in host halos, highlighting the importance of unifying peak-based cusp formation with merger-tree evolution and environmental dependence.
http://arxiv.org/abs/2601.18876 (pdf) (14)
Spectroscopic surveys now supply precise stellar label measurements such as element abundances for large samples of stars throughout the Milky Way. These element abundances are known to correlate with orbital actions or other dynamical invariants. We present a new data-driven method for empirically measuring the circular velocity curve of the Galaxy that uses element abundance gradients in the plane of radial kinematics. We use stellar surface abundances from the APOGEE survey combined with kinematic data from the Gaia mission. Our results confirm the ordered structure of the Milky Way disk in terms of average [Fe/H] and [Mg/Fe] abundance ratios, and suggest that ⟨[Fe/H]⟩ traces the radial position of stars in the disk, while ⟨[Mg/Fe]⟩ traces the orbital excursions around this radius. Our method uses the radial orbit structure in the Galaxy to enable an empirical measurement of the circular velocity curve, epicyclic and azimuthal frequencies, and kinematic gradients across the Milky Way disk. From these measurements, we infer a value of the circular velocity curve at the Solar radius of v_c,⊙ = 235.3^+2.8_-3.7 km s^-1 using the most constraining abundance ratio, [Mg/Fe]. We also measure the radial and azimuthal frequencies for a circular orbit at the solar radius, κ_0,R_⊙=36.9^+0.8_-1.0 km s^-1 kpc^-1 and Ω_0,R_⊙=28.5_-0.1^+0.4 km s^-1 kpc^-1, respectively. These values lead to an estimate of the Oort constants of A = 16.5^+0.1_-0.1 km s^-1 kpc^-1 and B=-11.9^+0.1_-0.3 km s^-1 kpc^-1. We measure the radial acceleration at the Solar radius to be (∂ Φ/∂ R)_⊙ = a_R_⊙=7.0^+0.2_-0.1 pc Myr^-2.
http://arxiv.org/abs/2601.19363 (pdf) (13)
Analyses of global magnetic fields in M dwarfs rely on many approximations regarding the derivation of average line profiles from spectropolarimetric data, interpreting them with analytical functions and modelling them using Zeeman Doppler imaging (ZDI). These assumptions have not been systematically tested. We assessed the accuracy of standard treatments of average polarisation profiles in M dwarfs and their interpretation with ZDI. We focused on the filling-factor approach, which attempts to represent coexisting global and small-scale fields. We performed polarised radiative transfer calculations across the near-infrared spectrum of a typical M dwarf. From these theoretical spectra, we derived mean Stokes profiles and approximated them with different line-synthesis methods. To test the recovery of global fields, we performed ZDI inversions using simulated Stokes V observations for low- and high-activity cases. The analytical approximation of mean polarisation profiles reproduces Stokes I and V only for fields up to 1 kG and fails for linear polarisation. ZDI with single-line analytical Stokes V profiles is adequate for weakly magnetic M dwarfs with fields below a few hundred gauss. However, combined with the filling-factor formalism, this traditional modelling approach produces unphysical local fields and distorted global geometries for active M dwarfs with multi-kilogauss fields. These issues are mitigated using a new mapping technique based on theoretical Stokes profiles that account for both global and randomly distributed small-scale fields. Our study reveals fundamental limitations of current ZDI analyses of active M dwarfs and questions the reliability of some published maps. (abridged)
http://arxiv.org/abs/2601.18878 (pdf) (13)
The FLAMES survey “Na-O anti-correlation and HB" uncovered the modern standard for globular clusters (GCs), that is their ubiquitous multiple stellar populations (MPs) distinct by the abundance of proton-capture elements. That survey can still be mined to extract a wealth of data. We derive new abundances of Mg, Si, Ca, and Ti for 948, 954, 1542, and 1350 red giant branch stars in 16 GCs, both formed in situ or accreted in the Milky Way. The program GCs cover the metallicity range from [Fe/H]=-2.35 dex to [Fe/H]=-0.74 dex. Both the halo and disc GCs show a clear overabundance of alpha-elements with the modulation in Mg and Si due to the MPs phenomenon in different clusters. We found star to star variations in Si abundance correlated to changes in Na in more than half of our sample, implying that temperatures in excess of about 65 MK were achieved in the polluters responsible for the enrichment. We confirm with an enlarged sample the previous result that significant variations in Mg are observed in GCs that are metal-poor, massive or both. Evidence of excess of Ca with respect to reference unpolluted field stars are found in NGC 6752 and NGC 7078, indicating the action of proton-capture reactions at very high temperature regime in these GCs. These excesses fit very well in a previously found relation as a function of a combination of cluster mass and metallicity shown by other typical signatures of MPs. At odds with previous results based on the Si abundance from APOGEE, we found that the average abundance of alpha-elements is not an efficient discriminating factor between in situ and accreted GCs.
http://arxiv.org/abs/2601.19450 (pdf) (13)
The redshifted 21 cm signal from the cosmic dawn and Epoch of Reionization (EoR) encodes important information about both astrophysical processes and primordial physics, such as inflation. In this work, we use morphological statistics to explore the sensitivity of the 21 cm signal to inflationary features and EoR dynamics simultaneously. Focusing on primordial features from particle production during inflation we generate semi-numerical simulations of the 21 cm signal across redshifts 5 < z < 35, incorporating these features. Using Minkowski Functionals (MFs), we analyze the morphology of 21 cm fields: density, neutral hydrogen fraction, spin temperature, and brightness temperature. We demonstrate that MFs are highly sensitive to both the amplitude and scale of primordial features, capturing rich morphological information. In particular, we show that MFs can robustly identify inflationary features and distinguish them from the standard model. We further explore various EoR scenarios, and demonstrate that combining MFs across redshifts can disentangle the signatures of primordial features from EoR effects. This approach opens new avenues for probing inflation with upcoming 21 cm surveys.
http://arxiv.org/abs/2601.19018 (pdf) (13)
Calcium-strong transients (CaSTs) are a subclass of faint and rapidly evolving supernovae (SNe) that exhibit strong calcium features and notably weak oxygen features. The small but growing population of CaSTs exhibits some aspects similar to thermonuclear supernovae and others that are similar to massive star core-collapse events, leading to intriguing questions on the physical origins of CaSTs. SN 2025coe is one of the most nearby CaSTs discovered to date, and our coordinated multi-wavelength observations obtained days to weeks post-explosion reveal new insights on these enigmatic transients. With the most robust NIR spectroscopic time-series of a CaST collected to date, SN 2025coe shows spectral signatures characteristic of Type Ib SNe (SNe Ib, i.e. He-rich stripped-envelope SNe). SN 2025coe is the third X-ray detected CaST and our analysis of the Swift X-ray data suggest interaction with 0.12 ± 0.11 M_⊙ of circumstellar material (CSM) extending to at least 2 × 10^15cm (∼ 30,000 R_⊙), while our analysis of the 1-240 GHz radio non-detections gives an outer radius of that CSM of at most ∼ 4× 10^15 cm. This inferred nearby high-density CSM extending out to 3± 1 ×10^15 cm is similar to that seen in the other two X-ray detected CaSTs, and its presence suggests that either intensive mass-loss or some polluting mechanism may be a common feature of this subclass. Our work also expands upon recent studies on the optical properties of SN 2025coe and explores our current understanding of different progenitor systems that could possibly produce CaSTs.
http://arxiv.org/abs/2601.18893 (pdf) (12)
The radio galaxy M87 is well known for its jet, which features a series of bright knots observable from radio to X-ray wavelengths. The most famous of these, HST-1, exhibits superluminal motion, and our analysis of Chandra data <cit.> reveals a correlation between the X-ray flux of HST-1 and its separation from the core. This correlation likely arises from moving shocks in the jet, allowing measurement of the internal structure of HST-1 in the X-ray band. To follow up on these results, we use observations from the Hubble Space Telescope Advanced Camera for Surveys HRC/WFC/SBC channel and the Wide Field Camera 3 (WFC3)&aposs UVIS to analyze the image and flux variability of HST-1. Our analysis includes 245 ACS and 120 WFC3 observations from 2002-2022, with a total exposure time of ∼345 ks. We study the brightness profile of the optical jet and measure the relative separation between the core and HST-1 for comparison to the X-ray. We find that the X-ray and the UV/optical emission could arise from physically distinct regions. The measured proper motion of the knot HST-1 is 1.04±0.04 c from 2002-2005 and around 2.1±0.05 c from 2005-2022. We discuss the acceleration of the jet and the flaring synchrotron emission from HST-1 from optical to X-rays.
http://arxiv.org/abs/2601.19353 (pdf) (12)
During the course of its mission, ESA&aposs Gaia spacecraft has generated a map of the stars of the Galaxy of exquisite detail. While in its L2 orbit, the satellite has been exposed to high energy cosmic rays and solar particles, that caused permanent damage to its CCDs. The main effect of radiation damage on Gaia data is the distortion of its images and spectra, caused by the CCDs charge transfer inefficiency (CTI) during the readout process, that, if not taken into account, can result in inaccurate measurements of a star&aposs location and flux. In this work, the impact of CTI in the serial readout direction, larger than in the parallel due to the presence of CCDs manufacturing defects, has been analysed and modelled. A pixel-based, physically motivated CTI model, CtiPixel, has been developed to characterise the damage in Gaia CCDs. The model has been calibrated using dedicated serial CTI diagnostic data, taken every 3-4 months over the course of the mission. The model is shown to be a good representation of the observed signatures of CTI in the calibration datasets, and its parameters reveal significant insights into the nature of the CCD defects generated by space irradiation. The evolution of the damage in the serial direction shows a general small linear increase over time, with sudden step changes after strong solar flares and coronal mass ejections directed towards Earth. The serial CTI showed a further step increase as a consequence of the engineering CCD annealing experiment carried out after the completion of Gaia science observations.
http://arxiv.org/abs/2601.18910 (pdf) (12)
We evaluate uncertainty calculations in the calstis pipeline for data in the low-count regime. Due to the low dark rate and read-noise free nature of MAMA detectors, observations of UV-dim sources can result in exposures with 0 or 1 counts in some pixels. In this regime, the "root-N" approximation widely used to calculate uncertainties breaks down, and one must compute Poisson confidence intervals for more accurate uncertainty calculations. The CalCOS pipeline was updated in 2020 to account for these low-count uncertainties. Here, we assess how STIS observations are currently affected by this phenomenon, describe a new Jupyter notebook exploring the issue, and introduce a new utility, stistools.poisson_err, to manually calculate Poisson confidence intervals for 1D STIS spectra. Additionally, we describe a related software bug in the stistools.inttag utility, which splits TIME-TAG data into sub-exposures. This newly fixed bug serves as a useful case-study for the proper use of Poisson confidence intervals.
http://arxiv.org/abs/2601.19584 (pdf) (11)
The migration of low-mass planets is tightly controlled by the torques exerted by both gas and solids in their natal disks. While canonical models assume a solar solid-to-gas mass ratio (epsilon=0.01) and neglect the back-reaction of solid component of the disk, recent work suggests that enhanced metallicity can radically alter these torques. We quantify how elevated metallicities (epsilon=0.03 and epsilon=0.1) modify the gas and solid torques, test widely used linear scaling prescriptions, and identify the regimes where solid back-reaction becomes decisive. We performed global, 2D hydrodynamic simulations that treat solid material as a pressureless fluid fully coupled to the gas through drag and include the reciprocal back-reaction force. The planet was maintained on a fixed circular orbit, thus we computed static torques. The Stokes number was varied from 0.01 to 10, three surface-density slopes (p=0.5, 1.0, and 1.5) and three accretion efficiencies (eta=0, 10, and 100%) were explored. Torques, obtained by rescaling canonical epsilon=0.01 results, were compared with direct simulations. Solid torques scale linearly with epsilon, but gas torques deviate by 50-100% and can even reverse sign for St<=1 in epsilon=0.1 disks. These are due to strong, feedback-driven, asymmetric gas perturbations in the co-orbital region, amplified by rapid planetary accretion. Solid back-reaction in high-metallicity environments can dominate the migration torque budget of low-mass planets. Simple metallicity rescalings are therefore unreliable for St<=2, implying that precise migration tracks - particularly in metal-rich disks – require simulations that fully couple solid and gas dynamics. These results highlight metallicity as a key parameter in shaping the early orbital architecture of planetary systems.
http://arxiv.org/abs/2601.19268 (pdf) (11)
Despite the assumption that an ideal FLRW observer is not dependent on the local environment, observations are biased by the positions of the observers due to the matter correlations in the large-scale structure (LSS) of the universe. The variation of the mass distribution of the LSS of the universe implies that observers residing in different locations may suffer bias in their measurements when they look at the images of distant galaxies. Here, we assess the influence of the local environment on weak gravitational lensing (WL) shear statistics in the context of relativistic N-body code, . We derive numerical constraints on the cosmological parameters from the WL shear angular power spectrum and comment on the local environment&aposs influence on WL shear. We find tighter constraints on the parameter Ω_m above redshift z = 0.2, which implies over this redshift the local environment&aposs impact is minor. We also investigate the bispectrum and conclude that on average the impact of the local environment on f_NL (a measure of non-Gaussianities) is minimal and consistent with zero effect. However, we find that within the assembly of all possible observers/locations, there will also be a few that could infer the parameter f_NL of the order 10. These results could thus be used to estimate the uncertainty in the inference of cosmological parameters such as f_NL based on WL shear bispectrum and thus may have implications for future surveys requiring precision at the percent level.
http://arxiv.org/abs/2601.19412 (pdf) (10)
Quasars accretion models predict a broad range of optical and ultraviolet properties that depend primarily on black hole mass and accretion rate. Yet, most optically selected luminous quasars display strikingly similar continuum spectra. We show that this uniformity can be explained by a nearly constant luminosity to mass (Eddington) ratio, L_EDD and by thermal emission from a standard, optically thick, geometrically thin accretion disc. A standard disk with an Eddington ratio L_EDD=0.1 reproduces both the black hole mass/luminosity distribution of Sloan Digital Sky Survey (SDSS) quasars and their principal continuum properties. In this framework, the spectral energy distribution peaks beyond the observable ultraviolet range for nearly all sources. We show that the few quasars, expected to be cold enough to shift the peak into the observable region, indeed show this behaviour. This scenario is further supported by an analysis of the relation between the luminosity of the main broad emission lines and the continuum luminosity (i.e. the Baldwin effect). We find that 1) the observed slopes of the line to continuum relations match the expectations from the standard disk model, if we assume that the line emission is a good proxy of the ionizing luminosity; 2) the dispersions of the line-continuum luminosity relations are very small (as small as 0.13 dex), suggesting that the physics of the disk-broad line region system is dominated by only one parameter (the black hole mass) with a nearly constant Eddington ratio. Finally, we notice that our hypothesis of constant L_EDD=0.1 provides a black hole mass estimate (based on the observed luminosity) with a smaller error than the virial estimate.
http://arxiv.org/abs/2601.19135 (pdf) (10)
Particle acceleration in relativistic collisionless shocks remains an open problem in high-energy astrophysics. Particle-in-cell (PIC) simulations predict that electron acceleration in weakly magnetized shocks proceeds via small-angle scattering, leading to a maximum electron energy significantly below the Bohm limit. This upper bound manifests observationally as a characteristic synchrotron cutoff, providing a direct probe of the underlying acceleration physics. Gamma-ray burst (GRB) afterglows offer an exceptional laboratory for testing these predictions. Here, we model the spectral evolution of GRB afterglows during the relativistic deceleration phase, incorporating PIC-motivated acceleration prescriptions and self-consistently computing synchrotron and synchrotron self-Compton emission. We find that low-energy bursts in low-density environments, typical of short GRBs, exhibit a pronounced synchrotron cutoff in the GeV band within minutes to hours after the trigger. Applying our framework to GRB 190114C and GRB 130427A, we find that current observations are insufficient to discriminate between PIC-motivated acceleration and the Bohm limit, primarily due to large uncertainties in the Fermi-LAT band. Nevertheless, future MeV-TeV afterglow observations can break model degeneracies and place substantially tighter constraints on particle acceleration in relativistic shocks.
http://arxiv.org/abs/2601.19531 (pdf) (9)
Only a handful of gas giant planets orbiting white dwarfs are known. It remains unclear whether this paucity reflects observational challenges or the consequences of stellar evolution. We aim to carry out population synthesis of substellar objects around white dwarfs to predict the fraction and properties of white dwarfs hosting substellar companions. We generated a representative population of white-dwarf progenitors with substellar companion and used the stellar-evolution codes MESA and SSE with standard prescriptions for mass loss and stellar tides to predict the resulting population of white dwarfs and their companions. We find that the predicted fraction of white dwarfs hosting substellar companions in the Milky Way is, independent of uncertainties related to initial distributions, stellar tides, or stellar mass loss during the asymptotic giant branch, below 3%. The occurrence rate peaks at relatively low-mass ( 0.53 Msun to 0.66 Msun) white dwarfs and relatively young ( 1-6 Gyr) systems, where it exceeds 3%. The semimajor axes of the surviving companions range from 3-24 au. We estimate that 95% of the predicted companions are gas-giant planets. Owing to the strong dependence of companion occurrence on the metallicity of the white dwarf progenitor, the assumed age-metallicity relation strongly affects the predictions. Based on recent estimates of the local age-metallicity relation, we estimate that the fraction of white dwarfs with companions close to the Sun might reach 8%. If the planetary and brown dwarf companion distributions derived from intermediate-mass giant stars through radial velocity surveys reflect the characteristics of the true population, less than 3% of white dwarfs host substellar companions. This most likely represents an upper limit on possible detections because a significant number of companions might not be detectable with current facilities.
http://arxiv.org/abs/2601.19382 (pdf) (8)
Context: The Necklace nebula is a bipolar, post-common-envelope planetary nebula, the central star of which has been shown to have a dwarf carbon star companion. Aims: We aim to understand the origins of the Necklace and its dwarf carbon central star. Methods: We study the carbon abundance of the nebula through far ultraviolet spectroscopy obtained with the Hubble Space Telescope. Furthermore, through simultaneous modelling of multiband light and velocity curves, we attempt to constrain the parameters of the central star system. Results: Puzzlingly, we find that the region of the inner nebula observed with the Hubble Space Telescope is seemingly not carbon-rich, at odds with the dwarf carbon star nature of the companion of the central star. The initial mass of the nebular progenitor was likely very close to the limit to become carbon-rich, perhaps experiencing a very late thermal pulse. The dwarf carbon star companion is found to be significantly inflated with respect to that expected for an isolated main sequence star of the same mass. Conclusions: The properties of the central binary are consistent with the progenitor having become carbon-rich and its companion having accreted a significant amount of that carbon-enriched material. However, it is unclear how this evolutionary hypothesis can be reconciled with the inner nebula potentially being carbon poor.
http://arxiv.org/abs/2601.19386 (pdf) (8)
We investigate the possibility of constraining primordial black holes (PBHs) with masses M_PBH≳ 10^15 g through Galactic diffuse synchrotron emissions. Due to Hawking radiation, these types of PBHs are expected to be stable sources of cosmic-ray (CR) electrons and positrons with energies below 𝒪(10 MeV). In many CR propagation models with diffusive re-acceleration characterized by a significant Alfvén velocity V_a∼𝒪(10) km/s, the energies of the evaporated electrons/positrons can be further enhanced to 𝒪(100) MeV through their scattering with the Galactic random magnetic fields. Consequently, the observation of Galactic synchrotron emissions at frequencies above ∼ 20 MHz can provide useful constraints on the abundance of PBHs. Using the AMS-02 and Voyager-1 data on the boron-to-carbon nuclei flux ratio, we confirm that a significant Alfvén velocity V_a ∼ 20 km/s is favored in several benchmark diffusive re-acceleration models. We show that, in this scenario, the observed low-frequency synchrotron emissions (from 22 MHz to 1.4 GHz) can provide stringent constraints on PBH abundance. The obtained conservative constraints are stronger than those derived from the Voyager-1 all-electron (electron plus positron) data by more than one order of magnitude for M_PBH≳ 1× 10^16 g, and also stronger than our previous constraints derived from the AMS-02 positron data for M_PBH≳ 2× 10^16 g.
http://arxiv.org/abs/2601.19188 (pdf) (8)
We have developed X-ray telescope mirrors using an original electroforming replication technique established through the fabrication of millimeter-aperture, ultra-short-focal-length nanofocusing mirrors for synchrotron X-ray microscopy. This paper presents detailed results of X-ray illumination tests of a 60-mm-diameter, full-circumference, double-reflection monolithic electroformed nickel mirror and its Mirror Module Assembly (MMA). The experiments were conducted at the 1-km beamline BL29XUL at SPring-8. To simulate a parallel X-ray beam from celestial sources, we constructed a dedicated evaluation system, the High-Brilliance X-ray Kilometer-long Large-Area Expanded-beam Evaluation System (HBX-KLAEES). Owing to the high photon flux and the quasi-point-like source with a small divergence provided by HBX-KLAEES, the imaging performance was evaluated with high fidelity, resolving both the sharp core and large-angle components of the Point Spread Function (PSF). The results show an extremely sharp core with a Full Width at Half Maximum (FWHM) of 0.7 arcsec and a Half Power Diameter (HPD) of 14 arcsec, even after integration into the MMA. In addition, a positive correlation was found between angular resolution and axial figure error in both the primary and secondary mirror sections, indicating that axial figure errors contribute to image degradation. Based on these results, the MMA was selected as one of the hard X-ray optics for the FOXSI-4 sounding rocket experiment, which performs high-resolution soft and hard X-ray imaging spectroscopy of solar flares and was successfully launched. These results demonstrate the potential for further improvements in angular resolution and the development of high-resolution, ultra-short focal length X-ray optics for small satellites, including CubeSats.
http://arxiv.org/abs/2601.19123 (pdf) (8)
Extreme mass-ratio inspirals (EMRIs) in relativistic accretion discs are a key science target for the upcoming LISA mission. Existing models of disc-EMRI interactions typically rely on crude dynamical friction or Newtonian planetary migration prescriptions, which fail to capture the relativistic fluid response induced by the binary potential. In this work we address this gap by providing the relativistic calculation. We apply standard methods from self-force theory, black hole perturbation theory, and relativistic stellar perturbation theory to perform the full fluid calculation of the relativistic analogue of planetary migration for the first time. We calculate the response of a fluid in the perturbing potential of an EMRI consistently incorporating pressure effects. Using a master enthalpy-like variable and linearised fluid theory, we reconstruct the fluid perturbations and relativistic spiral arm structure for a range of spin values in the Kerr geometry. We conclude by deriving a relativistic torque-balance equation that enables computation and comparison of local torques with advected angular momentum through the disc. This opens a promising route towards establishing torque-balance relations between integrated disc torques arising from fluid perturbations and the forces acting on EMRIs embedded in matter.
http://arxiv.org/abs/2601.18905 (pdf) (8)
We present 1066 speckle measurements of M dwarf multiples observed over 2021-2024, all taken with HRCam on the Southern Astrophysical Research 4.1 m telescope. Among these, 900 observations resolve companions in 212 pairs, with separations spanning 17 milliarcseconds to 3.4 arcsec and brightness differences ranging from 0 to 4.9 magnitudes in the I filter. We have characterized the orbits of 54 of these companions, spanning periods of 0.67-30 yr, by combining our data with literature astrometry, radial velocities, and, in four cases, Hipparcos-Gaia accelerations. Among the orbits presented here are 28 that are the first-ever such characterizations for their systems, and 26 that revise previously-published orbits, thus providing a significant update to the observed dynamics of M dwarfs in the solar neighborhood. From these orbits, we provide new and updated dynamical total masses for these systems, precise to 0.7-7% in nearly all cases. Future mass derivations for components in these systems will contribute to efforts in refining the mass-luminosity relation for the smallest stars, and will enhance investigations of age, magnetism, and metallicity effects on luminosities at a given mass.
http://arxiv.org/abs/2601.19463 (pdf) (7)
We discuss key aspects of the nature of radiation from global strings and its impact on the relic axion density. Using a simple model we demonstrate the dependence on the spectrum of radiation emitted by strings. We then study the radiation emitted by perturbed straight strings paying particular attention to the difference between the overall phase of the field and the small perturbations about the string solution which are the axions. We find that a significant correction is required to be sure that one is analyzing the axions and not the self-field of the string. Typically this requires one to excise a sizeable region around the string - something which is not usually done in the case of numerical field theory simulations of string networks. We have measured the spectrum of radiation from these strings and find that it is compatible with an exponential, as predicted by the Nambu-like Kalb-Ramond action, and in particular is not a “hard” spectrum often found in string network simulations. We conclude by attempting to assess the uncertainties on relic density and find that this leads to a range of possible axion masses when compared to the measured density from the Cosmic Microwave Background, albeit that they are typically higher than what is predicted by the Initial Misalignment Mechanism. If the decay is via a “soft spectrum” from loops produced close to the backreaction scale we find that m_a≈ 160 μ eV and a detection frequency f≈ 38 GHz. If axions are emitted directly by the string network, and we use emission spectra reported in field theory simulations, then m_a≈ 4 μ eV and f≈ 1 GHz, however this increases to m_a ≈ 125 μ eV and f≈ 30 GHz using our spectra for the case of an oscillating string. In all scenarios there are significant remaining uncertainties that we delineate.
http://arxiv.org/abs/2601.19374 (pdf) (7)
Motivated by recent experimental refinements of stellar reaction rates, we establish a non-perturbative Green&aposs function formalism based on the exact solution of the Dyson equation for sub-barrier proton-nucleus resonant scattering. By utilizing bare Green&aposs functions to map the quantum tunneling problem onto a scattering formalism, we demonstrate that the summation of infinite quantum paths recovers the exact tunneling coefficients, enabling an analytical solution of the Dyson equation where the strong nuclear force is modeled as a surface delta-shell impurity embedded within the Coulomb field. Applying this framework to the astrophysically relevant p + ^7Li, p + ^14N, and p + ^23Na systems, we achieve precise agreement with experimental resonance energies while revealing a fundamental physical distinction in resonance formation. The heavier ^23Na system is identified as a saturated state, residing on a geometric plateau where the resonance energy becomes insensitive to the interaction strength; our calculated value of 2.11 MeV aligns remarkably well with the experimental level of 2.08 MeV. In contrast, the lighter ^7Li and ^14N systems emerge as threshold states in a weak-coupling window, where the resonance energy is highly sensitive to the potential parameters and is sustained near the continuum edge. In this regime, our model yields energies of 0.489 MeV and 1.067 MeV, closely reproducing the experimental benchmarks of 0.441 MeV and 1.058 MeV, respectively. We demonstrate that these threshold states are characterized by a significant enhancement of the resonant cross-section, driven by the inverse relationship between the tunneling width and the spectral density peak.
http://arxiv.org/abs/2601.19328 (pdf) (7)
The study of astronomical transients at high energies provides insights into some of the most extreme physical events in the universe; however, carrying out their detection and fast follow-up studies are limited by instrumental constraints. Search for Transient Object in New observations using Known Sources (STONKS) is a near-real-time transient detection system for XMM-Newton offering the capability to detect transients in XMM-Newton observations at fainter fluxes than can be achieved with wide survey instruments. We present the transients detected with the STONKS pipeline found in an XMM-Newton multi-year heritage survey of the Galactic plane to identify and classify highly variable X-ray sources that have recently been reported in this region. We examined the alerts created by the STONKS pipeline from over 200 XMM observations of the Galactic plane, screening for instrumental effects. The 78 alerts associated with 70 real astrophysical sources were then subjected to further temporal and spectral analysis. From the 70 sources we identified, we were able to classify 32 with a high degree of confidence, including 7 X-ray binaries, 1 γ-Cas analogue, and 1 magnetar candidate. Of the 70 sources, 23 were detected for the first time in X-rays. This systematic analysis of publicly available data has shown the value and potential of STONKS in the application to XMM-Newton observations. It will enable the community to detect transient and highly variable sources at fainter fluxes than with any other X-ray transient detection systems.
http://arxiv.org/abs/2601.19696 (pdf) (7)
Weak-gravitational lensing tends to wash out scale and time-dependent features of the clustering of matter, such as the Baryonic Acoustic Oscillations (BAO) which appear in the form of wiggles in the matter power spectrum but that disappear in the analogous lensing C_ℓ. This is a direct consequence of lensing being a projected effect. In this paper, we demonstrate how the noise complexity – often deemed "erasing the signal" – induced by a particular de-projection technique, the Bernardeau-Nishimichi-Taruya (BNT) transform arXiv:1312.0430, can be used to extract the BAO signal and non-gaussian aperture-mass-like properties at chosen physical scales. We take into account parts of the data vectors that should effectively be without cosmological signature and also introduce an additional re-weighting designed to specifically highlight clustering features – both at the probe (summary statistics) or map (amplitude of the field) level. We thus demonstrate why weak-gravitational lensing by the large-scale structure of the Universe, though only in a tomographic setting, does not erase scale and time-dependent features of the dynamics of matter, while providing a tool to effectively extract them from actual galaxy-shapes measurements.