Peer Revenue An exchange of ideas

Papers

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

Physics Astrophysics

astro-ph Watch Astrophysics (sign in to use)

Monday 2 February

Astrophysics of Galaxies

The Rhythm of the ISM: Tracing the Timescales of Gas Evolution and Star Formation across Galactic Environments

We investigate the physical origin of the star formation scaling relations between the gas depletion time, the star-forming gas mass fraction, and the gas surface density, $\Sigma_{\rm gas}$, on kiloparsec scales, all of which are the key ingredients of the observed Kennicutt-Schmidt relation. To elucidate these trends, we employ an analytical framework that explicitly connects these kiloparsec-scale properties to the timescales governing the rapid, continuous ISM gas cycle on the scales of individual star-forming regions, including the formation, dispersal, and local depletion of star-forming gas. Using a suite of idealized disk galaxy simulations spanning a range of environments from dwarf and Milky Way-mass systems to a gas-rich starburst analog, we measure the timescales of the gas cycle and relate them to the dynamical and turbulent properties of the interstellar medium (ISM). We find that star-forming regions form on a timescale close to the vertical turbulent crossing time of the galactic disk, $\sim$3-30 Myr, which decreases at higher $\Sigma_{\rm gas}$ due to the increase in turbulent velocities in the ISM and the decrease in the disk thickness. In contrast, the local star formation and dispersal of such gas are set by the local conditions. Specifically, the local depletion time, $\sim$200-2000 Myr, is decreasing at higher $\Sigma_{\rm gas}$, as star-forming gas becomes denser and more efficient in forming stars. The lifetime of such gas is very short, $\sim$0.4-1 Myr, and only weakly increases with $\Sigma_{\rm gas}$. Together, our results demonstrate how the star formation properties of galaxies on kiloparsec scales emerge directly from the interplay between the galaxy-scale dynamics, ISM turbulence, and the state of star-forming gas.
Read the full abstract

High Energy Astrophysical Phenomena

Radiation-hydrodynamics of star-disc collisions for quasi-periodic eruptions

Quasi-periodic eruptions (QPEs) are recently discovered transients of unknown nature occurring near supermassive black holes, which feature bright X-ray bursts separated by hours to days. A promising model for QPEs is the star-disc collisions model, where a star repeatedly interacts with an accretion disc around a black hole, creating shocks that expel dense outflows of gas from which radiation emerges. We investigate the dynamics of the star-disc collisions, the properties of the outflows, and the resulting radiation signatures. Our study focuses on the generic case where the star remains unperturbed by the collision and the stellar crossing time through the disc is sufficiently long for shocked gas to flow around the star. We performed a three-dimensional (3D) radiation-hydrodynamics simulation of the star-disc collision. The star was modeled as a solid, spherical body, and the interaction was simulated for a small, local section of the accretion disc. We found that star-disc collisions generate a nearly paraboloidal bow shock. The heating of gas is not confined to the column of gas directly ahead of the star but also extends laterally as the shock front expands sideways while traveling with the star. As the star crosses the disc, it injects momentum preferentially along its direction of motion, leading to an asymmetric redistribution of energy and momentum. As a result, two outflows emerge on opposite sides of the disc with different properties: the forward outflow expands faster, contains more mass, carries more energy, and is about twice as luminous as the backward outflow. Our findings suggest that the asymmetry in outflow properties and luminosity arises naturally from the collision dynamics, offering a possible explanation for the alternating strong-weak flare patterns observed in several QPE sources.
Read the full abstract

All authors Taj Jankovič, Clément Bonnerot, Sergey Karpov, Aleksej Jurca

High Energy Astrophysical Phenomena · Cross-listed from hep-ph

The ubiquitous flavor pendulum

A system of classical interacting spins can develop collective instabilities which, in the nonlinear regime, mimic the motion of a gyroscopic pendulum. Known as the flavor pendulum, this behavior appears among the collective modes of a dense neutrino plasma after a strong reduction of phase space through symmetry assumptions. It has been identified in homogeneous slow and fast flavor systems and, most recently, in single-wave solutions of the fast system. We explain the reasons for its ubiquitous appearance. We show that a system of three classical spins must always be pendular, or only two in the presence of an external field. Furthermore, such a system always defines a continuum of vectors with time-independent length. If these are identified as interacting spins, they immediately lead to the continuum cases of slow and fast flavor pendula. As another new insight, any of these spins can be chosen as the pendulum, periodically exchanging flavor with the rest of the system.
Read the full abstract

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

Dynamical response of twin stars to perturbations

If a strong first-order phase transition takes place at sufficiently high rest-mass densities in the equation of state (EOS) modelling compact stars, a new branch will appear in the mass-radius sequence of stable equilibria. This branch will be populated by stars comprising a quark-matter core and a hadronic-matter envelope, i.e., hybrid stars, which represent “twin-star” solutions to equilibria having the same mass but a fully hadronic EOS. While both branches are stable to linear perturbations, it is unclear which of the twin solutions is the “favoured” one, that is, which of the two configurations is expected to be found in nature. We assess this point by performing a large campaign of general-relativistic simulations aimed at assessing the response of compact stars on the two branches to perturbations of various strength. In this way, we find that, independently of whether the stars populate the hadronic or the twin branch, their response is characterised by a critical-perturbation strength such that the star will oscillate on the original branch for subcritical perturbations and migrate to the neighbouring branch for supercritical perturbations while conserving rest-mass. Because the critical values are different for stars with the same rest-mass but sitting on either branch, it is possible to define as favoured the part of the branch that has the largest critical perturbation, thus correcting the common wisdom that stellar models on the twin branch are the favoured ones. Interestingly, we show that the binding energies on the two branches can be used to deduce without simulations which of the stellar configurations is more likely to be found in nature.
Read the full abstract

Astrophysics of Galaxies

WLM: Dynamics of an isolated Dwarf Irregular Galaxy Under Ram Pressure in the Local Group

WLM is an archetypal dwarf irregular galaxy that has not experienced interactions with major Local Group galaxies within the past 8 Gyr. It has recently been shown that WLM is losing its gas due to ram pressure forces exerted by the surrounding intergalactic medium (IGM). In this work, we explore how ram pressure may also affect the WLM gas kinematics, and we show that its dynamics is especially perturbed at its outskirts, explaining the asymmetric rotation between the approaching and receding sides. Moreover, we have been able to decompose WLM in two main components, a compact one with a solid-body rotation that resembles a bar-like structure, and a more extended one with a characteristic double-horn profile suggesting an edge-on disk. The former is relatively unaffected by ram pressure while the latter has its dynamics considerably affected by ram pressure. This study shows that mass estimates of a dwarf galaxy like WLM should account for a full modeling of its dynamical components, especially accounting for its asymmetric rotation curve.
Read the full abstract

All authors Neel Kolhe, Francois Hammer, Yanbin Yang, Brenda Namumba, Laurent Chemin, Philippe Amram, Roger Ianjamasimanana, Claude Carignan

High Energy Astrophysical Phenomena

Gravitational-Wave Signals for Supernova Explosions of Three-Dimensional Progenitors

Core-collapse supernovae (SNe) are sources of gravitational waves (GWs) produced by hydrodynamical instabilities and highly time-dependent anisotropies of the neutrino radiation. In this work we analyze both contributions to the GW signal for two state-of-the-art three-dimensional (3D) SN models computed with the Prometheus-Vertex neutrino-hydrodynamics code. In contrast to the far majority of models analyzed for GWs so far, our core-collapse simulations were started with 12.28 M_sun (18.88 M_sun) progenitors, whose final hour (7 min) of convective oxygen-shell burning was computed in 3D and featured a vigorous oxygen-neon shell merger. The corresponding large-scale asymmetries in the oxygen layer are conducive to buoyancy-aided neutrino-driven explosions. The models were continuously evolved in 3D from the pre-collapse evolution until 5.11 s (1.68 s) after the core bounce. The GW signals result from the well-known dynamical phenomena in the SN core such as prompt postshock convection, neutrino-driven convection, the standing accretion shock instability, proto-neutron star oscillations, and anisotropic ejecta expansion. They do not exhibit any new or specific features that can be unambiguously connected to the powerful pre-collapse activity in the progenitors, but we identify interesting differences compared to results in the literature. We also discuss measurement prospects by interferometers, confirming that GW signals from future Galactic SNe will be detectable with existing and next-generation experiments working in the frequency range f ~ 1-2000 Hz.
Read the full abstract

All authors Alessandro Lella, Giuseppe Lucente, Daniel Kresse, Robert Glas, H.-Thomas Janka, Alessandro Mirizzi

Earth and Planetary Astrophysics

Saturn's Evolutionary History and Seismology: Survival of Deep Stably Stratified Regions in Evolutionary Models of Saturn Consistent with Ring Seismology

With recent advances in the modeling of the solar system giant planets, rapid progress has been made in understanding the remaining questions pertaining to their formation and evolution. However, this progress has largely neglected the significant constraints on the interior of Saturn's structure imposed by the observed oscillation frequencies in its rings. Here, we study initial conditions for Saturn's evolution that, after $4.56\;\mathrm{Gyr}$ of evolution, give rise to planetary structures admitting oscillation frequencies consistent with those observed via Saturn's ring seismology. Restricting our attention to models without compact rocky cores, we achieve simultaneous good agreement with most observed properties of Saturn at the level of current evolutionary models and with key frequencies in the observed oscillation spectrum. Our preliminary work suggests that Saturn's interior stably stratified region may be moderately less extended ($\sim 0.4$–$0.5R_{\rm Sat}$) than previously thought, which is important for reconciling the seismic constraints with evolutionary models. We also tentatively find that the deep helium gradients inferred by previous, static structural modelling of Saturn's ring seismology may not be required to reproduce the observed seismology data.
Read the full abstract

All authors Yubo Su, Janosz W. Dewberry, Roberto Tejada Arevalo, Ankan Sur, Adam Burrows

Cosmology and Nongalactic Astrophysics · Cross-listed from gr-qc

The emergent Big Bang scenario

This paper proposes a new avenue for understanding the cosmological singularity. The standard cosmological model contains a generic initial singularity usually referred to as the big bang. Herein, we present a novel idea to extend the description of our Universe beyond this limit. The proposal relies on rewriting physics in a purely Riemannian, i.e. locally Euclidean, four-dimensional space and the emergence of Lorentzian patches owing to the interaction of all matter fields to a clock field that is responsible for a signature change. If our Universe is contained within one of these patches, the initial singularity is replaced by a smooth boundary on which the signature of the physical metric flips. In this paper, we first define the model and draw the necessary conditions on its arbitrary functions for solutions to exist. Next, we prove the existence of solutions that lead to an emergent universe with a primordial (almost) de Sitter phase. To finish, we discuss the consequences of this construction for the universe on scales much larger than our observable Universe: a large “Euclidean sea” in which Lorentzian islands locally emerge and host an expanding universe potentially similar to ours. While speculative, this scenario has specific features that can be tested, and the present paper sets the basis for further phenomenological investigations.
Read the full abstract

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

Probing the Charged Hayward Black Hole in Dark Matter and String Cloud Environments through Shadow, Geodesics, and Quasinormal Spectrum

We construct a charged Bardeen black hole (BH) surrounded by perfect fluid dark matter (PFDM) and coupled to a cloud of strings (CS). The metric function combines the magnetic monopole charge from nonlinear electrodynamics, the PFDM logarithmic correction, and the CS parameter that renders the spacetime asymptotically non-flat. We analyze the horizon structure, identifying parameter ranges yielding non-extremal BHs, extremal configurations, and naked singularities. The null geodesics, photon sphere radius, and shadow are computed, revealing that both CS and PFDM enlarge the shadow. For neutral particle dynamics, we derive the specific energy, angular momentum, and innermost stable circular orbit location. Quasiperiodic oscillations (QPOs) are examined through the azimuthal, radial, and vertical epicyclic frequencies, where notably the azimuthal frequency is independent of the CS parameter. Scalar field perturbations governed by the Klein-Gordon equation yield an effective potential whose peak decreases with both parameters, yet the transmission and reflection probabilities respond oppositely to CS and PFDM variations. The greybody factor bounds are obtained using semi-analytical methods. Our results demonstrate that the distinct effects of $\alpha$ and $\beta$ on various observables could allow independent constraints on these parameters through shadow measurements, QPO timing, and gravitational wave ringdown observations.
Read the full abstract

Cosmology and Nongalactic Astrophysics

Non-Singular Bouncing cosmology from Phantom Scalar-Gauss-Bonnet Coupling: Reconstruction with Observational Insights

We examine non-singular bounce cosmology within the framework of a phantom scalar field coupled to the Gauss-Bonnet term in both non-viscous and bulk-viscous cases. Using the scale factor ansatz $\alpha(t)=\left(\frac{\alpha}{\eta}+t^2\right)^{\frac{1}{2 \eta}}$, we reconstruct the scalar field potential $V(t)$, and observe a smooth potential well centered at the bounce point. The resulting energy density, pressure, and equation-of-state parameter show NEC violation necessary for successful bounce, while viscosity controls post-bounce dynamics with a positive and smooth squared speed of sound. In contrast, for the non-viscous model, sharp divergences occur just at the bounce and continues to be negative in the expanding phase, which in turn emphasises the stabilising role of dissipative effects. The energy condition analysis indicates a temporary NEC and SEC violation in the viscous scenario, whereas its persistent violation within the non-viscous model suggests a continuous accelerated expansion. Observational viability is found through Bayesian MCMC fitting in regards to the Pantheon+ supernova data, with best-fit parameters providing a reduced chi-squared of $\chi_{red}^2 =0.995$ while the inflation observables derived from the reconstructed potential place our model predictions inside $68\%$ CL Planck 2018 confidence contours. Our findings suggest that bounce cosmologies could offer a physically reasonable and observationally acceptable alternative or pre-inflationary scenario, while highlighting the role that viscosity could play for a stable and smooth cosmological evolution.
Read the full abstract

Astrophysics of Galaxies

DeepDive: Tracing the early quenching pathways of massive quiescent galaxies at $z>3$ from their star-formation histories and chemical abundances

We investigate the chemical abundances and star-formation histories (SFH) of ten massive ($\mathrm{log}_{10} (M_{\star}/\mathrm{M}_{\odot})>10.5$) quiescent galaxies at $3<z<4$ using deep, medium-resolution spectroscopic data obtained as part of the JWST DeepDive Cycle 2 GO program. Our DeepDive sample demonstrates early formation and quenching times inferred from spectro-photometric fitting, with most galaxies having formed 50% of their stellar mass by $z \sim 5$, and quenching by $z \sim 4$, showing good agreement across the various SFH parameterizations explored in this work. Though they differ slightly between SFH parameterizations, the inferred formation timescales for the DeepDive sample span both rapid ($\lesssim$ 100 Myr) and more extended ($\gtrsim$ 200 Myr) episodes, corresponding to star formation occurring over a few to several dynamical times given their compact sizes and high densities at $z\sim3-4$. On average, massive quiescent galaxies at $3<z<4$ are $\alpha$-enhanced ($\langle [\alpha/\mathrm{Fe}]\rangle$= $0.22^{+0.22}_{-0.17}$), although there is strong diversity ($\sim0.3$ dex in scatter) among individual [$\alpha$/Fe] values. Our results for $\alpha$-enhancement are consistent with lower-redshift studies, implying weak evolution in [$\alpha$/Fe] from $z \sim 4$ to $z\sim 1$. The SFH timescales associated with the low [$\alpha$/Fe] measurements suggest longer formation timescales, potentially pointing to earlier enrichment by Type Ia supernovae, or metals preferentially being removed via outflows driven either by powerful early active galactic nuclei or supernovae. Overall, this work represents the first, statistically representative combined study of the star-formation histories and chemical abundances of massive quiescent galaxies at $z>3$.
Read the full abstract

All authors Massissilia L. Hamadouche, Katherine E. Whitaker, Francesco Valentino, Jacqueline Antwi-Danso, Kei Ito, Aliza Beverage, Pengpei Zhu, Gabriel Brammer, Vasily Kokorev, Gabriella de Lucia, William M. Baker, Marion Farcy, Anna Gallazzi, Steven Gillman, Rashmi Gottumukkala, Michaela Hirschmann, Christian Kragh Jespersen, Takumi Kakimoto, Minju M. Lee, Masato Onodera, Rhythm Shimakawa, Masayuki Tanaka, John R. Weaver, Po-Feng Wu

Cosmology and Nongalactic Astrophysics

The FLAMINGO Project: Exploring the X-ray–cosmic-shear cross-correlation as a probe of large-scale structure

Baryonic feedback processes associated with galaxy formation directly influence the large-scale structure by redistributing gas. Recent measurements of the kinetic Sunyaev-Zel'dovich (kSZ) effect and eROSITA stacks of X-ray emission from optically-selected galaxy clusters suggest that feedback from Active Galactic Nuclei (AGN) is more efficient at expelling gas from low-mass clusters than previously thought. The measurement of the cross-correlation between cosmic shear and diffuse X-ray emission provides a new probe of the distribution of gas in groups and clusters. We use the FLAMINGO cosmological, hydrodynamical simulations to examine the X-ray–cosmic-shear cross-correlation. The cross-correlation is most sensitive to the distribution of gas in haloes with masses $10^{14}\leq M_{200\mathrm{c}}/\mathrm{M}_{\odot}\leq10^{15}$. It is sensitive to the strength of feedback, but the effects of variations in cosmology and baryonic physics are largely degenerate. We compare the FLAMINGO predictions with the cross-correlation between cosmic shear from the Dark Energy Survey and the ROSAT all-sky X-ray map. We account for AGN that would remain unresolved by ROSAT using either the direct FLAMINGO predictions or by abundance matching to the observed (extrapolated) AGN luminosity function. We find that the fiducial FLAMINGO model is in excellent agreement with the ROSAT–DES-Y3 cross-correlation, while models assuming weaker feedback and models with the significantly stronger feedback suggested by kSZ and eROSITA stacks are both incompatible with the observed X-ray–cosmic-shear cross-correlation.
Read the full abstract

All authors William McDonald, Joop Schaye, Konrad Kuijken, John Helly, Joey Braspenning, Matthieu Schaller

2602.02484 · 2 Feb · Accepted for publication in MNRAS. The revised version accounts for attenuation of the X-ray flux from unresolved AGN, reducing their impact and yielding stronger constraints on baryonic feedback strength. The X-ray-cosmic-shear cross-correlations predicted by FLAMINGO's strongest baryonic feedback variations are now incompatible with the observations; the fiducial model is preferred instead

Instrumentation and Methods for Astrophysics

X-ray characterization of fully-depleted p-channel Skipper-CCDs for the DarkNESS mission

The Dark matter Nanosatellite Equipped with Skipper Sensors (DarkNESS) mission is a 6U CubeSat designed to search for X-ray lines from decaying dark matter using Skipper-CCDs. Thick, fully-depleted p-channel Skipper-CCDs provide low readout noise and high quantum efficiency for 1-10 keV X-rays, but their X-ray performance has not yet been demonstrated in the space environment. DarkNESS will operate in low-Earth orbit, where trapped protons induce displacement damage in the sensor that increases charge-transfer inefficiency and degrades the X-ray energy resolution. This work measures the X-ray line response of Skipper-CCDs before and after proton irradiation and quantifies the associated degradation. A sensor was exposed to 217 MeV protons at a fluence of 8.4 x 10^10 protons cm^-2, corresponding to a displacement-damage dose more than an order of magnitude above the three-year expectation for representative mid-inclination and Sun-synchronous low-Earth orbits. A 55Fe source was used to compare the energy resolution of the beam-exposed quadrant to adjacent unexposed quadrants and a non-irradiated reference sensor. These measurements provide a quantitative assessment of radiation-induced spectral degradation in Skipper-CCDs and enable an estimate of the end-of-life X-ray energy resolution expected for DarkNESS operation in low-Earth orbit.
Read the full abstract

All authors Phoenix Alpine, Ana M. Botti, Brenda A. Cervantes-Vergara, Claudio R. Chavez, Fernando Chierchie, Alex Drlica-Wagner, Juan Estrada, Erez Etzion, Michael Lembeck, Pilar López Maggi, Joseph Noonan, Brandon Roach, Nathan Saffold, Javier Tiffenberg

Solar and Stellar Astrophysics

Painting a Family Portrait of the Yellow Super- and Hypergiants in the Milky Way I. Constraining the Distances and Luminosities

Context. Distances to evolved massive stars in the Milky Way are not well constrained by Gaia parallaxes due to their brightness and variability. This makes it difficult to determine their fundamental stellar parameters, such as radius or luminosity, and infer their evolutionary states. Aims. We aim to improve the distance estimates of Yellow Hypergiants (YHGs) and Yellow Supergiants (YSGs) by identifying possible cluster and association memberships. Using these distances, we derived updated luminosities and revised their positions in the Hertzsprung-Russell diagram. Methods. We compiled from the literature a sample of 35 luminous yellow massive stars (YHGs and the most luminous YSGs). We used Gaia DR3 astrometry to identify possible membership in clusters and OB associations. We derived distances by combining the parallaxes of nearby co-moving stars. We independently validated these distances by comparing the stellar radial velocities to the Galactic H I kinematic map. We combined angular diameters and effective temperature values from the literature with the new distances to estimate luminosities. Results. We improved the distance estimates for 28 of the 35 stars through association with co-moving stellar groups. For an additional six stars, we provided distance estimates based on the H I kinematic map. For one star, the distance remains unclear. Most YSGs are members of young stellar populations, while the environments of the YHGs are more diverse, and for some of them, their origin populations remain unclear. We derived updated luminosities for a subset of 20 stars. Most YHGs have luminosities above log L/L = 5.4, while YSGs occupy a wider range of luminosities and the most luminous YSGs have luminosities similar to YHGs.
Read the full abstract

All authors A. Kasikov, A. Mehner, I. Kolka, A. Aret

Instrumentation and Methods for Astrophysics

Depth Calibration of Double-sided Strip Germanium Detectors for the Compton Spectrometer and Imager Satellite

Double-sided strip high-purity germanium detectors with three-dimensional position reconstruction capability have been developed over three decades, with space-based applications in high-energy astrophysics and heliophysics. Position resolution in three dimensions is key to reconstruction of Compton scattering events, including for the upcoming Compton Spectrometer and Imager (COSI) satellite mission. Two-dimensional position reconstruction is enabled by segmentation of the two detector faces into orthogonal strip contacts, enabling a pixelized analysis. The depth of an interaction cannot be measured directly but must be inferred from the charge collection time difference between the two faces of the detector. Here, we demonstrate for the first time the depth calibration of a detector with the COSI satellite geometry read out using an application specific integrated circuit (ASIC) developed for the COSI mission. In this work, we map collection time difference to depth using the Julia-based simulation package SolidStateDetectors$.$jl and validate it with comparison to the timing distributions observed in data. We also use simulations and data to demonstrate the depth resolution on a per-pixel basis, with >90% of pixels having <0.9 mm (FWHM) resolution at 59.5 keV and <0.6 mm (FWHM) resolution at 122.1 keV.
Read the full abstract

All authors Field R. Rogers, Sean N. Pike, Samer Alnussirat, Robin Anthony-Petersen, Steven E. Boggs, Felix Hagemann, Sophia E. Haight, Alyson Joens, Carolyn Kierans, Alexander Lowell, Brent Mochizuki, Albert Y. Shih, Clio Sleator, John A. Tomsick, Andreas Zoglauer

Solar and Stellar Astrophysics

Nonlinear interaction between dynamo-generated magnetic fields, mean flows and internal gravity waves in stellar stably-stratified layers: From 3D to 1D

Magnetic fields have been constrained at the surface of several massive and intermediate-mass stars, but their origin and properties in deep stellar radiative interiors are still debated, despite recent detections in the core of some red giant stars. Therefore, the modelling of AM transport in stellar radiative layers only relies on theoretical and numerical estimates of magnetic fields. Recent 3D numerical simulations show that a dynamo could occur in deep radiative regions. A realistic setup for understanding AM transport in such layers thus requires to take into account the mutual interactions of IGW and dynamo-generated magnetic field. We model the dynamics induced by IGW and dynamo in rotating radiative stellar layers using a simple description applicable to various evolutionary stages. As dynamo action and the propagation of IGW are 3D processes that have characteristic timescales short compared to periods associated with structural evolution of stars, we propose a mean-field 1D model by taking advantage of the dynamo coefficients computed from 3D spherical simulations. In this model, the necessary mean shear flow to trigger the dynamo results from the dissipation of monochromatic IGW generated in existing adjacent convective layers, which are expected to drive the formation of an oscillating rotational shear layer, the so-called Shear Layer Oscillation (SLO). In turn, magnetic effects can act on the mean flow through the Lorentz force. We show that the inclusion of magnetic fields adds up to the already very complex nonlinear problem and gives rise to the emergence of new dynamical regimes. Particularly, the fast SLO generated very close to the place where IGW are generated is perturbed by magnetic fields. This dynamical change can filter the wave energy spectrum transmitted towards further layers, with potential influence on the long-term evolution of the inner rotation.
Read the full abstract

All authors Florentin Daniel, Ludovic Petitdemange, Charly Pinçon, Kévin Belkacem, Bruno Longo, Christophe Gissinger

Cosmology and Nongalactic Astrophysics

Mitigating half-wave plate systematics at the map-making level: calibration requirements for LiteBIRD

Although half-wave plates (HWPs) are becoming a popular choice of polarization modulators for cosmic microwave background (CMB) experiments, their non-idealities can introduce systematic effects that should be carefully characterized and mitigated. One possible mitigation strategy is to incorporate information about the non-idealities at the map-making level, which helps to reduce the HWP-induced distortions of the reconstructed CMB. Nevertheless, the non-idealities can only be known with finite precision. In this paper we investigate the consequences of discrepancies between their true frequency profiles and those assumed by the map-maker. We present an end-to-end framework, including a blind component-separation step, and use it to translate these discrepancies into a bias on the tensor-to-scalar ratio, $r$, for the LiteBIRD satellite mission. We subsequently derive realistic and conservative measurement requirements for accurately characterizing the HWP non-idealities to ensure they do not compromise LiteBIRD's ambitious scientific goals. We find that the obtained results are robust against sky models with varying complexity.
Read the full abstract

All authors N. Raffuzzi, A. Carones, M. Monelli, S. Giardiello, L. Pagano, Y. Sakurai, H. Ishino, E. Allys, A. Anand, J. Aumont, A. J. Banday, G. Barbieri Ripamonti, R. B. Barreiro, N. Bartolo, S. Basak, A. Basyrov, A. Besnard, M. Bortolami, T. Brinckmann, F. Cacciotti, E. Calabrese, P. Campeti, F. Carralot, F. J. Casas, J. Chandran, K. Cheung, M. Citran, L. Clermont, F. Columbro, A. Coppolecchia, F. Cuttaia, P. de Bernardis, T. de Haan, M. De Lucia, P. Diego-Palazuelos, H. K. Eriksen, F. Finelli, C. Franceschet, U. Fuskeland, G. Galloni, M. Galloway, M. Gerbino, M. Gervasi, T. Ghigna, C. Gimeno-Amo, A. Gruppuso, M. Hazumi, S. Henrot-Versillé, L. T. Hergt, E. Hivon, K. Kohri, L. Lamagna, M. Lattanzi, C. Leloup, F. Levrier, A. I. Lonappan, M. López-Caniego, G. Luzzi, J. Macias-Perez, V. Maranchery, E. Martínez-González, S. Masi, S. Matarrese, T. Matsumura, S. Micheli, M. Migliaccio, G. Morgante, L. Mousset, R. Nagata, T. Namikawa, P. Natoli, A. Novelli, F. Noviello, A. Occhiuzzi, A. Paiella, D. Paoletti, G. Pascual-Cisneros, G. Patanchon, F. Piacentini, G. Piccirilli, M. Pinchera, G. Polenta, L. Porcelli, M. Remazeilles, A. Ritacco, M. Ruiz-Granda, L. Salvati, J. Sanghavi, V. Sauvage, D. Scott, M. Shiraishi, G. Signorelli, R. M. Sullivan, Y. Takase, L. Terenzi, M. Tomasi, M. Tristram, L. Vacher, B. van Tent, P. Vielva, S. Vinzl, I. K. Wehus, G. Weymann-Despres, E. J. Wollack

Earth and Planetary Astrophysics

What to make of the Earth's curiously intermediate land fraction?

Approximately two-thirds of the Earth, the only known inhabited planet, is covered in ocean. Why not 0.01% or 99.99%? It has been previously suggested that this may represent a certain degree of fine-tuning, and thus perhaps observers are a-priori more likely to develop on those rare worlds with nearly equal land-ocean ratios, such as our own. In this work, we take the single datum of the Earth and then use Bayesian inference to compare four models for the probability distribution of a planet becoming inhabited by observers as a function of land-fraction, $f$, which we classify as i) land-centric ii) ocean-centric iii) equi-centric and iv) indifference. We find that no model is strongly favoured over the others, but that 1) the land-centric model is disfavoured over all others, and, 2) the equi-centric model is favoured over all competitors. Further, we show that more extreme models with heavy tail-weighting are strongly disfavoured even when conditioned upon the Earth alone. For example, a land-centric model where the median planet has $f=0.82$ (or greater) is in strong tension with our existence. Finally, we consider the potential addition of more data via Mars or exoplanets. Should paleo-Mars have once harboured life and had $f<0.20$, then this would strongly favour the ocean-centric model for life, over a land-centric hypothesis. We show that strong evidence for/against the equi-centric model versus its competitors would likely require at least a dozen inhabited exoplanets, offering a well-motivated sample size for future experiments.
Read the full abstract

High Energy Astrophysical Phenomena

Evaluating the Contribution of Active Galactic Nuclei to the Diffuse High-Energy Neutrino Flux

The detection of high-energy neutrinos from NGC 1068 and TXS-0506+56 suggests that active galactic nuclei (AGN) may contribute significantly to the the diffuse neutrino flux measured by IceCube. Using 10 years of publicly available IceCube data, we performed a systematic population analysis of X-ray-bright and gamma-ray-bright AGN to evaluate the extent to which this diffuse flux could originate from these sources. We find that gamma-ray-bright blazars can account for no more than 16% of IceCube's total diffuse flux. Although we find no evidence of neutrino emission from gamma-ray-bright, non-blazar AGN, we cannot exclude the possibility that these sources contribute significantly to the diffuse flux. In contrast, we report (pre-trials) evidence of neutrino emission from several nearby, X-ray-bright, Seyfert-type AGN, including \mbox{NGC 1068} ($4.9\sigma$), SWIFT J1041.4-1740 ($2.6\sigma$), SWIFT J0202.4+6824A/B ($2.6\sigma$), SWIFT J0744.0+2914 (2.6$\sigma$), NGC 4151 ($2.5\sigma$), and NGC 3079 ($2.5\sigma$). Although not fully conclusive, these results suggest that IceCube may be detecting neutrinos from a larger population of Seyfert galaxies. The fact that these sources are not gamma-ray bright indicates that their neutrino production must be taking place in optically thick environments, such as in the coronae surrounding these galaxies' supermassive black holes. We also identify a $4.2\sigma$ correlation between the neutrinos detected by IceCube and members of the Swift-BAT catalog of X-ray-bright AGN, although this correlation is dominated by NGC 1068. We estimate that this class of sources contributes between 11.2% and the entirety of IceCube's total diffuse neutrino flux. These results strengthen the emerging case for the prevalence of gamma-ray-obscured AGN as significant sources of high-energy neutrinos.
Read the full abstract

High Energy Astrophysical Phenomena

A Unified Framework for 10 TeV to EeV Diffuse Neutrino Sky and KM3-230213A

Establishing a unified framework that simultaneously accounts for the wideband diffuse neutrino flux and the physical origin of individual ultra-high-energy (UHE) neutrino detections, including KM3-230213A, remains a pressing challenge in multi-messenger astrophysics. In intrinsically low-luminosity gamma-ray bursts (LL~GRBs) driven by shock breakouts (SBOs), the evolving physical conditions naturally produce a multicomponent neutrino flux extending from 10 TeV to the EeV scale. By integrating prompt and afterglow phases within a unified framework grounded in multiwavelength observations of representative events, we show that LL GRB population accounts for this broadband neutrino emission through a characteristic two-hump spectrum. In this framework, the prompt emission from GRB~060218-like events accounts for $\gtrsim 10\%$ of the diffuse flux at 100~TeV, while GRB~100316D-like afterglow configuration predicts a distinct flux peak near $10^{-9}\rm~GeV~cm^{-2}~s^{-1}~sr^{-1}$ at 100~PeV. This two-hump spectrum provides a high-energy component flux consistent with the 220 PeV KM3-230213A event, while the low-energy component contributes non-trivially to the observed diffuse neutrinos and supports the lack of individual low-energy counterparts. Furthermore, we utilize Fermi-LAT gamma-ray upper limits to place constraints on the source distance and luminosity of the event, assuming a GRB 100316D-like afterglow configuration. Ultimately, this framework identifies SBO-like LL~GRBs as a unifying origin for these phenomena, providing a physical link across the 10 TeV to EeV neutrino sky that is testable by next-generation observatories, including GRAND, IceCube-Gen2, and RNO-G.
Read the full abstract

High Energy Astrophysical Phenomena

Multi-Messenger Modeling of Low-Luminosity Gamma-Ray Bursts

Low-luminosity gamma-ray bursts (LL GRBs), a subclass of the most powerful transients in the Universe, remain promising sources of high-energy astrophysical neutrinos, despite strong IceCube constraints on typical long GRBs. In this work, a novel approach is introduced to study a sample of seven LL~GRBs with their multi-wavelength observations to investigate leptohadronic processes during their prompt emission phases. The relative energy densities in magnetic fields, non-thermal electrons, and protons are constrained, with the latter defining the cosmic-ray (CR) loading factor. Our results suggest that LL~GRBs exhibit diverse emission processes, as confirmed by a machine-learning analysis of the fitted parameters. Across the seven LL~GRBs, we find the posterior medians of the CR loading factor in the range of $\xi_p \sim 0.2$–$1.6$. GRB~060218 and GRB~100316D, the lowest-luminosity bursts ($L_{\gamma, \rm iso} \sim 10^{46}$-$10^{47}\rm~erg~s^{-1}$) consistent with the shock-breakout (SBO) scenario, yield the highest CR loading factor and therefore are expected to produce neutrinos more efficiently. Our model predicts the expected number of neutrino signals that are consistent with current limits but would be detectable with next-generation neutrino observatories. These results strengthen the case for LL~GRBs as promising sources of high-energy astrophysical neutrinos and motivate real-time searches for coincident LL~GRB and neutrino events. Next-generation X-ray and MeV facilities will be critical for identifying more LL~GRBs and strengthening their role in multi-messenger astrophysics.
Read the full abstract

Cosmology and Nongalactic Astrophysics

Reconstructing the largest scales of the Universe with field-level inference applied to the Quaia Quasar Catalogue

The recently released Quaia quasar catalogue, with its broad redshift range and all-sky coverage, enables unprecedented three-dimensional reconstructions of matter across cosmic time. In this work, we apply the field-level inference algorithm BORG to the Quaia catalogues to reconstruct the initial conditions and present-day matter distribution of the Universe. We employ a physics-based forward model of large-scale structure using Lagrangian perturbation theory, incorporating light-cone effects, redshift-space distortions, quasar bias, and survey selection effects. This approach enables a detailed and physically motivated inference of the three-dimensional density field and initial conditions over the entire cosmic volume considered. We analyse both the G < 20.0 (Quaia Clean) and G < 20.5 (Quaia Deep) samples, where G denotes the Gaia broad optical-band magnitude, imposing conservative sky cuts to ensure robustness against foreground contamination. The resulting reconstructions span a comoving volume of (10h^{-1} Gpc)^3 with a maximum spatial resolution of 39.1 h^{-1}Mpc, making this the largest field-level reconstruction of the observable Universe in terms of comoving volume to date. We validate our reconstructions through a range of internal and external consistency checks, including the cross-correlation of the inferred density fields with Planck CMB lensing, where we detect a signal at ~4\sigma significance. Beyond delivering high-fidelity data products, including posterior maps of initial conditions, present-day dark matter, and velocity fields, this work establishes a framework for exploiting quasar surveys in field-level cosmology.
Read the full abstract

All authors Adam Andrews, Arthur Loureiro, Jens Jasche, Stuart McAlpine, Guilhem Lavaux, Florent Leclercq

Cosmology and Nongalactic Astrophysics

Negentropy as Diagnostic of Cosmic Density Fields and Dynamical Dark Energy Models

We employ negentropy ($J$), defined as the difference between the information content of a non-Gaussian probability distribution and a Gaussian with identical variance, as an information-theoretic probe of non-Gaussianity in the cosmic density field. We quantify its sensitivity to dynamical dark energy by studying the evolution of $J(a)$ and its derivatives $\Gamma_1(a)$ and $\Gamma_2(a)$ across three parameterisation schemes: CPL, JBP, and BA. We determine the characteristic redshift $z_{NG}$, marking the epoch of maximal non-Gaussian structure formation, and the turnaround redshift $z_{TA}$, when information production transitions due to dark-energy domination, finding $z_{NG}\sim0.81$ and $z_{TA}\sim0.18$ for $\Lambda$CDM. Our diagnostics clearly discriminate between thawing and freezing quintessence models and phantom dark energy at low redshifts. Thawing models show small departures from $\Lambda$CDM, freezing models display higher $z_{TA}$, while phantom models exhibit lower $z_{TA}$, reflecting late-time evolution. We provide a practical prescription for measuring negentropy from discrete galaxy distributions, establishing a framework that can be applied to simulations and observations. This information-theoretic approach offers a robust and complementary tool for probing dark energy dynamics, enabling sensitive discrimination between evolving and cosmological-constant scenarios.
Read the full abstract

Astrophysics of Galaxies

PRISMS. UNCOVER-26185, a metal-poor SFG at z=10.05 with no evidence for a X-ray-luminous AGN

This work presents the first results of the PRImordial galaxy Survey with MIRI Spectroscopy (PRISMS), a JWST cycle 4 program (PID 8051) aimed at the characterization of a relatively large sample of ten galaxies about 500 Myr after the Big Bang. Here, we present deep (13.9 hours) spectroscopy with the MIRI LRS of the lensed galaxy UNCOVER-26185 at a redshift of z=10.054. It is a faint UV galaxy (UV absolut magnitude of -18.83 mag) previously identified as a X-ray luminous AGN. MIRI LRS detects the H$\beta$+[OIII]4960,5008 complex and H$\alpha$ emission line with a significance of 10$\sigma$ and 8$\sigma$, respectively, as well as the optical continuum emission at rest-frame 0.45 $\mu$m and 0.57 $\mu$m with a signal-to-noise ratio of 6-7. The UV-to-optical spectral energy distribution, combining continuum and emission lines, is compatible with: (i) a low stellar (A$_V$= 0.2) and nebular (A$_V$=0.0) extinction, (ii) a SFH composed by a young (7 Myr) starburst and an intermediate-age (65 Myr) stellar population, and (iii) a total stellar mass of 1.7$\times$10$^{8}$ M$_{\odot}$. The H$\alpha$-derived star-formation rate is 1.3 M$_{\odot}$ yr$^{-1}$. The low optical emission line ratios locate UNCOVER-26185 as the most metal-poor (Z = 0.04 Z$_{\odot}$), and as outlier with the lowest ionization (logU=-2.5) galaxy identified so far at redshifts above 9. With no evidence of an active galactic nuclei in the rest-frame UV-to-optical spectrum, UNCOVER-26185 has the properties of a metal-poor, main-sequence star-forming galaxy at redshift 10, with ISM and ionization properties very different than those of the already studied UV-bright galaxies at redshifts beyond 10. PRISMS is starting to explore the population of intermediate-UV luminosity galaxies at z=10, covering UV absolute magnitudes in the range of -17.9 to -20.5, fainter than those of UV-bright galaxies studied so far.
Read the full abstract

All authors J. Álvarez-Márquez, L. Colina, A. Crespo-Gomez, S. Kendrew, J. Zavala, R. Marques-Chaves, C. Prieto-Jiménez, Abdurro'uf, C. Blanco-Prieto, L. A. Boogaard, M. Castellano, A. Fontana, Y. Fudamoto, S. Fujimoto, M. García-Marín, Y. Harikane, S. Harish, T. Hashimoto, T. Hsiao, E. Iani, A. K. Inoue, D. Langeroodi, R. Lin, J. Melinder, L. Napolitano, G. Ostlin, P. G. Pérez-González, P. Rinaldi, B. Rodríguez Del Pino, P. Santini, Y. Sugahara, T. Treu, A. Varo-O'ferral, G. Wright

Astrophysics of Galaxies

PRISMS. U37126, a very blue, ISM-naked starburst at z=10.255 with nearly 100% Lyman continuum escape fraction

We present very deep (~11h) JWST/MIRI low-resolution spectroscopy of the rest-frame optical emission of U37126, a UV-bright (M_UV ~ -20), mildly lensed ($\mu\simeq 2.2$) galaxy at z=10.255. The continuum emission is well detected in both NIRSpec and MIRI spectra, yet no nebular recombination or metal emission lines are observed (EW(Hbeta+[OIII])<300A and EW(Halpha)<400A, at 3sigma). Combined with the exceptionally blue UV continuum slope, beta_UV ~ -2.9, and weak/flat Balmer break, these constraints indicate a stellar population dominated by very young and massive stars with a strongly suppressed nebular contribution. Comparisons with synthetic stellar population models indicate that U37126 requires both a very high ionizing photon production efficiency, log(Xi_ion / Hz erg^-1) ~ 25.75, and a nearly unit LyC escape fraction, of fesc>86% (3sigma) based on Halpha flux limit and fesc=0.94+/-0.06 derived independently from SED fitting. The best-fit SED yields a (de-lensed) stellar mass of Mstar ~ 10^7.8 Msun and a star-formation rate of SFR~10Msun/yr (sSFR~160 Gyr^-1), that along with its very compact size, reff~61pc, yields very high stellar mass and star-formation-rate surface densities, Sigma_M ~ 3x10^3 Msun/pc^2 and Sigma_SFR ~ 400 Msun/yr/kpc^2. Together with the lack of detectable nebular emission, these properties suggest that U37126 is undergoing an “ISM-naked” starburst phase, possibly driven by an extremely efficient gas-to-star conversion followed by strong feedback that has cleared the remaining gas from its stellar core, allowing most LyC photons to escape. Finally, we show that even a small fraction of galaxies like U37126 (~ 3%-6%), with extreme LyC production and escape, could contribute disproportionately (~ 50%-100%) to the ionizing photon budget during cosmic reionization.
Read the full abstract

All authors R. Marques-Chaves, J. Álvarez-Márquez, L. Colina, S. Kendrew, Abdurro'uf, C. Blanco-Prieto, L. A. Boogaard, M. Castellano, K. I. Caputi, A. Crespo-Gomez, A. Fontana, Y. Fudamoto, S. Fujimoto, M. García-Marín, Y. Harikane, S. Harish, T. Hashimoto, T. Hsiao, E. Iani, A. K. Inoue, D. Langeroodi, R. Lin, J. Melinder, L. Napolitano, G. Ostlin, P. G. Pérez-González, C. Prieto-Jiménez, P. Rinaldi, B. Rodríguez Del Pino, P. Santini, Y. Sugahara, A. Varo-O'ferral, G. Wright, J. Zavala

Instrumentation and Methods for Astrophysics

Surface Interactions in Photon Monte Carlo Simulations

We implement a comprehensive simulation of photon surface interactions using a Monte Carlo approach. This is effective in simulating the interaction of light with telescope mirrors and lenses. We use a full electromagnetic solution to simulate the wavelength and angular dependence at surfaces. This includes bare interfaces, monolayer interfaces, protected layer coatings, and multilayer coatings. We handle special cases when multilayer data is incomplete or when there is photo-conversion in the interface as with sensors. We implement interactions with surface micro-roughness and predict the corresponding angular distribution using a Monte Carlo implementation of the Harvey-Shack scatter theory for a microroughness power spectrum. Finally, we simulate surface interaction with contamination from dust or condensation using Mie scattering applied efficiently to individual contaminants. The combination of these implementations can efficiently simulate rough to polished surfaces of arbitrary materials that are fully cleaned or dusty. The observational consequences includes complex wavelength and spatial-dependent photometric errors, the dominant effect of the wings of point-spread-functions, dust rings, and wavelength and angle-dependent throughput losses. We find agreement with the point-spread-function wings of WIYN ODI observations of bright stars and properties of dust rings, and demonstrate the ability to disentangle mirror microroughness from the turbulence PSF patterns. The comprehensive numerical implementation of surface interactions has wide applicability in non-astronomical applications as well.
Read the full abstract

All authors J. R. Peterson, D. Valls-Gabaud, A. Dutta, C. Kim, G. Sembroski

Solar and Stellar Astrophysics

Evolution of Coronal Mass Ejections in Different Data-Driven Solar Wind Conditions

Numerical models of the solar wind and coronal mass ejections (CMEs) utilize photospheric magnetic field observations to prescribe the inner boundary conditions for the plasma solutions. These magnetic field data are available to the community through various observational instruments, prepared via different methodologies and/or flux-transport models. The solar wind solution driven by these maps provides the ambient plasma environment into which CMEs travel, coupling, and interacting with the surrounding plasma and governing the CME evolution and propagation in the solar corona and inner heliosphere. In this work, we use different input magnetic field maps for the same time period to drive the global Alfven Wave Solar atmosphere Model (AWSoM). We obtain the ambient solar wind conditions and compare the plasma properties and magnetic morphology in the coronal domain to study the influence of the input maps. To understand how the resulting coronal solutions impact CMEs, we launch eruptions described by analytical flux ropes into these data-driven solutions and compare their evolution in the coronal domain (up to 24 solar radii radially). The CMEs achieve varying speeds, deceleration rates, propagation directions, mass and energies while coupling with the background solar wind. We quantify these differences to show that the different input driving maps can significantly impact the simulated CME propagation in the solar wind plasma. This also highlights the importance of understanding the uncertainties associated with data-driven modeling that become increasingly important in operational models and space weather prediction.
Read the full abstract

All authors Nishtha Sachdeva, Zhenguang Huang, Gabor Toth, Hongfan Chen, Ward B. Manchester, Bart van der Holst

Earth and Planetary Astrophysics

Helium escape in context: Comparative signatures of four close-in exoplanets

Observations of escaping atmospheres on close-in exoplanets show a wide range in the strength and morphology of He I 10830 A and H I absorption. Scaling relations attempt to link the He I signal to XUV irradiation, mass loss, and bulk planetary parameters. We test these relations with a comparative analysis of HD209458b, HD189733b, HD149026b, and GJ1214b using a 1D hydrodynamic, multi-species, full-atmosphere escape model. For the benchmark HD209458b, our previously validated solution reproduces the observed He I and Ha transit depths without imposing composition constraints. HD189733b exhibits comparable He I depths, but the broadest reported profiles require ~12 km/s of additional non-thermal broadening, whereas more recent measurements are narrower, consistent with our predictions. For HD149026b, despite similar system properties, our model shows that higher gravity suppresses escape and enhances diffusive separation, depleting helium at high altitudes and yielding extremely weak He I absorption. For the sub-Neptune GJ1214b, H/He-only models overestimate He I absorption; including H2 and its ions (H2+, H3+, HeH+) lowers the escape rate and modifies the ion/electron balance, reducing the metastable helium densities. Compared against scaling relations, HD189733b observations and our HD149026b prediction fall below the trend, whereas some observations of HD209458b and GJ1214b are consistent; however, the observed transit depths are variable. Across all targets, we find diffusive separation of helium and hydrogen, which may explain why sub-solar He/H ratios are often required in simplified models. We conclude that interpreting He I and Ha absorption requires first-principles models that include self-consistent temperature and velocity profiles, multi-species transport, and molecular chemistry.
Read the full abstract

All authors Anna Ruth Taylor, Tommi T. Koskinen, Chenliang Huang, Anthony Arfaux, Panayotis Lavvas

Instrumentation and Methods for Astrophysics

The Nulling Interferometry Cryogenic Experiment (NICE): Architecture, requirements, and preliminary warm precursor results

The success of the Large Interferometer For Exoplanets (LIFE) space mission depends on measuring the faint mid-infrared emission spectra of exoplanets while suppressing the glare of a host star. This requires an instrument capable of high-contrast nulling interferometry with exceptional sensitivity. While previous testbeds have proven the principle of deep, stable nulls, they have not combined high contrast with the high throughput and cryogenic operation required for LIFE. Here, we present the architecture of the Nulling Interferometry Cryogenic Experiment (NICE), a mid-infrared nulling testbed, to increase the technological readiness of LIFE. We derive the laboratory requirements necessary to validate the LIFE beam combiner and present the optical design of NICE. Finally, we report results from the ambient \enquote{Warm Bench} precursor, which has successfully demonstrated the required null depth ($< 10^{-5}$) using a polarized narrowband 4.7 um source, and the required throughput (> 17%) using one of the two nulling channels.
Read the full abstract

All authors Thomas Birbacher, Jonah T. Hansen, Felix A. Dannert, Germain Garreau, Adrian M. Glauser, Ryan Meierhofer, Julio Pino Jiménez, Mohanakrishna Ranganathan, Sascha P. Quanz

High Energy Astrophysical Phenomena

Search for ultra-high energy neutrons from Galactic sources with the Pierre Auger Observatory

Deflections in the propagation of charged ultra-high-energy cosmic rays (UHECRs) caused by magnetic fields make the identification of their sources challenging. On the other hand, the arrival directions at Earth of neutrons point directly to their origin. The emission of UHECRs from a source is expected to be accompanied by the production of neutrons in its vicinity through interactions with ambient matter and radiation. Since free neutrons travel a mean distance $d/\text{kpc}=9.2(E/\text{EeV})$ before decaying, a neutron flux in the EeV range could be detected on Earth from sources of UHECRs in our Galaxy. Using cosmic-ray data from the Phase\,I of the Surface Detector of the Pierre Auger Observatory, we search for neutron fluxes from Galactic candidate sources. We select more than 1000 objects of astrophysical interest, stacking them into target sets. The targets all have declinations within the exposure of the Observatory, ranging from $-90^\circ$ up to $+45^\circ$ for energies above 1 EeV (and up to $+20^\circ$ for energies down to 0.1 EeV). Given that a neutron air shower is indistinguishable from a proton one, there is a significant background due to cosmic rays. A neutron flux from the direction of a candidate source would be identified by a celestial density of events that significantly exceeds the expected density of cosmic rays for that direction. No significant excess is found at any tested target direction, and an upper limit on the neutron flux is calculated for each candidate source.
Read the full abstract

All authors The Pierre Auger Collaboration: A. Abdul Halim, P. Abreu, M. Aglietta, I. Allekotte, K. Almeida Cheminant, A. Almela, R. Aloisio, J. Alvarez-Muñiz, A. Ambrosone, J. Ammerman Yebra, L. Anchordoqui, B. Andrada, L. Andrade Dourado, L. Apollonio, C. Aramo, E. Arnone, J.C. Arteaga Velázquez, P. Assis, G. Avila, E. Avocone, A. Bakalova, A. Baluta, F. Barbato, A. Bartz Mocellin, C. Berat, M.E. Bertaina, M. Bianciotto, P.L. Biermann, V. Binet, K. Bismark, T. Bister, J. Biteau, J. Blazek, J. Blümer, M. Boháčová, D. Boncioli, C. Bonifazi, N. Borodai, J. Brack, P.G. Brichetto Orchera, A. Bueno, S. Buitink, M. Büsken, A. Bwembya, K.S. Caballero-Mora, S. Cabana-Freire, L. Caccianiga, F. Campuzano, J. Caraça-Valente, R. Caruso, A. Castellina, F. Catalani, G. Cataldi, L. Cazon, M. Cerda, B. Čermáková, A. Cermenati, K. Cerny, J.A. Chinellato, J. Chudoba, L. Chytka, R.W. Clay, A.C. Cobos Cerutti, R. Colalillo, R. Conceição, G. Consolati, M. Conte, F. Convenga, D. Correia dos Santos, P.J. Costa, C.E. Covault, M. Cristinziani, C.S. Cruz Sanchez, S. Dasso, K. Daumiller, B.R. Dawson, R.M. de Almeida, E.-T. de Boone, B. de Errico, J. de Jesús, S.J. de Jong, J.R.T. de Mello Neto, I. De Mitri, D. de Oliveira Franco, F. de Palma, V. de Souza, E. De Vito, A. Del Popolo, O. Deligny, N. Denner, K. Denner Syrokvas, L. Deval, A. di Matteo, C. Dobrigkeit, J.C. D'Olivo, L.M. Domingues Mendes, Y. Dominguez Ballesteros, Q. Dorosti, R.C. dos Anjos, J. Ebr, F. Ellwanger, R. Engel, I. Epicoco, M. Erdmann, A. Etchegoyen, C. Evoli, H. Falcke, G. Farrar, A.C. Fauth, T. Fehler, F. Feldbusch, A. Fernandes, M. Fernández Alonso, B. Fick, J.M. Figueira, P. Filip, A. Filipčič, T. Fitoussi, B. Flaggs, T. Fodran, A. Franco, M. Freitas, T. Fujii, A. Fuster, C. Galea, B. García, C. Gaudu, P.L. Ghia, U. Giaccari, C. Glaser, F. Gobbi, F. Gollan, G. Golup, M. Gómez Berisso, P.F. Gómez Vitale, J.P. Gongora, J.M. González, N. González, D. Góra, A. Gorgi, M. Gottowik, F. Guarino, G.P. Guedes, L. Gülzow, S. Hahn, P. Hamal, M.R. Hampel, P. Hansen, V.M. Harvey, A. Haungs, T. Hebbeker, C. Hojvat, J.R. Hörandel, P. Horvath, M. Hrabovský, T. Huege, A. Insolia, P.G. Isar, M. Ismaiel, P. Janecek, V. Jilek, K.-H. Kampert, B. Keilhauer, A. Khakurdikar, V.V. Kizakke Covilakam, H.O. Klages, M. Kleifges, J. Köhler, F. Krieger, M. Kubatova, N. Kunka, B.L. Lago, N. Langner, N. Leal, M.A. Leigui de Oliveira, Y. Lema-Capeans, A. Letessier-Selvon, I. Lhenry-Yvon, L. Lopes, J.P. Lundquist, M. Mallamaci, S. Mancuso, D. Mandat, P. Mantsch, F.M. Mariani, A.G. Mariazzi, I.C. Mariş, G. Marsella, D. Martello, S. Martinelli, M.A. Martins, H.-J. Mathes, J. Matthews, G. Matthiae, E. Mayotte, S. Mayotte, P.O. Mazur, G. Medina-Tanco, J. Meinert, D. Melo, A. Menshikov, C. Merx, S. Michal, M.I. Micheletti, L. Miramonti, M. Mogarkar, S. Mollerach, F. Montanet, L. Morejon, K. Mulrey, R. Mussa, W.M. Namasaka, S. Negi, L. Nellen, K. Nguyen, G. Nicora, M. Niechciol, D. Nitz, D. Nosek, A. Novikov, V. Novotny, L. Nožka, A. Nucita, L.A. Núñez, S.E. Nuza, J. Ochoa, M. Olegario, C. Oliveira, L. Östman, M. Palatka, J. Pallotta, S. Panja, G. Parente, T. Paulsen, J. Pawlowsky, M. Pech, J. Pękala, R. Pelayo, V. Pelgrims, E.E. Pereira Martins, C. Pérez Bertolli, L. Perrone, S. Petrera, C. Petrucci, T. Pierog, M. Pimenta, M. Platino, B. Pont, M. Pourmohammad Shahvar, P. Privitera, C. Priyadarshi, M. Prouza, K. Pytel, S. Querchfeld, J. Rautenberg, D. Ravignani, J.V. Reginatto Akim, A. Reuzki, J. Ridky, F. Riehn, M. Risse, V. Rizi, E. Rodriguez, G. Rodriguez Fernandez, J. Rodriguez Rojo, S. Rossoni, M. Roth, E. Roulet, A.C. Rovero, A. Saftoiu, M. Saharan, F. Salamida, H. Salazar, G. Salina, P. Sampathkumar, N. San Martin, J.D. Sanabria Gomez, F. Sánchez, E. Santos, F. Sarazin, R. Sarmento, R. Sato, P. Savina, V. Scherini, H. Schieler, M. Schimassek, M. Schimp, D. Schmidt, O. Scholten, H. Schoorlemmer, P. Schovánek, F.G. Schröder, J. Schulte, T. Schulz, S.J. Sciutto, M. Scornavacche, A. Sedoski, S. Sehgal, S.U. Shivashankara, G. Sigl, K. Simkova, F. Simon, R. Šmída, S. Soares Sippert, P. Sommers, R. Squartini, M. Stadelmaier, S. Stanič, J. Stasielak, P. Stassi, S. Strähnz, M. Straub, T. Suomijärvi, A.D. Supanitsky, Z. Svozilikova, Z. Szadkowski, F. Tairli, M. Tambone, A. Tapia, C. Taricco, C. Timmermans, O. Tkachenko, P. Tobiska, C.J. Todero Peixoto, B. Tomé, A. Travaini, P. Travnicek, C. Trimarelli, M. Tueros, M. Unger, R. Uzeiroska, L. Vaclavek, M. Vacula, I. Vaiman, J.F. Valdés Galicia, L. Valore, P. van Dillen, E. Varela, V. Vašíčková, A. Vásquez-Ramírez, D. Veberič, I.D. Vergara Quispe, S. Verpoest, V. Verzi, J. Vicha, S. Vorobiov, J.B. Vuta, C. Watanabe, A.A. Watson, A. Weindl, M. Weitz, L. Wiencke, H. Wilczyński, B. Wundheiler, B. Yue, A. Yushkov, E. Zas, D. Zavrtanik, M. Zavrtanik