Neutrino oscillation constraints on U(1)′ models: from non-standard interactions to long-range forces

We quantify the effect of gauge bosons from a weakly coupled lepton flavor dependent U(1)′ interaction on the matter background in the evolution of solar, atmospheric, reactor and long-baseline accelerator neutrinos in the global analysis of oscillation data. The analysis is performed for interaction lengths ranging from the Sun-Earth distance to effective contact neutrino interactions. We survey ∼ 10000 set of models characterized by the six relevant fermion U(1)′ charges and find that in all cases, constraints on the coupling and mass of the Z′ can be derived. We also find that about 5% of the U(1)′ model charges lead to a viable LMA-D solution but this is only possible in the contact interaction limit. We explicitly quantify the constraints for a variety of models including $$ \mathrm{U}{(1)}_{B-3{L}_e} $$ U 1 B − 3 L e , $$ \mathrm{U}{(1)}_{B-3{L}_{\mu }} $$ U 1 B − 3 L μ , $$ \mathrm{U}{(1)}_{B-3{L}_{\tau }} $$ U 1 B − 3 L τ , $$ \mathrm{U}{(1)}_{B-\frac{3}{2}\left({L}_{\mu }+{L}_{\tau}\right)} $$ U 1 B − 3 2 L μ + L τ , $$ \mathrm{U}{(1)}_{L_e-{L}_{\mu }} $$ U 1 L e − L μ , $$ \mathrm{U}{(1)}_{L_e-{L}_{\tau }} $$ U 1 L e − L τ , $$ \mathrm{U}{(1)}_{L_e-\frac{1}{2}\left({L}_{\mu }+{L}_{\tau}\right)} $$ U 1 L e − 1 2 L μ + L τ . We compare the constraints imposed by our oscillation analysis with the strongest bounds from fifth force searches, violation of equivalence principle as well as bounds from scattering experiments and white dwarf cooling. Our results show that generically, the oscillation analysis improves over the existing bounds from gravity tests for Z′ lighter than ∼ 10−8→ 10−11 eV depending on the specific couplings. In the contact interaction limit, we find that for most models listed above there are values of g′ and MZ′ for which the oscillation analysis provides constraints beyond those imposed by laboratory experiments. Finally we illustrate the range of Z′ and couplings leading to a viable LMA-D solution for two sets of models.