A light DM model for large B → K + invisible and K → π + invisible decays and its implications for Bₛ-B̄ₛ mixing and neutron EDM

We study the implications for Bₛ - B̄ₛ B s - B ¯ s mixing and the neutron electric dipole moment (EDM) in a light dark matter model with sizable invisible rare meson decays to accommodate the recent possible deviations from Standard Model (SM) predictions observed in B⁺→ K⁺ν ν ̄ B + → K + ν ν ¯ by Belle II and K⁺→ π ⁺ν ν ̄ K + → π + ν ν ¯ by NA62. Given that the neutrinos in these decays escape detection, they can be replaced by other invisible final states. Based on effective operator analysis, it has been proposed that branching ratios for B⁺→ K⁺ + invisible B + → K + + invisible and K⁺→ π ⁺ + invisible K + → π + + invisible can naturally be larger than the SM predictions due to the emission of light dark matter pairs. We demonstrate that this scenario can be realized within a UV-complete two-Higgs-doublet model (2HDM) where neutral Higgs bosons mediating dark matter interactions induce significant low-energy effects especially for Bₛ-B̄ₛ B s - B ¯ s mixing and neutron EDM. Within the allowed parameter space, we find non-negligible contributions to Bₛ - B̄ₛ B s - B ¯ s mixing. For the neutron EDM, there is a cancellation due to the exchange of neutral spin-zero particles, but QCD renormalization group evolution will lift this cancellation which in fact is generally true for any neutral Higgs contribution. However, we demonstrate that such a cancellation does not occur for charged scalar contributions. Ultimately, the allowed CP-violating phases in the Yukawa sector can generate a neutron EDM at a level consistent with the current experimental bound.