Papers

Pre-inflationary axion cosmologies admit a broad range of initial velocities and angles that generate a landscape of misalignment histories, whose nonlinear fate and small-scale structure remain uncertain. We use lattice simulations to follow the QCD axion in the kinetic misalignment scenario through the QCD crossover into the non-relativistic regime, towards matter-radiation equality. In the meV mass range, this landscape contains a continuum of ”summit-landing" solutions sensitive to both the initial angle and velocity. Nonlinear fragmentation reduces the comoving axion number and modifies the relation between the misalignment history and present relic abundance. For sufficiently large axion masses, strong self-interactions drive the small-scale spectrum towards the attractor found for post-inflationary axions, so that distinct early histories can lead to similar late-time structures.

We apply the two-particle irreducible (2PI) formalism as a framework for a consistent thermal resummation in studies of cosmological phase transitions. Considering a model with two mixing real scalar fields, we work within the Hartree approximation and renormalize the 2PI effective action, while introducing a connection to physical parameters. This yields the Hartree-resummed finite-temperature effective potential, which is valid for all temperatures and avoids the limitations of conventional methods based on the high-temperature approximation. With this potential, we study one- and two-step transitions within the model, and compare our results with those obtained using resummation schemes widely employed in the literature. Finally, we evaluate the gravitational wave spectrum generated from first-order phase transitions, demonstrating the impact of the choice of resummation scheme on the predicted spectrum.

We investigate an analytically tractable sector of the Extended Cuscuton model, a non-dynamical dark-energy realization within the framework of viable Horndeski gravity. We focus on four benchmark submodels and constrain them with current background probes, namely cosmic chronometers, Type-Ia supernovae, and BAO, while imposing theoretical viability, Lunar Laser Ranging, and Big Bang Nucleosynthesis bounds. We then forecast third-generation bright-standard-siren constraints with Einstein Telescope and Cosmic Explorer networks, considering prompt-emission, afterglow, and kilonova counterparts. Current data already restrict the viable parameter space to small departures from $\Lambda$CDM and do not remove the calibration-driven offset between the CC+SN and CC+BAO determinations of $H_0$. In principle, future bright sirens substantially sharpen the constraints, especially for kilonova catalogues and extended detector networks. Across the forecast configurations, the relative uncertainty on $H_0$ remains below $13.18\%$ and can reach $0.21\%$ in the most constraining cases, while $\Omega_\Lambda$ is recovered at the percent level in the best cases. These results show that third-generation standard sirens can provide a precise complementary test of non-dynamical dark energy beyond $\Lambda$CDM.

The next generation of gravitational wave detectors, such as the Einstein Telescope (ET), will observe orders of magnitude more binary black hole mergers than current facilities. Most of these events will lack an electromagnetic counterpart, also known as dark siren events, yet will still enable percent-level cosmological constraints. However, the likelihood traditionally used in Hierarchical Bayesian Inference (HBI) for population-level analyses becomes computationally prohibitive as the size of dark siren catalogues and population parameters grow. In this work we compare HBI against simulation-based inference (SBI) as a scalable alternative for cosmological population inference in the ET era. Studying a proof-of-concept example of a mock ET inference, we build a catalogue of $O(10^4)$ binary black hole events, then perform inference on the Hubble constant $H_0$ and matter density $\Omega_m$ in a flat $\Lambda$CDM cosmology, using both a hierarchical analytical likelihood and Marginal Neural Ratio Estimation (MNRE). We find excellent agreement between the two approaches, with SBI reproducing the HBI posteriors to high accuracy, while requiring orders of magnitude less computation once the simulation and training cost is amortized. We further demonstrate that SBI extends straightforwardly to a joint cosmology-plus-astrophysics analysis, simultaneously constraining $(H_0,\Omega_m)$ together with the parameters of the star formation rate density, at negligible additional cost compared to the significant increase in complexity such an extension would require within the HBI framework. Our results indicate that SBI is a promising and scalable tool for population inference with third-generation GW detectors.

$\Lambda$CDM provides a remarkably successful description of the Universe, yet it suffers from persistent discrepancies, with one of the most prominent being the $\sigma_8$ tension between early-Universe Cosmic Microwave Background (CMB) measurements and late-time structure growth probes. To address this mismatch, we investigate generalized interacting dark sector models governed by a quintessence field coupled to dark matter via conformal and disformal transformations. By generically engineering a quintessence potential that ensures the combined evolution of dark matter and dark energy exactly mimics $\Lambda$CDM, our models satisfy all background constraints while signaling modified gravity effects within the linear perturbation sector. We perform a statistical analysis using late-time structure growth data anchored to Planck 2018 constraints. Furthermore, we address the initial conditions problem inherent to purely disformal models, demonstrating the necessity of a non-zero initial kinetic energy to initiate dynamical evolution. Our results indicate that purely conformal couplings efficiently transfer energy from dark matter to the quintessence field, suppressing late-time structure formation and successfully accommodating CMB data within late-time modified structure growth. This purely conformal limit yields a strong statistical preference over the $\Lambda$CDM baseline according to both the Akaike and Bayesian Information Criteria. Conversely, while non-zero disformal couplings introduce a kinetic friction that dampens the dark sector energy exchange, their expanded parameter space incurs heavy statistical penalties without providing a decisive improvement over the purely conformal scenario.

The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) has begun its ten-year survey of the entire visible southern hemisphere. To ensure robust cosmological measurements, computationally inexpensive investigations of modeling choices must be made to gauge the performance of proposed cosmological analyses. In this paper, we introduce the $\texttt{Augur}$ tool of the Dark Energy Science Collaboration (DESC), which provides Fisher forecasts for cosmological inference for the LSST using software frameworks designed for DESC science. We test the pipeline by comparing it to forecasts produced by external code and direct sampling of the posterior via nested sampling methods, finding good agreement between all methods. We additionally investigate a range of modeling and hyperparameter choices for a 3$\times$2pt investigation in harmonic space, providing users with diagnostics to obtain reliable forecasts. $\texttt{Augur}$ will be continually updated to be compatible with the other tools in the DESC software ecosystem as additional probes and functionality become available.

Gravitational-wave (GW) cosmology is currently bottlenecked by the scarcity of electromagnetic counterparts for bright sirens and the systematic uncertainties of galaxy catalogs for dark sirens. We propose a purely gravitational resolution using the Integrated Cosmological Memory (ICM)-the cumulative GW strain encoded in the spacetime geometry of an expanding Universe. While the GW transient emitted by the source provides the luminosity distance, the ICM accumulates a mathematically distinct integral of the cosmic expansion history. We demonstrate that extracting both observables from a single binary merger completely breaks the distance-redshift degeneracy within the gravitational sector. This establishes a novel, catalog-free dark siren framework for third-generation GW detector networks. Crucially, the resulting constraints on late-time dark energy are only weakly sensitive to the local expansion rate, providing a robust cosmological probe that can potentially mitigate the impact of the H0 tension.

Departures from standard slow-roll (SR) inflation have attracted increasing interest in recent years. In particular, scenarios that strongly enhance the power spectrum of the curvature perturbation are often proposed as a mechanism for producing primordial black holes, which could account for part or the totality of dark matter. An ongoing debate in these models is whether non-linear interactions of cosmological perturbations on small scales, namely loop corrections, can become sufficiently large to backreact on CMB scales. In this work, we adopt a non-linear framework, dropping spatial-gradient terms, to study loop corrections at super-Hubble scales. Our focus is to clarify the role of spatial-diffeomorphism invariance, especially dilatation invariance, which is the symmetry under overall rescaling of spatial coordinates, in demonstrating the suppression of loop corrections to the CMB power spectrum. This analysis is valid for any inflaton potential, both with smooth and sharp transitions, and at any loop order, as long as the CMB scales cross the horizon during an SR phase. We compare this result with the different explanations proposed in the literature using various gauge choices, and show our analysis is consistent both in the comoving and in the spatially-flat gauges. In particular, we show that in the spatially-flat gauge, the cubic interaction contains a term that diverges in the sharp transition limit. However, this contribution is exactly cancelled by the quartic interaction.

The Large-Scale Structure (LSS) of the Universe forms a complex network of nodes, filaments, sheets, and voids known as the Cosmic Web. As current and upcoming galaxy surveys increasingly probe the quasi-linear and non-linear regimes of structure formation, understanding its geometry and dynamics is essential for precision cosmology. In particular, cosmic filaments, which channel matter across the web, are central to these dynamical processes. This thesis develops numerical tools to study the non-linear cosmic web. First, the Sahyadri suite of high-resolution cosmological $N$-body simulations is introduced, providing a framework for precision studies of LSS and its cosmological dependence. A calibration framework for filament reconstruction is then developed using controlled filament realizations, enabling systematic investigation of reconstruction biases. The effects of filament curvature and reconstruction noise on inferred filament properties are quantified, and a novel Fourier-space smoothing approach is introduced to improve profile recovery. The thesis further presents Skeletor, a Voronoi-based filament finder that identifies filamentary structures directly from discrete tracers while explicitly incorporating the hierarchical nature of the cosmic web. A novel framework for classifying sub-filamentary structure is developed. Applying these tools to cosmological simulations reveals distinct properties of filament substructure and provides a detailed view of filament phase space, including coherent inflows, multistreaming, and caustic-like features. Together, these developments provide a framework for studying the geometry, hierarchy, and dynamics of the non-linear cosmic web. More broadly, they contribute to the ongoing effort to build a physically motivated understanding of the non-linear cosmic web beyond traditional measures of clustering.

Recent James Webb Space Telescope (JWST) observations favor an earlier and more efficient reionization history, leading to Thomson-scattering optical depths larger than those inferred from the cosmic microwave background (CMB). We investigate whether this tension can be alleviated within the framework of self-interaction (SI) gravity by modifying the cosmological expansion history while retaining the standard astrophysical description of reionization. The SI gravity parameters are constrained through a Bayesian MCMC analysis of the Union3 Type Ia supernova and DESI DR2 Baryon Acoustic Oscillations (BAO) data. The resulting expansion history is then used to predict the ionization history and Thomson optical depth. We find that the predicted optical depth decreases from $\tau_{\rm{CMB}}\simeq0.076$ in $\Lambda$CDM to $\tau_{\rm{CMB}}\simeq0.061$, consistent with the Planck PR4 measurement within $1\sigma$, thereby substantially reducing the optical-depth tension.

Recent James Webb Space Telescope (JWST) observations suggest an unexpectedly high abundance of luminous galaxies at $z \gtrsim 10$, challenging traditional models of star formation. The 21-cm signal of neutral hydrogen offers a complementary, volume-averaged probe of these early epochs. In this work, we present a major extension to the explicitly photon-conserving semi-numerical framework SCRIPT, enabling the self-consistent calculation of spin temperature fluctuations driven by inhomogeneous Lyman-$\alpha$ coupling and X-ray heating during the Cosmic Dawn. Using this framework, we compare a JWST-informed fiducial galaxy model, featuring a redshift-evolving star-formation efficiency, against a static baseline model. We demonstrate that the global 21-cm signal is highly degenerate: a static galaxy population can replicate the absorption trough of an evolving population if a steep, mass-dependent X-ray efficiency is invoked to shift intergalactic medium heating toward highly abundant, low-mass halos. However, this degeneracy is definitively broken by the 21-cm power spectrum. Because spatial fluctuations retain a structural memory of the clustering bias of the source halos, the tuned static model predicts a substantially lower power spectrum amplitude across all characteristic peaks. We conclude that combining global and fluctuation measurements is essential to robustly constrain the redshift evolution of high-redshift galaxies.

Redshift information is central to nearly every extragalactic and cosmological application of sky surveys, yet only a small fraction of cataloged sources, and none of the photons forming diffuse backgrounds, have spectroscopic redshifts. Clustering-based redshift inference estimates the redshift distribution of an arbitrary dataset through spatial cross-correlation with a reference sample of known redshifts. It relies only on positional information, making it applicable to any tracer of large-scale structure, including source populations and diffuse intensity maps. Its broader adoption, however, has been limited by technical barriers: assembling and characterizing spectroscopic references, expensive pair-counting computations, theoretical corrections, and systematic control. To remove these barriers, we introduce Tomographer, an end-to-end clustering-redshift framework. The key design feature is a set of precomputed "activation maps" encoding the spatial pair information of 3 million SDSS spectroscopic galaxies and quasars up to $z\sim4$ over $10{,}000\,{\rm deg}^2$. This eliminates user-end pair counting, reducing computational scaling from $O(N\log N)$ to $O(1)$ map multiplications. Given a source catalog or intensity map, Tomographer returns the bias-weighted redshift distribution, $b(z)\,{\rm d}N/{\rm d}z(z)$ or $b(z)\,{\rm d}I/{\rm d}z(z)$. We validate the framework against samples with known redshifts, demonstrate accurate uncertainty estimates, and show robustness to survey footprint, spatially varying selection functions, beam smoothing, and foreground contamination. We showcase applications to source catalogs selected by flux, color, photometric redshift, morphology, or variability, and to intensity maps from radio to X-rays. Future Tomographer releases will incorporate additional wide-field spectroscopic reference samples as they become available.

These lecture notes provide an introduction to the generation of initial conditions for cosmological N-body simulations. Starting from the definition and properties of Gaussian random fields, we discuss their role in cosmology and the efficient generation of such fields using Fourier methods. The Vlasov-Poisson system is introduced as the governing framework for cold collisionless matter, and its solution via characteristics and Lagrangian perturbation theory (LPT) is detailed. We discuss the use of LPT for initializing N-body simulations, emphasizing the importance of high-order LPT and late-time starts to minimize truncation and discreteness errors. Finally, we discuss time integration schemes, including PT-informed integrators, and their role in accurately evolving the system. These notes aim to bridge the gap between theoretical perturbation methods and practical simulation techniques.

Fast radio burst (FRB) dispersion measures trace free electron column densities and will soon offer a unique probe of low-$z$ baryons. Such measurements can constrain cosmology directly, with the free electrons serving as a new tracer of large-scale structure, and probe the baryonic feedback of galaxies, a leading systematic for weak lensing surveys such as LSST and Euclid. We prepare for both science cases using effective field theory (EFT) and hydrodynamical simulations. We construct the one-loop EFT description of the free-electron auto-spectrum $P_{ee}$ and electron-galaxy cross-spectrum $P_{eg}$, placing FRB dispersion clustering on the same theoretical footing as spectroscopic galaxy analyses, and quantify the FRB densities at which this modeling is useful. We also investigate suitable galaxy samples for cross-correlations, finding that current spectroscopic catalogs provide appropriate redshift range and sufficient density. We validate the model against the FLAMINGO simulations, jointly fitting $P_{ee}$, $P_{eg}$, and $P_{gg}$ for DESI-like samples at $z=0.2$ and 0.5. The model describes all three spectra to $k\sim0.2\,h\,{\rm Mpc}^{-1}$, with an electron linear bias $b_{e,1}\simeq0.92$, higher-order biases consistent with zero, and all parameters stable across feedback variants. Together with the near-perfect electron-matter correlation $r_{em}\simeq1$, this establishes free electrons as nearly unbiased, feedback-robust tracers of matter, supporting a key assumption of FRB-based feedback constraints. These properties make electron clustering an ideal application for Hybrid Effective Field Theory (HEFT), which would extend the modeling reach by a further factor of 2–3. The low-$z$ electron spectrum becomes signal-dominated beyond the linear regime within the first few years of next-generation surveys; the models developed here will be necessary on these timescales.

We study the minimal classically scale-invariant extension of the Standard Model, containing a single $\mathbb Z_2$-stabilised real scalar singlet whose mass is not an independent input but is generated dynamically through the quantum effective potential and linked to radiative electroweak symmetry breaking through the Higgs portal. We construct the full two-field one-loop effective potential and introduce a modified on-shell renormalisation scheme that fixes the electroweak vacuum, the Higgs mass, the vanishing Higgs–singlet mixing and the singlet curvature at the physical point, yielding predictions stable under renormalisation-scale variation. Since thermal freeze-out is excluded by direct-detection limits, we identify a highly predictive freeze-in realisation of this model. Imposing perturbativity, vacuum stability and the observed relic abundance leads to a freeze-in solution with a dark-matter mass of around $2~{\rm MeV}$. The predicted electron-scattering cross section for this solution lies far below the current sensitivity of DAMIC-M. Current direct-detection experiments therefore do not constrain this scenario. The minimal scale-invariant singlet model thus provides a robust and highly predictive framework connecting radiative electroweak symmetry breaking and freeze-in dark-matter genesis.

We develop fCAMIRA (forced-mode CAMIRA), a tool for optical cluster confirmation, and apply it to a sample of 129 weak-lensing (WL) shear-selected galaxy clusters identified in aperture-mass maps obtained from the Hyper Suprime-Cam Subaru Strategic Program Three-Year (HSC-SSP Y3) weak-lensing data. fCAMIRA is built upon the CAMIRA cluster-finding algorithm and relies on a red-sequence (RS) galaxy model that is calibrated in a data-driven way. The RS model adopts the metallicity-luminosity relation measured in this work using X-ray-selected clusters up to redshift $z\approx1.3$, followed by the calibration of color offsets using large spectroscopic samples. With the RS model, we build two types of galaxy richness maps, one obtained with a spatial filter matched to a typical cluster size of $R=0.8\,h^{-1}\,\mathrm{Mpc}$ and the other obtained with a fixed angular-size filter identical to that used in constructing the WL aperture-mass maps. The fCAMIRA algorithm utilizes these two richness maps, identifies all optical counterpart candidates along the line of sight of each shear-selected cluster, and measures the cluster photometric redshift from the highest-ranked counterpart. The ranking is determined by the fractional lensing contribution $f_{\mathrm{lens}}$ of each candidate. Using available spectroscopic cluster redshifts, we quantify the mean bias and scatter in the cluster photometric redshifts at levels of approximately 0.005 and 0.008, respectively, demonstrating excellent photometric-redshift performance. We compare the fCAMIRA photometric redshifts with estimates from direct positional cross-matching and find that approximately 8% of the total sample exhibits redshift discrepancies greater than 0.15. This outlier fraction is primarily attributed to projection effects, leading to the misidentification of the optical counterparts. (abridged)

Gravitational waves originating from unresolved sources, or stochastic gravitational wave backgrounds (SGWBs), carry precious information about the physics underlying their diverse sources, and represent an important target for future experimental searches. Here, we model the SGWB due to core collapse supernovae. Using an extensive collection of state-of-the-art, three-dimensional, multi-second supernova simulations, we characterize the two main components of this background: one at $f \sim 10^{-2} - 1$ Hz, due to anisotropic neutrino emission (the gravitational wave memory); the other at $f \gtrsim 10^{2}$ Hz, from the near-core matter dynamics. We find that the memory component offers the best prospects of detection, as its characteristic peak at $f\sim 0.1$ Hz is within the reach of future space-born detectors. At the peak, the energy density might be comparable to that of backgrounds from slow roll inflation and from possible cosmological relics, thus potentially impacting searches of these important signals.

We compute the energy and angular momentum radiated by compact binaries through gravitational waves, scalar and vector fields, using the field-theoretic approach. We observe that each emitted graviton carries away an angular momentum of $\dot{J}/\dot {N}\simeq 2\hbar$, not only for elliptical orbits (as has already been pointed out by Page [arXiv:2409.00305]) but also for hyperbolic orbits. We extend the analysis to scalar and vector radiation and find that the angular momentum carried by each emitted quantum of radiation is approximately $\hbar$, i.e., $\dot{J}/\dot{N}\simeq \hbar$ for both types of radiation. We demonstrate that the radiated energy and angular momentum of binary systems can be determined from the evolution of the orbital angular frequency and the eccentricity of the quasi-Kepler orbit. However, a gauge-invariant decomposition of the angular momentum carried by the emitted particles into the spin and orbital components is not possible.

We revisit aspects of particle production during cosmological first-order phase transitions, and show that the expansion of true-vacuum bubbles can be a copious source of particle production. Specifically, for a massive spin-0 field linearly coupled to the bubble profile, spherical bubbles expanding even at constant radial velocity efficiently produce particles until the local-rest-frame radius of curvature of the bubble wall exceeds the particle Compton wavelength – contrary to the expectation that walls moving at constant speed cannot radiate. We compute the momentum spectrum of the produced particles for both accelerated and constant-velocity expansion histories, identify the regimes of coherent and incoherent production, quantify the validity of the perturbative treatment, and show that this mechanism can parametrically dominate other production processes of feebly coupled particles, including freeze-in. Our results apply to many models of light, feebly coupled particles studied in the literature, providing new sources of dark radiation and dark matter. In particular, a first-order deconfinement-confinement transition in a hidden Yang-Mills sector can abundantly produce axion-like particles, leading to dark radiation and/or dark matter signatures over large regions of parameter space.

We study effective reheating in a scalar–Gauss–Bonnet inflationary model with a phase-space-dependent coupling $\mu(\phi,X)$, in which a compact field-space feature is combined with a bounded kinetic gate. The modified inflationary background determines the pivot-scale quantities and the effective energy density at the end of inflation. These quantities are then used to derive the reheating duration $N_{\rm re}$ and temperature $T_{\rm re}$ through the thermal-history matching relation. We first perform two fixed-pivot reference scans by varying the overall Gauss–Bonnet strength $\lambda_{\rm GB}$ and the kinetic parameter $g_{_X}$ separately. For the reference parameter choices, increasing either parameter increases $N_{\rm re}$ and decreases $T_{\rm re}$ for the selected reheating equations of state. Additional benchmark calculations clarify how these variations depend on the dynamical regime of the model. In the $\lambda_{\rm GB}$ scan, the increase in $N_{\rm re}$ and the decrease in $T_{\rm re}$ persist, although both variations become strongly suppressed when the coupling is more localized or when the end of inflation is controlled more strongly by the E-model potential. In the $g_{_X}$ scan, stronger field-space localization and kinetic saturation can instead lead to a slight decrease in $N_{\rm re}$ and an increase in $T_{\rm re}$ as $g_{_X}$ is increased. When the bounded kinetic contribution is considered together with a weaker overall Gauss–Bonnet interaction, the resulting changes in the reheating quantities become nearly negligible. The fixed-pivot predictions of the representative and alternative benchmarks are compared with CMB constraints.These reheating constraints are then discussed for four representative values of the effective equation-of-state parameter, $\overline{w}_{\rm re}=-1/3,0,2/3,$ and $1$.

We study the strong backreaction regime in axion inflation, in which the friction from Abelian gauge fields generated via the pseudo-scalar interaction $\frac{\beta}{4 M_p} \phi F {\tilde F}$ is comparable to Hubble friction. The non-linear dynamics associated with the gauge fields render this regime phenomenologically particularly interesting, but also notoriously difficult to study with perturbative methods. Quantifying the perturbative control through the direct impact on the spectral backreaction using (i) a first-order gradient-expansion formalism including axion gradients, and (ii) a one-loop in-in calculation, we discover an extended mild backreaction regime at large couplings $\beta$, with a relatively large, nearly constant particle production parameter $\xi$. In passing, we point out that CMB non-Gaussianity bounds impose a general upper limit on the product of axion gauge field coupling and Hubble parameter during inflation, $\beta H/M_p < 1.4\cdot 10^{-3}$, independently of the choice of the axion potential.

Ray et al. [1] have claimed that the acceleration of the Hubble expansion rate inferred from Type Ia supernovae in the Pantheon+ catalogue "is equally strong in the CMB dipole ($q_m = + 1.45$) and anti-dipole hemispheres ($+ 1.55$), directly contradicting the SRS26 anisotropy argument". Here SRS26 refers to our analysis [2] which showed that "locally $q_0$ has a strong dipole anisotropy aligned approximately with the bulk flow, and only a small monopole component remains at distances exceeding a few hundred Mpc". We find that these authors confused Equatorial with Galactic coordinates, thereby getting the CMB dipole direction wrong. Hence their conclusion is baseless.

We present a class of cosmological scenarios in which a preferred spatial direction in the matter sector coexists with an exactly homogeneous and isotropic Friedmann–Lemaître–Robertson–Walker (FLRW) geometry. The Cosmological Principle is realized on shell, with the vector-field equations enforcing a vanishing momentum density and suitable interactions eliminating the anisotropic stress. Consequently, the construction yields an entire FLRW branch without tuning the matter-field initial conditions. The preferred direction reemerges in perturbations, producing direction-dependent propagation and mixing among scalar, vector, and tensor modes already at linear order. In particular, the mixing opens a linear channel through which perturbations in the scalar sector can in principle source gravitational waves. Our construction thus reveals a new route by which preferred-direction physics can leave observable cosmological signatures while remaining hidden in the background geometry.

The aim of the present work is two-fold: {\it(i)} Compute the drift of comoving extragalactic sources in axial Bianchi IX universes, in particular in those satisfying Einstein's equations with positive cosmological constant and comoving dust. {\it(ii)} Linearize this drift (and for comparison the previously calculated Lemaître-Hubble diagram) simultaneously in small anisotropy and in small positive curvature. We find that the drift in this linear approximation coincides with the drift in axial Bianchi I universes. This is not true for Lemaître-Hubble diagrams.

In a previous paper, we described a `standard speed-gun' (SSG) distance that uses the speed of light to standardize a ruler under the assumption that the radio variability seen in blazars is causally limited. The apparent size is then measured with Very Long Baseline Interferometry in order to derive the angular diameter distance. A key limitation of this method is that it requires knowledge of the relativistic Doppler factor. Previously, we estimated the distance to the bright radio source, 3C 84 at the center of the Perseus cluster assuming a Doppler factor of \delta ~ 1. In this paper, we aim to describe how a detected counter-jet and approaching jet proper motions can be used to remove the need for knowledge of the Doppler factor when measuring cosmological distances in this way. Under the assumption of a disk (or spherical) geometry and parameterizing the relationship between the physical emitting region and the variability timescale via a causality correction factor (kappa), we estimate a refined angular diameter distance to 3C 84 (z=0.0178) with statistical errors. Assuming kappa=1, we derive distances of D_A,disk = 78.9(-9.8+11.0) Mpc (or D_A,sphere = 71.2(-8.8+9.7) Mpc). Comparing these results to literature benchmarks, we find that the spherical assumption yields a distance consistent with local Type Ia supernovae calibrated to the SH0ES H0, while a disk-like geometry aligns with expectations from a lower H0 cosmology. Ultimately, this demonstrates that utilizing jet and counter-jet kinematics successfully removes the Doppler-factor dependence from the standard speed-gun method, providing a viable independent distance estimate once the geometric structure of the jet is resolved.

We investigate recent reports of a suppression in the growth rate of cosmic structures inferred from analyses of the $[f\sigma_8](z)$ dataset. To address this issue, we explore the hypothesis that the evolution of matter clustering is more accurately described within the framework of scale-dependent modified gravity. We perform a joint analysis of $[f\sigma_8](z)$, cosmic chronometer $H(z)$ measurements, luminosity distance data, and CMB observations using Markov Chain Monte Carlo techniques to constrain the parameters of a scale-dependent cosmological model and investigate its impact on the evolution of $[f\sigma_8](z)$. Our results indicate that the suppression of the growth rate of large-scale structures is more pronounced during the matter-dominated era than in the dark-energy-dominated epoch. We find evidence for scale-dependent growth at a statistical significance of $2.2\, \sigma$. In addition, we constrain the $S_8$ parameter and find it to be consistent with the value inferred from the CMB observations of the Planck Collaboration. Overall, our analysis shows that $k$-dependent growth models provide a viable explanation for the observed clustering of matter without exacerbating the current cosmological tensions.

BAO and CMB data are somewhat discrepant when interpreted in the context of $\Lambda$cdm, discrepancies that show up as a `matter density deficit' and as a `CMB lensing excess'. One possible resolution is an increased optical depth to scattering off of free electrons in the post-recombination universe, $\tau$, a possibility raised by Sailer et al. 2025 and Jhaveri et al. 2025. Since Planck measurements of the low-$\ell$ polarization `reionization bump' already constrain $\tau$ from standard stellar-driven reionization at $z<10$, we investigate additional optical depth sourced by transient or partial reionization at higher redshift from exotic processes. For specificity, we explore the impact of Hawking radiation from a monochromatic spectrum of primordial black holes, retaining the high-$\ell$ $TT/TE/EE$ data that constrain such histories and varying the reionization redshift jointly. We find that the CMB data do not significantly prefer these additional signals: the boost is at most $\Delta\tau \simeq 0.008$, well short of the $\Delta\tau \simeq 0.03$ that would completely eliminate the moderate discrepancy. The matter density deficit and the lensing excess are not significantly eased: we explain why, tracing it to compensation from the reionization redshift and the residual PBH signal at $\ell > 30$.

We revisit metastable dark energy (DE) models described by a radioactive-like decay law. We consider three scenarios: an effective, exponentially decaying DE component; decay of DE into non-baryonic dark matter (DM); and decay of DE into dark radiation (DR). We constrain the metastable DE models using DESI DR2 baryon acoustic oscillation (BAO) data, Type Ia supernovae (SNIa), cosmic microwave background (CMB) observations, and, for the first time, the current available DESI DR1 full-shape (FS) clustering measurements. The BAO+SNIa combinations show a mild preference for positive $\Gamma/H_0$, with a deviation from the $\Lambda$CDM limit at the $\gtrsim 2\sigma$ level. This corresponds to a decaying DE density and an effective quintessence-like behaviour at low redshift. Once CMB information from either Planck or P-ACT is included, however, the constraints become statistically consistent with $\Gamma/H_0=0$. The FS measurements probe the growth sector and help distinguish the interacting DM-DE behaviour of Model 2 from the effective decaying-DE response of Model 1 and the weaker DR-induced response of Model 3. For CMB+FS+DES-Dovekie, Model 1 shows a slight deviation from $\Gamma/H_0=0$ at the $\gtrsim 2\sigma$ level, while Models 2 and 3 remain consistent with the $\Lambda$CDM limit within $2\sigma$. Overall, metastable DE remains phenomenologically viable: current data allow late-time dynamics but do not provide decisive evidence for a nonzero decay rate. These results motivate extending our analysis to the upcoming DESI DR2 FS data to obtain tighter constraints on metastable dynamics.

Cosmological collider signals are usually discussed for isolated massive particles, whose exchange produces characteristic logarithmic oscillations in primordial correlators. In this work, we study how this signal is modified when the exchanged states form a continuous mass spectrum. We first develop a spectral representation for inflationary correlators mediated by a continuum field. In the soft limit, the non-analytic part of the seed function is expressed as a Fourier–Laplace transform of the spectral weight with respect to logarithmic momentum variables. This representation shows that a continuum superposes clock signals with different frequencies and can dephase the sharp oscillation associated with a single massive particle. We then realize this mechanism in an RS2-like inflationary braneworld, where the inflaton is localized on a de Sitter brane in a five-dimensional AdS bulk. The tensor sector contains a localized massless graviton and a continuum of Kaluza-Klein gravitons starting at \(m=3H/2\). We derive the KK wavefunctions and identify the brane spectral weight fixed by the continuum wavefunction on the brane. Applying this weight to the inflaton four-point function, we find a smooth seed function rather than a persistent logarithmic clock oscillation. This behavior follows from the fact that the continuum starts at zero clock frequency, while the KK spectral weight vanishes near threshold. Our results provide a concrete higher-dimensional example in which deviations from the standard cosmological collider signal encode the structure of a continuum high-energy spectrum.

Core-collapse supernovae are among the most promising yet still undetected sources of gravitational waves. A future detection would provide a direct view of the physical processes occurring deep inside a collapsing star. In this work, we investigate how networks of current and future gravitational-wave detectors can constrain the properties of rapidly rotating core-collapse supernovae using their characteristic core-bounce and early post-bounce signals. Using deep-learning techniques, we estimate the peak frequency, rotation rate, and signal amplitude from noisy detector data and compare the performance of different detector-network configurations. We find that detector networks improve both parameter recovery and sky coverage. For current-generation networks, estimation of the peak frequency is possible out to about 30 kpc, while the rotation rate and signal amplitude remain recoverable out to distances exceeding 100 kpc. Third-generation observatories extend these distances by nearly an order of magnitude.