Two-dimensional dilaton gravity provides a valuable framework to study the dynamics of quantum black holes. These models are often coupled to conformal scalar fields, which capture essential quantum effects such as the trace anomaly, while remaining analytically tractable. From the viewpoint of two-dimensional quantum field theory, unitary theories require a positive central charge. However, theories with a total negative central charge naturally arise from the contribution of the Faddeev-Popov ghosts to the effective action. Recent analyses of the Callan-Giddins-Harvey-Strominger (CGHS) model with a Russo-Susskin-Thorlacius (RST) counterterm have shown that a negative central charge can remove curvature singularities in the backreacted geometry. In this work, we argue that singularity resolution arises from the negative central charge itself, rather than the particular dynamics of a given model. To support this, we present analogous results in spherically reduced Einstein gravity.

This document, written by early career researchers (ECRs) in particle physics, aims to represent the perspectives of the European ECR community and serves as input for the 2025--2026 update of the European Strategy for Particle Physics. With input from a community-wide survey, it highlights key challenges faced by ECRs -- career stability, funding access and long-term research opportunities -- while proposing policy recommendations and targeted initiatives. It underscores the importance of practices fostering diverse, equitable, inclusive and healthy workplaces, as well as of stronger ECR communities, and highlights how effective communication and interdisciplinary collaborations reinforce the societal relevance of particle physics and promote continued support for large-scale and long-term projects. Finally, the future of both collider and beyond-collider experiments is addressed, emphasising the critical role of ECRs in shaping future projects. The ECR contribution is formed of two parts: the ten-page executive summary submitted as input to the European Strategy for Particle Physics Update and, as backup document, this extended white paper providing additional context.

Using lattice field simulations of the Abelian-Higgs model, we characterize the simultaneous emission of (scalar and gauge) particles and gravitational waves (GWs) by local string loops. We use {\it network} loops created in a phase transition, and {\it artificial} loops formed by either crossing straight-boosted or curved-static infinite strings. Loops decay via both particle and GW emission, on time scales $\Delta t_{\rm dec} \propto L^p$, where $L$ is the loop length. For particle production, we find $p \simeq 2$ for artificial loops and $p \simeq 1$ for network loops, whilst for GW emission, we find $p \simeq 1$ for all loops. We find that below a critical length, artificial loops decay primarily through particle production, whilst for larger loops GW emission dominates. However, for network loops, which represent more realistic configurations, particle emission always dominates, as supported by our data with length-to-core ratios up to $L/r_\text{c} \lesssim 6000$. Our results indicate that the GW background from a local string network should be greatly suppressed compared to estimations that ignore particle emission.

We search for dark matter (DM) annihilating subhalos of the Milky Way halo among the Fermi Large Area Telescope (LAT) unassociated sources. We construct, for the first time, a statistical model of the unassociated sources at latitudes above 10 degrees. The latter is built as a combination of both DM annihilation subhalos as well as Galactic and extragalactic astrophysical components. The astrophysical components are constructed based on distributions of associated sources, while the distribution of DM subhalos is derived from Monte Carlo simulations. In this model we take into account the differences in the distributions of associated and unassociated sources including both covariate and prior probability shifts (both being forms of ``dataset shifts''). Previous searches of DM subhalos were based on classify-and-count strategies, while the approach adopted in this work is based on quantification learning, which allows one to determine a well-defined statistical interpretation of the contribution of a population of DM subhalos to the unassociated Fermi-LAT sources. In the $b\bar{b}$ annihilation channel and for a range of DM masses from 10 GeV to 1 TeV, we don't find a significant contribution from DM subhalos and derive a statistical 95% confidence upper limit on the DM annihilation cross section in this channel. While the derived limits are consistent with previous classify-and-count approaches, our generative statistical model opens new avenues for population studies of Fermi-LAT sources and, more generally, for searches of anomalies on top of backgrounds in presence of statistical and systematic uncertainties.

We study in detail the fully inhomogeneous non-linear dynamics of axion inflation, identifying three regimes: weak-, mild-, and strong-backreaction, depending on the duration of inflation. We use lattice techniques that explicitly preserve gauge invariance and shift symmetry, and which we validate against other computational methods of the linear dynamics and of the homogeneous backreaction regime. Notably, we demonstrate that the latter fails to accurately describe the truly local dynamics of strong backreaction. We investigate the convergence of simulations of local backreaction, determining the requirements to achieve an accurate description of the dynamics, and providing useful parametrizations of the delay of the end of inflation. Additionally, we identify key features emerging from a proper local treatment of strong backreaction: the dominance of magnetic energy against the electric counterpart, the excitation of the longitudinal mode, and the generation of a scale-dependent chiral (im)balance. Our results underscore the necessity to accurately capture the local nature of the non-linear dynamics of the system, in order to correctly assess phenomenological predictions, such as e.g. the production of gravitational waves and primordial black holes.

We present the full Hubble diagram of photometrically-classified Type Ia supernovae (SNe Ia) from the Dark Energy Survey supernova program (DES-SN). DES-SN discovered more than 20,000 SN candidates and obtained spectroscopic redshifts of 7,000 host galaxies. Based on the light-curve quality, we select 1635 photometrically-identified SNe Ia with spectroscopic redshift 0.10< z <1.13, which is the largest sample of supernovae from any single survey and increases the number of known z>0.5 supernovae by a factor of five. In a companion paper, we present cosmological results of the DES-SN sample combined with 194 spectroscopically-classified SNe Ia at low redshift as an anchor for cosmological fits. Here we present extensive modeling of this combined sample and validate the entire analysis pipeline used to derive distances. We show that the statistical and systematic uncertainties on cosmological parameters are $\sigma_{\Omega_M,{\rm stat+sys}}^{\Lambda{\rm CDM}}=$0.017 in a flat $\Lambda$CDM model, and $(\sigma_{\Omega_M},\sigma_w)_{\rm stat+sys}^{w{\rm CDM}}=$(0.082, 0.152) in a flat $w$CDM model. Combining the DES SN data with the highly complementary CMB measurements by Planck Collaboration (2020) reduces uncertainties on cosmological parameters by a factor of 4. In all cases, statistical uncertainties dominate over systematics. We show that uncertainties due to photometric classification make up less than 10% of the total systematic uncertainty budget. This result sets the stage for the next generation of SN cosmology surveys such as the Vera C. Rubin Observatory's Legacy Survey of Space and Time.

Assuming a minimal $\Lambda$CDM cosmology with three massive neutrinos, the joint analysis of Planck cosmic microwave background data, DESI baryon acoustic oscillations, and distance moduli measurements of Type Ia supernovae from the Pantheon+ sample sets an upper bound on the total neutrino mass, $\sum m_\nu \lesssim 0.06$-$0.07$ eV, that lies barely above the lower limit from oscillation experiments. These constraints are mainly driven by mild differences in the inferred values of the matter density parameter across different probes that can be alleviated by introducing additional background-level degrees of freedom (e.g., by dynamical dark energy models). However, in this work we explore an alternative possibility. Since both $\Omega_\mathrm{m}$ and massive neutrinos critically influence the growth of cosmic structures, we test whether the neutrino mass tension may originate from the way matter clusters, rather than from a breakdown of the $\Lambda$CDM expansion history. To this end, we introduce the growth index $\gamma$, which characterizes the rate at which matter perturbations grow. Deviations from the standard $\Lambda$CDM value ($\gamma \simeq 0.55$) can capture a broad class of models, including non-minimal dark sector physics and modified gravity. We show that allowing $\gamma$ to vary significantly relaxes the neutrino mass bounds to $\sum m_\nu \lesssim 0.13$-$0.2$ eV, removing any tension with terrestrial constraints without altering the inferred value of $\Omega_\mathrm{m}$. However, this comes at the cost of departing from standard growth predictions: to have $\sum m_\nu \gtrsim 0.06$ eV one needs $\gamma > 0.55$, and we find a consistent preference for $\gamma > 0.55$ at the level of $\sim 2\sigma$. This preference increases to $\sim 2.5$-$3\sigma$ when a physically motivated prior $\sum m_\nu \ge 0.06$ eV from oscillation experiments is imposed.

The effect of meson-exchange currents on charged-current quasielastic neutrino scattering with single-nucleon emission is computed and analyzed within the relativistic Fermi gas model. This contribution arises primarily from the interference between one-body and two-body currents, where the two-body operator excites a 1p1h state in the presence of a second, spectator nucleon. The results obtained show a reduction of the vector, axial and vector-axial transverse response functions and, consequently, a decrease in the total neutrino cross section. In addition to a comparison with the non-relativistic limit, other models are also explored, such as the relativistic mean field model for nuclear matter and the superscaling analysis with relativistic effective mass, both of which yield qualitatively similar results.

While Primordial Black Holes (PBHs) with masses $M_{\rm PBH} \gtrsim 10^{-11} \, M_\odot$ cannot comprise the entirety of dark matter, the existence of even a small population of these objects can have profound astrophysical consequences. A sub-dominant population of PBHs will efficiently accrete dark matter particles before matter-radiation equality, giving rise to high-density dark matter spikes. We consider here the scenario in which dark matter is comprised primarily of Weakly Interacting Massive Particles (WIMPs) with a small sub-dominant contribution coming from PBHs, and revisit the constraints on the annihilation of WIMPs in these spikes using observations of the isotropic gamma-ray background (IGRB) and the Cosmic Microwave Background (CMB), for a range of WIMP masses, annihilation channels, cross sections, and PBH mass functions. We find that the constraints derived using the IGRB have been significantly overestimated (in some cases by many orders of magnitude), and that limits obtained using observations of the CMB are typically stronger than, or comparable to, those coming from the IGRB. Importantly, we show that $\sim \mathcal{O}(M_\odot)$ PBHs can still contribute significantly to the dark matter density for sufficiently low WIMP masses and p-wave annihilation cross sections.

The temperature and polarization of the cosmic microwave background (CMB), as measured today, may offer key insights into the topology of the early universe prior to inflation, for example, by discriminating between flat and warped geometries. In this paper, we focus on a Kaluza-Klein model with an extra spatial dimension that compactifies at the Grand Unified Theory (GUT) epoch, subject to mixed Neumann/Dirichlet boundary conditions at fixed points. As a consequence, a set of infrared cutoffs naturally emerges in both the scalar and tensor spectra, leading to observable consequences in the CMB. We examine in detail the possible signatures of such a topology, particularly in relation to the even-odd parity imbalance already reported by the COBE, WMAP and Planck missions in the temperature angular correlations at large scales. Furthermore, we extend our analysis to the existing Planck E-mode polarization data, and to the high-precision B-mode polarization measurements expected from the forthcoming LiteBIRD mission.

Next-generation neutrinoless double beta decay experiments aim for half-life sensitivities of ~$10^{27}$ yr, requiring suppressing backgrounds to <1 count/tonne/yr. For this, any extra background rejection handle, beyond excellent energy resolution and the use of extremely radiopure materials, is of utmost importance. The NEXT experiment exploits differences in the spatial ionization patterns of double beta decay and single-electron events to discriminate signal from background. While the former display two Bragg peak dense ionization regions at the opposite ends of the track, the latter typically have only one such feature. Thus, comparing the energies at the track extremes provides an additional rejection tool. The unique combination of the topology-based background discrimination and excellent energy resolution (1% FWHM at the Q-value of the decay) is the distinguishing feature of NEXT. Previous studies demonstrated a topological background rejection factor of ~5 when reconstructing electron-positron pairs in the $^{208}$Tl 1.6 MeV double escape peak (with Compton events as background), recorded in the NEXT-White demonstrator at the Laboratorio Subterráneo de Canfranc, with 72% signal efficiency. This was recently improved through the use of a deep convolutional neural network to yield a background rejection factor of ~10 with 65% signal efficiency. Here, we present a new reconstruction method, based on the Richardson-Lucy deconvolution algorithm, which allows reversing the blurring induced by electron diffusion and electroluminescence light production in the NEXT TPC. The new method yields highly refined 3D images of reconstructed events, and, as a result, significantly improves the topological background discrimination. When applied to real-data 1.6 MeV $e^-e^+$ pairs, it leads to a background rejection factor of 27 at 57% signal efficiency.

Distillation is a quark-smearing method for the construction of a broad class of hadron operators useful in lattice QCD computations and defined via a projection operator into a vector space of smooth gauge-covariant fields. A new orthonormal basis for this space is constructed which builds in locality. This basis is useful for the construction of stochastic methods to estimate the correlation functions computed in Monte Carlo calculations relevant for hadronic physics.

On February 13th, 2023, the KM3NeT/ARCA telescope detected a neutrino candidate with an estimated energy in the hundreds of PeVs. In this article, the observation of this ultra-high-energy neutrino is discussed in light of null observations above tens of PeV from the IceCube and Pierre Auger observatories. Performing a joint fit of all experiments under the assumption of an isotropic $E^{-2}$ flux, the best-fit single-flavour flux normalisation is $E^2 \Phi^{\rm 1f}_{\nu + \bar \nu} = 7.5 \times 10^{-10}~{\rm GeV cm^{-2} s^{-1} sr^{-1}}$ in the 90% energy range of the KM3NeT event. Furthermore, the ultra-high-energy data are then fit together with the IceCube measurements at lower energies, either with a single power law or with a broken power law, allowing for the presence of a new component in the spectrum. The joint fit including non-observations by other experiments in the ultra-high-energy region shows a slight preference for a break in the PeV regime if the ``High-Energy Starting Events'' sample is included, and no such preference for the other two IceCube samples investigated. A stronger preference for a break appears if only the KM3NeT data is considered in the ultra-high-energy region, though the flux resulting from such a fit would be inconsistent with null observations from IceCube and Pierre Auger. In all cases, the observed tension between KM3NeT and other datasets is of the order of $2.5\sigma-3\sigma$, and increased statistics are required to resolve this apparent tension and better characterise the neutrino landscape at ultra-high energies.

The annihilation of dark matter (DM) particles accumulated in the Sun could produce a flux of neutrinos, which is potentially detectable with neutrino detectors/telescopes and the DM elastic scattering cross section can be constrained. Although the process of DM capture in astrophysical objects like the Sun is commonly assumed to be due to interactions only with nucleons, there are scenarios in which tree-level DM couplings to quarks are absent, and even if loop-induced interactions with nucleons are allowed, scatterings off electrons could be the dominant capture mechanism. We consider this possibility and study in detail all the ingredients necessary to compute the neutrino production rates from DM annihilations in the Sun (capture, annihilation and evaporation rates) for velocity-independent and isotropic, velocity-dependent and isotropic and momentum-dependent scattering cross sections for DM interactions with electrons and compare them with the results obtained for the case of interactions with nucleons. Moreover, we improve the usual calculations in a number of ways and provide analytical expressions in three appendices. Interestingly, we find that the evaporation mass in the case of interactions with electrons could be below the GeV range, depending on the high-velocity tail of the DM distribution in the Sun, which would open a new mass window for searching for this type of scenarios.

We explore the use of symbolic regression to derive compact analytical expressions for angular observables relevant to electroweak boson production at the Large Hadron Collider (LHC). Focusing on the angular coefficients that govern the decay distributions of $W$ and $Z$ bosons, we investigate whether symbolic models can well approximate these quantities, typically computed via computationally costly numerical procedures, with high fidelity and interpretability. Using the PySR package, we first validate the approach in controlled settings, namely in angular distributions in lepton-lepton collisions in QED and in leading-order Drell-Yan production at the LHC. We then apply symbolic regression to extract closed-form expressions for the angular coefficients $A_i$ as functions of transverse momentum, rapidity, and invariant mass, using next-to-leading order simulations of $pp \to \ell^+\ell^-$ events. Our results demonstrate that symbolic regression can produce accurate and generalisable expressions that match Monte Carlo predictions within uncertainties, while preserving interpretability and providing insight into the kinematic dependence of angular observables.

Neutrino masses could originate in seesaw models testable at colliders, with light mediators and an approximate lepton number symmetry. The minimal model of this type contains two quasi-degenerate Majorana fermions forming a pseudo-Dirac pair. An important question is to what extent future colliders will have sensitivity to the splitting between the Majorana components, since this quantity signals the breaking of lepton number and is connected to the light neutrino masses. We consider the production of these neutral heavy leptons at the ILC, where their displaced decays provide a golden signal: a forward-backward charge asymmetry, which depends crucially on the mass splitting between the two Majorana components. We show that this observable can constrain the mass splitting to values much lower than current bounds from neutrinoless double beta decay and natural loop corrections.

In the path towards a muon collider with center of mass energy of 10 TeV or more, a stage at 3 TeV emerges as an appealing option. Reviewing the physics potential of such muon collider is the main purpose of this document. In order to outline the progression of the physics performances across the stages, a few sensitivity projections for higher energy are also presented. There are many opportunities for probing new physics at a 3 TeV muon collider. Some of them are in common with the extensively documented physics case of the CLIC 3 TeV energy stage, and include measuring the Higgs trilinear coupling and testing the possible composite nature of the Higgs boson and of the top quark at the 20 TeV scale. Other opportunities are unique of a 3 TeV muon collider, and stem from the fact that muons are collided rather than electrons. This is exemplified by studying the potential to explore the microscopic origin of the current $g$-2 and $B$-physics anomalies, which are both related with muons.

We explore whether Neural Networks (NNs) can {\it discover} the presence of symmetries as they learn to perform a task. For this, we train hundreds of NNs on a {\it decoy task} based on well-controlled Physics templates, where no information on symmetry is provided. We use the output from the last hidden layer of all these NNs, projected to fewer dimensions, as the input for a symmetry classification task, and show that information on symmetry had indeed been identified by the original NN without guidance. As an interdisciplinary application of this procedure, we identify the presence and level of symmetry in artistic paintings from different styles such as those of Picasso, Pollock and Van Gogh.

The conversion of gravitational to electromagnetic waves in the presence of background magnetic fields is known as the inverse Gertsenshtein effect, analogous to the Primakoff effect for axions. Rephrasing this conversion as a classical electrodynamics problem in the far-field regime of a magnetized region, we derive the angular distribution of the intensity and polarization of the emitted electromagnetic waves. We discuss the interplay of the internal structure of the magnetic field, the polarization of the gravitational wave and the scattering angle, demonstrating for example that a dipolar field can convert an unpolarized stochastic gravitational wave background into polarized electromagnetic emission, with peak emission intensity along the equator. We moreover outline how to incorporate medium effects in this framework, necessary for a realistic 3D description of gravitational wave to photon conversion in the magnetosphere of neutron stars.