2026-03-02 | | Total: 75
We report the first unambiguous detection of refractive substructure in an active galactic nucleus (AGN) using ground-based Very Long Baseline Interferometry (VLBI). Our analysis of TXS 2005+403 - observed at 1--5 GHz along a line of sight through the Cygnus region - reveals clear signatures of turbulence-induced substructure on long baselines that cannot be explained by the smooth scatter-broadened profile from diffractive effects alone. This signal persists across multiple observations spanning 2010-2019, demonstrating stable scattering properties along the line of sight. The combination of high flux density, compact intrinsic structure, and strong scattering establishes TXS 2005+403 as an exceptional laboratory for probing Galactic turbulence. This detection demonstrates that AGNs can serve as cosmic lighthouses illuminating interstellar plasma across the sky, complementing pulsar scintillation studies and informing scattering mitigation for millimeter-wavelength imaging of Sagittarius A*.
Intense solar eruptions occasionally accelerate solar energetic particles (SEPs) and can trigger ground-level enhancements (GLEs). Among the 77 known GLEs, the first four GLEs, #1 -- 4 in the 1940s took place before the advent of the standard neutron monitors and were missing from the International GLE Database. This data gap challenged their quantification. To overcome this difficulty, we systematically gathered, digitised, and quantified contemporaneous cosmic-ray records pertaining to these GLEs. These data allow us to reconstruct the temporal evolution, with the 1 -- 15 min resolutions, of these GLEs, and broaden their geographical coverage to a global scale. GLEs #1 and #3 exhibited gradual increases in their rise times, measured at 45 +/- 15 and 105 +/- 15 min, respectively. In contrast, GLEs #2 and #4 both exhibited abrupt increases of 15 +/- 15 min. We also compared integral ionisation increase on the standard ionisation chambers and their local geomagnetic cutoff rigidities Pc to qualitatively compare these GLE's spectral hardness: Our result indicates that their spectra are extremely hard for GLEs #2 and #4 and mildly hard for GLEs #1 and #3. GLE #3 showed the greatest integral ionisations for polar detectors among them.
The RADiation-hard Electron Monitor (RADEM) aboard the Jupiter Icy Moons Explorer (JUICE), launched on 14 April 2023, measures high-energy protons and electrons during the cruise phase and will continue throughout the nominal mission. Initial in-flight observations could not be explained by pre-flight ground calibration, motivating an in-flight calibration campaign. We calibrated the RADEM sensors using galactic cosmic rays by progressively increasing detector thresholds, thereby modifying their response to high-energy particles. Threshold-dependent in-flight count rates were compared with theoretical expectations derived from the Badhwar-O'Neill 2020 galactic cosmic ray model and corresponding response functions. These results were used to derive new in-flight calibration coefficients and to develop a flux reconstruction algorithm based on the bow-tie method. In several cases, the in-flight calibration slopes differ by up to an order of magnitude from ground calibration values. Proton fluxes from solar energetic particle events reconstructed with this method agree within a factor of two with independent measurements from the Solar and Heliospheric Observatory. These results demonstrate that RADEM provides accurate proton flux measurements in interplanetary space and is well suited for both single-spacecraft analyses and coordinated multi-mission studies of solar energetic particles. While electrons were clearly detected during the JUICE Lunar-Earth gravity assist, reliable reconstruction of their fluxes requires further analysis.
We examined the long-term behavior of the superbias calibration frames for the Advanced Camera for Surveys Wide Field Channel (ACS/WFC) aboard the Hubble Space Telescope (HST). Superbias frames are used to remove detector-level bias structure from science images and are currently generated after an anneal and delivered monthly. The primary goal of this study was to determine whether the frequency of superbias generation could be reduced without compromising calibration quality, potentially aligning with the Wide Field Camera 3 UVIS (WFC3/UVIS) approach of generating only one superbias per year. We analyzed superbias frames produced from 2007 through 2024 to investigate whether these calibration products have changed significantly over time, and whether the frequency of superbias generation and delivery could be safely reduced without loss of calibration accuracy. In addition to visual inspections and pixel-level comparisons, we employed Principal Component Analysis (PCA) to evaluate whether any long-term, global structure exists beneath the apparent noise in these frames. Our findings show that the superbias structure has remained fairly stable post-Servicing Mission 4 (SM4), a 15-year period, and no significant or unexpected global trends or systematic shifts were detected. However, due to unstable hot columns and increasing readout dark observed in ACS/WFC data, it is likely that these calibrations still benefit from more frequent superbias updates than the annual cadence adopted for WFC3/UVIS.
Being the most massive known young stellar cluster in the Milky Way, Westerlund 1 (Wd1) constitutes an ideal benchmark for understanding the evolution of massive stars. However, the cluster age remains highly controversial (~4-10 Myr), hindering the use of Wd1 as a reference for massive star evolution. One of the main issues is high foreground extinction, which has so far prevented the detection of the main sequence. Using infrared spectroscopy we seek to detect the cluster's main sequence for the first time, to characterise the Hertzsprung-Russell diagram, and to use the cluster's turn-off to obtain a robust age estimate. We obtained multi-epoch, near-infrared VLT/KMOS spectroscopic observations of Wd1 to map its population of massive stars. The spectra of ~110 members were analysed with CMFGEN models to derive stellar parameters, populate the cluster Hertzsprung-Russell diagram, and compare it with isochrones from evolutionary models. Our observations returned 47 new spectroscopically identified cluster members, with spectral types O9-B1 III-V. The cluster turn-off indicates an age of 5.5+/-1.0 Myr at a distance of 4.23+0.23-0.21 kpc, displaying a moderate degree of coevality. We demonstrate that our estimate of the age of Wd1 is robust against reasonable changes in the distance and extinction law, and the adopted rotational velocity and metallicity of the stellar isochrones. We further find that ~65% of the OB stars with multi-epoch coverage exhibit radial-velocity variability. Infrared observations of the unevolved stellar population support a single episode of star formation with an age of ~5.5 Myr, reinforcing its potential as a benchmark for massive star evolution and providing a reference sample for future binary population studies.
We present a broadband study of the WLQ SDSS J101353.45+492758.1, which displays a nearly featureless UV-optical spectrum with only a weak Mg II line alongside an exceptionally low X-ray flux. We model its spectral energy distribution using the relativistic thin-disk model kerrbb with a power law, and the multicomponent AGN model relagn, a physically motivated extension of agnsed incorporating warm and hot Comptonizing regions. Our fits constrain the black hole mass, accretion rate, X-ray loudness, and coronal energetics. Both approaches yield consistent BH masses of M_{BH} \approx 2 \times 10^{9} M_\odot and an Eddington accretion rate of \dot m \approx 0.1. The relagn fit including a warm Comptonizing region provides a significantly improved representation of the UV-soft X-ray continuum. The warm corona, characterized by kTe \simeq 0.20 keV, Γ \simeq 3.8, and an optical depth τ \simeq 7.26, extends to \sim 34 R_g. The hot corona appears compact and energetically suppressed, leading to an intrinsically weak X ray output with log(L_{X}/L_{bol}) \simeq -4.29, among the lowest reported for WLQs. The α_{ox} \sim 2.06 indicates the source to be in high/soft AGN spectral state. The combination of a luminous, standard disk and extremely weak hot corona suggests that this quasar hosts a highly inefficient inner coronal region. This explains its X-ray faintness and extreme deficit of high-ionization emission lines. The source may represent an AGN analog in "ultrasoft" accretion state, or a system in which the ionizing continuum is suppressed by a compact or quenched corona. Our study suggests that the source is not accreting at high Eddington ratio, highlighting the physical diversity of WLQs, and supports the view that geometric and radiative effects jointly shape their extreme spectral properties.
One key piece of evidence for dark matter is the rotation-curve problem: the disagreement between measured galactic rotation curves and their luminous mass. A novel solution to this problem is presented here, in a model that predicts observed Doppler-shifted spectra based only on the luminous matter estimates and one free model parameter. This model is applied to fit the rotation curves of the SPARC sample of 175 galaxies, yielding mass-to-light ratios, goodness of fit measurements, and the free parameter. The model's average reduced chi square compares favorably with the dark matter model for the same data, and more galaxies are successfully fit by this model. The model provides a useful formulation linking luminous matter to the observed rotation curves, with the dark matter contribution to galaxies encoded in two transformation terms of the luminous mass. It also offers a lower-parameter characterization of the rotation curve problem, and a power law relationship between the model's free parameter and galactic photometric quantities is observed, potentially removing the need for the free parameter.
Shock waves generated during cluster mergers offer a powerful probe of how large-scale structure grows and evolves in the Universe. As part of the MeerKAT-South Pole Telescope (SPT) survey, we report the discovery of a single arc-like radio relic in the galaxy cluster Abell 4067 ($z=0.099$), one of the lowest-mass clusters known to host such a structure. MeerKAT UHF-band (0.58--1.09 GHz) observations reveal a relic with a largest linear size of $\sim 1.48 \pm 0.02$ Mpc, located at a projected distance of 0.95 Mpc from the cluster centre. XMM-Newton X-ray data show that the relic's position and orientation relative to the intracluster medium (ICM) elongation are consistent with a merger-driven shock-wave scenario. The relic has an estimated radio power of $3.10 \pm 0.03 \times 10^{24}$ W Hz$^{-1}$ at 150 MHz. When placed in the $P_{150\,\mathrm{MHz}}$--$M_{500}$ scaling relation, the Abell 4067 relic appears less luminous compared to relics in more massive clusters, suggesting an association with weak merger shocks. This finding supports the idea that relics in low-mass clusters may form through less energetic merger events, leading to weak merger shocks. This is further supported by the absence of a detectable central radio halo in Abell 4067, reinforcing the idea that luminous radio halos are not a universal outcome of cluster mergers and highlighting the role of cluster mass, merger energetics, and evolutionary stage in shaping diffuse radio emission in the intracluster medium.
We forecast that the Terahertz Intensity Mapper (TIM) cross-correlated with Euclid's Fornax deep field (EDF-F), TIM$\times$EDF-F, will detect the [CII]-galaxy cross-power spectrum at a median redshift of 1.1 with $\gtrsim 7 σ$ confidence. The Poisson component of the cross-power spectrum at $0.1 \leq k \leq 10$ hMpc$^{-1}$ (i.e. cross-shot noise) will be detected at $\gtrsim 3 σ$ in 4 bins spanning $0.5 < z< 1.7$. This measurement will constrain the mean [CII] specific intensity over half of cosmic history and assess the degree to which Euclid-selected galaxies account for the [CII] intensity observed by TIM. We find that TIM can detect the cross-power spectrum across a wide range of [CII] intensity models.
Using the POLISH2 polarimeter at the Lick Observatory Shane 3-m and Nickel 1-m telescopes, we discover rotation phase-locked variations in linear polarization to be common among asteroids and a NEO in a clear, 383 to 720 nm bandpass. Essentially all bodies in our eight-year study harbor statistically significant, repeatable linear polarization variations at the $0.01\% = 100$ ppm level or above (1 Ceres, 2 Pallas, 3 Juno, 4 Vesta, 6 Hebe, 7 Iris, 12 Victoria, 15 Eunomia, 16 Psyche, 132 Aethra, 216 Kleopatra, and 65803 Didymos). Since polarimetry is a differential technique, such variations cannot be due to shape changes and must be caused by heterogeneity in surface albedo and/or composition. While (4) Vesta has long been known to exhibit large, repeatable polarization variations across its surface, we discover the variations on (6) Hebe, (12) Victoria, and (65803) Didymos to be 1.5 to 3.5 times as large. We hypothesize that the \textit{DART} impact with Dimorphos blanketed Didymos with depolarizing ejecta, which suggests pristine variations across Didymos to have been even larger. As the only NEO in this study with data quality sufficient to investigate polarization variations, Didymos' huge variations suggest they may be common among NEOs. We also discover optical circular polarization to be enhanced for low-albedo, M type asteroids, which is correlated with their large radar albedos. Thus, we present optical circular polarimetry as an alternative method for the identification of metalliferous bodies.
When in situ solar energetic electron (SEE) events are closely associated with nonthermal flares, the escaping electron population is frequently observed to be much smaller than the nonthermal-radiation-emitting population near the solar surface. If a single accelerated population drives both signatures, the physical mechanism causing this severe deficit of upward-propagating electrons remains poorly understood.Focusing on one of the 2022 November 10--12 SEE events associated with recurrent solar jets and interplanetary type III radio bursts, we present a new, combined microwave--X-ray analysis using the Expanded Owens Valley Solar Array (EOVSA) and the Spectrometer/Telescope for Imaging X-rays (STIX) aboard Solar Orbiter. This synergy enables, for the first time for such an event, spatially resolved diagnostics over a broad energy spectrum of the near-Sun energetic electrons, complemented by in situ measurements made by spacecraft at multiple heliocentric longitudes and distances. Consistent with earlier results based on in situ and X-ray data, our results show that only 0.1--1\% of energetic electrons escape into interplanetary space. Crucially, the new microwave spectral imaging analysis suggests that energetic electrons are strongly concentrated in a compact region just above a mini-flare arcade at the base of the jet spire, and that their number density decreases by at least two orders of magnitude in the direction of the jet spire away from this region. This steep gradient, revealed by the microwave diagnostics, points to efficient local acceleration and trapping in the region analogous to the above-the-looptop ``magnetic bottle'' region in major eruptive flares, allowing only a small fraction of electrons to access open magnetic field lines and enter interplanetary space.
We present a detailed study of HUDFJ0332.4-2746.6 (HUDF46), a $z \approx 1.84$ overdensity in the Hubble Ultra Deep Field, previously identified with HST as a proto-cluster. JWST/NIRISS spectroscopy expands its confirmed membership from 18 to 37 galaxies, while deep HST/ACS, JWST/NIRCam, and JWST/MIRI imaging provide a comprehensive multiwavelength view from the rest-frame UV to the mid-infrared. This dataset probes the population across three dex in stellar mass ($M_\bigstar \approx 10^{7.5\text{--}10.5}\,M_\odot$), delivering the first direct view of a young cluster down to such low-$M_\bigstar$ at $z\gtrsim1$. Assuming virialization, we derive a velocity dispersion of $σ\approx 670\pm 91\,\mathrm{km\,s^{-1}}$ and a halo mass of $M_{200} \approx (1.2\pm0.2) \times 10^{14}\,M_\odot$, in agreement with X-ray constraints from deep {\it Chandra} data. Despite residing in a massive halo likely in the hot-halo regime, the population is overwhelmingly star-forming, with no established red sequence and no extended X-ray emission from a hot intracluster medium. HUDF46 members have stellar and structural properties nearly indistinguishable from coeval field galaxies, and the structure hosts only one AGN candidate, found in its brightest galaxy, which lies at the cluster center. Overall, HUDF46 appears to be in a transitional phase prior to the onset of environmental quenching, making its galaxy population a key benchmark for tracing the processes that will later build a passive population and shape the assembly of massive clusters at later cosmic times.
Going beyond second order in weak lensing (WL) statistics is known to break degeneracies among cosmological parameters. We take a step further here, investigating whether higher-order statistics (HOS) in weak lensing can disentangle among General Relativity (GR) and modified gravity (MG), also taking into account the presence of massive neutrinos. To this end, we rely on mock convergence maps obtained from GR and $f(R)$ gravity N - body simulations, and we look for MG signatures in a wide set of higher-order WL probes. We rely on different metrics to quantify the discriminatory power of each probe, also varying the measurement setup. We find out that WL HOS can indeed disentangle MG and GR also in the presence of massive neutrinos.
Across the peak of cosmic star formation at $z\sim1-2$, inflow, processing, and feedback drive rapid changes in the spatial distribution and chemical composition of baryons in galaxies and surrounding reservoirs; this baryon cycle can be tomographically mapped by line intensity mapping (LIM) of atomic hydrogen, ionized carbon, and carbon monoxide. We present a simulation-based forecasting framework for detecting auto- and cross-power spectra between spectroscopic surveys of four such tracers at $z\sim0.5-1.7$ mapping the same deep field - TIM, EoRSpec/FYST, MeerKAT, & Euclid. We forward-model 3-D distributions for these tracers from magnetohydrodynamic simulations, directly capturing the two-halo, one-halo, and shot statistics without relying on analytical decompositions. We further detail a signal-to-noise formalism, tailored to LIM surveys with highly anisotropic geometries and Fourier-space coverage. We demonstrate that galaxy cross-correlations will be the dominant discovery channel for current-generation surveys. These instruments will detect the auto-spectra for CO and HI 21 cm and the CO $\times$ 21 cm cross-spectrum at modest S/N $\sim 1-10$, while placing upper limits on the [CII]-158$μ$m signals. [CII], CO, and HI LIM will be $\sim3-30\times$ ($0.5-1.5$ dex) more sensitive to cross-correlation with the Euclid survey, however, than their respective auto-correlations, constraining all three models of line emission at high significance (S/N $\sim 10-40$) within this decade. Finally, we formulate a staged instrumental trajectory with planned or reasonable improvements, including the as-proposed SKA-Mid. We forecast advancing the per-$k$-mode sensitivities of each auto-, galaxy-line, and line-line spectrum by several orders of magnitude, enabling new percent- and sub-percent level constraints on cosmology and the redshift evolution of star formation and the baryon cycle.
The evolution of the global solar magnetic field directly impacts the interplanetary magnetic field (IMF). During the solar maximum of Cycle 24, the monthly averaged IMF strength doubled over five Carrington rotations in late 2014. To understand the physical origin of this increase, we investigate the temporal evolution of open magnetic flux resulting from the emergence and decay of bipolar magnetic regions (BMRs). Using surface flux transport and potential field source surface models, we simulated how BMR characteristics, spatial distributions, and interaction with background magnetic fields affect open flux evolution. Our simulation confirmed that the relative configuration of BMRs can either inhibit open flux expansion via closed loops or promote it through favorable connections. The increase in open flux is primarily driven by the equatorial dipole component, which is enhanced by differential rotation acting on tilted BMRs. These behaviors suggest that large open field structures develop from equatorial dipole components formed by these stretched BMRs. We attribute the rapid IMF increase in 2014 (Carrington rotations 2152-2157) to the combination of the following three factors: (1) a specific sunspot configuration that facilitated the expansion of the southern coronal hole; (2) the emergence of a giant sunspot group (active region 12192) with high magnetic intensity; and (3) the diffusion of these regions, which reinforced the global magnetic field. These results imply that rapid open flux variations during solar maximum are governed not only by the characteristics of emerging BMRs but also by their interaction with pre-existing large coronal holes.
The reionization process is driven by ionizing photons from dwarf galaxies in halos with virial temperature $T_{\rm vir} \gtrsim 10^4$ K, while minihalos whose $T_{\rm vir}\lesssim 10^4$ K consume ionizing photons and have negative contributions to reionization. Since ionizing sources and minihalos have different clustering characteristics, not only the reionization history, but also the morphology of the ionization field, is sensitive to the small-scale power spectrum. If the power spectrum at small scales is enhanced compared with the standard six-parameter $Λ$CDM model, then both the sources and sinks of ionizing photons would be boosted and the net impact depends on the competition between them. Therefore, the 21 cm signal that can probe the morphology of the ionization field will be a useful tool for detecting the small-scale power spectrum. Using the power spectrum proposed by Cielo et al. (2025) (C25) as a demonstration, we investigate the influence of small-scale power spectrum enhancement on the ionization field and the 21 cm signal. We find that for the C25 model, even under the constraints of observed UV luminosity functions for high-$z$ galaxies and reionization history, the 21~cm power spectrum and the bubble size distribution could be still significantly different from the regular $Λ$CDM model. The upcoming SKA-low AA* telescope, and a further imaging telescope, have the potential to detect the small-scale power spectrum more deeply.
Relativistic jets from Active Galactic Nuclei (AGN) are highly energetic and emit radiation across a wide range of frequencies. Despite several observational studies, their particle composition still remains a key open question. The detection of high-energy neutrinos from blazar sources such as TXS 0506+056 has highlighted the plausibility of hadronic/lepto-hadronic models for AGN jets. To understand the origin of high-energy neutrinos from such sources, it is imperative to capture the complex interplay between the jet dynamics, their composition, and the mechanism of particle acceleration and cooling in relativistic jets. In this pilot study, we have coupled a numerical multi-zone framework for lepto-hadronic modeling, with 3D relativistic magneto-hydrodynamic simulations of AGN jets, including external photon fields. Our framework provides synthetic multi-wavelength and neutrino flux by spatially sampling the simulated jet into multiple zones. We investigate the implications of such a framework in exploring the different intrinsic and extrinsic pathways for proton-enrichment in jets. Essentially, we find that for low proton-to-electron number density ratios, producing a substantial jet neutrino flux, requires the underlying proton energy distribution to have a relatively flat spectrum with a power-law index of $\simeq 2.0$. We further find that while intrinsic shocks triggered by kink-instabilities in the jet can accelerate electrons to high energies, they may not be sufficient to produce such flat particle energy distributions for the chosen set of parsec-scale jet parameters. Finally, to produce a significant jet neutrino emission, our simulations suggest the need to consider particle acceleration mechanisms through alternative pathways, either internal or external.
The origin of the X-ray plateau phase observed in a large fraction of gamma-ray burst afterglows remains debated. We present a novel analytic framework for the hydrodynamics of ultra-relativistic, radially stratified outflows interacting with an external medium. By explicitly accounting for a continuous distribution of Lorentz factors within the ejecta, we derive analytic expressions describing the evolution of a long-lived, mildly relativistic reverse shock and determine its crossing time. Then, we compute the resulting synchrotron emission from both the forward and reverse shocks. The forward shock naturally produces a shallow, long-lasting X-ray decay consistent with the observed properties of X-ray plateaus (including the Dainotti relation). We further show that reproducing the observed plateau durations requires the stratified ejecta to extend to Lorentz factors of order $\gtrsim 100$, consistent with the ultra-relativistic outflow that powers the prompt $γ$-ray emission. The reverse shock generates a long-lived millimeter emission component that outshines the forward shock emission at these wavelengths. Both the plateau and reverse shock emission terminate smoothly once the slowest ejecta are processed, marking a transition to the standard Blandford-McKee self-similar evolution without requiring late-time energy injection or an additional emission component. Such stratified outflows are expected on physical grounds, as the ultra-relativistic ejecta responsible for the prompt $γ$-ray emission are unlikely to be launched with a single Lorentz factor. This model provides a unified picture in which the same outflow powers the prompt emission, the X-ray plateau, and the subsequent afterglow evolution.
We present observations of the classical T Tauri star DO Tau collected with the near-infrared SPIRou spectropolarimeter and precision velocimeter at the Canada-France-Hawaii Telescope from early 2020 to late 2025. Circularly polarized Zeeman signatures were clearly detected at most epochs in the atomic spectral lines of DO Tau, yielding longitudinal magnetic fields of up to 280 G modulated with a period of 5.128+-0.002 d which we identified as the rotation period of DO Tau. Applying Zeeman-Doppler imaging to the SPIRou data recorded in 2021, 2024 and 2025, we found that DO Tau hosts an unusual large-scale magnetic field that is weaker, less poloidal, more inclined to the rotation axis, and varies more rapidly with time than those of previously studied T Tauri stars, possibly as a result of intense accretion between the inner disk and the stellar surface. The dipole component of this large-scale field of about 0.2-0.3 kG even flipped polarity toward the end of our observing campaign, making DO Tau the first T Tauri star for which a magnetic polarity reversal is reported. The magnetospheric gap surrounding the central star was quite compact, extending to ~1.6 Rstar (0.014 au) as a result of the strong accretion rate (log Mdot = -7.7 Msun/yr), with the inner accretion disk being warped by the tilted stellar magnetic field. Radial velocity variations suggest the presence of a close-in planet of a few Mjup or a density structure in the inner accretion disk at an orbital period of 21 d (corresponding to 0.12 au), which might be linked to the wiggle in the jet axis of DO Tau.
Recent high-cadence observations by Subaru-HSC have identified a population of ultrashort-timescale microlensing events, providing a compelling window for planet-mass primordial black holes (PBHs) to constitute the entirety of dark matter. In this Letter, we demonstrate that this PBH population and the nanohertz stochastic gravitational-wave (GW) background reported by pulsar timing arrays (PTAs) can be naturally unified by a single primordial origin: a broad, nearly-flat enhancement of the curvature power spectrum with an amplitude of $O(10^{-2})$. The resulting PBH mass function spans the planet-to-solar mass range, while remaining consistent with all current observational constraints. This unified PBH--induced-GW framework makes concrete multi-messenger predictions, which can be decisively scrutinized by forthcoming microlensing surveys, next-generation PTAs, space-borne interferometers, precision astrometry, and laser ranging experiments.
Wide binaries (WBs) offer a unique opportunity to test gravity in the low-acceleration regime, where modifications such as Milgromian dynamics (MOND) predict measurable deviations from Newtonian gravity. We construct a rigorous framework for conducting the wide binary test (WBT), emphasizing high quality sample selection, filtering of poor astrometric solutions, contamination mitigation, and uncertainty propagation. We show that undetected close binaries, chance alignments, and improper treatment of projection effects can mimic MOND-like signals. We introduce a checklist of best practices to identify and avoid these pitfalls. Applying this framework to Gaia DR3 data, we compile a high-purity sample of WBs within 130 pc with projected separations of 1 - 30 kAU, spanning the transition between the Newtonian and MOND regimes. We find that the scaled relative velocity distribution of wide binaries does not exhibit the 20% enhancement expected from MOND and is consistent with Newtonian gravity across all separations. A meta-analysis of previous WBTs shows that apparent MOND signals diminish as methodological rigour improves. We conclude that when stringent quality controls are applied, there is no observational evidence for MOND-induced velocity boosts in wide binaries. Our results place strong empirical constraints on modified gravity theories operating between a0/10 and 200 a0, where a0 is the MOND acceleration scale. Across this range of internal accelerations, Newtonian gravity is up to 1500x more likely than MOND for our cleanest sample.
This work investigates symbolic regression (SR) as an interpretable alternative to black-box machine learning for the classification of stars, galaxies, and quasars in the Sloan Digital Sky Survey Data Release 17 (SDSS DR17). We conduct a systematic comparative study of four state-of-the-art SR frameworks: {\tt PySR}, Exhaustive Symbolic Regression ({\tt ESR}) with MDL-based selection, Physical Symbolic Optimization ({\tt PhySO}) using deep reinforcement learning, and Multi-View Symbolic Regression ({\tt MvSR}). By deriving compact analytic functions (complexity $\leq 10$) on a representative training subset and subsequently evaluating them via a 5-fold stratified cross-validation protocol on 100,000 spectroscopically confirmed objects, we map spectroscopic redshift ($z$) to continuous classification scores. Our results demonstrate that these low-complexity expressions achieve high predictive reliability, with {\tt MvSR} reaching a Cohen's Kappa of 0.8948 and {\tt PhySO} achieving exceptional parametric stability ($σ< 0.002$). We show that these models not only match the performance of traditional baselines but also provide a transparent, mathematically concise characterization of the astrophysical boundaries separating galactic and extragalactic populations.
The ionization mechanisms of low-ionization nuclear emission-line regions (LINERs), which are common in the local Universe, have been debated for decades. Our nearest large neighbor, M31, is classified as a LINER based on its optical emission line properties within the central kpc. In this work, we present a detailed photoionization modeling of the circumnuclear ionized gas in M31, explicitly tailored to its well-constrained physical conditions, including the absence of ongoing star formation and a currently inactive active galactic nucleus (AGN). Using spatially resolved CFHT/SITELLE observations, we find that photoionization by hot, evolved low-mass stars distributed throughout the bulge can roughly reproduce the observed radial intensity profiles of Hα, H\b{eta}, and [NII]. However, these models fail to match the observed [OIII] emission, producing radial profiles and [O III]/H\b{eta} ratios that are significantly steeper than observed. This discrepancy indicates a deficit of high-energy ionizing photons in standard stellar photoionization models, even with extended ionizing sources. We explore whether this tension can be alleviated by invoking either a bulge-filling, low-density ionized medium surrounding a denser Hα-emitting disk, or enhanced AGN activity in the recent past. While both scenarios can partially increase the [O III] emission, neither provides a fully satisfactory explanation under physically plausible conditions. Together with our earlier results for M81, these findings underscore persistent challenges in explaining LINER-like emission solely through conventional photoionization mechanisms.
We use the galaxy ultraviolet luminosity function measurements at $z=6-9$ to constrain modification to standard inflationary power spectrum. These observables are sensitive to the matter power spectrum which itself depends on inflationary initial conditions. We consider specific models where a bump or oscillatory features are introduced to the standard power law inflation spectrum. We find that the galaxy luminosity observables can probe such modifications at wavenumbers $0.5\lesssim k \lesssim 20$ Mpc$^{-1}$. We obtain upper limits on the amplitude of bump-like features at the mentioned wavenumbers. We obtain constraints which are similar to previous constraints on these models using measurements of optical depth of reionization. However, the galaxy luminosity functions are a more direct probe for these type of models and, therefore, can complement indirect constraints coming from measurements of IGM properties.
The Milky Way's nuclear star cluster (NSC) is a unique laboratory to study the formation and evolution of dense stellar systems around a supermassive black hole. Previous work suggests that most stars in the NSC are old; however, the detailed age and metallicity distributions remain uncertain. We constrain the star formation history (SFH) and metallicity of a poorly explored region located $\sim$3 pc from SagittariusA*. We analyse VLT/NACO imaging in an intermediate-band filter centred at 2.24 $μ$m, complemented by $H$-band data. We construct completeness-corrected $K$-band luminosity functions (LFs), clearly identifying the Red Clump and Red Giant Branch Bumps. The SFH is derived by fitting cumulative LFs with MIST, PARSEC, and BaSTI models spanning a wide range of ages and metallicities, using Monte Carlo sampling to estimate uncertainties. Metallicity constraints are refined using spectroscopic measurements from the literature. The stellar population is predominantly old and metal-rich: $75.6 \pm 9.5$% of the stellar mass formed $\gtrsim 10$ Gyr ago, with median [M/H] $\sim +0.35$. An intermediate-age component at 2-3 Gyr contributes $20.8 \pm 8.7$%, while minor populations are present at $\sim$400 Myr ($0.9 \pm 0.8$%) and 20 Myr ($3.6 \pm 1.4$%), the latter representing a small but non-negligible young population. Systematic uncertainties from stellar models, binning, photometric range, unresolved binaries, and filter choice are assessed. These results indicate early dominant formation, a significant 2-3 Gyr episode, and minor recent activity, consistent with spectroscopic measurements and with properties of the inner NSC and nuclear stellar disc.