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)

Friday 20 February

Cosmology and Nongalactic Astrophysics

Cosmological Constraints on Temperature-Dependent Interaction between Dark Matter and Neutrinos

We study the influence of the temperature-dependent interaction between dark matter (DM) and neutrinos on the measurement of cosmological parameters. We pay attention to the neutrino mass effects, so that the derivation of Boltzmann equations needs to specify the concrete form of interaction. We work in a model in which the DM-neutrino scatterings are induced by a dimension-six operator, and present the details for deriving the full Boltzmann hierarchy for DM and neutrinos, including a novel method to obtain the fluid approximation for modes entering the horizon. It is shown that our interaction can induce the dark acoustic oscillation in the DM-neutrino fluid, leaving distinct signatures on the CMB and matter power spectra. By using the latest CMB and BAO datasets from Planck, DESI and ACT, the constraint on today's DM-neutrino interaction parameter for the normal neutrino mass ordering reaches $u^0_{\chi-\nu} \lesssim {\cal O}(10^{-13})$, nearly nine orders stronger than that for temperature-independent case in the literature. This can be understood by noting that the scattering cross section increases nearly quadratically with cosmological temperature in the early universe, leading to enhanced effects. We have investigated alternative scenarios with different neutrino mass assumptions. In particular, models with degenerate neutrino masses give rise to weaker constraint of $u^0_{\chi-\nu} \lesssim {\cal O}(10^{-11})$, showing the importance to incorporate the realistic neutrino mass ordering in the fits. Finally, when employing the logarithmic flat prior for $u^0_{\chi-\nu}$, we have shown hints to a nonzero interaction at $95\%$ CL by combining Planck, DESI and ACT data.
Read the full abstract

Astrophysics of Galaxies

Single- and double-headed odd radio circles in the LOFAR surveys

Deep radio surveys are now producing catalogs with millions of radio sources. Radio sources can have complex morphologies that depend on both the production mechanisms and the local environment. Recently, an unusual class of circular radio sources named "odd radio circles" (ORCs) were identified. They have sizes of about 1 arcmin, a circular/elliptical shape and appear edge-brightened. Subsequent observations have suggested that this class may comprise a variety of sources. Despite various attempts to explain them, their origin remains unclear. The main goal of this work is to increase the number of known ORCs and derive common characteristics that can help identify the origin of these sources. We searched the LOFAR Two Metre Sky Survey (LoTSS) Data Release 3 (DR3) at 144 MHz for ORCs using a combination of parameter filtering on catalog entries and visual inspection. We then identified possible optical counterparts and derived information such as redshift, physical size, and spectral index using further radio data at 54 and 1400 MHz. We isolated 18 sources with ORC structures. Four of these are double-headed ORCs (ORCs with two rings), and two are new discoveries. We also found five new single-headed ORCs and nine candidate ORCs. With this work we significantly expand the population of known ORCs. Our findings confirm that ORCs are a rare and heterogenous population of radio sources. We confirm the association with large ellipticals in most cases, and we note a relation between the ORCs' physical sizes and their integrated spectral indexes, with small ORCs avoiding steep spectra.
Read the full abstract

All authors F. De Gasperin, M. Brüggen, T. Pasini, H. Andernach, A. Ivleva, C. Spingola, L. M. Böss, J. Callingham, M. Hardcastle, G. Di Gennaro

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

GR-Athena++: Binary Neutron Star Merger Simulations with Neutrino Transport

We present general-relativistic radiation magnetohydrodynamics simulations of binary neutron star mergers performed with GR-Athena++. Neutrino transport is treated using a moment-based, energy-integrated scheme (M1), augmented by neutrino number density evolution (N0). Our implementation is validated through an extensive suite of standard tests and demonstrated to perform robustly under adaptive mesh refinement. As a first application, we simulate the gravitational collapse of a uniformly rotating, magnetized neutron star, demonstrating stable radiation evolution through apparent-horizon formation using a novel excision technique based on the tapering of state vector evolution inside the horizon. To further test robustness in highly dynamic environments, we apply our code to two demanding binary neutron star merger scenarios. We investigate a long-lived remnant with the DD2 equation of state, evolved with full general-relativistic magnetohydrodynamics and M1 neutrino transport. Following this, a gravitational collapse scenario with the SFHo equation of state is explored. We showcase long-term stable evolution on neutrino cooling time-scales, demonstrating robust handling of excision and stable evolution of the post-collapse accretion phase in three-dimensional mergers with magnetic fields and neutrino radiation.
Read the full abstract

All authors Boris Daszuta, Sebastiano Bernuzzi, Maximilian Jacobi, Eduardo M. Gutiérrez, Peter Hammond, William Cook, David Radice

Instrumentation and Methods for Astrophysics

Chromaticity-Optimized Antenna Design and Bayesian Foreground Validation for the CANTAR Global 21 cm Experiment

Detecting the global 21 cm signal from the epoch of reionization remains a major observational challenge due to bright foregrounds and instrumental systematics. As part of the Colombian Antarctic Telescopes for 21 cm Absorption during Reionization (CANTAR) initiative, we present a simulation and analysis framework to evaluate antenna chromaticity, optimize instrument design, and assess site suitability for global 21 cm experiments. Using frequency-dependent beam models and Haslam-based sky maps, we compute dynamic spectra for the EDGES blade dipole and a set of dipole and novel monopole antennas optimized via particle swarm optimization. The optimized designs exhibit improved spectral smoothness compared to EDGES, particularly in the 70-120 MHz range. We also evaluate latitude-dependent sky brightness and identify mid-latitude sites (-40{\deg} to +5{\deg}) as optimal for foreground suppression. We apply Bayesian inference together with posterior predictive model validation to the publicly released EDGES data, assessing statistical consistency rather than hypothesis testing or model comparison. We find that physically motivated foreground and ionospheric models are statistically consistent with the data only when a 21 cm absorption feature is excluded. From the validated posterior, we generate a statistically validated ensemble of foreground corrections for use in beam-sky simulations. These results support a two-phase strategy for CANTAR: Antarctic deployments for calibration and testing, and future science operations at mid-latitude sites. Our framework provides a validated path toward robust foreground modeling, antenna design, and systematics control for global 21 cm signal detection.
Read the full abstract

All authors Michelle Mora, German Chaparro, Juan D. Guerrero, Catalina Alzate, Juan P. Urrego, Jimena Giraldo, Pablo Cuartas-Restrepo, Julian Rodriguez-Ferreira, Oscar Restrepo

Cosmology and Nongalactic Astrophysics

CMB anisotropies from cosmic (super)strings in light of ACT DR6

We present updated constraints on cosmic string and superstring parameters derived from Cosmic Microwave Background (CMB) anisotropies. The constraints are obtained via Markov Chain Monte Carlo (MCMC) analyses of the full Planck temperature and polarization data combined with the Atacama Cosmology Telescope (ACT) Data Release 6 (DR6). For ordinary cosmic strings, we constrain the string tension $G\mu$, the string wiggliness parameter $\alpha$, and the self-chopping efficiency $\tilde{c}$. For cosmic superstrings, we constrain the fundamental string tension $G\mu_F$, the string coupling $g_s$, and a parameter $w$ describing the volume of the compact extra dimensions. In both cases, we find significantly tighter bounds on the string tension compared to previous analyses, obtaining $2\sigma$ upper limits of $G\mu < 3.66\times10^{-8}$ and $G\mu_F < 1.38\times10^{-8}$. We also discuss the significant prior-dependence of these results. The computational pipeline used in this work, including a modified version of CAMB capable of computing CMB anisotropies sourced by any active network described via unequal-time correlators, is released publicly as CAMBactive \cite{Raidal_CAMBactive_CAMB_extension_2026}.
Read the full abstract

All authors Juhan Raidal, Anastasios Avgoustidis, Edmund Copeland, Adam Moss

Instrumentation and Methods for Astrophysics

Mathematical derivation and verification of the amplitude of LISA's interferometric signals on an ultra-stable interferometer testbed

The Laser Interferometer Space Antenna (LISA) mission aims to detect gravitational waves by interferometrically measuring the change of separation between free-falling test masses (TMs). LISA's interferometers must deliver pm/rtHz sensitivity while accommodating beam tilts up to 1 mrad at the photodiodes, which degrade the interferometric amplitude and increase the induced readout noise coupling. This paper uses an analytical framework developed by the authors in a previous work, based on minimal and justified approximations, that relates beam tilt to the resulting heterodyne signal amplitude in a generic two-beam interferometer with circular-area photodiodes (PDs). A set of interferometric topologies is analyzed, all of high relevance for LISA. We derive the exact amplitude response for an infinite detector and a closed-form approximation for finite detectors, and we validate both against numerical simulations and experimental measurements on an ultra-stable LISA-representative testbed. We then use this model to quantify the phase-noise amplification arising from reduced signal-to-noise ratio (SNR) under tilt, showing that curvature mismatches between the interfering beams substantially enhance this effect. Finally, we introduce a compact function that captures the angular dependence of correlated and uncorrelated phase noises in quadrant photodiode (QPD)-based readouts. Here, a new noise feature, caused by wavefront curvature mismatch, is predicted and measured for the first time. These results indicate that controlling wavefront curvature mismatch in the test mass interferometer (TMI) is essential to limit excess phase noise. The models and results derived in this paper, although originating in the context of LISA, are general and can be applied to any interferometric topology undergoing tilts with pivot on the detector plane.
Read the full abstract

All authors Alvise Pizzella, Lennart Wissel, Miguel Dovale-Alvarez, Pablo Martinez Cano, Rodrigo Garcia Alvarez, Christoph Bode, Juan Jose Esteban Delgado, Gerhard Heinzel

Cosmology and Nongalactic Astrophysics

First Dark Photon Search Results from the Dandelion Experiment

This paper presents the first results from the Dandelion experiment, a directional detection, which searches for 1 meV dark photon dark matter. We use a spherical mirror to convert dark photons into standard millimeter-wavelength photons that can then be detected with an array of 221 Kinetic Inductance Detectors (KIDs) cooled down to 150 mK within the KISS-NIKA camera ( (KIDs Interferometric Spectral Surveyor) and operating between 150 and 350 GHz. We used 1480 minutes of data to search for the signal of dark photons in the KID detectors, which is expected to be modulated due to the Earth's rotation. Our main challenge was to deal with a large background from room temperature and stray-light fluctuations. We used a de-correlation analysis to remove these background fluctuations. Templates of the background fluctuations were constructed from a Principal Component Analysis decomposition of detector measurements outside the expected Field of View trajectory of dark photons. We found that the dark photon signal was consistent with zero, giving a new upper limit on the dark photon's kinetic mixing, $\chi$, with masses between 0.6 meV and 1.4 meV. These are the first constraints on dark photons as a dark matter candidate using an array of KIDs at millimeter wavelength.
Read the full abstract

All authors I. Ourahou, S. Savorgnano, C. Beaufort, M. Bastero-Gil, J. Bounmy, A. Catalano, J. Macias-Perez, D. Santos, C. Smith, F. Naraghi, D. Tourres, F. Vezzu

Astrophysics of Galaxies

The GUAPOS project – VII: Physical structure and molecular environment of the G31.41+0.31 HII region

Ionised regions around OB-type stars are formed at an early stage of their evolution and are important to investigate the formation process of these objects. However, so far only few observations of their physical structure and interaction with the parental molecular cloud have been made. The high resolution and sensitivity of new instruments such as ALMA and the upgraded VLA allow us to fill this gap in our knowledge. We investigate the well known core-halo ultracompact HII region G31.41+0.31 and the surrounding molecular clump with the aim to determine the density and temperature of both the ionised and neutral gas, and possibly obtain a 3D picture of their spacial distribution. We take advantage of the full-band frequency coverage at 3 mm obtained with ALMA for the GUAPOS project to image the emission of a plethora of hydrogen recombination lines towards the G31.41+0.31 HII region as well as several molecular transitions which are tracers of medium-density ($\sim$$10^4$–$10^6$ cm$^{-3}$) gas. The line data are complemented by continuum measurements obtained with the VLA at 1 cm and 7 mm. By fitting these lines also using a model that takes into account non-LTE effects we can investigate the density and temperature structure and the velocity field of the region. Our findings, based on a model fit accounting for non-LTE effects, indicate that the electron temperature of the HII region is mostly spanning a range between 5000 and 6000 K, while the density varies between 2500 and 7500 cm$^{-3}$. All in all, the distribution of these parameters as well as the corresponding velocity field hint at a cometary shaped HII region expanding away from the observer to the NW. The molecular gas appears to be still infalling towards the peak of the UC HII region, and its density and temperature are consistent with pressure confinement of the ionised gas to the SE.
Read the full abstract

All authors R. Cesaroni, M. T. Beltrán, V. M. Rivilla, Á. Sánchez-Monge, L. Colzi, F. Fontani, Á. López-Gallifa, A. Lorenzani, C. Mininni

Earth and Planetary Astrophysics

The Dispersed Matter Planet Project Sample – Detection limits, Occurrence Rates and New Planets

DMPP is a radial-velocity survey that aims to detect planets around stars exhibiting anomalous activity signatures, consistent with the presence of close-in evaporating planets. Here, we report the discovery of 7 new planetary signals in 5 different systems: DMPP-2c & d, HD67200/DMPP-6b & c, HD118006/DMPP-7b, HD191122/DMPP-8b, and HD200133/DMPP-9b. We update the orbital parameters of the DMPP-1, DMPP-2, and DMPP-3 systems, along with those of the planetary systems orbiting HD181433, HD39194, and HD89839. We derive detection limits for all 24 targets in our sample with adequate observational coverage, and test the DMPP hypothesis by calculating the occurrence rates for planets in this configuration. We find that the occurrence rates of planets in our sample with orbital periods shorter than $50~\mathrm{d}$ and masses in the range $3$-$10$ M$_\oplus$ are $83.0^{+27.1}_{-24.4}\%$, for $10$-$30$ M$_\oplus$ are $27.0^{+15.0}_{-11.2}\%$, and for $30$-$100$ M$_\oplus$ are $13.9^{+11.8}_{-7.5}\%$. This is significantly higher than the occurrence rates reported by other radial velocity surveys, providing strong support for the DMPP hypothesis.
Read the full abstract

All authors Matthew R. Standing, John R. Barnes, Carole A. Haswell, Adam T. Stevenson, João P. Faria, Erwan Quintin, Zachary O. B. Ross, Luca Fossati, James S. Jenkins, Douglas Alves, Daniel Staab

High Energy Astrophysical Phenomena

A geometric physics-informed machine learning inference for the neutron star maximum mass and the inverse problem

The existence of a distinct mass boundary between the heaviest neutron stars and the lightest black holes remains in question. It is an artefact of our ignorance of the properties of matter at supra-nuclear densities, which exist in the cores of neutron stars. The study addresses these problems with a physics-informed machine learning approach, guided by astrophysical observations. The Transformer model is trained on an agnostically generated ensemble of equations of state. Two geometric parameters are defined on the mass-radius sequence of a neutron star–the front bending and the back bending. The transformer provides a two-step solution: first, the model predicts the maximum mass and radius using the bending parameters. Second, it predicts the square of the sound speed profile, completing the inverse mapping. The prediction is that massive neutron stars form when the sound speed peaks at low density, leading to strong back-bending and an early phase transition to quark matter. Massive stars favour a stiff equation of state at low density, and the density of matter at the star's core is sufficiently small. The maximum mass for a neutron star predicted by the astrophysical constrained transformer model is $2.477$ solar masses, and a minimum radius of about $11.498$ km for a neutron star of $1.4$ solar masses.
Read the full abstract

Instrumentation and Methods for Astrophysics

Euclid: An automated system to match Rubin transient alerts to Euclid observations

The Vera C. Rubin observatory is expected to produce 10 million transient alerts per night in ugrizy filters, whilst Euclid is a visible to near-infrared space telescope engaged in a wide field survey. We present a prototype system to automatically match the transient alerts from Rubin to Euclid observations. The system produces joint light-curves containing both visible and near-infrared photometry, and joint image cutouts. Using Zwicky Transient Facility alerts as a proxy for Rubin, we demonstrate the system in use in cases where Euclid did and did not detect the transient and highlight the value that can be added in each case. For transients detected by Euclid these benefits include identifying the supernovae (SNe) in observations taken prior to ground-based detection, thereby better constraining the explosion time, such as SN 2024pvw detected ~3 d prior to ground based detections. In cases where Euclid did not detect the transient, we demonstrate the benefit of adding Euclid observations to improve host morphology measurements and associations.
Read the full abstract

All authors C. Duffy, I. M. Hook, C. M. Gutierrez, K. Paterson, V. Petrecca, T. J. Moriya, F. Poidevin, R. Kotak, B. Altieri, A. Amara, S. Andreon, N. Auricchio, C. Baccigalupi, M. Baldi, A. Balestra, S. Bardelli, P. Battaglia, A. Biviano, E. Branchini, M. Brescia, S. Camera, G. Cañas-Herrera, V. Capobianco, C. Carbone, J. Carretero, S. Casas, M. Castellano, G. Castignani, S. Cavuoti, K. C. Chambers, A. Cimatti, C. Colodro-Conde, G. Congedo, C. J. Conselice, L. Conversi, Y. Copin, F. Courbin, H. M. Courtois, M. Cropper, J.-C. Cuillandre, H. Degaudenzi, G. De Lucia, H. Dole, F. Dubath, X. Dupac, S. Dusini, S. Escoffier, M. Farina, R. Farinelli, S. Ferriol, M. Frailis, E. Franceschi, M. Fumana, S. Galeotta, K. George, W. Gillard, B. Gillis, C. Giocoli, J. Gracia-Carpio, A. Grazian, F. Grupp, S. V. H. Haugan, M. S. Holliman, W. Holmes, F. Hormuth, A. Hornstrup, K. Jahnke, M. Jhabvala, S. Kermiche, A. Kiessling, R. Kohley, B. Kubik, M. Kunz, H. Kurki-Suonio, A. M. C. Le Brun, S. Ligori, P. B. Lilje, V. Lindholm, I. Lloro, G. Mainetti, D. Maino, E. Maiorano, O. Mansutti, S. Marcin, O. Marggraf, M. Martinelli, N. Martinet, F. Marulli, R. J. Massey, E. Medinaceli, S. Mei, Y. Mellier, M. Meneghetti, E. Merlin, G. Meylan, A. Mora, M. Moresco, L. Moscardini, R. Nakajima, C. Neissner, S.-M. Niemi, C. Padilla, S. Paltani, F. Pasian, K. Pedersen, V. Pettorino, S. Pires, G. Polenta, M. Poncet, L. A. Popa, F. Raison, A. Renzi, J. Rhodes, G. Riccio, E. Romelli, M. Roncarelli, C. Rosset, R. Saglia, Z. Sakr, D. Sapone, B. Sartoris, M. Schirmer, P. Schneider, A. Secroun, G. Seidel, S. Serrano, E. Sihvola, P. Simon, C. Sirignano, G. Sirri, J. Skottfelt, L. Stanco, J. Steinwagner, P. Tallada-Crespí, A. N. Taylor, I. Tereno, N. Tessore, S. Toft, R. Toledo-Moreo, F. Torradeflot, I. Tutusaus, L. Valenziano, J. Valiviita, T. Vassallo, Y. Wang, J. Weller, G. Zamorani, E. Zucca, J. García-Bellido, E. Jullo, J. Martín-Fleitas, A. A. Nucita, V. Scottez

Solar and Stellar Astrophysics

Power-law Indices of EUV Intensity Power Spectrum in Flaring Coronal Active Regions

Solar intensity power spectra are usually characterised by coloured noise, with the spectral energy following a segmented power-law function of frequency, $S(f)\propto f^{-\alpha}$, over different frequency ranges. Typically, the power-law index exceeds 1 in the low-frequency part ($\alpha_\mathrm{lf}$) and is around 0 at high frequencies ($\alpha_\mathrm{hf}$). This work investigates the spatial and temporal evolution of the power-law indices of coronal EUV intensity power spectra in flare-hosting active regions. The spatial distribution of the power-law index in the low-frequency domain ($\alpha_\mathrm{lf}$) closely mirrors EUV intensity images, indicating that $\alpha_\mathrm{lf}$ can reveal the dynamics of coronal plasma structures. Temporally, $\alpha_\mathrm{lf}$ remains stable in quiescent active regions, but it exhibits significant variability before the flare onset. Motivated by this behaviour, we analysed 14 flare events, quantifying the temporal variation of the indices $\alpha_\mathrm{lf}$ and $\alpha_\mathrm{hf}$ as potential flare precursors. In all flare events considered, notable deviations of $\alpha_\mathrm{lf}$ beyond a defined threshold consistently occurred at the flare site within a few minutes before the flare. In some cases, the change in the value of $\alpha_\mathrm{lf} - \alpha_\mathrm{hf}$ was detected within 30–90\, minutes before the flare. This proof-of-concept study suggests that the temporal variation of the power-law indices in coronal EUV intensity power spectra could potentially serve as short-term precursors of solar flares, which needs to be validated on a larger flare sample.
Read the full abstract

Astrophysics of Galaxies

MIGHTEE HI observations of low surface brightness and ultra-diffuse galaxies in the XMM-LSS field

Untargeted neutral hydrogen (HI) surveys are well suited to identifying low surface brightness galaxies (LSBGs) that are gas rich, and they offer a complementary view to optically selected populations. We examined the LSBG population as identified via stellar and gaseous content using the MIGHTEE HI XMM-LSS early science data and the publicly available catalogs of optically identified LSBGs. There is currently little overlap between these datasets, with only three galaxies commonly detected. We performed surface brightness photometry of selected MIGHTEE HI detections to find 29 LSBGs, and 26 of these meet the size requirement (R_eff > 1.5 kpc) to be ultra-diffuse galaxies (UDGs). Furthermore, we extracted HI spectra at the location of all optically identified galaxies, placing upper limits on the HI-to-stellar mass ratio in these systems. While the HI-identified population overall tends toward bluer colors, the HI-identified and the optically selected samples mostly overlap in mean effective surface brightness, effective radii, and color. Although it is not straightforward to discern why the HI-identified LSBGs were missed in optical searches, this work highlights the utility of HI surveys in finding these faint systems. The HI-identified LSBGs are gas rich compared to the general HI-selected population. Furthermore, three out of four HI-selected UDGs with available kinematics show no systematic offset from the baryonic Tully-Fisher relation, although we are biased away from sources with low rotational velocities due to the low spectral resolution of the data. This work demonstrates the utility of HI observations for finding and characterizing the low surface brightness Universe.
Read the full abstract

All authors Elizabeth A. K. Adams, Barbara Šiljeg, Anastasia Ponomareva, Natasha Maddox, Pavel Enrique Mancera Piña, Maarten Baes, Bradley Frank, Marcin Glowacki, Matt Jarvis, Sambatriniaina Rajohnson, Gauri Sharma

Solar and Stellar Astrophysics

A Photometric and Spectroscopic investigation of 11 TESS eclipsing contact binaries

By cross-matching the eclipsing binary catalog provided by Prsa et al. (2022) with LAMOST medium resolution spectra, we obtained 11 targets. Combining light and radial velocity curves analysis, we have derived accurate physical parameters for these 11 targets. The results indicate that there are 3 deep contact binaries, 3 moderate ones, and 5 shallow ones. Among them, 3 targets exhibit the O'Connell effect, which is attributed to the presence of star-spot on the component's surface. One target is a low-mass ratio deep contact binary and may be contact binary merging candidates. The evolutionary status of these 11 targets was studied using the mass-luminosity and mass-radius relation diagrams. Based on the O-C (Observed minus Calculated) analysis of 10 targets, we found that the orbital periods of 5 contact binaries show a long-term decreasing trend, likely due to the combined effects of mass transfer between the two components and loss of angular momentum. Meanwhile, the orbital periods of the other 4 stars are continuously increasing, which is attributed to mass transfer. Besides, the O-C curves of 3 targets show clear periodic changes, which might result from the Applegate mechanism or the light travel time effect.
Read the full abstract

All authors Yani Guo, Kai Li, Yanke Tang, Xiang Gao, Qiqi Xia, Liheng Wang, Meng Guo

High Energy Astrophysical Phenomena

Updated Air-Shower $X_{\rm max}$ Moment Parametrizations for UHECR Composition with Latest Hadronic Interaction Models

The mass composition of ultra-high-energy cosmic rays (UHECRs) is commonly inferred from the first two moments of the depth of shower maximum, $X_{\rm max}$, measured by fluorescence and hybrid detectors. Such analyses require fast and accurate mappings between the moments of $X_{\rm max}$ and those of the logarithmic mass, $\ln A$, based on realistic air-shower simulations. In this work we provide updated parametrizations of the $X_{\rm max}$ moments and distributions for air showers initiated by nuclei from proton to iron, simulated with CONEX for three state-of-the-art hadronic interaction models: Epos LHC-R, Sibyll 2.3e, and QGSJet-III-01. We parametrize the mean depth $\langle X_{\rm max}\rangle$ and the variance $\sigma^2(X_{\rm max})$ as functions of energy and mass. For the variance we compare a second-order polynomial model with an exponential model. In addition, we model the full $X_{\rm max}$ distributions with a three-parameter generalized Gumbel function. The Gumbel parameters are fitted using an unbinned likelihood and are validated by comparing the implied mean and variance with the raw CONEX samples and with the moment parametrizations. Across the full energy range considered, residuals between the parametrizations (or the Gumbel representation) and the simulations are at the level of a few g cm$^{-2}$ for the mean and a few (g cm$^{-2}$)$^2$ for the variance, making these parametrizations suitable for precision UHECR composition studies and forward-folding analyses of $X_{\rm max}$ distributions.
Read the full abstract

All authors Carmelo Evoli, Igor Vaiman, Sergio Petrera, Francesco Salamida

Earth and Planetary Astrophysics

Testing the Icy Pebble Accretion Hypothesis with Primordial Main Belt Asteroids

Large main-belt asteroids (diameter $D \gtrsim 120\ \mathrm{km}$) exhibit a surface composition gradient as a function of heliocentric distance, ranging from anhydrous bodies to those rich in hydrated and, possibly, ammoniated materials. Their primordial nature holds key clues to the evolution of the Solar System. It has been suggested that the volatile-rich bodies formed in the outer Solar System and were implanted into the main belt. Alternatively, volatiles may have been delivered via inward-drifting icy pebbles in the protosolar disk. Here, we examine whether in-situ formed rocky embryos can acquire volatiles through pebble accretion as the snowline migrated inward. With the turbulence strength of the disk, radial pebble flux, and the dimensionless stopping time of pebbles (St) as parameters, we calculate the growth of large asteroids. The results are then compared with mass and compositional constraints based on asteroid observations. We find that a moderate pebble flux ($\lesssim18~M_\oplus / \text{Myr}$) is required to enable volatile delivery while preventing the largest asteroids from becoming more massive than Ceres. Water accretion is feasible with $\mathrm{St} \sim 10^{-3}$ ($\sim 1$ mm). However, only the largest asteroids (D > 200 km) can accumulate sufficient ammonia under such conditions. For most asteroids with D between 100 and 200 km, ammonia ice accretion requires $\mathrm{St} \sim 10^{-4}$ ($\sim 100\,\mu$m). Such small particle sizes may pose both theoretical and observational challenges. Thus, we propose that the intermediate-sized, potentially ammonia-bearing asteroids serve as a record of the Solar System's dynamic evolution.
Read the full abstract

Earth and Planetary Astrophysics

The GAPS Programme at the TNG: LXX. TOI-5734b: A hot sub-Neptune orbiting a relatively young K dwarf with an Earth-like density

Increasing interest in young exoplanets is leading to a growing effort to understand the formation and evolutionary processes responsible for their different architectures. One interesting target is TOI-5734, a relatively young K3-K4 dwarf star ($500_{-150}^{+300}$ Myr) showing a transiting candidate in photometric observations followed up with high-resolution spectroscopic data. Using Transiting Exoplanet Survey Satellite (TESS) photometry and High Accuracy Radial velocity Planet Searcher for the Northern hemisphere (HARPS-N) radial-velocity (RV) data, we aim to validate the presence of the companion TOI-5734b, measure its planetary mass, size, and its orbital parameters after having characterised its host star. We then aim to study its possible planetary composition and atmospheric evolution. By simultaneously modelling photometry and high-cadence RVs, we measured the radius, mass, and density of TOI-5734b precisely. In particular, we employed Gaussian processes (GPs) with a flexible kernel to discriminate between the stellar activity of the young host and planetary signals. We confirmed the planetary nature of TOI-5734b and measured its orbital period ($P_{\rm b}\sim6.18$ d), radius ($R_{\rm b} = 2.10^{+0.12}_{-0.12}$ $R_\oplus$), and mass ($M_{\rm b}=9.1^{+2.6}_{-2.6}$ $M_\oplus$). By measuring its density ($\rho_{\rm b}=0.98_{-0.30}^{+0.36}$ $\rho_\oplus$), we infer that TOI-5734b is close to having a rocky composition and an almost completely depleted primary envelope. Our results point toward the possibility of considering the target for atmospheric studies with present and future ground- and space-based facilities.
Read the full abstract

All authors S. Filomeno, T. Trifonov, M. Damasso, M. Baratella, S. Benatti, K. Biazzo, K.A. Collins, R. Cosentino, S. Desidera, C. Di Maio, D. Locci, A. Maggio, L. Mancini, S. Messina, L. Naponiello, D. Nardiello, K. G. Stassun, F. Amadori, S. Antoniucci, F. Biassoni, A.S. Bonomo, L. Cabona, C.A. Clark, M. Gonzalez, A. F. Lanza, D.W. Latham, V. Lorenzi, L. Malavolta, G. Mantovan, M. Pedani, G. Scandariato, A. Shporer, A. Sozzetti, T. Zingales, D. Ciardi, M. Everett, M.B. Lund, H. Osborn, I. A. Strakhov, J. Villaseñor, D. Watanabe, A. Youngblood, R. Zambelli

High Energy Astrophysical Phenomena

TROYE: Modeling Dynamic Phase Transitions in Gravitational Waves from Neutron Star-Black Hole Mergers

The Equation of State (EoS) of dense nuclear matter remains one of the most compelling open questions in high-energy astrophysics. While static EoS models are increasingly well-constrained by observations of binary neutron star (BNS) inspirals, the possibility of a dynamic phase transition occurring during the coalescence has been thus far deferred from standard gravitational-wave (GW) analyses. In this work, we investigate the detectability of such a phase transition, manifesting as a macroscopic shift in the tidal deformability parameter $\Lambda$, using GWs from Neutron Star-Black Hole (NSBH) coalescences. We argue that NSBH systems serve as a cleaner laboratory for this phenomenology than BNS systems due to the absence of the $\tilde{\Lambda}(\Lambda_1,\Lambda_2)$ degeneracy, allowing for the isolation of single-body tidal evolution. We introduce a phenomenological waveform model, TROYE (Transitional Representation Of varYing Equation-of-state), which stitches together two waveform approximants in the time domain to simulate a smooth but rapid transition between two equations of state during the late inspiral. We perform a comprehensive Bayesian injection and recovery campaign on 100 simulated events using the bilby inference library. Our results demonstrate that a phase transition corresponding to a tidal shift of $|\Delta\Lambda| \gtrsim 400$ is detectable with Advanced LIGO design sensitivity, yielding decisive statistical evidence ($\ln B > 5$). We further identify a "V-shape" asymmetry in detectability, where "softening" transitions (decreasing $\Lambda$) are systematically easier to detect than "stiffening" ones due to the specific phase evolution of the tidal sector. Finally, we present "stress tests" showing that the transition remains recoverable even when marginalized over uncertainties in the stitching time and binary mass ratio.
Read the full abstract

Earth and Planetary Astrophysics

Inhomogeneous magnetic coupling in exoplanets: the stop & go of WASP-18 b's atmospheric flows

Early studies of ionization in hot Jupiter atmospheres suggest that magnetic coupling can shape their dynamics. These effects may be most pronounced in ultra-hot Jupiters that sustain global magnetic fields. WASP-18 b hosts an ionized dayside atmosphere extending deep enough to be strongly influenced by magnetic forces. Phase curve observations suggest effective magnetic drag, yet its impact on the atmospheric circulation remains poorly constrained. This work explores how magnetic drag in an inhomogeneously ionized atmosphere shapes local and global dynamics to provide a pathway to constrain the planet's magnetic field strength. An analytical parameterization for anisotropic magnetic drag, including both Pedersen and Hall drag components, and associated frictional heating in the globally neutral atmosphere, is implemented in the 3D General Circulation Model ExoRad to study WASP-18 b's atmosphere. Climate characteristics are compared for different drag formulations to assess whether anisotropic physics is required to capture magnetic coupling effects. Anisotropic magnetic drag and frictional heating, both set by local ionization, strongly affect wind strength and direction in the upper atmosphere, modify the day-night circulation, and produce observable temperature asymmetries. They enhance the evening-morning terminator temperature difference near 0.1 bar and generate two off-equator hotspots with reduced eastward shift. The terminator regions are particularly sensitive to how magnetic drag is modeled. Anisotropic magnetic drag damps and redirects dayside-to-nightside winds, partially decoupling the equatorial flow at the morning terminator while maintaining the nightside jet. Locally varying drag forces and frictional heating create asymmetric temperature patterns manifesting as primary and secondary hotspot regions.
Read the full abstract

Solar and Stellar Astrophysics

3D radiative transfer modeling of scattering polarization with partial frequency redistribution I. Verification and disk-center results for the solar Ca I 4227 Å line

Several strong solar resonance lines show observable linear scattering polarization signals, holding a great potential for investigating the magnetism of the outer solar atmosphere. Accurately modeling these signals requires solving the radiative transfer (RT) problem for polarized radiation in comprehensive 3D models of the solar atmosphere, in non-local thermodynamic equilibrium, accounting for partial frequency redistribution (PRD) effects. This problem has so far been computationally inaccessible. We present the first scientific application of TRIP, a novel software for the massively parallel solution of the 3D non-LTE RT problem for polarized radiation, including scattering polarization and PRD. We aim to verify the code and explore the combined action of PRD and the 3D structure of the solar atmosphere on scattering polarization. We run TRIP to synthesize the Stokes profiles of the Ca I line at 4227 Å in a 3D model of the solar atmosphere extracted from a radiation magneto-hydrodynamic simulation. We efficiently solve the resulting large-scale problem, with up to $4 \times 10^{10}$ degrees of freedom, with a state-of-the-art preconditioned Krylov method, using up to 20 thousand parallel CPUs. After including verification tests, we find that the joint impact of PRD effects and the detailed 3D structure of the atmospheric model produce disk-center scattering polarization signals in the line wings. These signals are sensitive to the magnetic field, via magneto-optical effects, and to bulk velocity gradients. We also show that the CRD approximation underestimates the amplitude of disk-center line-core signals. This achievement represents a crucial step forward for diagnosing the magnetism of the solar chromosphere and transition region through the quantitative comparisons of synthetic and observational data.
Read the full abstract

All authors Pietro Benedusi, Simone Riva, Tanausú del Pino Alemán, Gioele Janett, Fabio Riva, Jirí Štepán, Rolf Krause, Javier Trujillo Bueno, Luca Belluzzi

Solar and Stellar Astrophysics

The attempted polarity reversal and evolving magnetic environment of AD Leo

In the past two decades, the observed large-scale magnetic field of the active M dwarf star AD Leo has evolved from strongly to mildly negative, raising a suspicion that it might switch polarity. Although magnetic field reversals are observed every 11 years for the Sun, such reversals are poorly understood for M dwarfs. Further, no reversals have been observed for fast-rotating M dwarfs. We examine the properties of AD Leo's large-scale magnetic field and investigate how its evolution affects the space weather environment. We analysed spectropolarimetric data collected by ESPaDOnS and SPIRou in late-2022 and early-2023. With the optical and near-infrared data we computed the longitudinal magnetic field, and with the near-infrared data reconstructed the large-scale magnetic field using Zeeman-Doppler imaging. Using five magnetograms, from 2019 to 2023, we simulated three-dimensional Alfven wave-driven stellar winds using the space weather code SWMF. Although we see an evolution of the large-scale magnetic field of AD Leo, we find no polarity reversal. Rather, we see a restoration of the field to a simpler configuration with consistently negative values for the longitudinal magnetic field strength. Our new large-scale field reconstruction for AD Leo is characterised by a highly axisymmetric, poloidal-dipolar field with an increased mean large-scale field strength. SWMF simulations find the stellar mass loss rates to be, on average, an order of magnitude greater than that of the Sun. Additionally, we find that the habitable zone resides beyond the Alfven surface. Hypothetical magnetised habitable zone planets (with planetary field strengths greater than 0.34 G) would likely be shielded from the incident wind and atmospheric erosion would be negligible. Further, we find variable conditions across each epoch due to the evolving axisymmetry of the stellar large-scale magnetic field.
Read the full abstract

All authors K.G. Smith, D. Evensberget, S. Bellotti, J. Morin, A.A. Vidotto, B.J.S. Pope

Solar and Stellar Astrophysics

Modeling brightness temperature of sunspots using ALMA single-dish observations

We model sunspot brightness temperature in the wavelength range 0.3-10mm using the Atacama Large Millimeter/submillimeter Array (ALMA) single-dish measurements together with other instruments. Our modeling idea starts from an existing model, the 1D semi-empirical Avrett-Tian-Landi-Curdt-Wülser (ATLCW) model for sunspots. The original model is then modified in density and temperature with a constant factor over all heights. A least-square fit in both density and temperature to the observed values yields an improved model for sunspot brightness temperature. Observations show that sunspots are darker than QS at wavelengths below 2 mm and brighter beyond, while fitting of the ATLCW sunspot model imply that sunspots should be darker than QS over the entireALMArange. The best-fitted sunspot brightness temperature profiles resulted in 3.6-4.0 times lower sunspot density compared to the original sunspot ATLCW model, or 1.4-1.8 times higher density compared to QS. The temperature came out around 1.2 times higher than the original ATLCW predictions for sunspots, or 1.2 times lower than QS temperature. Even with improved density and temperature parameters, the remaining discrepancy between the bestfitted models and observations indicated that additional improvements are needed in the ATLCW model and modeling assumptions.
Read the full abstract

All authors F. Matković, R. Brajša, A. O. Benz, H.-G. Ludwig, C. L. Selhorst, I. Skokić, D. Sudar, A. Hanslmeier

Astrophysics of Galaxies

Investigating the role of mergers in galaxy assembly in the early Universe (z > 5)

Galaxy mergers play a crucial role in shaping the morphology, the star formation, and the mass growth of galaxies across cosmic time. While mergers have been extensively investigated in the local Universe, the evolution of their frequency and physical properties in the early Universe has yet to be fully understood. We investigate the role of mergers in a large spectroscopic sample of 1233 galaxies in the range 5<z<14 with good detection (S/N-pixel > 3) in JWST imaging, covering six different extragalactic fields. We identify mergers from rest-frame optical disturbances in F444W, using a combination of Gini, M-20, and Asymmetry parameters. We find a morphological merger fraction f_m that does not strongly evolve with redshift from z=0 to z ~ 8. The average f_m of our primary major merger condition (Gini+0.14xM-20 > 0.33, A>0.35) is ~ 5 %, which increases to ~13 % for major+minor merger tracers. Accounting for the evolving observability timescale of each tracer, we find that the merger rate is strongly increasing from z=1 to 7 by more than 1 dex, averaging ~ 2 merger/galaxy/Gyr at 5<z<10 for major mergers (in agreement with photometric pair studies), and a factor of 3 higher for minor+major mergers. We also perform SED modeling using available HST+JWST photometry to infer stellar masses and SFRs, using a non parametric star-formation history. We find that mergers at z > 5 have a significant impact, although significantly lower than at z<1, on the SFR of galaxies. When averaged over 10 Myr (comparable to the observability timescale of morphological disturbances), their SFRs are a factor of 1.7 higher than a mass and redshift matched sample of non-mergers, suggesting that mergers trigger new star-formation through short-lived powerful bursty episodes. Despite this, mergers contribute only by 5% - 10% to the mass build-up of galaxies in the redshift range explored.
Read the full abstract

All authors A. Calabrò, L. Pentericci, M. Llerena, S. Rossi, L. Napolitano, D. Bevacqua, M. Giavalisco, R. Somerville, G. Gandolfi, E. Daddi, M. Dickinson, S. Finkelstein, A. Fontana, M. Hirschmann, J. S. Kartaltepe, D. Kocevski, A. Koekemoer, H. Leung, R.A. Lucas, A. Taylor, R. Tripodi, X. Wang, L.Y.A. Yung

Earth and Planetary Astrophysics

Probabilistic Methods for Initial Orbit Determination and Orbit Determination in Cislunar Space

In orbital mechanics, Gauss's method for orbit determination (OD) is a popular, minimal assumption solution for obtaining the initial state estimate of a passing resident space object (RSO). Since much of the cislunar domain relies on three-body dynamics, a key assumption of Gauss's method is rendered incompatible, creating a need for a new, minimal assumption method for initial orbit determination (IOD). In this work, we present a framework for short and long term probabilistic target tracking in cislunar space which produces an initial state estimate with as few assumptions as possible. Specifically, we propose an IOD method involving the kinematic fitting of several series of noisy, consecutive ground-based observations. Once a probabilistic initial state estimate in the form of a particle cloud is formed, we apply the powerful Particle Gaussian Mixture (PGM) Filter to reduce the uncertainty of our state estimate over time. This combined IOD/OD framework is demonstrated for several classes of trajectories in cislunar space and compared to better-known filtering frameworks.
Read the full abstract

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

Dark Photon mediated Inelastic Dark Matter in Cosmology, Astrophysics and Colliders

We explore the phenomenology of Dark Photon iDM ($A^\prime$IDM) where the Standard Model (SM) is extended by a dark sector containing an additional $U(1)_D$ gauge symmetry with coupling $\alpha_D$ under which all SM particles are neutral, and that couples to the SM hypercharge gauge boson through a kinetic mixing parameter $\epsilon$. The model contains two Majorana states $\chi_1$ and $\chi_2$ with $\delta=M_{\chi_2}-M_{\chi_1}>0$ and $\chi_1$ the dark matter candidate, and a dark photon $A^{\prime}$ with mass $M_{A^{\prime}}$. Our analysis represents an integration of existing ones, where only specific benchmarks of the A'IDM scenario have been discussed. We perform a complete scan of the parameters of the model ($M_{\chi_1}$, $\delta$, $M_{A^{\prime}}$, $\epsilon$, $\alpha_D$) and discuss the $\chi_1$ relic abundance, its direct and indirect searches, as well as potential signals from astrophysics and accelerators. Our scan shows that $\alpha_D=\alpha_{EM}$ is not disfavored, as some previous analyses, limited to specific benchmarks, may suggest. We also find that when the $\chi_1$ relic density matches observation direct and indirect searches are not kinematically accessible. On the other hand we find that the projected luminosity of FASER, a detector searching for Long Lived Particles (LLP) decay at the LHC, can probe or rule out the parameters space of the model for $M_{\chi_1}\lesssim$ 7 GeV, 100 MeV$\lesssim\delta\lesssim$300 MeV and $M_{A^\prime}\lesssim$25 GeV. This range of parameter could be significantly extended by the FASER 2 upgrade proposed for the High-Luminosity phase at the LHC. The parameter space probed by LLP seaches partially overlaps with that probed by $\chi_1$ capture in neutron stars.
Read the full abstract

Astrophysics of Galaxies

MATLAS-42, A Globular Cluster-Rich Ultra-Diffuse Galaxy That Diverges from the "Failed Galaxy” Formation Pathway

To date, there has been significant interest in globular cluster (GC)-rich ultra-diffuse galaxies (UDGs) and the evidence that they have formed via an unexpected, “failed galaxy” formation pathway. The majority of the evidence for “failed galaxy” UDGs originates from spectroscopic observations targeting passive GC-rich UDGs, with a focus on those residing in galaxy clusters. In this work, we study the gas-rich, GC-rich group UDG MATLAS-42 and derive its stellar population properties using the Keck Cosmic Web Imager. We measure a redshift for the galaxy ($V_{\rm R, \star}=2433\pm8$~km s$^{-1}$), confirming the previous assumptions that it is both part of the NGC~502 group and has an associated HI-reservoir ($V_{\rm R,HI}=2423\pm 15$~km s$^{-1}$). We measure integrated stellar populations and find the galaxy to be both young (mass-weighted age $=3.2^{+2.6}_{-1.5}$Gyr) and of average-to-low metallicity ($[M/H]=-1.19^{+0.42}_{-0.30}$ dex). When considering these properties in the context of the galaxy's formation, we note it likely does not follow the “failed galaxy” formation pathway commonly attributed to GC-rich, cluster UDGs, as it has experienced recent star formation. At most it started failed, however, it has recently rejuvenated its star formation. Finally, we build a toy model of the passive evolution of this galaxy, finding that its relative GC-richness (i.e., $M_{\rm GC}/M_\star$) will likely decrease with time as GCs slowly evaporate/disrupt to contribute to the stellar mass of the galaxy. Due to this, we hypothesise that it is likely not a low redshift analogue of the progenitor to a “failed galaxy” UDGs.
Read the full abstract

All authors Jonah S. Gannon, Duncan A. Forbes, Francine R. Marleau, Anna Ferré-Mateu, Aaron J. Romanowsky, Maria Luisa Buzzo, Jean P. Brodie

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

Decay of a multi-axionic SU(N) symmetric color aether in the early Universe as an origin of emergence of a many-component dark matter

We establish a new SU(N) symmetric model of interaction between the fields of four types: the multiplet of vector fields, which describes the so-called color aether, the multiplet of pseudoscalar fields, which is associated with the multi-component cosmic dark matter, the gauge and gravitational fields. The extended Lagrangian of the model contains a new constructive element, which is based on the covariant SU(N) symmetric divergence of the multiplet of vector fields; this new element, being the multiplet of scalars from the point of view of spacetime transformations and the color vector from the point of view of the SU(N) group space, gives us the possibility to formulate properly the multi-axionic extension of the Peccei-Quinn theory. The hypothesis of a spontaneous polarization of the multi-axionic color aether in the early Universe is presented. The set of self-consistent master equations of the model is derived. An application to cosmology is considered: the obtained master equations are solved for the truncated test model based on the Bianchi-I spacetime platform.
Read the full abstract

Astrophysics of Galaxies

A Comprehensive Catalog of Radio Sources and Rotation Measures in the Perseus Molecular Cloud from Very Large Array Observations

We present a comprehensive radio polarization study of the Perseus molecular cloud using wideband L-band observations from the Karl G. Jansky Very Large Array. Our survey covers $\sim13.8$~deg$^2$ with a mean Stokes~$I$ sensitivity of $\sim80~\mu$Jy~beam$^{-1}$, enabling the detection of 1410 compact radio sources. From this population, we construct a catalog of source properties, including positions, integrated flux densities, and spectral indices measured across nine spectral windows. The majority of sources exhibit negative spectral indices, consistent with non-thermal synchrotron emission. Using RM Synthesis and RM CLEAN techniques, we detect 205 polarized background sources above an $8\sigma$ threshold. This corresponds to a sampling density of $\sim14.8$~deg$^{-2}$, representing more than an order-of-magnitude increase compared to previous NVSS-based measurements. The resulting rotation measures exhibit coherent large-scale variations across the surveyed region, with additional small-scale structure superimposed. The enhanced sensitivity, frequency coverage, and sampling density of our observations enable a substantially improved mapping of the line-of-sight magnetic field component toward the Perseus molecular cloud compared to previous surveys.
Read the full abstract

All authors Haleh Hajizadeh, Jeroen Stil, René Plume, Mehrnoosh Tahani, Preshanth Jagannathan

Astrophysics of Galaxies

A Bayesian Exploration of the Mass of Ursa Major III: Kinematics, Rotation and their influence on the Mass to Light Ratio

We investigate the kinematics of the potential ultra-faint dwarf galaxy (UFD) UMa III/U1 using Bayesian inference to search for the signal of any potential intrinsic rotation. The magnitude of rotation is relevant to estimating the total mass of UMa III/U1, which is critical in determining whether or not UMa III/U1 is in fact a UFD, or possibly a star cluster home to a significant binary fraction. A non-rotating model and a rotational model are fitted for the current total population of member stars of UMa III/U1, finding that a non-rotating model was preferred by a factor of $\sim 5-12 \times$. This was repeated on a reduced population of UMa III/U1, where potential contaminant stars were removed. A similar preference for non-rotation was found for these reduced populations. We calculate a lower-bound rotational mass estimate for UMa III/U1 and a corresponding lower bound mass-to-light ratio of $ 734.4^{+339.0}_{-176.2} \mathrm{M_\odot} / \mathrm{L_\odot} $ for the total population. We conclude that UMa III/U1 still remains an ambiguous object with viable arguments for both the UFD and self-gravitating star cluster scenarios, however under both, UMa III/U1 is unlikely to be supported by rotational pressure.
Read the full abstract

Instrumentation and Methods for Astrophysics

dfcosmic: A Python package for cosmic ray removal

Astronomical images often show sharp features that are caused by cosmic ray (CR) hits, hot pixels, or non-Gaussian noise. L.A.Cosmic (van Dokkum 2001) is a widely used edge detection algorithm that identifies and replaces such features. Here we describe dfcosmic, a direct python port of L.A.Cosmic utilizing PyTorch and C++ to enable efficient performance on both CPUs and GPUs. The code was developed for the MOTHRA array, which is projected to produce more than 1000 large format CMOS images every 15 minutes. Compared to previous python implementations, dfcosmic achieves a speed gain of at least 20%.
Read the full abstract

All authors Carter Lee Rhea, Pieter van Dokkum, Steven R. Janssens, Imad Pasha, Roberto Abraham, William P. Bowman, Deborah Lokhorst, Seery Chen