RG evolution and effect of intermediate new physics on $ΔB=2$ six-quark operators

Abstract: The recent identification of possible 11 neutron-antineutron ($n$-$\bar{n}$) oscillation candidate events at Super-Kamiokande has renewed the interest in $\Delta B = 2$ transitions. In this work, we analyze the Renormalization Group (RG) running of mass dimension-9 six-quark operators, in $\bar{MS}$ scheme, that generate processes like $nn\to \pi^0\pi^0$, deuteron decay, $n$-$\bar{n}$ oscillations etc, evolving them from the electroweak scale to baryon number violating scale ($\mathcal{O}(10^3~\text{TeV})$). Our goal is to systematically account for the influence of potential new physics at intermediate energies ($\gtrsim \mathcal{O}(10~ \text{TeV})$), especially given the fact that {\it Large Hadron Collider} has not ruled out new physics beyond $\sim 10~\text{TeV}$. To comprehensively investigate their influence, we consider two scenarios: (i) a minimal setup with only Standard Model degrees of freedom up to the high scale at $\mathcal{O}(10^3~\text{TeV})$, and (ii) an extended framework involving scalar and vector bosons above $\sim 10~\text{TeV}$ up till BNV scale. To facilitate further studies, we also provide a Python script that performs RG evolution of the BNV Wilson coefficients in the presence of generic bosonic new physics at any intermediate energy scale. It can be modified easily to meet the needs of the user to investigate the running of the BNV Wilson coefficients. We then compare the result with the experimental bound from the neutron-antineutron oscillation process and constrain the scale of baryon number violating new physics.