Precision $e^+e^-$ Hemisphere Masses in the Dijet Region with Power Corrections

We derive high-precision results for the $e^+e^-$ heavy jet mass (HJM) $d \sigma/d \rho$ and dihemisphere mass (DHM) $d^2\sigma/(d s_1 d s_2)$ distributions, for $s_1\sim s_2$, in the dijet region. New results include: i) the N$^3$LL resummation for HJM of large logarithms $\ln^n(\rho)$ at small $\rho$ including the exact two-loop non-global hemisphere soft function, the 4-loop cusp anomalous dimension and the 3-loop hard and jet functions, ii) N$^3$LL results for DHM with resummation of logarithms $\ln(s_{1,2}/Q^2)$ when there is no large separation between $s_1$ and $s_2$, iii) profile functions for HJM to give results simultaneously valid in the peak and tail regions, iv) a complete two-dimensional basis of non-perturbative functions which can be used for double differential observables, that are needed for both HJM and DHM in the peak region, and v) an implementation of renormalon subtractions for large-angle soft radiation to ${\cal O}(\alpha_s^3)$ together with a resummation of the additional large $\ln(Q\rho/\Lambda_{QCD})$ logarithms. Here $Q$ is the $e^+e^-$ center-of-mass energy. Our resummation results are combined with known fixed-order ${\cal O}(\alpha_s^3)$ results and we discuss the convergence and remaining perturbative uncertainty in the cross section. We also prove that, at order $1/Q$, the first moment of the HJM distribution involves an additional non-perturbative parameter compared to the power correction that shifts the tail of the spectrum (where $1\gg \rho\gg \Lambda_{QCD}/Q$). This differs from thrust where a single non-perturbative parameter at order $1/Q$ describes both the first moment and the tail, and it disfavors models of power corrections employing a single non-perturbative parameter, such as the low-scale effective coupling model. In this paper we focus only on the dijet region, not the far-tail distribution for $\rho \gtrsim 0.2$.