We provide the most complete analysis so far of quasinormal modes of rotating black holes in a general higher-derivative extension of Einstein's theory. By finding the corrections to the Teukolsky equation and expressing them in a simple form, we are able to apply a generalized continued fraction method that allows us to find the quasinormal mode frequencies including overtones. We obtain the leading-order corrections to the Kerr quasinormal mode frequencies of all the $(l,m,n)$ modes with $l=2,3,4$, $-l\le m\le l$ and $n=0,1,2$, and express them as a function of the black hole spin $\chi$ using polynomial fits. We estimate that our results remain accurate up to spins between $\chi\sim 0.7$ and $\chi\sim 0.95$, depending on the mode. We report that overtones are overall more sensitive to corrections, which is expected from recent literature on this topic. We also discuss the limit of validity of the linear corrections to the quasinormal mode frequencies by estimating the size of nonlinear effects in the higher-derivative couplings. All our results are publicly available in an online repository.

The Euclid mission aims to measure the positions, shapes, and redshifts of over a billion galaxies to provide unprecedented constraints on the nature of dark matter and dark energy. Achieving this goal requires a continuous reassessment of the mission's scientific performance, particularly in terms of its ability to constrain cosmological parameters, as our understanding of how to model large-scale structure observables improves. In this study, we present the first scientific forecasts using CLOE (Cosmology Likelihood for Observables in Euclid), a dedicated Euclid cosmological pipeline developed to support this endeavour. Using advanced Bayesian inference techniques applied to synthetic Euclid-like data, we sample the posterior distribution of cosmological and nuisance parameters across a variety of cosmological models and Euclid primary probes: cosmic shear, angular photometric galaxy clustering, galaxy-galaxy lensing, and spectroscopic galaxy clustering. We validate the capability of CLOE to produce reliable cosmological forecasts, showcasing Euclid's potential to achieve a figure of merit for the dark energy parameters $w_0$ and $w_a$ exceeding 400 when combining all primary probes. Furthermore, we illustrate the behaviour of the posterior probability distribution of the parameters of interest given different priors and scale cuts. Finally, we emphasise the importance of addressing computational challenges, proposing further exploration of innovative data science techniques to efficiently navigate the Euclid high-dimensional parameter space in upcoming cosmological data releases.

We study the projected clustering of photometric luminous red galaxies from the DESI Legacy Survey, combining their angular power spectrum, bispectrum, and cross-correlation with maps of the CMB lensing convergence from the Planck satellite. We employ a perturbative bias expansion in Eulerian space to describe the clustering of galaxies, modelling the power spectrum and bispectrum at one-loop and tree level, respectively. This allows us to use the power spectrum to self-consistently calibrate the perturbative bias parameters. We validate this model against an $N$-body simulation, and show that it can be used up to scales of at least $k_{\rm max}^P\simeq 0.2\,h{\rm Mpc}^{-1}$ and $k_{\rm max}^B\simeq 0.08\,h{\rm Mpc}^{-1}$, saturating the information recovered from the data. We obtain constraints on the amplitude of matter fluctuations $\sigma_8=0.761\pm 0.020$ and the non-relativistic matter fraction $\Omega_m=0.307\pm 0.015$, as well as the combination $S_8\equiv\sigma_8\sqrt{\Omega_m/0.3}=0.769 \pm 0.020$. Including the galaxy bispectrum leads to a $10$-$20\%$ improvement on the cosmological constraints, which are also in good agreement with previous analyses of the same data, and in mild tension with Planck at the $\sim2.5\sigma$ level. This tension is largely present in the standard two-point function dataset, and the addition of the bispectrum increases it slightly, marginally shifting $\sigma_8$ downwards and $\Omega_m$ upwards. Finally, using the bispectrum allows for a substantially more precise measurement of the bias parameters of this sample, which are in reasonable agreement with existing coevolution relations.

The idea that new physics could take the form of feebly interacting particles (FIPs) - particles with a mass below the electroweak scale, but which may have evaded detection due to their tiny couplings or very long lifetime - has gained a lot of traction in the last decade, and numerous experiments have been proposed to search for such particles. It is important, and now very timely, to consistently compare the potential of these experiments for exploring the parameter space of various well-motivated FIPs. The present paper addresses this pressing issue by presenting an open-source tool to estimate the sensitivity of many experiments - located at Fermilab or the CERN's SPS, LHC, and FCC-hh - to various models of FIPs in a unified way: the Mathematica-based code SensCalc.

Recent strong-field regime tests of gravity are so far in agreement with general relativity. In particular, astrophysical black holes appear all to be consistent with the Kerr spacetime, but the statistical error on current observations allows for small yet detectable deviations from this description. Here we study superradiance of scalar and electromagnetic test fields around the Kerr-like Konoplya--Zhidenko black hole and we observe that for large values of the deformation parameter superradiance is highly suppressed with respect to the Kerr case. Surprisingly, there exists a range of small values of the deformation parameter for which the maximum amplification factor is larger than the Kerr one. We also provide a first result about the superradiant instability of these non-Kerr spacetimes against massive scalar fields.

The detection of primordial $B$-mode polarisation of the Cosmic Microwave Background (CMB) is a major observational goal in modern Cosmology, offering a potential window into inflationary physics through the measurement of the tensor-to-scalar ratio $r$. However, the presence of Galactic foregrounds poses significant challenges, possibly biasing the $r$ estimate. In this study we explore the viability of using Minkowski functionals (MFs) as a robustness test to validate a potential $r$ detection by identifying non-Gaussian features associated with foregrounds contamination. To do so, we simulate sky maps as observed by a LiteBIRD-like CMB experiment, with realistic instrumental and foregrounds modelling. The CMB $B$-mode signal is recovered through blind component separation algorithms, and the obtained (biased) value of $r$ is used to generate Gaussian realisation of CMB signal. Their MFs are then compared with those computed on maps contaminated by foreground residual left by component separation, looking for a detection of non-Gaussianity. Our results demonstrate that, with the experimental configuration considered here, MFs can not be reliably adopted as a robustness test of an eventual $r$ detection, as we find that in the majority of the cases MFs are not able to raise significant warnings about the non-Gaussianity induced by the presence of foreground residuals. In the most realistic and refined scenario we adopted, the test is able to flag non-Gaussianity in $\sim 26\%$ of the simulations, meaning that there is no warning on the biased tensor-to-scalar ratio in $\sim 74\%$ of cases. These results suggest that more advanced statistics than MFs must be considered to look for non-Gaussian signatures of foregrounds, in order to be able to perform reliable null tests in future CMB missions.

We perform an analysis of the full shapes of Lyman-$\alpha$ (Ly$\alpha$) forest correlation functions measured from the first data release (DR1) of the Dark Energy Spectroscopic Instrument (DESI). Our analysis focuses on measuring the Alcock-Paczynski (AP) effect and the cosmic growth rate times the amplitude of matter fluctuations in spheres of $8$ $h^{-1}\text{Mpc}$, $f\sigma_8$. We validate our measurements using two different sets of mocks, a series of data splits, and a large set of analysis variations, which were first performed blinded. Our analysis constrains the ratio $D_M/D_H(z_\mathrm{eff})=4.525\pm0.071$, where $D_H=c/H(z)$ is the Hubble distance, $D_M$ is the transverse comoving distance, and the effective redshift is $z_\mathrm{eff}=2.33$. This is a factor of $2.4$ tighter than the Baryon Acoustic Oscillation (BAO) constraint from the same data. When combining with Ly$\alpha$ BAO constraints from DESI DR2, we obtain the ratios $D_H(z_\mathrm{eff})/r_d=8.646\pm0.077$ and $D_M(z_\mathrm{eff})/r_d=38.90\pm0.38$, where $r_d$ is the sound horizon at the drag epoch. We also measure $f\sigma_8(z_\mathrm{eff}) = 0.37\; ^{+0.055}_{-0.065} \,(\mathrm{stat})\, \pm 0.033 \,(\mathrm{sys})$, but we do not use it for cosmological inference due to difficulties in its validation with mocks. In $\Lambda$CDM, our measurements are consistent with both cosmic microwave background (CMB) and galaxy clustering constraints. Using a nucleosynthesis prior but no CMB anisotropy information, we measure the Hubble constant to be $H_0 = 68.3\pm 1.6\;\,{\rm km\,s^{-1}\,Mpc^{-1}}$ within $\Lambda$CDM. Finally, we show that Ly$\alpha$ forest AP measurements can help improve constraints on the dark energy equation of state, and are expected to play an important role in upcoming DESI analyses.

We constrain cosmological parameters from a joint cosmic shear analysis of peak-counts and the two-point shear correlation functions, as measured from the Dark Energy Survey (DES-Y1). We find the structure growth parameter $S_8\equiv \sigma_8\sqrt{\Omega_{\rm m}/0.3} = 0.766^{+0.033}_{-0.038}$, which at 4.8% precision, provides one of the tightest constraints on $S_8$ from the DES-Y1 weak lensing data. In our simulation-based method we determine the expected DES-Y1 peak-count signal for a range of cosmologies sampled in four $w$CDM parameters ($\Omega_{\rm m}$, $\sigma_8$, $h$, $w_0$). We also determine the joint covariance matrix with over 1000 realisations at our fiducial cosmology. With mock DES-Y1 data we calibrate the impact of photometric redshift and shear calibration uncertainty on the peak-count, marginalising over these uncertainties in our cosmological analysis. Using dedicated training samples we show that our measurements are unaffected by mass resolution limits in the simulation, and that our constraints are robust against uncertainty in the effect of baryon feedback. Accurate modelling for the impact of intrinsic alignments on the tomographic peak-count remains a challenge, currently limiting our exploitation of cross-correlated peak counts between high and low redshift bins. We demonstrate that once calibrated, a fully tomographic joint peak-count and correlation functions analysis has the potential to reach a 3% precision on $S_8$ for DES-Y1. Our methodology can be adopted to model any statistic that is sensitive to the non-Gaussian information encoded in the shear field. In order to accelerate the development of these beyond-two-point cosmic shear studies, our simulations are made available to the community upon request.

[Abridged] Galaxy clusters are the most massive gravitationally-bound systems in the universe and are widely considered to be an effective cosmological probe. We propose the first Machine Learning method using galaxy cluster properties to derive unbiased constraints on a set of cosmological parameters, including Omega_m, sigma_8, Omega_b, and h_0. We train the machine learning model with mock catalogs including "measured" quantities from Magneticum multi-cosmology hydrodynamical simulations, like gas mass, gas bolometric luminosity, gas temperature, stellar mass, cluster radius, total mass, velocity dispersion, and redshift, and correctly predict all parameters with uncertainties of the order of ~14% for Omega_m, ~8% for sigma_8, ~6% for Omega_b, and ~3% for h_0. This first test is exceptionally promising, as it shows that machine learning can efficiently map the correlations in the multi-dimensional space of the observed quantities to the cosmological parameter space and narrow down the probability that a given sample belongs to a given cosmological parameter combination. In the future, these ML tools can be applied to cluster samples with multi-wavelength observations from surveys like LSST, CSST, Euclid, Roman in optical and near-infrared bands, and eROSITA in X-rays, to constrain both the cosmology and the effect of the baryonic feedback.

We present a cosmological analysis using the second and third moments of the weak lensing mass (convergence) maps from the first three years of data (Y3) data of the Dark Energy Survey (DES). The survey spans an effective area of 4139 square degrees and uses the images of over 100 million galaxies to reconstruct the convergence field. The second moment of the convergence as a function of smoothing scale contains information similar to standard shear 2-point statistics. The third moment, or the skewness, contains additional non-Gaussian information. The data is analysed in the context of the $\Lambda$CDM model, varying 5 cosmological parameters and 19 nuisance parameters modelling astrophysical and measurement systematics. Our modelling of the observables is completely analytical, and has been tested with simulations in our previous methodology study. We obtain a 1.7\% measurement of the amplitude of fluctuations parameter $S_8\equiv \sigma_8 (\Omega_m/0.3)^{0.5} = 0.784\pm 0.013$. The measurements are shown to be internally consistent across redshift bins, angular scales, and between second and third moments. In particular, the measured third moment is consistent with the expectation of gravitational clustering under the $\Lambda$CDM model. The addition of the third moment improves the constraints on $S_8$ and $\Omega_{\rm m}$ by $\sim$15\% and $\sim$25\% compared to an analysis that only uses second moments. We compare our results with {\it Planck} constraints from the Cosmic Microwave Background (CMB), finding a $2.2$ \textendash $2.8\sigma$ tension in the full parameter space, depending on the combination of moments considered. The third moment independently is in $2.8\sigma$ tension with {\it Planck}, and thus provides a cross-check on analyses of 2-point correlations.

Characterizing the evolution of the star-forming main sequence (SFMS) at high redshift is crucial to contextualize the observed extreme properties of galaxies in the early Universe. We present an analysis of the SFMS and its scatter in the THESAN-ZOOM simulations, where we find a redshift evolution of the SFMS normalization scaling as $\propto (1+z)^{2.64\pm0.03}$, significantly stronger than is typically inferred from observations. We can reproduce the flatter observed evolution by filtering out weakly star-forming galaxies, implying that current observational fits are biased due to a missing population of lulling galaxies or overestimated star-formation rates. We also explore star-formation variability using the scatter of galaxies around the SFMS ($\sigma_{\mathrm{MS}}$). At the population level, the scatter around the SFMS increases with cosmic time, driven by the increased importance of long-term environmental effects in regulating star formation at later times. To study short-term star-formation variability, or ''burstiness'', we isolate the scatter on timescales shorter than 50 Myr. The short-term scatter is larger at higher redshift, indicating that star formation is indeed more bursty in the early Universe. We identify two starburst modes: (i) externally driven, where rapid large-scale inflows trigger and fuel prolonged, extreme star formation episodes, and (ii) internally driven, where cyclical ejection and re-accretion of the interstellar medium in low-mass galaxies drive bursts, even under relatively steady large-scale inflow. Both modes occur at all redshifts, but the increased burstiness of galaxies at higher redshift is due to the increasing prevalence of the more extreme external mode of star formation.

We introduce an extension of the evolution mapping framework to cosmological models that include massive neutrinos. The original evolution mapping framework exploits a degeneracy in the linear matter power spectrum when expressed in ${\rm Mpc}$ units, which compresses its dependence on cosmological parameters into those that affect its shape and a single extra parameter $\sigma_{12}$, defined as the RMS linear variance in spheres of radius $12 {\rm Mpc}$. We show that by promoting the scalar amplitude of fluctuations, $A_{\rm s}$, to a shape parameter, we can additionally describe the suppression due to massive neutrinos at any redshift to sub-0.01\% accuracy across a wide range of masses and for different numbers of mass eigenstates. This methodology has been integrated into the public COMET package, enhancing its ability to emulate predictions of state-of-the-art perturbative models for galaxy clustering, such as the effective field theory (EFT) model. Additionally, the updated software now accommodates a broader cosmological parameter space for the emulator, enables the simultaneous generation of multiple predictions to reduce computation time, and incorporates analytic marginalisation over nuisance parameters to expedite posterior estimation. Finally, we explore the impact of different infrared resummation techniques on galaxy power spectrum multipoles, demonstrating that any discrepancies can be mitigated by EFT counterterms without impacting the cosmological parameters.

Our peculiar velocity imprints a dipole on galaxy density maps derived from redshift surveys. The dipole gives rise to an oscillatory signal in the multipole moments of the observed power spectrum which we indicate as the finger-of-the-observer (FOTO) effect. Using a suite of large mock catalogues mimicking ongoing and future $\textrm{H}\alpha$- and $\textrm{H}\scriptstyle\mathrm{I}$-selected surveys, we demonstrate that the oscillatory features can be measured with a signal-to-noise ratio of up to 7 (depending on the sky area coverage and provided that observational systematics are kept under control on large scales). We also show that the FOTO effect cannot be erased by correcting the individual galaxy redshifts. On the contrary, by leveraging the power of the redshift corrections, we propose a novel method to determine both the magnitude and the direction of our peculiar velocity. After applying this technique to our mock catalogues, we conclude that it can be used to either test the kinematic interpretation of the temperature dipole in the cosmic microwave background or to extract cosmological information such as the matter density parameter and the equation of state of dark energy.

We study the consistency of Scalar Gauss-Bonnet Gravity, a generalization of General Relativity where black holes can develop non-trivial hair by the action of a coupling $F(\Phi){\cal G}$ between a function of a scalar field and the Gauss-Bonnet invariant of the space-time. When properly normalized, interactions induced by this term are weighted by a cut-off, and take the form of an Effective Field Theory expansion. By invoking the existence of a Lorentz invariant, causal, local, and unitary UV completion of the theory, we derive positivity bounds for $n$-to-$n$ scattering amplitudes including exchange of dynamical gravitons. These constrain the value of all even derivatives of the function $F(\Phi)$, and are highly restrictive. They require some of the scales of the theory to be of Planckian order, and rule out most of the models used in the literature for black hole scalarization.

Combining multiple observational probes is a powerful technique to provide robust and precise constraints on cosmological parameters. In this letter, we present the first joint analysis of cluster abundances and auto/cross correlations of three cosmic tracer fields measured from the first year data of the Dark Energy Survey: galaxy density, weak gravitational lensing shear, and cluster density split by optical richness. From a joint analysis of cluster abundances, three cluster cross-correlations, and auto correlations of galaxy density, we obtain $\Omega_{\rm{m}}=0.305^{+0.055}_{-0.038}$ and $\sigma_8=0.783^{+0.064}_{-0.054}$. This result is consistent with constraints from the DES-Y1 galaxy clustering and weak lensing two-point correlation functions for the flat $\nu\Lambda$CDM model. We thus combine cluster abundances and all two-point correlations from three cosmic tracer fields and find improved constraints on cosmological parameters as well as on the cluster observable--mass scaling relation. This analysis is an important advance in both optical cluster cosmology and multi-probe analyses of upcoming wide imaging surveys.

Searches for variations of fundamental constants require a comprehensive understanding of measurement errors. This paper examines a source of error that is usually overlooked: the impact of continuum placement error. We investigate the problem using a high resolution, high signal to noise spectrum of the white dwarf G191$-$B2B. Narrow photospheric absorption lines allow us to search for new physics in the presence of a gravitational field approximately $10^4$ times that on Earth. Modelling photospheric lines requires knowing the underlying spectral continuum level. We describe the development of a fully automated, objective, and reproducible continuum estimation method. Measurements of the fine structure constant are produced using several continuum models. The results show that continuum placement variations result in small systematic shifts in the centroids of narrow photospheric absorption lines which impact significantly on fine structure constant measurements. This effect should therefore be included in the error budgets of future measurements. Our results suggest that continuum placement variations should be investigated in other contexts, including fine structure constant measurements in stars other than white dwarfs. The analysis presented here is based on NiV absorption lines in the photosphere of G191$-$B2B. Curiously, the inferred measurement of the fine structure constant obtained in this paper using NiV (the least negative of our measurements is $\Delta\alpha/\alpha = -1.462 \pm 1.121 \times 10^{-5}$) is inconsistent with the most recent previous G191$-$B2B photospheric measurement using FeV ($\Delta\alpha/\alpha = 6.36 \pm 0.35_{stat} \pm 1.84_{sys} \times 10^{-5}$). Given both measurements are derived from the same spectrum, we presume (but in this work are unable to check) that this 3.2$\sigma$ difference results from unknown laboratory wavelength systematics.