Subleading Logarithmic QED Initial State Corrections to $e^+e^- \rightarrow \gamma^*/{Z^{0}}^*$ to $O(\alpha^6 L^5)$

Using the method of massive operator matrix elements, we calculate the subleading QED initial state radiative corrections to the process $e^+e^- \rightarrow \gamma^*/Z^*$ for the first three logarithmic contributions from $O(\alpha^3 L^3), O(\alpha^3 L^2), O(\alpha^3 L)$ to $O(\alpha^5 L^5), O(\alpha^5 L^4), O(\alpha^5 L^3)$ and compare their effects to the leading contribution $O(\alpha^6 L^6)$ and one more subleading term $O(\alpha^6 L^5)$. The calculation is performed in the limit of large center of mass energies squared $s \gg m_e^2$. These terms supplement the known corrections to $O(\alpha^2)$, which were completed recently. Given the high precision at future colliders operating at very large luminosity, these corrections are important for concise theoretical predictions. The present calculation needs the calculation of one more two--loop massive operator matrix element in QED. The radiators are obtained as solutions of the associated Callen--Symanzik equations in the massive case. The radiators can be expressed in terms of harmonic polylogarithms to weight {\sf w = 6} of argument $z$ and $(1-z)$ and in Mellin $N$ space by generalized harmonic sums. Numerical results are presented on the position of the $Z$ peak and corrections to the $Z$ width, $\Gamma_Z$. The corrections calculated result into a final theoretical accuracy for $\delta M_Z$ and $\delta \Gamma_Z$ which is estimated to be of O(30 keV) at an anticipated systematic accuracy at the FCC\_ee of \sim 100 keV. This precision cannot be reached, however, by including only the corrections up to $O(\alpha^3)$.