The scintillation time response of liquid argon has a key role in the discrimination of electronic backgrounds in dark matter search experiments. However, its extraordinary rejection power can be affected by various detector effects such as the delayed light emission of TetraPhenyl Butadiene, the most commonly used wavelength shifter, and the electric drift field applied in Time Projection Chambers. In this work, we characterized the TetraPhenyl Butadiene delayed response and the dependence of the pulse shape discrimination on the electric field, exploiting the data acquired with the ARIS, a small-scale single-phase liquid argon detector exposed to monochromatic neutron and gamma sources at the ALTO facility of IJC Lab in Orsay.

We present results from the search for an isotropic gravitational-wave background using Advanced LIGO and Advanced Virgo data from O1 through O4a, the first part of the fourth observing run. This background is the accumulated signal from unresolved sources throughout cosmic history and encodes information about the merger history of compact binaries throughout the Universe, as well as exotic physics and potentially primordial processes from the early cosmos. Our cross-correlation analysis reveals no statistically significant background signal, enabling us to constrain several theoretical scenarios. For compact binary coalescences which approximately follow a 2/3 power-law spectrum, we constrain the fractional energy density to $\Omega_{\rm GW}(25{\rm Hz})\leq 2.0\times 10^{-9}$ (95% cred.), a factor of 1.7 improvement over previous results. Scale-invariant backgrounds are constrained to $\Omega_{\rm GW}(25{\rm Hz})\leq 2.8\times 10^{-9}$, representing a 2.1x sensitivity gain. We also place new limits on gravity theories predicting non-standard polarization modes and confirm that terrestrial magnetic noise sources remain below detection threshold. Combining these spectral limits with population models for GWTC-4, the latest gravitational-wave event catalog, we find our constraints remain above predicted merger backgrounds but are approaching detectability. The joint analysis combining the background limits shown here with the GWTC-4 catalog enables improved inference of the binary black hole merger rate evolution across cosmic time. Employing GWTC-4 inference results and standard modeling choices, we estimate that the total background arising from compact binary coalescences is $\Omega_{\rm CBC}(25{\rm Hz})={0.9^{+1.1}_{-0.5}\times 10^{-9}}$ at 90% confidence, where the largest contribution is due to binary black holes only, $\Omega_{\rm BBH}(25{\rm Hz})=0.8^{+1.1}_{-0.5}\times 10^{-9}$.

The large-scale analysis task of deciphering gravitational wave signals in the LISA data stream will be difficult, requiring a large amount of computational resources and extensive development of computational methods. Its high dimensionality, multiple model types, and complicated noise profile require a global fit to all parameters and input models simultaneously. In this work, we detail our global fit algorithm, called ``Erebor,'' designed to accomplish this challenging task. It is capable of analysing current state-of-the-art datasets and then growing into the future as more pieces of the pipeline are completed and added. We describe our pipeline strategy, the algorithmic setup, and the results from our analysis of the LDC2A Sangria dataset, which contains Massive Black Hole Binaries, compact Galactic Binaries, and a parameterized noise spectrum whose parameters are unknown to the user. The Erebor algorithm includes three unique and very useful contributions: GPU acceleration for enhanced computational efficiency; ensemble MCMC sampling with multiple MCMC walkers per temperature for better mixing and parallelized sample creation; and special online updates to reversible-jump (or trans-dimensional) sampling distributions to ensure sampler mixing and accurate initial estimates for detectable sources in the data. We recover posterior distributions for all 15 (6) of the injected MBHBs in the LDC2A training (hidden) dataset. We catalog $\sim12000$ Galactic Binaries ($\sim8000$ as high confidence detections) for both the training and hidden datasets. All of the sources and their posterior distributions are provided in publicly available catalogs.

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 perform a frequentist analysis using the standard profile likelihood method for clustering measurements from Data Release 1 of the Dark Energy Spectroscopic Instrument (DESI). While Bayesian inferences for Effective Field Theory models of galaxy clustering can be highly sensitive to the choice of priors for extended cosmological models, frequentist inferences are not susceptible to such effects. We compare Bayesian and frequentist constraints for the parameter set $\{\sigma_8, H_0, \Omega_{\rm{m}}, w_0, w_a\}$ when fitting to the full-shape of the power spectrum multipoles, the post-reconstruction Baryon Acoustic Oscillation (BAO) measurements, as well as external datasets from the CMB and type Ia supernovae measurements. Bayesian prior effects are very significant for the $w_0w_a$CDM model; while the $1 \sigma$ frequentist confidence intervals encompass the maximum a posteriori (MAP), the Bayesian credible intervals almost always exclude the maximum likelihood estimate (MLE) and the MAP - indicating strong prior volume projection effects - unless supernovae data are included. We observe limited prior effects for the $\Lambda$CDM model, due to the reduced number of parameters. When DESI full-shape and BAO data are jointly fit, we obtain the following $1\sigma$ frequentist confidence intervals for $\Lambda$CDM ($w_0w_a$CDM): $\sigma_8 = 0.867^{+0.048}_{-0.041} , \ H_0 = 68.91^{+0.80}_{-0.79} \ \rm{km \ s^{-1}Mpc^{-1}} , \ \Omega_{\rm{m}} = 0.3038\pm0.0110$ ($\sigma_8 = 0.793^{+0.069}_{-0.048} , \ H_0 = 64.9^{+4.8}_{-2.8} \ \rm{km \ s^{-1}Mpc^{-1}} , \ \Omega_{\rm{m}} = 0.369^{+0.029}_{-0.059}$ , $w_0 = -0.24^{+0.17}_{-0.64}$ , $w_a = -2.5^{+1.9}_{}$), corresponding to 0.7$\sigma$, 0.3$\sigma$, 0.7$\sigma$ (1.9$\sigma$, 3.4$\sigma$, 5.6$\sigma$, 5.5$\sigma$, 5.6$\sigma$) shifts between the MLE relative to the Bayesian posterior mean for $\Lambda$CDM ($w_0w_a$CDM) respectively.

A holographic model of QCD in the limit of large number of colors, $N_c$, and massless fermion flavors, $N_f$, but constant ratio $x_f=N_f/N_c$ is analyzed at finite temperature and chemical potential. The five dimensional gravity model contains three bulk fields: a scalar dilaton sourcing ${\rm Tr}F^2$, a scalar tachyon dual to $\bar qq$ and a 4-vector dual to the baryon current $\bar q \gamma^{\mu} q$. The main result is the $\mu,T$ phase diagram of the holographic theory. A first order deconfining transition along $T_h(\mu)$ and a chiral transition at $T_\chi(\mu)>T_h(\mu)$ are found. The chiral transition is of second order for all $\mu$. The dependence of thermodynamical quantities including the speed of sound and susceptibilities on the chemical potential and temperature is computed. A new quantum critical regime is found at zero temperature and finite chemical potential. It is controlled by an AdS$_2\times R^3$ geometry and displays semi-local criticality.

Gravitational lensing allows the detection of binary black holes (BBH) at cosmological distances with chirp masses that appear to be enhanced by $1+z$ in the range $195\%$), from which we infer a factor $\simeq 5$ higher intrinsic rate of NSBH events than BBH events, reflecting a higher proportion of neutron stars formed by early star formation. We predict a distinctive locus for lensed NSBH events in the observed binary mass plane, spanning $1

The Laser Interferometer Space Antenna (LISA) is expected to detect thousands of individually resolved gravitational wave sources, overlapping in time and frequency, on top of unresolved astrophysical and/or primordial backgrounds. Disentangling resolved sources from backgrounds and extracting their parameters in a computationally intensive "global fit" is normally regarded as a necessary step toward reconstructing the properties of the underlying astrophysical populations. Here, we show that it is possible to infer the properties of the most numerous population of LISA sources - Galactic double white dwarfs - directly from the frequency (or, equivalently, time) strain series, by using a simulation-based approach that bypasses the global fit entirely. By training a normalizing flow on a custom-designed compression of simulated LISA frequency series from the Galactic double white dwarf population, we demonstrate how to infer the posterior distribution of population parameters (e.g., mass function, frequency, and spatial distributions). This allows for extracting information on the population parameters from both resolved and unresolved sources simultaneously and in a computationally efficient manner. Our approach to target population properties directly can be readily extended to other source classes (e.g., massive and stellar-mass black holes, extreme mass ratio inspirals), provided fast simulations are available, and to scenarios involving non-Gaussian or non-stationary noise (e.g., data gaps).

Extreme mass ratio inspirals (EMRIs) are among the primary targets for the Laser Interferometer Space Antenna (LISA). The extreme mass ratios of these systems result in relatively weak GW signals, that can be individually resolved only for cosmologically nearby sources (up to $z\approx2$). The incoherent piling up of the signal emitted by unresolved EMRIs generate a confusion noise, that can be formally treated as a stochastic GW background (GWB). In this paper, we estimate the level of this background considering a collection of astrophysically motivated EMRI models, spanning the range of uncertainties affecting EMRI formation. To this end, we employed the innovative Augmented Analytic Kludge waveforms and used the full LISA response function. For each model, we compute the GWB SNR and the number of resolvable sources. Compared to simplified computations of the EMRI signals from the literature, we find that for a given model the GWB SNR is lower by a factor of $\approx 2$ whereas the number of resolvable sources drops by a factor 3-to-5. Nonetheless, the vast majority of the models result in potentially detectable GWB which can also significantly contribute to the overall LISA noise budget in the 1-10 mHz frequency range.

Blazars dominate the extragalactic $\gamma$-ray sky and show pronounced flares. Using public Fermi-LAT light curves for 732 blazars with secure redshifts, I implement an automated pipeline to identify and characterize $\gamma$-ray bursts from blazars (GRBBLs). Each event is modeled with an exponential rise/decay profile, and spectral variability is quantified via a constant fit. From 679 high-quality GRBBLs, I apply extreme deconvolution for unsupervised classification. The GRBBL population is remarkably homogeneous; the most robust split is in achromatic vs. chromatic events, with significant overlap. Removing spectral information yields a luminosity-driven classification in type-1 and type-2 GRBBLs, although this classification is not identified in all tests. This study establishes GRBBL population studies as a tool to study blazars. As a by-product of this project I identify a correlation between peak luminosity and timescales in GRBBLs.

We compare the performance of the flat-sky approximation and Limber approximation for the clustering analysis of the photometric galaxy catalogue of Euclid. We study a 6 bin configuration representing the first data release (DR1) and a 13 bin configuration representative of the third and final data release (DR3). We find that the Limber approximation is sufficiently accurate for the analysis of the wide bins of DR1. Contrarily, the 13 bins of DR3 cannot be modelled accurately with the Limber approximation. Instead, the flat-sky approximation is accurate to below $5\%$ in recovering the angular power spectra of galaxy number counts in both cases and can be used to simplify the computation of the full power spectrum in harmonic space for the data analysis of DR3.

This paper presents a novel coherent multiband analysis framework for characterizing stellar- and intermediate-mass binary black holes using LISA and next-generation ground-based detectors (ET and CE), leveraging the latest developments in the \texttt{PyCBC} pipeline. Given the population parameters inferred from LVK results and LISA's sensitivity limits at high frequencies, most stellar-mass binary black holes would likely have SNRs below 5 in LISA, but the most state-of-the-art multiband parameter estimation methods, such as those using ET and CE posteriors as priors for LISA, typically struggle to analyze sources with a LISA SNR less than 5. We present a novel coherent multiband parameter estimation method that directly calculates a joint likelihood, which is highly efficient; this efficiency is enabled by multiband marginalization of the extrinsic parameter space, implemented using importance sampling, which can work robustly even when the LISA SNR is as low as 3. Having an SNR of $\sim 3$ allows LISA to contribute nearly double the number of multiband sources. Even if LISA only observes for one year, most of the multiband detector-frame chirp mass's 90\% credible interval (less than $10^{-4} \mathrm{M}_\odot$) is still better than that of the most accurately measured events for ET+2CE network in 7.5 years of observation, by at least one order of magnitude. For the first time, we show efficient multiband Bayesian parameter estimation results on the population scale, which paves the way for large-scale astrophysical tests using multibanding.

We derive the energy fluxes radiated by compact binary systems, including "hairy" black holes, in Einstein-Maxwell-dilaton theories, for circular orbits and at quadrupolar order. This enables to include in their resummed effective-one-body (EOB) dynamics the effect of the radiation reaction force at the origin of their inspiral and merger. We also exhibit typical examples of the resulting tensor and scalar waveforms.

Primordial Black Holes (PBH) could dominate in the early universe and, evaporating before Big bang Nucleosynthesis, can provide new freeze in mechanism of dark matter (DM) production. The proposed scenario is considered for two possible mechanisms of PBH formation and the corresponding continuous PBH mass spectra so that the effect of non-single PBH mass spectrum is taken into account in the results of PBH evaporation, by which PBH dominance in the early universe ends. We specify the conditions under which the proposed scenario can explain production of dark matter in very early Universe.

We construct a class of bottom-up holographic models with physics comparable to the one expected from QCD in the Veneziano limit of large N_f and N_c with fixed x = N_f/N_c. The models capture the holographic dynamics of the dilaton (dual to the YM coupling) and a tachyon (dual to the chiral condensate), and are parametrized by the real parameter x, which can take values within the range 0 <= x < 11/2. We analyze the saddle point solutions, and draw the phase diagram at zero temperature and density. The backreaction of flavor on the glue is fully included. We find the conformal window for x >= x_c, and the QCD-like phase with chiral symmetry breaking at x < x_c, where the critical value x_c lies close to four. We also find Miransky scaling as x -> x_c as well as Efimov-like saddle points. By calculating the holographic beta-functions, we demonstrate the "walking" behavior of the coupling in the region near and below x_c.

We propose a framework for Lorentz-invariant Lagrangian field theories where Ostrogradsky's scalar ghosts could be absent. A key ingredient is the generalized Kronecker delta. The general Lagrangians are reformulated in the language of differential forms. The absence of higher order equations of motion for the scalar modes stems from the basic fact that every exact form is closed. The well-established Lagrangian theories for spin-0, spin-1, p-form, spin-2 fields have natural formulations in this framework. We also propose novel building blocks for Lagrangian field theories. Some of them are novel nonlinear derivative terms for spin-2 fields. It is nontrivial that Ostrogradsky's scalar ghosts are absent in these fully nonlinear theories.

Generic solutions are studied in Einstein-scalar gravity in an ansatz that can interpolate between de Sitter and Anti-de Sitter regimes. The scalar potential is arbitrary. All solutions are determined by their end-points in the scalar field space. All such end-points are classified. This provides a complete classification and characterization of the full space of regular solutions. It is shown that there are no regular (Centaur) solutions that interpolate between an AdS boundary and a dS interior, within our ansatz, when $d>2$. This no-go theorem persists in the presence of multiple scalar fields with a non-trivial field space metric. The Gubser classification of regular solutions is also upgraded to include cases that are not Lorentz invariant and do not contain AdS boundaries.

A methodology to provide the polarization angle requirements for different sets of detectors, at a given frequency of a CMB polarization experiment, is presented. The uncertainties in the polarization angle of each detector set are related to a given bias on the tensor-to-scalar ratio $r$ parameter. The approach is grounded in using a linear combination of the detector sets to obtain the CMB polarization signal. In addition, assuming that the uncertainties on the polarization angle are in the small angle limit (lower than a few degrees), it is possible to derive analytic expressions to establish the requirements. The methodology also accounts for possible correlations among detectors, that may originate from the optics, wafers, etc. The approach is applied to the LiteBIRD space mission. We show that, for the most restrictive case (i.e., full correlation of the polarization angle systematics among detector sets), the requirements on the polarization angle uncertainties are of around 1 arcmin at the most sensitive frequency bands (i.e., $\approx 150$ GHz) and of few tens of arcmin at the lowest (i.e., $\approx 40$ GHz) and highest (i.e., $\approx 400$ GHz) observational bands. Conversely, for the least restrictive case (i.e., no correlation of the polarization angle systematics among detector sets), the requirements are $\approx 5$ times less restrictive than for the previous scenario. At the global and the telescope levels, polarization angle knowledge of a few arcmins is sufficient for correlated global systematic errors and can be relaxed by a factor of two for fully uncorrelated errors in detector polarization angle. The reported uncertainty levels are needed in order to have the bias on $r$ due to systematics below the limit established by the LiteBIRD collaboration.

Liquid argon is a bright scintillator with potent particle identification properties, making it an attractive target for direct-detection dark matter searches. The DarkSide-50 dark matter search here reports the first WIMP search results obtained using a target of low-radioactivity argon. DarkSide-50 is a dark matter detector, using two-phase liquid argon time projection chamber, located at the Laboratori Nazionali del Gran Sasso. The underground argon is shown to contain Ar-39 at a level reduced by a factor (1.4 +- 0.2) x 10^3 relative to atmospheric argon. We report a background-free null result from (2616 +- 43) kg d of data, accumulated over 70.9 live-days. When combined with our previous search using an atmospheric argon, the 90 % C.L. upper limit on the WIMP-nucleon spin-independent cross section based on zero events found in the WIMP search regions, is 2.0 x 10^-44 cm^2 (8.6 x 10^-44 cm^2, 8.0 x 10^-43 cm^2) for a WIMP mass of 100 GeV/c^2 (1 TeV/c^2 , 10 TeV/c^2).

We use the effective field theory of dark energy to explore the space of modified gravity models which are capable of driving the present cosmic acceleration. We identify five universal functions of cosmic time that are enough to describe a wide range of theories containing a single scalar degree of freedom in addition to the metric. The first function (the effective equation of state) uniquely controls the expansion history of the universe. The remaining four functions appear in the linear cosmological perturbation equations, but only three of them regulate the growth history of large scale structures. We propose a specific parameterization of such functions in terms of characteristic coefficients that serve as coordinates in the space of modified gravity theories and can be effectively constrained by the next generation of cosmological experiments. We address in full generality the problem of the soundness of the theory against ghost-like and gradient instabilities and show how the space of non-pathological models shrinks when a more negative equation of state parameter is considered. This analysis allows us to locate a large class of stable theories that violate the null energy condition (i.e. super-acceleration models) and to recover, as particular subsets, various models considered so far. Finally, under the assumption that the true underlying cosmological model is the $\Lambda$ Cold Dark Matter ($\Lambda$CDM) scenario, and relying on the figure of merit of EUCLID-like observations, we demonstrate that the theoretical requirement of stability significantly narrows the empirical likelihood, increasing the discriminatory power of data. We also find that the vast majority of these non-pathological theories generating the same expansion history as the $\Lambda$CDM model predict a different, lower, growth rate of cosmic structures.