Papers

We present simulations of the microwave sky at 2.5 arcsecond resolution over 100 square degrees, generated from existing full-sky, lower-resolution simulations, and use them to demonstrate extragalactic foreground removal for a CMB-HD survey. Our cleaning method detects and removes the cosmic infrared background and radio galaxies, yielding source catalogs that are 95% complete down to flux limits of 0.008 and 0.04 mJy at 90 and 148 GHz, respectively. It also identifies and removes galaxy clusters via the thermal Sunyaev-Zel'dovich effect, producing a cluster sample that is nearly complete above $M_{500c} = 5 \times 10^{13} M_\odot$. After cleaning, the residual foreground-plus-noise power spectrum of the coadded 90 and 148 GHz temperature map is 50% higher than previous idealized estimates, increasing cosmological parameter uncertainties for an 11-parameter $\Lambda$CDM + $N_\mathrm{eff}$ + $\sum m_\nu$ + $T_\mathrm{AGN}$ + $A_\mathrm{kSZ}$ + $n_\mathrm{kSZ}$ model by less than 7%. The small impact on cosmological parameters reflects the strong constraining power of CMB polarization, the temperature-polarization cross spectra, and CMB lensing spectra reconstructed from polarization-only estimators, all of which are minimally affected by extragalactic foregrounds. In particular, the survey remains a sensitive probe of light thermal relic particles, achieving $\sigma(N_\mathrm{eff}) = 0.0154$, which can exclude any new species ($\Delta N_\mathrm{eff} \ge 0.027$) with at least 90% confidence. Our simulation and foreground-removal codes are publicly available and should aid the development of analysis pipelines for ultradeep, ultrahigh-resolution microwave surveys.

The "Little red dots" (LRDs) are a population of accreting supermassive black holes (SMBHs) in the early Universe which often exhibit undermassive or even undetectable stellar hosts. Their early emergence, high space density, and extremely large black-hole-to-stellar mass ratios pose a serious challenge to conventional seeding scenarios that rely on baryon for both the formation and growth of black holes. Here we demonstrate that the above anomalies can be naturally resolved if dark matter is self-interacting. We apply a fully relativistic, non-equilibrium halo-evolution model, first developed in our earlier work, to trace the complete gravothermal evolution of self-interacting dark matter (SIDM) halos, from the initial collapse into BH seeds to the subsequent accretion of dark matter. We find that in highly concentrated halos assembled before reionization, gravothermal collapse efficiently produces stellar-mass black-hole seeds within a few hundred million years. Remarkably, and contrary to standard expectations for dark-matter accretion, heat conduction in SIDM then sustains a prolonged super-Bondi inflow that drives these seeds to supermassive scale by the LRD epoch, without baryonic assistance. The halo conditions required for completing these processes, together with the probability of avoiding major mergers that disrupt gravothermal evolution, result in an SMBH population consistent with the observed abundance and redshift distribution of LRDs. Our findings establish a pathway in which SMBHs are seeded and assembled primarily from dark matter, well before substantial galaxies form around them, thereby offering both a compelling physical explanation for LRDs and a new observational probe of dark-matter microphysics.

We introduce a model of dark matter charged under a new confining gauge group akin to the Standard Model's quantum chromodynamics (QCD). The fermionic sector contains both heavy and light quark-like fermions charged under the dark QCD. As the Universe cools, the gauge force confines and dark radiation transitions from free particles to low-mass dark pions. This form of dark radiation behaves as an imperfect fluid with nonzero viscosity. If the confinement occurs at the MeV scale, the model can alleviate the tension between early- and late-time measurements of the Hubble parameter, with a Bayesian posterior of $H_0 = 71.60\pm 0.63$ km/s/Mpc and a $3\sigma$ residual tension across Planck, ACT, BBN, Pantheon+, DESI DR2, and SH0ES. In the late Universe, the heavy dark matter can evade existing constraints on dark matter self-interactions while containing a rich phenomenology that may be accessible in future observations.

The cosmological constant $\Lambda$, and its generalization as Dark Energy, makes up roughly 70% of the present-epoch universe and has been part of the dominant $\Lambda$CDM cosmological paradigm for more than a quarter of a century. It has been put on trial by recent measurements of galaxy clustering, Type Ia supernovae, and the cosmic microwave background. In this Perspective we evaluate the strength of the evidence for a constant $\Lambda$, versus Dynamical Dark Energy (DDE) in the light of recent developments. We also assess what the future holds both for our Universe and our understanding of its most mysterious component.

We investigate how neutrino-mass constraints from cosmology depend on the assumed thermal history of the universe. Photon injection after neutrino decoupling would decrease the neutrino abundance inferred from the temperature of the cosmic microwave background (CMB), potentially loosening the upper limit on their masses. We first evaluate how the cosmological neutrino-mass bound is altered by the decay of massive particles into photons and dark radiation after Big Bang nucleosynthesis. To test the robustness of this constraint more generally, we also assess the impact of varying the temperature of the cosmic neutrino background without assuming a specific physical mechanism. We perform MCMC analyses of both frameworks with primary CMB observations from Planck, CMB lensing measurements from Planck and ACT, and baryon acoustic oscillation data from DESI. For the degenerate mass ordering, the $95\%$ credible limit tightens from $\sum m_\nu<0.0691$ eV in a standard thermal history to $\sum m_\nu<0.0652$ eV if the massive particles decay only into photons, while the injection of dark radiation in addition to photons slightly relaxes this limit to $\sum m_\nu<0.0710$ eV. The same pattern holds for the normal and inverted orderings, and the decay scenario shifts the bound on the sum of the neutrino masses by at most $0.004$ eV for fixed mass splittings. Allowing model-agnostic changes in the neutrino-to-photon ratio yields a $95\%$ credible limit of $\sum m_\nu<0.0724$ eV for the degenerate ordering, indicating that the stringency of our neutrino-mass bounds is not driven by constraints on the decay scenario. We find that the neutrino temperature and the sum of the neutrino masses are positively correlated, which implies that reducing the pre-recombination radiation density will only worsen the emerging tension between cosmological bounds on the neutrino masses and the measured mass splittings.

We confront a class of dark-energy equations of state emerging from group field theory (GFT) quantum gravity with DESI Data Release 2 baryon acoustic oscillations and Pantheon+ type-Ia supernovae. We introduce sampling parametrisations that replace microscopic initial-condition parameters by combinations more directly measured by background probes. The GFT solutions separate into logarithmic, power-law and oscillatory branches, determined by the microscopic interaction parameter $m$. The logarithmic branch is constrained to lie extremely close to a cosmological constant, while the power-law branch permits a small phantom deviation. Without perturbative-theory priors, oscillatory solutions can reproduce the mild preference of the distance data for a dip in $w(z)$ near $z\simeq0.5$–$1$. Profile-likelihood constraints favour $m\sim-2$ and $m\sim-5$ from BAO and supernovae, shifting towards $m\sim-3.5$ when CMB information is included. Conservative perturbative priors strongly suppress these deviations from $\Lambda$CDM. The quantum-gravity scale $z_q$, related to the average number of quantum gravity atoms, remains unconstrained, although its role in the time evolution makes higher-redshift observations a promising route to probing it. We further find that strong projection effects highlight the importance of performing likelihood profiling alongside our marginal posterior constraints. Our results provide a first direct test of GFT-motivated dynamical dark energy and demonstrate the potential for cosmological observations to inform quantum-gravity model building.

The ultraviolet continuum from a high-redshift accreting black hole is routinely used to infer its mass. GN-z11 offers a sharp test of whether that inference is secure. A thin-disk fit to its continuum gives an Eddington mass of $M_{\rm E}=1.12\times10^{7}\,\Msun$, an order of magnitude above the $\log(\Mbh/\Msun)=6.2\pm0.3$ inferred from broad N~iv. We show that this apparent tension is not imposed by the ultraviolet slope alone. We derive a closed-form criterion comparing the radius where observed photons are produced with the radius where a supercritical flow departs from efficient solution. For GN-z11, the criterion places the supercritical transition inside the ultraviolet-emitting region. Composite disks retaining the outer solution confirm this numerically, requiring $\Mdot_{\rm out}=3.1$–$4.4\,\Msun\,\mathrm{yr^{-1}}$ and placing $\Rsph=(1.2$–$1.8)\times10^{3}\,\rg$ inside $\Ruv=(2.6$–$2.9)\times10^{3}\,\rg$. Replacing the inner $\sim10^{3}\rg$ by a supercritical flow shifts the fitted slope by only $\Delta\betaUV=+0.08$ to $+0.11$, comparable to uncertainties from standard spectral modelling. Reversing the radial ordering requires the effective transition to move outward by a factor $1.7$–$2.1$. Within our composite models, the measured slope constrains black-hole mass only to around $10^{6}\,\Msun$, an order of magnitude below the continuum Eddington mass. Continuum Eddington masses should therefore be treated as model-dependent bounds that assume global thin-disk efficiency. The same hierarchy appears in three additional JWST sources under the disk-dominated interpretation, where separation masses are only $4$–$6$% of continuum Eddington masses. This distinction can substantially weaken the seed-mass and early-growth demands inferred from ultraviolet continua at cosmic dawn for the first massive black holes.

Thermal conduction can act as a heating mechanism for cool-core clusters, yet the same transport must erode any temperature discontinuity it crosses. Ophiuchus permits both effects to be tested in the same atmosphere. Its temperature rises from about $1$ keV in the innermost core to $9$ keV at $r\sim30$ kpc, making conductive heating especially favorable, while two sharply resolved cold fronts independently constrain the same transport. \textcolor{black}{We construct a steady radial energy balance from deprojected Chandra density and temperature profiles, with XRISM-based limits on turbulent heating and a subsonic inflow contribution.} Closing the balance with conduction requires only sub-Spitzer coefficients — \Ftwentyfive\ at 25 kpc and \Finner\ and \Fouter\ at the inner and outer cold fronts — so conduction can energetically supply the missing heat. Yet the same coefficients, applied normal to the cold-front surfaces, would broaden the fronts to their observed width limits within \Tinner\ and \Touter\ Myr, short compared to the characteristic sloshing timescales that generate and sustain such interfaces. For a representative 200 Myr residence time, the radial conductivity demanded by the energy balance exceeds the cross-front conductivity permitted by front survival by roughly a factor of three, with the required anisotropy exceeding unity in 99.9% and 98.8% of accepted models. The result leaves two broad regimes. Thermal conduction is either subdominant in Ophiuchus or contributes substantially but is strongly anisotropic, with heat transport suppressed across the cold fronts relative to the radial direction. Together they disfavor locally isotropic conduction as a dominant heating mechanism in cool cores generally.

We present the first application of neural-network-based baryon acoustic oscillation (BAO) reconstruction to real galaxy survey data, restoring the acoustic signature damped by nonlinear structure growth. {NERV} (Neural-network Enhanced Reconstruction of the UniVerse) augments standard reconstruction with a convolutional neural network trained on cubic $N$-body simulations, and explicitly accounts for realistic observational effects including the curved-sky geometry, the redshift-dependent selection function, and finite survey boundaries, by tessellating the survey volume into local patches. We validate the method on the MultiDark-Patchy mock catalogs, recovering unbiased BAO dilation parameters. Applied to the BOSS DR12 galaxy sample, NERV improves the precision of the BAO distance measurements significantly. These results establish neural reconstruction as a practical component of BAO analyses for ongoing surveys such as DESI, with the potential to substantially tighten constraints on the cosmic expansion history and the nature of dark energy.

When propagating through the universe, gamma rays at very-high energy (VHE, E > 100 GeV) and ultra-high energy cosmic rays (UHECRs, E > 1 EeV) can interact with the optical, infrared and microwave photon fields that permeate the universe. These interactions result in a characteristic absorption imprint in the spectra of extragalactic gamma-ray sources at VHE, and in a change in the mass-composition of UHECRs. The study of both VHE gamma rays and UHECRs therefore requires precise knowledge of the intensity of the cosmic photon fields. In this work, we explore the impact of the current uncertainties on the optical and infrared photon fields on the propagation of astroparticles. We restrict the range of available models to those best matching the recent measurements, and compare the different reconstructions resulting from these models. We find that the knowledge on the cosmic background light is no longer the dominant source of uncertainties in understanding the phenomenology of astroparticle sources in the low redshift universe (z < 0.1), enabling robust spectral and composition inference with the next generation of both gamma-ray and UHECR measurements.

We study a curvature-coupled dark energy model that can modify cosmological evolution both before recombination and during the late-time accelerated era. The model belongs to the class of scalar-tensor theories, in which a quintessence field is non-minimally coupled to the Ricci scalar. We specify the model through a shifted quartic coupling, $f(\varphi)=\alpha(\varphi^2-\varphi_{\rm today}^2)^2$, and an inverse power-law potential, $V(\varphi)=\Lambda\varphi^{-\sigma}$, where $\varphi_{\rm today}$ is a constant fixed by requiring the effective Planck mass to recover its present-day normalization. We implement the model in a modified version of $\mathtt{hi\_class}$ and compute its background and linear cosmological predictions. This specific form of non-minimal coupling allows an effective crossing of the phantom divide at late times while naturally suppressing deviations from standard gravity today and satisfying local gravity constraints. At the same time, the scalar field can modify the expansion history before recombination, shifting the acoustic scale in the direction required to alleviate the $H_0$ tension, while remaining dynamically relevant as dark energy at low redshift. For the parameter choices studied here, we find tens-of-percent deviations from $\Lambda$CDM in the expansion history, matter clustering, and metric-potential spectra. The modified evolution of the gravitational potentials leaves characteristic signatures in relativistic observables, with weak-lensing power suppressed by $\mathcal{O}(20$-$40\%)$ at low multipoles and order-unity changes in the late Integrated Sachs-Wolfe signal. Together, these results reveal a broad, scale-dependent phenomenology that motivates both a full parameter-space analysis and an extension of the present linear treatment to a dedicated non-linear $N$-body implementation.

The kinetic Sunyaev-Zeldovich (kSZ) effect is a sensitive probe of the free-electron momentum field, tapping into spatial correlations relevant for both large-scale cosmology and smaller-scale astrophysics. Precision kSZ measurements will provide a robust channel for constraining inflationary physics and equivalence principle violation, while simultaneously characterizing the small-scale gas distribution and associated baryonic feedback. However, robust kSZ measurement for cosmological applications is challenging, as the primary source of large-scale information is convolved with unknown small-scale dynamics in a degenerate combination. In this paper we address this challenge using the Effective Field Theory (EFT) framework that provides a principled way of treating the effects of small-scale dynamics on large scale-structure. Extracting the kSZ signal can be reformulated as a measurement of the one-loop galaxy-galaxy-electron momentum bispectrum. We compute here a pure dark matter version of this bispectrum in EFT and compare it to simulations. We derive the EFT counterterms that capture small-scale backreaction and find that two of them acquire a transverse direction. One transverse counterterm gets generated entirely by the vorticity of the velocity field. Our computation of the bispectrum dipole agrees with simulation data to 5% up to k_max = 0.23 h/Mpc, which is significantly larger than the reach of tree level theory, which breaks down beyond k_max = 0.07 h/Mpc (at redshift z = 0.5). For the bispectrum computed with the projected momentum field, relevant to measuring the kSZ field, we find a similar k_max, and make the first detection of the transverse counterterms from simulations. Our results suggest that one-loop bispectrum computations can play an important role in advancing kSZ cosmology.

Many well-motivated dark matter (DM) scenarios, including dark photon DM, axion-like particles, sterile neutrinos, and primordial black holes, can deposit energy into ordinary matter and significantly raise the thermal energy of astrophysical gas. We present a simple, model-independent implementation of DM-induced heating in the RAMSES hydrodynamic simulation code, parameterized by an effective local heat-transfer rate $\Gamma_{\rm heat}$. Using simulations of an isolated gas-rich dwarf galaxy with properties similar to WLM, we find that DM heating substantially alters the thermal structure of the circumgalactic medium (CGM) while leaving the star formation rate and the dense ISM nearly unchanged. For $\Gamma_{\rm heat} = 10^{-27}$ to $10^{-25}\ {\rm s}^{-1}$, the neutral hydrogen column at 10 kpc is suppressed by 0.6 to 3.9 dex relative to the control simulation. Similarly, the C II and Si II column densities decrease by 4.1 and 5.5 dex, and by 2.2 and 3.0 dex for C IV and Si IV. These signatures are accessible to quasar absorption-line spectroscopy, so CGM observations of nearby field dwarf galaxies provide a promising avenue for constraining dark matter models that inject heat into ordinary matter.

The recent observation of a single event at nuclear recoil energy $248\pm23\pm23$ keV by the LZ collaboration has prompted dark matter (DM) model builders to explore new physics explanations. While most of the existing literature invokes an endothermic scattering by a TeV-scale DM to explain the event, the favored DM interaction strength is in tension with the indirect detection constraints for most of the models. To evade this, this work explores exothermic DM scattering in a hadrophilic $q_1-q_2$ gauge extension setup with a new $Z'$ boson as a mediator to the dark sector. This can naturally escape the indirect searches in the absence of tree-level coupling of the mediator to neutrinos, $W,\tau$, and $b$ quark. We consider freeze-in production at a low reheating temperature to populate DM in the early universe, with a comparatively smaller coupling than the thermal freeze-out scenario. We identify the $Z'$ parameter space that can explain the observed relic density and LZ event simultaneously. We also discuss the relevant constraints coming from flavor mixing and FCNC.

Photometric galaxy clustering is a major probe of large-scale structure whose interpretation relies on accurate characterization of the galaxy redshift distribution, $n(z)$. Standard analyses assume the $n(z)$ of a selected lens sample to be isotropic, whereas observational systematics can imprint spatial variations. We quantify how this anisotropy biases angular galaxy clustering, internal consistency, and cosmological inference. We develop an analytical formalism for clustering under spatially varying $n(z)$ and forward-model realistic observational systematics, galaxy selection, and tomographic binning. For a fiducial magnitude-limited sample, modeling $w(\theta)$ with a single global $\bar n(z)$ underestimates the clustering signal by up to $\sim10\%$ at small scales in the highest-redshift bin, decreasing towards lower redshift. In a $2\times2$pt analysis, cosmological parameter shifts remain below $1\sigma$ for KiDS-like and DES-like setups. For an LSST-like setup retaining KiDS-like selection and observing-condition variation, $\sigma_8$ is biased low by $1.7\sigma$ and $\Omega_m$ high by $1.4\sigma$, while the high-redshift galaxy bias shifts by up to $\sim4\sigma$. The anisotropy also affects whether the probes are consistent enough to combine: a calibrated posterior-predictive internal-consistency test finds the contaminated angular clustering improbable given the joint cosmic-shear and galaxy–galaxy-lensing constraints (median calibrated $\tilde p\approx8\times10^{-4}$ under LSST-like noise). Anisotropic redshift distributions are therefore a potentially significant systematic for future photometric clustering analyses. Our forward model can be used to mitigate it given a realistic selection function and galaxy population model.

Dark-matter self-interactions can drive subhalos through gravothermal core collapse, but conventional lensing observables vary continuously with uncertain inner profiles. We show that gravitational-wave image multiplicity provides a discrete alternative. Near a macro-critical minimum, increasing subhalo central convergence opens two caustics in sequence, converting one detectable waveform into three and then five stationary images. Using weighted SASHIMI subhalos and an ET+CE selection model, we forecast almost 2 systems with at least three copies detectable in five years if dense remnants survive, compared with the prediction of 0.36 in 5yrs for standard CDM. The much more diagnostic second case is rare: the corresponding 5yrs yield with at least four detectable copies is 0.149, versus $1.55\times10^{-6}$ for CDM. According to the Poisson statistics, this gives only a 13.8% probability of at least one event detected in 5 years. Thus this is a rare-event, exposure-limited test rather than a guaranteed detection: a resolved second pair would be difficult to reconcile with the fiducial single-halo CDM population, whereas a five-year null result would not by itself constrain remnant survival strongly.

The cosmic web governs the growth of dark matter halos through anisotropic matter flows, yet the role of local velocity-field dynamics in halo mass assembly remains poorly understood. The velocity field carries non-linear information about matter transport that is not fully captured by density alone, particularly within filaments. We investigate this connection through velocity divergence and vorticity using dark matter halos in the IllustrisTNG TNG50-1-Dark simulation at $z=0$. Cosmic filaments are extracted with the 1-DREAM framework, and we statistically study how halo mass relates to divergence, vorticity, and radial distance from filament spines across three overdensity regimes. The velocity field correlates strongly with the cosmic web, with converging flows tracing filaments and vorticity arising primarily in overdense, non-linear regions. Within filaments, both divergence and vorticity show enhanced kinematic activity at intermediate distances from the spine, with amplitudes increasing with overdensity. We confirm a clear mass segregation, with more massive halos preferentially located near filament centers. Most notably, halos with masses $M \geq 10^{12} M_\odot$ are found almost exclusively in regions of low absolute divergence and low vorticity, while highly dynamical environments are dominated by low-mass halos. We further find that halos residing in kinematically calm regions assembled later and maintain higher recent accretion rates than those in dynamically active environments, linking present-day velocity-field properties to halo growth histories. These results demonstrate that the local kinematic state of the cosmic flow encodes information about halo mass and assembly beyond that contained in density alone, providing a complementary dynamical description of halo growth within the cosmic web.

Observations of the Large Scale Structure sit on the light cone. Cosmological statistics must then be projected from 3D to angular space. For the bispectrum and higher orders, the polyspectra crucially depend on the angles between Fourier modes in addition to their moduli. This leads to delicate cancellations in the squeezed limit and make Limber's approximation fail. This is problematic for predicting the Super-Sample Covariance (SSC), with no known full-sky computation of the angle-dependent terms. Here we develop a new analytical approach to the angular bispectrum, writing Fourier vectors as derivatives that we project on the sphere and tackle with spin-weighted spherical harmonics. We find bispectrum equations free of cancellations and amenable to beyond-Limber numerical methods. Limber's approximation becomes well-behaved and works in the squeezed limit and any other configuration. We then derive the consequence for the covariance between 1-point (e.g. cluster or peak counts) and 2-point statistics (e.g. galaxy clustering, weak lensing or 21cm), which relates to the squeezed bispectrum. The results split into intra-survey covariance and SSC, where we find significant differences with past literature. Our results are in line with the 2+1 geometry of the light cone, having e.g. temporal derivatives in addition to spatial ones. We also find that the tidal term is degenerate with the simpler growth-only term. The total SSC contains a wealth of new terms including wide-angle effects and responses to the background velocity. However we show numerically that it is well approximated by a single term with a density-only response. This ultimately makes our results simpler than past literature. This full-sky result should be more representative than the past flat-sky limit for large survey areas, thus being the one relevant for coming observations with e.g. Euclid, LSST and Roman.

Halo assembly bias describes the dependence of halo clustering on properties beyond halo mass and provides a means of connecting halo assembly history to the large-scale matter distribution. We develop and test a relative clustering methodology for measuring differences in the large-scale bias of two halo populations using the ratio of their cluster–galaxy pair counts. In the linear regime, this ratio can be related to the difference in cluster bias, while the galaxy autocorrelation independently constrains the galaxy bias. We apply the method to a volume-limited X-ray selected galaxy clusters from the eROSITA eRASS1 catalogue, dividing them into high- and low-temperature subsamples with matched distributions in X-ray luminosity, redshift, and exposure. We first verify the validity of the linear approximation, finding $\xi(r)\ll1.0$ over $4$–$10\,h^{-1}\mathrm{Mpc}$. The joint clustering analysis yields $b_{\rm g}=1.610^{+0.018}_{-0.019}$ and a bias difference $\Delta b=0.8\pm1.1$. DES Y3 weak gravitational lensing measurements provide an independent test of the halo masses, yielding $\log(M_{200\mathrm{m}}/h^{-1}M_\odot)=14.33^{+0.07}_{-0.07}$ and $14.33^{+0.10}_{-0.09}$ for the high- and low-temperature samples, respectively, with a mass difference of $\Delta\log M_{200\mathrm{m}}=-0.01\pm0.12$. We therefore find no statistically significant evidence for assembly bias at the $1\sigma$ level, but our measurement is consistent with the bias difference measured for analogous halo populations in the FLAMINGO simulation. The result therefore does not rule out assembly bias at the level expected from the simulation, while demonstrating that the relative clustering methodology provides a framework for future assembly-bias measurements with larger X-ray cluster and galaxy surveys.

We discuss the effectiveness of the recently proposed quantity $J_{s}$ in suppressing the contribution of a string's self-field to its spectrum, which is used in field theory simulations to track the emission of axions by decaying cosmic strings. We compute the contribution of the self-field to the spectrum of this quantity and to the usually computed spectrum of $\phi\partial_{t}\alpha$ for an infinitely long straight string. Although we demonstrate that the $J_s$ approach is a substantial improvement, we also point out that this highly symmetric model doesn't capture the full contribution of the self-field to the spectrum emitted by strings of a more complex shape, which are typical in network simulations. We then illustrate this point numerically using a simulation of a sinusoidally perturbed straight string. For this simple configuration, we managed to separate the contribution from the self-field and the radiation to the spectrum of $J_{s}$, using the self-field subtraction method. This allows us to show that the spectrum of $J_{s}$ is still dominated by the $n=1$ mode of the string's oscillation, which can be attributed to the self-field and is largely suppressed when the self-field is removed. We also demonstrate that this mode is predominantly sourced by variations along the direction parallel to the string, which is missing in the unconnected segment model, used by arXiv:2512.13653 to claim the effectiveness of $J_{s}$ in suppressing the self-field's contribution to the spectrum.

Temporal gaps in Swift-XRT X-ray afterglows of gamma-ray bursts (GRBs) limit how precisely the standard empirical plateau function constrains three parameters: the plateau end time $T_a$, the plateau flux $F_a$, and the post-plateau decay index $\alpha$. These parameters underpin the luminosity-time and fundamental-plane correlations of GRB plateaus so that smaller measurement uncertainties can improve the statistical power of GRBs as cosmological probes. As the fifth study in a series on light-curve reconstruction, this work benchmarks four models on a parent sample of 545 Swift-XRT GRBs: a Physics-Informed Neural Network (PINN) with seven afterglow-motivated priors, ReFANN, a Siamese dual-branch network with Monte Carlo dropout, and Polynomial Quantile Regression (PQR). All four reduce the fractional uncertainties in $\log T_a$, $\log F_a$, and $\alpha$ relative to the original observations. The PINN with a broken power-law prior gives the largest reductions ($\sim$41-49%), at a higher outlier rate ($\sim$15-20%), while ReFANN, Siamese network and PQR give more modest reductions ($\sim$19-27%) with outlier fractions $\lesssim$4%, across a wider range of light-curve morphologies. A reduced-$\chi^2$ prior-selection scheme recovers a morphological classification consistent with independent labelling, though it does not match the best single-prior reduction. Together, these results provide a systematic benchmark of physics-informed and data-driven reconstruction strategies for irregularly sampled GRB afterglows.

FOSSIL (FTS fOr CMB Spectral diStortIon expLoration) is a proposed ESA M8 mission tailored to measure the spectral distortions of the Cosmic Microwave Background (CMB) with a sensitivity three orders of magnitude beyond the COBE/FIRAS legacy measurement. Achieving this sensitivity demands an extraordinarily challenging cryogenic architecture: the scientific instrument must be maintained at the lowest achievable temperature (4.5 K), while the detector focal plane assembly operates at 50 mK. This paper presents the thermal architecture of the FOSSIL payload and instrument, from the spacecraft service module at 293 K down to the sub-kelvin detector stage. The thermal design draws on heritage from the Planck and ARIEL missions and relies on a staged passive cooling chain comprising a multi-layer insulation blanket, three V-groove radiators (operating at approximately 130 K, 90 K, and 50 K), and a 25 K actively cooled shield fed by an ESA-provided 4 K cryocooler chain. Sub-kelvin temperatures are achieved via a multi-stage adiabatic demagnetisation refrigerator (ADR) system developed for NewAthena/X-IFU, providing continuous cooling at 1.8 K and 350 mK, and 50 mK with an 80% duty cycle. The Focal Plane Assembly (FPA), which houses four Kinetic Inductance Detector (KID) arrays at 50 mK, is thermally isolated from the 4.5~K bench via a carbon-fibre reinforced polymer (CFRP) hexapod structure with staged heat interception. We present the steady-state thermal budget across all stages, demonstrating comfortable margins at every temperature level, and discuss key thermal design drivers on the Blackbody Internal Reference (BBIR).

We investigate whether the recently reported high-energy nuclear recoil event by the LUX-ZEPLIN (LZ) collaboration can be interpreted as a signal of freeze-in dark matter (DM) in a non-standard cosmological history. We consider DM production from the Standard Model thermal bath during an era of modified expansion preceding Big Bang nucleosynthesis. Requiring consistency with the observed DM abundance, cosmological and collider constraints, we identify viable parameter regions that can accommodate the LZ event. Our results demonstrate the potential of direct detection experiments to probe non-standard cosmological histories and, more broadly, alternative theories of gravity.

The LUX-ZEPLIN (LZ) experiment has reported a high-energy nuclear recoil candidate, LZ230616, at $E_R = 248 \pm 23_{\rm stat} \pm 23_{\rm sys}\,$keV, which is difficult to reconcile with elastic scattering of halo dark matter (DM). Endothermic inelastic scattering offers a natural explanation, but it rests on a nucleon coupling that is off-diagonal rather than diagonal. We show that this feature arises automatically within a radiative Dirac neutrino mass framework. The Standard Model is extended by three right-handed neutrinos (RHNs), a pair of vector-like neutral fermions for each generations, and a scalar sector comprising an inert doublet and a real singlet, governed by a $\mathcal{Z}_2 \times \mathcal{Z}_4$ symmetry. The $\mathcal{Z}_2$ symmetry stabilizes the DM, while the $\mathcal{Z}_4$ symmetry forbids the tree-level Dirac Yukawa coupling together with all renormalizable Majorana mass terms, thereby allowing the generation of Dirac neutrino masses at the one-loop level through a softly broken scalar trilinear coupling `$\kappa$'. Since the neutral dark-sector fields can be expressed in terms of real scalar mass eigenstates, the $Z$ boson couples to them purely off-diagonally and elastic $Z$-mediated scattering is absent identically rather than simply suppressed. The singlet–doublet mixing induced by `$\kappa$' governs both the inelastic rate and the one-loop Dirac neutrino mass, which therefore vanish as $\kappa \to 0$. We find that DM masses in the few hundred GeV to TeV range, with splittings of $\mathcal{O}(340-360)\,$ keV, simultaneously reproduce the observed relic abundance, account for the LZ event and yield neutrino masses of the correct order, while respecting elastic direct-detection limits and all relevant theoretical and experimental constraints. A substantial part of the surviving parameter space lies within reach of DARWIN.

We report new optical spectroscopic results for the gamma-ray blazar PKS~0903$-$57. The optical properties of this source have historically been complicated to establish, due to a bright foreground star separated from the blazar by only 0.67$^{\prime\prime}$. During its strong high-energy (HE) and very high-energy (VHE) gamma-ray flaring period between 2020 and 2022, the Southern African Large Telescope and Very Large Telescope (VLT) took spectra of its optical counterpart. Photometric monitoring with the Rapid Eye Mount (REM) telescope complemented the observations. The observational results were inconclusive. A major improvement was achieved with a repeat of the VLT observations in 2024 using a 0.5$^{\prime\prime}$ slit under excellent seeing conditions. This led to the detection of five narrow emission features in the spectra, yielding a robust redshift measurement of $z$ = 0.2621 $\pm$ 0.0006. We also detected a broad symmetric H$\alpha$ emission component with a full width at half maximum of 4020 $\pm$ 30 km/s in the spectra. Non-thermal jet emission was found to dominate the continuum, supporting the classification of PKS~0903$-$57 as a flat-spectrum radio quasar. Another interesting observation was a significant velocity displacement of $\sim$1500~km/s between the broad H$\alpha$ line and its corresponding narrow line. Such a displacement may be associated with unusual dynamics in the broad-line region or accretion flow, a system containing two supermassive black holes, or a recoiling black hole produced following a supermassive black-hole merger. Building on these results, we carried out additional VLT, REM, and Swift observations towards the end of 2024 and during 2025. We present in this paper the preliminary results from the analysis of these new data.

We forecast constraints on the parameters of the coupled dark energy model using DESI and Euclid data with the one-loop galaxy power spectrum. We investigate the distinguishability of our model from the zero-coupling scenario at the $1\sigma$ level and explore how the parameter constraints depend on the fiducial coupling $\beta_{\rm fid}$, the fixed potential slope parameter $\nu$, the maximum wavenumber $k_{\max}$, and different prior choices. We find that the inclusion of mildly non-linear scales improves the constraints on the coupling by roughly a factor of five. Then, we address the question of which is the minimum value of $\beta$ that can be distinguished from zero at $1\sigma$. We find that for our reference case $k_{\rm max}=0.2 h/$Mpc, $\beta=0.15$ lies roughly $1\sigma$ above zero. In the most optimistic case with $k_{\rm max}=0.3 h/$Mpc and including a Planck prior on $\Omega_{m0}$, this value can be reduced to $0.05$. These values are substantially larger than the current constraints on $\beta$, but the latter have been obtained assuming $\beta$ to be constant from at least the decoupling epoch to today, while we only employ late-time data. We conclude therefore that only models that allow for time-varying couplings can be detected with late-time datasets.

We present the first application of neural field-level baryon-acoustic oscillation (BAO) reconstruction to real spectroscopic survey data. We develop Linear Field Transformer (LiFT), a 3D vision transformer that takes as input the observed galaxy field, its standard reconstruction, and a set of context channels encoding local line of sight, survey coverage, and redshift, and learns to correct standard reconstruction toward the linear density field. We construct a forward-modeling pipeline to produce mock lightcones similar to the DESI Data Release 1 (DR1) luminous red galaxy (LRG) sample, and train LiFT on these. We validate LiFT on held-out simulations, as well as on additional mocks which differ in gravity solver, halo finder, HOD, cosmology, and fiber-assignment history, as well as on mocks analyzed with a distorted distance-redshift relation; ultimately, we find unbiased dilation parameters with consistently tighter constraints than standard reconstruction. Applied to the DESI DR1 LRGs, LiFT improves errors on $\alpha_{\rm iso}$ by 13%, 24%, and 30% and on $\alpha_{\rm AP}$ by 5%, 22%, and 30% relative to the DESI DR1 standard reconstruction analysis in the LRG1, LRG2, and LRG3 bins respectively; meanwhile, our DR1 central values remain consistent with DR1 and DR2. This equates to a factor of 1.2, 1.7 and 2.0 increase in Figure of Merit (or effective survey volume) if one were to only use standard reconstruction. Ultimately, these results establish LiFT as a validated, survey-ready tool for current and upcoming galaxy surveys.

A cosmological source can have identically zero comoving energy density and still gravitate through its pressure. This differs from an energy density passing through zero at an isolated time. Maintaining zero density with nonzero pressure requires energy exchange in a nonstatic FLRW universe. We construct proportional-pressure and constant-pressure backgrounds, including steady-state, bouncing, recollapsing, and $\Lambda$CDM expansion histories, and relate them to creation pressure, bulk viscosity, and time-dependent vacuum descriptions. Scalar fields can realize zero density along a trajectory; imposing it as an identity on a timelike $P(X,\phi)$ domain selects $P=A(\phi)\sqrt{X}$, locally a signed potential-free cuscuton for nonzero $A$. We formulate linear perturbations without dividing by the vanishing background density. For a nondegenerate $P(X,\phi)$ scalar with a timelike field gradient, zero density and negative pressure cannot coexist with positive kinetic and gradient coefficients when the interaction leaves its principal kinetic structure unchanged. Two models couple the scalar to number-conserving dust through a field-dependent mass. The proportional-pressure trajectory is a homogeneous saddle and has no growing power-law dust mode on its accelerating branch in the matter-sourced quasistatic subhorizon approximation. A distinct constant-pressure model reproduces the $\Lambda$CDM expansion history, but in the same approximation the dust mode that initially grows during matter domination reaches a maximum at low redshift and then decays. These restrictions apply to the stated scalar class and models. Exact Bianchi I solutions demonstrate separately conserved zero-density sources supported by anisotropic stress. An inverse-pressure construction distinguishes the density of an interacting constituent from an inferred dark-energy density that may cross zero.

The LUX-ZEPLIN experiment has reported a single nuclear recoil candidate at $E_R = 248 \pm 23(\mathrm{stat}) \pm 23(\mathrm{sys})$~keV, disfavouring the background-only hypothesis at a global significance of $2.6\sigma$. The difficulty such an event poses is not the recoil energy itself but the absence of any accompanying excess at low energy as elastic scattering of halo dark matter yields a monotonically falling spectrum, and supplying the required momentum transfer $q \simeq 246$~MeV already demands $m_\chi \gtrsim 79$~GeV. We confront the event with the two kinematically distinct mechanisms that evade this limitation, endothermic inelastic dark matter, in which a mass splitting $\delta \sim \mathcal{O}(100)$~keV forbids low-energy recoils and boosted dark matter, in which a light relativistic flux supplies the momentum, treating both with the same model-independent scalar–scalar, pseudoscalar–scalar and pseudoscalar–pseudoscalar effective operators. The two scenarios prove spectrally distinguishable. The inelastic spectra sit near the observed energy for every operator, placing only $6$–$17\%$ of events below 150~keV, whereas the boosted spectra depend critically on the operator: the scalar and pseudoscalar–scalar interactions place $99\%$ and $92\%$ of their events below 150~keV, while the pseudoscalar–pseudoscalar interaction places $75\%$ above it. Momentum dependence is thus essential to the boosted interpretation, but in the inelastic case it trades against the splitting, the preferred $\delta$ decreasing monotonically from $\mathcal{O}_{ss}$ to $\mathcal{O}_{ps}$ to $\mathcal{O}_{pp}$. Because the scenarios differ across the whole high-energy window, a handful of additional events would separate them, placing the question within reach of the full LZ exposure.

The tip of the red giant branch (TRGB) is a relative distance measure that can be used, as part of a cosmic distance ladder, to infer the Hubble constant. We present a Bayesian framework for TRGB inference that incorporates selection effects and propagates per-star flux uncertainties, modelling the observed catalogue as a realisation of an inhomogeneous Poisson point process which includes the RGB and contamination from the asymptotic giant branch (AGB). We derive expressions for the predicted TRGB uncertainty in terms of the photometric noise, stellar density fainter than the tip, and the ratio of the AGB to AGB+RGB number densities, $r$. We validate these results using synthetic data from a simplified luminosity function. Applying the full model to synthetic catalogues, we find that the dependence on $r$ is the strongest of the three, and that reliable inference requires $r\lesssim0.5$; otherwise, the posterior can be multimodal, and the TRGB cannot reliably be identified. We then fit data from 11 Hubble Space Telescope fields of the calibrator galaxy NGC 4258. We find intrinsic field-to-field scatter of $\tau=0.067^{+0.023}_{-0.017}$ mag (beyond the measurement noise), and a galaxy-level, extinction-corrected TRGB magnitude of $25.324 \pm 0.026$ in the F814W filter. Combining this with the maser distance modulus to NGC 4258 of $29.397 \pm 0.032$ mag, we determine the absolute magnitude of the TRGB to be $M_{\mathrm{F814W}}^{\mathrm{TRGB}} = -4.073 \pm 0.026$ (tip) $\pm 0.032$ (dist), consistent with, but brighter than, other recent analyses of the same data.