Preprints @ Department of Physics

The Circular Electron-Positron Collider (CEPC), a proposed next-generation Higgs factory, provides new opportunities to explore physics beyond the Standard Model (SM). With its clean electron-positron collision environment and the ability to collect large samples of Higgs, W, and Z bosons, the CEPC enables precision measurements and searches for new physics. This white paper outlines the CEPC's discovery potential, including studies of exotic decays of the Higgs, Z, and top quarks, dark matter and dark sector phenomena, long-lived particles, supersymmetry, and neutrino-related signatures. Advanced detector technologies and reconstruction techniques, such as one-to-one correspondence reconstruction and jet origin identification, significantly improve sensitivity to rare and weakly interacting processes. The CEPC is particularly well suited to probe the electroweak phase transition and test models of electroweak baryogenesis and dark sector interactions. In addition, global fit analyses highlight the CEPC's complementary role in constraining a wide range of new physics scenarios. These features position the CEPC as a powerful tool for exploring the next frontier in fundamental particle physics in the post-Higgs discovery era.

In collider physics, experiments are often based on counting the numbers of events in bins of a histogram. We present a new way to build and analyze statistical models that describe these experiments, based on the probabilistic programming language Stan and the HistFactory specification. A command-line tool transpiles HistFactory models into Stan code and data files. Because Stan is an imperative language, it enables richer and more detailed modeling, as modeling choices defined in the HistFactory declarative specification can be tweaked and adapted in the Stan language. Stan was constructed with automatic differentiation, allowing modern computational algorithms for sampling and optimization.

We rely on Monte Carlo (MC) simulations to interpret searches for new physics at the Large Hadron Collider (LHC) and elsewhere. These simulations result in noisy and approximate estimators of selection efficiencies and likelihoods. In this context we pioneer an exact-approximate computational method - exact-approximate Markov Chain Monte Carlo - that returns exact inferences despite noisy simulations. To do so, we introduce an unbiased estimator for a Poisson likelihood. We demonstrate the new estimator and new techniques in examples based on a search for neutralinos and charginos at the LHC using a simplified model. We find attractive performance characteristics - exact inferences are obtained for a similar computational cost to approximate ones from existing methods and inferences are robust with respect to the number of events generated per point.

In recent years, the prospect of detecting gravitational waves sourced from a strongly first-order cosmological phase transition has emerged as one of the most exciting frontiers of gravitational wave astronomy. Cosmological phase transitions are an essential ingredient in the Standard Model of particle cosmology, and help explain the mechanism for creation of matter in the early Universe, provide insights into fundamental theories of physics, and shed light on the nature of dark matter. This underscores the significance of developing robust end-to-end tools for determining the resulting gravitational waves from these phase transitions. In this article we present PhaseTracer2, an improved version of the C++ software package PhaseTracer, designed for mapping cosmological phases and transitions in Standard Model extensions of multiple scalar fields. Building on the robust framework of its predecessor, PhaseTracer2 extends its capabilities by including new features crucial for a more comprehensive analysis of cosmological phase transitions. It can calculate more complex properties, such as the bounce action through the path deformation method or an interface with BubbleProfiler, thermodynamic parameters, and gravitational wave spectra. Its applicability has also been broadened via incorporating the dimensionally reduced effective potential for models obtained from DRalgo, as well as calculations in the MSbar and OS-like renormalisation schemes. This modular, flexible, and practical upgrade retains the speed and stability of the original PhaseTracer, while significantly expanding its utility.

For several decades now, Bayesian inference techniques have been applied to theories of particle physics, cosmology and astrophysics to obtain the probability density functions of their free parameters. In this study, we review and compare a wide range of Markov Chain Monte Carlo (MCMC) and nested sampling techniques to determine their relative efficacy on functions that resemble those encountered most frequently in the particle astrophysics literature. Our first series of tests explores a series of high-dimensional analytic test functions that exemplify particular challenges, for example highly multimodal posteriors or posteriors with curving degeneracies. We then investigate two real physics examples, the first being a global fit of the $\Lambda$CDM model using cosmic microwave background data from the Planck experiment, and the second being a global fit of the Minimal Supersymmetric Standard Model using a wide variety of collider and astrophysics data. We show that several examples widely thought to be most easily solved using nested sampling approaches can in fact be more efficiently solved using modern MCMC algorithms, but the details of the implementation matter. Furthermore, we also provide a series of useful insights for practitioners of particle astrophysics and cosmology.

Building on previous work, we extend an algebraic realisation of three fermion generations within the complex Clifford algebra $\mathbb{C}\ell(8)$ by incorporating a $U(1)_{em}$ gauge symmetry. The algebra $\mathbb{C}\ell(8)$ corresponds to the algebra of complex linear maps from the (complexification of the) Cayley-Dickson algebra of sedenions, $\mathbb{S}$, to itself. Previous work represented three generations of fermions with $SU(3)_C$ colour symmetry permuted by an $S_3$ symmetry of order-three, but failed to include a $U(1)$ generator that assigns the correct electric charge to all states. Furthermore, the three generations suffered from a degree of linear dependence between states. By generalising the embedding of the discrete group $S_3$, corresponding to automorphisms of $\mathbb{S}$, into $\mathbb{C}\ell(8)$, we include an $S_3$-invariant $U(1)$ that correctly assigns electric charge. First-generation states are represented in terms of two even $\mathbb{C}\ell(8)$ semi-spinors, obtained from two minimal left ideals, related to each other via the order-two $S_3$ symmetry. The remaining two generations are obtained by applying the $S_3$ symmetry of order-three to the first generation. In this model, the gauge symmetries, $SU(3)_C\times U(1)_{em}$, are $S_3$-invariant and preserve the semi-spinors. As a result of the generalised embedding of the $S_3$ automorphisms of $\mathbb{S}$ into $\mathbb{C}\ell(8)$, the three generations are now linearly independent.

We uncover two precise interpretations of traditional electroweak fine-tuning (FT) measures that were historically missed. (i) a statistical interpretation: the traditional FT measure shows the change in plausibility of a model in which a parameter was exchanged for the $Z$ boson mass relative to an untuned model in light of the $Z$ boson mass measurement. (ii) an information-theoretic interpretation: the traditional FT measure shows the exponential of the extra information, measured in nats, relative to an untuned model that you must supply about a parameter in order to fit the $Z$ mass. We derive the mathematical results underlying these interpretations, and explain them using examples from weak scale supersymmetry. These new interpretations allow us to rigorously define FT in particle physics and beyond; shed fresh light on the status of extensions to the Standard Model; and lastly, allow us to precisely reinterpret historical and recent studies using traditional FT measures.

We compare the structures and methods in the theory of causal fermion systems with approaches to fundamental physics based on division algebras, in particular the octonions. We find that octonions and, more generally, tensor products of division algebras come up naturally to describe the symmetries of the vacuum configuration of a causal fermion system. This is achieved by associating the real and imaginary octonion basis elements with the neutrino and charged sectors of the vacuum fermionic projector, respectively. Conversely, causal fermion systems provide octonionic theories with spacetime structures and dynamical equations via the causal action principle. In this way, octonionic theories and causal fermion systems complement each other..

The Bayes factor surface is a new way to present results from experimental searches for new physics. Searches are regularly expressed in terms of phenomenological parameters - such as the mass and cross-section of a weakly interacting massive particle. Bayes factor surfaces indicate the strength of evidence for or against models relative to the background only model in terms of the phenomenological parameters that they predict. They provide a clear and direct measure of evidence, may be easily reinterpreted, but do not depend on choices of prior or parameterization. We demonstrate the Bayes factor surface with examples from dark matter, cosmology, and collider physics.

When exploring equations of nonlinear electrodynamics in effective medium formed by mutually parallel external electric and magnetic fields, we come to special static axial-symmetric solutions of two types. The first are comprised of fields referred to as electric and magnetic responses to a point-like electric charge when placed into the medium. In electric case, this is a field determined by the induced charge density. In magnetic case, this is a field carrying no magnetic charge and determined by an induced current. Fields of second type require presence of pseudoscalar constants for their existence. These are singular on the axis drawn along the external fields. In electric case this is a field of an inhomogeneously charged infinitely thin thread. In magnetic case this is the magnetic monopole with the Dirac string supported by solenoidal current. In both cases the necessary pseudoscalar constant is supplied by field derivatives of nonlinear Lagrangian taken on external fields. There is also a magnetic thread solution dual to electric thread with null total magnetic charge.

Effective field theory of low-energy exitations-magnons that describes antiferromagnets is mapped into scalar electrodynamics of a charged scalar field interacting with an external electromagnetic potential. In the presence of a constant inhomogeneous external magnetic field the latter problem is technically reduced to the problem of charged-particle creation from the vacuum by an electric potential step (x step). Magnetic moment plays here the role of the electric charge, and magnons and antimagnons differ from each other in the sign of the magnetic moment. In the framework of such a consideration, it is important to take into account the vacuum instability (the Schwinger effect) under the magnon-antimagnon production on magnetic field inhomogeneities (an analog of pair creation from the vacuum by electric-like fields). We demonstrate how to use the strong field QED with x steps developed by the authors (SPG and DMG) to study the magnon-antimagnon pair production on magnetic field inhomogeneities. Characteristics of the vacuum instability obtained for some magnetic steps that allows exact solving the Klein-Gordon equation are presented. In particular, we consider examples of magnetic steps with very sharp field derivatives that correspond to a regularization of the Klein step. In the case of smooth-gradient steps, we describe an universal behavior of the flux density of created magnon pairs. We also note that since the low-energy magnons are bosons with small effective mass, then for the first time maybe the opportunity will arise to observe the Schwinger effect in the case of the Bose statistics, in particular, the bosonic Klein effect in laboratory conditions. Moreover, it turns out that in the case of the Bose statistics appears a new mechanism for amplifying the effect of the pair creation, which we call statistically assisted Schwinger effect.

The BNL and FNAL measurements of the anomalous magnetic moment of the muon disagree with the Standard Model (SM) prediction by more than $4\sigma$. The hadronic vacuum polarization (HVP) contributions are the dominant source of uncertainty in the SM prediction. There are, however, tensions between different estimates of the HVP contributions, including data-driven estimates based on measurements of the $R$-ratio. To investigate that tension, we modeled the unknown $R$-ratio as a function of CM energy with a treed Gaussian process (TGP). This is a principled and general method grounded in data-science that allows complete uncertainty quantification and automatically balances over- and under-fitting to noisy data. Our tool yields exploratory results are similar to previous ones and we find no indication that the $R$-ratio was previously mismodeled. Whilst we advance some aspects of modeling the $R$-ratio and develop new tools for doing so, a competitive estimate of the HVP contributions requires domain-specific expertise and a carefully curated database of measurements.

Several pulsar timing array collaborations recently reported evidence of a stochastic gravitational wave background (SGWB) at nHz frequencies. Whilst the SGWB could originate from the merger of supermassive black holes, it could be a signature of new physics near the 100 MeV scale. Supercooled first-order phase transitions (FOPTs) that end at the 100 MeV scale are intriguing explanations, because they could connect the nHz signal to new physics at the electroweak scale or beyond. Here, however, we provide a clear demonstration that it is not simple to create a nHz signal from a supercooled phase transition, due to two crucial issues that could rule out many proposed supercooled explanations and should be checked. As an example, we use a model based on non-linearly realized electroweak symmetry that has been cited as evidence for a supercooled explanation. First, we show that a FOPT cannot complete for the required transition temperature of around 100 MeV. Such supercooling implies a period of vacuum domination that hinders bubble percolation and transition completion. Second, we show that even if completion is not required or if this constraint is evaded, the Universe typically reheats to the scale of any physics driving the FOPT. The hierarchy between the transition and reheating temperature makes it challenging to compute the spectrum of the SGWB.

Gravitational waves (GWs) were recently detected for the first time. This revolutionary discovery opens a new way of learning about particle physics through GWs from first-order phase transitions (FOPTs) in the early Universe. FOPTs could occur when new fundamental symmetries are spontaneously broken down to the Standard Model and are a vital ingredient in solutions of the matter anti-matter asymmetry problem. The purpose of our work is to review the path from a particle physics model to GWs, which contains many specialized parts, so here we provide a timely review of all the required steps, including: (i) building a finite-temperature effective potential in a particle physics model and checking for FOPTs; (ii) computing transition rates; (iii) analyzing the dynamics of bubbles of true vacuum expanding in a thermal plasma; (iv) characterizing a transition using thermal parameters; and, finally, (v) making predictions for GW spectra using the latest simulations and theoretical results and considering the detectability of predicted spectra at future GW detectors. For each step we emphasize the subtleties, advantages and drawbacks of different methods, discuss open questions and review the state-of-art approaches available in the literature. This provides everything a particle physicist needs to begin exploring GW phenomenology.

Using the GAMBIT global fitting framework, we constrain the MSSM with an eV-scale gravitino as the lightest supersymmetric particle, and the six electroweakinos (neutralinos and charginos) as the only other light new states. We combine 15 ATLAS and 12 CMS searches at 13\,TeV, along with a large collection of ATLAS and CMS measurements of Standard Model signatures. This model, which we refer to as the $\tilde G$-EWMSSM, exhibits quite varied collider phenomenology due to its many permitted electroweakino production processes and decay modes. Characteristic $\tilde G$-EWMSSM signal events have two or more Standard Model bosons and missing energy due to the escaping gravitinos. While much of the $\tilde G$-EWMSSM parameter space is excluded, we find several viable parameter regions that predict phenomenologically rich scenarios with multiple neutralinos and charginos within the kinematic reach of the LHC during Run 3, or the High Luminosity LHC. In particular, we identify scenarios with Higgsino-dominated electroweakinos as light as 140 GeV that are consistent with our combined set of collider searches and measurements. The full set of $\tilde G$-EWMSSM parameter samples and GAMBIT input files generated for this work is available via Zenodo.

We address the problem of the electromagnetic radiation produced by charge distributions in the framework of a semiclassical approach proposed in the work by Bagrov, Gitman, Shishmarev and Farias [J. Synchrotron Rad. (2020). 27, 902-911]. In this approach, currents, generating the radiation are considered classically, while the quantum nature of the radiation is kept exactly. Quantum states of the electromagnetic field are solutions of Schr\"odinger's equation and relevant quantities to the problem are evaluated with the aid of transition probabilities. This construction allows us to introduce the quantum transition time in physical quantities and assess its role in radiation problems by classical currents. We study radiated electromagnetic energies in detail and present a definition for the rate at which radiation is emitted from sources. In calculating the total energy and rate radiated by a pointlike charged particle accelerated by a constant and uniform electric field, we discover that our results are compatible with results obtained by other authors in the framework of the classical radiation theory under an appropriate limit. We also perform numerical and asymptotic analysis of the results.

We explore the origins of parameters in adimensional theories -- fundamental theories with no classical massive scales. If the parameters originate as draws from a distribution, it should be possible to write a distribution for them that doesn't depend on or introduce any massive scales. These distributions are the invariant distributions for the renormalization group (RG). If there exist RG invariant combinations of parameters, the RG invariant distributions are specified up to arbitrary functions of the RG invariants. If such distributions can be constructed, adimensional theories could predict the values of their parameters through distributions that are constrained by the RG. If they can't be constructed, the parameters must originate in a different way. We demonstrate the RG invariant distributions in QCD, the Coleman-Weinberg model and a totally asymptotically free theory.

Nested sampling (NS) is a popular algorithm for Bayesian computation. We investigate statistical errors in NS both analytically and numerically. We show two analytic results. First, we show that the leading terms in Skilling's expression using information theory match the leading terms in Keeton's expression from an analysis of moments. This approximate agreement was previously only known numerically and was somewhat mysterious. Second, we show that the uncertainty in single NS runs approximately equals the standard deviation in repeated NS runs. Whilst intuitive, this was previously taken for granted. We close by investigating our results and their assumptions in several numerical examples, including cases in which NS uncertainties increase without bound.

We investigate the extent to which perturbative calculations of the electroweak phase transition are arbitrary and uncertain, owing to their gauge, renormalisation scale and scheme dependence, as well as treatments of the Goldstone catastrophe and daisy diagrams. Using the complete parameter space of the Standard Model extended by a real scalar singlet with a $\mathbb{Z}_2$ symmetry as a test, we explore the properties of the electroweak phase transition in general $R_\xi$ and covariant gauges, OS and $\overline{\text{MS}}$ renormalisation schemes, and for common treatments of the Goldstone catastrophe and daisy diagrams. Reassuringly, we find that different renormalisation schemes and different treatments of the Goldstone catastrophe and daisy diagrams typically lead to only modest changes in predictions for the critical temperature and strength of the phase transition. On the other hand, the gauge and renormalisation scale dependence may be significant, and often impact the existence of the phase transition altogether.

We review Skilling's nested sampling (NS) algorithm for Bayesian inference and more broadly multi-dimensional integration. After recapitulating the principles of NS, we survey developments in implementing efficient NS algorithms in practice in high-dimensions, including methods for sampling from the so-called constrained prior. We outline the ways in which NS may be applied and describe the application of NS in three scientific fields in which the algorithm has proved to be useful: cosmology, gravitational-wave astronomy, and materials science. We close by making recommendations for best practice when using NS and by summarizing potential limitations and optimizations of NS.