The recent observation of possible 11 neutron-antineutron ( \(n \text- \overline{n }\) ) oscillation candidates with an expected background of 9.3 ± 2.7 events at Super-Kamiokande has renewed the interest in ∆B = 2 transitions. In this work, we analyze the Renormalization Group (RG) running of mass dimension-9 six-quark operators, in \(\overline{M }S\) scheme, that generate processes like nn → π0π0, deuteron decay, \(n \text- \overline{n }\) oscillations etc, evolving them from the electroweak scale to baryon number violation (BNV) scale ( \(\mathcal{O}\) (103 TeV)). Our goal is to systematically account for the influence of potential new physics at intermediate energies (≳ \(\mathcal{O}\) (10 TeV)), especially given the fact that Large Hadron Collider has not ruled out new physics beyond ~ 10 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}\) (103 TeV), and (ii) an extended framework involving scalar and vector bosons above ~ 10 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. The code and instructions can be downloaded from https://github.com/rp-winter/bnv-running.
