Deep inelastic scattering Created 2026-10-05 Updated 2026-10-18
Scattering that resolves the constituents of a hadron at large spacelike momentum transfer, with an inclusive sum over final hadron states.
In Electron scattering from a proton, write the reaction as : the outgoing Electron is measured, while denotes all allowed final hadron states. Large momentum transfer supplies a short-distance probe; a large final invariant mass separates the deep continuum from elastic scattering and individual particle resonances. The useful surprise is that an extended proton can have a response characteristic of much smaller constituents. Deep inelastic scattering turns measurements outside the proton into constraints on its internal parton distribution functions.
The SLAC-MIT deep inelastic electron scattering experiment connects this definition to an actual measurement. Its evidence came from distributions of outgoing Electron energies and scattering angles, rather than a photograph of a quark. Bjorken scaling concerns deep-inelastic structure functions after known scattering factors have been accounted for; it does not say that the raw differential scattering cross-section is constant.
Work in natural units with , neglect the electron mass, and put the initial proton at rest with mass . Let and be the initial and final Electron four-momenta, , and . Define the laboratory energy transfer . This has units of energy; the invariant has units of energy squared. The Bjorken scaling variable and the final hadron invariant mass obey
These follow by squaring and , using and . Solving the Bjorken x definition gives the outgoing energy and then the scattering angle:
Alternatively, substituting before solving gives . A chosen must satisfy and . In addition, must exceed the lightest allowed inelastic hadron threshold; positive alone is insufficient.
For an original numerical example, take , , and . Then , , , and . Varying and at this same Bjorken x keeps every setting physical and gives , above the proton-pion threshold and the stricter chosen continuum selection . These are calculated settings, not historical measurements.
Let a projectile have mass , incident speed , and incident kinetic energy . The stationary target has mass . In a head-on elastic collision, the signed outgoing projectile velocity and target velocity obey momentum conservation and conservation of energy:
Factor the second equation and eliminate . The nontrivial collision branch has , hence
This is the maximum over scattering angles, not the energy of every recoil. To see the bound without assuming a head-on event, set the recoil momentum magnitude to and its laboratory angle to the incident direction to . Conservation of energy with outgoing projectile momentum gives
Thus for a nonzero recoil, and the endpoint occurs at .
For an original numerical example, stipulate and , with the proton mass. Approximate hydrogen, helium and nitrogen target masses by , and respectively. Their maximum fractions and recoil energies are
These stipulated values illustrate the approximation ; they are not measurements from 1932. With , the incident ratio gives , consistent with the nonrelativistic limit used here. A million electronvolts is an energy unit, whereas a mass expressed in requires the factor when inserted into an energy equation.
The hydrogen endpoint permits essentially complete transfer because the projectile and target masses match. For a very heavy target, , the fraction approaches , which is small: the target can take substantial momentum with little kinetic energy. The exact inequality follows from , which implies .
Figure 1.
Recoil energy transfer and the two-target constraint on an unseen projectile's mass
. Details.
Consider a nonrelativistic particle of mass undergoing an elastic collision with a stationary free target of mass , where . Assume the outgoing direction is uniformly distributed over solid angle in the center of mass frame; this assumption concerns the scattering angle distribution and need not hold for every target or energy. Let be the cosine of that frame's scattering angle. Transforming the outgoing velocity back to the laboratory gives
The uniform distribution of on makes uniform on . For , integration yields
For equal masses, take the limit : . The expected value of the fractional energy after a collision is instead . These are different averages: averaging a logarithm does not give the logarithm of an average.
For independent collisions with the same distribution, the logarithms of successive energy ratios add. Thus
Here is the geometric-mean energy, and is a characteristic collision count for reaching an energy scale . It is not an exact first-passage expected value. Individual histories fluctuate, and neither the time between collisions nor absorption is specified by this model.
In an original chosen example, set and . Approximate hydrogen by and carbon by . Then
These are illustrative model values, not measured collision counts from 1934. The energy scale marks where neglected thermal motion matters; the calculation does not assert continued cooling to zero energy.
For an ultrarelativistic Electron scattered elastically by a static, spherical, spin-zero atomic nucleus, a plane-wave Born approximation retains a single electromagnetic interaction. The Fourier transform of its Coulomb potential is proportional to . Relative to the same point charge, the scattering amplitude is therefore multiplied by and the differential scattering cross-section by :
Here is the incident total Electron energy, the scattering angle, a solid angle element, and the fine-structure constant. The cosine factor already accounts for the projectile's spin; a spinless projectile would give a different point-charge factor. The formula neglects nuclear recoil, electron rest-mass corrections, radiative corrections and Coulomb distortion in electron scattering. The latter can be substantial at large : electron waves then require treatment in the target's Coulomb potential, rather than undistorted plane waves.
Nonzero target spin can introduce magnetic and additional electric contributions, so division by this point-charge expression does not generally isolate one charge form factor for a proton or deuteron. Even for spin zero, finite energy resolution, recoil and Coulomb distortion in electron scattering can change a measured minimum. A zero of the uniform-sphere model is a clear illustration of amplitude cancellation, not a universal prediction of a zero experimental counting rate.
Assume a real, central, short-range potential and elastic quantum scattering. At large radius, normalize the stationary wavefunction as
where and is the scattering angle. Equal incoming and outgoing speeds give the differential scattering cross-section .
The partial-wave scattering expansion separates the spherical harmonics. The incident plane wave has expansion , with Spherical Bessel functions and Legendre polynomials . The regular radial solution outside the potential has asymptotic form . Its real scattering phase shift changes the ratio of outgoing to incoming waves to the Partial-wave S-matrix . Matching outgoing coefficients yields
Here labels the orbital angular momentum, with squared magnitude . The assumption of purely elastic scattering is what makes every have unit modulus.
Integrate the squared partial-wave expansion of a scattering amplitude over the solid angle. Legendre polynomial orthogonality eliminates terms with unequal angular-momentum indices, giving
At zero scattering angle, and . Consequently
This is the optical theorem. In the present elastic problem the total and elastic cross sections coincide. For an absorptive interaction the inclusive total still obeys the optical theorem, but need not equal the integral of the elastic differential cross section.
In the first Born approximation, let be the transferred wave number. The angular integral in the three-dimensional Fourier transform is , so
Writing and , with value one at zero, the delta-shell Born amplitude gives
For the scattering is approximately isotropic, with total cross section . For it has a forward lobe of angular width of order and oscillatory zeros where . At a fixed nonzero scattering angle, its envelope decreases as . The right-hand panel of the figure shows this angular dependence.
For real the first-order scattering amplitude is real, so inserting it literally in the optical theorem would predict zero total cross section, contrary to its nonzero squared modulus. The first Born amplitude alone does not satisfy the optical theorem at order . This is an order-of-approximation issue: the squared first-order amplitude is already second order, and the second Born term supplies the required second-order imaginary forward amplitude. At first order both sides vanish. A weak-shell regime is needed for the Born truncation itself to be accurate.
A movable magnetic spectrometer selected scattered Electrons by momentum and scattering angle. Magnetic fields dispersed momenta at the focal plane. With relativistic Electrons, the inferred momentum determines outgoing energy . Commissioning included steering the beam through the spectrometer at zero angle to check its optics and accepted solid angle. Counters sorted events into momentum and angle bins. The acceptance entered the conversion from counts to differential scattering cross-sections.[ref]
Figure 1. . Original schematic, not an engineering drawing. The magnetic spectrometer measures the outgoing Electron; the inclusive final hadron state is unobserved. Details.
Radiative corrections matter because photon emission can move events from an elastic scattering peak into lower outgoing-energy bins. The analyses subtracted the elastic radiative tail before unfolding inelastic radiation. Independent MIT and SLAC procedures provided a check. An uncorrected low- tail alone therefore cannot establish deep inelastic scattering.[ref]
Alongside the 1927 Davisson–Germer experiment, George Paget Thomson developed an independent transmission test. Alexander Reid first used celluloid; Thomson used thin metal films. A narrow Electron beam crossed a film about thick, reaching a photographic plate perpendicular to the beam. Rings appeared. Accelerating voltages of – set the de Broglie wavelengths; metal structures known from X-ray diffraction allowed quantitative comparison. A downstream magnetic field deflected ring-forming Electrons like those passing through holes, identifying the carriers as Electrons rather than X-rays. Celluloid's uncertain structure limited its quantitative interpretation.[ref]
The Davisson–Germer experiment uses crystal reflection and angular electric current scans; this transmission camera records a spatial pattern behind the specimen. Both probe interference in electron diffraction, with different specimen preparations and readouts. A developed photographic plate records accumulated detections, not a direct image of the wavefunction.
Why rings? In an ideal film containing many small crystals with uniformly distributed orientations, a family of lattice planes with spacing can satisfy Bragg's law at every azimuth about the incident beam. First-order elastic scattering then forms a cone of opening half-angle , where . A flat detector a distance downstream cuts this cone in a circle:
The last expression requires small scattering angles. It explains a useful control: hold the specimen and camera distance fixed, increase accelerating voltage, and the ring radii should shrink with the electron accelerating-voltage wavelength. It predicts positions, not intensities: preferred crystal orientations, film thickness and repeated quantum scattering can alter a real pattern.
For a chosen example, take and . The exact relativistic electron accelerating-voltage wavelength gives at and at . The corresponding exact ring radii are and . These are calculated examples with chosen spacings and camera distance, not Thomson's measurements. The higher voltage lies outside the range quoted from his lecture, deliberately illustrating a prediction rather than reconstructing his data.
Figure 1.
Calculated electron diffraction rings shrink as accelerating voltage increases
. Original transmission geometry and calculated first-order ring positions for chosen plane spacings and at camera distance . Both detector views have identical scales. Exact relativistic de Broglie wavelengths and determine the radii; ring brightness and width are drawing choices. The geometry sketch exaggerates the scattering angles. No experimental data or borrowed diffraction photograph is plotted. Source. Details.
Consider one chosen model at a fixed total energy, with two open scattering channels, no absorption and no omitted outgoing modes. Channel 1 leaves the target unchanged; channel 2 excites it. Let and . For the same projectile mass, conservation of energy then requires a target excitation . These are hypothetical parameters, not a fit to a particular atom. Each incoming or outgoing mode is normalized to unit probability current; the outgoing modes will have the angular profiles specified below.
Choose the S-matrix
This phenomenological unitary matrix specifies the channel conversion; it is not derived here from an interaction potential energy. For incidence in channel 1, its first column gives the outgoing flux amplitudes
To recover ordinary wave coefficients relative to an incident coefficient one, undo the velocity normalization:
Thus is allowed: channel 2 travels more slowly. Its outgoing probability current fraction is , not . Renormalizing the unweighted squares would give and , the wrong probabilities. More generally the ordinary coefficient map is , with ; it obeys . The conserved quantity is the squared coefficient norm weighted by velocity, which represents probability current.
Now specify an original angular realization of these two outgoing modes, with scattering angle and no azimuth dependence:
The first normalization follows from , where . These angular profiles define the modes of this model; they are not asserted to follow from Born's approximation or an atomic potential energy. The joint directional probability densities are
Figure 1.
Complex angular amplitudes become channel-resolved scattering intensity
. Original two-channel model. Left: arrows in the complex plane show , the angular factor of . Right: the joint directional probability densities , per solid angle, retain the channel flux weights and . The shaded interval is the cone ; its probability requires the measure , not the area under a curve against degrees. These are chosen amplitudes, not historical measurements. Source. Details.
Suppose a detector accepts the cone , with full azimuth and unit efficiency. Its angular fractions within the channels are
Multiplying the angular fractions by the outgoing flux weights gives
For independent incident events, the expected values are therefore and accepted detections, compared with and over all directions. The remaining expected events fall outside the cone. The chosen incident number makes the expectations integers; actual counts fluctuate. Conditional on a detection inside the cone, the channel probabilities are and , rather than and . Angular acceptance changes the observed channel proportions even with equal detector efficiencies.
Figure 2.
Velocity weighting and angular acceptance determine expected channel counts
. Original model with incident events. Left: raw squared wave coefficients and become conserved outgoing flux fractions and after multiplying by . Right: all-direction expectations and become and inside the cone. Bars show calculated expected values, not a generated sample or experimental counts. Source. Details.
The final atomic states are orthogonal, so ignoring the atomic readout adds the two channel intensities; it does not produce . Coherent paths that reach the same final state can instead interfere, as the existing probability amplitude and quantum superposition examples demonstrate. The distinction is physical distinguishability of final states, not whether the detector operator chooses to display their labels. Flux normalization and angular integration make the scattering interpretation useful: they turn complex outgoing waves into conserved event weights and then into the statistics of a specified detector.