We present the results of Phase I of an ongoing review of Monte Carlo tools relevant for low-energy hadronic cross sections. This includes a detailed comparison of Monte Carlo codes for electron-positron scattering into a muon pair, pion pair, and electron pair, for scan and radiative-return experiments. After discussing the various approaches that are used and effects that are included, we show differential cross sections obtained with AfkQed, BabaYaga@NLO, KKMC, MCGPJ, McMule, Phokhara, and Sherpa, for scenarios that are inspired by experiments providing input for the dispersive evaluation of the hadronic vacuum polarisation.

This paper derives a next-to-leading power (NLP) soft theorem for multi-photon emission to all orders in the electromagnetic coupling constant, generalising the leading-power theorem of Yennie, Frautschi, and Suura. Working in the QED version of heavy-quark effective theory, multi-emission amplitudes are shown to reduce to single- and double-radiation contributions only. Single soft-photon emission, in turn, is described by the recent all-order extension of the Low-Burnett-Kroll theorem, where the tree-level formula is supplemented with a one-loop exact soft function. The same approach is used in this article to prove that the genuine double-emission contribution is tree-level exact. As a validation and a first non-trivial application of the multi-photon theorem, the real-real-virtual electron-line corrections to muon-electron scattering are calculated at NLP in the soft limit.

We provide an overview of the status of Monte-Carlo event generators for high-energy particle physics. Guided by the experimental needs and requirements, we highlight areas of active development, and opportunities for future improvements. Particular emphasis is given to physics models and algorithms that are employed across a variety of experiments. These common themes in event generator development lead to a more comprehensive understanding of physics at the highest energies and intensities, and allow models to be tested against a wealth of data that have been accumulated over the past decades. A cohesive approach to event generator development will allow these models to be further improved and systematic uncertainties to be reduced, directly contributing to future experimental success. Event generators are part of a much larger ecosystem of computational tools. They typically involve a number of unknown model parameters that must be tuned to experimental data, while maintaining the integrity of the underlying physics models. Making both these data, and the analyses with which they have been obtained accessible to future users is an essential aspect of open science and data preservation. It ensures the consistency of physics models across a variety of experiments.

This report presents a short summary of the activities of the "Standard Model" working group for the "Physics at TeV Colliders" workshop (Les Houches, France, 12-30 June, 2023).

We study the soft limit of one-photon radiation at next-to-leading power (NLP) in the framework of heavy-quark effective theory (HQET) to all orders in perturbation theory. We establish the soft theorem that for unpolarised scattering the radiative contribution up to NLP is entirely determined by the non-radiative amplitude. This generalises the Low-Burnett-Kroll (LBK) theorem for QED to all orders. All hard matching corrections can be calculated by applying the LBK differential operator to the non-radiative amplitude. The virtual corrections in the effective theory vanish beyond one loop, resulting in a one-loop exact soft function. As a first, non-trivial application we calculate the real-virtual-virtual electron-line corrections to muon-electron scattering at NLP in the soft limit.

We present an improved theoretical prediction of the positron energy spectrum for the polarised Michel decay \mu^+\to e^+ \nu_e\bar{\nu}_\mu μ+→e+νeν‾μ . In addition to the full next-to-next-to-leading order correction of order \alpha^2 α2 in the electromagnetic coupling, we include logarithmically enhanced terms at even higher orders. Logarithms due to collinear emission are included at next-to-leading accuracy up to order \alpha^4 α4 . At the endpoint of the Michel spectrum, soft photon emission results in large logarithms that are resummed up to next-to-next-to-leading logarithmic accuracy. We apply our results in the context of the MEG II and Mu3e experiments to estimate the impact of the theory error on the branching ratio sensitivity for the lepton-flavour-violating decay \mu^+\to e^+ X μ+→e+X of a muon into an axion-like particle X X .

We present the complete next-to-next-to-leading order (NNLO) pure pointlike QED corrections to lepton-proton scattering, including three-photon-exchange contributions, and investigate their impact in the case of the MUSE experiment. These corrections are computed with no approximation regarding the energy of the emitted photons and taking into account lepton-mass effects. We contrast the NNLO QED corrections to known next-to-leading order corrections, where we include the elastic two-photon exchange (TPE) through a simple hadronic model calculation with a dipole ansatz for the proton electromagnetic form factors. We show that, in the low-momentum-transfer region accessed by the MUSE experiment, the improvement due to more sophisticated treatments of the TPE, including inelastic TPE, is of similar if not smaller size than some of the NNLO QED corrections. Hence, the latter have to be included in a precision determination of the low-energy proton structure from scattering data, in particular for electron-proton scattering. For muon-proton scattering, the NNLO QED corrections are considerably smaller.

Pulak Banerjee, Tim Engel, Nicolas Schalch, Adrian Signer, Yannick Ulrich a Paul Scherrer Institut, CH-5232 Villigen PSI, Switzerland b Physik-Institut, Universität Zürich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland c Albert Einstein Center for Fundamental Physics, Institut für Theoretische Physik, Universität Bern, Sidlerstrasse 5, CH-3012 Bern, Switzerland d Institute for Particle Physics Phenomenology, University of Durham, South Road, Durham DH1 3LE, United Kingdom

This white paper concerns theoretical and phenomenological aspects relevant to the physics of future $e^+e^-$ colliders, in particular regarding initial-state QED radiation. The contributions each contain key technical aspects, and are formulated in a pedagogical manner so as to render them accessible also to those who are not directly working on these and immediately-related topics. This should help both experts and non-experts understand the theoretical challenges that we shall face at future $e^+e^-$ colliders. Specifically, this paper contains descriptions of the treatment of initial state radiation from several Monte Carlo collaborations, as well as contributions that explain a number of more theoretical developments with promise of future phenomenological impact.

We present two novel results about the universal structure of radiative QED amplitudes in the soft and in the collinear limit. On the one hand, we extend the well-known Low-Burnett-Kroll theorem to the one-loop level and give the explicit relation between the radiative and non-radiative amplitude at subleading power in the soft limit. On the other hand, we consider a factorisation formula at leading power in the limit where the emitted photon becomes collinear to a light fermion and provide the corresponding one-loop splitting function. In addition to being interesting in their own right these findings are particularly relevant in the context of fully-differential higher-order QED calculations. One of the main challenges in this regard is the numerical stability of radiative contributions in the soft and collinear regions. The results presented here allow for a stabilisation of realvirtual amplitudes in these delicate phase-space regions by switching to the corresponding approximation without the need of explicit computations.

McMule is a generic framework for higher-order QED calculations of scattering and decay processes involving leptons. In this dataset we include all data used in McMule publications as well as the analysis code used. An up-to-date version can be found here.

A critical subject in fully differential QED calculations originates from numerical instabilities due to small fermion masses that act as regulators of collinear singularities. At next-to-next-to-leading order (NNLO) a major challenge is therefore to find a stable implementation of numerically delicate real-virtual matrix elements. In the case of Bhabha scattering this has so far prevented the development of a fixed-order Monte Carlo at NNLO accuracy. In this paper we present a new method for stabilising the real-virtual matrix element. It is based on the expansion for soft photon energies including the non-universal subleading term calculated with the method of regions. We have applied this method to Bhabha scattering to obtain a stable and efficient implementation within the McMule framework. We therefore present for the first time fully differential results for the photonic NNLO corrections to Bhabha scattering.

We review the current status of the theory predictions for elastic $$\mu $$-e scattering, describing the recent activities and future plans of the theory initiative related to the proposed MUonE experiment.

We review the current status of the theory predictions for elastic $\mu$-$e$ scattering, describing the recent activities and future plans of the theory initiative related to the proposed MUonE experiment.

McMule is a framework for fully differential higher-order QED calculations of scattering and decay processes involving leptons. It keeps finite lepton masses, which regularises collinear singularities. Soft singularities are treated with dimensional regularisation and using FKS$^\ell$ subtraction. We describe the implementation of the framework in Fortran 95, list the processes that are currently implemented, and give instructions on how to run the code. In addition, we present new phenomenological results for muon-electron scattering and lepton-proton scattering, including the dominant NNLO corrections. While the applications presented focus on MUonE, MUSE, and P2, the code can be used for a large number of planned and running experiments.

We present the heavy-to-light form factors with two different non-vanishing masses at next-to-next-to-leading order and study its expansion in the small mass. The leading term of this small-mass expansion leads to a factorized expression for the form factor. The presence of a second mass results in a new feature, in that the soft contribution develops a factorization anomaly. This cancels with the corresponding anomaly in the collinear contribution. With the generalized factorization presented here, it is possible to obtain the leading small-mass terms for processes with large masses, such as muon-electron scattering, from the corresponding massless amplitude and the soft contribution.