Peer Revenue Publish or profit

Papers

Fresh from arXiv, typeset properly. Highlighted words trade on the market.

Physics Astrophysics

astro-ph Watch Astrophysics (sign in to use)

Tuesday 14 July

Earth and Planetary Astrophysics

Plasma chemistry and electron-moderated pathways in substellar atmospheres: a new perspective on the L/T transition

The long-standing puzzle of the CO/CH$_{4}$ transition in brown dwarfs endures. Although the bulk spectral evolution across an atmosphere can be accounted for through thermal equilibrium cloud models, the behaviour in the NIR remains unaccounted for, indicating that additional, non-thermal processes may influence atmospheric chemistry alongside conventional pathways. We explore cloud-driven electrical activation, where low-energy sparks to full lightning discharges, unlocks non-equilibrium reaction pathways inaccessible under thermal conditions alone. To quantify this, the aim of this paper is to model the electron-moderated atmospheric chemistry with SPARCKS, a bespoke zero-dimensional code that solves the coupled set of particle balance equations for substellar plasma activation and reaction kinetics, focusing on the key CO-CH$_{4}$ electron-moderated chemistry across the parameter ranges $T_{\rm gas} \in [700, 1600]$ K} and $T_{e} \in [2, 5]$ eV. We simulate a 1 microsecond pulse, representing a short dart-stepped leader; and, two pulsed systems with $(t_{\rm on}, t_{\rm off})$ = ($10^{-8}$ s, $10^{-6}$ s) and $(10^{-9}$ s, $10^{-9}$ s), representing small-scale inter-grain discharges, consistent with a typical characteristic substellar atmosphere. Our results show that even modest, physically plausible energies can strongly perturb atmospheric composition: an electron energy of 3.0 eV is sufficient to halve the CH$_{4}$/CO ratio in our sample atmosphere within one microsecond. Beyond the CO-CH$_{4}$ system, electron-moderated plasma chemistry exerts a far-reaching influence on substellar atmospheric composition.
Read the full abstract

High Energy Astrophysical Phenomena

Identifying potentially missed extended sources in the Fermi-LAT 4FGL Catalog using clustering analysis

Many Fermi Large Area Telescope (LAT) sources remain unassociated. Some may be extended sources represented by multiple catalog entries, as found for HESS J1813-178. We search for clusters of unassociated Fermi-LAT sources and test whether single extended-source models describe them better than multiple catalog sources. We applied DBSCAN to 4FGL sources with a linking scale of $0.3^\circ$ over 5 GeV-1 TeV. Each cluster contains at least one unassociated source and up to one extended or point-like associated source from selected categories, including pulsars, pulsar wind nebulae and supernova remnants. Using Fermipy, we compared extended- and multiple-source models and characterized each candidate spectrally and morphologically, focusing on the Galactic plane. We identified 48 clusters containing 124 sources, each with at least one unassociated source. For all 8 clusters passing our quality selection, an extended-source model is statistically preferred. At linking scales of $0.4^\circ$ and $0.5^\circ$, all 8 clusters retain their core sources. Cross-matches with the Second Fermi Galactic Extended Sources Catalog (2FGES) and the HESS Galactic Plane Survey (HGPS) showed that 5 clusters overlap known extended sources. The other 3 clusters, including one cluster with a morphology dependent on the interstellar emission model, have no counterpart in either catalog and are new extended-source candidates. Spatial clustering combined with likelihood-based model comparison can uncover extended sources missed in the Fermi-LAT catalog and complements existing searches.
Read the full abstract

All authors Giovanni Cozzolongo, Alison M. W. Mitchell, Samuel T. Spencer, Dmitry V. Malyshev, Tim Unbehaun

High Energy Astrophysical Phenomena · Cross-listed from gr-qc

Evaluating KAGRA upgrade scenarios for multimessenger observations of binary neutron stars

Binary neutron star mergers are key targets for multimessenger astronomy, motivating future upgrades of gravitational-wave detectors. For KAGRA, both broadband sensitivity improvements that increase the binary neutron star detection range, and high-frequency optimizations targeting neutron-star physics are under consideration. We present a computationally efficient framework to evaluate the multimessenger performance of detector upgrades by combining Fisher-matrix estimates of localization area and localization volume with detector duty factors and binary neutron star merger rates. We apply this framework to proposed KAGRA upgrade scenarios within the LIGO-Virgo-KAGRA network. For identical sources, the high-frequency upgrade improves sky localization by about 20% compared with the broadband option. However, when detection rates are taken into account, the broadband upgrade yields a larger number of well-localized events. Despite its shorter binary neutron star range than the other detectors, the inclusion of KAGRA increases the number of events localized within $10^3$ Mpc$^3$ volume by about 60%. These results provide a quantitative framework for evaluating future detector upgrades from the perspective of multimessenger observations.
Read the full abstract

All authors Yuta Michimura, Soichiro Morisaki, Kenta Hotokezaka, Masaomi Tanaka

Instrumentation and Methods for Astrophysics

A genetic algorithm for peer-review panel composition

The composition of scientific review panels is a constrained optimization problem in which a finite pool of experts must be distributed among multiple panels while balancing scientific expertise and demographic diversity. As the number of possible panel configurations grows very rapidly with the number of reviewers, exhaustive searches rapidly become computationally impractical. In this paper I present a genetic algorithm designed to optimize panel composition for the European Southern Observatory (ESO) proposal evaluation process. Panel assignments are represented through a chromosome-based encoding. Candidate solutions are evaluated using a fitness function based on four imbalance indicators: scientific expertise, gender, country affiliation, and professional seniority. The method is tested using real reviewer data from the ESO proposal handling system. The results show that the genetic algorithm rapidly identifies panel configurations with substantially lower imbalance than those obtained from random assignments and progressively improves the quality of the overall population. Beyond producing a single optimized configuration, the approach generates a set of high-quality panel realizations that can subsequently be filtered according to additional operational constraints not explicitly included in the fitness function. Although developed for ESO, the methodology is general and applicable to a wide range of panel-based peer-review systems.
Read the full abstract

Cosmology and Nongalactic Astrophysics

Exploring Hu-Sawicki-like modified gravity with Genetic Algorithms

We investigate whether viable Hu-Sawicki-like $f(R)$ models can produce deviations from $\Lambda\mathrm{CDM}$ that can be tested against current background cosmological data. We adopt a machine-learning approach based on Genetic Algorithms (GA) to reconstruct analytical perturbations around the Hu-Sawicki class of models. We develop a pipeline that interfaces the GATO GA code with the CANDI cosmology code. Each $f(R)$ function generated by the GA is first tested against theoretical viability conditions, including stability, the recovery of a standard matter-dominated epoch, the General Relativity limit, and chameleon screening mechanism. Viable candidates are then passed to CANDI to reconstruct the corresponding background cosmology and are tested against DESI DR2 BAO measurements and the Pantheon+ Type Ia supernova catalogue. %\newline The deviations we find are largest at late times, where the lower curvature makes modified-gravity effects more relevant, and are rapidly suppressed at higher redshift, in agreement with the imposed matching to the matter-dominated era. To further quantify deviations from the standard cosmological model, we compute the $Om(z)$ diagnostic. It shows only a very small departure from the constant $\Lambda\mathrm{CDM}$ behaviour. The effective dark energy equation of state associated with the reconstructed $f(R)$ function also evolves only weakly, showing a mild transition from an effective quintessence-like nature to an effective phantom-like regime. Overall, our results indicate that, within perturbations around the Hu-Sawicki class of models, current background data allow only limited deviations from $\Lambda\mathrm{CDM}$.
Read the full abstract

All authors Chiara De Leo, Elisa Fazzari, Matteo Martinelli, Savvas Nesseris

High Energy Astrophysical Phenomena

An internal shock model calibrated with real gamma-ray burst light curves using a genetic algorithm

The origin of gamma-ray burst (GRB) prompt emission remains an open question. The internal shock (IS) model is a leading scenario for converting relativistic ejecta kinetic energy into gamma rays, but its parameters have not yet been fully calibrated against observed GRB light curves (LCs) to reproduce their diversity. We adopt a machine-learning framework to optimise the IS model by comparing simulated and observed LC properties from three GRB catalogues (Swift/BAT, Fermi/GBM, CGRO/BATSE). Assuming a redshift-dependent GRB formation rate, we employ a genetic algorithm to minimise a loss function based on six independent metrics capturing both average behaviours and statistical distributions. The optimised model reproduces several key observational properties, including the average post-peak temporal profile, autocorrelation function, and the distributions of duration, signal-to-noise ratio, number of peaks, peak flux, and fluence. We also derive constraints on the central engine activity: (i) the number of emitted shells is well described by a generalised Zipf distribution, analogous to the Gutenberg-Richter law for earthquakes, and (ii) the rest-frame shell-emission times follow a negative exponential distribution, indicating a stochastic process with a constant ejection probability. This calibrated IS model provides a physically grounded framework for interpreting GRB variability and predicting GRB populations detectable by future missions.
Read the full abstract

All authors Manuele Maistrello, Cristiano Guidorzi, Shiho Kobayashi, Romain Maccary

High Energy Astrophysical Phenomena

Constraining the High-Density Equation of State with Present and Future NICER Observations Using Physics-Informed Regularized Machine Learning

The precise mass and radius measurements of neutron stars by NICER have significantly advanced our ability to constrain the properties of matter at supranuclear densities. In this work, we develop a physics-informed regularized conditional Invertible Neural Network (cINN) that bijectively maps mass–radius posterior distributions directly onto the corresponding central energy density and pressure, eliminating the need for explicit high-dimensional parameter sampling. The physics-informed regularisation guarantees that all inferred solutions satisfy causality and thermodynamic stability, ensuring physically consistent predictions without explicit forward modelling. We demonstrate that the framework accurately reconstructs central EoS posteriors for NICER-like observations while preserving the mapping between macroscopic stellar observables and the microscopic properties of dense matter. Exploiting the computational efficiency of the cINN, we perform a systematic optimisation study of 62,400 simulated mass–radius observations to identify the most informative targets for constraining the high-density EoS. We find that the constraining power depends strongly on the location of the observation in the mass–radius plane, with an optimal strategy that alternates between compact high-mass stars and extended intermediate-mass stars, reducing the uncertainty in the inferred EoS by up to $\sim 9\%-10\%$ relative to the current NICER baseline. These results establish physics-informed invertible neural networks as a powerful framework for rapid, physically consistent inference of dense-matter properties from present and future multi-messenger observations.
Read the full abstract

Astrophysics of Galaxies

Perturber-Driven Dynamics of Supermassive Black Hole Binaries in Galaxy Merger

The orbital eccentricity of massive black hole binaries (MBHBs) at binary formation shapes the stochastic gravitational-wave background (GWB) detectable by pulsar timing arrays (PTAs). Previous $N$-body simulations show large run-to-run scatter in this quantity, dominated by Poisson noise, raising the question of whether physical substructure adds genuine astrophysical stochasticity. We test this with high-resolution re-simulations of a major merger from IllustrisTNG100-1, evolved with the Griffin $N$-body code. A no-perturber control is compared with two matched suites in which $f_{\mathrm{target}}=0.1$ of the primary bulge mass is redistributed into equal-mass perturbers of $10^7,M_\odot$ ($\mu_{\mathrm{p}}\approx3.2\times10^{-3}$) and $10^8,M_\odot$ ($\mu_{\mathrm{p}}\approx3.2\times10^{-2}$), with four realisations per scenario. The control gives $\sigma_e\approx0.11$, consistent with the Poisson noise floor at this resolution. The $10^7,M_\odot$ case gives $\sigma_e\approx0.115$, indistinguishable from the control, whereas the $10^8,M_\odot$ case gives $\sigma_e\approx0.26$, a factor of $2.4$ above the floor, although statistically marginal given only four realisations. This excess scatter coincides with larger event-aligned residuals in orbital energy and angular momentum and stronger torque spikes, consistent with near-impulsive perturber–MBHB encounters. In binary–single scattering theory, the transition is set by the perturber–MBHB mass ratio $\mu_{\mathrm{p}}$: the $10^7,M_\odot$ case remains diffusive, whereas the $10^8,M_\odot$ case approaches the near-impulsive regime. Because the expected perturber population in massive ellipticals lies mostly below this regime, perturber-driven eccentricity randomisation is unlikely to affect GWB-relevant MBHB mergers.
Read the full abstract

All authors Julian Chan, Alessia Gualandris, Walter Dehnen, Justin I. Read

Astrophysics of Galaxies

Radar-Chart Analysis of Star-Formation Quenching Stages Across Circular Velocity Curve Classes in Nearby CALIFA Galaxies

We analyze the circular velocity curves of 215 non-active galaxies from the Calar Alto Legacy Integral Field spectroscopy Area survey. Using radar-chart analysis, we compare two classification schemes: the circular velocity curve classification and the quenching classification. We find a systematic progression from slow-rising to flat, round-peaked, and sharp-peaked curves that mirrors the transition from star-forming to fully retired systems. Star-forming galaxies are mainly slow-rising and flat, whereas fully retired galaxies are predominantly round-peaked and sharp-peaked, with intermediate stages showing a gradual shift. Overall, advanced quenching stages are associated with more centrally peaked curves and higher outer velocities (i.e., velocities along the flat part of the curve), indicating increasingly concentrated central mass distributions and more massive disks of the galaxies.
Read the full abstract

High Energy Astrophysical Phenomena

Coupled Shock Cooling and Radioactive Heating in the Type IIb Supernova SN 2024aecx: An Extended Envelope and Rapid Optical Decline

SN~2024aecx is a nearby, rapidly evolving stripped-envelope supernova with a prominent double-peaked ultraviolet–optical light curve. We model its multiband evolution with an extended version of TransFit, in which the early shock-cooling emission and the subsequent radioactive heating are treated within a single time-dependent radiative diffusion calculation. To describe the stratified ejecta expected for a Type~IIb progenitor, we adopt a compact inner ejecta connected to a dilute extended outer envelope and fit the outer density slope directly from the early light curve. The model reproduces the short-lived first peak, the rise to the radioactive main peak, and the overall multiband evolution. We infer an effective outer radius of $R_0=109.6^{+6.6}_{-3.5}\,R_\odot$, an ejecta mass of $M_{\rm ej}=2.14^{+0.21}_{-0.19}\,M_\odot$, a nickel mass of $M_{\rm Ni}=0.050\pm0.002\,M_\odot$, and a steep outer density slope of $n_{\rm out}=13.33^{+0.11}_{-0.12}$. The steep outer profile favors a low-mass extended envelope, while the low ejecta mass explains the rapid evolution of the main peak. However, a control model with standard $\gamma$-ray leakage fades too slowly after maximum. We therefore introduce an effective optical-output factor to quantify the additional late-time suppression of the ultraviolet–optical luminosity. These results support the shock-cooling plus radioactive-heating interpretation of SN~2024aecx, but show that its rapid optical decline requires physics beyond the simplest radioactive-diffusion prescription.
Read the full abstract

All authors Na Wei, Yu-Hao Zhang, Liang-Duan Liu, Guang-Lei Wu, Jing-Yao Li, Yun-Wei Yu, Ning-Chen Sun

Solar and Stellar Astrophysics

Symbiotic binaries in the Gaia data. II. Symbiotic candidates from general variability classification in DR3

We investigate the reliability of the symbiotic star class in the general variability classification of Gaia DR3 and search for new genuine symbiotic systems among the objects assigned to this category. The sample contains 649 sources, including 246 previously known symbiotic stars, 61 literature candidates, and 339 new candidates proposed by the Gaia variability pipeline. Diagnostics based on the Gaia colour-magnitude diagram, near-infrared photometry, and the pseudo-equivalent width of H$\alpha$ indicate that a large fraction of the new candidates are likely contaminants, predominantly pulsating red giants. To quantify the contamination, we constructed a Random Forest classifier trained on confirmed symbiotic stars and on Mira and semi-regular variables, using Gaia photometry, variability parameters, H$\alpha$ measurements from XP spectra, and infrared colours. The classifier reaches a balanced accuracy of $\approx$0.94 and efficiently separates most symbiotic binaries from single evolved stars, leaving only eight strong candidates among the 339 newly proposed objects. Follow-up spectroscopy confirms three new symbiotic stars through the presence of high-excitation emission lines, while several additional objects remain possible symbiotics. Our results show that the Gaia DR3 variability classification efficiently recovers known symbiotic stars but has low purity due to overlap with pulsating red giants. The small number of newly confirmed systems implies that the discrepancy between predicted and observed Galactic symbiotic populations remains unresolved, although Gaia provides a powerful basis for future searches combining variability, spectroscopic indicators, and multi-wavelength data.
Read the full abstract

All authors Jaroslav Merc, Lionel Mulato, Joanna Mikołajewska, Paul G. Beck, Ana Escorza, Michael Abdul-Masih, Stéphane Charbonnel, Olivier Garde, Pascal Le Dû, Thomas Petit

Instrumentation and Methods for Astrophysics

Status and future development of the COSMOCal Project for absolute CMB polarization calibration

Cosmic Microwave Background (CMB) polarization measurements are pushing instrumental sensitivities to levels where calibration systematics become a dominant limitation. The large dynamic range between cosmological and Galactic emission prevents future experiments from relying only on diffuse sky measurements or standard celestial calibrators. To address this challenge, the COSmological Microwave Observations Calibrator (COSMOCal) project proposes an artificial calibration source deployed as a guest payload on a geostationary satellite, scheduled for launch by the Eutelsat group by 2030. This source will provide stable, well-characterized polarized microwave signals accessible to multiple ground-based observatories. In this work, we present the status of the project, the updated development timeline, and the refined scientific and technical requirements, defined with the observatories that plan to use this calibration source. Furthermore, we investigate the interplay between instrumental systematics and component separation in the presence of complex models of interstellar dust emission. We discuss in this paper how this can impact the recovery of the primordial signal, and whether residual calibration errors can degrade the performance of foreground cleaning algorithms.
Read the full abstract

All authors Silvia Micheli, Alessia Ritacco, Alessandro Carones, Jonathan Aumont, Stefano Berta, Ludovico Bizzarri, François Boulanger, Andrea Catalano, Francesco Cuttaia, Anne Denis, François-Xavier Désert, Anne Laure Fontana, David González-Ovejero, Benedetta Kalemi, Nicoletta Krachmalnicoff, Samuel Leclercq, Thibaut Louis, Juan-Francisco Macías-Pérez, Bruno Maffei, Tania Mcnamara, Marina Migliaccio, Ludovic Montier, Pascal Morfin, Michel Moulin, Louise Mousset, Matteo Murgia, Ioannis Myserlis, Federico Nati, Pierluigi Ortu, Michel Pérault, Giampaolo Pisano, Tonino Pisanu, Nicolas Ponthieu, Sofia Savorgnano, Luca Terenzi, Jeanne Treuttel, Léo Vacher, Christophe Vescovi, Mario Zannoni

Astrophysics of Galaxies

Tidal Grinding of Dwarf Galaxies in Cluster Environments

Dwarf elliptical galaxies (dEs) dominate galaxy clusters and provide key constraints on environmentally driven galaxy evolution. Here we examine whether the projected shapes of dEs retain information about their accretion and transformation histories using a homogeneous sample of 1,108 bright (m_g < 19 mag) dEs in the Virgo cluster. Based on the axis-ratio (b/a), we define flat (< 0.70) and round (> 0.74) subsamples and compare their spatial and kinematic properties. We find that flat dEs are distributed more uniformly across the cluster, whereas round dEs preferentially occupy regions of stronger tidal fields around massive (M_* > 10^{10} M_sun) galaxies. Within the central 5^\circ x 5^\circ region around the Virgo central galaxy (M87), 149 dEs have spectroscopic radial velocities compiled from public archives. In this region, the two shape classes also exhibit clear kinematic segregation: flat dEs have systematically larger line-of-sight velocity offsets from the cluster mean (median $\Delta v = 654 km/s$), whereas round dEs have smaller offsets (median $\Delta v = 414 km/s$), as expected for a more dynamically relaxed population. Flat dEs are consistent with a population that has experienced weaker tidal processing and consequently retains more flattened morphologies. By contrast, round dEs are consistent with prolonged tidal processing (“tidal grinding”) that may have transformed initially flattened systems into rounder spheroids. However, projection contamination of the round subsample may have introduced some uncertainty in the interpretation of intrinsic galaxy shapes.
Read the full abstract

All authors Sanjaya Paudel, Suk-Jin Yoon, Tek Prasad Adhikari, Eun-Taek Gim, Myung-Hun Kim, Inhyuk Park, Nau Raj Pokhrel

Instrumentation and Methods for Astrophysics

ngVLTI: Cost and Feasibility of a Four-Telescope VLTI Expansion for Milliarcsecond-Scale Imaging

We assess the cost and technical feasibility of extending the Very Large Telescope Interferometer (VLTI) through the addition of four 8 m-class Unit Telescopes (UTs), an upgrade we refer to as the new generation VLTI (ngVLTI). Such an upgrade would provide a dense and homogeneous uv coverage with baselines up to 220 m, enabling true imaging at milliarcsecond angular resolution across science cases ranging from Solar System bodies to distant active galactic nuclei. Because each image is reconstructed at a single wavelength, repeating the reconstruction across the spectral channels of the instrument would deliver spectral-imaging cubes – a qualitatively new capability for the VLTI, bringing it close to the imaging power of ALMA at near-infrared wavelengths. Motivated by this scientific case, we examine a compact telescope concept consisting of a fast ($f/0.64$), segmented, parabolic primary mirror feeding a subterranean coudé focus compatible with the existing VLTI infrastructure. We summarise the optical design, which achieves diffraction-limited performance over a 1 arcmin field of view and a well-matched reimaged pupil, discuss the mechanical trade-offs behind the choice of a 60-segment, 1.2 m primary – a mirror mass of ~13.5 tons and an altitude moving mass of ~50 tons – and quantify the gravitational flexure of the telescope structure and its resulting optical sensitivity as a function of pointing elevation. We then present a back-of-the-envelope cost estimate of order 80 MEuro per telescope (2026 prices), broken down into the segmented primary, adaptive secondary, coudé train, mount, enclosure, and ancillary instrumentation. The proposed upgrade appears both cost-competitive and technically achievable, offering a long-term perspective for Paranal Observatory in the ELT era.
Read the full abstract

All authors Sylvestre Lacour, Johan Kosmalski, Michael Mueller, Célia Desgrange, Oscar Carrión-González, Mathias Nowak, Elsa Huby, Julien Woillez, Jens Kammerer, Jean-Baptiste Le Bouquin, Jason Spyromilio

Solar and Stellar Astrophysics

Symbiotic binaries in the Gaia data. I. Known symbiotics in DR3 and FPR

Symbiotic stars are long-period interacting binaries composed of an evolved giant and a hot compact companion. Their complex spectra and variability make them both astrophysically valuable and observationally challenging. We investigate how known symbiotic stars are represented in Gaia DR3 and the Focused Product Release (FPR), assess the reliability of these data, and evaluate their usefulness for candidate analysis and searches for new systems. We crossmatched Gaia DR3 and the FPR with confirmed symbiotic stars from the New Online Database of Symbiotic Variables and examined their astrometric, photometric, and spectroscopic data, along with derived products such as astrophysical parameters, orbital solutions, variability properties, and emission-line classifications. Astrometric data reliably constrain the position of symbiotics in the color-magnitude diagram, aiding searches for new systems, while RUWE is generally not a reliable indicator of their binarity. Most symbiotics are variable in Gaia photometry. Mean and epoch radial velocities, as well as inferred orbital solutions, are broadly consistent with the literature, and we provide the first tentative orbital solutions for two systems. Effective temperatures and metallicities are unreliable due to contamination from nebular continuum and strong emission lines. H$\alpha$ emission is detected in nearly all symbiotics, making it a robust diagnostic. Additionally, resolved companions were identified for two systems, and one previously confirmed symbiotic star was reclassified as non-symbiotic. Gaia DR3 provides a rich dataset for the study of symbiotic binaries. The forthcoming DR4 promises a major leap forward with its longer time baseline, new data products, and epoch data.
Read the full abstract

High Energy Astrophysical Phenomena · Cross-listed from nucl-th

Thermal conductivity of dense $np\Lambda$ matter in neutron star cores

The possible presence of hyperons in the cores of massive neutron stars has important implications for their microscopic transport properties and thermal evolution. Despite recent progress in modelling core transport, a dedicated analysis regarding the thermal conductivity of a specific \(np\Lambda\) mixture remains absent from the existing literature. In this work, we investigate the thermal conductivity of dense, \(\beta \)-equilibrated \(np\Lambda\) matter within the framework of the variational linearized Boltzmann kinetic approach, employing the density-dependent DDME2 equation of state across a baryon density range of \((0.5\text{--}4.5)\,n_0\). We find that neutrons still dominate thermal transport, while the onset of \(\Lambda \) hyperons induces only a remarkably small reduction in neutron conductivity compared to pure nucleonic matter. These results suggests that the core thermal relaxation timescale remains practically unaltered in the presence of $\Lambda$ hyperons.
Read the full abstract

Astrophysics of Galaxies

Discovery of a star sensitive to the spin of Sgr A*

Residing in the center of the Milky Way, Sgr A* is the closest massive black hole (MBH). Its vicinity has allowed measuring individual stellar orbits around it. The stars act as test particles and probe the gravitational potential around the $4.3 \times 10^6 M_\odot$ MBH. These observations have determined the central mass to sub-percent precision, and the mildly relativistic motions of stars have given access to the dominant relativistic corrections, the gravitational redshift, the transverse Doppler effect, and the prograde precession imposed by the Schwarzschild metric nature of the potential. These effects are of order $\beta^2 = (v/c)^2$ (for velocity $v$ and speed of light $c$). The Kerr metric for a rotating black hole leads to corrections of order $\beta^3$. Here, we report the discovery of a faint main-sequence star ($m_K = 19.3$), S301, on a 8.7-year orbit and with small enough a pericenter distance, such that the star's peak velocity reaches $25000\,$km/s. Within the measurement capabilities of current near-infrared interferometry and future spectroscopy on an extremely large telescope, S301's motion is directly sensitive to the spin of Sgr A*. The high eccentricity of S301 suggests that it is the captured component of a binary that was torn apart via the Hills mechanism.
Read the full abstract

All authors K. Abd El Dayem, R. Abuter, N. Aimar, P. Amaro-Seoane, A. Berdeu, J.-P. Berger, G. Bourdarot, W. Brandner, A. Burkert, D. Calderon, C. Correia, J. Cuadra, R. Davies, D. Defrere, L. Delit, A. Drescher, F. Eisenhauer, L. Esteras Otal, M. Fabricius, H. Feuchtgruber, N.M. Foerster Schreiber, A. Foschi, P. Garcia, R. Garcia Lopez, A. Generozov, R. Genzel, S. Gillessen, F. Gonte, X. Haubois, S.F. Hoenig, M. Houlle, S. Joharle, A. Kaufer, J. Kammerer, P. Kervella, J. Kolb, L. Kreidberg, L. Labadie, S. Lacour, O. Lai, R. Laugier, J.-B. Le Bouquin, J. Leftley, B. Lopez, D. Lutz, F. Mang, A. Merand, F. Millour, M. Montarges, N. Morujao, H. Nowacki, M. Nowak, S. Oberti, J. Osorno, T. Ott, T. Paumard, C. Paladini, H.B. Perets, K. Perraut, G. Perrin, R. Petrov, P.O. Petrucci, T. Piran, N. Pourre, S. Rabien, D.C. Ribeiro, S. Robbe-Dubois, M. Sadun Bordoni, J. Sanchez Bermudez, D. Santos, R. Sari, J. Sauter, S. Scheithauer, J. Scigliuto, J. Shangguan, T.T. Shimizu, F. Soulez, J. Stadler, C. Straubmeier, E. Sturm, M. Subroweit, C. Sykes, L.J. Tacconi, P. Thevenet, I. Urso, F. Vincent, J. Woillez, G. Zins

Earth and Planetary Astrophysics

Direct Retrieval of Protoplanetary Disk Dust Properties using Auto-differentiable Gaussian Processes and Its Application to the HD 169142 Disk

Retrieving dust properties in protoplanetary disks, including the surface density distribution, temperature, and grain size distribution, is a fundamental task in observational studies of planet formation. While multi-wavelength analysis of the spectral energy distribution (SED) using interferometers such as the Atacama Large Millimeter/submillimeter Array (ALMA) is a powerful diagnostic tool, traditional methods are often hindered by strong biases arising from a limited imaging beamsize. In this paper, we present a new retrieval framework for dust disk properties designed to overcome this challenge. We assume that the underlying physical structures are expressed as sample paths from Gaussian processes, compute the radial intensity distributions at observed wavelengths, and produce one-dimensional visibility models. The models are compared with the observed data and the posterior distributions are sampled via the Markov-Chain Monte-Carlo method. The whole procedure is implemented in JAX, which enables end-to-end auto-differentiation and significantly accelerates the inference. We validate our methodology using mock datasets, and find that the results are not strongly biased and better reproduce the input profiles. We also demonstrate its capabilities through an application to ALMA Band 3, 6, and 9 observations of the HD 169142 disk, revealing a new complex structure. Our developed code is publicly available as a Python module, FRAP (Flexible Radial Analysis of Protoplanetary disks). This framework provides a next-generation infrastructure for disk SED modeling, enabling high-precision studies of the physical environments in which planets form.
Read the full abstract

All authors Tomohiro C. Yoshida, Enrique Macías, Giovanni Rosotti, Stefano Facchini, Elena Viscardi, Kiyoaki Doi

Instrumentation and Methods for Astrophysics

PSRDISP: A novel approach to modeling dispersive processes in single-pulsar noise analysis using epoch-wise dispersion measures

We present PSRDISP, a novel approach to modeling deterministic and stochastic dispersive processes in pulsar timing datasets using high-precision epoch-wise dispersion measure (DM) estimates, with a Gaussian Process-based approach. Unlike the conventional single-pulsar noise analysis methodology, which is applied to frequency-resolved times of arrival (ToAs) of pulses, this technique is applied to epoch-wise DMs which are derived from these ToAs. It can also be applied to wideband DMs measured simultaneously with wideband ToAs. Therefore, this framework provides a paradigm-agnostic approach to characterise single-pulsar dispersive processes. This method is expected to minimise the impact of achromatic red noise processes while characterising these dispersive effects. We substantiate the discussed technique with representative examples using simulated narrowband and wideband datasets with realistic noise injections. We found the recovery to be in close agreement with the injections, and agnostic to the estimation technique. Our method applies to pulsar timing experiments where precise, epoch-wise DM estimates are possible, such as the Indian Pulsar Timing Array. This technique can serve as a powerful diagnostic tool for validating single-pulsar noise analyses, which is crucial for precision pulsar timing experiments, such as Pulsar Timing Arrays.
Read the full abstract

Solar and Stellar Astrophysics

Mean-field interpretation of star-in-a-box simulations of red giants

Context: The origin of magnetic fields and the dynamo mechanism in red giants are still not fully understood. Aims: We aim to interpret the dynamo behaviour of global 3D simulations of red giants using mean-field dynamo models. Methods: We use mean-field models constrained by the differential rotation profile extracted from 3D simulations. We perform $\alpha^2$, $\alpha\Omega$, and $\alpha^2\Omega$ mean-field dynamo simulations, varying the strength of the $\alpha$ effect and the differential rotation. Results: The mean-field models can reproduce the growth rate of several 3D runs. The morphology of the large-scale magnetic field is better reproduced for the slowly rotating 3D cases than for the rapidly rotating ones. Faster rotation enhances dynamo action and it also modifies the dominant mode of the dynamo. Rapidly rotating 3D runs produce predominantly non-axisymmetric equatorial dipoles instead of axisymmetric fields at slower rotation. Mean-field models are supercritical even in the absence of differential rotation, indicating $\alpha^2$ dynamo action. By contrast, models using the $\alpha\Omega$ approximation require sufficiently strong differential rotation to become supercritical. Conclusions: Our results suggest that the magnetic field in red giants requires dynamo action. In our mean-field runs, differential rotation speeds up the magnetic decay, disfavouring the idea of a persistent fossil magnetic field in red giants.
Read the full abstract

Astrophysics of Galaxies

The Three Hundred project: Low Gas Fraction Galaxy Clusters properties and their environment

Galaxy cluster samples based on X-ray and Sunyaev-Zel'dovich (SZ) observations are affected by selection biases. These catalogs preferentially include systems with high gas content and surface brightness. Excluding objects with depleted gas content, low-gas-fraction clusters (LGFCs), could lead to an incomplete sampling. We aim to investigate the abundance and the properties of the LGFCs population using The Three Hundred hydrodynamical simulations, focusing on the Gadget-X code. In particular, we study outliers in the $f_{\mathrm{g},500} - M_{500}$ relation, environmental influences, and their behavior in key scaling relations, with a focus on the Compton-Y observable. We analyze a sample of $N_{\mathrm{tot}} = 9858$ simulated objects from The Three Hundred, in the redshift band $z \in [0;0.817]$. LGFCs are selected statistically as outliers of the $f_{\mathrm{g},500}-M_{500}$ relation. To analyze environmental effects, we compare the gas density and temperature radial profiles of LGFCs against the No-LGFCs population. Finally, we study how the temperature, entropy, and spherical Compton parameter scaling relations are affected by the inclusion of LGFCs. We find that LGFCs are preferentially found at the low-mass end and their abundance increases toward low redshift. Radial profiles of LGFCs show lower gas concentrations in the core regions and higher temperatures, suggesting a more diffuse and heated ICM. This behavior is also reflected in the entropy scaling relation, where LGFCs are extreme positive outliers. Contrary to observations, the $Y_{\mathrm{sph},500}$ values of LGFCs show no significant deviation from the general population. Nevertheless, we cannot rule out that these differences are partly driven by the mass incompleteness at the low-mass end and the environmental bias of our simulated sample.
Read the full abstract

All authors Francesco Guidi, Marco De Petris, Sara Santoni, Raphael Wicker, Antonio Ferragamo, Hippolyte Froget, Stefano Andreon, Weiguang Cui, Elena Rasia, Gustavo Yepes

High Energy Astrophysical Phenomena

Machine Learning Detection of Non-Axisymmetric Fast Flavor Instabilities in Compact Objects

Neutrinos in dense astrophysical environments such as core-collapse supernovae (CCSNe) and neutron star mergers (NSMs) can undergo FFCs, which could develop on extremely small scales. A necessary condition for the occurrence of FFCs is the presence of a zero crossing in the electron lepton number (ELN) angular distribution of neutrinos. In this work, we explore machine learning (ML) approaches to detect non-axisymmetric ELN crossings in these environments, based on input features of the $\nu_e$ and $\bar\nu_e$ zeroth and first angular moments. Overall, the ML models demonstrate relatively good generalizability for most of the unseen test datasets generated by various methods that do not assume the same underlying angular distributions as used in the training set. Interestingly, while the model's performance is mediocre for an axisymmetric distribution dataset derived by solving the discretized Boltzmann transport equation under 1D CCSN background, imposing an artificial non-axisymmetry substantially improves the performance. We also find that for the flavor-equilibrated angular distributions, although our ML model trained based solely on ELN inputs performs poorly when the true crossings depend on the post-equilibrated angular distributions of heavy lepton neutrinos and antineutrinos, which become different, it delivers strong performance in detecting ELN crossings when the heavy-lepton neutrino and antineutrino distributions are artificially removed. This highlights the need for additional input features to further improve the model. This is a crucial step toward successfully integrating FFCs into large-scale CCSN and NSM simulations.
Read the full abstract

All authors Madhurima Chakraborty, Sajad Abbar, Meng-Ru Wu, Francois Foucart, Zewei Xiong

High Energy Astrophysical Phenomena

Photometric and Spectroscopic Studies of Type Ic Supernovae SN 2020akf and SN 2021mxx

We present a comprehensive analysis of optical photometry and medium-resolution spectroscopy for Type Ic supernovae (SNe) SN 2020akf and SN 2021mxx. Our study covers the evolution of SN 2020akf from -5 to 109 days and SN 2021mxx from -3 to 115 days relative to their B-band maximum. Their peak quasi-bolometric luminosities are estimated at logL_bol = 42.37 +- 0.02 and 42.21 +- 0.02 erg/s, respectively. The velocities of the ejecta, derived from the Fe II 5169 A line at maximum light, are approximately 15000 km/s for SN~2020akf and about 10000 km/s for SN~2021mxx, which are consistent with those observed in other Type~Ic SNe. Using the semianalytical Arnett model, we estimate that SN 2020akf has a kinetic energy of E_k = 6.33 (+0.68)(-0.62) x 10^51 erg and an ejected mass of M_ej = 4.71 (+0.50)(-0.46) M_sun, while for SN 2021mxx, we get E_k = 0.54 (+0.08)(-0.12) x 10^51 erg and M_ej = 0.90 (+0.14)(-0.18),M_sun. The mass of Ni56 synthesized in the explosion is estimated at 0.10 +- 0.02 M_sun for SN 2020akf and 0.05 +- 0.01 M_sun for SN 2021mxx. The metallicities of the host galaxies near the SN regions are ~0.81 Z_sun for SN 2020akf and ~0.76 Z_sun for SN 2021mxx, where Z_sun indicates solar metallicity. Our detailed analysis suggests that SN 2020akf falls into the category of a transitional Type Ic SN, having spectral properties between normal and broad-line Type Ic SNe, while SN 2021mxx is a normal Type Ic SN with an extremely low value of the M_ej / E_k ratio.
Read the full abstract

All authors Mulchand Kurre, K. R. Sahu, Shrutika Tiwari, Yogita Patel, N. K. Chakradhari, D. K. Sahu, Rishabh Singh Teja, G. C. Anupama, Mridweeka Singh, Anirban Dutta, Anjasha Gangopadhyay

Instrumentation and Methods for Astrophysics

Estimation of the laser guide star uplink tip-tilt using aperture size diversity

Laser guide star (LGS) adaptive optics cannot directly measure tip-tilt (TT), forcing reliance on natural guide stars and limiting sky coverage. We propose estimating the uplink TT from telemetry acquired at the laser launch telescope alone, operated in a monostatic configuration as both emitter and receiver. Extracting TT over concentric disks of different diameters within the receiving pupil yields signals mixing uplink and downlink contributions in different proportions; this aperture size diversity, combined with an LMMSE estimator, disentangles the uplink component. Simulations show a residual error of 24 mas for a single turbulent layer and 34 mas with two layers.
Read the full abstract

All authors Pierre Jouve, Thierry Fusco, Carlos M. Correia, Cédric Taissir Heritier, Jean-François Sauvage, Benoit Neichel

Earth and Planetary Astrophysics

Discrimination of Tectonic, Swarm and Impact-related Marsquakes using Spectral Characteristics

The NASA InSight mission observed over 2000 marsquakes in the course of its three year mission. These quakes varied in magnitude between 1.5 and 4.5, as well as in spectral content. We present a simple framework to describe the spectral characteristics of all observed marsquakes, based on source process; propagation through the mantle or crust; and local, receiver-side amplification. We assign to each quake an objective measure of its amplitude, as well as the spectral decay created by the duration of the rupture and the dampening of high frequencies due to visco-elastic attenuation. Together, this allows us to obtain characteristic patterns of the whole marsquake dataset, e.g. in terms of event magnitudes, source size, and - for quakes caused by meteoritic impacts - crater size. We show that a significant fraction of all marsquakes - the high-frequency quakes - form a swarm that is likely not caused by tectonic processes in rocks. Our analysis allows separation of the whole marsquake catalogue into three event classes, of tectonic quakes, meteoritic impacts, and swarm events. We finally conclude that the largest marsquake, S1222a, most likely belongs to the group of meteoritic impacts.
Read the full abstract

All authors Simon C. Stähler, Savas Ceylan, Nikolaj Dahmen, John Clinton, Domenico Giardini

Instrumentation and Methods for Astrophysics

Scaling laws for multi-object spectrographs: empirical relationships and the photonic advantage of CAWSMOS

Conventional multi-object spectrographs (MOS) incur masses and costs that scale unfavourably with channel count~$N_\mathrm{ch}$ and resolving power~$\mathcal{R}$. A dataset of eight fibre-fed MOS-VIS instruments (1995–2024) is compiled and analysed. Ordinary least-squares fitting yields a mass function $M \propto N_\mathrm{ch}^{0.78\pm0.28}$; the resolving-power exponent is empirically consistent with zero ($\gamma = 0.20\pm0.37$) but is predicted analytically to lie in $\mathcal{R}^{0.8\text{--}2.0}$. Leave-one-out cross-validation confirms sub-linear channel-count scaling ($\beta\in[0.59,1.13]$, mean 0.79), with HERMES and AAOmega – sharing the same telescope and $N_\mathrm{ch}=392$ across a factor of~7 in $\mathcal{R}$ at essentially equal estimated mass – providing a direct empirical confirmation that $\gamma\approx 0$. Cost per channel decreases as $N_\mathrm{ch}^{-0.22\text{ to }-0.53}$ for all assumed cost–mass exponents $k\in[0.6,1.0]$, demonstrating a genuine economy of scale. The photonic integrated-circuit (PIC) approach embodied by the arrayed waveguide grating (AWG)-based PAWS demonstrator and the proposed CAWSMOS instrument decouples dispersive element size from $N_\mathrm{ch}$ and $\mathcal{R}$, with chip-area scaling as $\mathcal{R}^{1\text{--}1.5}$ instead of $\mathcal{R}^{3}$. Applied to the Wide-field Spectroscopic Telescope (WST), the model predicts $\sim$140 tonnes for a conventional NIR spectrograph system against $\sim$200 kg for a CAWSMOS equivalent, with projected costs of EUR 300–700M conventional versus EUR 40–120M photonic ($6$–$15\times$ reduction). AWG-based photonics are identified as a strategially important enabling technology for WST.
Read the full abstract

Cosmology and Nongalactic Astrophysics · Cross-listed from hep-ph

Neutrino superradiance constraint on asteroid-mass PBH Dark Matter and beyond

Primordial Black Holes (PBHs) are an attractive candidate for Dark Matter (DM) and there has been extensive experimental efforts to look for them. The asteroid-mass window, $M_{\text{PBH}} \sim 10^{17} - 10^{23}$ g, is particularly interesting, since PBHs in this range may still constitute all of DM. In this work, we study a scenario in which rotating PBHs are surrounded by superradiantly produced boson clouds that emit an approximately steady and nearly monochromatic flux of neutrinos in the few MeV range. We compute both the Galactic and extragalactic neutrino fluxes from such PBH populations and compare them with existing low-energy antineutrino limits from Borexino, KamLAND, and Super-Kamiokande. For scalar bosons, we find that these neutrino searches can strongly constrain a significant part of the asteroid-mass window and extend to somewhat larger masses. For instance, for rapidly rotating PBHs with spin ã = 0.9 and gravitational fine-structure coupling $\alpha_g=0.25$, the strongest bound on dark matter fraction reaches approximately $f_{\text{PBH}} \sim 10^{-7}$ around $M_{\text{PBH}} \sim 2 \times 10^{22}$ g for a Yukawa coupling $g_{\nu \phi} = 10^{-4}$. These constraints can be significantly stronger than existing microlensing limits in the same mass range. Our results provide a complementary neutrino probe of PBH DM, distinct from previous neutrino constraints based mainly on Hawking evaporation from lighter PBHs.
Read the full abstract

Earth and Planetary Astrophysics

A Fast Concentric-disk Contour Integration Method For Microlensing Limb-darkening Effect

Incorporating limb darkening is a computationally demanding step in contour-integration-based microlensing modeling. Conventional concentric-ring integration is incompatible with high-order quadrature schemes, limiting efficiency. We develop a new concentric-disk method, which reformulates the limb-darkening integral and enables the application of high-order adaptive quadrature, significantly accelerating the computation. While mainly demonstrated using a linear limb-darkening profile, the method readily extends to more general limb-darkening profiles. As a source effect, it applies to lens systems of any complexity. The method achieves a convergence rate scaling faster than $N_{\rm{uni}}^{-4}$ with the number of uniform-source magnification evaluations $N_{\rm{uni}}$, a significant improvement over the concentric-ring $N_{\rm{uni}}^{-2}$ scaling. For a relative accuracy of $10^{-6}$, the concentric-disk approach typically requires only $30\%$ or less of the computational cost of traditional algorithms. This method has been implemented in the binary-lens contour integration code Twinkle, providing an efficient and precise tool for analyzing current and future high-precision microlensing observations.
Read the full abstract

Astrophysics of Galaxies

Accretion-disk formation around orbiting stellar black holes in gaseous star clusters

We consider low-mass black holes (BHs) moving in regular orbits in the cores of non-rotating gaseous star clusters, representative of proto-stellar clusters or the centers of protogalaxies. We argue that as the BH's sphere of influence – the Bondi sphere – is advected along the BH trajectory, the transverse velocity shear between the inner and outer hemispheres injects angular momentum, driving the formation of an accretion disk. Coriolis forces oppose angular momentum injection, delaying disk formation but not preventing it. The disk lies in the BH orbital plane and is counter-rotating with respect to the orbital BH motion. We verify this picture with 2D and 3D hydrodynamic simulations in the non-inertial frame of the orbiting BH. We find that the disk-formation timescale following a disruption event is of order the Bondi crossing timescale, $\tau_{\rm d} \sim R_{\rm B}/V_{\bullet}$, and that the disk radius is of order $R_{\rm d} \sim \omega_{\bullet}^2 R_{\rm B}^4/ G m_{\bullet} $, set by the circularization radius of gas captured in the Bondi sphere. For the case of a BH with $m_{\bullet} = 50\,{\rm M}_\odot$, inside the core of a typical compact proto-stellar cluster, these values read $\tau_{\rm d} \sim 0.1 P_{\rm orbit}$ and $R_{\rm d} \sim 10^{-3} R_{\rm B}$.
Read the full abstract

Cosmology and Nongalactic Astrophysics · Cross-listed from hep-ph

Measuring Ultralight-Axion Coherence with Galaxy Polarization Correlations

Ultralight axion-like particles coupled to photons rotate the linear polarization of distant sources. We propose using the three-dimensional two-point correlation of galaxy polarization-rotation angles to measure not only the amplitude of this birefringence field but also its spatial coherence scale. For a nonrelativistic ALP component with an isotropic Gaussian velocity distribution, the equal-time field correlation has an $e^{-1}$ scale $L_{\rm G}^{\rm phys}=\sqrt{6}/(m_a v_a)$, where $v_a$ is the three-dimensional rms ALP velocity dispersion. A detected turnover in the galaxy-pair correlation therefore measures the characteristic momentum scale $m_a v_a$, while the correlation amplitude constrains $g_{a\gamma}\sqrt{\Omega_a/\Omega_{\rm DM}}$. For a fiducial survey with $10^6$ polarized galaxies over a quarter of the sky to $z=2$ and effective per-galaxy scatter $\sigma_{\rm tot}=10^\circ$, we find $5\sigma$ sensitivity to sub-degree correlated rotations over $m_a(10^{-3}c/v_a)\sim10^{-29}$–$10^{-27}\,{\rm eV}$, with $\sigma[\log_{10}(m_a v_a)]\lesssim0.1$ for detected signals across much of this range. This provides a geometric late-time probe complementary to CMB birefringence and structure-formation constraints.
Read the full abstract

2607.12446 · 14 Jul · The author is withdrawing this manuscript because its key predictions, namely a 5{\sigma} detection and a 0.1 dex precision constraint on ultralight axion coherence, cannot be adequately substantiated under realistic observational conditions. These predictions rely on idealized assumptions, and the calculations do not provide sufficient support for the claimed observational sensitivity