Neutral-current NSIs: Beyond Standard Oscillations

• Neutral-current NSIs are effective four-fermion operators that modify neutrino propagation in matter through new vector interactions.
• They are probed through complementary methods including oscillation experiments, CEνNS scattering, and collider missing-energy searches.
• Light mediator models and loop-induced effects offer viable frameworks, influencing both neutrino oscillation degeneracies and direct-detection signatures.

Searching arXiv for recent and foundational papers on neutral-current NSIs to support the encyclopedia entry. Neutral-current non-standard interactions (NSIs) are effective beyond-Standard-Model neutrino interactions that modify neutral-current amplitudes involving neutrinos and matter fermions. In the conventional oscillation framework, they are encoded as four-fermion operators of the form $$\mathcal{L}_{\rm NSI} = -2 \sqrt{2} G_F \varepsilon_{\alpha\beta}^{fP} \left( \overline{\nu_\alpha} \gamma^\mu P_L \nu_\beta \right) \left( \overline{f} \gamma_\mu P_C f \right)$$, with $f=e,u,d$, and they affect neutrino propagation through coherent forward scattering in matter (Ohlsson, 2012). In this standard sense, only the vector combination is relevant for matter effects in unpolarized media, and the effective matter Hamiltonian takes the form $$H_{\rm mat}=\sqrt{2} G_F N_e \left( \begin{matrix} 1+\epsilon_{ee} & \epsilon_{e\mu} & \epsilon_{e\tau} \cr \epsilon_{e\mu}^* & \epsilon_{\mu\mu} & \epsilon_{\mu\tau}\cr \epsilon_{e\tau}^*& \epsilon_{\mu\tau}^*&\epsilon_{\tau\tau}\end{matrix}\right)$$ (Bakhti et al., 2016). The same quark-level couplings also appear in neutral-current scattering, including coherent elastic neutrino–nucleus scattering (CE$\nu$NS), deep-inelastic scattering, and collider missing-energy signatures, so neutral-current NSIs sit at the interface of oscillation phenomenology, neutrino scattering, and new-physics model building (Dev et al., 2019).

1. Effective operator framework and matter-potential formulation

The conventional neutral-current NSI operator is written as

$$\mathcal{L}_{\mathrm{NC-NSI}} = -2\sqrt{2}G_F \varepsilon_{\alpha\beta}^{fC} \bigl(\overline{\nu_\alpha}\gamma^\mu P_L \nu_\beta\bigr) \bigl(\overline{f}\gamma_\mu P_C f\bigr),$$

with $\alpha,\beta=e,\mu,\tau$, $f=e,u,d$, and chirality $C=L,R$ (Chatterjee et al., 2020). Hermiticity requires $$\varepsilon_{\beta\alpha}^{fC}=(\varepsilon_{\alpha\beta}^{fC})^*$$ (Agarwalla et al., 2016). For oscillations in matter, the relevant effective quantities are matter-composition-weighted combinations,

$$\varepsilon_{\alpha\beta} = \sum_{f,P}\varepsilon_{\alpha\beta}^{fP}\,\frac{n_f}{n_e},$$

or equivalently $f=e,u,d$0 for approximately neutral, isoscalar Earth matter (0907.0097). In the same approximation, $f=e,u,d$1, $f=e,u,d$2, and $f=e,u,d$3 (Chatterjee et al., 2020).

The propagation Hamiltonian is the sum of the vacuum term and the matter term,

$f=e,u,d$4

with $f=e,u,d$5 (Chatterjee et al., 2020). Only differences of diagonal NSI parameters are physical in oscillations, because adding a term proportional to the identity does not affect flavor transitions; accordingly, analyses commonly set $f=e,u,d$6 and interpret $f=e,u,d$7 and $f=e,u,d$8 relative to it (Bakhti et al., 2016). Off-diagonal entries such as $f=e,u,d$9, $$H_{\rm mat}=\sqrt{2} G_F N_e \left( \begin{matrix} 1+\epsilon_{ee} & \epsilon_{e\mu} & \epsilon_{e\tau} \cr \epsilon_{e\mu}^* & \epsilon_{\mu\mu} & \epsilon_{\mu\tau}\cr \epsilon_{e\tau}^*& \epsilon_{\mu\tau}^*&\epsilon_{\tau\tau}\end{matrix}\right)$$0, and $$H_{\rm mat}=\sqrt{2} G_F N_e \left( \begin{matrix} 1+\epsilon_{ee} & \epsilon_{e\mu} & \epsilon_{e\tau} \cr \epsilon_{e\mu}^* & \epsilon_{\mu\mu} & \epsilon_{\mu\tau}\cr \epsilon_{e\tau}^*& \epsilon_{\mu\tau}^*&\epsilon_{\tau\tau}\end{matrix}\right)$$1 are generally complex and introduce new phases beyond the standard Dirac phase (Ohlsson, 2012).

A recurring point in the literature is that the phenomenologically relevant NC NSI coefficients are environment-dependent effective combinations, not fundamental couplings by themselves. The same microscopic theory therefore maps differently onto Earth-like and solar-like matter (Ohlsson, 2012). The review literature also stresses that the operators written in this form are low-energy effective interactions; in the heavy-mediator limit they are understood as arising after integrating out new degrees of freedom, whereas light-mediator realizations require separate treatment (Dev et al., 2019).

2. Flavor structure, diagonal ambiguities, and standard bounds

The flavor structure of neutral-current NSIs divides naturally into diagonal non-universal terms and off-diagonal flavor-changing terms. Diagonal entries $$H_{\rm mat}=\sqrt{2} G_F N_e \left( \begin{matrix} 1+\epsilon_{ee} & \epsilon_{e\mu} & \epsilon_{e\tau} \cr \epsilon_{e\mu}^* & \epsilon_{\mu\mu} & \epsilon_{\mu\tau}\cr \epsilon_{e\tau}^*& \epsilon_{\mu\tau}^*&\epsilon_{\tau\tau}\end{matrix}\right)$$2, $$H_{\rm mat}=\sqrt{2} G_F N_e \left( \begin{matrix} 1+\epsilon_{ee} & \epsilon_{e\mu} & \epsilon_{e\tau} \cr \epsilon_{e\mu}^* & \epsilon_{\mu\mu} & \epsilon_{\mu\tau}\cr \epsilon_{e\tau}^*& \epsilon_{\mu\tau}^*&\epsilon_{\tau\tau}\end{matrix}\right)$$3, and $$H_{\rm mat}=\sqrt{2} G_F N_e \left( \begin{matrix} 1+\epsilon_{ee} & \epsilon_{e\mu} & \epsilon_{e\tau} \cr \epsilon_{e\mu}^* & \epsilon_{\mu\mu} & \epsilon_{\mu\tau}\cr \epsilon_{e\tau}^*& \epsilon_{\mu\tau}^*&\epsilon_{\tau\tau}\end{matrix}\right)$$4 modify flavor-dependent forward scattering, while off-diagonal entries $$H_{\rm mat}=\sqrt{2} G_F N_e \left( \begin{matrix} 1+\epsilon_{ee} & \epsilon_{e\mu} & \epsilon_{e\tau} \cr \epsilon_{e\mu}^* & \epsilon_{\mu\mu} & \epsilon_{\mu\tau}\cr \epsilon_{e\tau}^*& \epsilon_{\mu\tau}^*&\epsilon_{\tau\tau}\end{matrix}\right)$$5, $$H_{\rm mat}=\sqrt{2} G_F N_e \left( \begin{matrix} 1+\epsilon_{ee} & \epsilon_{e\mu} & \epsilon_{e\tau} \cr \epsilon_{e\mu}^* & \epsilon_{\mu\mu} & \epsilon_{\mu\tau}\cr \epsilon_{e\tau}^*& \epsilon_{\mu\tau}^*&\epsilon_{\tau\tau}\end{matrix}\right)$$6, and $$H_{\rm mat}=\sqrt{2} G_F N_e \left( \begin{matrix} 1+\epsilon_{ee} & \epsilon_{e\mu} & \epsilon_{e\tau} \cr \epsilon_{e\mu}^* & \epsilon_{\mu\mu} & \epsilon_{\mu\tau}\cr \epsilon_{e\tau}^*& \epsilon_{\mu\tau}^*&\epsilon_{\tau\tau}\end{matrix}\right)$$7 induce flavor-changing matter effects (0907.0097). Because only diagonal differences are observable, bounds on diagonal matter NSIs are intrinsically less direct than bounds on off-diagonal ones (Ohlsson, 2012).

A useful compiled benchmark for effective matter NSIs in Earth-like matter is

$$H_{\rm mat}=\sqrt{2} G_F N_e \left( \begin{matrix} 1+\epsilon_{ee} & \epsilon_{e\mu} & \epsilon_{e\tau} \cr \epsilon_{e\mu}^* & \epsilon_{\mu\mu} & \epsilon_{\mu\tau}\cr \epsilon_{e\tau}^*& \epsilon_{\mu\tau}^*&\epsilon_{\tau\tau}\end{matrix}\right)$$8

and for solar-like matter

$$H_{\rm mat}=\sqrt{2} G_F N_e \left( \begin{matrix} 1+\epsilon_{ee} & \epsilon_{e\mu} & \epsilon_{e\tau} \cr \epsilon_{e\mu}^* & \epsilon_{\mu\mu} & \epsilon_{\mu\tau}\cr \epsilon_{e\tau}^*& \epsilon_{\mu\tau}^*&\epsilon_{\tau\tau}\end{matrix}\right)$$9

with the explicit caveat that these are approximate effective combinations built from bounds obtained assuming one nonzero NSI coupling at a time (0907.0097). Interpreted entry by entry, the strongest bounds apply to $\nu$0-sector parameters such as $\nu$1, while $\nu$2 and especially $\nu$3 can remain comparatively large (0907.0097).

The same hierarchy is echoed in the broader review literature. Robustly constrained matter NSIs are mainly those involving $\nu$4-flavor, especially $\nu$5, and to a lesser extent $\nu$6 and $\nu$7, whereas weakly constrained matter NSIs remain in the $\nu$8- and especially $\nu$9-diagonal sectors (Ohlsson, 2012). A further atmospheric relation,

$$\mathcal{L}_{\mathrm{NC-NSI}} = -2\sqrt{2}G_F \varepsilon_{\alpha\beta}^{fC} \bigl(\overline{\nu_\alpha}\gamma^\mu P_L \nu_\beta\bigr) \bigl(\overline{f}\gamma_\mu P_C f\bigr),$$0

illustrates that matter-NSI bounds are often shaped by correlations rather than by independent one-parameter exclusions (0907.0097).

A common misconception is to treat “NSI” as automatically meaning the usual vector matter operator. One study explicitly broadens the term to the full set of dimension-6 neutrino–quark four-fermion operators built from scalar, pseudo-scalar, vector, axial-vector, and tensor bilinears, but also notes that only the vector-current interaction corresponds most directly to the usual neutral-current NSI of oscillation and CE$$\mathcal{L}_{\mathrm{NC-NSI}} = -2\sqrt{2}G_F \varepsilon_{\alpha\beta}^{fC} \bigl(\overline{\nu_\alpha}\gamma^\mu P_L \nu_\beta\bigr) \bigl(\overline{f}\gamma_\mu P_C f\bigr),$$1NS phenomenology (Chao et al., 2019). This distinction matters because the scalar, axial-vector, pseudo-scalar, and tensor interactions belong to a broader class of exotic neutrino interactions even when they are discussed under the NSI label (Chao et al., 2019).

3. Propagation effects in oscillation experiments

Long-baseline, atmospheric, and solar oscillation experiments probe neutral-current NSIs through propagation in matter. The qualitative reason is that NSIs modify the flavor structure of the matter potential and therefore alter the interference pattern between vacuum oscillations and matter effects (Dev et al., 2019). In long-baseline accelerator appearance channels, the key parameters are usually the $$\mathcal{L}_{\mathrm{NC-NSI}} = -2\sqrt{2}G_F \varepsilon_{\alpha\beta}^{fC} \bigl(\overline{\nu_\alpha}\gamma^\mu P_L \nu_\beta\bigr) \bigl(\overline{f}\gamma_\mu P_C f\bigr),$$2-flavor off-diagonal entries $$\mathcal{L}_{\mathrm{NC-NSI}} = -2\sqrt{2}G_F \varepsilon_{\alpha\beta}^{fC} \bigl(\overline{\nu_\alpha}\gamma^\mu P_L \nu_\beta\bigr) \bigl(\overline{f}\gamma_\mu P_C f\bigr),$$3 and $$\mathcal{L}_{\mathrm{NC-NSI}} = -2\sqrt{2}G_F \varepsilon_{\alpha\beta}^{fC} \bigl(\overline{\nu_\alpha}\gamma^\mu P_L \nu_\beta\bigr) \bigl(\overline{f}\gamma_\mu P_C f\bigr),$$4, because the $$\mathcal{L}_{\mathrm{NC-NSI}} = -2\sqrt{2}G_F \varepsilon_{\alpha\beta}^{fC} \bigl(\overline{\nu_\alpha}\gamma^\mu P_L \nu_\beta\bigr) \bigl(\overline{f}\gamma_\mu P_C f\bigr),$$5 channel is directly sensitive to them (Chatterjee et al., 2020).

For $$\mathcal{L}_{\mathrm{NC-NSI}} = -2\sqrt{2}G_F \varepsilon_{\alpha\beta}^{fC} \bigl(\overline{\nu_\alpha}\gamma^\mu P_L \nu_\beta\bigr) \bigl(\overline{f}\gamma_\mu P_C f\bigr),$$6, the appearance probability can be written approximately as

$$\mathcal{L}_{\mathrm{NC-NSI}} = -2\sqrt{2}G_F \varepsilon_{\alpha\beta}^{fC} \bigl(\overline{\nu_\alpha}\gamma^\mu P_L \nu_\beta\bigr) \bigl(\overline{f}\gamma_\mu P_C f\bigr),$$7

with a standard atmospheric term $$\mathcal{L}_{\mathrm{NC-NSI}} = -2\sqrt{2}G_F \varepsilon_{\alpha\beta}^{fC} \bigl(\overline{\nu_\alpha}\gamma^\mu P_L \nu_\beta\bigr) \bigl(\overline{f}\gamma_\mu P_C f\bigr),$$8, a standard interference term $$\mathcal{L}_{\mathrm{NC-NSI}} = -2\sqrt{2}G_F \varepsilon_{\alpha\beta}^{fC} \bigl(\overline{\nu_\alpha}\gamma^\mu P_L \nu_\beta\bigr) \bigl(\overline{f}\gamma_\mu P_C f\bigr),$$9, and an NSI interference term

$\alpha,\beta=e,\mu,\tau$0

where $\alpha,\beta=e,\mu,\tau$1 (Chatterjee et al., 2020). Because $\alpha,\beta=e,\mu,\tau$2, matter-sensitive experiments such as NO$\alpha,\beta=e,\mu,\tau$3A and DUNE are substantially more affected than shorter-baseline experiments such as T2K (Chatterjee et al., 2020). This is the basis of the NO$\alpha,\beta=e,\mu,\tau$4A–T2K NSI explanation: with $\alpha,\beta=e,\mu,\tau$5 or $\alpha,\beta=e,\mu,\tau$6 and phases near $\alpha,\beta=e,\mu,\tau$7, the discrepancy in the inferred $\alpha,\beta=e,\mu,\tau$8 between the two experiments can be alleviated in normal ordering (Chatterjee et al., 2020).

The same mechanism produces degeneracies in other parameters. At DUNE, the appearance probability acquires an additional NSI-induced interference term carrying the phase $\alpha,\beta=e,\mu,\tau$9, and for $f=e,u,d$0 or $f=e,u,d$1, the discovery potential of the octant of $f=e,u,d$2 can be completely lost for unfavorable combinations of $f=e,u,d$3 and $f=e,u,d$4 (Agarwalla et al., 2016). At the DUNE first oscillation maximum, the standard octant-sensitive splitting $f=e,u,d$5, the standard interference amplitude $f=e,u,d$6, and the NSI interference amplitude $f=e,u,d$7 are numerically comparable when $f=e,u,d$8, which is the quantitative reason the octant sensitivity can collapse (Agarwalla et al., 2016).

Shorter-baseline, lower-energy experiments can instead act as degeneracy breakers. For MOMENT, the inequalities

$f=e,u,d$9

make both standard and non-standard matter effects small, so the determination of $C=L,R$0 is comparatively robust against propagation NSIs (Bakhti et al., 2016). Combining MOMENT with T2K and NO$C=L,R$1A removes wrong-$C=L,R$2 solutions and improves some NSI bounds, including $C=L,R$3 at $C=L,R$4 in the combined analysis (Bakhti et al., 2016).

The general review literature consistently emphasizes that NC NSIs can mimic or obscure CP violation, weaken mass-ordering and octant sensitivity, and generate generalized degeneracies in long-baseline fits (Dev et al., 2019). A plausible implication is that robust extraction of standard oscillation parameters requires external information from scattering, atmospheric, and matter-insensitive measurements.

4. Neutral-current scattering, CE$C=L,R$5NS, and direct-detection implications

Neutral-current NSIs also modify scattering cross sections. In CE$C=L,R$6NS, the differential cross section is

$C=L,R$7

with the weak charge

$C=L,R$8

for flavor $C=L,R$9 (Khan et al., 2021). This expression makes clear that flavor-diagonal amplitudes interfere with the Standard Model, whereas flavor-changing detection amplitudes contribute incoherently in the summed final state (Khan et al., 2021).

A particularly specific result concerns CP-violating phases in CE$$\varepsilon_{\beta\alpha}^{fC}=(\varepsilon_{\alpha\beta}^{fC})^*$$0NS. Observable phase dependence requires simultaneously nonzero $$\varepsilon_{\beta\alpha}^{fC}=(\varepsilon_{\alpha\beta}^{fC})^*$$1 and $$\varepsilon_{\beta\alpha}^{fC}=(\varepsilon_{\alpha\beta}^{fC})^*$$2 for the same $$\varepsilon_{\beta\alpha}^{fC}=(\varepsilon_{\alpha\beta}^{fC})^*$$3, because the rate depends only on the relative phase

$$\varepsilon_{\beta\alpha}^{fC}=(\varepsilon_{\alpha\beta}^{fC})^*$$4

through an interference term proportional to $$\varepsilon_{\beta\alpha}^{fC}=(\varepsilon_{\alpha\beta}^{fC})^*$$5 (Khan et al., 2021). A single isolated off-diagonal NSI parameter has no observable phase dependence in CE$$\varepsilon_{\beta\alpha}^{fC}=(\varepsilon_{\alpha\beta}^{fC})^*$$6NS (Khan et al., 2021). This sharply contrasts with oscillation propagation, where a single complex off-diagonal matter NSI can already affect the inferred CP structure.

Neutral-current scattering also extends beyond CE$$\varepsilon_{\beta\alpha}^{fC}=(\varepsilon_{\alpha\beta}^{fC})^*$$7NS. A DUNE-like experiment using NC deep-inelastic scattering can probe axial NC NSI of neutrinos with $$\varepsilon_{\beta\alpha}^{fC}=(\varepsilon_{\alpha\beta}^{fC})^*$$8, $$\varepsilon_{\beta\alpha}^{fC}=(\varepsilon_{\alpha\beta}^{fC})^*$$9, and $$\varepsilon_{\alpha\beta} = \sum_{f,P}\varepsilon_{\alpha\beta}^{fP}\,\frac{n_f}{n_e},$$0 quarks, which can hide from oscillations and CE$$\varepsilon_{\alpha\beta} = \sum_{f,P}\varepsilon_{\alpha\beta}^{fP}\,\frac{n_f}{n_e},$$1NS because those standard probes are dominantly sensitive to vector couplings (Abbaslu et al., 2023). Using both near and far detectors, such a DUNE-like setup can significantly improve present bounds on all the flavor elements studied there, and for the first time the couplings to strange quarks become accessible in a direct way (Abbaslu et al., 2023).

NSIs also propagate into dark-matter direct-detection phenomenology through the neutrino background. In the effective-operator treatment of neutrino–quark interactions, vector-current and scalar-current NSIs can significantly change the neutrino floor, pseudo-scalar NSIs can raise it by about $$\varepsilon_{\alpha\beta} = \sum_{f,P}\varepsilon_{\alpha\beta}^{fP}\,\frac{n_f}{n_e},$$2, while axial-vector and tensor NSIs alter it only at the $$\varepsilon_{\alpha\beta} = \sum_{f,P}\varepsilon_{\alpha\beta}^{fP}\,\frac{n_f}{n_e},$$3 level (Chao et al., 2019). The reason is structural: vector NSIs use the coherent $$\varepsilon_{\alpha\beta} = \sum_{f,P}\varepsilon_{\alpha\beta}^{fP}\,\frac{n_f}{n_e},$$4 response and interfere with the Standard Model, scalar NSIs are coherently enhanced but non-interfering, and axial/tensor effects depend on spin-sensitive responses that are small or vanish for some targets (Chao et al., 2019). This has led to the explicit suggestion that dark-matter direct-detection experiments should be combined with CE$$\varepsilon_{\alpha\beta} = \sum_{f,P}\varepsilon_{\alpha\beta}^{fP}\,\frac{n_f}{n_e},$$5NS measurements to test whether a recoil excess is due to exotic neutral-current neutrino interactions rather than WIMP scattering (Chao et al., 2019).

5. Collider probes and long-baseline neutral-current event samples

Collider missing-transverse-momentum searches provide a complementary high-energy probe of neutrino–quark NC NSIs. In a standard NC NSI parameterization,

$$\varepsilon_{\alpha\beta} = \sum_{f,P}\varepsilon_{\alpha\beta}^{fP}\,\frac{n_f}{n_e},$$6

the LHC is effectively flavor blind because final-state neutrino flavor is not identified (Liu et al., 2020). Recasting 13 TeV ATLAS and CMS mono-object data with NLO QCD matched to parton showering yields a final EFT-style collider bound $$\varepsilon_{\alpha\beta} = \sum_{f,P}\varepsilon_{\alpha\beta}^{fP}\,\frac{n_f}{n_e},$$7 at $$\varepsilon_{\alpha\beta} = \sum_{f,P}\varepsilon_{\alpha\beta}^{fP}\,\frac{n_f}{n_e},$$8 CL for $$\varepsilon_{\alpha\beta} = \sum_{f,P}\varepsilon_{\alpha\beta}^{fP}\,\frac{n_f}{n_e},$$9, while in a simplified $f=e,u,d$00 model with $f=e,u,d$01 the combined bounds are $f=e,u,d$02 for $f=e,u,d$03 TeV and $f=e,u,d$04 for $f=e,u,d$05 TeV (Liu et al., 2020). The paper emphasizes that the collider limits are especially useful for heavy mediators or contact-like interactions and that the EFT approximation requires the new-physics scale to exceed the typical LHC partonic energy (Liu et al., 2020).

Within oscillation experiments themselves, neutral-current event samples have emerged as a direct NSI probe. One recent line of work shows that NC event rates at long-baseline experiments become oscillation sensitive once NSIs make NC cross sections flavor dependent (Gehrlein et al., 2024). Existing NOvA NC data then yield one-parameter-at-a-time 90% C.L. bounds such as

$f=e,u,d$06

for axial-vector NSIs under the isospin-conserving assumption, together with competitive vector constraints and a disfavoring of large diagonal vector NSI associated with the LMA-Dark region (Gehrlein et al., 2024). The paper emphasizes that these NC long-baseline data are highly complementary to SNO because SNO is not sensitive to the isospin-conserving axial direction $f=e,u,d$07 in the same way (Gehrlein et al., 2024).

An even more specific complementarity has been developed in terms of isoscalar and isovector quark combinations. Charged-current oscillation analyses constrain NSIs through matter effects in propagation and probe almost exclusively the isoscalar combination of up- and down-quark couplings, while the orthogonal isovector combination is suppressed by a factor of $f=e,u,d$08 in Earth matter (Gehrlein et al., 17 Apr 2026). NC event samples, by contrast, depend on modified NC cross sections and the flavor composition of the beam at near and far detectors, and therefore probe both isoscalar and isovector couplings with comparable weight (Gehrlein et al., 17 Apr 2026). Using NOvA data and DUNE projections, this framework yields the first bounded long-baseline constraints on isovector NSI and resolves the individual quark couplings only when CC and NC information are combined (Gehrlein et al., 17 Apr 2026).

6. UV completions, loop generation, and theoretical caveats

Large neutral-current NSIs are nontrivial to realize in realistic models. In the low-energy EFT, one expects $f=e,u,d$09, so a mediator scale $f=e,u,d$10 suggests $f=e,u,d$11, whereas $f=e,u,d$12 suggests $f=e,u,d$13 (Ohlsson, 2012). This scaling already indicates that order-one matter NSIs are difficult in heavy-mediator theories unless special structures are present.

The most developed successful constructions therefore exploit light mediators. One class of viable models introduces a $f=e,u,d$14 gauge boson $f=e,u,d$15 with mass $f=e,u,d$16 and tiny couplings $f=e,u,d$17, so that forward scattering in matter can be large while finite-momentum-transfer scattering bounds are suppressed (Farzan, 2016). In one model aimed at the LMA-Dark pattern, the generated quark NSIs satisfy

$f=e,u,d$18

with $f=e,u,d$19 and negligible off-diagonal entries (Farzan, 2016). A second model generates diagonal and off-diagonal NSI through mixing of active neutrinos with a new Dirac fermion $f=e,u,d$20, leading to

$f=e,u,d$21

and allowing complex off-diagonal matter NSIs relevant for CP-violating propagation effects (Farzan, 2016).

A related light-$f=e,u,d$22 framework focused on lepton-flavor-violating matter NSIs shows that sizable propagation NSIs are realistically possible only in the $f=e,u,d$23-$f=e,u,d$24 sector, because $f=e,u,d$25 and $f=e,u,d$26 cases are excluded by charged-lepton flavor-violating limits such as $f=e,u,d$27 and $f=e,u,d$28 (Farzan et al., 2015). In the viable $f=e,u,d$29-$f=e,u,d$30 case, the model predicts correlated propagation parameters

$f=e,u,d$31

with

$f=e,u,d$32

so that long-baseline data can in principle determine both the overall strength and the flavor angle (Farzan et al., 2015).

Loop generation provides an additional route. A one-loop extension of the NSI formalism shows that effective matter NSIs arise from loop-corrected $f=e,u,d$33 couplings plus dimension-six four-fermion operators,

$f=e,u,d$34

with gauge invariance enforcing finiteness through cancellations among wave-function and vertex corrections (Bellazzini et al., 2010). In the MSSM example developed there, the largest loop-induced effects found numerically are mostly source/detector rather than matter NSIs, but the formalism makes clear how genuine neutral-current matter NSIs emerge at one loop (Bellazzini et al., 2010).

Several caveats recur across the literature. The heavy-mediator EFT description can fail for light mediators or large momentum transfer (Liu et al., 2020). Oscillation matter effects constrain vector combinations, not axial ones (Abbaslu et al., 2023). Many quoted bounds assume one nonzero NSI parameter at a time, and correlated multiparameter fits can be substantially weaker (0907.0097). A plausible implication is that any single quoted interval must be interpreted together with its flavor, mediator, and operator assumptions.

7. Open issues, complementarity, and current status

The current status of neutral-current NSIs is defined less by a single bound than by a network of complementary probes. Oscillation experiments constrain effective matter combinations and are especially sensitive to propagation phases and degeneracies (Dev et al., 2019). CE$f=e,u,d$35NS constrains the same quark-level vector couplings through scattering and is directly relevant to oscillation-motivated scenarios such as the NO$f=e,u,d$36A–T2K tension and the LMA-Dark region (Khan et al., 2021). Collider missing-energy searches constrain heavy-mediator neutrino–quark NSIs in a way that is flavor blind but sensitive to absolute coupling magnitudes (Liu et al., 2020). Long-baseline NC event samples now probe both vector and axial structures, including directions that are nearly invisible to CC-only oscillation fits (Gehrlein et al., 2024).

The literature repeatedly emphasizes complementarity rather than replacement. One status report frames this broadly in terms of oscillation data, scattering data, colliders, and model-building constraints all being necessary to map the parameter space (Dev et al., 2019). More recent long-baseline work sharpens that point by showing that CC analyses probe almost only the isoscalar quark direction through matter effects, whereas NC event analyses access both isoscalar and isovector directions through cross sections (Gehrlein et al., 17 Apr 2026). This suggests that “neutral-current NSI” is best understood not as a single observable effect, but as a family of operator deformations whose phenomenology depends strongly on whether one studies propagation, scattering, or both.

Several controversies or common misunderstandings are now clear. Neutral-current NSI in the standard oscillation sense means the vector neutrino–matter interaction entering the propagation Hamiltonian, not the entire space of scalar, pseudo-scalar, axial, or tensor neutrino operators (Chao et al., 2019). Large diagonal effective NSIs are not automatically excluded in a model-independent sense, although they are difficult to embed in gauge-invariant UV completions (Ohlsson, 2012). Apparent NSI-like effects in oscillation experiments can also arise from source and detector corrections or non-unitarity, so identifying a propagation NSI requires correlating oscillation and scattering information (Bellazzini et al., 2010).

The most conservative synthesis is that neutral-current NSIs remain one of the standard EFT descriptions of sub-leading new physics in neutrino flavor transitions, but their interpretation is highly assumption dependent. Vector matter NSIs are tightly constrained in many directions yet still phenomenologically relevant because of degeneracies in long-baseline fits and because diagonal sectors remain comparatively weak (Dev et al., 2019). Axial neutral-current NSIs can evade the standard oscillation and CE$f=e,u,d$37NS probes and therefore motivate dedicated NC scattering analyses at DUNE-like facilities (Abbaslu et al., 2023). Light-mediator constructions show that large observable matter NSIs are not automatically ruled out, but they require specific flavor structures and are tied to broader laboratory, astrophysical, and cosmological signatures (Farzan, 2016).