Direct measurement of alpha_QED(mZ) at the FCC-ee

When the measurements from the FCC-ee become available, an improved determination of the standard-model"input"parameters will be needed to fully exploit the new precision data towards either constraining or fitting the parameters of beyond-the-standard-model theories. Among these input parameters is the electromagnetic coupling constant estimated at the Z mass scale, alpha_QED(mZ). The measurement of the muon forward-backward asymmetry at the FCC-ee, just below and just above the Z pole, can be used to make a direct determination of alpha_QED(mZ) with an accuracy deemed adequate for an optimal use of the FCC-ee precision data.


Introduction
The design study of the Future Circular Colliders (FCC) in a 100-km ring in the Geneva area has started at CERN at the beginning of 2014, as an option for post-LHC particle accelerators. The study has an emphasis on proton-proton and electron-positron high-energy frontier machines [1]. In the current plans, the first step of the FCC physics programme would exploit a high-luminosity e + e − collider called FCC-ee, with centre-of-mass energies ranging from the Z pole to the tt threshold and beyond. A first look at the physics case of the FCC-ee can be found in Ref. [2].
In this first look, an estimate of the achievable precision on a number of Z-pole observables was inferred and used in a global electroweak fit to set constraints on weakly-coupled new physics up to a scale of 100 TeV [3]. These constraints were obtained under two assumptions: (i) the precision of the pertaining theoretical calculations will match the expected experimental accuracy by the time of the FCC-ee startup; and (ii) the determination of standard-model input parameters -four masses: m Z , m W , m top , m Higgs ; and three coupling constants: α s (m 2 Z ), G F , α QED (m 2 Z ) -will improve in order not to be the limiting factors to the constraining power of the fit. The determinations of the Higgs boson mass from the LHC data [4] and of the Fermi constant from the muon lifetime measurement [5] are already sufficient for this purpose. It is argued in Refs. [2,6] that the FCC-ee can adequately improve the determination of the other three masses and of the strong coupling constant by one order of magnitude or more: the experimental precision targets for the FCC-ee are 100 keV for the Z-boson mass, 500 keV for the W-boson mass, 10 MeV for the top-quark mass, and 0.0001 for the strong coupling constant. (The FCC-ee also aims at reducing the Higgs boson mass uncertainty down to 8 MeV.) No mention was made, however, of a way to improve the determination of the electromagnetic coupling constant evaluated at the Z mass, and it was simply assumed that a factor 5 improvement with respect to today's uncertainty -down to 2 × 10 −5 -could be achieved by the time of the FCC-ee startup. Today, α QED (m 2 Z ) is determined from α QED (0) (itself known with an accuracy of 10 −10 ) with the running coupling constant formula: had (m 2 Z ) .
(1.1) Its uncertainty is dominated by the experimental determination of the hadronic vacuum polarization, ∆α (5) had (m 2 Z ), obtained from the dispersion integral: where R γ (s) is the hadronic cross section σ 0 (e + e − → γ * → hadrons) at a given centreof-mass energy √ s, normalized to the muon pair cross section at the same centre-of-mass energy. At small values of √ s, typically up to 5 GeV, and in the Υ resonance region from 9.6 to 13 GeV, the evaluation of the dispersion integral relies on the measurements made with low-energy e + e − data accumulated by the KLOE, CMD-2/SND, BaBar, Belle, CLEO and BES experiments. The most recent re-evaluation [7,8] gives ∆α (5) had (m 2 Z ) = (275.7 ± 1.0) × 10 −4 , which leads to corresponding to a relative uncertainty on the electromagnetic coupling constant, ∆α/α, of 1.1 × 10 −4 . It is hoped that future low-energy e + e − data collected by the BES III and VEPP-2000 colliders will improve this figure to 5 × 10 −5 or better [9]. In this study, it is shown that the FCC-ee can provide another way of determining the electromagnetic coupling constant with a similar or better accuracy, from the precise measurement of muon forward-backward asymmetry, A µµ FB , just above and just below the Z peak, as part of the resonance scan. This method does not rely on the experimental determination of the vacuum polarization ∆α (5) had . Here, the point is not to extrapolate α QED (m 2 Z ) from α QED (0), but to provide a direct evaluation of α QED at √ s m Z , hence with totally different theoretical and experimental uncertainties. This measurement would in turn be combined with other determinations for an even smaller uncertainty.
This letter is organized as follows. In Section 2, the reasons for the choice of A µµ FB as an observable sensitive to α QED are given, and the sensitivity is determined as a function of the centre-of-mass energy. The optimal centre-of-mass energies, as well as the integrated luminosities and running time needed to achieve a statistical uncertainty of a few 10 −5 are determined in Section 3. Possible systematic uncertainties are discussed and evaluated in Section 4.
2 The muon forward-backward asymmetry and the electromagnetic coupling constant At the FCC-ee, the muon pair production proceeds via the graph depicted in Fig. 1 through either a Z or a γ exchange. At tree level, the cross section σ µµ therefore contains three terms: (i) the γ-exchange term squared, proportional to α 2 QED (s); (ii) the Z-exchange term squared, proportional to G 2 F (where G F is the Fermi constant); and (iii) the γ-Z interference term, proportional to α QED (s) × G F . These three terms are denoted G, Z, and I in the following. Their expressions as a function of the centre-of-mass energy √ s can be found in Ref. [10] and reported below.
with the following definitions: (2.4) and where θ W is the effective Weinberg angle (sin 2 θ W 0.2315).
An absolute measurement of the µ + µ − production cross section σ µµ = Z + I + G is therefore a priori sensitive to α QED through the interference term and the γ-exchange term. The cross section and the three contributing terms are displayed in Fig. 2 as a function of the centre-of-mass energy √ s, with the inclusion of initial state radiation (ISR). In this figure, the effective collision energy after ISR, denoted √ s , is required to satisfy s 0.99s. The importance of such a requirement on s , together with the way to control it experimentally, is discussed in Section 4.3.2.
At a given √ s, a small variation ∆α of the electromagnetic coupling constant translates to a variation ∆σ µµ of the cross section : As is well visible in Fig. 2, the interference term can be neglected in the above equation.
As a consequence, if the cross section can be measured with a precision ∆σ µµ , the relative precision on the electromagnetic coupling constant amounts to The target statistical precision of 2 × 10 −5 on α QED can therefore be achieved with more than 10 9 µ + µ − events and at centre-of-mass energy where the Z contribution to the cross section is much smaller than the photon contribution. These two conditions call for a centre-of-mass energy smaller than 70 GeV, where the cross section is both large and dominated by the photon contribution. Beside the fact that this centre-of-mass energy is not in the current core programme of the FCC-ee and that the needed integrated luminosity of 50 ab −1 would require at least a year of running at this energy in the most favourable conditions, the measurement itself poses a number of intrinsic difficulties. Indeed, the absolute measurement of a cross section with a precision of a few 10 −5 requires the selection efficiency, the detector acceptance, and the integrated luminosity to be known with this precision or better. Even if not impossible to meet, these requirements are exceedingly challenging in the extraction of α QED from this method with the needed precision.
The muon forward-backward asymmetry, A µµ FB , defined as where σ F(B) µµ is the µ + µ − cross section for events with the µ − direction in the forward (backward) hemisphere with respect to the e − -beam direction, hence with σ F µµ + σ B µµ = σ µµ , solves most of these obstacles. Indeed, it is a self-normalized quantity, which thus does not  Cross section for the e + e − → µ + µ − process (red curve) and the three contributions, calculated from the analytical expressions of Ref. [10]: pure γ-exchange term (blue curve); pure Z-exchange term (green curve); and the absolute value of the γ-Z interference term (black curve). The initial-state radiation is included, and s /s is required to exceed 0.99.
need the measurement of the integrated luminosity. Moreover, most uncertainties on the selection efficiency and the detector acceptance simply cancel in the ratio. This observable is therefore a good candidate for a measurement with an exquisite precision. At lowest order, and if the terms proportional to m 2 µ /m 2 Z ∼ 10 −6 are neglected, the angular distribution of the µ − from the e + e − → µ + µ − production can be written in the following way [11]: where G 1 (s) and G 3 (s) can be expressed as a function of G, Z and I as follows: After integration over the muon polar angle θ, the forward-backward asymmetry therefore amounts to: (2.10) The variation of A µµ FB as a function of the centre-of-mass energy, as obtained from Eq. 2.10, is shown in Fig. 3. In the above expressions, the photon-exchange term is totally symmetric, hence is absent from the numerator. Because v 4 /a 4 3 × 10 −5 , the Z-exchange term contribution to the asymmetry is minute, except at the Z pole where the interference term vanishes and the asymmetry is small: A µµ FB,0 = (3/4) × 4v 2 a 2 /(a 2 + v 2 ) 2 0.016. The interference term, on the other hand, is almost 100% anti-symmetric and contributes mostly to the numerator. (The contribution of the interference term to the denominator, i.e., to the total cross section, can be neglected as shown in Fig. 2.) The off-peak muon forward-backward asymmetry can therefore be expressed as follows: At a given √ s, a small variation ∆α of the electromagnetic coupling constant translates to a variation ∆A µµ F B of the muon forward-backward asymmetry: In first approximation, the asymmetry is therefore not sensitive to α QED when the Zand photon-exchange terms are equal, i.e., at √ s = 78 and 112 GeV (Fig. 2), where the asymmetry is maximal (Fig. 3). Similarly, the sensitivity to the electromagnetic coupling constant vanishes in the immediate vicinity of the Z pole. The red curve of Fig. 4 shows the variation of A µµ FB for a relative change of α QED by +1.1 × 10 −4 , as a function of √ s. In other words, the red curves displays the absolute precision with which A µµ FB must be measured to start improving the accuracy on α QED (m 2 Z ) with respect to today's determination.  to α QED when the red curve crosses the black horizontal line. The blue area represents the absolute statisitical uncertainty with which the muon forward-backward asymmetry can be measured at the FCC-ee in one year of data taking at any given centre-of-mass energy.
For a positive variation of ∆α, the sign of ∆A µµ FB , i.e., the sign of A µµ , changes at each of these centre-of mass energies: it is positive below 78 GeV, where the asymmetry is negative and the Z contribution is smaller than the photon contribution, becomes negative between 78 GeV and the Z pole, where the Z contribution dominates, then positive again from the Z pole all the way to 112 GeV because the asymmetry becomes positive, and negative for larger centre-of-mass energies where the photon contribution takes over. This interesting property, in particular the sign change around the Z pole, is fully exploited in Section 4. Written the other way around and in a perhaps more useful manner for the following, the relative precision on the electromagnetic coupling constant amounts where the approximation in the last term of the equality is valid off the Z peak.

Statistical power of the method
The optimal centre-of-mass energies are those which minimize the statistical uncertainty on α QED (s). For a given integrated luminosity L, the statistical uncertainty on the forwardbackward asymmetry amounts to The target luminosities for the FCC-ee in a configuration with four interaction points are 215 × 10 34 cm −2 s −1 per interaction point at the Z pole and 38 × 10 34 cm −2 s −1 per interaction point at the WW pair production threshold [12]. With 10 7 effective seconds per year, the total integrated luminosity is therefore expected to be 86 ab −1 / year at the Z pole and 15.2 ab −1 / year at the WW threshold. Between these two points, the variation of the luminosity with the centre-of-mass energy is assumed to follow a simple power law: The very large Z pole luminosity is achieved by colliding about 60,000 bunches of electrons and positrons, which fill the entirety of the 400 MHz RF buckets available over 100 km. It also corresponds to a time between two bunch crossings of 5 ns, which is close to the minimum value acceptable today for the experiments. With a constant number of bunches, the luminosity was therefore conservatively assumed to linearly decrease with the centre-of-mass energy (and reach 0. for √ s = 0.), leading to the profile of  Fig. 4, for a one-year running at any given centre-of-mass energy. An improvement on the determination of α QED (s) is possible wherever the red curve lies outside the blue area, and is largest when the absolute value of the ratio between the red and blue curves is maximum.
The corresponding relative accuracy for the α QED (s) determination is shown in Fig. 6. The best accuracy of ∼ 3 × 10 −5 is obtained for one year of running either just below or just above the Z pole, specifically at √ s − ∼ 87.9 GeV and √ s + ∼ 94.3 GeV.
The value of the electromagnetic coupling constant extracted from the muon forwardbackward asymmetry measured at either energy, α − ≡ α QED (s − ) and α + ≡ α QED (s + ), are then extrapolated towards a determination of α 0 ≡ α QED (m 2 Z ) with the running coupling constant expression around the Z pole, valid at all orders in the leading-log approximation: where β is proportional to the well-known QED β-function. In the standard model and at the lowest QED/QCD order, it reads β 0 = f Q 2 f /3π, where the sum runs over all active fermions at the Z pole (f = e, µ, τ , d, u, s, c b) and Q f is the fermion electric charge in unit of e. The standard model extrapolation correction from α ± to α 0 therefore amounts to −0.033 from the measurement below the Z pole, and +0.030 from the measurement above the Z pole, corresponding to a relative correction of ±2.5 × 10 −4 in both cases, i.e., an order of magnitude larger than the targeted uncertainty on α 0 . While this correction is known with an excellent precision in the standard model -the QED β-function is now known with QED corrections up to five loops and QCD corrections up to four loops [13,14] -, it is certainly preferable to remove this model dependence (and the residual theory uncertainty) from the determination of α 0 . The dual measurements of α − and of α + solve this issue and yields the straightforward combination: without any model dependence related to the running of the electromagnetic constant. This combination of a measurement below the Z peak and a measurement above the Z peak has other advantages, the most important of which is the cancellation to a large extent of many systematic uncertainties, as explained in the next section. With this weighted average, the targeted precision of 2 × 10 −5 can be obtained from one year at 87.9 GeV and one year  Figure 6. Relative statistical uncertainty for the α QED (s) determination from a measurement of the muon forward-backward asymmetry at the FCC-ee, with a one-year running at any given centre-of-mass energy. The best accuracy is obtained for one year of running either just below or just above the Z pole. at 94.3 GeV with the sole measurement of the muon forward-backward asymmetry. The running time at each energy can be reduced to six months -as is assumed in the following -if additional measurements are considered, e.g., the tau forward-backward asymmetry and, possibly, the electron forward-backward asymmetry.

Systematic uncertainties
For the forward-backward asymmetry measurement to be relevant in the determination of the electromagnetic coupling constant, all possible systematic uncertainties must be kept well below the statistical uncertainty aimed at in Section 3. Systematic uncertainties may be of experimental, parametric, and theoretical origin, and are studied in turn below.

Beam energy calibration
Circular e + e − colliders have the unique feature of providing the possibility to measure the beam energy with the resonant depolarization method [15] with an outstanding accuracy.
At the FCC-ee [16], this accuracy has been estimated [2,17] to be of the order of 50 keV around the Z pole, of which 45 (23) keV are (un)correlated between all energy points, corresponding to a total relative uncertainty of 10 −6 . The derivative of the muon forwardbackward asymmetry with respect to the centre-of-mass energy, however, is largest around the Z pole, as can be seen from Fig. 3. It is therefore important to check that this expected precision is indeed sufficient.
At √ s − and √ s + , the photon contribution is only 5% of the total cross section and varies slowly with the centre-of-mass energy around the Z pole (Fig. 2): this contribution can be considered as a second order effect in the uncertainty evaluation. Equations 2.11 and 2.13 therefore simplify to (4.1) The dependence of I and Z on s and m Z is given at the beginning of Section 2. The forward-backward asymmetry dependence on s and m Z in the vicinity of the Z pole is simply which in turn simplifies to, at √ s ± , after neglecting the much smaller term proportional to (σ Σ /Σ) 2 . Numerically, the relative uncertainties on α ± , or equivalently on 1/α ± , arising from the beam energy measurement both amount to 1.4 × 10 −5 and are uncorrelated. The uncertainty on the coefficient ξ (±0.00001) was found to have a totally negligible contribution (±2 × 10 −9 ) to the relative uncertainty on α 0 . Only the (uncorrelated) errors on α − and α + contribute. As a consequence, the relative uncertainty on α QED (m 2 Z ) arising from the beam energy measurement amounts to (4.5)

Beam energy spread
At the FCC-ee, the relative beam energy spread δ for centre-of-mass energies around the Z pole is expected [12] to be of the order of 0.12%, i.e., two orders of magnitude larger than the accuracy of the (average) beam energy measurement. The relative centre-of-mass energy spread δ is √ 2 times smaller, i.e., of the order of 0.08%. The shift ∆A µµ FB between the predicted asymmetry and its measured value at a centre-of-mass energy √ s ± is therefore which yields, with the functional form of A µµ FB (s) given in Eq. 4.2 expanded around s ± : i.e., numerically under the reasonable assumption that the beam energy spread values are similar at √ s ± and m Z . The relative changes of A FB (s ± ) are of the order of the statistical uncertainty, and larger than the uncertainty originating from the beam energy measurement. These changes are, however, of opposite sign, and lead to a remarkable cancellation by more than one order of magnitude in the determination of α 0 . Indeed, the combination of Eqs. 2.13 and 3.3 leads to the following estimate of the bias on α 0 : The uncertainty on this small bias (which is to be corrected for) depends on the accuracy with which the beam energy spread in known. For example, the measurement of bunch length from the distribution of the µ + µ − event primary vertices determined directly by the FCC-ee experiments would allow a precise determination of the beam energy spread. A precision of 2.5% could be reached with this method at LEP [18], yielding a negligible uncertainty on the α QED (m 2 Z ) determination.

Muon identification efficiency and detector acceptance
In Eq. 2.10, the asymmetry is determined under the assumption of a 100% muon identification efficiency and a 4π detector acceptance. This equation is still valid for a smaller efficiency, with the condition that it is independent of the muon polar angle. If instead the identification efficiency times the detector acceptance is a non-trivial function of the polar angle, ε(cos θ), the measured muon angular distribution gets modified accordingly, and so does the measured forward-backward asymmetry. This issue can be solved experimentally with the observation [19] that a e + e − → µ + µ − event contains not only a negative muon but also a positive muon, the measured angular distributions of which are given by Eqs. 2.8 and 2.10 modified with ε(cos θ): where ε(cos θ) is assumed to be independent on the muon electric charge. This very reasonable assumption can be verified with an adequate accuracy from the 3 × 10 11 Z → µ + µ − events collected at √ s = 91.2 GeV, by tagging one of the two muons in each event, and probing the other to determine ε(cos θ) separately for positive and negative muons. The ratio of the difference to the sum of the numbers of negative and positive muons detected a given cos θ bin, N − (cos θ) and N + (cos θ), therefore amounts to which allows the muon forward-backward asymmetry to be determined in each bin as follows: an expression from which ε(cos θ) has simplified away in the ratio, hence without any impact on the measurement uncertainty. The muon forward-backward asymmetry for the complete event sample is then obtained by the statistically-weighted average of the bin-bybin determination over the detector acceptance. Any systematic effect related to the choice of the bin size -related for example to the muon angular resolution -can be eliminated by the use of an unbinned likelihood instead.

Charge inversion
The asymmetry of the sample of events where both muon charges are wrongly measured equals −A µµ FB . The relative change of the asymmetry arising from double charge inversion therefore amounts to ∆A µµ where f ± is the probability for a muon to be measured with the wrong charge sign. For this effect to be relevant (i.e., larger than 2 × 10 −5 ), f ± would need to exceed 0.5% -a typical value for LEP detectors. If FCC-ee detectors were similar to LEP detectors, f ± would be measured with an outstanding precision from the several 10 7 million µ ± µ ± events collected at √ s ± , thus allowing the effect to be corrected with no additional uncertainty on the asymmetry. On the other hand, the next generation of detectors for FCC-ee is likely to provide a track momentum resolution better than that of the LEP detectors by up to one order of magnitude, reducing f ± to ridiculously small values, with no sizable effect on the asymmetry.

Background from tau-pair production
The background from e + e − → τ + τ − , where the two taus decay into µν e ν µ has a cross section of the order of 3% of the e + e − → µ + µ − cross section. It can be greatly reduced by cuts on the muon impact parameters, on the angle between the two muons, and on the muon momenta. The very small residual contribution from this process is however not an issue, as the tau forward-backward asymmetry is expected to be identical to the muon forward-backward asymmetry. No additional uncertainty is therefore expected from this source.

Parametric uncertainties
The cross section given in Section 2 depends solely on four parameters -beside α QED (s)namley the Fermi constant G F , the Z boson mass and width, m Z and Γ Z , and the Weinberg angle, sin 2 θ W . The precision with which these parameters are known is the source of additional uncertainties for the muon forward-backward asymmetry, and in turn, on the electromagnetic coupling constant. These uncertainties are examined in turn below.

Z mass
The uncertainty on the Z mass is fully correlated to the uncertainty on the beam energy. Its effect on the forward backward asymmetry is already taken into account in Section 4.1.1.

Z width
The Z width simplifies away from the ratio given in Eq. 4.1, which contains only the dominant contributions from I and Z to the asymmetry. To exhibit the sub-leading dependence on Γ Z , it is necessary to go back to the more complete expression given in Eq. 2.11, which contains also the G term. The uncertainty arising from the knowledge of Γ Z is therefore not expected to be dominant. Straightforward error propagation yields (4.14) The uncertainty on Γ Z is dominated by the energy calibration error and amounts to σ D = 46 keV, i.e., to about 2 × 10 −5 Γ Z . At √ s ± , the photon contribution G is about 5% of the Z contribution Z, itself about 50% of its pole value. As a consequence, the relative uncertainties on α ± are equal and amount to 10 −6 with a 100% correlation. The relative uncertainty on α QED (m 2 Z ) arising from the Z width is therefore at the same level of 10 −6 .

Weinberg angle
Only the terms proportional to I and Z in the complete asymmetry expression (Eqs. 2.9 and 2.10) vary with sin 2 θ W , through the vectorial coupling v. In the vicinity of the Z pole, the small photon contribution can anyway be neglected, and the asymmetry expression simplifies to (4.15) The derivative of A µµ FB (s) with respect to sin 2 θ W can be obtained analytically, yielding in a straightforward manner The propagation to α 0 from Eq. 4.9 leads to a partial cancellation by almost one order of magnitude: For the current precision of the effective Weinberg angle determination, of the order 1.6 × 10 −4 , this uncertainty on α QED is large and amounts to 1.4 × 10 −4 . At the FCC-ee, however, the measurement of the asymmetry at the Z pole (insensitive to the electromagnetic coupling constant) can be used to improve the precision of the effective Weinberg angle [20] by a factor 30 to 6 × 10 −6 -an uncertainty dominated by the absolute calibration of the beam energy -thus reducing the uncertainty on α QED to 5 × 10 −6 .

Theoretical uncertainties
Theoretical uncertainties on the muon forward-backward asymmetry arise from the lack of higher orders in the calculation of the muon angular distribution. The dominant higherorder effects on the muon angles originate from QED corrections, namely (i) initial-state radiation (ISR), i.e., the emission of one or several photons by the incoming beams; (ii) final state radiation (FSR), i.e., the emission of photons from the outgoing muons; and (iii) the interference between ISR and FSR (IFI), which becomes significant when the muons are produced at small angles with respect to the beam axis. The effect of these three QED corrections on the muon angular distributions and on the muon forward-backward asymmetry have been studied analytically in Ref. [22] with a complete O(α) calculation and soft-photon exponentiation, and in a more pragmatic way by the OPAL experiment in Ref. [23]. Their conclusions are summarized and the relevant effects of the O(α 2 ) corrections are examined in this section. Other electroweak corrections are discussed at the end of the section.

Final-state radiation (FSR)
Final-state radiation is mostly collinear and is symmetric around the muon directions, at all orders in α. The effect on the forward-backward asymmetry is therefore expected to be unmeasurably small [23]. It was checked in Ref. [22] that the effect is rigorously 0 at order α (with soft-photon exponentiation) if no cut is applied on the final-state photon energy, and vanishingly small if a tight upper cut is applied to the final-state photon energy, typically of the order of α π Ecut √ s , i.e., ∼ 2 × 10 −6 for E cut ∼ 100 MeV. The theoretical uncertainty on this effect due to the O(α 2 ) corrections is typically one-to-two orders of magnitude smaller than that. It is therefore ignored in the following.

Initial-state radiation (ISR)
Initial-state radiation corrections are known up to order O(α 2 ) with soft-photon exponentiation [24]. Unlike FSR, ISR has a macroscopic influence on the forward-backward asymmetry. Photons from ISR are emitted mostly along the beam axis, with a twofold consequence: (i) the centre-of-mass frame of the muon pair therefore acquires a longitudinal boost, which modifies the angular distribution of both muons in a non trivial way; and (ii) the effective centre-of-mass energy of the collision is reduced to s and x ± = E γ ± / √ s are the fractional radiated energies by the e ± beams, which modifies the asymmetric term of the cross section through A µµ FB (s ). As A µµ FB (s ) varies quite fast with √ s , as displayed in Fig. 3 and expressed in Eq. 4.15, a large, negative, variation of the measured asymmetry is indeed to be expected.
When only one ISR photon is radiated by one of the two beams, the effects can be largely mitigated. In the vast majority of the cases, the photon is radiated exactly along the beam axis. The polar angles of the outgoing muons, denoted θ ± , suffice in that case to determine the effective centre-of-mass energy √ s : as well as the µ + direction in the centre-of-mass frame of the muon pair In this simplest configuration, the use of cos θ * entirely corrects for the effect of the longitudinal boost on the angular distribution, and the forward-backward asymmetry dependence on s /s can be studied explicitly. Furthermore, the events relevant for the determination of α ± can be selected by requiring s /s to be close to unity. If the initial-state photon is radiated with a finite angle with respect to the beam axis, however, Eqs. 4.19 and 4.20 no longer hold, but the corresponding events can be rejected by requiring the two muons to be back-to-back in the plane transverse to the beam axis. In rare cases, both beams can radiate photons, which render these two equations only approximate, and may still create a bias in the forward-backward asymmetry. To determine the effect of this approximation, large samples of µ + µ − events were generated at √ s ± . The simulation of ISR was performed with the REMT package [25] modified to include O(α 2 ) correcctions with soft-photon exponentiation, and the possibility to radiate up to two photons. For the reasons just explained, only events with s /s in excess of 0.999 and an acoplanarity angle between the two muons smaller than 0.35 mrad were considered. These two cuts typically select about 80% of the cross section in Fig. 2, and tremendously increase the purity towards events without ISR. The blue histograms in Figs. 7 show, for √ s = √ s − and √ s + , the relative biasses on A µµ FB (s ± ) with respect to the standard model prediction, as a function of 1 − s /s and for a perfect muon angular resolution, σ θ = σ φ = 0.
Events with only one ISR photon would lead to a blue straight line at ∆A FB /A FB ≡ 0.0, as s /s can be exactly determined in that case from Eq. 4.19. The possibility to radiate photons from the two beams, however, induces a visible systematic effect on the measured asymmetry, up to 0.1%, as soon as the energy of these photons is non-zero (i.e., in all but the first bin of the blue histogram). The sign of this effect can be understood as follows: the emission of two photons in opposite directions reduces the effective boost of the µ + µ − pair, causing the value of s , as determined from Eq. 4.19, to be larger than the true effective centre-of-mass energy. The corresponding event migration towards the left of the blue histogram therefore tends to reduce the value of the forward-backward asymmetry (hence a relative increase below the Z peak and a relative decrease above). The effect in the first bin is much smaller than in the other bins because it contains the vast majority of the events, hence is little affected by the migration from bins with a much smaller number of events. A non-perfect muon angular resolution also affects the determination of s /s from Eq. 4.19 for all events (with or without ISR), hence has a non-trivial effect on the measured asymmetry, as shown in the red histograms of Figs. 7 for a uniform angular resolution σ θ = σ φ = 0.1 mrad. Because a worse angular resolution would lead to a proportionally larger systematic bias, the measurement of α QED can be used as a benchmark to define the muon reconstruction performance of the FCC-ee detectors. (As an example, the track angular resolution of the pixel detector of SiD [26], with a hit resolution of 5 µm, a thickness of 1.5% of a radiation length, a radius of 6 cm, and a length of 34 cm is better than 0.1 mrad over the whole acceptance.) For most of the events, without significant ISR photons, the measured angle between the two muons can only decrease from its maximal value of π rad, yielding a smaller value of s from Eq. 4.19. Part of these events therefore migrate to the other bins of the distribution, and therefore increase significantly the value of the measured forward-backward asymmetry in these bins, by up to 0.4% (with a relative variation negative below the Z peak, and positive above). The residual migration from events with initial state radiation, but with two wrongly measured back-to-back muons, has the opposite effect in the first bin of the distribution, albeit with a much smaller amplitude, because of the much larger number of events originally in this bin.
These two systematic biasses are much larger, by two orders of magnitude, than the target precision with which the forward-backward asymmetry needs to be measured. These biasses can be corrected for if (i) the energy and angular distributions of initial-state radiation can be predicted with an accuracy better than 1%, which is probably the case already today; and if (ii) the muon angular resolution can be mapped with a precision of a few per mil over the whole detector acceptance, which is probably feasible with the large samples of K 0 s 's, Λ's and even J/ψ's expected at the FCC-ee. The predicted relative biasses, however, appear to be quasi-"universal", in the sense that they are similar in amplitude below and above the Z peak, albeit in opposite directions. The combination of the two measurements towards a determination of α QED (m 2 Z ) with Eq. 4.9 exhibits an almost perfect cancellation in all bins, as displayed in Fig. 8 as a function of 1 − s /s, with the same vertical scale as in Figs. 7. When integrated over all bins, the total relative bias on α QED (m 2 Z ) amounts to −8 × 10 −6 , i.e., well below the target statistical precision of a few 3 × 10 −5 . The aforementioned theoretical knowledge of initial-state radiation and the in-situ determination of the angular resolution would allow this residual bias to be predicted and corrected for, with a precision at least an order of magnitude better. 10 × Relative bias and statistical uncertainty on the electromagnetic coupling constant estimated at the Z mass scale, as a function 1 − s /s , from measurements of the muon forwardbackward asymmetry at √ s ± , with a perfect muon angular resolution (blue histogram) and with σ θ = σ φ = 0.1 mrad (red histogram).

Interference between initial-and final-state radiation (IFI)
While initial-state radiation does not change the functional form of the muon angular distribution, the interference between initial-state and final-state occurs preferably when the final state muons are close to the initial state electrons, hence does affect their distribution in the forward and backward directions beyond the usual (1 + cos 2 θ * ) + 8/3A FB cos θ * formula. It is shown in Ref. [23] that the muon angular distribution is be modified by a multiplicative factor with a characteristic logarithmic dependence on cos θ * , and can be parameterized as in presence of a tight muon acoplanarity cut as suggested in the previous section. The multiplicative factor contains an additional asymmetric term, which enhances the integrated muon forward-backward asymmetry. The tight cut on s /s aimed at rejecting ISR also reduces IFI in similar proportions. To mitigate the very small residual effect of IFI on the angular distribution, the specific shape of the additional contribution can be fitted away, as was done at LEP and with the benefit of the much larger data samples expected at the FCC-ee. On the other hand, this additional contribution appears to be "universal", (i.e., with an amplitude that depends only on s /s, similarly to what is observed for ISR), hence cancels out in the determination of α QED (m 2 Z ) from a combination of the measurements at the two centre-of-mass energies, with no loss of statistical power.

Other electroweak higher-order corrections
Other electroweak corrections have so far "only" been computed off-peak with complete one-loop calculation [27]. One-loop box corrections lead to relative changes of −9 × 10 −4 at √ s − and −8 × 10 −4 at √ s + from the improved Born approximation of A µµ FB − A µµ FB,0 , hence to a shift of α QED (m Z ) at the per-mil level. A shift of similar size arises from one-loop vertex corrections. The theoretical uncertainty arising from the missing higher orders in the asymmetry calculation, estimated to be at the level of a few 10 −4 [28], was adequate at the time of LEP but is insufficient today to match the precision offered by the FCC-ee.
An order of magnitude improvement would be achievable today, with proven techniques, by including the dominant two-loop and leading three-loop corrections, and would represent a major breakthrough towards the FCC-ee targets. Meeting these targets might require a complete three-loop calculation, including three-loop box corrections, perhaps a serious challenge with the current techniques, and definitely beyond the scope of the present work. It is not unlikely, however, that a large part of these missing corrections affect in the same way the asymmetry at 87.9 GeV and the asymmetry at 94.3 GeV. If it were the case, the α QED (m 2 Z ) determination would enjoy a cancellation similar to the that observed for QED corrections, which could suffice even without a complete three-loop calculation.

Conclusions and outlook
In this paper, it has been shown that the measurement of the muon forward-backward asymmetry at the FCC-ee, with six months of data taking just below ( √ s = 87.9 GeV) and just above ( √ s = 94.3 GeV) the Z peak, as part of the Z resonance scan, would open the opportunity of a direct measurement of the electromagnetic constant α QED (m 2 Z ), with a relative statistical uncertainty of the order of 3 × 10 −5 .
A comprehensive list of sources for experimental, parametric, theoretical systematic uncertainties has been examined. Most of these uncertainties have been shown to be under control at the level of 10 −5 or below, as summarized in Table 1. A significant fraction of those benefits from a delicate cancellation between the two asymmetry measurements. The knowledge of the beam energy, both on-and off-peak, turns out to be the dominant contribution, albeit still well below the targeted statistical power of the method.

Type
Source Uncertainty  Table 1. Summary of relative statistical, experimental, parametric and theoretical uncertainties to the direct determination of the electromagnetic coupling constant at the FCC-ee, with a one-year running period equally shared between centre-of-mass energies of 87.9 and 94.3 GeV, corresponding to an integrated luminosity of 85 ab −1 .
The fantastic integrated luminosity and the unique beam-energy determination are the key breakthrough advantages of the FCC-ee in the perspective of a precise determination of the electromagnetic coupling constant. Today, the only obstacle towards this measurement -beside the construction of the collider and the delivery of the target luminosities -stems from the lack of higher orders in the determination of the electroweak corrections to the forward-backward asymmetry prediction in the standard model. With the full one-loop calculation presently available for these corrections, a relative uncertainty on A µµ FB of the order of a few 10 −4 is estimated. An improvement deemed adequate to match the FCC-ee experimental precision might require a calculation beyond two loops, which may be beyond the current state of the art, but is possibly within reach on the time scale required by the FCC-ee.
A consistent international programme for present and future young theorists must therefore be set up towards significant precision improvements in the prediction of all electroweak precision observables, in order to reap the rewards potentially offered by the FCCee.

Interference between Initial and Final state radiation (IFI)
In Section 4.3.3 of Ref. [1], it is stated that the tight cut on s /s aimed at rejecting ISR also reduces IFI in similar proportions. This statement is incomplete. While rejecting indeed events with energetic photons in similar proportions with or without IFI, a tight cut on s /s has no effect on soft photons, which in turn have a macroscopic effect (several per-cent increase) on the numerical value of the muon forward-backward asymmetry when IFI is included. As pointed out in Ref. [3], an important cancellation, already exploited in Ref. [1] for ISR, still occurs for IFI between the two centre-of-mass energies √ s − and √ s + .
The precision with which this cancellation can be predicted is being evaluated and is the topic of two forthcoming publications [4,5].