Next-to-leading-order QCD matching for ∆F = 2 processes in scalar leptoquark models

Leptoquarks provide viable solutions to the flavour anomalies, i.e. they can explain the tensions between the measurements and the Standard Model predictions of the anomalous magnetic moment of the muon as well as b → sℓ$^{+}$ℓ$^{−}$ and b → cτν processes. However, LQs also contribute to other flavour observables, such as ∆F = 2 processes, at the loop-level. In particular, B$_{s}$−$ \overline{B} _{s}$ mixing provides a crucial bound in setups addressing b → cτν data, often excluding a big portion of the parameter space that could otherwise account for it. In this article, we first derive the complete leading order matching, including all five scalar LQ representations, for D$^{0}$−$ \overline{D} ^{0}$, K$^{0}$−$ \overline{K} ^{0}$ and B$_{s,d}$−$ \overline{B} _{s,d}$ mixing (at the dimension-six level). We then calculate the next-to-leading order α$_{s}$ matching corrections to these ∆F = 2 processes in generic scalar leptoquark models. We find that the two-loop corrections increase the effects in ∆F = 2 processes by ~ 5-10% and significantly reduce the matching scale uncertainty.