New articles on Astrophysics


[1] 2507.08048

A candidate field for deep imaging of the Epoch of Reionization observed with MWA

Deep imaging of structures from the Cosmic Dawn (CD) and the Epoch of Reionization (EoR) in five targeted fields is one of the highest priority scientific objectives for the Square Kilometre Array (SKA). Selecting 'quiet' fields, which allow deep imaging, is critical for future SKA CD/EoR observations. Pre-observations using existing radio facilities will help estimate the computational capabilities required for optimal data quality and refine data reduction techniques. In this study, we utilize data from the Murchison Widefield Array (MWA) Phase II extended array for a selected field to study the properties of foregrounds. We conduct deep imaging across two frequency bands: 72-103 MHz and 200-231 MHz. We identify up to 2,576 radio sources within a 5-degree radius of the image center (at RA (J2000) $8^h$ , Dec (J2000) 5°), achieving approximately 80% completeness at 7.7 mJy and 90% at 10.4 mJy for 216 MHz, with a total integration time of 4.43 hours and an average RMS of 1.80 mJy. Additionally, we apply a foreground removal algorithm using Principal Component Analysis (PCA) and calculate the angular power spectra of the residual images. Our results indicate that nearly all resolved radio sources can be successfully removed using PCA, leading to a reduction in foreground power. However, the angular power spectra of the residual map remains over an order of magnitude higher than the theoretically predicted CD/EoR 21 cm signal. Further improvements in data reduction and foreground subtraction techniques will be necessary to enhance these results.


[2] 2507.08056

DESI and fast radio burst used to constrain modified theories of gravity

This paper is devoted to study the viability of DESI and fast radio burst to constraint free parameters of modified theories of gravity. Thus, we present a model supported in $f(R)$ gravity involving a function of the Ricci scalar named Starobinsky-type with the peculiarity that the Non-commutative essence is intrinsic to the coefficients. Additionally, to understand the dynamics within a flat Friedmann-Lemaître-Robertson-Walker Universe, we explore the possibility of deriving a Friedman equation (in measure) that results from an adequate mathematical treatment. As we mention previously, to test the outlined model, a Monte Carlo Markov Chain analysis is implemented, using Cosmic Chronometers, type Ia supernovae, Hydrogen II galaxies, Intermediate-luminosity quasars, Baryon Acoustic oscillations and Fast Radio Bursts data, to constraint the free parameters of the model and presenting $H(z)$, $q(z)$ and $\omega_{eff}(z)$. The final results are compared with the $\Lambda$CDM model and it is presented a robust discussion about the viability of DESI and fast radio burst to constraint free parameters in specific of a modified theory of gravity.


[3] 2507.08070

TAMBO: A Deep-Valley Neutrino Observatory

Although the field of neutrino astronomy has blossomed in the last decade, physicists have struggled to fully map the high-energy neutrino sky. TAMBO, a mountain-based neutrino observatory, aims to solve that issue -- and find clues of new physics along the way.


[4] 2507.08074

Exploring substructures in the Milky Way halo Neural networks applied to Gaia and APOGEE DR 17

The identification of stellar structures in the Galactic halo, including stellar streams and merger remnants, often relies on the dynamics of their constituent stars. However, this approach has limitations due to the complex dynamical interactions between these structures and their environment. Perturbations such as tidal forces exerted by the Milky Way, the potential escape of stars, and passages through the Galactic plane can result in the loss of dynamical coherence of stars in these structures. Consequently, relying solely on dynamics may be insufficient for detecting such disrupted or dispersed remnants. We combine chemistry and dynamics, integrated through a system of neural networks, to develop a clustering method for identifying accreted structures in the Galactic halo. We developed an integrated approach combining Siamese neural networks (SNNs), graph neural networks (GNNs), autoencoders, and the OPTICS algorithm to create a comprehensive procedure named CREEK. This method is designed to uncover stellar structures in the Galactic halo. Initially, CREEK was trained on known globular clusters (GCs) and then applied to the dataset to identify stellar streams. CREEK successfully recovered 80% of the GCs present in the APOGEE dataset, re-identified several known stellar streams, and identified a potential new stream. Additionally, within highly populated stellar structures, CREEK can identify substructures that exhibit distinct chemical compositions and orbital energies. This approach provides an objective data-driven method for selecting stars associated with streams and stellar structures in general.


[5] 2507.08075

Fluctuations in Hill's equation parameters and application to cosmic reheating

Cosmic inflation provides a compelling framework for explaining several observed features of our Universe, but its viability depends on an efficient reheating phase that converts the inflaton's energy into Standard Model particles. This conversion often proceeds through non-perturbative mechanisms such as parametric resonance, which is described by Hill's equation. In this work, we investigate how stochastic fluctuations in the parameters of Hill's equation can influence particle production during reheating. We show that such fluctuations can arise from couplings to light scalar fields, and can significantly alter the stability bands in the resonance structure, thereby enhancing the growth of fluctuations and broadening the region of efficient energy transfer. Using random matrix theory and stochastic differential equations, we decompose the particle growth rate into deterministic and noise-induced components and demonstrate analytically and numerically that even modest noise leads to substantial particle production in otherwise stable regimes. These results suggest that stochastic effects can robustly enhance the efficacy of reheating across a wide swath of parameter space, with implications for early Universe cosmology, UV completions involving multiple scalar fields, and the resolution of the cosmological moduli problem.


[6] 2507.08077

MIRACLE II: Unveiling the multi-phase gas interplay in the circumnuclear region of NGC 1365 via multi-cloud modeling

We present a multi-phase study of the gas in the circumnuclear region (~1.1x1.0 kpc^2) of the nearby Seyfert 1.8 galaxy NGC 1365, observed in the context of the Mid-IR Activity of Circumnuclear Line Emission (MIRACLE) program. We combined spatially resolved spectroscopic observations from JWST/MIRI, VLT/MUSE, and ALMA to investigate the ionized atomic gas and the warm and cold molecular phases. MIRI data revealed over 40 mid-IR emission lines from ionized and warm molecular gas. Moment maps show that both cold and warm molecular gas follow the rotation of the stellar disk along the circumnuclear ring. The ionized gas displays flux and kinematic patterns that depend on ionization potential (IP): low-IP species (<25 eV) trace the disk, while higher-IP lines (up to ~120 eV) trace outflowing material. The [O III]5700 and [Ne V]14 lines both trace the southeast nuclear outflow cone. Additionally, [Ne V]14 detects the northwest counter-cone, obscured in the optical and thus invisible in [O III]5700. Mid-IR diagnostics, unlike optical ones, clearly reveal the AGN as the primary ionization source in the nucleus. Emission from high-IP species is spatially coincident with the ionization cones and not with star-forming regions. Using the [Ne V]24/[Ne V]14 ratio, we derive an electron density of (750+-440) cm^(-3), in agreement with values from the [S II] optical doublet. For the first time, we apply a fully self-consistent approach combining advanced photoionization and kinematic models (HOMERUN+MOKA3D) to constrain intrinsic outflow properties, overcoming the limitations of simplified classical methods. Exploiting the synergy of JWST/MIRI and VLT/MUSE, HOMERUN reproduces fluxes of over 60 emission lines from optical to mid-IR, disentangling AGN and star formation contributions and yielding robust estimates of outflow mass, geometry, and energetics.


[7] 2507.08078

SN 2023ixf: The Closest Supernova of the Decade

Supernova 2023ixf occurred on May 18th 2023 in the nearby galaxy Messier 101 (D ~ 6.85 Mpc), making it the closest supernova in the last decade. Following its discovery, astronomers around the world rushed to observe the explosion across the electromagnetic spectrum in order to uncover its early-time properties. Based on multi-wavelength analysis during its first year post-explosion, supernova 2023ixf is a type II supernova that interacted with dense, confined circumstellar material in its local environment -- this material being lost from its red supergiant progenitor in the final years before explosion. In this article, we will review the findings of >80 studies already published on this incredible event as well as explore how the synthesis of SN 2023ixf observations across the electromagnetic spectrum can be used to constrain both type II supernova explosion physics in addition to the the uncertain mass loss histories of red supergiant stars in their final years.


[8] 2507.08079

Dissecting NGC 3132: Tracing the mass-loss history of the southern ring planetary nebula

Central to our understanding of stellar evolution and its impact on processes in our Galaxy and across the Universe is the study of mass loss. While the general framework is well established, recent JWST observations of objects like NGC 3132 have revealed intricate nebular structures, suggesting complex mass-loss processes likely driven by multiple star system at its core. These findings pose new challenges for the currently available investigation tools. The primary goal of this study is the first detailed comparison of the physical properties and chemical composition obtained for NGC 3132, based on the latest detailed 3D model and observations from MUSE, JWST and Spitzer. We evaluate the reliability of the traditional empirical method and photoionization model for abundances estimations, both based on the same available high-quality, spatially resolved observations. We find that the model and empirical method yield consistent results for the integrated total properties such as Te, ne and chemical abundances. However, when applied to simulated observations from the model, the empirical method fails to recover the model input abundances, providing only an approximate estimate. This discrepancy arises in part from the loss of information when summing fluxes over regions which have complex ionisation structures. This discrepancy in the case of oxygen has been estimated to be up to 35%. Moreover, the latest IR data reveal a spatial correlation between H2, c(Hb) as well as the [8.0]/[4.5] IRAC ratio. Finally, new clumps are discovered in [Ni II] 7378 Å, [Fe II] 8617 Åand [Fe III] 5270 Åemission lines.


[9] 2507.08082

Tidal Disruption Event Demographics in Supermassive Black Hole Binaries Over Cosmic Times

Tidal disruption events (TDEs) offer a unique probe of supermassive black hole (SMBH) demographics, but their observed rates remain difficult to reconcile with standard single-SMBH models. In this work, we use simulations of SMBH binaries, including the combined effects of eccentric Kozai-Lidov oscillations and two-body relaxation, to explore how TDE rates scale with SMBH mass and redshift. We find that binary systems exhibit increasing TDE rates with mass, in contrast to the declining trend expected for single SMBHs. These binary-driven rates match those observed in post-starburst galaxies, suggesting that a subset of TDE hosts may contain SMBH binaries. TDE light curves in some massive galaxies exhibit unexpectedly short durations, suggesting that the disrupting SMBH may be less massive than implied by host galaxy scaling relations, consistent with disruptions by the less massive black hole in a binary. By convolving our mass-dependent rates with the SMBH mass function, we predict redshift-dependent TDE rates, which we show can be used to constrain the supermassive black hole binary fraction. Our results provide a testable framework for interpreting TDE demographics in upcoming wide-field surveys such as LSST and Roman.


[10] 2507.08084

Reverberation Mapping Data of NGC 5548 Imply a Multicomponent Broad-Line Region

Line broadening and variability are observational hallmarks of active galactic nuclei which allow us to measure supermassive black hole masses as well as constrain the geometry and kinematics of the emitting gas, with the most precise measurements requiring a degree of modeling. Two popular models of the broad-line region describe the emitting gas as either a distribution of puffed up clouds or a thin disk with strong velocity gradients. As we show in this work, key features in the reverberation mapping dataset obtained by the AGN STORM team in NGC 5548 cannot be accounted for by either simple model. In several emission lines the observed broad-line profile has a single peak yet the delay profile has a distinct double peak, strongly motivating a BLR with emission from multiple components. We demonstrate a few possibilities that may alleviate the tension and better represent the true nature of the broad-line emitting gas in NGC 5548 and beyond.


[11] 2507.08093

Observation and Modeling of Small Spatial Structures of Solar Radio Noise Storms using the uGMRT

One of the most commonly observed solar radio sources in the metric and decametric wavelengths is the solar noise storm. These are generally associated with active regions and are believed to be powered by the plasma emission mechanism. Since plasma emission emits primarily at the fundamental and harmonic of the local plasma frequency, it is significantly affected by density inhomogeneities in the solar corona. The source can become significantly scatter-broadened due to the multi-path propagation caused by refraction from the density inhomogeneities. Past observational and theoretical estimates suggest some minimum observable source size in the solar corona. The details of this limit, however, depends on the modeling approach and details of the coronal turbulence model chosen. Hence pushing the minimum observable source size to smaller values can help constrain the plasma environment of the observed sources. In this work, we for the first time, use data from the upgraded Giant Metrewave Radio Telescope in the 250--500 MHz band, to determine multiple instances of very small-scale structures in the noise storms. We also find that these structures are stable over timescales of 15--30 minutes. By comparing the past observations of Type III radio bursts and noise storms, we hypothesize that the primary reason behind the detection of these small sources in noise storm is due to the local environment of the noise storm. We also build an illustrative model and propose some conditions under which the minimum observable source size predicted by theoretical models, can be lowered significantly.


[12] 2507.08097

Measuring the Neutrino Flux in Segments along the Galactic Plane with IceCube

Gamma-ray emission from the plane of the Milky Way is understood as partly originating from the interaction of cosmic rays with the interstellar medium. The same interaction is expected to produce a corresponding flux of neutrinos. In 2023, IceCube reported the first observation of this galactic neutrino flux at 4.5$\sigma$ confidence level. The analysis relied on neutrino flux predictions - based on gamma ray observations - to model the expected neutrino emission from the galactic plane. Three signal hypotheses describing different possible spatial and energy distributions were tested, where the single free parameter in each test was the normalization of the neutrino flux. We present first results of an analysis that can improve the characterization of Galactic neutrino emission by dividing the galactic plane into segments in galactic longitude. An unbinned maximum-likelihood analysis is used that can fit the spectral index and the flux normalization separately in each segment. While gamma ray telescopes can not differentiate between hadronic and leptonic emission, neutrino production must come from hadronic processes. Measuring a spectral index can further help to understand the contribution of unresolved neutrino sources inside the galactic plane. This work uses a full-sky cascade dataset and provides model-independent insight into the variation of the neutrino flux and energy distribution from different regions of the galactic plane.


[13] 2507.08098

JWST Discovery of Warm Dust in the Circumgalactic Medium of the Makani Galaxy

We report the detection of near- and mid-infrared emission from polycyclic aromatic hydrocarbons (PAHs) out to ~ 35 kpc in the Makani Galaxy, a compact massive galaxy with a record-breaking 100-kpc scale starburst-driven wind at redshift z = 0.459. The NIRCam and MIRI observations with JWST take advantage of a coincidental match between the PAH spectral features at 3.3, 7.7, and (11.3 + 12.2) microns in Makani and the bandpasses of the MIRI and NIRCam filters. The warm dust is not only detected in the cool-gas tracers of the galactic wind associated with the more recent (7 Myr) starburst episode, but also in the outer warm-ionized gas wind produced by the older (0.4 Gyr) episode. The presence of PAHs in the outer wind indicates that the PAHs have survived the long (R/v ~ 10^8 yrs) journey to the halo despite the harsh environment of the galactic wind. The measured F1800W/F1130W flux ratios in the unresolved nucleus, inner halo (R = 10 - 20 kpc), and outer halo (R = 20 - 35 kpc), tracers of the PAH (11.3 + 12.2)/7.7 ratios, indicate decreasing starlight intensity incident on the PAHs, decreasing PAH sizes, and increasing PAH ionization fractions with increasing distance from the nucleus. These data provide the strongest evidence to date that the ejected dust of galactic winds survives the long journey to the CGM, but is eroded along the way.


[14] 2507.08111

The Kinematic Age of 3I/ATLAS and its Implications for Early Planet Formation

The recent discovery of the third interstellar object (3I/ATLAS) expands the known census from two to three and significantly improves statistical inferences regarding the underlying galactic population. 3I/ATLAS exhibits detectable activity that may increase as the object approaches perihelion. In this paper, we argue that the cometary activity likely significantly contributes to 3I/ATLAS's brightness, since the nuclear size of 3I/ATLAS assuming an asteroidal reflectance implies an untenable interstellar object mass per star. 3I/ATLAS exhibits an excess velocity of $v_\infty=58$ km/s relative to the Sun which is significantly larger than that of 1I/`Oumuamua or 2I/Borisov. This velocity implies that 3I/ATLAS is relatively old in comparison to previous interstellar objects. Here, we calculate the posterior distribution of ages implied by the kinematics of the interstellar objects and find that 3I/ATLAS is likely $\sim3-11$ Gyr old, assuming that the interstellar object and stellar age-velocity relations are equivalent. We also calculate the distribution of host star metallicities and find that 3I/ATLAS likely originated in a lower-metallicity system than the previous interstellar objects. These results show that interstellar object formation is likely efficient at low metallicities and early in the history of the Galaxy. Finally, we estimate and approximate interstellar object formation rate throughout Galactic history implied by these three objects.


[15] 2507.08117

Magnetic properties of orphan penumbrae

Orphan penumbrae (OPU) are features resembling sunspot penumbrae, but are not connected to an umbra. Here we compare OPUs and sunspot penumbrae, including their filaments. We also identify and describe the main mechanisms for the formation of OPUs and we characterise their decay process. Our study is based on spectropolarimetric inversions of active regions observed with the Hinode spectropolarimeter. We manually identified 80 individual OPUs, allowing us to study them statistically. In addition, we analysed the time-evolution of selected OPUs using data provided by the Helioseismic and Magnetic Imager. Orphan penumbrae display a broad range of shapes, associated with typically $\Omega$-shaped magnetic field configurations, where opposite polarity fields predominate at the two ends of the OPU. In addition, the properties of the OPU filaments are remarkably uniform between different OPUs, resembling the ones in sunspot penumbrae. Most OPUs form by either a patch of a penumbra separating from a sunspot, or by new magnetic flux emerging close to the polarity inversion line of an active region. We observe chromospheric fibrils above almost all OPUs in Hinode H$\alpha$ images, indicating that a part of the magnetic field of the OPUs extends to the chromosphere. Our results show that OPU filaments can form given a broad range of boundary conditions for the magnetic field.


[16] 2507.08122

The Blue Lurker WOCS 14020: A Long-Period Post-Common-Envelope Binary in M67 Originating from a Merger in a Triple System

We present Hubble Space Telescope far-ultraviolet (FUV) spectra of a blue-lurker$-$white-dwarf (BL-WD) binary system in the 4 Gyr open cluster M67. We fit the FUV spectrum of the WD, determining it is a C/O WD with a mass of $0.72^{+0.05}_{-0.04}$ M$_\odot$ and a cooling age of $\sim400$ Myr. This requires a WD progenitor of $\sim3$ M$_\odot$, significantly larger than the current cluster turnoff mass of 1.3 M$_\odot$. We suggest the WD progenitor star formed several hundred Myr ago via the merger of two stars near the turnoff of the cluster. In this scenario, the original progenitor system was a hierarchical triple consisting of a close, near-equal-mass inner binary, with a tertiary companion with an orbit of a few thousand days. The WD is descended from the merged inner binary, and the original tertiary is now the observed BL. The likely formation scenario involves a common envelope while the WD progenitor is on the AGB, and thus the observed orbital period of 359 days requires an efficient common envelope ejection. The rapid rotation of the BL indicates it accreted some material during its evolution, perhaps via a wind prior to the common envelope. This system will likely undergo a second common envelope in the future, and thus could result in a short-period double WD binary or merger of a 0.72 M$_\odot$ C/O WD and a 0.38 $M_\odot$ Helium WD, making this a potential progenitor of an interesting transient such as a sub-Chandrasekhar Type Ia supernova.


[17] 2507.08126

SILCC -- IX. The multi-phase interstellar medium at low metallicity

The gas-phase metallicity affects heating and cooling processes in the star-forming galactic interstellar medium (ISM) as well as ionising luminosities, wind strengths, and lifetimes of massive stars. To investigate its impact, we conduct magnetohydrodynamic simulations of the ISM using the FLASH code as part of the SILCC project. The simulations assume a gas surface density of 10 M$_\odot$ pc$^{-2}$ and span metallicities from 1/50 Z$_\odot$ to 1 Z$_\odot$. We include non-equilibrium thermo-chemistry, a space- and time-variable far-UV background and cosmic ray ionisation rate, metal-dependent stellar tracks, the formation of HII regions, stellar winds, type II supernovae, and cosmic ray injection and transport. With the metallicity decreasing over the investigated range, the star formation rate decreases by more than a factor of ten, the mass fraction of cold gas decreases from 60% to 2.3%, while the volume filling fraction of the warm gas increases from 20% to 80%. Furthermore, the fraction of H$_\mathrm{2}$ in the densest regions drops by a factor of four, and the dense ISM fragments into approximately five times fewer structures at the lowest metallicity. Outflow mass loading factors remain largely unchanged, with values close to unity, except for a significant decline at the lowest metallicity. Including the major processes that regulate ISM properties, this study highlights the strong impact of gas phase metallicity on the star-forming ISM.


[18] 2507.08131

Results from IceCube Searches for High-energy Neutrinos Coincident with Gravitational-Wave Alerts in LVK O4

Mergers of compact objects, binary black holes and mergers including at least one neutron star, are a predicted source of high-energy neutrinos. These astrophysical events are now routinely detected through observation of their gravitational wave signature and, at least in one instance, their electromagnetic counterparts were also detected. Particles accelerated during the coalescence of compact objects may also interact to produce high-energy neutrinos, which have yet to be detected, but observations are ongoing. The LIGO-Virgo-KAGRA Collaboration publicly releases information on candidate gravitational wave events from compact binary coalescences in low latency during the current observing run (O4). To aid the electromagnetic follow-up, using data from the IceCube Neutrino Observatory, we search, in real time, for neutrinos spatially and temporally coincident with these gravitational wave candidate events using a time window of 1000 seconds centered on the merger time. We use two methods, both of which have been previously used to search for neutrino emission from gravitational-wave transients: an unbinned maximum likelihood analysis applied to significant alerts and a Bayesian analysis with astrophysical priors, applied to both significant and low-significance alerts. In addition, we search for long-duration neutrino emission up to 14 days after the merging of binaries containing a neutron star. We report analysis results determined in real time for these searches, and set upper limits on both flux and isotropic-equivalent energy emitted in neutrinos.


[19] 2507.08132

Enhancing searches for astrophysical neutrino sources in IceCube with machine learning and improved spatial modeling

Searches for astrophysical neutrino sources in IceCube rely on an unbinned likelihood that consists of an energy and spatial component. Accurate modeling of the detector, ice, and spatial distributions leads to improved directional and energy reconstructions, resulting in increased sensitivity. In this work, we utilize our best knowledge of the detector ice properties and detector calibrations to reconstruct in-ice particle showers. The spatial component of the likelihood is parameterized either by a 2D Gaussian or a von Mises Fisher (vMF) distribution at small and large angular uncertainties, respectively. Here, we use a gradient-boosted decision tree with a vMF spatial likelihood loss function, reparameterized through two coordinate transformations, to predict per-event point spread functions (PSF). Additionally, we discuss the search for PeV cosmic ray sources using the IceCube Multi-Flavor Astrophysical Neutrino (ICEMAN) sample. Our search contains both an analysis of individual neutrino sources coincident with greater than 100 TeV gamma-ray sources and also a stacking analysis. We outline the prospects for extended neutrino emission originating from the Cygnus Cocoon region.


[20] 2507.08133

newASTROGAM -- The New MeV to GeV Gamma-ray Observatory

newASTROGAM is a breakthrough mission concept for the study of the non-thermal Universe from space with gamma rays in the energy range from 15 keV to 3 GeV. It is based on advanced space-proven detector technologies, which will achieve unprecedented sensitivity, angular and energy resolution combined with polarimetric capability. Since the MeV gamma-ray energy range is the most under-explored electromagnetic window to the Universe, a mission in this energy range can for the first time sensitively address fundamental astrophysics questions connected to the physics of compact objects and merger events, jets and their environments, supernovae and the origin of the elements, potentially constrain the nature of dark matter and many more science objectives. The mission will detect and follow-up many of the key sources of multi-messenger astronomy in the 2040s. newASTROGAM provides an unprecedentedly broad energy coverage from keV to GeV energies. The payload concept consists of a Silicon tracker combined with a crystal calorimeter. Both detectors are surrounded by an anti-coincidence detector to reject charged cosmic rays. In addition, a thin X-ray coded mask provides very good imaging capabilities. Such a mission can uniquely detect gamma rays via the photoelectric effect, Compton scattering and electron-positron pair production. newASTROGAM is proposed to the ESA call for medium-class mission ideas (M8).


[21] 2507.08146

Universal energy limits of radiation belts in planetary and brown dwarf magnetospheric systems

Radiation belts are regions of magnetically trapped particle radiation found around all of the sufficiently magnetized planets in the Solar System and recently also observed around brown dwarfs, yet despite their ubiquity, there is not yet a general theory or model to predict the uppermost energy limits that any particular magnetospheric system's radiation belts can attain. By considering only the most fundamental loss processes, a model and corresponding theory are developed that successfully bound and explain the maximum observed energies of all documented radiation belt systems. Interestingly, this approach yields a relatively simple function for the uppermost energy limit that depends on only the surface magnetic field strength of the system. The model predicts an energy limit for all radiation belt systems that asymptotes at 7 +/- 2 TeV (for protons and electrons), offering intriguing new insight on potential sources of galactic cosmic rays. This model is also applied to an exoplanetary system, demonstrating that the planet is likely a synchrotron emitter and showcasing the model's use for identifying candidate targets for synchrotron-emitting astrophysical systems and revealing details critical to habitability at those remote worlds.


[22] 2507.08152

Analysis of the Diffuse Astrophysical Flux from the Galactic Plane with 12.1 Years of IceCube Starting Tracks, Throughgoing Tracks and Cascades

The IceCube Neutrino Observatory has measured an isotropic astrophysical neutrino flux through various detection channels for over 12 years. IceCube has also detected neutrino emission from the Galactic plane at the 4.5$\sigma$ significance level compared to a background-only hypothesis, testing three models of Galactic diffuse emission: Fermi-LAT $\pi^0$, KRA$_{\gamma}^{5}$, KRA$_{\gamma}^{50}$. We present an analysis combining 3 detection channels: throughgoing tracks, starting tracks and cascades. The throughgoing track sample is restricted to the northern sky to reduce atmospheric backgrounds, while the starting track and cascade samples reduce the atmospheric neutrino backgrounds in the southern sky by vetoing accompanying muons. We will use this combination of event samples from 12.1 years of data to measure the Galactic neutrino spectrum in the TeV to PeV energy range and independently for multiple galactic regions in a model independent procedure. We will simultaneously measure the isotropic cosmic neutrino flux.


[23] 2507.08172

A group of merging galaxies falling onto Abell 2142

Galaxy clusters produce a very hostile environment to galaxies, whose gas gets stripped by ram-pressure, suffer galaxy interactions and witness quenching of their star formation. Clusters, like Abell 2142, grow not only through galaxy accretion but also through group infall. Our goal is to study the physical and dynamical state of the most conspicuous infalling group, on a filament projected at 1.3 Mpc from the Abell 2142 centre. The galaxy group is the leading edge of a spectacular trailing 700-kpc-long X-ray tail of hot gas stripped by ram-pressure. The infalling galaxies are not quenched yet, and are ideal objects to study the transformation processes due to the cluster environment. We use integral field spectroscopy from MaNGA to derive stellar and gas kinematics, and MegaCam for ugr photometry. Stellar populations (with age and metallicity) are obtained through full-spectrum fitting using Nburst. The gas kinematics and excitation are derived from the line emission of H$\alpha$, [NII], [OIII] and H$\beta$. The group contains four galaxies, of which two are merging and partly superposing on the line of sight. With a simple parametric model for each velocity field, we succeed in disentangling the contribution of each galaxy and derive their physical state and kinematics. The galaxies are perturbed, and intra-group gas is observed as tidal tails and loops. They are mainly disks in rotation, although some regions reveal elevated dispersion, typical of out-of-equilibrium gas. All galaxies show sustained star formation, with a global star formation rate of 45 M$_\odot$/yr. We conclude that the long X-ray tail must have come from the hot intra-group medium, present before the group infall, and does not correspond to the ram-pressure stripping of the galaxy gas. The galaxy interactions within the group are still enhancing the star formation, from the disks that are still rich in dense gas.


[24] 2507.08176

Development of a laser stabilized on an ultra-stable silicon cryogenic Fabry-Perot cavity for dark matter detection

Astrophysical observations suggest the existence of an unknown kind of matter in the Universe, in the frame of the $\Lambda$CDM model. The research field of dark matter covers an energy scale going from massive objects to ultra-light scalar fields, which are the focus of the present work. It is supposed that ultra-light scalar fields affect the length of objects, whereas the speed of light stays unchanged. It follows that Fabry-Perot cavities are ideal tools for ultra-light dark matter detection since the fluctuations in the length of a cavity can be detected on the frequency of the laser stabilized to it. At FEMTO-ST, we have set up an ultra-stable silicon cavity suitable for a test of detection of ultra-light dark matter in an energy range close to 10$^{-10}$ eV. Our 14 cm cavity is composed of two mirrors optically bonded to an ultra-rigid spacer, with each element made in single-crystal of silicon, and cooled at 17 K in order to cancel the first order thermal expansion coefficient of the silicon spacer. The projected fractional frequency stability of the laser is $3 \times 10^{-17}$, mainly limited by the thermal noise of the amorphous dielectric reflective coatings. To reach this remarkable stability, several effects have to be reduced below the thermal noise limit. While the contribution of the residual amplitude modulation is now acceptable, we are currently implementing a laser power lock with residual fluctuations lower than 3~nW and a piezoelectric-based servo loop to actively reduce the vibration noise that has to be inferior to $-110 \textrm{~dB(m s}^{-2})^2/\textrm{Hz}$ at 1 Hz. Here, we present both the status of the development of our ultra-stable laser and the mechanical response of the cavity in the presence of ultra-light dark matter.


[25] 2507.08179

BASS XLVIII: [Ne v] λ3427 Emission in Powerful Nearby Active Galactic Nuclei

We investigate the high-ionization, narrow [Ne v] $\lambda$3427 emission line in a sample of over 340 ultrahard X-ray (14-195 keV) selected Active Galactic Nuclei (AGN) drawn from the BASS project. The analysis includes measurements in individual and stacked spectra, and considers several key AGN properties such as X-ray luminosity, supermassive black hole (SMBH) mass, Eddington ratios, and line-of-sight column density. The [Ne v] $\lambda$3427 line is robustly detected in ~43% (146/341) of the AGN in our sample, with no significant trends between the detection rate and key AGN/SMBH properties. In particular, the detection rate remains high even at the highest levels of obscuration (>70% for log[N_H/cm^-2] > 23). On the other hand, even some of our highest signal-to-noise spectra (S/N > 50) lack a robust [Ne v] detection. The typical (median) scaling ratios between [Ne v] line emission and (ultra-)hard X-ray emission in our sample are log L[Ne v]/L(14-150 keV) = -3.75 and log L[Ne v]/L(2-10 keV) = -3.36. The scatter on these scaling ratios, of ~0.5 dex, is comparable to, and indeed smaller than, what is found for other commonly used tracers of AGN radiative outputs (e.g., [O III] $\lambda$5007). Otherwise, we find no significant relations between the (relative) strength of [Ne v] and the basic AGN/SMBH properties under study, in contrast with simple expectations from models of SMBH accretion flows. Our results reaffirm the usability of [Ne v] as an AGN tracer even in highly obscured systems, including dual AGN and high redshift sources.


[26] 2507.08201

Detecting Intermediate-mass Black Holes Using Miniature Pulsar Timing Arrays in Globular Clusters

Theoretical models predict that intermediate-mass black holes (IMBHs) exist in globular clusters (GCs), but observational evidence remains elusive. Millisecond pulsars (MSPs), which are abundant in GCs and have served as precise probes for gravitational waves (GWs), offer a unique opportunity to detect potential IMBH binaries in GCs. Here, we consider the possibility of using multiple MSPs in a GC to form a miniature pulsar timing array (PTA), so as to take advantage of their correlated timing residuals to search for potential IMBH binaries in the same cluster. Our semi-analytical calculations reveal that nearby IMBH binaries around MSPs in GCs could induce microsecond-level timing residuals. In GCs like $\omega$ Centauri and M15, favorable configurations are found which could lead to the detection of binaries with mass ratios $q\gtrsim0.1$ and orbital periods of a few days. We estimate that future higher-precision timing programs could achieve $100$-nanosecond sensitivity, substantially expanding the searchable parameter space and establishing mini-PTAs as powerful detectors of IMBHs.


[27] 2507.08219

GW231123: a Binary Black Hole Merger with Total Mass 190-265 $M_{\odot}$

On 2023 November 23 the two LIGO observatories both detected GW231123, a gravitational-wave signal consistent with the merger of two black holes with masses $137^{+22}_{-17}\, M_\odot$ and $103^{+20}_{-52}\, M_\odot$ (90\% credible intervals), at luminosity distance 0.7-4.1 Gpc and redshift of $0.39^{+0.27}_{-0.24}$, and a network signal-to-noise ratio of $\sim$22.5. Both black holes exhibit high spins, $0.9^{+0.10}_{-0.19}$ and $0.80^{+0.20}_{-0.51}$ respectively. A massive black hole remnant is supported by an independent ringdown analysis. Some properties of GW231123 are subject to large systematic uncertainties, as indicated by differences in inferred parameters between signal models. The primary black hole lies within or above the theorized mass gap where black holes between 60-130 $M_\odot$ should be rare due to pair instability mechanisms, while the secondary spans the gap. The observation of GW231123 therefore suggests the formation of black holes from channels beyond standard stellar collapse, and that intermediate-mass black holes of mass $\sim$200 $M_\odot$ form through gravitational-wave driven mergers.


[28] 2507.08237

Evidence of .Ia Supernova Detonations in 3D Hydrodynamical Simulations of Double Degenerate Mergers

We report detailed 3D simulations of 1.1 $\mathrm{M_{\odot}}$ Oxygen-Neon (ONe) white dwarfs (WDs) merging with a 0.35 $\mathrm{M_{\odot}}$ helium WD, conducted with the moving-mesh hydrodynamic code AREPO. The simulations utilise self-consistent chemical profiles for the primary WD which were generated by a stellar evolution code incorporating the effects of semi-degenerate carbon burning. We find that a helium detonation is ignited at the base of the helium layer, starting a thermonuclear runaway which encircles the WD and ejects material as a sub-luminous supernovae. Our canonical simulation, (C-120), ejects 0.103 $\mathrm{M_{\odot}}$ of primarily $^{4}\mathrm{He}$, $^{28}\mathrm{Si}$, and $^{32}\mathrm{S}$, after which the primary begins accreting again from the surviving secondary. Our results depend qualitatively on the "inspiral time" simulation parameter, which describes the length of a period of accelerated angular momentum loss. For example, the binary does not survive when inspiral time is too long. We compare the results using our self-consistent chemical profiles to a constant-composition WD structure and find the same explosion pattern when the inspiral time is short. However, we are able to obtain a typical Type Ia supernova (SN Ia) which destroys the primary by using the constant-composition structure and long inspiral. The shock-convergence in this simulation follows the "x-scissor mechanism" described by Gronow et al. in 2020, and causes a secondary detonation due to higher temperatures and higher number density of $^{12}\mathrm{C}$ at the convergence site. These results highlight the potential for unrealistic outcomes when conducting simulations that incorporate unrealistically large enhancements in angular momentum losses and (or) non-realistic chemical structures for the primary WD.


[29] 2507.08242

Time Variation in the TeV Cosmic Ray Anisotropy with IceCube and Energy Dependence of the Solar Dipole

There is an observed anisotropy in the arrival direction distribution of cosmic rays in the TeV-PeV regime with variations on the scale of one part in a thousand. While the origin of this anisotropy is an open question, a possible factor is cosmic-ray interactions with interstellar and heliospheric magnetic fields. These magnetic fields may change over time - for example, due to changes in solar activity throughout its 11-year solar cycle. The cosmic-ray anisotropy can reflect these time-dependent magnetic fields. In addition to these speculative sources, there are several known sources of time variation in this anisotropy, such as the Compton-Getting Effect from the Earth's orbital motion. We discuss a preliminary study with limited statistics of time variation undertaken by the IceCube Neutrino Observatory, including a measurement of the Compton-Getting Effect as well as a general, model-independent search for other time variations. Further, we use the Compton-Getting Effect to present a preliminary measurement of the cosmic-ray spectral index as a function of energy below the knee.


[30] 2507.08245

The JWST Spectroscopic Properties of Galaxies at $z=9-14$

We characterize the JWST spectra of $61$ galaxies at $z=9-14$, including $30$ newly-confirmed galaxies. We directly compare the $z>9$ spectroscopic properties against $401$ galaxies at $69$ galaxies. Our dataset indicates that the emission lines undergo a marked change at $z>9$, with extremely large CIII], H$\beta$, and H$\gamma$ EWs becoming $2-3\times$ more common at $z>9$ relative to $69$ galaxies. These results may reflect a larger dispersion in UV luminosity at fixed halo mass and larger baryon accretion rates at $z>9$, although other physical effects may also contribute. We suggest that the shift in star formation conditions explains the prevalence of extreme nebular spectra that have been detected at $z>9$, with hard ionizing sources and nitrogen-enhancements becoming more typical at the highest redshifts. Finally, we identify five $z>9$ spectroscopically confirmed galaxies with red UV colors ($\beta\gtrsim-1.5$), either revealing a small population with moderate dust attenuation ($\tau_V=0.23-0.35$) or very high density nebular-dominated galaxies with hot stellar populations.


[31] 2507.08246

JVLA Measurement of Grain Size in the Compact Dust Ring around Class I Protostar WL 17

The maximum grain size in protoplanetary disks is a critical parameter for planet formation, as the efficiency of mechanisms like streaming instability and pebble accretion depend on grain size. Even young class 0/I objects, such as HL Tau, show substructures in their disks, indicating the potential for early planet formation. In this study, we investigated the grain size in the dust surrounding the class I object WL 17 using the Karl G. Jansky Very Large Array. Observations were conducted across seven frequency bands (Q, Ka, K, Ku, X, C, and S bands) ranging from 2 to 48 GHz, corresponding to wavelengths of 15 cm to 6.3 mm, with a spatial resolution exceeding 0\farcs5. While the ring structure at 0\farcs1 of WL 17 remains unresolved in our data, its emission is clearly detected at all observed frequencies, except at 2 GHz. To estimate the maximum grain size ($a_{\rm max}$) within the ring, we compared the observed spectral energy distribution (SED) with theoretical SEDs calculated for various $a_{\rm max}$ values using radiative transfer models. Assuming the dust opacity follows the DSHARP model, our analysis suggests that certain structures internal to the ring achieved a maximum grain size of approximately 4.2 mm. Additionally, we discuss the gravitational stability of the ring and the potential planetary core mass that could form through pebble accretion within the structure.


[32] 2507.08256

Dynamo Confinement of a Radiatively Spreading Solar Tachocline Revealed by Self-consistent Global Simulations

The helioseismically observed solar tachocline is a thin internal boundary layer of shear that separates the rigidly-rotating solar radiative zone from the differentially-rotating convective zone and is believed to play a central role in the 22-year solar dynamo cycle. The observed thinness of the tachocline has long been a mystery, given the expected tendency of such shear to undergo radiative spreading. Radiative spreading is the process by which the meridional circulation and angular velocity burrow into a stably-stratified fluid owing to the mitigating effect of radiative thermal diffusion. A confinement mechanism is thus required to keep the tachocline so thin. In previous work using global dynamo simulations, we achieved a statistically-stationary confined tachocline where the confinement mechanism was derived from the Maxwell stress arising from a dynamo-generated nonaxisymmetric poloidal magnetic field. However, the parameters chosen meant that the tachocline was confined against viscous spread instead of radiative spread. Here, we show that the previously identified dynamo confinement mechanism still succeeds in a simulation that lies in the more solar-like radiative spreading regime. In particular, a nonaxisymmetric, quasicyclic dynamo develops in the convective zone and overshoot layer, penetrates into the radiative zone via a novel type of skin effect, and creates a Maxwell stress that confines the tachocline over many magnetic cycles. To the best of our knowledge, this is the first fully self-consistent rendering of a confined tachocline in a global numerical simulation in the parameter regime appropriate to the Sun.


[33] 2507.08257

SDSS-V LVM: Collisionless Shocks in the Supernova Remnant RCW86

The supernova remnant (SNR) RCW86 is among the few SNRs with Balmer-emission lines containing broad and narrow spectral components that trace fast, non-radiative shocks in partially-ionized gas.\ These are invaluable laboratories for collisionless shock physics, especially for poorly-understood phenomena like electron-ion equilibration, and shock precursors. Here we present the first $\sim$0.3 pc spatial scale integral field unit (IFU) observations of the southwestern RCW86 shock, obtained as part of the Sloan Digital Sky Survey-V Local Volume Mapper (SDSS-V LVM). The forward shock, clearly visible as thin filaments in narrowband images, have broad H$\alpha$ components, indicating shock velocities varying from 500--900 km/s in the south to 1000--1500 km/s in the north. The varying velocity widths and broad-to-narrow intensity ratios show that electrons and ions have lower equilibration ($T_e/T_p \rightarrow 0.1$) in faster ($>$800 km/s) shocks, in line with previous studies. The broad components are generally redshifted from the narrow components by $\lesssim$100 km/s, likely due to shock-obliquity or non-Maxwellian post-shock distributions. We observe high extinction-corrected Balmer-decrements of 3--5 in the narrow components, indicating that conversion of Ly$\beta$ photons to H$\alpha$ is more efficient than Ly$\gamma$ to H$\beta$. Broad HeII$\lambda$4686 was marginally ($\gtrsim$2$\sigma$) detected in the southern shock, meaning the shock is impacting gas with high ($>$30--100\%) neutral fraction. We also find the first evidence of an intermediate H$\alpha$ component in RCW86, with $\Delta$V(FWHM) = 193--207 km/s, likely due to a neutral precursor. We also briefly discuss the southwestern radiative shock, and lay out the exciting future of studying astrophysical shocks with LVM.


[34] 2507.08266

Evaluating the Limits of Rotation Period Recovery through Gyrochronology Criteria

Contamination from nearby sources often compromises stellar rotation periods derived from photometric light curves, particularly in data with large pixel scales such as TESS. This problem is compounded when both the target and contaminant are intrinsically variable, a scenario that challenges deblending algorithms, which often assume constant contaminants. We assess the reliability of rotation period detections using wide binary systems, whose components share a common age and rotational history. By applying gyrochronology constraints, we identify period combinations that yield consistent ages between components, helping to isolate true rotation signals. Simulating blends with degraded Kepler data, our method recovers correct rotation periods with an 88\% success rate for periods $<12$ days, where TESS detections are most reliable. Applying this framework to nearly 300 wide binaries observed by TESS, we find that despite significant contamination, a subset of pairs shows consistent gyrochronological ages. We establish a practical detection threshold for TESS blended observations, finding that periods shorter than $\sim8$ days are reliably recovered, while those longer than $\sim10$ days become significantly more challenging and often remain unresolved. As expected, rotation periods are more often recovered when the highest-amplitude periodogram peak is linked to the brighter star and the second to the dimmer star, although many cases deviate from this pattern, indicating it cannot always be assumed. Our results highlight the limitations of standard deblending methods and demonstrate that astrophysical constraints, such as gyrochronology, provide a valuable tool for extracting reliable rotation periods from complex photometric blends.


[35] 2507.08275

Constraining decaying dark matter models with gravitational lensing and cosmic voids

Despite overwhelming observational evidence for dark matter, we still have no evidence of direct detection. Consequently, our knowledge about dark matter is limited, for example, we do not know if dark matter is a stable particle or if it decays. Without a theoretical particle model, the parameter space of possible decay models is highly variable, and astrophysical or cosmological means of indirectly constraining the phenomenological models are required. This paper investigates a scenario in which a dark matter decays and the dark daughter particle moves with respect to the comoving mother particle. The model is parameterised by the decay rate and the injection velocity of the dark matter particles, which can be converted to the mass ratio. In previous work, a simpler model was used to investigate the evolution of cosmic voids typified as regions with low content of galaxies and non-baryonic matter. It was found that the growth of S-type voids is modified by the dark matter decay, leading to imprints at the present day. Here we extend our study and improve the method used to model the decay. We also study the gravitational weak-lensing signal that will be able to detect or constrain the parameter space of decaying dark matter. The results of this study suggest that future weak lensing surveys may provide unique probes of the phenomological parameters of dark matter.


[36] 2507.08298

Obscured and unobscured X-ray AGNs I: Host galaxy properties

Active galactic nuclei (AGN) play a crucial role in galaxy evolution by influencing the observational properties of their host galaxies. We investigate the host galaxy properties of X-ray selected AGNs, focusing on differences between obscured and unobscured AGNs, and between high-(log([OIII]5007/H$\beta$}$>$0.5) and low-excitation sources (log([OIII]5007/H$\beta$$<$0.5). We selected a sample of AGNs from the spectroscopic zCOSMOS survey with 0.5 $< zsp <$ 0.9 based on the Mass-Excitation (MEx) diagram and X-ray emission. AGNs were classified as obscured or unobscured using hydrogen column density, and as high- or low-excitation based on the [OIII]5007/H$\beta$ ratio. We analysed various AGN properties, including the hardness ratio, X-ray luminosity, emission line ratios such as the ionisation-level sensitive parameter O32=log([OIII]5007/[OII]3727), and the metallicity sensitive parameter R23=log(([OIII]5007+[OII]3727/H$\beta$), and the specific black hole accretion rate. Unobscured AGNs exhibit a more evident correlation between the [OIII]5007/[OII]3727 ionisation ratio and X-ray luminosity than obscured AGNs, while high-excitation obscured AGNs reach, on average, higher X-ray luminosities. Furthermore, high-excitation AGNs typically show high values of R23, suggesting low metallicities, similar to that observed in high-redshift galaxies (4 $< z <$ 6). We find a positive correlation between the parameters $\lambda$sBHAR and N$_H$, R23 and O32 parameters. The correlation suggests that AGNs with a high specific accretion rate have not only a higher production of high-energy photons, which ionise the surrounding medium more intensely, but are also usually associated with environments less enriched in heavy elements. These results provide insights into the complex interplay between AGN activity, host galaxy properties, and the role of obscuration in shaping galaxy evolution.


[37] 2507.08304

Einstein Probe discovery of the short period intermediate polar EP J115415.8-501810

The X-ray transient source EP240309a/EP\,J115415.8$-$501810 was first detected by the Wide-Field X-ray Telescope (WXT) on board Einstein Probe (EP) during the commissioning phase. Subsequent optical observations confirmed it as a Cataclysmic Variable of the intermediate polar type with a 238.2\,s spinning white dwarf in a $\sim$3.76\,hr orbit. We report on the source discovery and follow-up studies made with the Follow-up X-ray Telescope (FXT) of EP. A periodic variation of 231\,s is detected in the 0.3$-$2\,keV band, while no obvious pulsation appears in the 2$-$10\,keV band. The spectral analysis shows that the X-ray emission could be described by an absorbed bremsstrahlung model with $kT$\textgreater\,11\,keV. The partial covering absorption, with an hydrogen column density $N_H$ = 2.0$\times 10^{22}\,\rm cm^{-2}$ and covering fraction around 0.9, is much larger than the interstellar absorption along the line of sight. According to the distance $d = 309.5$\,pc obtained from Gaia parallax, we estimate that the luminosity of this source in the 0.3$-$10\,keV range is $\sim 2\times10^{32}$\,erg\,s$^{-1}$. In addition, phase-resolved spectral analysis reveals that the detected periodic variation is mainly caused by the change in the absorption column density. In this scenario the spin modulation arises due to absorption from the pre-shock accretion flow of the X-ray emitting pole, while the optical radiation is modulated at the orbital side band ($\omega_{\rm spin} - \Omega_{\rm orbit}$) due to reprocessing in regions within the binary system. Due to its unusual transient behaviour for an intermediate polar, we have also searched for radio signals similar to those observed in the new class of long period transients. We derived upper limits with ASKAP (200--300\,$\mu$Jy\,beam$^{-1}$ between 800--1500 MHz) and MWA (40--90\,mJy\,beam$^{-1}$ between 80--300 MHz).


[38] 2507.08323

On-axis afocal telescopes as framework for Cubesat based astronomical imagers and slit-less spectrographs

Cubesats present unique opportunities for observational astronomy in the modern era. They are useful in observing difficult-to-access wavelength regions and long-term monitoring of interesting astronomical sources. However, conventional telescope designs are not necessarily the best fit for restricted envelope of a Cubesat. Additionally, fine-pointing stability on these platforms is difficult due to the low mass of the spacecraft and special allocations within the optical design are needed to achieve stable pointing. We propose afocal telescope designs as the framework to realise imagers and low-resolution spectrographs on Cubesat platforms. These designs help reduce the number of components in the optical chain and aim to improve throughput and sensitivity compared to conventional designs. Additionally, they also provide a fine steering mechanism within a collimated beam section. Fine beam steering within the collimated beam section avoids issues of image degradation due to out-of-plane rotation of the image plane or offset in the rotation axis of the mirror. This permits the use of simple and mostly off-the-shelf tip-tilt mirrors for beam steering. The designs discussed here also allow for a standard telescope design to be used in many instrument types; thus reducing the complexity as well as the development time and cost. The optical design, performance and SNR estimations of these designs along with some interesting science cases are discussed. A number of practical aspects in implementation such as guiding, tolerancing, choice of detectors, vibration analysis and laboratory test setups are also presented.


[39] 2507.08372

Spectroscopic Observations of Four Candidates for Blue Large-Amplitude Pulsators. No BLAPs at High Galactic Latitudes

We have obtained spectroscopic observations for four short-period variable objects detected in ZTF, Gaia, and Pan-STARRS (ZGP) data and classified as Blue Large-Amplitude Pulsators (BLAPs) in McWhirter and Lam (2022): ZGP-BLAP-03, ZGP-BLAP-04, ZGP-BLAP-10, and ZGP-BLAP-15. The variables have periods between 46 and 56 min, full amplitudes of 0.13-0.22 mag in the $r$ band, and light curve shapes typical for radially pulsating stars. Three of them were found at high galactic latitudes (|$b$|>30 deg). We have identified object ZGP-BLAP-03 as an early F-type star, while objects ZGP-BLAP-04 and ZGP-BLAP-15 as low-metallicity late A-type stars. These are the three objects found at high galactic latitudes and located several kiloparsecs from the Sun. Thus, they are SX Phoenicis-type variable stars residing in the Galactic halo. In the case of low-latitude object ZGP-BLAP-10, we report the presence of helium lines in its spectrum and atmospheric parameters in agreement with known BLAPs. This and other results indicate that BLAPs are absent in metal-poor environments.


[40] 2507.08405

How a Klein-Nishina Modified Eddington limited accretion explains rapid black hole growth in the early universe

The discovery of quasars and their supermassive black holes (SMBHs) over $10^{9} M_{\odot}$ merely hundreds of millions of years after the Big Bang generates tension with the idea of Eddington-limited accretion and pressures the community into exploring the concept of massive black hole seeds and/or super-Eddington accretion. The observation that many black holes have reached supermassive status while obeying the Eddington limit is puzzling as accretion models are not spherically symmetric. We address this issue by illustrating the physics behind a picture of inner disk accretion involving a geometrically thick, hot quasi-spherical flow and argue that such an inner region provides the radiation that instantiates the Eddington limit. Given the energetics of the inner disk edge, we show how the characteristic electron cross-section drops below its Thomson value, allowing black holes to grow rapidly despite being Eddington-limited. Indeed, after implementing a modified cross-section calculated via the Klein-Nishina Formula, we find that SMBH formation time drops by up to $47\%$. In this context, we show how a $10^{9} M_{\odot}$ black hole can form from a seed $10 M_{\odot}$ black hole within $500$ Myr by way of accretion and mergers. While our picture is over-simplified and contrived in a number of ways that we discuss, we suggest that our scenario is interesting in that it offers a solution to two issues at the intersection of astrophysics and cosmology, namely the reason the Eddington limit is obeyed and how some black holes have grown rapidly despite that limit.


[41] 2507.08413

NuAncestor: an artificial satellite-borne star for accurate frequency calibration of ground-based EPRV spectrographs

The accuracy of state-of-the-art Extreme Precision Radial Velocity (EPRV) spectrographs depends on the access to extremely precise and stable wavelength calibration sources. There are several available calibration sources (e.g., emission lamps, laser frequency combs, reference cavities) that can be used to calibrate an astronomical spectrograph. However, the calibration as it is currently performed is always 'local'. In the proposed talk we will present the NuAncestor concept that proposes an accurate (absolute) and common wavelength calibration for astronomical high-resolution, high-precision spectrographs by embarking an optical frequency comb on-board a satellite equipped with an actively pointing telescope and precision orbitography. This calibration satellite shall be available and serve EPRV spectrographs in all major observatories around the world.


[42] 2507.08415

The Optical Sensor for IceCube-Gen2

An innovative optical module (OM) with segmented light-sensitive area has been developed for IceCube-Gen2 that will take neutrino astronomy at the South Pole to the next level. It builds on the successful features of the mDOM and D-Egg modules of IceCube Upgrade while adapting to the smaller borehole diameter of IceCube-Gen2. The newly developed OM, which is being tested in IceCube Upgrade, serves as a prototype for the planned mass production of about 10,000 OMs for IceCube-Gen2. To simplify the assembly process, important changes were made to the design, in particular to integrate the new gel pad concept. This replaces the 3D-printed support structure of the mDOM while maintaining through total internal reflection the increased light collection efficiency of the reflector rings. In addition, the design features local generation of high voltage for each photomultiplier tube (PMT) via a Cockcroft-Walton circuit and the full digitization of the signal on each PMT base with a sampling rate of 60 MSpS. This significantly reduces the complexity of the mainboard so that it fits into the limited space available. This article describes the development status and presents the performance of the first prototypes.


[43] 2507.08421

Astrophysical neutrino flux measurement and search for tau neutrino induced cascades with 11 years of IceCube data

IceCube measured the diffuse astrophysical neutrino flux for all flavors up to PeV energies. The high energy (TeV-PeV) IceCube cascade sample is particularly effective at selecting electron and tau neutrinos. We present the results of Single Power Law (SPL) and Broken Power Law (BPL) flux measurements based on 11 years of cascade data. From this cascade sample, we study the identification of high energy (~100 TeV) tau neutrinos by detecting a double cascade signature produced by a charged-current neutrino interaction and the subsequent decay of the tau lepton. A Boosted Decision Tree (BDT) is employed to search for the double cascade signature, achieving a significantly improved signal-to-background ratio of 9:1 compared to previous analyses. The selected sample has a tau-neutrino purity of approximately 90% and a weighted mean reconstruction error on the tau decay length of about 4m. We present sensitivities for a maximum likelihood fit of the flavor composition and constraints on the astrophysical neutrino flavor ratios using this sample in combination with IceCube's northern muon neutrino-induced track sample.


[44] 2507.08430

A co-deployed dust-logging instrument for the IceCube Upgrade and IceCube-Gen2

A precise understanding of the optical properties of the instrumented Antarctic ice sheet is crucial to the performance of optical Cherenkov telescopes such as the IceCube Neutrino Observatory and its planned successor, IceCube-Gen2. One complication arising from the large envisioned footprint of IceCube-Gen2 is the larger impact of the so-called ice tilt. It describes the undulation of ice layers of constant optical properties within the detector. In this contribution, we will describe the project to build a co-deployed laser dust logger. This is a device to measure the stratigraphy of impurities in the ice to derive the ice tilt. It consists of a light source that will be co-deployed with the photosensor modules, meaning it is part of the deployment string and operated during the deployment of the detector. The newly developed device will be tested during the deployment of the IceCube Upgrade in the 2025/26 austral summer to pave the way for IceCube-Gen2.


[45] 2507.08436

The Rosetta Stone Project. I. A suite of radiative magnetohydrodynamics simulations of high-mass star-forming clumps

Context. Star formation and, in particular, high-mass star formation are key astrophysical processes that are far from being fully understood. Unfortunately, progress in these fields is slow because observations are hard to interpret as they cannot be directly compared to numerical simulations. Synthetic observations are therefore necessary to better constrain the models. Aims. With the Rosetta Stone project, we aim to develop an end-to-end pipeline to compare star formation simulations with observations as accurately as possible in order to study the evolution from clumps scales to stars. Methods. Using the adaptive mesh-refinement code RAMSES, we computed a first grid of model of star-forming clumps to develop our pipeline and explore the impact of the clump initial conditions on their evolution. The main purpose of this set of simulations is to be converted into synthetic observations to enable a direct comparison with real star-forming clumps observed with Herschel and ALMA. Results. The Rosetta Stone simulations presented here provide a catalog available for full post-processing and subsequent comparison with observations (RS1). Among all the parameters explored here, the strength of the magnetic field has the strongest influence on the clump evolution (fragmentation, star formation, global collapse) at both large and small scales. Numerical parameters such as the resolution per Jeans length or the threshold for accretion onto sink particles affects the formation of low-mass sinks. Finally, the widely used L/M ratio is found to be a good indicator of the clump evolutionary state regardless of its initial condition, but this could change when more feedback processes (jets, HII regions) are included. Conclusions. We now have a new suite of simulations of star-forming clumps that is available for full post-processing and subsequent comparison with the observations,


[46] 2507.08437

Plasma instability in the front of ejected energetic electrons and Type III solar radiobursts

Type III radio bursts are a signature of the flux of near-relativistic electrons ejected during solar flares. These bursts are frequently observed by spacecraft such as the Parker Solar Probe. It is traditionally believed that these electron beams generate Langmuir waves through the two-stream instability, which are then converted into electromagnetic waves. In this study, we revise that model by examining how the electron distribution becomes truncated due to the "time-of-flight" effect as the beam travels through a randomly inhomogeneous, and gently varying solar-wind plasma. Rather than the two-stream instability, this truncation destabilizes the distribution and leads to the generation of Langmuir waves via a linear instability; we confine our analysis to this linear regime and do not take into account the back reaction of the generated Langmuir waves on the electron distribution, which is nonlinear. The instability grows until slower electrons arrive and dampen the waves. Our qualitative analysis shows that the resulting wave intensity growth and decay closely match the intensity-time profile of observed Type III radio bursts at the fundamental frequency, supporting this modified theory.


[47] 2507.08452

Ion acceleration in Fermi-LAT behind-the-limb solar flares: the role of coronal shock waves

We investigate the relationship between the gamma-ray emission measured with Fermi-LAT and radio signatures of coronal shock waves in four behind-the-limb (BTL) solar flares. All events were associated with metric type II radio burst. Both start and end times of the radio bursts were synchronized with the gamma-ray emission. The type II bursts associated with the BTL gamma-ray flares had higher speeds and lower formation heights than those of an average sample. These findings support the notion that the highly relativistic ions that produce the gamma-rays in BTL flares are accelerated at CME-driven propagating coronal shock waves rather than in large-scale coronal loops.


[48] 2507.08464

An Earth-Sized Planet in a 5.4h Orbit Around a Nearby K dwarf

We present the discovery and confirmation of the ultra-short period (USP) planet TOI-2431 b orbiting a nearby ($d\sim36$ pc) late K star ($T_{\mathrm{eff}}$ = $4109 \pm 28 \, {\rm K}$) using observations from the Transiting Exoplanet Survey Satellite (TESS), precise radial velocities with the NEID and the Habitable-zone Planet Finder (HPF) spectrographs, as well as ground-based high contrast imaging from NESSI. TOI-2431 b has a period of 5 hours and 22 minutes, making it one of the shortest-period exoplanets known to date. TOI-2431 b has a radius of $1.536 \pm 0.033\, \rm{R_\oplus}$, and a mass of $6.2 \pm 1.2\, \rm{M_\oplus}$, suggesting it has a density compatible with an Earth-like composition and, due to its high irradiation, is likely a 'lava-world' with a $T_{\mathrm{eq}}$ = $2063 \pm 30 \, {\rm K}$. We estimate that the current orbital period is only 30% larger than the Roche-limit orbital period, and that it has an expected orbital decay timescale of only $\sim$31 Myr. Finally, due to the brightness of the host star ($V = 10.9$, $K = 7.6$), TOI-2431 b has a high Emission Spectroscopy Metric of 27, making it one of the best USP systems for atmospheric phase-curve analysis.


[49] 2507.08469

Ion-rich acceleration during an eruptive flux rope event in a multiple null-point configuration

We report on the $\gamma$-ray emission above 100~MeV from the GOES M3.3 flare SOL2012-06-03. The hard X-ray (HXR) and microwave emissions have typical time profiles with a fast rise to a well-defined peak followed by a slower decay. The $>$100~MeV emission during the prompt phase displayed a double-peaked temporal structure with the first peak following the HXR and microwaves, and the second one, about three times stronger, occurring $17 \pm 2$ seconds later. The time profiles seem to indicate two separate acceleration mechanisms at work, where the second $\gamma$-ray peak reveals a potentially pure or at least largely dominant ion acceleration. The Atmospheric Imaging Assembly imaging shows a bright elliptical ribbon and a transient brightening in the north-western (NW) region. Nonlinear force-free extrapolations at the time of the impulsive peaks show closed field lines connecting the NW region to the south-eastern part of the ribbon and the magnetic topology revealed clusters of nulls. These observations suggest a spine-and-fan geometry, and based on these observations we interpret the second $\gamma$-ray peak as being due to the predominant acceleration of ions in a region with multiple null points. The $>$100 MeV emission from this flare also exhibits a delayed phase with an exponential decay of roughly 350 seconds. We find that the delayed emission is consistent with ions being trapped in a closed flux tube with gradual escape via their loss cone to the chromosphere.


[50] 2507.08479

Beyond LambdaCDM: How the Hubble tension challenges early universe physics

Differences in the values of the Hubble constant obtained from the local universe and the early universe have resulted in a significant tension. This tension signifies that our understanding of cosmology (physical processes and/or cosmological data) is incomplete. Some of the suggested solutions include physics of the early Universe. In this paper we aim to investigate common features of various early universe solutions to the Hubble constant tension. The physics of the early universe affects the size of the sound horizon which is probed with the Cosmic Microwave Background (CMB) data. Within the standard model, the size of the horizon (within limits of current measurements) is affected by processes that could occur between (approximately) 1 day after the Big Bang and the last scattering instant. We focus on simple extensions incorporating Early Dark Energy (EDE) and show how such a model affects the inferred values of the Hubble constant. We compare this model to LambdaCDM models using MCMC analysis, likelihoods over the parameter space and Bayesian evidence. The MCMC analysis shows that EDE leads to a decrease in the size of the sound horizon that is consistent with H0 = 73.56 km/s/Mpc but we also show that MCMC analysis favours increasing redshift and proportion of EDE. The Bayesian evidence favours our EDE model for very narrow, finely-tuned parameter space. The LambdaCDM model used for comparison has good evidence across a wide parameter space. We interpret this as an indication that more sophisticated models are required. We conclude that if the Hubble tension were to be related to the physics of the early universe, EDE could be used as a window to explore conditions of the early universe and extend our understanding of that era.


[51] 2507.08500

Determination of the energy scale of cosmic ray measurements using the Auger Engineering Radio Array

The accurate determination of the absolute energy scale in cosmic ray measurements is both a challenging and fundamentally important task. We present how measurements of radio pulses from extensive air showers with the Auger Engineering Radio Array, combined with per-event simulations of radio emission using the CoREAS extension of CORSIKA, allow us to determine the energy scale of cosmic rays between $3\cdot 10^{17}$\,eV and several $10^{18}$\,eV. Our analysis accounts for many factors, each of which is controlled on the 5\% level or better. The absolute calibration of the antennas and the entire analog signal chain builds on a Galactic calibration in combination with a detailed understanding of the antenna-gain patterns. Additional key elements include compensation for temperature-dependent signal amplification, continuous detector health monitoring, an active veto for thunderstorm conditions, an unbiased event reconstruction, and per-event atmospheric modeling in the simulations. The analysis benefits from a high-statistics dataset of over 800 measured cosmic ray showers. We describe our analysis method, perform multiple cross-checks, and evaluate systematic uncertainties. We find that absolute energies determined with AERA are 12\% higher than those established with the Auger Fluorescence Detector, a result well in agreement within systematic uncertainties and thus a strong independent confirmation of the absolute energy scale of the Pierre Auger Observatory.


[52] 2507.08504

Expanding the Pierre Auger Observatory Open Data program

Since 2021, the Open Data Portal has provided access to the Pierre Auger Observatory's data for both the scientific community and the general public. The data release process has been in place since the Observatory's foundation. It continues to be strengthened as outlined in the approved policy and the Observatory's Data Management Plan. More than 80 000 cosmic-ray events above $10^{17}$ eV, detected with the surface and fluorescence detectors, have been released at various levels, from calibrated traces to high-level reconstruction parameters. Additionally, atmospheric data and low-energy particle counting rates have been made available for space weather studies. The Collaboration is committed to releasing FAIR (Findable, Accessible, Interoperable, and Reusable) data, along with accompanying software and detailed documentation, enabling users to perform their own queries and analyses for both research and educational purposes. These datasets have already served as a basis for several scientific papers and have been widely used in various outreach activities. After 20 years of stable data acquisition, the Pierre Auger Collaboration will disclose 30% of the cosmic ray events above $2.5 \cdot 10^{18}$ eV collected with the main surface detector array between 2004 and 2022, corresponding to an exposure of about 24 000 km$^2$ sr yr, together with events detected with the fluorescence detector and used for energy calibration. This release will provide an unprecedented public dataset for ultra-high-energy cosmic rays, enabling in-depth studies of their properties. Together with the published catalog of the 100 most energetic events recorded, this initiative represents the Pierre Auger Collaboration's strong commitment to distributed and collective knowledge, sharing progress with the entire scientific community.


[53] 2507.08509

Investigating field burial by magnetically confined accretion mounds on Neutron Stars

We explore the problem of magnetic confinement of accreted matter forming an accretion mound near the magnetic poles of a neutron star. We calculate the magnetic field geometry of the accreted mound by solving the magnetostatic Grad Shafranov (GS) equation in radially stretched spherical coordinates with high resolution and an extended domain. In this work, we propose a new physically motivated multipolar current free boundary condition at the outer radial boundary. We have evaluated a large suite of GS solutions for different neutron star magnetic fields and mound configurations. We find that with sufficient resolution, the ring-shaped mound profiles spread latitudinally on the neutron star surface, towards the equator, with a potential decline in dipole moment at outer radii, demonstrating the onset of field burial. A higher latitudinal spread towards the equator leads to more effective magnetic field burial. Along with the ring-shaped mound profile on a hard crust majorly used in this work, we also model mounds formed on a pre-existing ocean, which is more physically motivated. Additionally, we explore different GS solutions for a quadru-dipolar surface magnetic field. We find that such configurations lead to asymmetric polar mounds. We discuss the validity of such solutions for different relative strengths of the quadrupole and dipolar components.


[54] 2507.08510

Recent Highlights from the Auger Engineering Radio Array

The Auger Engineering Radio Array (AERA) consists of 153 autonomous antenna stations deployed over 17 km^2 to measure the radio emission from extensive air showers initiated by cosmic rays with energies between 0.1 and 10 EeV in the 30 to 80 MHz frequency band. It operates in coincidence with the other detectors of the Pierre Auger Observatory particularly the Surface Detector (SD) and the Fluorescence Detector (FD). As the largest cosmic-ray radio detector worldwide before the recent Observatory upgrade AugerPrime, AERA has played a pioneering role in the development of radio technique for cosmic rays, providing complementary measurements and serving as a testbed for ideas that motivated the build-up of a new Radio Detector (RD) as part of AugerPrime. We report on measurements of the depth of the shower maximum using the radio footprint, demonstrating compatibility and competitive resolution with established FD results. An absolute calibration of AERA is achieved by monitoring the sidereal modulation of the diffuse Galactic radio emission for nearly a decade, confirming the long-term stability of a radio detector with no significant aging effects observed. This stability suggests that radio detectors could also be used to monitor potential aging effects in other detector systems. Additionally, we investigate the muon content of inclined air showers using hybrid SD-AERA events. Our results indicate that the muon content in measured data is consistent with expectations for iron nuclei as predicted by current-generation hadronic interaction models, confirming the well-known muon deficit for the first time with radio data. These findings reinforce the value of radio detection for cosmic-ray studies and provide a foundation for the next generation of analyses with the AugerPrime RD.


[55] 2507.08511

Geometric Considerations in Hot Jupiter Magnetic Drag Models

Magnetic fields are expected to impact the atmospheric dynamics of hot and ultra-hot Jupiters due to their increased ionization fractions, compared to that of cooler exoplanets, but our ability to model these magnetic processes is limited by the different coupling regimes between the day and night sides of the planets. One common approach is to approximate the magnetic interactions as a drag acting on the atmosphere. In this work, we examine, within the context of this drag approximation, the impact of including vertical and meridional drag, in addition to zonal drag, from a background dipole magnetic field on the flows in hot Jupiter atmospheres as well as a relaxation of the assumption of solely meridional currents and demonstrate that the inclusion of meridional and vertical drag can limit flows over the poles in hotter atmospheres, something not seen in models that only consider zonal drag, and the assumption of only meridional currents results in an underestimation of the equatorial drag in all cases examined.


[56] 2507.08519

Limits on WIMP-Scattering Cross Sections using Solar Neutrinos with Ten Years of IceCube Data

Although dark matter (DM) comprises 84\% of the matter content of the Universe, its nature remains unknown. One broad class of particle DM motivated by extensions of the Standard Model (SM) is weakly interacting massive particles (WIMPs). Generically, WIMPs will scatter off nuclei in large celestial bodies such as the Sun, thus becoming gravitationally bound. Subsequently, WIMPs can annihilate to stable SM particles, ultimately releasing most of their energy as high-energy neutrinos which escape from the Sun. Thus, an excess of neutrinos from the Sun's direction would be evidence for WIMPs. The IceCube Neutrino Observatory is well-suited to such searches since it is sensitive to WIMPs with masses in the region preferred by supersymmetric extensions of the SM. I will present the results of IceCube's most recent solar WIMP search, which includes all neutrino flavors, covers the WIMP mass range from 20 GeV to 10 TeV, and has world-leading sensitivity over this entire range for most channels considered.


[57] 2507.08524

Simulation and Performance Studies for the Tau Air-Shower Mountain-Based Observatory

While IceCube's detection of astrophysical neutrinos at energies up to a few PeV has opened a new window to our Universe, much remains to be discovered regarding these neutrinos' origin and nature. In particular, the difficulty of differentiating electron- and tau-neutrino charged-current (CC) events limits our ability to measure precisely the flavor ratio of this flux. The Tau Air-Shower Mountain-Based Observatory (TAMBO) is a next-generation neutrino observatory capable of producing a high-purity sample of tau-neutrino CC events in the energy range from 1 PeV--100 PeV, i.e. just above the IceCube measurements. An array of water Cherenkov tanks and plastic scintillators deployed in the Colca Canyon will observe the air-shower produced when a tau lepton, produced in a tau-neutrino CC interaction, emerges from the opposite face and decays in the air. In this contribution, I will present the performance studies for TAMBO -- including the expected rates, effective areas, and discrimination potential -- as well as the simulation on which these studies are based.


[58] 2507.08526

Anisotropic Diffusion of $e^\pm$ in Pulsar Halos over Multiple Coherence of Magnetic Field

The slow particle diffusion in pulsar halos, inferred from TeV gamma-ray surface brightness profiles, is attributed to cross-field diffusion under the anisotropic diffusion model. This model assumes sub-Alfvénic interstellar turbulence in the surrounding medium of the pulsar and a rough alignment of the line-of-sight of observers towards the pulsars with the local mean magnetic field direction in the halo. In this model, the expected morphology of a pulsar halo is highly dependent on the properties of the interstellar magnetic field. We investigate the anisotropic diffusion of electron-positron pairs across multiple coherence of magnetic fields in pulsar halos in this work. We focus particularly on their influences on the predicted gamma-ray surface brightness profile and the asymmetry of the halo's morphology, as well as the observational expectations by the Large High Altitude Air Shower Observatory (LHAASO). Our results indicate that the requirement of a specific magnetic field geometry can be alleviated when accounting for a limited (and realistic) coherence length of the magnetic field in the model. Also, the halo's morphology may appear less asymmetric, especially after being smoothed by the point spread function of instruments. It largely relaxes the tension between the asymmetric morphology of halos predicted by the model and lack of apparent asymmetric halos detected so far. Our findings demonstrate the important influence of the coherence length of interstellar magnetic field on the distribution of particles around their accelerators, and the consequence on the measured source morphology.


[59] 2507.08532

A thermonuclear supernova interacting with hydrogen- and helium-deficient circumstellar material. SN 2020aeuh as a SN Ia-CSM-C/O?

Identifying the progenitors of thermonuclear supernovae (Type Ia supernovae; SNe Ia) remains a key objective in contemporary astronomy. The rare subclass of SNe Ia that interacts with circumstellar material (Type Ia-CSM) allows for studies of the progenitor's environment before explosion, and generally favours single-degenerate progenitor channels. The case of SN Ia-CSM PTF11kx clearly connected thermonuclear explosions with hydrogen-rich CSM-interacting events, and the more recent SN 2020eyj connected SNe Ia with helum-rich companion progenitors. Here we present a study of SN 2020aeuh, a Type Ia-CSM with delayed interaction. We analyse photometric and spectroscopic data that monitor the evolution of SN 2020aeuh and compare its properties with those of peculiar SNe Ia and core-collapse SNe. At early times, the evolution of SN 2020aeuh resembles a slightly overluminous SN Ia. Later, the interaction-dominated spectra develop the same pseudocontinuum seen in Type Ia-CSM PTF11kx and SN 2020eyj. However, the later-time spectra of SN 2020aeuh lack hydrogen and helium narrow lines. Instead, a few narrow lines could be attributed to carbon and oxygen. We fit the pseudobolometric light curve with a CSM-interaction mode, yielding a CSM mass of 1-2 M$_{\odot}$. We propose that SN 2020aeuh was a Type Ia supernova that eventually interacted with a dense medium which was deficient in both hydrogen and helium. Whereas previous SNe Ia-CSM constitute our best evidence for nondegenerate companion progenitors, the CSM around SN 2020aeuh is more difficult to understand. We include a hydrodynamical simulation for a double-degenerate system to showcase how the dynamical evolution of such a progenitor scenario could produce the CSM observed around SN 2020aeuh. It is clear that SN 2020aeuh challenges current models for stellar evolution leading up to a SN Ia explosion.


[60] 2507.08534

A survey for radio pulsars and transients in the 10 pc region around Sgr A*

Here we report on a new survey for pulsars and transients in the 10 pc region around Sgr A* using the Effelsberg radio telescope at frequencies between 4 to 8 GHz. Our calibrated full-Stokes data were searched for pulsars and transients using PulsarX, TransientX and PRESTO. Polarisation information is used in the scoring of the candidates. Our periodicity acceleration and jerk searches allowed us to maintain good sensitivity towards binary pulsars in $\gtrsim$ 10-hr orbits. In addition we performed a dedicated search in linear polarisation for slow transients. While our searches yielded no new discovery beyond the redetection of the magnetar SGR J1745-2900, we report on a faint single pulse candidate in addition to several weak periodicity search candidates. After thoroughly assessing our survey's sensitivity, we determined that it is still not sensitive to a population of millisecond pulsars. Next generation radio interferometers can overcome the limitations of traditional single-dish pulsar searches of the Galactic Centre.


[61] 2507.08544

Update on the intermediate arrival-direction analyses of the Pierre Auger Observatory

We present an update on the arrival-direction analyses conducted on intermediate angular scales using the complete Phase I data of the Pierre Auger Observatory up to the end of $2022$ with a total exposure of $135{,}000\,\text{km}^2\,\text{sr}\,\text{yr}$. We show the arrival-direction distribution of the ultra-high-energy cosmic rays along the supergalactic plane above $20\,\text{EeV}$, and an update in the search for magnetically-induced signatures in the arrival directions. Furthermore, we present the potential of introducing estimators for the rigidity ordering of the events to enhance arrival-direction analyses on small to intermediate angular scales. To achieve this, we take advantage of two estimators working on the response of the surface detector: an analytical fit based on air-shower universality and a deep neural network.


[62] 2507.08545

Euclid: Early Release Observations. A combined strong and weak lensing solution for Abell 2390 beyond its virial radius

Euclid is presently mapping the distribution of matter in the Universe in detail via the weak lensing (WL) signature of billions of distant galaxies. The WL signal is most prominent around galaxy clusters, and can extend up to distances well beyond their virial radius, thus constraining their total mass. Near the centre of clusters, where contamination by member galaxies is an issue, the WL data can be complemented with strong lensing (SL) data which can diminish the uncertainty due to the mass-sheet degeneracy and provide high-resolution information about the distribution of matter in the centre of clusters. Here we present a joint SL and WL analysis of the Euclid Early Release Observations of the cluster Abell 2390 at z=0.228. Thanks to Euclid's wide field of view of 0.5 deg$^$2, combined with its angular resolution in the visible band of 0."13 and sampling of 0."1 per pixel, we constrain the density profile in a wide range of radii, 30 kpc < r < 2000 kpc, from the inner region near the brightest cluster galaxy to beyond the virial radius of the cluster. We find consistency with earlier X-ray results based on assumptions of hydrostatic equilibrium, thus indirectly confirming the nearly relaxed state of this cluster. We also find consistency with previous results based on weak lensing data and ground-based observations of this cluster. From the combined SL+WL profile, we derive the values of the viral mass $M_{200} = (1.48 \pm 0.29)\times10^{15}\, \Msun$, and virial radius $r_{200} =(2.05\pm0.13 \, {\rm Mpc}$), with error bars representing one standard deviation. The profile is well described by an NFW model with concentration c=6.5 and a small-scale radius of 230 kpc in the 30\,kpc $< r <$ 2000\,kpc range that is best constrained by SL and WL data. Abell 2390 is the first of many examples where Euclid data will play a crucial role in providing masses for clusters.


[63] 2507.08556

Event reconstruction with the Radio detector of the Pierre Auger Observatory

The surface detector of the Pierre Auger Observatory has recently been upgraded with the addition of radio antennas, forming the radio detector (RD). This contribution outlines the standard methods for reconstructing extensive air showers using the RD, along with recent developments. The reconstruction pipeline is based on a robust understanding of the detector itself. The entire instrument, including the antenna pattern and analog chain, has been meticulously characterized within the Offline software framework, based on measurements in the laboratory as well as in the field. To ensure data integrity, stations identified as unreliable through monitoring are excluded before event reconstruction. Absolute calibration is achieved at the 5 percent level by analyzing the diffuse galactic radio emission. Next, the electric field that induced voltages in the antenna is calculated by unfolding the antenna response pattern. Key observables, such as the energy fluence (the energy deposited in the ground per unit area) and the arrival time of the pulse, are then determined. With these quantities, shower parameters can be reconstructed with very good accuracy in two chi-square-minimization fits: one to determine the shower's arrival direction via a spherical wavefront fit (predicted within 0.2 degree), and the other to estimate the distance to the shower maximum and the electromagnetic cascade energy (predicted within 5 percent) using a lateral density function


[64] 2507.08561

Tracing Missing Baryons in the Cosmic Filaments with tSZ and CMB-Lensing Stacking

We investigate the distribution of missing baryons in cosmic filaments by stacking approximately 31,300 filaments across the northern and southern SDSS sky regions using Planck Compton-$y$ and CMB lensing maps. Filaments are identified using the DisPerSE algorithm applied to the SDSS LOWZ-CMASS galaxy samples, selecting structures with lengths between 30-100 cMpc and redshifts in the range $0.2 < z < 0.6$. Radial profiles are extracted out to 25 cMpc from the filament spines, and galaxy clusters with halo masses above $3 \times 10^{13}$ solar masses are masked to reduce contamination. We detect the thermal Sunyaev-Zel'dovich (tSZ) signal at $7.82\sigma$ and the CMB lensing signal at $7.78\sigma$. The stacked profiles are corrected by a geometric bias related to filament inclination with respect to the line of sight, and modeled assuming isothermal, cylindrically symmetric gas distributions. We explore different gas and matter density distributions, focusing on $\beta$-models with $(\alpha, \beta) = (2, 2/3)$ or $(1, 1)$. By jointly fitting the Compton-$y$ and lensing convergence profiles, we constrain the central electron overdensity and temperature to be $\delta = 5.90^{+4.18}_{-1.98}$ and $T_e = 2.71^{+0.63}_{-0.51} \times 10^6$ K for the standard $\beta$-model. These results suggest that the filamentary WHIM in our selected long filaments contributes a significant baryon fraction of $0.124^{+0.020}_{-0.021} \times$ {\Omega}_b to the cosmic baryon budget.


[65] 2507.08564

A composition-informed search for large-scale anisotropy with the Pierre Auger Observatory

The large-scale dipolar structure in the arrival directions of ultra-high-energy cosmic rays with energies above $8\,$EeV observed by the Pierre Auger Collaboration is a well-established anisotropy measurement. This anisotropy is understood to be of extragalactic origin, as the maximum of the dipolar component is located ${\sim}115^\circ$ away from the Galactic Center. Cosmic rays interact with background radiation and magnetized regions on their path from their sources to Earth. These interactions, which depend on the cosmic-ray energy, charge and mass composition, give rise to different horizons and deflections that are expected to lead to different anisotropies in the arrival directions of cosmic rays at Earth. The Auger Collaboration has determined that the mass composition of cosmic rays at ultra-high energies is mixed, becoming increasingly heavier as a function of energy. Thus, different dipole amplitudes are expected to be measured at a given energy when separating the data into composition-distinct subsets of lighter and heavier elements. In this contribution, we investigate the composition signature on the large-scale anisotropy taking advantage of composition estimators obtained from the data gathered with the surface detector. A way of probing for composition signatures in anisotropy patterns is then to divide the data into subsets of ``lighter'' and ``heavier'' elements per energy bin. In a simulation library, we evaluate the possibility of measuring a separation in total dipole amplitude between two such populations of the measured dataset under a source-agnostic model. We present the results using two different composition estimators, one based on air-shower universality and one inferred with deep learning.


[66] 2507.08569

Prospects for sub-GeV astrophysical neutrino detection with IceCube

The IceCube Neutrino Observatory is currently the largest and most sensitive detector for astrophysical neutrinos and has pioneered the field of high-energy neutrino astronomy. Despite being designed with the primary goal of identifying astrophysical TeV neutrinos and their corresponding sources, recent studies, utilising the DeepCore subdetector, have shown IceCube's proficiency in being sensitive to astrophysical neutrinos at GeV energies. Currently, there is a gap in sensitivity between the supernova detection system at MeV energies and the lowest-energy triggering events around 1 GeV. In this contribution, we present the ongoing efforts to cover this gap and increase the sensitivity of IceCube to sub-GeV astrophysical neutrinos. Despite high background rates, we show how the complimentary use of manifold and supervised machine learning can make IceCube sensitive to neutrinos from transient sources down to energies of 100 MeV.


[67] 2507.08573

Measurement of the cosmic-ray energy spectrum above 2.5 EeV using 19 years of operation of the Pierre Auger Observatory

We present the spectrum of cosmic rays with energies above 2.5 EeV measured at the Pierre Auger Observatory after 19 years of operation, covering the period before the AugerPrime upgrade. Two independent event sets from the surface array of 1500 m-spaced detectors are combined, yielding a total exposure of approximately 100,000 km$^2$ sr yr. The first set includes events with zenith angles less than 60$^\circ$, while the second consists of events between 60$^\circ$ and 80$^\circ$, for which azimuthal asymmetries must be accounted for in the energy estimator. The threshold energy is chosen to ensure a trigger efficiency of the surface detector greater than 97%, thus minimizing composition biases. The energy scale is determined using high-quality fluorescence measurements, providing calorimetric estimates without reliance on simulations. A statistically successful combination is achieved within the uncorrelated systematic uncertainties of the individual spectra. All spectra are consistent when analyzing potential declination dependences, except for a mild modulation expected from the previously reported dipolar anisotropy. In particular, this statement applies to the northernmost declination band [+25$^\circ$,+45$^\circ$], where only contribute events with zenith angles between 60$^\circ$ and 80$^\circ$. Beyond the firmly established ankle and suppression spectral features, the combined spectrum across declinations $-90^\circ$ to +45$^\circ$ provides high-precision measurements of the instep feature with more than 5$\sigma$ confidence.


[68] 2507.08582

Observation of the Shadows of the Moon and Sun Using the Pierre Auger Observatory at an Average Energy of $7\times10^{17}\,$eV

The interaction of cosmic rays with celestial bodies such as the Moon or the Sun produces a shadow in the arrival direction distribution of the cosmic rays reaching the Earth. Such deficits from an isotropic flux have been observed by astroparticle observatories below energies of $10^{15}\,$eV. Above this energy, measurements were limited due to the low number of events as a result of the steeply falling cosmic-ray flux with energy. With more than 10.6 million events recorded during 20 years of operation of the Pierre Auger Observatory, we report the first observation of the shadow of the Moon at an average energy of $7\times10^{17}\,$eV with a maximum significance above 3$\sigma$. The shadow is an end-to-end check that the celestial directions are correctly reconstructed from the air shower data, and it is used here to derive the effective angular resolution for this dataset. Additionally, we present the results of a similar study on the shadow of the Sun.


[69] 2507.08583

Classical Fluid Analogies for Schrdinger-Newton Systems

The Schrödinger-Poisson formalism has found a number of applications in cosmology, particularly in describing the growth by gravitational instability of large-scale structure in a universe dominated by ultra-light scalar particles. Here we investigate the extent to which the behaviour of this and the more general case of a Schrödinger-Newton system, can be described in terms of classical fluid concepts such as viscosity and pressure. We also explore whether such systems can be described by a pseudo-Reynolds number as for classical viscous fluids. The conclusion we reach is that this is indeed possible, but with important restrictions to ensure physical consistency.


[70] 2507.08589

A new type-1 transition in NGC 1346 revealed by (spectro)polarimetry

We present the first photo- and spectropolarimetric observations of the changing-look active galactic nucleus NGC~1346 taken with the Perkins telescope (November 2022, March 2025) and VLT/FORS2 (August and September 2024), which have revealed a new spectral transition. We find that NGC~1346 has reverted to its former (2001) type-1 spectral state after having been a type-2 active galactic nucleus from 2004 to 2022. More and more intense broad H$\beta$ (3720 $\pm$ 310 km.s$^{-1}$ full width at half maximum) and H$\alpha$ (2985 $\pm$ 230 km.s$^{-1}$) emission lines appeared between August and September 2024, suggesting that the broad-line region (BLR) is increasingly irradiated. The strong blueshift of the Balmer lines ($\sim$ 1600 km.s$^{-1}$) indicates that the gas is moving toward the observer, possibly due to a radiation-driven outflow. The strong difference between 2022's ($\sim$ 3\%) and 2024-2025's de-biased polarization ($\le$ 0.5\%), the blueing of the spectra between August and September 2024, and the slow rise in the integrated flux between 2022 and 2025 argue against asymmetric or temporary obscuration by a cloud passing in front of the line of sight. Two options, either a re-illumination of the BLR or a tidal disruption event, can explain the observed properties of NGC~1346. Timely follow-up observations taken while the phenomenon is still ongoing are needed to determine which of the two solutions is the correct one.


[71] 2507.08592

XMM-Newton follow-up of two eROSITA X-ray binary candidates

We report on the follow-up observations with XMM-Newton of two X-ray binary candidates identified in the first eROSITA all-sky survey data (eRASS1), 1eRASS J061330.8+160440 and 1eRASS J161201.9-464622. Based on the obtained results, in particular, the observed X-ray spectra and lack of pulsations, as well as properties of the identified optical counterparts, we conclude that both candidates are unlikely to be XRBs. Based on LAMOST optical spectroscopy and SED fit results for 1eRASS J061330.8+160440 we classify it as an M0 chromospherically active subgiant star. ZTF and TESS photometry reveal highly significant period for this object of 7.189 days, which likely attributed to starspot(s). On the other hand, SALT follow-up spectroscopy of 1eRASS J161201.9-464622 solidly classifies this source as a bright novalike cataclysmic variable (CV), the second discovered with eROSITA. A persistent 4.802 h signal is found across all three available TESS observations, and is tentatively identified as the orbital period of the binary. Follow-up high-speed photometry and time-resolved spectroscopy are required to confirm the derived orbital modulation.


[72] 2507.08593

Spectra for the Vacuum Cherenkov Effect in Astrophysical Electromagnetic Cascades with Lorentz Invariance Violation

Lorentz invariance violation is a feature of several quantum gravity models in which Lorentz symmetry is broken at high energies, possibly leading to changes in particle behavior and interactions. In this work, we investigate vacuum Cherenkov radiation, a reaction in which an electron spontaneously emits a photon. This process, forbidden when Lorentz symmetry is unbroken, is a phenomenological consequence of some quantum gravity models. We derive, for the first time, the spectra for the vacuum Cherenkov reaction, and confirm our results numerically. These results can be used to derive limits on Lorentz invariance violation.


[73] 2507.08604

UNIONS Optical Identifications for VLASS Radio Sources in the Euclid Sky (UNVEIL) I. A Catalog of ~146,000 Radio Galaxies up to z~5

We present results of optical identifications for Very Large Array Sky Survey (VLASS) radio sources using the Ultraviolet Near Infrared Optical Northern Survey (UNIONS). A cross-match between UNIONS and VLASS Epoch 2 catalogs yields 146,212 radio galaxies down to $r=24.5$ mag over a wide area of $\sim4,200\ \mathrm{deg^2}$. We perform g-dropout selections and find $>200$ sources at $z\sim4$ optimistically. Of 63,019 sources with valid photo-$z$, 8,692 are at $z_\mathrm{photo}>1$, and 1,171 are at $z_\mathrm{photo}>2$. Based on spectral luminosities at 1.4 GHz using the valid photo-$z$, we identify $\sim49,000$ radio-loud AGNs (RLAGNs) with $\textit{L}_\mathrm{1.4GHz}>10^{24}$ W/Hz, and all radio galaxies at $z>1$ are RL. Adopting radio loudness instead, 138,266 out of 146,212 UNIONS-VLASS radio galaxies are RL. Thus, our catalog greatly increases the number counts of RLAGNs at $z>1$. We further cross-match the UNIONS-VLASS catalog with LOw-Frequency ARray Two-metre Sky Survey (LoTSS) at 144 MHz and Faint Images of the Radio Sky at Twenty-cm (FIRST) at 1.4 GHz, yielding 101,671 UNIONS-VLASS-LoTSS, 79,638 -FIRST, and 64,672 -LoTSS-FIRST sources, respectively. This multifrequency radio dataset reveals sources of various spectral shapes, including the steep spectrum of aged populations, the peaked spectrum of young populations, and the upturned spectrum potentially associated with transient sources. The UNIONS-VLASS radio galaxies will be covered by the Euclid wide survey, offering legacy values to benefit multi-faceted studies related to RLAGNs and beyond.


[74] 2507.08612

Nature or Nurture: LMC-like Dust in the Solar Metallicity Galaxy M31

Using the {\it Hubble Space Telescope}/Space Telescope Imaging Spectrograph, ultraviolet (UV) extinction curves have been measured in M31 along thirteen new sightlines, increasing the M31 sample to seventeen. This sample covers a wide area of M31 having galactocentric distances of 5 to 16 kpc, enabling the analysis of UV extinction curve variations over a large region of an external galaxy similar to the Milky Way with global galactic characteristics such as metallicity for the first time. No correlation is found between the extinction parameters and galactocentric distance which might be expected if there is a radial metallicity gradient in M31. Most of the new UV extinction curves presented here are significantly different from the average extinction curves of the Milky Way, LMC, and SMC, but the average M31 extinction curve is similar to the average extinction curve in the 30-Dor region of the LMC. The wide range of extinction curves seen in each individual Local Group galaxy suggests that global galactic properties such as metallicity may be less important than the local environmental conditions such as density, UV radiation field, and shocks along each sightline. The combined behavior of the Milky Way, LMC, SMC, and now M31 UV extinction curves supports the idea that there is a family of curves in the Local Group with overlapping dust grain properties between different galaxies.


[75] 2507.08614

First images of Antares photosphere from spectropolarimetry

Antares is the closest red supergiant (RSG) to Earth. The discovery of linear polarization in the atomic lines of the star opens the path to produce direct images of the photosphere, hence probing the dynamics at the surface. We analyze this linear polarization signals following the same scheme as has been previously done for the RSG Betelgeuse, and find that they are comparable in all its details. This allows us to use the same models for the analysis of these polarization signals in both stars. We found that as in Betelgeuse, the linear polarization signal of Antares is due to the depolarization of the continuum combined with brightness inhomogeneities. This allows us to produce images of the photosphere of star. We show that in Antares, convective cells can last several months and occupy roughly 30\% of the stellar radius for the largest ones.


[76] 2507.08633

Metal-THINGS: gas metallicity gradients in nearby galaxies

This paper explores the gas metallicity gradients in a sample of 25 nearby galaxies using new Integral Field Spectroscopy observations from the Metal-THINGS survey. We derive and study the resolved diffuse ionised gas content, Baldwin, Phillips and Terlevich diagrams and gas metallicities for our entire sample, at spatial resolutions of 40-300 pc. Gas metallicity gradients are studied as a function of the galaxy's stellar mass, H I gas fraction, diffuse ionised gas content, and using different parametric length scales for normalisation. The metallicity gradients are analysed using Bayesian statistics based on data from the Metal-THINGS survey. Bayesian MCMC models are developed to explore how metallicity gradients vary with a galaxy's mass and how they correlate with properties such as the stellar mass or the atomic gas fraction. For our sample, we find that the metallicity typically decreases with galactic radius, consistent with inside-out galaxy growth. We find a trend dependent on the stellar mass, with a break at log(M_star/M_sun)=9.5, and another between the metallicity gradients and the atomic gas fraction (f_g,HI) of a galaxy at fg,HI=0.75, indicating relatively shallower gradients for lower gas fractions. We find that normalisation using NUV-band effective radii are preferable for galaxies with a higher atomic gas content and lower stellar masses, while r-band radii are better suited for those with lower atomic gas fractions and more massive ones. Our results highlight a strong connection between gas content, stellar mass, and metallicity gradients. The breaks at log(M_star/M_sun)=9.5 and fg,HI=0.75 mark shifts in chemical enrichment behaviour, with low-mass galaxies showing greater sensitivity to gas processes. Overall, this points to gas accretion and removal as key drivers of chemical evolution in low-mass systems.


[77] 2507.08643

Precision spectral measurements of Chromium and Titanium from 10 to 250 GeV$/n$ and sub-Iron to Iron ratio with the Calorimetric Electron Telescope on the ISS

The Calorimetric Electron Telescope (CALET), in operation on the International Space Station since 2015, collected a large sample of cosmic-ray (CR) iron and sub-iron events over a wide energy interval. In this paper we report an update of our previous measurement of the iron flux and we present - for the first time - a high statistics measurement of the spectra of two sub-iron elements Cr and Ti in the energy interval from 10 to 250 GeV/n. The analyses are based on 8 years of data. Differently from older generations of cosmic-ray instruments which, in most cases, could not resolve individual sub-iron elements, CALET can identify each nuclear species from proton to nickel (and beyond) with a measurement of their electric charge. Thanks to the improvement in statistics and a more refined assessment of systematic uncertainties, the iron spectral shape is better resolved, at high energy, than in our previous paper and we report its flux ratio to chromium and titanium. The measured fluxes of Cr and Ti show energy dependences compatible with a single power law with spectral indices $-2.74 \pm 0.06$ and $-2.88 \pm 0.06$, respectively.


[78] 2507.08647

Growth of Cosmic Structures in generalized mass-to-horizon relation Entropic Cosmology

We investigate the growth of cosmic structures in the thermodynamically consistent generalised mass-to-horizon entropic cosmology (MHEC). For the Bekenstein case the entropic energy density augments the Friedmann equations without modifying the Hawking temperature and automatically satisfies the Clausius relation, thereby avoiding the inconsistencies that afflicted earlier entropy-force models. We then derive the linear perturbation equations, emphasising the distinction between a fully perturbed interaction term and the common approximation in which the perturbation is neglected. Numerical solutions show that fully perturbed follows the LCDM matter-growth history within the current growth uncertainties. Our results demonstrate that MHEC matches both background and growth probes as well as LCDM without extra free parameters, providing a viable entropic explanation for recent accelerated expansion of the Universe.


[79] 2507.08651

Chemical abundances and small-scale gas kinematics in the local star-forming galaxy NGC 2366

We present a detailed spectroscopic study of ionised gas in the nearby (~3.3Mpc) dwarf galaxy NGC 2366, a local analogue of Green Pea galaxies, based on observations with the SCORPIO-2 instrument at the Russian 6-m BTA telescope. Using scanning Fabry-Perot interferometry and long-slit spectroscopy, we examine the gas kinematics, excitation mechanisms, and chemical abundances across the disc of NGC 2366, including its prominent starburst region Mrk 71 and the companion region NGC 2363. We identified 20 regions with locally elevated Ha velocity dispersion, only four of which correspond to known high-energy sources. We argue that one of the remaining objects can be a previously unidentified Wolf-Rayet star and two - supernova remnants. For 15 HII regions, we measure electron temperatures, oxygen and nitrogen abundances via the `direct' Te method, with 12 + log(O/H) ranging from 7.6 to 8.0 in most of the regions. We show that Mrk 71 has higher oxygen abundance compared to the other \HII regions in the galaxy, contrary to the previous indirect estimates suggesting flat gradient throughout the galaxy. Together with the localized spatial variations of metallicity in the area, it is indicative of metal enrichment by the outflow from the super star cluster in the centre of Mrk 71.


[80] 2507.08654

Spectral analysis of two directly imaged benchmark L dwarf companions at the stellar-substellar boundary

We used multiple epochs of high-contrast imaging spectrophotometric observations to determine the atmospheric characteristics and thermal evolution of two previously detected benchmark L dwarf companions, HD 112863 B and HD 206505 B. We analyzed IRDIS and IFS data from VLT/SPHERE of each companion, both of which have dynamical masses near the stellar-substellar boundary. We compared each companion with empirical spectral standards, as well as constrained their physical properties through atmospheric model fits. From these analyses, we estimate that HD 112863 B is spectral type $\rm{L}3\pm1$ and that HD 206505 B is spectral type $\rm{L}2\pm1$. Using the BT-Settl atmospheric model grids, we find a bimodal solution for the atmospheric model fit of HD 112863 B, such that $T_{\rm{eff}}=1757^{+37}_{-36}$ K or $2002^{+23}_{-24}$ K and $\log{g}=4.973^{+0.057}_{-0.063}$ or $5.253^{+0.037}_{-0.033}$, while for HD 206505 B, $T_{\rm{eff}}=1754^{+13}_{-13}$ K and $\log{g}=4.919^{+0.031}_{-0.029}$. Comparing the bolometric luminosities of both companions with evolutionary models imply that both companions are likely above the hydrogen burning limit.


[81] 2507.08666

The IceCube Collaboration -- Contributions to the 39th International Cosmic Ray Conference (ICRC2025)

The IceCube Observatory at the South Pole has been operating in its full configuration since May 2011 with a duty cycle of about 99%. Its main component consists of a cubic-kilometer array of optical sensors deployed deep in the Glacial ice designed for the detection of high-energy astrophysical neutrinos. A surface array for cosmic ray air shower detection, IceTop, and a denser inner subdetector, DeepCore, significantly enhance the capabilities of the observatory, making it a multipurpose facility. This list of contributions to the 39th International Cosmic Ray Conference in Geneva, Switzerland (July 15-24, 2025) summarizes the latest results from IceCube covering a broad set of key questions in physics and astrophysics. The papers in this index are grouped topically to highlight IceCube contributions related to high-energy neutrino and multi-messenger astrophysics, atmospheric fluxes, cosmic-ray physics, low-energy neutrino transients, physics beyond the Standard Model, detector calibration and event reconstruction, and the status and performance of the IceCube Upgrade, a dense sensor infill complemented by calibration devices to be deployed by the end of 2025. Contributions related to IceCube-Gen2, the planned future extension of IceCube, are available in a separate collection.


[82] 2507.08667

The IceCube-Gen2 Collaboration -- Contributions to the 39th International Cosmic Ray Conference (ICRC2025)

IceCube-Gen2 is a planned next-generation neutrino observatory at the South Pole that builds upon the successful design of IceCube. Integrating two complementary detection technologies for neutrinos, optical and radio Cherenkov emission, in combination with a surface array for cosmic-ray air shower detection, IceCube-Gen2 will cover a broad neutrino energy range from MeV to EeV. This index of contributions to the 39th International Cosmic Ray Conference in Geneva, Switzerland (July 15-24, 2025) describes research and development efforts for IceCube-Gen2. Included are summaries of the design, status, and sensitivity of the IceCube-Gen2 optical, surface, and radio components; performance studies of next-generation surface detectors and in-ice optical sensors; advanced reconstruction techniques of cosmic-ray air showers and neutrino events; sustainability and environmental impact; and sensitivity studies of astrophysical neutrino fluxes and cosmic-ray physics. Contributions related to IceCube and the scheduled IceCube Upgrade are available in a separate collection.


[83] 2507.08680

The role of dissipation distance on reconnection-driven multi-messenger signals from blazar jets

Blazars are characterized by relativistic jets that are closely aligned with our line of sight. This results in relativistic beaming, making blazars among the most luminous extragalactic sources across the electromagnetic spectrum, from radio waves to gamma-rays and, potentially, in high-energy neutrinos. We present a comprehensive study of multi-messenger emission from blazar jets powered by magnetic reconnection occurring at varying distances from the supermassive black hole (SMBH). By generalizing previous models, we explore how the emission characteristics depend self-consistently on the spatial evolution of key jet properties, including magnetization, bulk Lorentz factor, and external photon fields (accretion disc, broad-line region, and dusty torus). Using numerical simulations, we examined the impact of the initial jet magnetization, particle acceleration efficiency, jet-to-accretion power ratio, and mass accretion rate on the broadband photon spectra and neutrino emission. Our findings reveal distinct emission regimes characterized by different dominant radiative processes: synchrotron and synchrotron self-Compton dominate closer to the SMBH where magnetization is high, while external Compton (EC) processes become significant near the broad-line region (BLR). Neutrino production efficiency is highest upstream of the BLR, driven by enhanced photon target densities from synchrotron and external photons available for photopion interactions, whereas the proton particle distribution is hard. Our model predictions are compared with observations of gamma-ray luminosities and synchrotron peak energies of Fermi-detected blazars, highlighting magnetic reconnection as a potential mechanism driving both electromagnetic and neutrino emissions in astrophysical jets.


[84] 2507.08687

Interpreting Hubble tension with a cascade decaying dark matter sector

Hubble tension can be alleviated by altering either early- or late-time $\Lambda$CDM. So far, extensive studies have shown that only early dark energy or ad hoc combinations of those two-fold effects can reduce the tension to $3\sigma$ or lower. In this work, we improve the later solution by considering a cascade decaying dark matter sector, where a parent particle produces relativistic dark matter in the early Universe and the dark matter subsequently decays to neutrinos in the late-time Universe. We parametrize this model with four model parameters, carry out a Markov Chain Monte Carlo fit to Planck 2018+BAO+Pantheon data sets, and compare it to Big Bang Nucleosynthesis limits, neutrino flux bounds, and Planck and SDSS data about matter power spectra. Within parameter regions reducing the Hubble tension to $3.2\sigma$ and compatible with the existing constraints, the main phenomenological features include a larger $S_{8}\sim 0.84$, smaller $n_{s}\sim 0.95$ and larger matter power spectra relative to the $\Lambda$CDM, which are left for future tests.


[85] 2507.08691

Advances in reconstructing the cosmic-ray energy spectrum with IceTop

The IceTop array at the surface of the IceCube Neutrino Observatory measures extensive air showers produced by cosmic-ray particles with energies from PeV up to EeV, covering the transition region from galactic to extragalactic sources. This contribution presents significant improvements that will enhance the measurement of the cosmic-ray energy spectrum. (I) To analyze more than a decade of data with increasing snow overburdens on the detector, an improved method to handle the time-dependent attenuation of the detector signals was developed. (II) New analysis cuts have been developed to increase the measured event rate while improving the reconstruction quality. (III)~A~new reconstruction method that separately fits the electromagnetic and muonic components, while incorporating data from the deep in-ice detector, enables the reconstruction of air showers with cores outside the IceTop array, up to zenith angles of 50°. These improvements will significantly increase the event statistics and extend the IceTop spectrum towards higher energies.


[86] 2507.08695

The optimal Padé polynomial for reconstruction of luminosity distance based on 10-fold cross-validation

The cosmography known as the Padé polynomials has been widely used in the reconstruction of luminosity distance, and the orders of Padé polynomials influence the reconstructed result derived from Padé approximation. In this paper, we present a more general scheme of selecting optimal Padé polynomial for reconstruction of luminosity distance based on 10-fold cross-validation. Then the proposed scheme is applied to Pantheon+ dataset. The numerical results clearly indicate that the proposed procedure has a remarkable ability to distinguish Padé approximations with different orders for the reconstruction of the luminosity distance. We conclude that the (2,1) Padé approximation is the optimal approach that can well explain Pantheon+ data at low and high red-shifts. Future applications of this scheme could help choose the optimal model that is more suitable for cosmological observation data at hand and gain a deeper understanding of the universe.


[87] 2507.08706

IceCube searches for GeV neutrino counterparts associated with high-energy starting events

The origin of the astrophysical neutrino flux observed by IceCube is largely unknown. To help decipher its astrophysical origin, we propose an IceCube analysis that conducts follow-up searches for GeV neutrinos associated with neutrino events above 60 TeV, which are known to have a high probability to be of astrophysical origin. It could not only identify a new component of the astrophysical neutrino flux, but also characterize how its spectrum extrapolates from GeV to PeV energies. This would in turn give valuable insights into the internal processes of neutrino sources. Astrophysical transients, such as collapsars, have been proposed as sources of time-correlated GeV- and high-energy neutrinos. Conducting this search in such short time scales allows for a substantial reduction in the dominant background rate for GeV neutrino candidate events. We introduce the statistical method and sensitivity of this search as well as dedicated data quality checks. None of the searches yield statistically significant results, and we present the first limits on GeV neutrino emission associated to VHE neutrino events at short time scales.


[88] 2507.08709

Unambiguous GeV gamma-rays from molecular clouds illuminated by particles diffusing from the adjacent supernova remnant G335.2+0.1 that is confined in an expanding bubble

We report the detection of GeV gamma-ray emission likely associated with supernova remnant (SNR) G335.2+0.1 and the finding of a molecular cloud ($\sim20$--$30^\prime$ in angular size) that is very likely in physical contact with the SNR and responsible for the gamma-ray emission. Using the 16.8 yr Fermi-LAT data, an extended emission, with a significance of 13.5 $\sigma$ and a radius 0.24° in 0.2--500 GeV in the uniform-disk model, was found to the adjacent east of the SNR. With archival Mopra CO-line data, a large molecular clump at local-standard-of-rest velocity $\sim-48$ to $-43$ km s$^{-1}$ was revealed appearing coincident with the gamma-ray source. The SNR was found located in a cavity encircled by a 'C'-shaped ring-like molecular shell at $-45$ to $-43$ km s$^{-1}$. This morphological agreement, together with the position-velociy diagrams made along lines cutting across the cavity, suggests that the SNR was evolving in the expanding molecular bubble created by the stellar wind of the progenitor with a mass $\gtrsim 15 M_{\mathrm{sun}}$. The giant molecular cloud, visible at around $-46$ km s$^{-1}$, and the associated SNR are thus estimated to lie at a kinematic distance of 3.1 kpc, with the HI absorption taken into account. We suggest that the SNR has entered the radiative phase after the blastwave recently struck the cavity wall. With the evolutionary scenario of the SNR, we demonstrate that the gamma-ray emission reported here can be naturally interpreted by the hadronic interaction between the accelerated protons that escaped from the SNR shock and the eastern large molecular clump.


[89] 2507.08723

Monitoring and performance of AugerPrime

The upgrade of Pierre Auger Observatory, AugerPrime, is a multi-hybrid system designed to improve the sensitivity and precision of ultra-high-energy cosmic ray measurements. It includes scintillator detectors positioned both atop the enhanced Water-Cherenkov detectors and buried nearby for direct muon measurements, along with radio and fluorescence detectors. In this contribution, we present an overview of the monitoring tools developed for all the components of AugerPrime, focusing on real-time performance assessment and long-term stability metrics. By continuously tracking key parameters, we can identify potential issues early, enabling timely interventions and improving overall data quality. These strategies are crucial for maintaining the long-term reliability of the measurements taken at the Auger Observatory and providing high-quality data for cosmic ray research in the coming decades.


[90] 2507.08727

Constraining the survival of HCN during cometary impacts

Cometary impacts have been invoked as an atmosphere-independent method of stockpiling hydrogen cyanide (HCN), a key prebiotic feedstock molecule, into environments favourable for the onset of prebiotic chemistry on the early Earth. This work revisits the prospects for cometary delivery of HCN through new impacts simulations of idealised cometary bodies using the shock physics code iSALE combined with simple chemical modelling. Using temperature and pressure profiles for material within spherical, non-porous comets with a high resolution of Lagrangian tracer particles, we assess the survival rate of HCN across a range of impact velocities, sizes and angles, assuming both steady state and equilibrium chemistry. We find that HCN survival is extremely limited at impact velocities above the escape velocity of the Earth, unless the impact occurs at extreme obliquity ($\theta \sim 15^\circ$). We present a parametrisation of the survival of HCN as a function of impact velocity, angle, and cometary diameter, which provides an upper limit to survival in more realistic scenarios to aid with future studies investigating the role of comets in the origins of life. Although successful HCN delivery may be possible in our idealised model, we neglect to consider the effect of atmospheric passage and our results suggest that delivery alone is not likely to be sufficient for the onset of prebiotic chemistry.


[91] 2507.08731

Statistical Analysis of Early Spectra in Type II and IIb Supernovae

We present a comprehensive analysis of the early spectra of type II and type IIb supernovae (SNe) to explore their diversity and distinguishable characteristics. Using 866 publicly available spectra from 393 SNe, 407 from type IIb SNe (SNe IIb) and 459 from type II SNe (SNe II), we analysed H$\alpha$ and He~I 5876 A at early phases ($<40$ days from the explosion) to identify possible differences between these two SN types. By comparing the pseudo-equivalent width (pEW) and full width at half maximum (FWHM), we find that the strength of the absorption component of these lines serves as a quantitative discriminator, with SNe IIb exhibiting stronger lines at all times. The most significant differences emerge within the first 10-20 days. To assess the statistical significance of these differences, we apply statistical methods and machine-learning techniques. Population density evolution reveals a clear distinction in both pEW and FWHM. Quadratic Discriminant Analysis confirms distinct evolutionary patterns, particularly in pEW, while FWHM variations are less pronounced. A combination of t-distributed Stochastic Neighbour Embedding and Linear Discriminant Analysis effectively separates the two SN types. Additionally, a Random Forest Classifier demonstrates the robustness of pEW and FWHM as classification criteria, allowing for accurate classification of newly observed SNe II and IIb based on computed classification probabilities. Applying our method to low-resolution spectra obtained from the Zwicky Transient Facility Bright Transient, a magnitude-limited survey, we identified 34 misclassified SNe. This revision increases the estimated fraction of SNe IIb from 4.0% to 7.26%. This finding suggests that misclassification significantly impacts the estimated core-collapse SN rate. Our approach enhances classification accuracy and provides a valuable tool for future supernova studies.


[92] 2507.08747

Performance and extensions of the central data acquisition system for Phase II of the Pierre Auger Observatory

The Pierre Auger Observatory is a hybrid detector designed to observe and study ultra-high-energy particles of extraterrestrial origin. With its 27 fluorescence telescopes and over 1600 autonomously operating water-Cherenkov detectors spread over an area of 3000 km$^2$, it is world-leading in terms of exposure to cosmic rays and offers an unparalleled window into the physical processes that happen at energy scales unattainable by particle accelerators on Earth. Measurement information of candidate air-shower events from all associated detectors and telescopes is collected at a central data-acquisition system located in the nearby town of Malargüe, and processed for higher-level physics analysis. On top of this, data for monitoring the long-term stability and operation of the observatory is forwarded to the central server as well. In this work, we briefly review the central data-acquisition system of the Pierre Auger Observatory. We examine the rates, efficiencies, and purity of detected events for Phase II of the Pierre Auger Observatory and compare them to performance parameters during Phase I. We detail challenges in the event detection up until now and present recent changes in the central data acquisition system and local station software that aim to streamline the data acquisition chain.


[93] 2507.08748

Low-latency neutrino follow-up combining diverse IceCube selections

Neutrino observations are crucial for multi-messenger astronomy, but currently limited by effective area and atmospheric background. However, while other telescopes must point their limited field of view, IceCube can perform full-sky realtime follow-up of astrophysical transient sources with high uptime. For this, IceCube uses its Fast Response Analysis (FRA), which can provide results from an unbinned maximum likelihood analysis within tens of minutes of an astrophysical transient. Besides manually selected transient candidates, it also routinely scans areas of the sky compatible with gravitational wave alerts from LIGO/Virgo/KAGRA and the IceCube event singlets most likely to originate from an astrophysical source. Currently the analysis uses TeV muon neutrino candidate events whose signature permits especially precise angular reconstruction, selected and reconstructed at the South Pole and transmitted with low-latency over satellite. Recently, different event selections are also being included in IceCube analyses. These efforts include the follow-up of gravitational wave events with GeV neutrinos detected by IceCube-DeepCore and the observation of the Galactic plane with cascade events produced by all neutrino flavors. If made available on a day-scale latency, these event selections can also be used in FRAs. Moreso, multiple event samples can be combined in a Fast Response Analysis that is sensitive to a broader energy range of a neutrino transient spectrum and ensures the inclusion of all neutrino flavors. We present the analysis method and technical aspects of such an extension of the existing framework. This includes a proposed new pipeline allowing the inclusion of the more computationally-intensive reconstruction methods used by the aforementioned event selections. The extension is validated using example analyses implemented in this framework. (abbreviated)


[94] 2507.08753

Measuring the Astrophysical Galactic Plane Neutrino Flux and Searching for Galactic PeVatrons using the IceCube Multi-Flavor Astrophysical Neutrino Sample

The IceCube Neutrino Observatory has provided new insights into the high-energy universe, in particular, unveiling neutrinos from the galactic plane. However, galactic neutrino sources are still unresolved. The recent detection of multi-PeV photons by LHAASO from the Cygnus region highlights its potential as a galactic neutrino source. Additionally, LHAASO, HAWC, and HESS have reported over forty galactic gamma-ray sources with energies above 100 TeV. Detecting neutrinos correlated with high-energy gamma-ray sources would provide compelling evidence of cosmic-ray acceleration in these galactic sources. In this work, we compile a 12.3-year, full-sky, all-flavor dataset, the IceCube Multi-Flavor Astrophysics Neutrino sample (ICEMAN). ICEMAN is the combination of three largely independent neutrino samples of different event morphologies and builds upon the previous work of the DNN-based cascade sample, Enhanced Starting Track Event Selection, and the Northern Track sample. Recent improvements in ice modeling and detector calibration are also incorporated into the cascade reconstruction. In addition to revisiting the galactic plane, we adopt two different analysis methods to search for galactic PeVatrons. First, we use a template-based approach to probe the Cygnus Cocoon region. Second, we use a point source hypothesis to find correlations between IceCube neutrinos and gamma-ray sources detected at energies greater than 100 TeV.


[95] 2507.08760

Recent Advances in Understanding R-Process Nucleosynthesis in Metal-Poor Stars and Stellar Systems

The rapid neutron-capture process (r-process) is responsible for the creation of roughly half of the elements heavier than iron, including precious metals like silver, gold, and platinum, as well as radioactive elements such as thorium and uranium. Despite its importance, the nature of the astrophysical sites where the r-process occurs, and the detailed mechanisms of its formation, remain elusive. The key to resolving these mysteries lies in the study of chemical signatures preserved in ancient, metal-poor stars. In this review, we explore r-process nucleosynthesis, focusing on the sites, progenitors, and formation mechanisms. We discuss the role of potential astrophysical sites such as neutron star mergers, core-collapse supernovae, magneto-rotational supernovae, and collapsars, that can play a key role in producing the heavy elements. We also highlight the importance of studying these signatures through high-resolution spectroscopic surveys, stellar archaeology, and multi-messenger astronomy. Recent advancements, such as the gravitational wave event GW170817 and detection of the r-process in the ejecta of its associated kilonovae, have established neutron star mergers as one of the confirmed sites. However, questions remain regarding whether they are the only sites that could have contributed in early epochs or if additional sources are needed to explain the signatures of r-process found in the oldest stars. Additionally, there are strong indications pointing towards additional sources of r-process-rich nuclei in the context of Galactic evolutionary timescales. This review summarizes what has been learned so far, the challenges that remain, and the exciting prospects for future discoveries. The increasing synergy between observational facilities, computational models, and large-scale surveys is poised to transform our understanding of r-process nucleosynthesis in the coming years.


[96] 2507.08775

All-flavor Time-dependent Search for Transient Neutrino Sources

Transient sources are among the preferred candidates for the sources of high-energy neutrino emission. Intriguing examples so far include blazar flares and tidal disruption events coincident with IceCube neutrinos. Here, we report the first all-flavor, all-sky time-dependent search for neutrino sources by combining IceCube throughgoing tracks, starting tracks and cascades. Throughgoing tracks provide the best sensitivity in the Northern Sky, while cascades have worse angular resolution but yield better sensitivity in the Southern Sky than tracks. The relatively new starting tracks sample has reduced contamination from atmospheric muons. This analysis takes advantage of the strengths of each of the datasets, combining them for increased statistics and obtaining the best accessible all-sky sensitivity for transient searches. In this search, we look for unbound $E^{-\gamma}$ power-law sources, as well as $E^{-2}$ sources with low and high-energy exponential cutoffs, optimizing the sensitivity for the duration of the flares.


[97] 2507.08778

Mining the Alerts: A Preliminary Catalog of Compact Binaries from the Fourth Observing Run

We present a preliminary catalog of compact binary merger candidates from the ongoing fourth observing run (O4) of Advanced LIGO, Virgo, and KAGRA, based on an analysis of public alerts distributed through GraceDB as of May 2025. We developed and applied methods to estimate the source-frame chirp mass for each candidate by utilizing information from public data products, including source classification probabilities, sky localizations, and observatory status. Combining our O4 analysis with previous catalogs, we provide updated estimates for the local merger rate density. For sources with chirp mass characteristic of binary neutron stars ($[1, 1.5]\,M_\odot$), we find a rate of $56^{+99}_{-40}$ $\textrm{Gpc}^{-3}\,\textrm{yr}^{-1}$. For systems in the expected neutron star--black hole chirp mass range ($[1.5, 3.5]\,M_\odot$), the rate is $36^{+32}_{-20}$ $\textrm{Gpc}^{-3}\,\textrm{yr}^{-1}$, and for heavier binary black holes ($[3.5, 100]\,M_\odot$), we estimate a rate of $19^{+4}_{-2}$ $\textrm{Gpc}^{-3}\,\textrm{yr}^{-1}$. This work provides an early glimpse into the compact binary population being observed in O4; we identify a number of high-value candidates up to signal-to-noise $\sim 80$, which we expect to enable precision measurements in the future.


[98] 2507.08779

Advancements in the IceAct Energy Spectrum Analysis

The IceAct telescopes are Imaging Air Cherenkov telescopes installed as part of the IceCube Neutrino Observatory at the geographic South Pole. They consist of a 61 pixel camera and are small and robust to withstand the harsh environmental conditions. IceAct detects Cherenkov light produced by cosmic-ray particles with energies above approximately 10\,TeV interacting inside the atmosphere, which is complementary to the measurement of the air shower at the surface by IceTop and the high-energy muons in the deep ice. Two telescopes have been taking data since 2019 with a conservative estimated duty cycle of around 10\%. A graph neural network is used to reconstruct the basic air shower properties, like geometry and primary energy. This work focuses on the current progress in analyzing the energy spectrum of cosmic rays using IceAct data.


[99] 2507.08786

Non-relativistic effective theories for fields with general potentials and their implications for cosmology

Non-relativistic effective field theories (NREFTs) play a crucial role in various areas of physics, from cold atom experiments to cosmology. In this paper, we present a systematic framework for deriving NREFTs from relativistic theories with generic self-interactions. Our approach allows for (but is not limited to) non-power-law potentials (such as those arising from dilatons or axions) or potentials that are non-analytic around the classical vacuum (such as those with logarithmic radiative corrections). These are of theoretical and phenomenological interest but have largely been unexplored in the non-relativistic regime. NREFTs are typically viewed as approximations for systems with low velocities, weak couplings, and small field amplitudes. The latter assumption is relaxed in our approach, as long as the mass term remains dominant (ensuring the validity of the non-relativistic limit). Additionally, we establish an effective fluid description for the non-relativistic scalar field, identifying key quantities such as energy density, pressure, and sound speed. To enable cosmological applications, we extend our formalism to account for the expanding universe, providing a reliable tool for investigating ultra-light dark matter models with arbitrary self-interactions. Finally, we demonstrate the applicability of our NREFT in analyzing solitons, which is also relevant to cosmology for studying celestial objects such as boson stars and the cores of dark matter halos.


[100] 2507.08787

The THESAN-ZOOM project: Mystery N/O more -- uncovering the origin of peculiar chemical abundances and a not-so-fundamental metallicity relation at $3<z<12$

We present an analysis of metallicities and chemical abundances at $3-0.6$) without the need for exotic yields in our chemical network. Instead, bursty star formation naturally generates order-of-magnitude excursions in N/O on $\lesssim$100 Myr timescales due to temporally differential galactic winds; after a starburst, stellar feedback expels gas, leaving a large population of asymptotic-giant-branch stars to dominate the enrichment of the relatively low-mass interstellar medium. NRGs lie below the main sequence and typically exhibit $\mathrm{EW}[H$\beta$]\lesssim40$ Å, in apparent tension with observed high-EW NRGs. This tension is reconciled if observed NRGs are in the initial stages of a subsequent starburst, illuminating previously enriched gas, which is supported by the finding of high SFR surface density nitrogen-rich giant molecular clouds.


[101] 2507.08790

Triggering and quenching in the shadow of AGN: How does AGN proximity affect star formation in the EAGLE simulation?

Active galactic nuclei (AGN) are among the most energetic phenomena in the universe, capable of regulating star formation in galaxies via radiative and mechanical feedback. While AGN feedback is well studied in host galaxies, its influence on neighbouring galaxies within the same large-scale environment remains less understood. In this work, we use the EAGLE cosmological hydrodynamical simulation to examine how proximity to AGN affects star formation in nearby star-forming galaxies (SFGs) out to 2 Mpc. A control sample matched in stellar mass and local density allows us to isolate AGN-driven environmental effects, quantified through the $\Delta$SFR offset defined as the difference between a galaxy's star formation rate (SFR) and that of its matched controls. We find a significant and mass-dependent impact where $61\%$ of AGN-adjacent SFGs show suppressed star formation, while $39\%$ exhibit enhancement. Although the overall magnitude of these offsets is modest, the trends are statistically robust and consistent across the population. Suppression is prevalent in massive, gas-poor galaxies in high-mass halos, consistent with thermal feedback inhibiting gas cooling. Enhanced SFRs appear in low-mass, gas-rich galaxies, suggesting AGN-driven compression may locally trigger star formation. These trends emerge despite EAGLE implementing AGN feedback as a single-mode, stochastic thermal process, indicating that such models can reproduce both quenching and triggering regimes. Our findings demonstrate that AGN feedback extends well beyond host galaxies, impacting neighbouring systems in a non-uniform, distance-, mass-, and gas-dependent manner. This non-local AGN influence represents a critical, yet often overlooked, mechanism in shaping galaxy evolution within large-scale structures, bridging high-energy astrophysical processes and cosmic-scale environmental regulation.


[102] 2507.08792

The Effective Field Theory of Large Scale Structure for Mixed Dark Matter Scenarios

We initiate a systematic study of the perturbative nonlinear dynamics of cosmological fluctuations in dark sectors comprising a fraction of non-cold dark matter, for example ultra-light axions or light thermal relics. These mixed dark matter scenarios exhibit suppressed growth of perturbations below a characteristic, cosmologically relevant, scale associated with the microscopic nature of the non-cold species. As a consequence, the scale-free nonlinear solutions developed for pure cold dark matter and for massive neutrinos do not, in general, apply. We thus extend the Effective Field Theory of Large Scale Structure to model the coupled fluctuations of the cold and non-cold dark matter components, describing the latter as a perfect fluid with finite sound speed at linear level. We provide new analytical solutions wherever possible and devise an accurate and computationally tractable prescription for the evaluation of the one-loop galaxy power spectrum, which can be applied to probe mixed dark matter scenarios with current and upcoming galaxy survey data. As a first application of this framework, we derive updated constraints on the energy density in ultra-light axions using a combination of Planck and BOSS data. Our refined theoretical modeling leads to somewhat weaker bounds compared to previous analyses.


[103] 2507.08797

Multi-frequency analysis of the ALMA and VLA high resolution continuum observations of the substructured disc around CI Tau. Preference for sub-mm-sized low-porosity amorphous carbonaceous grains

(Abridged) We present high angular resolution and sensitivity ALMA 3.1 mm and VLA 9.1 mm observations of the disc around CI Tau. These new data were combined with similar-resolution archival ALMA 0.9 and 1.3 mm observations and new and archival VLA 7.1 mm, 2.0, 3.0, and 6.0 cm photometry to study the properties of dust in this system. At wavelengths <3.1 mm, CI Tau's continuum emission is very extended and highly substructured (with three gaps, four rings, and two additional gap-ring pairs identified by non-parametric visibility modelling). Instead, the VLA 9.1 mm data are dominated by a bright central component, only partially (< 50%) due to dust emission, surrounded by a marginally detected, faint, and smooth halo. We fitted the ALMA and VLA 9.1 mm data together, adopting a physical model that accounts for the effects of dust absorption and scattering. For our fiducial dust composition ("Ricci" opacities), we retrieved a flat maximum grain size distribution across the disc radius of $(7.1\pm0.8)\times10^{-2}$ cm, that we tentatively attributed to fragmentation of fragile dust or bouncing. We tested, for the first time, the dependence of our results on the adopted dust composition model to assess which mixture can best reproduce the observations. We found that the "Ricci" opacities work better than the traditionally adopted "DSHARP" ones, while graphite-rich mixtures perform significantly worse. We also show that, for our fiducial composition, the data prefer low-porosity (< 70%) grains, in contrast with claims of highly porous aggregates in younger sources, which we tentatively justified by time-dependent compaction. Our results are in line with constraints from disc population synthesis models and naturally arise from CI Tau's peculiar spectral behaviour, making this disc an ideal target for deeper cm-wavelength and dust polarisation follow-ups.


[104] 2507.08798

The Atacama Cosmology Telescope: High-redshift measurement of structure growth from the cross-correlation of Quaia quasars and CMB lensing from ACT DR6 and $\textit{Planck}$ PR4

We measure the amplitude of matter fluctuations over a wide range of redshifts by combining CMB lensing observations from ACT DR6 and $\textit{Planck}$ PR4 with the overdensity of quasars from Quaia, a $\textit{Gaia}$ and $\textit{unWISE}$ quasar catalog. Our analysis includes the CMB lensing power spectrum from ACT DR6, the auto-correlation of two Quaia quasar samples centered at $z \simeq 1.0$ and $z \simeq 2.1$, and their cross-correlations with CMB lensing from both ACT DR6 and $\textit{Planck}$ PR4. By performing a series of contamination and systematic null tests, we find no evidence for contamination in the lensing maps, contrary to what was suggested in previous Quaia cross-correlation analyses using $\textit{Planck}$ PR4 CMB lensing data. From the joint analysis of the quasar auto- and cross-correlations with CMB lensing, and including BOSS BAO data to break the degeneracy between $\Omega_m$ and $\sigma_8$, we obtain $\sigma_8 = 0.802^{+0.045}_{-0.057}$, consistent with $\Lambda$CDM predictions from $\textit{Planck}$ primary CMB measurements. Combining the CMB lensing auto-spectrum with the cross-correlation measurement improves the constraint on $\sigma_8$ by $12\%$ relative to the lensing auto-spectrum alone, yielding $\sigma_8 = 0.804 \pm 0.013$. This dataset combination also enables a reconstruction of structure growth across redshifts. We infer a $12\%$ constraint on the amplitude of matter fluctuations at $z > 3$, with a measurement at the median redshift of the signal of $\sigma_8(\tilde{z}=5.1) = 0.146^{+0.021}_{-0.014}$, consistent with $\textit{Planck}$ at the $1.4\sigma$ level. These results provide one of the highest redshift constraints on the growth of structure to date.


[105] 2507.08054

Gas of thin tubes of massless scalar field as a possible model of unified dark fluid

In this work, we investigated the influence of constants that characterize the motion of thin tubes of a massless scalar field (TToMSF) in a gas on the global properties of such a gas, in particular, on its ability to accelerated expansion. The peculiarity of the TToMSF gas is related to the fact that the parameter of its equation of state $w$ depends on the degree of rarefaction of the gas and can take values in a very wide range of values $-1< w <1$. In this case, fluctuations in the density of the TToMSF gas, caused by its gravitational interaction with non-uniformly distributed baryonic matter, should lead to the implementation of local equations of state with different values and different laws of change of the parameter $w$. The possibility of implementing local equations of state can alleviate the tensions associated with the observed non-uniformity of the distribution of matter and the growth of structures in the Universe.


[106] 2507.08095

A model-agnostic gravitational-wave background characterization algorithm

As ground-based gravitational-wave (GW) detectors improve in sensitivity, gravitational-wave background (GWB) signals will progressively become detectable. Currently, searches for the GWB model the signal as a power law, however deviations from this model will be relevant at increased sensitivity. Therefore, to prepare for the range of potentially detectable GWB signals, we propose an interpolation model implemented through a transdimensional reverse-jump Markov chain Monte Carlo (RJMCMC) algorithm. This interpolation model foregoes a specific physics-informed model (of which there are a great many) in favor of a flexible model that can accurately recover a broad range of potential signals. In this paper, we employ this framework for an array of GWB applications. We present three dimensionless fractional GW energy density injections and recoveries as examples of the capabilities of this spline interpolation model. We further demonstrate how our model can be implemented for hierarchical GW analysis on $\Omega_{\rm GW}$.


[107] 2507.08116

Singularities in loop quantum cosmology

Quantum effects are expected to modify the cosmological dynamics of the early universe while maintaining some (potentially discrete) notion of space-time structure. In one approach, loop quantum cosmology, current models are shown here to either be incompatible with a consistent space-time structure, or to have physical singularities. The latter happens in spite of a non-zero scale factor in the isotropic background dynamics. A new effective Friedmann equation shows that a bounce is obtained at sub-Planckian densities, preceded by a physical singularity at infinite scale factor that resembles a time-reversed big rip. The entire phase is accompanied by rapid changes of the Hubble radius. In addition, a new version of perturbative inhomogeneity in loop quantum cosmology is introduced that maintains a consistent space-time structure and has a non-singular background dynamics.


[108] 2507.08234

Maneuver Detection via a Confidence Dominance Maneuver Indicator

Accurate and efficient maneuver detection is critical for ensuring the safety and predictability of spacecraft trajectories. This paper presents a novel maneuver detection approach based on comparing the confidence levels associated with the orbital state estimation and the observation likelihood. First, a confidence-dominance maneuver indicator (CDMI) is proposed by setting a confidence level for the state estimation and computing the maximum likelihood of the observation and its confidence level. The CDMI then flag a maneuver when the observation's confidence level exceeds that of the state estimation, indicating that the observation is unlikely under the no-maneuver hypothesis while maintaining consistency with the prior state estimation confidence. To efficiently compute the maximum likelihood of the observation and obtain the CDMI, a recursive polynomial optimization method is developed, taking advantage of convex optimization and polynomial approximation. In addition, an integrated CDMI approach is developed to eliminate the need to manually select the state confidence level. The integrated CDMI approach maintains high detection accuracy while simultaneously providing an indication of maneuver likelihood, thereby enhancing robustness and practical applicability. The performance of the proposed CDMI-based maneuver detection approaches is evaluated against an optimal control distance metric and two mixture-based approaches. The simulation results demonstrate that the proposed integrated CDMI approach can achieve up to 99.33\% detection accuracy, at least 10% higher than the competing methods, while substantially reducing computational costs.


[109] 2507.08318

Improving gravitational wave search sensitivity with TIER: Trigger Inference using Extended strain Representation

We introduce a machine learning (ML) framework called $\texttt{TIER}$ for improving the sensitivity of gravitational wave search pipelines. Typically, search pipelines only use a small region of strain data in the vicinity of a candidate signal to construct the detection statistic. However, extended strain data ($\sim 10$ s) in the candidate's vicinity can also carry valuable complementary information. We show that this information can be efficiently captured by ML classifier models trained on sparse summary representation/features of the extended data. Our framework is easy to train and can be used with already existing candidates from any search pipeline, and without requiring expensive injection campaigns. Furthermore, the output of our model can be easily integrated into the detection statistic of a search pipeline. Using $\texttt{TIER}$ on triggers from the $\texttt{IAS-HM}$ pipeline, we find up to $\sim 20\%$ improvement in sensitive volume time in LIGO-Virgo-Kagra O3 data, with improvements concentrated in regions of high masses and unequal mass ratios. Applying our framework increases the significance of several near-threshold gravitational-wave candidates, especially in the pair-instability mass gap and intermediate-mass black hole (IMBH) ranges.


[110] 2507.08345

Defining eccentricity for spin-precessing binaries

Standardizing the definition of eccentricity is necessary for unambiguous inference of the orbital eccentricity of compact binaries from gravitational wave observations. In previous works, we proposed a definition of eccentricity for systems without spin-precession that relies solely on the gravitational waveform, is applicable to any waveform model, and has the correct Newtonian limit. In this work, we extend this definition to spin-precessing systems. This simple yet effective extension relies on first transforming the waveform from the inertial frame to the coprecessing frame, and then adopting an amplitude and a phase with reduced spin-induced effects. Our method includes a robust procedure for filtering out spin-induced modulations, which become non-negligible in the small eccentricity and large spin-precession regime. Finally, we apply our method to a set of Numerical Relativity and Effective One Body waveforms to showcase its robustness for generic eccentric spin-precessing binaries. We make our method public via Python implementation in \texttt{gw\_eccentricity}.


[111] 2507.08431

Null geodesics, causal structure, and matter accretion in Lorentzian-Euclidean black holes

Recently, we introduced the Lorentzian-Euclidean black hole, a static and spherically symmetric solution of vacuum Einstein equations that exhibits a change in metric signature across the event horizon. In this framework, the analysis of radial trajectories of freely falling bodies proves that the central singularity can be avoided via a mechanism we interpret as atemporality, which is responsible for the shift of the time variable from real to imaginary values. In this paper, we further explore this model by first examining the behavior of null geodesics. Our investigation requires a set of signature-adaptive coordinate changes that generalize the local Lorentz transformations underlying General Relativity. We find that photon orbits, like their massive counterparts, cannot traverse the event horizon, thereby strengthening the previous result on the impossibility to reach the $r=0$ singularity. Additionally, we discuss the causal structure of the spacetime, provide the corresponding Penrose diagram, and analyze the process of matter accretion in the outer region of the black hole.


[112] 2507.08737

Quantum production of gravitational waves after inflation

A variety of mechanisms in the early Universe lead to the generation of gravitational waves (GWs). We introduce here a novel source of GWs generated by vacuum fluctuations after inflation. Being gravitons minimally coupled particles, their quantum creation takes place during inflation, but is absent in an unperturbed Universe during the radiation-dominated epoch, since they behave as conformally invariant particles. However, the presence of inhomogeneities breaks the conformal flatness of the metric, allowing scalar metric perturbations to induce the quantum production of gravitons. We compute the resulting GW spectrum from this mechanism for different models of the primordial scalar power spectrum. We find that this GW signal peaks around the GHz frequency range, distinguishing it from other astrophysical and cosmological backgrounds and underscoring the need for detectors sensitive to these high frequencies.


[113] 2507.08739

Thermal Gravitons from Warm Inflation

In warm inflation (WI), the persistent thermal bath that is sustained by dissipative interactions with the inflaton field produces a stochastic background of gravitational waves (GWs). In this paper we study the production and evolution of these GWs. Specifically, we investigate the emission of thermal gravitons (gravitons emitted by a thermal bath) from particle scattering in the bath and the evolution of the corresponding GWs. We find that the bulk of thermal graviton production in WI occurs during the transition to radiation domination after inflation. Further, the energy density of thermal gravitons is enhanced by roughly one to two orders of magnitude compared to that in a radiation-dominated scenario with the same reheating temperature. We also calculate the spectrum of the resulting stochastic GW background and find that it has a distinctive shape, consisting of a peak at high frequencies ~100 GHz and an almost flat spectrum extending to low frequencies. The peak arises from emission of sub-horizon modes that follow the temperature of the bath. The flat part of the spectrum corresponds to the modes that exit the horizon during WI and re-enter during radiation domination. We show that the detection prospects for the high-frequency peak of the GW spectrum, while improved slightly compared to the radiation-dominated case, still remain challenging. The thermal spectrum's low-frequency plateau is typically subdominant to the amplitude of the standard vacuum tensor modes from inflation, although WI models can exist where the thermal graviton plateau surpasses the vacuum contribution without exceeding current observational limits on the tensor-to-scalar ratio. Furthermore, we calculate the thermal graviton contribution from WI to dark radiation and show that WI models are generally expected to satisfy current observational bounds, including those from the cosmic microwave background.


[114] 2507.08789

Potential science with GW250114 -- the loudest binary black hole merger detected to date

On January 14, 2025 the LIGO interferometers detected a gravitational wave from the merger of two black holes, GW250114. Using publicly available information, we estimate that the signal-to-noise ratio (SNR) of GW250114 was $\sim 80$. This would make it three to four times louder than any other gravitational wave detected to date. GW250114 therefore offers a unique opportunity to make precise measurements of its source parameters and to test general relativity. In anticipation of its public data release, we analyze a set of simulated signals that have parameters similar to what we estimate for GW250114 and explore what new insights may be gained from this significant event.


[115] 2312.06434

Measurement of the Weyl potential evolution from the first three years of Dark Energy Survey data

The Weyl potential, which is the sum of the spatial and temporal distortions of the Universe's geometry, provides a direct way of testing the theory of gravity and the validity of the $\Lambda$CDM (Lambda Cold Dark Matter) model. Here we present measurement of the Weyl potential at four redshifts bins using data from the first three years of observations of the Dark Energy Survey. We find that the measured Weyl potential is 2$\sigma$, respectively 2.8$\sigma$, below the $\Lambda$CDM predictions in the two lowest redshift bins. We show that these low values of the Weyl potential are at the origin of the tension between Cosmic Microwave Background measurements and weak lensing measurements, regarding the parameter $\sigma_8$ which quantifies the clustering of matter. Interestingly, we find that the tension remains if no information from the Cosmic Microwave Background is used. Dark Energy Survey data on their own prefer a high value of the primordial fluctuations, together with a slow evolution of the Weyl potential. An important feature of our method is that the measurements of the Weyl potential are model-independent and can therefore be confronted with any theory of gravity, allowing efficient tests of models beyond General Relativity.


[116] 2401.06558

Heat loss and internal dynamics of Venus from lithosphere strength

The absence of plate tectonics and the young surface age (0.3-1 billion years) of Venus have led to diverse geodynamic models for Venus. The energetics of the Venusian interior drive these models; however, the lack of direct constraints on surface heat flow hampers their quantitative assessment. Here we present a global heat flow model for Venus, as well as estimates of the total heat loss, obtained from an inversion of geophysical data, including lithospheric effective elastic thickness, crustal thickness, and radioactive heat production. Heat flow on Venus is lower and less geographically structured than on Earth, with an average of 31 mW m$^{-2}$, but with highs associated to rifts systems reaching values typical of active terrestrial areas. The obtained total heat loss is 11-17 TW, similar to estimates of the total radioactive heat production. Therefore, at present, Venus proportionally dissipates much less heat than Earth. Furthermore, the calculated crustal temperatures imply that crustal melting or eclogitization are not dominant processes in the Venusian crust.


[117] 2409.17860

Gamma-ray burst pulse structures and emission mechanisms

The prompt phase X- and $\gamma$-ray light curves of gamma-ray bursts (GRBs) exhibit erratic and complex behaviour, often with multiple pulses. The temporal shape of individual pulses is often modelled as 'fast rise exponential decay' (FRED). Here, we introduce a novel fitting function to quantify pulse asymmetry. We conduct a light curve and a time-resolved spectral analysis on 61 pulses from 22 GRBs detected by the Fermi Gamma-ray Burst Monitor. Contrary to previous claims, we find that only $\sim 50\%$ of pulse lightcurves in our sample show a FRED shape, while about $25\%$ have a symmetric lightcurve, and the other $25\%$ have a mixed shape. Furthermore, our analysis reveals a clear trend: in multi-pulse bursts, the initial pulse tends to exhibit the most symmetric light curve, while subsequent pulses become increasingly asymmetric, adopting a more FRED-like shape. Additionally, we correlate the temporal and spectral shapes of the pulses. By fitting the spectra with the classical "Band" function, we find a moderate positive Spearman correlation index of 0.23 between pulse asymmetry and the low-energy spectral index $\alpha_{\max}$ (the maximum value across all time bins covering an individual pulse). Thus, during GRB light curves, the pulses tend to get more asymmetric and spectrally softer with time. We interpret this as a transition in the dominant emission mechanism from photospheric (symmetric-like and hard) to non-thermal emission above the photosphere and show that this interpretation aligns with a GRB jet Lorentz factor of the order of a few 10s in many cases.


[118] 2409.18185

Long-lived neutron-star remnants from asymmetric binary neutron star mergers: element formation, kilonova signals and gravitational waves

We present 3D general-relativistic neutrino-radiation hydrodynamics simulations of two asymmetric binary neutron star mergers producing long-lived neutron stars remnants and spanning a fraction of their cooling time scale. The mergers are characterized by significant tidal disruption with neutron rich material forming a massive disc around the remnant. The latter develops one-armed dynamics that is imprinted in the emitted kilo-Hertz gravitational waves. Angular momentum transport to the disc is initially driven by spiral-density waves and enhanced by turbulent viscosity and neutrino heating on longer timescales. The mass outflows are composed by neutron-rich dynamical ejecta of mass ${\sim}10^{-3}-10^{-2}M_\odot$ followed by a persistent spiral-wave/neutrino-driven wind of ${\gtrsim}10^{-2}M_\odot$ with material spanning a wide range of electron fractions, ${\sim}0.1-0.55$. Dynamical ejecta (winds) have fast velocity tails up to ${\sim}0.8$ (${\sim}0.4$) c. The outflows are further evolved to days timescale using 2D ray-by-ray radiation-hydrodynamics simulations that include an online nuclear network. We find complete $r$-process yields and identify the production of $^{56}$Ni and the subsequent decay chain to $^{56}$Co and $^{56}$Fe. Synthetic kilonova light curves predict an extended (near-) infrared peak a few days postmerger originating from $r$-process in the neutron-rich/high-opacity ejecta and UV/optical peaks at a few hours (ten minutes) postmerger originating from weak $r$-process (free-neutron decay) in the faster ejecta components. Additionally, the fast tail of tidal origin generates kilonova afterglows potentially detectable in radio and X band on a few to ten years time scale. Quantitative effects originating from the tidal disruption merger dynamics are reflected in the multimessenger emissions.


[119] 2410.02857

Reconstructing Galaxy Cluster Mass Maps using Score-based Generative Modeling

We present a novel approach to reconstruct gas and dark matter projected density maps of galaxy clusters using score-based generative modeling. Our diffusion model takes in mock SZ and X-ray images as conditional inputs, and generates realizations of corresponding gas and dark matter maps by sampling from a learned data posterior. We train and validate the performance of our model by using mock data from a cosmological simulation. The model accurately reconstructs both the mean and spread of the radial density profiles in the spatial domain, indicating that the model is able to distinguish between clusters of different mass sizes. In the spectral domain, the model achieves close-to-unity values for the bias and cross-correlation coefficients, indicating that the model can accurately probe cluster structures on both large and small scales. Our experiments demonstrate the ability of score models to learn a strong, nonlinear, and unbiased mapping between input observables and fundamental density distributions of galaxy clusters. These diffusion models can be further fine-tuned and generalized to not only take in additional observables as inputs, but also real observations and predict unknown density distributions of galaxy clusters.


[120] 2412.14826

J-PLUS: Tomographic analysis of galaxy angular density and redshift fluctuations in Data Release 3. Constraints on photo-$z$ errors, linear bias, and peculiar velocities

The {\it Javalambre Photometric Local Universe Survey} (J-PLUS) is a {\it spectro-photometric} survey covering about 3,000~deg$^2$ in its third data release (DR3), and containing about 300,000 galaxies with high quality ({\it odds}$>0.8$) photometric redshifts (hereafter photo-$z$s). We use this galaxy sample to conduct a tomographic study of the counts and redshift angular fluctuations under Gaussian shells sampling the redshift range $z\in[0.05,0.25]$. We confront the angular power spectra of these observables measured under shells centered on 11 different redshifts with theoretical expectations derived from a linear Boltzmann code ({\tt ARFCAMB}). Overall we find that J-PLUS DR3 data are well reproduced by our linear, simplistic model. We obtain that counts (or density) angular fluctuations (hereafter ADF) are very sensitive to the linear galaxy bias $b_g(z)$, although weakly sensitive to radial peculiar velocities of the galaxy field, while suffering from systematics residuals for $z>0.15$. Angular redshift fluctuations (ARF), instead, show higher sensitivity to radial peculiar velocities and also higher sensitivity to the average uncertainty in photo-$z$s ($\sigma_{\rm Err}$), with no obvious impact from systematics. For $z<0.15$ both ADF and ARF agree on measuring a monotonically increasing linear bias varying from $b_g(z=0.05)\simeq 0.9\pm 0.06$ up to $b_g(z=0.15)\simeq 1.5\pm 0.05$, while, by first time, providing consistent measurements of $\sigma_{\rm Err}(z)\sim 0.014$ that are $\sim 40~\%$ higher than estimates from the photo-$z$ code {\tt LePhare}, ($\sigma_{\rm Err}^{\rm LePhare}=0.010$). As expected, this photo-$z$ uncertainty level prevents the detection of radial peculiar velocities in the modest volume sampled by J-PLUS DR3, although prospects for larger galaxy surveys of similar (and higher) photo-$z$ precision are promising.


[121] 2412.14827

The J-PLUS collaboration. Additive versus multiplicative systematics in surveys of the large scale structure of the Universe

Observational and/or astrophysical systematics modulating the observed number of luminous tracers can constitute a major limitation in the cosmological exploitation of surveys of the large scale structure of the universe. Part of this limitation arises on top of our ignorance on how such systematics actually impact the observed galaxy/quasar fields. In this work we develop a generic, hybrid model for an arbitrary number of systematics that may modulate observations in both an additive and a multiplicative way, after applying a nonlinear power law transformation. This model allows us devising a novel algorithm that addresses the identification and correction for either additive and/or multiplicative contaminants. We test this model on galaxy mocks and systematics templates inspired from data of the third data release of the {\it Javalambre Photometric Local Universe Survey} (J-PLUS). We find that our method clearly outperforms standard methods that assume either an additive or multiplicative character for all contaminants in scenarios where both characters are actually acting on the observed data. In simpler scenarios where only an additive or multiplicative imprint on observations is considered, our hybrid method does not lie far behind the corresponding simplified, additive/multiplicative methods. Nonetheless, in scenarios of mild/low impact of systematics, we find that our hybrid approach converges towards the standard method that assumes additive contamination, as predicted by our model describing systematics. Our methodology also allows for the estimation of biases induced by systematics residuals on different angular scales and under different observational configurations, although these predictions necessarily restrict to the subset of {\em known/identified} potential systematics, and say nothing about ``unknown unknowns" possibly impacting the data.


[122] 2412.17912

Gravitational Waves from Particles Produced from Bubble Collisions in First-Order Phase Transitions

We discuss a new source of gravitational waves (GWs) from first-order phase transitions. The collisions of bubbles of the new phase can efficiently produce particles that couple to the background field undergoing the transition, thereby transferring a significant fraction of the released vacuum energy into a distribution of inhomogeneous and dynamic particle populations that persist long after the bubbles have disappeared. We study the GWs produced by such particle distributions, showing that GWs arise from the quadrupolar anisotropy in the radiation emitted from the bubble collisions, and present a semi-analytical calculation of the two-point correlation function for the associated energy distributions. We find that this new contribution can qualitatively modify the overall GW signal from such phase transitions, creating a distinct shift in the spectral slope at low frequencies that could be observed by future GW experiments. It is therefore important to take this new contribution into account for any transition where the background field has significant self-coupling or couplings to other fields that could lead to efficient particle production at bubble collision.


[123] 2501.12138

Globular cluster formation from inertial inflows: accreting extremely massive stars as the origin of abundance anomalies

We use the inertial-inflow model of massive star formation to describe the formation of globular clusters (GCs) in turbulent molecular clouds. A key aspect of this model is that the maximum stellar mass scales linearly with cloud mass, such that extremely massive stars (EMSs, $10^{3-4}\,\msun$) form in massive GCs ($\gtrsim10^5\,\msun$). The total wind mass loss is dominated by accreting EMSs (aEMSs), whose wind mass-loss rates have become comparable to their accretion rates ($\gtrsim10^{-2}\,\msun\,\yr^{-1}$). These winds pollute the intra-cluster medium with hot-hydrogen burning yields during GC formation. We propose a parameterised model for the evolution of the stellar mass function during GC formation ($\sim 1-2\,\myr$), accounting for gas inflow, wind mass loss and mixing of aEMS yields with pristine gas that has initial proto-GC abundances. Low-mass stars ($\lesssim1\,\msun$) form continuously from this mixed gas and their abundances resemble observed abundance trends with GC mass and metallicity, specifically: (i) the helium spread in a typical GC is small ($\Delta Y \simeq 0.01$) and increases with GC mass; (ii) the fraction of polluted stars increases with GC mass and metallicity; (iii) the extent of the Mg-Al anticorrelations is more pronounced in metal-poor and massive GCs. We conclude that GCs formed with a population of EMSs from gas with surface densities $\gtrsim10^3\,\msun\,\pc^{-2}$ and that nitrogen-rich galaxies discovered by the James Webb Space Telescope ({\it JWST}) are dominated by EMS-rich GCs that formed in the earliest phases of galaxy formation. These EMSs may have left behind intermediate-mass black holes with masses above the pair-instability gap ($\gtrsim120\,\msun$) that could be found with ongoing gravitational wave experiments.


[124] 2501.14025

Cataclysmic Variables in Triples: Formation Models and New Discoveries

The formation of cataclysmic variables (CVs) has long been modeled as a product of common envelope evolution (CEE) in isolated binaries. However, a significant fraction of intermediate-mass stars -- the progenitors of the white dwarfs (WDs) in CVs -- are in triples. We therefore investigate the importance of triple star dynamics in CV formation. Using Gaia astrometry and existing CV catalogs, we construct a sample of $\sim50$ CVs in hierarchical triples within 1 kpc of the Sun, containing main-sequence (MS) and WD tertiaries at separations of 100 - 30,000 au. We infer that at least 10% of CVs host wide tertiaries. To interpret this discovery, we evolve a population of 2000 triples using detailed three-body simulations, 47 of which become CVs. We predict that 20% of CVs in triples form without ever experiencing CEE, where the WD and donor are brought together by the eccentric Kozai-Lidov (EKL) mechanism after the formation of the WD. These systems favor larger donor stars and longer birth orbital periods (8-20 hrs) than typical CVs. Among systems that do undergo CEE, about half would not have interacted without the presence of the tertiary. Triple formation channels both with and without CEE require initially wide inner orbits ($\gtrsim 1$ au), which in turn require larger tertiary separations to be stable. Consistent with this prediction, we find that the observed Gaia CV triples have wider separations on average than normal wide binaries selected in the same way. Our work underscores the importance of triples in shaping interacting binary populations including CVs, ultracompact binaries, and low-mass X-ray binaries.


[125] 2502.10758

Gravitational wave lensing: probing Fuzzy Dark Matter with LISA

Gravitational lensing is a universal phenomenon: it affects both gravitational waves (GWs) and electromagnetic signals travelling through the gravitational field of a massive object. In this work, we explore the prospects of observing lensed GW signals from the mergers of massive black holes, lensed by dark matter halos composed of Fuzzy Dark Matter (FDM), which form dense cores known as solitons. We focus on wave optics phenomena, where frequency-dependent signatures can be observed in the weak lensing regime (i.e. single-image). Our results show that lensing diffraction signatures differ for low-mass halos in FDM, and can reveal the presence of a solitonic core. Furthermore, we demonstrate that FDM and cold dark matter profiles can be distinguished in GW signals from binary massive black hole mergers, which will be observed by the Laser Interferometer Space Antenna (LISA) mission. However, the dense solitonic core does not substantially enhance the detectability of FDM halos at large source-lens offsets, relative to standard cold dark matter. Our analysis confirms FDM halos as a promising signature of dark matter on GW observations.


[126] 2502.12667

The Preference for Evolving Dark Energy from Cosmological Distance Measurements and Possible Signatures in the Growth Rate of Perturbations

In this study, we use a flexible parametrization of the equation of state of dark energy to explore its possible evolution with datasets from the Dark Energy Spectroscopic Instrument (DESI), Planck cosmic microwave background, and either the 5-year Dark Energy Survey (DES) or the Pantheon+ (PP) supernova (SN) compilation. This parametrization, called transitional dark energy, allows for rapid changes in the equation of state but also changes like that in the Chevallier-Polarski-Linder parametrization. We find a 3.8{\sigma} preference for evolving dark energy over {\Lambda}CDM with the DES dataset and a weaker 2.4{\sigma} preference when using the PP dataset. This corroborates the finding of the DESI Collaboration, who found that their baryon acoustic oscillation data preferred evolving dark energy when fit with the CPL parametrization of the equation of state. Our analysis reveals no significant outliers in the DESI data around the TDE best-fit, while the data is asymmetrically distributed around the {\Lambda}CDM best-fit model such that the measured distances are on average smaller. The DESI and SN data both prefer an expansion history that implies a higher dark energy density around z=0.5 than in the Planck-{\Lambda}CDM model, with the inferred equation of state being greater than -1 around z=0 and close to or below -1 at z>0.5. We show that when the expansion rate is greater than that in the Planck-{\Lambda}CDM model (around z=0.5), the growth rate calculated assuming General Relativity is suppressed relative to the Planck-{\Lambda}CDM model, and it rebounds as the expansion rate differences between the models become smaller closer to the present time. The resulting flattening of the $f\sigma_8(z)$ curve compared to the {\Lambda}CDM model could be an independent signature of the temporal evolution of dark energy.


[127] 2503.12372

How stable are the impact factors of major astrophysics journals

We examine the long-term trends in impact factor over the past decade (2015-2014) for seven flagship journals widely used in Astrophysics: ApJ, AJ, ApJL, ApJS,A &A, JCAP, and MNRAS. We check the variation of both the traditional impact factor as well as the median-based impact factor, which we had studied in a previous work. We find that ApJS exhibits the largest drop in impact factor from 14.9 to 9.0, while ApJ and MNRAS show a steady decline from 6.7 to 5.2 and 5.6 to 5.0, respectively. However, the impact factors for AJ and ApJL have shown a steady increase from 4.55 to 5.35 and 6.4 to 9.7, respectively. The median-based impact factors are much more stable showing variations of at most $\pm 1$ within the last decade.


[128] 2503.23856

Explodability criteria for the neutrino-driven supernova mechanism

Massive stars undergoing iron core-collapse at the end of their evolution terminate their lives either in successful or failed supernovae (SNe). The physics of core-collapse supernovae (CCSNe) is complex, and their understanding requires computationally expensive simulations. Using these to predict CCSN outcomes over large, densely sampled parameter spaces of SN progenitors, as is needed e.g. for population synthesis studies, is thus not feasible. To remedy this situation, we present explodability criteria that allow us to predict the final fates of stars by evaluating stellar structure variables at the onset of core-collapse. The criteria are calibrated to predictions of a semi-analytical SN model, evaluated over a set of $\sim$~3,900 heterogeneous stellar progenitors (single, binary-stripped and accretor stars). Over these, the criteria achieve an accuracy of >99\% agreement with the semi-analytical model. The criteria are tested on 29 state-of-the-art 3D CCSN simulation outcomes from two different groups. Furthermore, we find that all explodability proxies needed for our pre-SN structure-based criteria have two distinct peaks and intervening valleys as a function of the carbon-oxygen (CO) core mass $M_\mathrm{CO}$, which coincide with failed and successful SNe, respectively. The CO core masses of explodability peaks shift systematically with metallicity, $Z$, and with timing of hydrogen-rich envelope removal by binary mass transfer. With these, we identify critical values in $M_\mathrm{CO}$ that define windows over which black holes form by direct collapse and formulate a CCSN recipe based on $M_\mathrm{CO}$ and $Z$, applicable for rapid binary population synthesis and other studies. Our explodability formalism is consistent with observations of Type~IIP, IIb/Ib and Ic supernova progenitors and partially addresses the missing Red Supergiant Problem by direct black hole formation.


[129] 2504.06247

The Resolved Structure of a Low Metallicity Photodissociation Region

Photodissociation Regions (PDRs) are key to understanding the feedback processes that shape interstellar matter in galaxies. One important type of PDR is the interface between HII regions and molecular clouds, where far-ultraviolet (FUV) radiation from massive stars heats gas and dissociates molecules. Photochemical models predict that the C/CO transition occurs deeper in the PDR compared to the H/H$_2$ transition in low-metallicity environments, increasing the extent of CO-dark H$_2$ gas. This prediction has been difficult to test outside the Milky Way due to the lack of high spatial resolution observations tracing H$_2$ and CO. This study examines a low-metallicity PDR in the N13 region of the Small Magellanic Cloud (SMC) where we spatially resolve the ionization front, the H$_2$ dissociation front, and the C/CO transition using 12CO J=2-1, 3-2 and [CI] (1-0) observations from the Atacama Large Millimeter/sub-mm Array (ALMA) and near-infrared spectroscopy of the H$_2$ 2.12 1-0S(1) vibrational line, and H recombination lines from the James Webb Space Telescope (JWST). Our analysis shows that the separation between the H/H$_2$ and C/CO boundaries is approximately 0.043 $\pm$ 0.013(stat.) $\pm$ 0.0036(syst.) pc (equivalent to 0".146 $\pm$ 0".042(stat.) $\pm$ 0".012(syst.) at the SMC's distance of 62 kpc), defining the spatial extent of the CO-dark H$_2$ region. Compared to our plane-parallel PDR models, we find that a constant pressure model matches the observed structure better than a constant density one. Overall, we find that the PDR model does well at predicting the extent of the CO-dark H$_2$ layer in N13. This study represents the first resolved benchmark for low metallicity PDRs.


[130] 2504.09287

Looking into the Jet Cone of the Neutrino-Associated Very High Energy Blazar PKS 1424+240

The acceleration process of massive particles as well as the production of very high energy (VHE) photons and neutrinos remains a fundamental challenge in astrophysics. We investigate the parsec-scale jet structure and magnetic field of the blazar PKS 1424+240, that has been selected on the basis of strong VHE gamma-ray emission and identified with one of the highest peaks in the IceCube 9-year neutrino sky. We analyze 15 GHz VLBA observations of this BL Lac object, stacking 42 polarization-sensitive images collected in 2009-2025 to enhance the signal and reveal persistent parsec-scale structure. Our observations uncover a rare scenario. The object is viewed inside the jet cone, very close to the axis of its relativistic jet, with a viewing angle of <0.6 deg. This effectively maximizes Doppler boosting to values ~30 and enhances both electromagnetic and neutrino emission in the observer's direction. Based on polarimetric observations, we unambiguously detect a net toroidal component of the jet's magnetic field, indicating a current carrying jet flowing almost directly towards our line of sight. Blazars with very small jet viewing angles offer a solution to the longstanding mismatch between Doppler factors inferred from low VLBI apparent jet speed and those derived from VHE observations -- the so-called "Doppler factor crisis". We show that relativistic beaming plays the critical role in the gamma-ray and neutrino emission of blazars, with direct implications for models of their multi-messenger emission.


[131] 2504.12897

OntoPortal-Astro, a Semantic Artefact Catalogue for Astronomy

The astronomy communities are widely recognised as mature communities for their open science practices. However, while their data ecosystems are rather advanced and permit efficient data interoperability, there are still gaps between these ecosystems. Semantic artefacts (SAs) -- e.g., ontologies, thesauri, vocabularies or metadata schemas -- are a means to bridge that gap as they allow to semantically described the data and map the underlying concepts. The increasing use of SAs in astronomy presents challenges in description, selection, evaluation, trust, and mappings. The landscape remains fragmented, with SAs scattered across various registries in diverse formats and structures -- not yet fully developed or encoded with rich semantic web standards like OWL or SKOS -- and often with overlapping scopes. Enhancing data semantic interoperability requires common platforms to catalog, align, and facilitate the sharing of FAIR (Findable, Accessible, Interoperable and Reusable) SAs. In the frame of the FAIR-IMPACT project, we prototyped a SA catalogue for astronomy, heliophysics and planetary sciences. This exercise resulted in improved vocabulary and ontology management in the communities, and is now paving the way for better interdisciplinary data discovery and reuse. This article presents current practices in our discipline, reviews candidate SAs for such a catalogue, presents driving use cases and the perspective of a real production service for the astronomy community based on the OntoPortal technology, that will be called OntoPortal-Astro.


[132] 2504.14255

Inferring Dense Confined Circumstellar Medium around Supernova Progenitors via Long-term Hydrodynamical Evolution

Circumstellar interaction of supernova (SN) ejecta is an essential process in its evolution and observations of SNe have found the signature of circumstellar interaction both in the early and late evolutionary phase of SNe. In this Letter, we show that if the SN forward shock plunges into tenuous stellar wind from dense circumstellar medium (CSM) in the vicinity of the progenitor (i.e., confined CSM), the subsequent time evolutions of the SN-CSM interaction system deviates from the prediction of self-similar solution. In this case, after all of the confined CSM is swept up by the SN forward shock (roughly $10$ days after the explosion), the propagation of the shocked shell will be driven by the freely expanding ram pressure of the confined CSM component, instead of the SN ejecta. Meanwhile, the forward shock decelerates faster than the prediction of thin-shell approximation once the confined CSM component reaches homologous expansion. This lasts until the reverse shock in the confined CSM component reaches the head of the SN ejecta, leading to the restoration of the system into the evolutionary model without confined CSM, where the SN ejecta drives the expansion of the system. We also show that this peculiar evolution will be reflected in observational signatures originating from SN-CSM interaction, taking rapid decline and rebrightening of radio emission as examples. Our results shed light on the importance of taking into account the effect of initial SN-CSM interaction even when we focus on observational properties of SNe a few years after the explosion.


[133] 2505.15455

Magnetohydrodynamic turbulence and the associated spatial diffusion tensor of cosmic rays in dynamical galactic halos

A detailed understanding of cosmic-ray transport in galactic halos is essential for explaining various observations, such as the radio continuum measurements of synchrotron radiation from energetic electrons. Of central importance is the spatial diffusion tensor of cosmic rays, which can be computed in an ab~initio manner if the turbulence in the background medium is known. The study aimed to establish a suitable framework to compute the evolution of magnetohydrodynamic (MHD) turbulence and, hence, the diffusion tensor in the dynamical halos of galaxies. The Reynolds-averaged single-fluid MHD equations were solved numerically on a cylindrical grid, assuming axial symmetry and fixed boundary conditions on a central ellipsoid representing the galaxy. The physical properties of both large-scale MHD and small-scale turbulent quantities, including the coefficients of parallel and perpendicular diffusion, were extracted from the resulting steady state. Hydrodynamic validation revealed a persistent instability of the near-axis flow, which could be traced to the galaxy's gravitating mass. The results for the evolution of MHD turbulence in a galactic halo -- using parameters approximating those of NGC4631 -- enabled an ab initio computation of the spatial diffusion tensor of cosmic rays through the application of a state-of-the-art nonlinear theory. The simulation results reveal variation in turbulence quantities, (i.e., fluctuation energy, correlation scale, and cross helicity) in a dynamical halo. Furthermore, the corresponding diffusion tensor of cosmic rays exhibits significant variation throughout such a halo.


[134] 2505.19784

Do multifrequency polarimetric observations of BL Lac rule out a hadronic origin for its X-ray emission?

Recent multifrequency polarimetric observations of the eponymous blazar BL Lac reveal an extremely large degree of polarization in the optical band (average of $25\%$, reaching $45\%$), together with a small ($\lesssim 7\%$) degree of polarization in the X-ray band. This has been interpreted as evidence that the X-rays are produced through inverse Compton emission by relativistic electrons, thus ruling out alternative models based on hadronic processes. Here we revisit the observational evidence, interpreting it in a framework where the observed radiation is entirely produced through synchrotron emission. Electrons produce the radio-to-optical component and protons produce the X-rays and the gamma-rays. We determine the jet magnetic fields from an MHD model of magnetically dominated stationary axisymmetric outflows, and show that the X-ray emission from the protons is naturally less polarized than the optical emission from the electrons. The model parameters required to reproduce the multifrequency polarimetric observations are fully compatible with blazar jets.


[135] 2506.01733

Observation of the Crab Nebula with the Single-Mirror Small-Size Telescope stereoscopic system at low altitude

The Single-Mirror Small-Size Telescope (SST-1M) stereoscopic system is composed of two Imaging Atmospheric Cherenkov Telescopes (IACTs) designed for optimal performance for gamma-ray astronomy in the multi-TeV energy range. It features a 4-meter-diameter tessellated mirror dish and an innovative SiPM-based camera. Its optical system features a 4-m diameter spherical mirror dish based on the Davies-Cotton design, maintaining a good image quality over a large FoV while minimizing optical aberrations. In 2022, two SST-1M telescopes were installed at the Ondřejov Observatory, Czech Republic, at an altitude of 510 meters above sea level, and have been collecting data for commissioning and astronomical observations since then. We present the first SST-1M observations of the Crab Nebula, conducted between September 2023 and March 2024 in both mono and stereoscopic modes. During this observation period, 46 hours for the SST-1M-1 and 52 hours for the SST-1M-2 were collected for which 33 hours are in stereoscopic mode. We use the Crab Nebula observation to validate the expected performance of the instrument, as evaluated by Monte Carlo simulations carefully tuned to account for instrumental and atmospheric effects. We determined that the energy threshold at the analysis level for the zenith angles below $30^\circ$ is 1 TeV for mono mode and 1.3 TeV for stereo mode. The energy and angular resolutions are approximately 20% and $0.18^\circ$ for mono mode and 10% and $0.10^\circ$ for stereo mode, respectively. We present the off-axis performance of the instrument and a detailed study of systematic uncertainties. The results of a full simulation of the telescope and its camera is compared to the data for the first time, allowing a deep understanding of the SST-1M array performance.


[136] 2506.22272

A pipeline to search for signatures of line-of-sight acceleration in gravitational wave signals produced by compact binary coalescences

Compact binary coalescences (CBCs), such as merging binary black holes (BBHs), binary neutron stars (BNSs), or neutron star black holes (NSBHs), hosted by dense stellar environments, could produce gravitational waves (GWs) that contain signatures of line-of-sight acceleration (LOSA) imparted by the environment's gravitational potential. We calculate the Post-Newtonian (PN) corrections to the $(2,\,2)$ mode GW phase due to a finite LOSA, starting from the leading order at -4 PN below the quadrupole order, up to 3.5 PN above the quadrupole order. We do so for binaries whose component spins are aligned with the orbital angular momentum, as well as for binaries with non-zero tidal deformation. We implement these corrections into the LIGO-Virgo-Kagra (LVK) collaboration's flagship parameter estimation (PE) software \textsc{Bilby\_tgr}. We study the systematics associated with recovering LOSAs. We find that, when the injection and recovery waveform models are identical, LOSAs are recovered as expected. We test the robustness of the pipeline against waveforms with strong higher-mode signatures or signatures of beyond-general-relativistic (beyond-GR) effects, to delimit the range of applicability of our GR-consistent quasi-circular LOSA-corrected waveforms.


[137] 2506.22547

MOSDEF-3D: Keck/OSIRIS Maps of the Ionized ISM in $z \sim 2$ Galaxies

We present spatially-resolved rest-frame optical emission line maps of four galaxies at $z \sim 2$ observed with Keck/OSIRIS to study the physical conditions of the ISM at Cosmic Noon. Our analysis of strong emission line ratios in these galaxies reveals an offset from the local star-forming locus on the BPT diagram, but agrees with other star-forming galaxies at similar redshifts. Despite the offset towards higher [O III]$\lambda5008$/H$\beta$ and [N II]$\lambda6585$/H$\alpha$, these strong-line ratios remain consistent with or below the maximum starburst threshold even in the inner $\sim 1$ kpc region of the galaxies, providing no compelling evidence for central AGN activity. The galaxies also exhibit flat radial gas-phase metallicity gradients, consistent with previous studies of $z \sim 2$ galaxies and suggesting efficient radial mixing possibly driven by strong outflows from intense star formation. Overall, our results reveal the highly star-forming nature of these galaxies, with the potential to launch outflows that flatten metallicity gradients through significant radial gas mixing. Future observations with JWST/NIRSpec are crucial to detect fainter emission lines at higher spatial resolution to further constrain the physical processes and ionization mechanisms that shape the ISM during Cosmic Noon.


[138] 2507.01671

Constraints on Earth's Core-Mantle boundary from nutation

Periodic variations in the Sun and Moon's gravitational pull cause small changes in Earth's rotational axis direction called nutation. Nutation components in the retrograde quasi-diurnal frequency band measured in the terrestrial reference frame are amplified by resonance with the Free Core Nutation (FCN), a rotational mode of Earth's fluid core. Dissipative processes at the core-mantle boundary (CMB) dampen this resonance, contributing to the observed phase lag between tidal forcing and Earth's rotational response. This phase lag is commonly attributed to electromagnetic (EM) coupling between the core and the electrically conducting lower mantle. However, estimates of mantle conductivity and radial magnetic field strength at the CMB suggest these effects are insufficient. We show that the missing dissipation arises naturally from the excitation of internal waves in the fluid core by topographic features at the CMB. Adapting a theoretical framework originally developed for tidal flow over oceanic topography, we compute the form drag and associated power flux induced by CMB topography. Our results are consistent with a CMB topography characterized by a root mean square amplitude of a few kilometers. The model favors weak stratification at the top of the core, though stronger stratification remains compatible with increased topographic amplitude. This mechanism provides independent constraints on CMB topography and stratification, complementing seismological and magnetic observations. Its generality offers a new framework for probing deep-interior dynamics across terrestrial planets.


[139] 2507.02777

Advanced techniques of searching for flares of ultra-high-energy photons from point sources

Astrophysical flares are one of the possible prominent source classes of ultra-high-energy (UHE, $E > 10^{17}$ eV) cosmic rays, which can be detected by recording clusters of extensive air showers in arrays of detectors. The search for sources of neutral particles offers distinct advantages over searching for sources of charged particles, as the former traverse cosmic distances undeflected by magnetic fields. While no cosmic-ray photons exceeding $10^{17}$ eV have been definitively detected, identifying the clustering of events in cosmic-ray data would provide compelling evidence for their existence. We compare two analysis methods for detecting direction-time clustering in UHE extensive air showers: an approach that examines multiplets, and the stacking method, which analyzes sets of doublets that are not necessarily consecutive, thus making it sensitive to multiple flares. Both techniques combine time-clustering algorithms with unbinned likelihood study. Background events (initiated by hadrons) can be more efficiently distinguished from photon-induced events (signals) by using a photon tag that employs probability distribution functions to classify each event as more likely to be initiated by either a photon or a hadron. We demonstrate that these methods can effectively distinguish between events initiated by photons and those initiated by hadrons (background), and can accurately reproduce both the number of photon events within flares and their duration. We calculate the discovery potentials, i.e., the number of events required to identify a photon flare. The methods discussed can be used to search for cosmic ray sources and/or improve limits on the fluxes of UHE photons.


[140] 2507.04066

Metallicities of M Dwarf Planet Host Stars from Kepler, K2, and TESS observed by APOGEE: Trends with Exoplanetary Radii and Orbital Periods

One important property in studying the exoplanet population is the host star metallicity ([M/H]). In this study, we derived stellar metallicities and oxygen abundances for 45 M dwarf stars (with 65 exoplanets) using the near-infrared high-resolution spectra from the SDSS APOGEE survey and synthetic spectra computed in LTE. We also investigated the exoplanetary radii distribution for a larger sample of 246 exoplanets orbiting 188 M dwarf stars. The [M/H] versus [O/M] distribution obtained indicates that our sample is composed mainly of thin disk stars, which follow the behavior of the low-alpha sequence in the Milky Way thin disk. Small planets with radii smaller than 3R$_{\oplus}$ were found around stars with a range of metallicities (-0.6$<$[M/H]$<$+0.3), while larger planets of the sample orbit only stars with [M/H]$>=0.0$. These results indicate that while small planets can form in different environments, larger planets preferentially form in metal-rich protoplanetary disks. Exoplanets with P$_{\rm orb}<$4.3 days orbit on average more metal-rich stars than planets with longer periods. This threshold is smaller than that found for FGK stars (8--10 days) and might be related to M dwarfs having a smaller dust sublimation radius. The distribution of exoplanets with R$_{\rm p}>$4R$_{\oplus}$ shows a concentration at orbital periods between 2 and 5 days, which may result from inward orbital migration. There is also a different behavior between single-detected exoplanets and planets from multiplanetary systems, with the latter being found on average around more metal-poor stars, and with planetary radii roughly up to 3 R$_{\oplus}$.


[141] 2507.07346

CLASS_SZ II: Notes and Examples of Fast and Accurate Calculations of Halo Model, Large Scale Structure and Cosmic Microwave Background Observables

These notes are very much work-in-progress and simply intended to showcase, in various degrees of details (and rigour), some of the cosmology calculations that class_sz can do. We describe the class_sz code in C, Python and Jax. Based on the Boltzmann code class, it can compute a wide range of observables relevant to current and forthcoming CMB and Large Scale Structure surveys. This includes galaxy shear and clustering, CMB lensing, thermal and kinetic Sunyaev and Zeldovich observables, Cosmic Infrared Background, cross-correlations and three-point statistics. Calculations can be done either within the halo model or the linear bias model. For standard $\Lambda$CDM cosmology and extensions, class_sz uses high-accuracy cosmopower emulators of the CMB and matter power spectrum to accelerate calculations. With this, along with efficient numerical integration routines, most class_sz output can be obtained in less than 500 ms (CMB $C_\ell$'s or matter $P(k)$ take $\mathcal{O}(1\mathrm{ms})$), allowing for fast or ultra-fast parameter inference analyses. Parts of the calculations are "jaxified", so the software can be integrated into differentiable pipelines.


[142] 2507.07616

A forecasting framework for galactic cosmic ray flux in space weather applications

The intensity and energy spectrum of galactic cosmic rays in the heliosphere are significantly influenced by the 11-year solar cycle, a phenomenon known as solar modulation. Understanding this effect and its underlying physical mechanisms is essential for assessing radiation exposure and associated risks during space missions. Starting from a previously developed effective predictive model of solar modulation, validated using cosmic ray flux measurements from space-based detectors such as PAMELA and AMS-02, we build a generalizable forecasting strategy for the long-term evolution of cosmic ray fluxes. This strategy is based on identifying delayed cross-correlation relationships between solar proxies and the model's parameters. It integrates recent findings on time lags between cosmic ray fluxes and solar activity, and incorporates advanced time-series signal processing techniques. The framework not only performs well in reproducing observed data, but also shows strong potential for applications in space radiation monitoring and forecasting. By efficiently capturing the long-term variability of galactic cosmic rays, our approach contributes valuable insights for evaluating radiation risks, ultimately supporting safer and more effective space exploration.


[143] 2411.10216

On the dark matter origin of an LDMX signal

Fixed target experiments where beam electrons are focused upon a thin target have shown great potential for probing new physics, including the sub-GeV dark matter (DM) paradigm. However, a signal in future experiments such as the light dark matter experiment (LDMX) would require an independent validation to assert its DM origin. To this end, we propose to combine LDMX and next generation DM direct detection (DD) data in a four-step analysis strategy, which we here illustrate with Monte Carlo simulations. In the first step, the hypothetical LDMX signal (i.e. an excess in the final state electron energy and transverse momentum distributions) is $\textit{recorded}$. In the second step, a DM DD experiment operates with increasing exposure to test the DM origin of the LDMX signal. Here, LDMX and DD data are simulated. In the third step, a posterior probability density function (pdf) for the DM model parameters is extracted from the DD data, and used to $\textit{predict}$ the electron recoil energy and transverse momentum distributions at LDMX. In the last step, $\textit{predicted}$ and $\textit{recorded}$ electron recoil energy and transverse momentum distributions are compared in a chi-square test. We present the results of this comparison in terms of a threshold exposure that a DD experiment has to operate with to assert whether $\textit{predicted}$ and $\textit{recorded}$ distributions $\textit{can}$ be statistically dependent. We find that this threshold exposure grows with the DM particle mass, $m_\chi$. It varies from 0.012 kg-year for a DM mass of $m_\chi=4$ MeV to 1 kg-year for $m_\chi=25$ MeV, which is or will soon be within reach.


[144] 2412.02651

Costs of Bayesian Parameter Estimation in Third-Generation Gravitational Wave Detectors: a Review of Acceleration Methods

Bayesian inference with stochastic sampling has been widely used to obtain the properties of gravitational wave (GW) sources. Although computationally intensive, its cost remains manageable for current second-generation GW detectors because of the relatively low event rate and signal-to-noise ratio (SNR). The third-generation (3G) GW detectors are expected to detect hundreds of thousands of compact binary coalescence (CBC) events every year with substantially higher SNR and longer signal duration, presenting significant computational challenges. In this study, we systematically evaluate the computational costs of CBC source parameter estimation (PE) in the 3G era by modeling the PE time cost as a function of SNR and signal duration. We examine the standard PE method alongside acceleration methods including relative binning, multibanding, and reduced order quadrature. We predict that PE for a one-month-observation catalog with 3G detectors could require at least billions of CPU core hours with the standard PE method, whereas acceleration techniques can reduce this demand to less than millions of core hours, which is as high as the cost of analyzing GW events in the past 10 years. These findings highlight the necessity for more efficient PE methods to enable cost-effective and environmentally sustainable data analysis for 3G detectors. In addition, we assess the accuracy of accelerated PE methods, emphasizing the need for careful treatment in high-SNR scenarios.


[145] 2412.14259

Black hole spectroscopy and nonlinear echoes in Einstein-Maxwell-scalar theory

In the context of Einstein-Maxwell-scalar theory with a nonminimal coupling between the electromagnetic and scalar field, we study linear (non)radial perturbations and nonlinear radial dynamics of spherically symmetric black holes. In a certain region of the parameter space, this theory admits hairy black holes with a stable photon sphere. This has a counterpart in the effective potential of linear perturbations, featuring multiple maxima and minima. The corresponding quasinormal mode spectrum contains long-lived modes trapped in the potential cavity and the time-domain linear response displays echoes, as previously observed for horizonless compact objects. Interestingly, the black-hole dynamics in this theory can be studied at the nonlinear level. By performing fully-fledged 1+1 simulations, we show that echoes are present even when the nonlinearities are significant. To our knowledge, this is the first example of echoes appearing in a consistent theory beyond a linearized analysis. In a follow-up work we will study whether this feature is also present in the post-merger signal from black hole collisions in this theory.


[146] 2502.14960

Lepton Flavor Asymmetries: from the early Universe to BBN

Large primordial lepton flavor asymmetries with almost vanishing total baryon-minus-lepton number can evade the usual BBN and CMB constraints if neutrino oscillations lead to perfect flavor equilibration. Solving the momentum averaged quantum kinetic equations (QKEs) describing neutrino oscillations and interactions, we perform the first systematic investigation of this scenario, uncovering a rich flavor structure in stark contradiction to the assumption of simple flavor equilibration. We find (i) a particular direction in flavor space, $\Delta n_e \simeq -2/3 \, (-1) \Delta n_\mu$ for normal (inverted) neutrino mass hierarchy, in which the flavor equilibration is efficient and primordial asymmetries are essentially unconstrained, (ii) a minimal washout factor, $\Delta n_e^2|_\mathrm{BBN} \leq 0.03 \, (0.016) \sum_\alpha \Delta n_\alpha^2|_{\mathrm{ini}}$ yielding a conservative estimate for the allowed primordial asymmetries in a generic flavor direction, and (iii) particularly strong or weak washout if one of the initial flavor asymmetries vanishes due to non-adiabatic muon- or electron-driven MSW transitions. These results open up the possibility of a first-order QCD phase transition facilitated by large lepton asymmetries as well as baryogenesis from large and compensated $\Delta n_e =- \Delta n_\mu$ asymmetries. Our first-principles approach of deriving momentum averaged QKEs includes collision terms beyond the damping approximation, energy transfer between the neutrino and electron-photon plasma, and provides a fast and reliable way to investigate the impact of primordial lepton asymmetries at the time of BBN. We publicly release the Mathematica code COFLASY-M on GitHub which solves the QKEs numerically.


[147] 2504.16156

Extreme mass ratio inspirals around topological stars

We study a point scalar charge in circular orbit around a topological star, a regular, horizonless soliton emerging from dimensional compactification of Einstein-Maxwell theory in five dimensions, which could describe qualitative properties of microstate geometries for astrophysical black holes. This is the first step towards studying extreme mass-ratio inspirals around these objects. We show that when the particle probes the spacetime close to the object, the scalar-wave flux deviates significantly from the corresponding black hole case. Furthermore, as the topological star approaches the black-hole limit, the inspiral can resonantly excite its long-lived modes, resulting in sharp features in the emitted flux. Although such resonances are too narrow to produce detectable dephasing, we estimate that a year-long inspiral down to the innermost stable circular orbit could accumulate a significant dephasing for most configurations relative to the black hole case. While a full parameter-estimation analysis is needed, the generically large deviations are likely to be within the sensitivity reach of future space-based gravitational-wave detectors.


[148] 2505.08543

Tangent point position and atmosphere composition in limb scanning instruments

In the forward model for limb-scanning instruments, ray tracing must be accounted for because variations in air refractivity cause the lines of sight to bend from straight paths into curves. The tangent point of a line of sight is defined as the minimum height point. The lines of sight can be adjusted by varying the nadir angles of the instrument, which must be calibrated to account for Earth's ellipticity and the actual atmospheric conditions sampled by these paths. In this study, we investigate the relationship between the nadir angles, atmospheric properties, and tangent point positions, using a configuration relevant to the CAIRT instrument. Atmospheric data from reanalysis databases are utilized, and the lines of sight are determined by numerically solving the Eikonal equation. Our findings are compared to those of a previous study conducted for the MIPAS instrument, highlighting key differences.


[149] 2507.01895

Freeze-In Dark Matter and Leptogenesis: a $ψ'$SM route

We investigate the possibility of \emph{freeze-in} dark matter production and baryogenesis via leptogenesis in a $\psi'$SM model, which is an $E_6$ extension of the Standard Model, featuring a residual $U(1)_{\psi'}$ gauge symmetry. This symmetry arises from a linear combination of $U(1)\chi$ and $U(1)_{\psi}$, both of which are subgroups of the $E_6$. The spontaneous breaking of $U(1)_{\psi'}$ symmetry governs the dynamics of a singlet fermion, which naturally serves as a freeze-in dark matter candidate. The dark matter mass arises from dimension-five operators, and a discrete symmetry ensures their stability. We show that freeze-in production from scalar decay can yield the correct relic abundance for dark matter masses between a few MeV to a few hundred GeV. Simultaneously, heavy right-handed neutrinos generate light neutrino masses via the type-I seesaw and produce the observed baryon asymmetry via leptogenesis.