Unified framework for B-anomalies, muon g-2 and neutrino masses

We present a model of radiative neutrino masses which also resolves anomalies reported in B-meson decays, R-D(*),R- and R-D(*),R-, as well as in muon g- 2 measurement, Delta a(mu). Neutrino masses arise in the model through loop diagrams involving TeV-scale leptoquark (LQ) scalars R-2 and S-3. Fits to neutrino oscillation parameters are obtained satisfying all flavor constraints which also explain the anomalies in R-D(*),R-, RRD(*), and Delta a(mu) within 1 sigma. An isospin-3/2 Higgs quadruplet plays a crucial role in generating neutrino masses; we point out that the doubly-charged scalar contained therein can be produced in the decays of the S3 LQ, which enhances its reach to 1.1 (6.2) TeV at p root s = 14TeV high-luminosity LHC (root s = 100TeV FCC-hh). We also present flavor-dependent upper limits on the Yukawa couplings of the LQs to the first two family fermions, arising from non-resonant dilepton (pp -> l(+)l(-)) processes mediated by t-channel LQ exchange, which for 1TeV LQ mass, are found to be in the range (0.15 - 0.36). These limits preclude any explanation of R-D(*) through LQ-mediated B-meson decays involving nu(e) or nu(mu) in the final state. We also find that the same Yukawa couplings responsible for the chirally-enhanced contribution to Delta a(mu) give rise to new contributions to the SM Higgs decays to muon and tau pairs, with the modifications to the corresponding branching ratios being at (2-6)% level, which could be tested at future hadron colliders, such as HL-LHC and FCC-hh.