The excess in electron recoil events reported recently by the XENON1T experiment may be interpreted as evidence for a sizable transition magnetic moment $\mu_{\nu_e\nu_\mu}$ of Majorana neutrinos. We show the consistency of this scenario when a single component transition magnetic moment takes values $\mu_{\nu_e\nu_\mu} \in(1.65 - 3.42) \times 10^{-11} \mu_B$. Such a large value typically leads to unacceptably large neutrino masses. In this paper we show that new leptonic symmetries can solve this problem and demonstrate this with several examples. We first revive and then propose a simplified model based on $SU(2)_H$ horizontal symmetry. Owing to the difference in their Lorentz structures, in the $SU(2)_H$ symmetric limit, $m_\nu$ vanishes while $\mu_{\nu_e\nu_\mu}$ is nonzero. Our simplified model is based on an approximate $SU(2)_H$, which we also generalize to a three family $SU(3)_H$-symmetry. Collider and low energy tests of these models are analyzed. We have also analyzed implications of the XENON1T data for the Zee model and its extensions which naturally generate a large $\mu_{\nu_e\nu_\mu}$ with suppressed $m_\nu$ via a spin symmetry mechanism, but found that the induced $\mu_{\nu_e\nu_\mu}$ is not large enough to explain recent data. Finally, we suggest a mechanism to evade stringent astrophysical limits on neutrino magnetic moments arising from stellar evolution by inducing a medium-dependent mass for the neutrino.

We present a simple extension of the Standard Model that leads to renormalizable long-range vector-mediated neutrino self-interactions. This model can resolve the Hubble tension by delaying the onset of neutrino free-streaming during recombination, without conflicting with other measurements. The extended gauge, scalar and neutrino sectors lead to observable signatures, including invisible Higgs and $Z$ decays, thereby relating the Hubble tension to precision measurements at the LHC and future colliders. The model has a new neutrinophilic gauge boson with $m_{Z'}\sim\mathcal{O}(10~\mathrm{eV})$ and charged Higgses at a few $100~\mathrm{GeV}$. It requires hidden neutrinos with active-hidden mixing angles larger than $5\times10^{-4}$ and masses in the range $1\div300\mathrm{eV}$, which could also play a role for short baseline neutrino oscillation anomalies.

Ultra light particles $(m_a \sim 10^{-21}eV-10^{-22}eV)$ with axion-like couplings to other particles can be candidates for fuzzy dark matter (FDM) if the axion decay constant $f_a\sim 10^{17}GeV$. If a compact star is immersed in such a low mass axionic potential it develops a long range field outside the star. This axionic field is radiated away when the star is in a binary orbit. The orbital period of a compact binary decays mainly due to the gravitational wave radiation, which was confirmed first in the Hulse-Taylor binary pulsar. The orbital period can also decay by radiation of other light particles like axions and axion like particles(ALPs). For axionic radiation to take place, the orbital frequency of the periodic motion of the binary system should be greater than the mass of the scalar particle which can be radiated. This implies that, for most of the observed binaries, particles with mass $m_a< 10^{-19}eV$ can be radiated, which includes FDM particles. In this paper, we consider four compact binary systems: PSR J0348+0432, PSR J0737-3039, PSR J1738+0333, and PSR B1913+16 (Hulse Taylor Binary) and show that the observations of the decay in orbital period put the bound on axion decay constant, $f_a\lesssim \mathcal{O}(10^{11}GeV)$. This implies that Fuzzy Dark Matter cannot couple to gluons.

The standard model particles can be gauged in an anomaly free way by three possible gauge symmetries namely ${L_e-L_\mu}$, ${L_e-L_\tau}$, and ${L_\mu-L_\tau}$. Of these, ${L_e-L_\mu}$ and ${L_e-L_\tau}$ forces can mediate between the Sun and the planets and change planetary orbits. It is well known that a deviation from the $1/r^2$ Newtonian force can give rise to a perihelion advancement in the planetary orbit, for instance, as in the well known case of Einstein's gravity which was tested from the observation of the perihelion advancement of the Mercury. We consider the Yukawa potential of ${L_e-L_{\mu,\tau}}$ force which arises between the Sun and the planets if the mass of the gauge boson is $M_{Z^{\prime}}\leq \mathcal{O}(10^{-19})\rm {eV}$. We derive the formula for the perihelion advancement for such Yukawa type fifth force. We find that perihelion advancement is proportional to the square of the semi major axis of the orbit for the Yukawa potential, unlike GR, where it is largest for the nearest planet. We take the observational limits of all planets for which the perihelion advancement is measured and we obtain the gauge boson coupling $g$ in the range $10^{-18}$ to $10^{-16}$ for the mass range $10^{-22}$eV to $10^{-18}$eV. This mass range of gauge boson can be a possible candidate of fuzzy dark matter whose effect can therefore be observed in the precession measurement of the planetary orbits.

We propose a new way to probe non-standard interactions (NSI) of neutrinos with matter using the ultra-high energy (UHE) neutrino data at current and future neutrino telescopes. We consider the Zee model of radiative neutrino mass generation as a prototype, which allows two charged scalars -- one $SU(2)_L$-doublet and one singlet, both being leptophilic, to be as light as 100 GeV, thereby inducing potentially observable NSI with electrons. We show that these light charged Zee-scalars could give rise to a Glashow-like resonance feature in the UHE neutrino event spectrum at the IceCube neutrino observatory and its high-energy upgrade IceCube-Gen2, which can probe a sizable fraction of the allowed NSI parameter space.

The orbital period of a compact binary system decays mainly due to quadrupole gravitational radiation, which agrees with the observation to within one percent. Other types of radiation such as ultralight scalar or pseudoscalar radiation, massive vector boson radiation also contribute to the decay of orbital period as long as the mass of the emitted particle is less than the orbital frequency of the compact binary system. We obtain an expression of the energy loss due to the radiation of massive vector field from the neutron star-neutron star and neutron star-white dwarf binaries. Due to large chemical potential of the degenerate electrons, neutron stars have large muon charge. We derive the energy loss due to $U(1)_{L_\mu-L_\tau}$ gauge boson radiation from the binaries. For the radiation of vector boson, the mass is restricted by $M_{Z^\prime}<\Omega \simeq10^{-19}eV$ are the orbital frequencies of the compact star binaries. Using the formula of orbital period decay, we obtain constraints on the coupling constant of the gauge boson in the gauged $L_\mu-L_\tau$ theory for the four compact binary systems. For vector gauge boson muon coupling we find that for $M_{Z^\prime}<10^{-19}eV$, constraints on the coupling constant is $g<\mathcal{O}(10^{-20})$. We also obtain the exclusion plots of the massive vector proca field and the gauge field which can couple to muons.

This report summarizes the present status of neutrino non-standard interactions (NSI). After a brief overview, several aspects of NSIs are discussed, including connection to neutrino mass models, model-building and phenomenology of large NSI with both light and heavy mediators, NSI phenomenology in both short- and long-baseline neutrino oscillation experiments, neutrino cross-sections, complementarity of NSI with other low- and high-energy experiments, fits with neutrino oscillation and scattering data, DUNE sensitivity to NSI, effective field theory of NSI, as well as the relevance of NSI to dark matter and cosmology. We also discuss the open questions and interesting future directions that can be pursued by the community at large. This report is based on talks and discussions during the Neutrino Theory Network NSI workshop held at Washington University in St. Louis from May 29-31, 2019 (https://indico.cern.ch/event/812851/)

General relativity may be an effective theory with higher curvature correction terms in the action. Inclusion of these terms leads to exciting new possibilities, e.g., gravitational and electromagnetic perturbations following different geodesics; leading to a time delay. Such a time delay was observed between the gravitational wave event GW170817 and its electromagnetic counterpart GRB 170817A. We describe how this effect can be used to constrain the coupling of the higher curvature term. Our method is sufficiently general and applicable to any higher curvature theory.

A novel constraint on $f(R)$ theories of gravity is obtained from the gravitational wave signal emitted from the binary neutron star merger event GW170817. The $f(R)$ theories possess an additional massive scalar degree of freedom apart from the massless spin-2 modes. The corresponding scalar field contributes an additional attractive, short-ranged "fifth" force affecting the gravitational wave radiation process. We realize that chameleon screening is necessary to conform with the observation. A model independent bound $\vert f'(R_0) -1\vert < 3\times 10^{-3}$ has been obtained, where the prime denotes the derivative with respect $R$ and $R_0$ is the curvature of our Universe at present. Though we use the nonrelativistic approximations and obtain an order of magnitude estimate of the bound, it comes from direct observations of gravitational waves and thus it is worth noting. This bound is stronger/equivalent compared to some earlier other bounds such as from the Cassini mission in the Solar-System, Supernova monopole radiation, the observed CMB spectrum, galaxy cluster density profile, etc., although it is weaker than best current constraints ($\vert f'(R_0) -1\vert \lesssim 10^{-6}$) from cosmology. Using the bound obtained, we also constrain the parameter space in the $f(R)$ theories of dark energy like Hu-Sawicki, Starobinsky, and Tsujikawa models.

The observations of gravitational waves from the binary neutron star merger event GW170817 and the subsequent observation of its electromagnetic counterparts from the gamma-ray burst GRB 170817A provide us a significant opportunity to study theories of gravity beyond general relativity. An important outcome of these observations is that they constrain the difference between the speed of gravity and the speed of light to less than $10^{-15}c$. Also, the time delay between the arrivals of gravitational waves at different detectors constrains the speed of gravity at the Earth to be in the range $0.55c < v_{gw} < 1.42c$. We use these results to constrain a widely studied modified theory of gravity: Eddington-inspired Born-Infeld (EiBI) gravity. We show that, in EiBI theory, the speed of gravitational waves in matter deviates from $c$. From the time delay in arrival of gravitational wave signals at Earth-based detectors, we obtain the bound on the theory parameter $\kappa$ as $\vert\kappa\vert \lesssim 10^{21}\, m^2$. Similarly, from the time delay between the signals of GW170817 and GRB 170817A, in a background Friedmann-Robertson-Walker universe, we obtain $\vert \kappa \vert \lesssim 10^{37}\, m^2$. Although the bounds on $\kappa$ are weak compared to other earlier bounds from the study of neutron stars, stellar evolution, primordial nucleosynthesis, etc., our bounds are from the direct observations and thus worth noting.

Recently proposed Born-Infeld (BI) theories of gravity assume a constant BI parameter ($\kappa$). However, no clear consensus exists on the sign and value of $\kappa$. Recalling the Brans-Dicke (BD) approach, where a scalar field was used to generate the gravitational constant $G$, we suggest an extension of Born-Infeld gravity with a similar Brans-Dicke flavor. Thus, a new action, with $\kappa$ elevated to a spacetime dependent real scalar field, is proposed. We illustrate this new theory in a cosmological setting with pressureless dust and radiation as matter. Assuming a functional form of $\kappa(t)$, we numerically obtain the scale factor evolution and other details of the background cosmology. It is known that BI gravity differs from general relativity (GR) in the strong-field regime but reduces to GR for intermediate and weak fields. Our studies in cosmology demonstrate how, with this new, scalar-tensor BI gravity, deviations from GR as well as usual BI gravity, may arise in the weak-field regime too. For example, we note a late-time acceleration without any dark energy contribution. Apart from such qualitative differences, we note that fixing the sign and value of $\kappa$ is no longer a necessity in this theory, though the origin of the BD scalar does remain an open question.

Cosmology in Eddington-inspired Born-Infeld gravity is investigated using a scalar Born-Infeld field (e.g. tachyon condensate) as matter. In this way, both in the gravity and matter sectors we have Born-Infeld-like structures characterized by their actions and via two separate constants, $\kappa$ and $\alpha_T^2$ respectively. With a particular choice of the form of $\dot{\phi}$ (the time derivative of the Born-Infeld scalar), analytical cosmological solutions are found. Thereafter, we explore some of the unique features of the corresponding cosmological spacetimes. For $\kappa>0$, our solution has a de Sitter-like expansion both at early and late times, with an intermediate deceleration sandwiched between the accelerating phases. On the other hand, when $\kappa<0$, the initial de Sitter phase is replaced by a bounce. Our solutions, at late time, fit well with available supernova data-- a fact we demonstrate explicitly. The estimated properties of the Universe obtained from the fitting of the $\kappa>0$ solution, are as good as in $\Lambda$CDM cosmology. However, the $\kappa<0$ solution has to be discarded due to the occurrence of a bounce at an unacceptably low redshift.

A perfect fluid, spatially flat cosmology in a $f(T)$ model, derived from a recently proposed general Born-Infeld type theory of gravity is studied. Four dimensional cosmological solutions are obtained assuming the equation of state $p=\omega \rho$. For a positive value of $\lambda $ (a parameter in the theory) the solution is singular (of big-bang type) but may have accelerated expansion at an early stage. For $\lambda<0$ there exists a non-zero minimum scale factor and a finite maximum value of the energy density, but the curvature scalar diverges. Interestingly, for $\lambda <0$, the universe may undergo an eternal accelerated expansion with a de Sitter expansion phase at late times. We find these features without considering any extra matter field or even negative pressure. Fitting our model with Supernova data we find that the simplest dust model ($p=0$), with $\lambda >0$, is able to generate acceleration and fits well, although the resulting properties of the universe differ much from the known, present day, accepted values. The best fit model requires (with $\lambda > 0$) an additional component of the physical matter density, with a negative value of the equation of state parameter, along with dust. The $\lambda < 0$ solutions do not fit well with observations. Though these models do not explain the dark energy problem with consistency, their analysis does shed light on the plausibility of an alternative geometrical explanation.