Phenomenology of early universe, cosmic microwave background, cosmological parameters, primordial element abundances, extragalactic distance scale, large-scale structure of the universe.
Looking for a broader view? This category is part of:
The Hubble tension is shaped not only by shifts between early- and late-time parameter estimates, but also by the stiffness of the constraints that define them. In this work, we analyze this geometric structure in the wCDM model by separating the discrepancy into two components: a parameter displacement and a directional Fisher curvature. Within the local Gaussian approximation, the quadratic tension along a given direction factorizes into the squared shift and the combined directional curvature contributed by the datasets. Applying this framework to Planck, DESI DR2, and SH0ES, we show that extending \LambdaCDM to wCDM primarily reshapes the Fisher geometry of the CMB constraint rather than opening a genuinely new route to concordance. Allowing the dark-energy equation-of-state parameter w to vary suppresses the leading Planck Fisher eigenvalue to only \sim 2.7 % of its \LambdaCDM value, while producing only a modest rotation of the dominant acoustic-scale eigenmode. The net effect is a strong softening of the effective acoustic rigidity. At the same time, high-precision late-time data, especially from DESI DR2, inject substantial curvature along the expansion-rate direction. This added stiffness acts as a geometric wall, closing off phantom-like escape routes and sharply limiting tension relief within the extended parameter space. Our results indicate that changes in the inferred H_0 tension under model extension are best understood as a reconfiguration of the constraint manifold rather than as evidence for new physical agreement. The shift-curvature decomposition thus offers a simple, fast, and physically transparent way to diagnose cosmological tensions.
We present constraints on the Hubble constant ($H_0$) derived from the observed dispersion measure (DM) distribution of unlocalized fast radio bursts (FRBs). While localized FRBs with redshift measurements have been used to investigate the Hubble tension, their sample remains limited. Here we demonstrate that unlocalized FRBs -- which are far more numerous -- can independently constrain $H_0$ without requiring redshift information, as cosmic expansion imprints itself on their DM distribution. Analyzing a selected sample of 2124 unlocalized FRBs from the CHIME Catalog II, we obtain $H_0 = 73.8^{+14.0}_{-12.3}~\mathrm{km\,s^{-1}\,Mpc^{-1}}$ at the $1σ$ confidence level, corresponding to an uncertainty of about 18%. Breaking the degeneracy between $H_0$ and the characteristic cutoff energy $E_*$ of the FRB isotropic energy distribution would reduce this uncertainty to 9%. This work constitutes the first $H_0$ measurement derived solely from the DM distribution of unlocalized FRBs, highlighting their potential as a new cosmological probe. Future joint analyses with localized FRBs promise even tighter constraints.
We study the performance of the flat $Λ$CDM model and the dynamical dark energy parameterizations $w_0$CDM and $w_0w_a$CDM, in which the dark energy (DE) equation of state is either constant ($w=w_0$) or redshift-dependent [$w(z)=w_0+w_a z/(1+z)$], without and with a varying CMB lensing consistency parameter $A_L$, using combinations of Planck PR4 CMB data (PR4 and lensing), and a compilation of non-CMB data composed of baryon acoustic oscillation (BAO) data that do not include DESI BAO data, Pantheon+ type Ia supernova observations, Hubble parameter measurements $H(z)$, and growth rate $fσ_8$ data. We also compare results from earlier Planck PR3 data with those obtained using PR4 data in order to assess the stability of cosmological constraints. For the largest data combinations, PR3/PR4+lensing+non-CMB, the cosmological parameters inferred from PR3 and PR4 data are consistent, almost all differing by $1σ$ or less. For the $Λ$CDM$+A_L$ model, we have $A_L=1.087 \pm 0.035$ for PR3 and $A_L=1.053 \pm 0.034$ ($1.6σ$ above unity) for PR4, which indicates that the CMB lensing anomaly is reduced when PR4 data are used. For the $w_0 w_a$CDM parameterization, we find $w_0 = -0.863\pm0.060$ (quintessence-like) and $w_0+w_a=-1.37^{+0.19}_{-0.17}$ (phantom-like), suggesting that the current observations favor dynamical DE over a cosmological constant at about $1.8σ$. For the $w_0w_a$CDM$+A_L$ parameterization, we find $w_0=-0.877\pm 0.060$ and $w_0 + w_a =-1.29_{-0.17}^{+0.20}$, corresponding to a preference for dynamical DE over a cosmological constant of about $1.5σ$ and with $A_L = 1.042 \pm 0.037$ exceeding unity at $1.1σ$. These results indicate that while the PR4 data mildly favor a time-evolving DE, part of this preference may be associated with possible residual excess smoothing present in the Planck PR4 CMB anisotropy spectra (abridged).
Scalar-induced gravitational waves (SIGWs) provide a powerful probe of inflationary dynamics on scales far smaller than those accessible to the cosmic microwave background and large-scale structure. In scenarios with a transient ultra-slow-roll (USR) phase, the curvature power spectrum can be strongly enhanced on small scales, potentially generating an observable stochastic GW background. In this regime, scalar dynamics during inflation can become nonlinear, challenging the validity of standard perturbative predictions. Existing semi-analytical calculations of SIGWs rely on the linear evolution of inflation fluctuations. In this work, we compute SIGWs from USR inflation using lattice simulations. We evolve the inflaton field fully nonlinearly during inflation and extract the curvature perturbation nonperturbatively, then simulate its post-reheating horizon re-entry by evolving the Newtonian potential linearly while retaining the full non-Gaussian structure of the initial conditions for the primordial fluctuations in the tensor source. For moderate non-Gaussianity, the semi-analytical prediction captures the correct order of magnitude of the GW signal but receives important corrections. When inflationary non-Gaussianities are large, it can fail dramatically in both amplitude and spectral shape, independently of the overall size of the tensor power spectrum. Our results show that reliable predictions of SIGWs in such scenarios require nonperturbative control of the inflationary scalar dynamics. The code used for this work is available at https://github.com/caravangelo/inflation-easy.git.
At Stage-III sensitivities, cosmic shear $B$ modes unambiguously indicate systematic contamination and are often used to inform data selection and scale cuts for cosmological inference. We validate $B$ modes for the Ultraviolet Near-Infrared Optical Northern Survey (UNIONS)-3500 (2894 deg$^2$, $n_\mathrm{eff} \approx 5.0$ arcmin$^{-2}$) using three $E$/$B$-separable statistics: pure-mode correlation functions $ξ_\pm^{\mathrm{B}}(θ)$, Complete Orthogonal Sets of $E$/$B$-mode Integrals (COSEBI) $B$-mode amplitudes $B_n$, and harmonic-space power spectra $C_\ell^{BB}$. For each statistic, we compute probability-to-exceed (PTE) values over a two-dimensional grid of scale-cut boundaries; our adopted cuts lie in broad stable regions of acceptable PTE. $B$-mode detections and PTE failures on initial catalog versions led us to investigate galaxy size cuts and stellar halo masking. After cuts, all three statistics pass the null test (minimum PTE $= 0.18$). Before scale cuts, we measure an oscillatory COSEBI $B$-mode pattern consistent with repeating additive shear bias, a detector-level effect seen across multiple Stage-III surveys including CFHTLenS, which used the same MegaCam camera; scale cuts that exclude the charge-coupled device (CCD) angular scale suppress it. Although these statistics probe the same two-point shear field, scale cuts in one do not map exactly onto cuts in another, because their respective filter functions weight angular scales differently. The most conservative validation therefore requires scale and sample selections that pass null tests across all frameworks simultaneously, an approach that applies directly to Stage-IV surveys where systematic errors dominate.
The redshift drift provides a kinematic test of the cosmic expansion history through the slow time variation of the redshift of comoving sources. Motivated by the expected Sandage-Loeb measurements from future facilities, we investigate the drift within a cosmographic framework, modeling the Hubble rate through both a second-order Taylor expansion and a $(2,1)$ Padé approximant. We constrain the cosmographic parameters $(H_0,q_0,j_0)$ by combining Pantheon+ and SH0ES type Ia supernovae with gamma-ray bursts and then examine the impact of adding baryon acoustic oscillation measurements from the second DESI data release. The resulting constraints are used to construct a mock Sandage-Loeb catalog, after which the analyses are repeated including the simulated drift data. In this way, we assess the internal consistency of the reconstructed background rather than perform an independent forecast. Accordingly, we find that, for the SNeIa+GRB analysis, the Taylor reconstruction is compatible at the $1σ$ level with the $ω_0ω_1$CDM scenario, whereas the Padé parameterization improves the agreement of $q_0$ with the $Λ$CDM paradigm. Once DESI BAO data are included, the agreement with the reference background models weakens to the $2σ$ level. The addition of the mock Sandage-Loeb sample mainly tightens the bounds on $q_0$ and $j_0$, with moderate shifts in the central values. We finally compare the reconstructed redshift drift with the corresponding behavior predicted by the $Λ$CDM and $ω_0ω_1$CDM scenarios.
The potential association between gravitational waves (GWs) and fast radio bursts (FRBs) offers a unique multi-messenger probe for cosmology. In this paper, we develop a redshift-independent framework to constrain cosmological parameters using the luminosity distance - dispersion measure relation, accounting for realistic astrophysical uncertainties. We perform a comprehensive comparative analysis across different GWs detector sensitivities and modeling assumptions. Specifically, we investigate the performance of the current LIGO-Virgo (LV) network (at $z < 0.2$) versus the future Cosmic Explorer (CE). Our study further evaluates the impact of different dispersion measure (DM) distributions -- specifically the corrected Macquart's PDF (Zhuge+2025) and the log-normal distribution -- and explores the influence of including or excluding host galaxy DM contributions. Using realistic simulated observations, we find that while the current LV network lacks the precision to provide meaningful constraints, CE will enable high-precision cosmology. Even without spectroscopic redshifts, CE observations can effectively break parameter degeneracies and robustly constrain both cosmology and host galaxy parameters. These results highlight the necessity of next-generation detectors.
Dark-matter haloes do not form in isolation but within the surrounding cosmic web. By the time a halo begins to collapse, its larger-scale environment has typically collapsed along two axes, forming filaments that channel anisotropic infall toward the halo. In this work, we derive from first principles the characteristic Lagrangian scale ratio at which such an anisotropic tidal field most strongly influences halo formation. Specifically, we identify the inflection point of the conditional probability that the tidal field, smoothed on a scale Rsd, undergoes two-dimensional compression, given the presence of a density peak of rarity nu on a smaller scale Rpk. For a standard LambdaCDM cosmology, we find (Rsd/Rpk)infl = 2.2 + (nu-2.5) for Rpk corresponding to a tophat filter of 8Mpc/h. This result implies that the anisotropic tidal influence on a collapsing halo typically extends to 2-3 times the size of its Lagrangian patch. Recast as a function of formation redshift z, the characteristic filament scale around 2.5 sigma peaks can be approximated by Rsd(z) = 31 /(2+(1+z)**2)Mpc/h. We provide practical scaling laws for selecting dynamically relevant smoothing scales in large-scale surveys and for setting initial patch sizes in high-resolution zoom simulations.
The early dark energy resolution of Hubble tension seems to be suggesting a scale-invariant Harrison-Zeldovich spectrum of primordial scalar perturbation, i.e. $n_s=1$ ($|n_s-1|\sim {\cal O}(0.001)$) for $H_0\sim 73$km/s/Mpc. In this work, we propose a possibility to acquire $n_s=1$ in single field slow-roll models of inflation. In our consideration, the potential of inflaton during inflation still preserve the shape of well-known single field inflation models in deep slow-roll region, but inflation ends suddenly due to a large step of inflaton potential. In particular, we investigate the implication of our scheme for chaotic inflation and Starobinski inflation, and show how they can be compatible with the observation for $n_s=1$.
Investigating the properties of the first stars in the universe is essential, yet it remains an open question. One way to explore these stars is by examining their effects on the surrounding gas during the epoch of reionization. In this study, we investigate whether the 21-cm global signal can constrain the typical mass and star formation efficiency of first-generation stars. We perform semi-numerical simulations that include the escape fraction of ionizing photons, which depends on stellar and halo masses, as well as the heating structure surrounding a halo that hosts the first star, determined by radiation hydrodynamics (RHD) simulations. By applying Fisher analysis, while accounting for foreground emissions, we demonstrate that future observations with instruments such as the Radio Experiment for the Analysis of Cosmic Hydrogen (REACH) could provide meaningful constraints on these properties.
We carry out a test of the fundamental Etherington relation (cosmic distance duality relation) which relates the luminosity distance $D_{\rm L}$ and angular diameter distance $D_{\rm A}$ in metric theories of gravity. We use the latest measurements of the angular diameter distance as a function of redshift from the Dark Energy Spectroscopic Instrument Data Release 2 (DESI-DR2) and the luminosity distance from a variety of compilations of Supernovae of Type Ia (SNeIa). Our results indicate that these measurements are statistically consistent with the Etherington relation. In addition to providing a confirmation of the underlying assumptions of the Etherington relation, i.e., the metric nature of gravity, Lorentz invariance and photon number conservation, our results are also a stringent test of any residual systematic effects. We interpret the absence of evidence of any deviation from this relation to constrain the evolution of the absolute magnitude of SNeIa to $dM/dz = 0.07 \pm 0.07$ over and above the systematics that are already accounted for in the SNeIa analyses. We discuss how the Etherington relation can be used to constrain systematic parameters in the analyses of dynamical dark energy using geometric probes, to make it more robust against systematic effects.
Cosmological constraints on neutrino mass offer a promising avenue for advancing our understanding of both fundamental particle physics and the evolution of cosmic large-scale structure. To overcome challenges associated with traditional probes of neutrino mass, particularly degeneracies with other parameters, we consider topological summaries of the cosmic web based on the formalism of persistent homology. We introduce persistence strips, a novel representation that segments Betti curves by topological persistence, effectively condensing two-dimensional persistence diagrams into a set of one-dimensional curves. Applied to the FLAMINGO suite of cosmological simulations, these topological descriptors demonstrate pronounced sensitivity to neutrino mass. By constructing an emulator spanning a 10-dimensional $w_0 w_a\text{CDM} +ν$ cosmological parameter space that includes parameters degenerate with neutrino masses in conventional approaches, assuming a volume of only $(350 \, \mathrm{Mpc})^3$, we obtain neutrino mass constraints with an uncertainty of $0.05 \, \mathrm{eV}$ for the total matter field and $0.13 \, \mathrm{eV}$ for the dark matter-only field, with the strongest constraints coming from void topology. Persistence strips exhibit roughly twice the constraining power of unbinned Betti curves and, through their multi-scale, environment-dependent sensitivity, systematically break degeneracies between neutrino mass and other cosmological parameters. We pinpoint the precise physical origin of the signal, identifying two equally important contributions: sensitivity to the neutrino mass fraction, which is highest in underdense regions, and the impact of neutrinos on the distribution of dark matter. Our findings highlight the particular promise of applying topological statistics to weak lensing measurements, which directly probe the total matter distribution.
We discuss the amplitude and phase fluctuations of gravitational waves due to wave optics lensing in the presence of both a strong lens and cosmological weak lenses. By applying the geometric optics approximation to the strong lens and treating the weak lensing potential perturbatively, we obtain the amplification factor up to the second order in the weak lensing potential. Additionally, we establish a methodology to systematically evaluate the weak lensing effects based on diagrammatic rules. Based on the derived amplification factor, we evaluate the statistics of the fluctuations and demonstrate that the consistency relations originally established in the absence of a strong lens still hold in exactly the same form when a strong lens is present. The physical origin of these relations is also discussed. Furthermore, we demonstrate that for the mean of the weak lensing signal, both the magnification of the signal and the shift of the Fresnel scale to larger scales occur, consistent with the behavior observed in the variance.
We propose a new non-minimally coupled quintessence model to account for the late-time dark energy dynamics indicated by recent DESI measurements. Within this framework, the quintessence density begins to decrease only when it starts to dominate the universe, which naturally accounts for the late-time onset of dark energy weakening. The coupling also induces a sign change in the effective energy transfer between dark matter and dark energy during cosmic evolution. While the scalar field itself remains canonical and never crosses the phantom divide, the modified evolution of the dark matter density gives rise to an effective crossing behavior in the observationally inferred dark energy sector. Compared with both $Λ$CDM and $w_0w_a$CDM models, our model is favored more strongly by current cosmological data. This work may provide a promising avenue for understanding the observational late-time weakening of dark energy and the origin of its dynamics.
Recently, the Dark Energy Spectroscopic Instrument (DESI) collaboration has presented results indicating that dark energy may exhibit dynamical behavior. Calibrated gamma-ray burst (GRB) correlations can be employed to verify or reject a time-evolution of the dark energy (DE) equation of state, $ω(z)$, up to redshifts $z\sim 9$. We use the most updated catalog of GRBs fulfilling the Combo correlation and improve its calibration employing three catalogs of type Ia supernovae at redshifts $z\leq0.075$ and the Bézier interpolation of the Hubble rate, as an alternative to the cosmographic series that fails to be constraining at high redshifts. To test the evolution of $ω(z)$, we adopt a model-independent, redshift-binned DE parametrization. In both the calibration and the DE reconstruction analyses the impact of the spatial curvature on the results is explored. The calibrated Combo correlation yields a Hubble constant $H_0\sim70$ km/s/Mpc which alleviates the existing Hubble tension and is broadly consistent with current measurements, although the uncertainties prevent a high-precision measurement. Regarding the reconstruction of $ω(z)$ of DE, spatially curved scenarios are disfavored and, despite the apparent ''phantom'' behavior at $z\lesssim0.55$ due to the limited statistics caused by the shortage of nearby events, at $z>0.55$ the analysis provides statistically robust evidence in favor of the cosmological constant scenario. The Combo correlation alleviates the Hubble tension and shows no significant evidence in favor of dynamical DE. This suggests that GRBs, as distance indicators, are broadly consistent with the current cosmic distance ladder.
The work is devoted to the study of the possibilities of observational manifestations of traversable wormholes (WHs). The simplest binary system model consisting of a traversable WH candidate (black hole (BH), supermassive BH) and a companion star, whose motion is perturbed by a massive object (star) located on the other side of the wormhole throat, is considered. In the first case of supermassive BH as WH candidate the perturbing acceleration is analyzed and compared with a competing effect in the form of the stochastic influence of stars adjacent to the companion star. In the second case the features of the change in the radial velocity of the companion star in the model of a wide binary system with a WH are also analyzed in order to distinguish it from the following models: 1) a binary system with a BH, and 2) a triple system. For the observational accuracy in radial velocity expected in the near future, at the level of 1.5 km/s the radial velocity perturbations are detectable for all considered observation time spans. For a more realistic accuracy of 10 km/s, the spectral analysis methods become statistically significant after approximately 17 years of data accumulation. The application of spectral and non-parametric methods significantly decreases the required accumulation time compared to matched-filtering applied in isolation.
We investigate the imprints of the Generalized Uncertainty Principle on cosmological scales by using redshift-space distortion measurements in combination with background cosmological data to determine constraints on the deformation parameter $β$. We consider the modified Poisson bracket related to the existence of a minimal length, which leads to a modified Raychaudhuri equation for the standard $Λ$CDM model and gives rise to a phenomenological one-parameter dynamical dark energy scenario. Through this modification, we can reveal the effects of the minimal length on the late-time structure of the universe. We employ the $f$ and $fσ_8$ measurements of the growth rate combined with background data, including cosmic chronometers, baryon acoustic oscillations and Type Ia supernova observations. The observational constraints reveal a systematically negative value for the deformation parameter $β$, with the $Λ$CDM limit lying within the 95\% credible interval. When supernova data are included, the Akaike Information Criterion indicates weak-to-strong support in favour of the GUP-modified model depending on the SNIa catalogue, while the Bayesian evidence suggests a weak preference.
In arXiv:2507.07171, we demonstrate that the high-redshift supermassive black holes in the so-called "little red dots" discovered by James Webb Space Telescope (JWST) can be explained by the primordial black hole (PBH) clustering on small scales. In this paper, we present a comprehensive simulation of the successive PBH mergers within a cluster by solving the Smoluchowski coagulation equation. We derive the coagulation kernel considering both cases with and without the effects of mass segregation. Then we employ the Monte Carlo method to solve the equation, implementing the full-conditioning scheme using the discrete inverse transformation method. Our simulations determine the runaway timescales of clusters and the mass population evolution of PBHs across a wide range of cosmic redshifts, depending on the number of PBHs within the cluster and the associated density.
We present the first public cosmic-web environment catalog built on any DESI data release. Using ASTRA (Algorithm for Stochastic Topological RAnking), we classify each object in the DESI Early Data Release into void, sheet, filament, or knot by combining observed positions with matched random catalogs, without reconstructing a continuous density field. We apply this method to four DESI extragalactic tracers (BGS, LRG, ELG, and QSO) across the 20 EDR rosettes ($\sim 175$ deg$^2$ total), running 100 realizations per tracer-zone pair to derive per-object membership probabilities and classification entropies. We calibrate the classification thresholds using BGS as an anchor to match the volume-filling fractions reported for GAMA, and recover a physically consistent web morphology across all tracers. For BGS, the resulting web-type fractions and the environmental dependence of star formation are consistent with GAMA, COSMOS, and SDSS-based references, validating the method against established benchmarks. A normalized mutual information analysis on BGS reveals a clear dependence of the statistical associations between galaxy color, stellar mass, and specific star formation rate across environments. These results provide a new observational baseline for galaxy evolution studies with DESI. All data products and the open-source pipeline are publicly available.
In this work we account for this skewness in parameter inference by modelling the likelihood through an Edgeworth expansion which involves the complete skewness tensor, composed of 1-point, 2-point, and 3-point correlators. To simplify the calculations of this expansion we perform a change of basis which reduces the precision matrix to the identity. In this basis, the off-diagonal elements of the skewness tensor are consistent with zero, while the amplitude of its diagonal match the level expected for a Gaussian underlying field. We perform parameter inference with this likelihood model and find that including only the diagonal part of the skewness is sufficient, while incorporating the full skewness tensor injects noise without improving accuracy. Despite the estimated excess skewness in the original basis, the cosmological constraints remain effectively unchanged when adopting a Gaussian likelihood or considering the more complete Edgeworth expansion, with variations in the figure of merit of cosmological parameters between the two cases below $5\%$. This result remains unchanged against variations of the survey volume and geometry, scale-cut, and 2-point statistic (power spectrum or correlation function). Using $10\, 000$ cloned \Euclid large mocks based on realistic galaxy catalogues with characteristics close to future \Euclid data, we find no detectable excess skewness on intermediate scales, due to the level of shot noise expected for the \Euclid spectroscopic sample. We conclude that the Gaussian likelihood assumption is robust for \Euclid 2-point statistics analyses in both Fourier and configuration space.