research28 min studyLesson 919 in Physics

Hadronic neutrinos require targets, thresholds and decay

Compare pp and photohadronic kinematics, interaction efficiency and the detector’s energy-dependent response.

On this page

What you’ll understand

  • Distinguish fixed-target pp and photon-target pion thresholds.
  • Separate kinematic allowance from pion-production and decay efficiency.
  • Relate an emitted neutrino population to flavor propagation and effective-area detection.

Accelerating protons is only the first requirement for a hadronic neutrino source. The protons need suitable gas or photon targets, enough interaction time, and charged secondaries that decay before losing their energy. Photons and neutrinos emerging from the same hadronic cascade can then experience very different propagation. We will distinguish exact threshold kinematics, broad energy-partition estimates and detector response instead of assigning every detected neutrino one fixed parent-proton energy.

Gas targets and photon targets have different thresholds

For a proton striking a stationary proton, the lowest neutral-pion channel is pp→ppπ⁰. Its projectile kinetic threshold is about 280 MeV. Charged-pion neutrino channels have different final-state masses and slightly different thresholds; a neutral-pion threshold by itself is not a neutrino threshold. In photon-target production, evaluate the proton and photon energies in the same frame. A head-on ultrarelativistic proton colliding with a photon reaches a specified nucleon–pion final state when their energy product supplies its invariant mass.

The often-used delta-resonance approximation describes where a photohadronic cross section is enhanced, rather than the first allowed pion threshold. Broad photon spectra and collision angles require an interaction-rate integral; a monochromatic head-on estimate cannot replace it. Throughout, primes indicate a local plasma frame. Observed neutrino energies acquire the factor δ/(1+z) from that frame. A target photon energy measured by the observer must first be transformed before combining it with a comoving proton energy.

Tp,thrppπ0=2mπ0c2+mπ02c22mp,Ep′ϵ′≳[(mN+mπ)2−mp2]c44(head-on, Ep′≫mpc2)T_{p,\rm thr}^{pp\pi^0}=2m_{\pi^0}c^2+\frac{m_{\pi^0}^2c^2}{2m_p},\qquad E_p^\prime\epsilon^\prime\gtrsim\frac{[(m_N+m_\pi)^2-m_p^2]c^4}{4}\quad(\text{head-on},\ E_p^\prime\gg m_pc^2)
The pp target is stationary in the frame defining the projectile kinetic energy Tp. For pγ the final-state nucleon N and pion specify the channel; the displayed product is the ultrarelativistic head-on limit of invariant kinematics. Energies may be in GeV consistently, making the product GeV².↔ Scroll wide equations horizontally.

Production efficiency and decay cooling are independent tests

At high proton energy, an illustrative gas-target energy-loss time is the reciprocal of density times inelastic cross section times speed times inelasticity. With constant coefficients over a residence interval, the fractional proton energy lost is one minus an exponential. This energy-loss fraction is not the interaction probability, because inelasticity is smaller than one. Only a channel-dependent fraction of that lost energy enters charged pions. Their neutrino yield also depends on whether pions and muons decay before synchrotron, adiabatic or other losses modify them.

A rough characteristic neutrino energy of five percent of the parent proton energy assumes a representative pion energy partition and an unsuppressed decay chain. It is broad, energy dependent and inappropriate as an exact threshold relation. Neutral pion photons can pair-produce and cascade, while neutrinos usually escape the source with much less interaction. Consequently a gamma-ray upper limit is not automatically the same upper limit on emitted neutrino power. Conversely a leptonic gamma-ray source need not produce a corresponding hadronic neutrino flux.

tpp−1≃κppnpσppc,fpp=1−e−tres/tpp,Eν,obs∼δ1+z 0.05Ep′t_{pp}^{-1}\simeq\kappa_{pp}n_p\sigma_{pp}c,\quad f_{pp}=1-e^{-t_{\rm res}/t_{pp}},\qquad E_{\nu,\rm obs}\sim\frac{\delta}{1+z}\,0.05E_p^\prime
np is gas proton density in m⁻³, σpp an effective inelastic cross section in m², κpp a dimensionless mean inelasticity and tres,tpp in s, all in the same target frame. The final energy estimate applies well above threshold to a representative pion-decay yield, not each event.↔ Scroll wide equations horizontally.

Propagation and selection enter before a source is inferred

The conventional charged-pion decay chain gives approximate source flavor numbers one electron to two muon neutrinos when particles and antiparticles are combined and muons decay without major cooling. Over astrophysical baselines, averaged vacuum mixing gives a composition near one to one to one, not an exact universal ratio. Muon damping, neutron decay or other source conditions change it. Neutrinos travel nearly straight, but detection probability depends strongly on energy and direction, and Earth absorption becomes relevant at sufficiently high energy. A measured track energy is not a unique incoming-neutrino energy.

For a point-source number fluence per energy, the expected selected event count is its energy integral weighted by the direction-dependent effective area, summed over flavors. That response includes interaction, propagation through Earth and the chosen event cuts; it is not merely geometric detector area. The Markov lab illustrates only repeated removal: Reset and set both controls p=0.1,q=0; Step twice for two transitions, predicting surviving state-zero probability 0.81. Reset and set p=0.2,q=0; Step twice, predicting 0.64. There is no energy-dependent hadronic cascade, flavor oscillation or neutrino detector in this two-state analogy.

Compare two explicitly different pion thresholds

Use mₚc²=0.938 GeV, mₙc²=0.940 GeV, mπ⁰c²=0.135 GeV and mπ⁺c²=0.140 GeV. Find the fixed-target pp→ppπ⁰ kinetic threshold, then the head-on pγ→nπ⁺ proton energy for comoving ε′=1 keV. Use the ultrarelativistic product approximation for the latter.

  1. The pp neutral-pion threshold is 2×0.135+0.135²/(2×0.938)=0.270+0.00971≈0.2797 GeV of projectile kinetic energy.

  2. For the charged photopion channel, the product constant is [(0.940+0.140)²−0.938²]/4=(1.1664−0.879844)/4=0.071639 GeV².

  3. Since 1 keV=10⁻⁶ GeV, E′p,thr≈0.071639/10⁻⁶=7.16×10⁴ GeV=71.6 TeV. This is a charged-channel threshold in the specified frame, not a characteristic neutrino energy or a guaranteed efficient interaction.

The result

The selected pp neutral-pion threshold is about 280 MeV kinetic energy; the selected head-on charged photopion threshold is about 71.6 TeV for 1 keV targets.

Explore the model

Interactive lab

Interactive lab · markov

Check repeated removal without claiming a cascade

Reset and set BOTH controls p=0.1,q=0; Step twice to two transitions: predict P(state 0)=0.9²=0.81. Reset and set p=0.2,q=0; Step twice: predict 0.8²=0.64. State 1 is absorbing. This two-state exact probability model includes no hadronic energies, pion decays or detector response.

Check repeated removal without claiming a cascadeExact finite-step probabilities for a two-state chain starting in state 0, not sampled trajectories. Playback advances two transitions per second. When p = q = 0 every distribution is stationary; when p = q = 1 the chain alternates forever. A single zero transition creates an absorbing state. Current readings: P(state 1) at step 0: 0%; Stationary P(state 1): 60%; Chain behavior: Monotone convergence. All parameters are adjustable with the labeled controls below.Two-state Markov distribution after 0 steps: state zero probability 1, state one probability 0.Initial distribution (1, 0) · step 0p = 0.3q = 0.201100%0%1−p1−qP(0)P(1)Monotone convergence
P(state 0)P(state 1)Stationary distribution
P(state 1) at step 0
0%
Stationary P(state 1)
60%
Chain behavior
Monotone convergence
t = 0 s

Exact finite-step probabilities for a two-state chain starting in state 0, not sampled trajectories. Playback advances two transitions per second. When p = q = 0 every distribution is stationary; when p = q = 1 the chain alternates forever. A single zero transition creates an absorbing state. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forA removal probability must be specified separately from particle-production kinematics and energy partition.

Check your understanding

Why is the estimate Eν≈0.05Ep unsuitable as an exact event-by-event inference?

Try it yourself

For a hypothetical high-energy gas target take nₚ=10¹² m⁻³, σpp=4×10⁻³⁰ m², κpp=0.5 and tᵣₑₛ=10⁸ s. Use c=3×10⁸ m/s and an injected proton energy reservoir of 10⁴⁰ J. Suppose two thirds of lost energy enters charged pions and three quarters of charged-pion energy becomes all-flavor neutrinos, with no secondary cooling. Estimate fpp and neutrino energy.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

First calculate the energy-loss exponent κppnₚσppctᵣₑₛ. Keep the two stated energy fractions separate from the loss fraction.

Compare with the solution

The inverse energy-loss time is 0.5×10¹²×4×10⁻³⁰×3×10⁸=6×10⁻¹⁰ s⁻¹, so the exponent is 0.06 and fpp=1−e⁻⁰·⁰⁶≈0.0582. The assumed charged-pion and decay fractions multiply to (2/3)(3/4)=1/2. Thus the all-flavor neutrino energy is about 0.5×0.0582×10⁴⁰=2.91×10³⁸ J in this target frame. Different channel fractions or secondary cooling change that result.

Take it with you

  • Specify the final-state channel and frame before calculating a pion threshold.
  • Target optical depth, energy loss and decay efficiency are different factors.
  • Flavor propagation and effective detector area are needed to relate emission to selected events.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.