We consider the economical 331 model, based on $\beta=-1/\sqrt{3}$, with three $SU(3)$ triplets with a softly broken $\mathbb{Z}_2$ symmetry. The resulting scalar potential is commonly used in phenomenology. We systematically determine all the potential minima and obtain the conditions under which the electroweak vacuum is global with the help of orbit space methods. For the case the electroweak vacuum is not global, we calculate bounds on the scalar couplings from metastability. We find a parametrisation of the potential couplings in terms of physical quantities and use it to show the available parameter space.

Starobinsky inflation and non-minimally coupled Higgs inflation have been among the most favored models of the early universe, as their predictions for the scalar spectral index $n_s$ and tensor-to-scalar ratio $r$ fall comfortably within the constraints set by Planck and BICEP/Keck. However, new results from the Atacama Cosmology Telescope (ACT) suggest a preference for higher values of $n_s$, introducing tension with the simplest realizations of these models. In this work, being agnostic about the nature of the inflaton, we show that incorporating one-loop corrections to a Higgs-like inflationary scenario leads to a shift in the predicted value of $n_s$, which brings Higgs-like inflation into better agreement with ACT observations. Remarkably, we find that this can be achieved with non-minimal couplings \xi < 1, in contrast to the large values typically required in conventional Higgs inflation, thereby pushing any unitarity-violation scale above the Planck scale. The effect is even more significant when the model is formulated in the Palatini approach, where the modified field-space structure naturally enhances deviations from the metric case. These findings highlight the importance of quantum corrections and gravitational degrees of freedom in refining inflationary predictions in light of new data.

We analyze a model of quintessence governed by an exponential potential and non-minimally coupled to gravity, in light of recent datasets, including cosmic microwave background, baryon acoustic oscillations, and supernovae distance moduli observations. Mainly focusing on the Palatini formulation of gravity, a phase space analysis reveals the existence of a late-time stable de Sitter attractor as long as the non-minimal coupling constant is negative, regardless of the value of the slope of the exponential. Fitting to CMB+DESI+DESY5 data, we find strong evidence for our model over $\Lambda$CDM, with a Bayes factor $\log B = 5.52$. Furthermore, the data seem to prefer dynamical dark energy at $>3\sigma$ C.L. and a phantom crossing in the barotropic parameter of dark energy at $2-3\sigma$ C.L.. We find that the scalar field dynamics in the Palatini formalism provides marginally better agreement to the data compared to the metric formalism.

In this paper we review the physics opportunities at linear $e^+e^-$ colliders with a special focus on high centre-of-mass energies and beam polarisation, take a fresh look at the various accelerator technologies available or under development and, for the first time, discuss how a facility first equipped with a technology mature today could be upgraded with technologies of tomorrow to reach much higher energies and/or luminosities. In addition, we will discuss detectors and alternative collider modes, as well as opportunities for beyond-collider experiments and R\&D facilities as part of a linear collider facility (LCF). The material of this paper will support all plans for $e^+e^-$ linear colliders and additional opportunities they offer, independently of technology choice or proposed site, as well as R\&D for advanced accelerator technologies. This joint perspective on the physics goals, early technologies and upgrade strategies has been developed by the LCVision team based on an initial discussion at LCWS2024 in Tokyo and a follow-up at the LCVision Community Event at CERN in January 2025. It heavily builds on decades of achievements of the global linear collider community, in particular in the context of CLIC and ILC.

The inert doublet model is a well-motivated extension of the Standard Model that contains a dark matter candidate and modifies the dynamics of the electroweak symmetry breaking. In order to detail its phenomenology, we perform a comprehensive study of cosmic phase transitions and gravitational wave signals implied by the framework, accounting for the latest results of collider experiments. We require the neutral inert scalar to constitute, at least, a subdominant part of the observed dark matter abundance. While most of the phase transitions proceed through a single step, we identify regions of the parameter space where the electroweak vacuum is reached after multiple phase transitions. The resulting gravitational wave spectrum is generally dominated by single-step transitions and, in part of the parameter space, falls within the reach of future gravitational wave detectors such as LISA, BBO or DECIGO. We find that direct detection experiments efficiently probe the part of parameter space associated with multi-step phase transitions, which remain unconstrained only in the Higgs resonance region testable with future monojet searches. The implications of the new determination of the $W$ boson mass are also discussed.

We propose a mechanism to generate Primordial Black Holes (PBHs) which is independent of cosmological inflation and occurs slightly below the QCD phase transition. Our setup relies on the collapse of long-lived string-domain wall networks and is naturally realized in QCD axion models with domain wall number $N_{DW}>1$ and Peccei-Quinn symmetry broken after inflation. In our framework, dark matter is mostly composed of axions in the meV mass range along with a small fraction, $\Omega_{\text{PBH}} \gtrsim 10^{-6} \Omega_{\text{CDM}}$ of heavy $M \sim 10^4-10^7 M_\odot$ PBHs. The latter could play a role in alleviating some of the shortcomings of the $\Lambda$CDM model on sub-galactic scales. The scenario has distinct signatures in ongoing axion searches as well as gravitational wave observatories.

We study the formation of primordial black holes (PBHs) in strongly supercooled first-order phase transitions. The mechanism is based on the presence of remnants dominated by the false vacuum that scale slower with the expansion of the Universe than their surroundings where this energy was already converted into radiation. We compute the PBH formation from these remnants including the contribution from the false vacuum and the bubble walls, by estimating the collapse using the hoop conjecture and by considering both regions collapsing immediately when entering the horizon and sub-horizon regions that collapse as their compactness grows. We show that for exponential bubble nucleation rate, $\Gamma \propto e^{\beta t}$, the primordial black hole formation implies $\beta/H \gtrsim 3.8$, where $H$ denotes the Hubble rate, if the potential energy of the false vacuum is $\Delta V \lesssim (10^{12} {\rm GeV})^4$, as otherwise a too large abundance of long-lived PBHs forms. The observed dark matter abundance can be formed in asteroid mass PBHs if $\beta/H \simeq 3.8$ and $10^5 {\rm GeV} \lesssim \Delta V^{1/4} \lesssim 10^8 {\rm GeV}$. Finally, we consider also the effect of the second order correction to the exponential nucleation rate showing that the PBH abundance is mainly determined by the average radius of the true vacuum bubbles.

In light of the latest results from ACT observations we review a class of potentials labeled as fractional attractors, that can originate from Palatini gravity. We show that, for certain choices of the scalar potential $V(\phi)$, the fractional attractors predict both a spectral index $n_s$ and a tensor-to-scalar ratio $r$ that fall within the $1\sigma$ region of the combined ACT+Planck data for a wide range of parameters. We also provide a numerical fit for the parameter space of this models in the case of a simple quadratic and quartic fractional potential.

This note collects together useful unit conversions and numerical values from early universe cosmology. It is a quick reference that can be used to make easy order-of-magnitude estimates. Included are tables for unit conversions, the thermal history of the universe, and collected properties of astronomical objects. The note also introduces a modifiable Mathematica package NaturalUnits (newest version 2.0), which makes it easy to convert between natural and physical units.

In this letter, we provide a simple algorithm, anyPUB, to systematically derive the $2 \rightarrow 2$ scattering matrix in the high-energy limit for any kind of models, irrespective of their gauge group or their field representation. After computing the eigenvalues analytically and/or numerically from this matrix, we impose perturbative unitarity bounds on them. We tested our method on various models and validated the results against the literature. Finally, as a concrete application of our approach, we discuss the case of the minimal left-right symmetric model and derive, for the first time, the perturbative unitarity constraints in the Pati-Salam model.

We propose a threshold model of language competition which includes intermediate bilingual state. The model is based on the Minett-Wang model but through the introduction of thresholds in the language shift rates it incorporates the effects of memory and learning. The model is piecewise-linear, allowing the exact analytical treatment. We study the symmetric case where two competing languages are equivalent in terms of status and social pressure and provide a complete list of the various dynamical regimes. We also study several limiting regimes corresponding to asymmetric systems and characterize the full spectrum of possible asymptotic behaviors. Unlike the Minett-Wang model, which always predicts the extinction of one of the languages, the proposed new model exhibits a wide range of possible equilibrium scenarios, including equilibrium states of coexistence. Most commonly, in such coexistence regimes the minority language speakers are either completely monolingual or completely bilingual.

Efficient and accurate algorithms are necessary to reconstruct particles in the highly granular detectors anticipated at the High-Luminosity Large Hadron Collider and the Future Circular Collider. We study scalable machine learning models for event reconstruction in electron-positron collisions based on a full detector simulation. Particle-flow reconstruction can be formulated as a supervised learning task using tracks and calorimeter clusters. We compare a graph neural network and kernel-based transformer and demonstrate that we can avoid quadratic operations while achieving realistic reconstruction. We show that hyperparameter tuning significantly improves the performance of the models. The best graph neural network model shows improvement in the jet transverse momentum resolution by up to 50% compared to the rule-based algorithm. The resulting model is portable across Nvidia, AMD and Habana hardware. Accurate and fast machine-learning based reconstruction can significantly improve future measurements at colliders.

We revise primordial black holes (PBHs) production in the axion-curvaton model, in light of recent developments in the computation of their abundance accounting for non-gaussianities (NGs) in the curvature perturbation up to all orders. We find that NGs intrinsically generated in such scenarios have a relevant impact on the phenomenology associated to PBHs and, in particular, on the relation between the abundance and the signal of second-order gravitational waves. We show that this model could explain both the totality of dark matter in the asteroid mass range and the tentative signal reported by the NANOGrav and IPTA collaborations in the nano-Hz frequency range. En route, we provide a new, explicit computation of the power spectrum of curvature perturbations going beyond the sudden-decay approximation.

Thermal corrections in classically conformal models typically induce a strong first-order electroweak phase transition, thereby resulting in a stochastic gravitational wave background that could be detectable at gravitational wave observatories. After reviewing the basics of classically conformal scenarios, in this paper we investigate the phase transition dynamics in a thermal environment and the related gravitational wave phenomenology within the framework of scalar conformal extensions of the Standard Model. We find that minimal extensions involving only one additional scalar field struggle to reproduce the correct phase transition dynamics once thermal corrections are accounted for. Next-to-minimal models, instead, yield the desired electroweak symmetry breaking and typically result in a very strong gravitational wave signal.

Dark Matter (DM) can become captured, deposit annihilation energy, and hence increase the heat flow in exoplanets and brown dwarfs. Detecting such a DM-induced heating in a population of exoplanets in the inner kpc of the Milky Way thus provides potential sensitivity to the galactic DM halo parameters. We develop a Bayesian Hierarchical Model to investigate the feasibility of DM discovery with exoplanets and examine future prospects to recover the spatial distribution of DM in the Milky Way. We reconstruct from mock exoplanet datasets observable parameters such as exoplanet age, temperature, mass, and location, together with DM halo parameters, for representative choices of measurement uncertainty and the number of exoplanets detected. We find that detection of $\mathcal{O}(100)$ exoplanets in the inner Galaxy can yield quantitative information on the galactic DM density profile, under the assumption of 10% measurement uncertainty. Even as few as $\mathcal{O}(10)$ exoplanets can deliver meaningful sensitivities if the DM density and inner slope are sufficiently large.

We build models where Dark Matter candidates arise as composite states of a new confining gauge force, stable thanks to accidental symmetries. Restricting to renormalizable theories compatible with SU(5) unification, we find 13 models based on SU(N) gauge theories and 9 based on SO(N). We also describe other models that require non-renormalizable interactions. The two gauge groups lead to distinctive phenomenologies: SU(N) theories give complex DM, with potentially observable electric and magnetic dipole moments that lead to peculiar spin-independent cross sections; SO(N) theories give real DM, with challenging spin-dependent cross sections or inelastic scatterings. Models with Yukawa couplings also give rise to spin-independent direct detection mediated by the Higgs boson and to electric dipole moments for the electron. In some models DM has higher spin. Each model predicts a specific set of lighter composite scalars, possibly observable at colliders.

We analyze nucleon decay modes in a no-scale supersymmetric flipped SU(5) GUT model, and contrast them with the predictions for proton decays via dimension-6 operators in a standard unflipped supersymmetric SU(5) GUT model. We find that these GUT models make very different predictions for the ratios $\Gamma(p \to \pi^0 \mu^+)/\Gamma(p \to \pi^0 e^+)$, $\Gamma(p \to \pi^+ \bar{\nu})/\Gamma(p \to \pi^0 e^+)$, $\Gamma(p \to K^0 e^+)/\Gamma(p \to \pi^0 e^+)$ and $\Gamma(p \to K^0 \mu^+)/\Gamma(p \to \pi^0 \mu^+)$, and that predictions for the ratios $\Gamma(p \to \pi^0 \mu^+)/\Gamma(p \to \pi^0 e^+)$ and $\Gamma(p \to \pi^+ \bar{\nu})/\Gamma(p \to \pi^0 e^+)$ also differ in variants of the flipped SU(5) model with normal- or inverse-ordered light neutrino masses. Upcoming large neutrino experiments may have interesting opportunities to explore both GUT and flavour physics in proton and neutron decays.

A novel gauge-theoretic framework for spacetime and gravitation is proposed, in which a {\it Cartan khronon} field dynamically breaks the symmetry between space and time, enabling the emergence of temporality within a fundamentally Euclidean setting. Based on a $Spin(4)$ gauge structure, the theory provides a real-valued formulation of chiral spacetime, wherein the effects typically attributed to dark matter could instead be explained by the dynamics of gravitation. New results are presented with implications for a broad range of phenomena, including cosmology, large-scale structure, gravitational waves, black holes, and potential signatures accessible to laboratory experiments. By avoiding the pitfalls of complexification and reinterpreting chiral spacetime geometry through a real, dynamical, and gauge-theoretic lens, the {\it Cartan khronon} theory offers a fresh and compelling framework for revisiting the foundations of spacetime and gravitation.

The spin-1/2 alternating Heisenberg chain system Na$_3$Cu$_2$SbO$_6$ features two relevant exchange couplings: $J_{1a}$ within the structural Cu$_2$O$_6$ dimers and $J_{1b}$ between the dimers. Motivated by the controversially discussed nature of $J_{1a}$, we perform extensive density-functional-theory (DFT) calculations, including DFT+$U$ and hybrid functionals. Fits to the experimental magnetic susceptibility using high-temperature series expansions and quantum Monte Carlo simulations yield the optimal parameters $J_{1a}$ = $-$217 K and $J_{1b}$ = 174 K with the alternation ratio $\alpha = J_{1a}/J_{1b} \simeq$$-$1.25. For the closely related system Na$_2$Cu$_2$TeO$_6$, DFT yields substantially enhanced $J_{1b}$, but weaker $J_{1a}$. The comparative analysis renders the buckling of the chains as the key parameter altering the magnetic coupling regime. Numerical simulation of the dispersion relations of the alternating chain model clarify why both antiferromagnetic and ferrromagnetic $J_{1a}$ can reproduce the experimental magnetic susceptibility data.

The large-scale structure of the Universe is characterised by a web-like structure made of voids, sheets, filaments, and knots. The structure of this so-called cosmic web is dictated by the local velocity shear tensor. In particular, the local direction of a filament should be strongly aligned with e3, the eigenvector associated with the smallest eigenvalue of the tensor. That conjecture is tested here on the basis of a cosmological simulation. The cosmic web delineated by the halo distribution is probed by a marked point process with interactions (the Bisous model), detecting filaments directly from the halo distribution (P-web). The detected P-web filaments are found to be strongly aligned with the local e3: the alignment is within 30 degree for 80% of the elements. This indicates that large-scale filaments defined purely from the distribution of haloes carry more than just morphological information, although the Bisous model does not make any prior assumption on the underlying shear tensor. The P-web filaments are also compared to the structure revealed from the velocity shear tensor itself (V-web). In the densest regions, the P- and V-web filaments overlap well (90%), whereas in lower density regions, the P-web filaments preferentially mark sheets in the V-web.