Dihadron Angular Correlations in the $e^+e^-$ Collision

Abstract: The precision of fixed-order calculations on the dihadron production in electron-positron annihilation is paramount for probing QCD factorization and constraining non-perturbative inputs. This paper investigates the QCD corrections to the angular separation distribution $\theta_{12}$ between two observed hadrons, $H_1$ and $H_2$, in the process $e^+e^- \to H_1 H_2 + X$ up to $\mathcal{O}(\alpha_s^2)$, with particular emphasis on the intermediate region $\theta_{12} \in (0,\pi)$. The partonic processes at this accuracy consist of two sorts of contributions, the real-virtual and double-real corrections. Of them, the evaluation of four-body phase space integrals in the latter case is at the core of this study. To address them, we first employ the integration-by-parts (IBP) identities to reduce the number of independent integrals and then apply the differential equations (DE) method to recursively solve the resulting master integrals. In kinematic regions where the invariant mass of the unresolved partons vanishes, IBP coefficients can develop divergences. To this end, we resum higher-order terms in the dimensional regulator for each master integral based on the asymptotic behavior of the canonical DEs. After combining the real and virtual corrections with the counter terms from fragmentation function renormalization, we demonstrate that the pole terms in the final analytic expressions exactly cancel out in all partonic channels, thereby providing a non-trivial validation of collinear factorization at the next-to-leading order (NLO). Eventually, when presenting our analytic expressions of the finite partonic coefficients, we transform the transcendental functions resulting from the DE solutions into classical (poly)logarithmic functions, in order to facilitate the implementation in event generators.