Physics World Model — Modality Catalog
3 imaging modalities with descriptions, experimental setups, and reconstruction guidance.
Muon Tomography
Muon tomography uses naturally occurring cosmic-ray muons (mean energy ~4 GeV, flux ~1/cm2/min at sea level) to image the interior of large, dense objects by measuring the scattering angle of each muon as it traverses the object. High-Z materials (uranium, plutonium, lead) cause large-angle scattering that is readily distinguished from low-Z materials. Position-sensitive detectors (drift tubes, RPCs) above and below the object track each muon's trajectory. The scattering density is proportional to Z^2/A. Reconstruction uses the point-of-closest-approach (POCA) algorithm or maximum-likelihood/expectation-maximization (ML-EM). Long exposure times (minutes to hours) are needed due to the low natural muon flux. Applications include nuclear material detection and volcano interior imaging (muography).
Muon Tomography
Description
Muon tomography uses naturally occurring cosmic-ray muons (mean energy ~4 GeV, flux ~1/cm2/min at sea level) to image the interior of large, dense objects by measuring the scattering angle of each muon as it traverses the object. High-Z materials (uranium, plutonium, lead) cause large-angle scattering that is readily distinguished from low-Z materials. Position-sensitive detectors (drift tubes, RPCs) above and below the object track each muon's trajectory. The scattering density is proportional to Z^2/A. Reconstruction uses the point-of-closest-approach (POCA) algorithm or maximum-likelihood/expectation-maximization (ML-EM). Long exposure times (minutes to hours) are needed due to the low natural muon flux. Applications include nuclear material detection and volcano interior imaging (muography).
Principle
Muon tomography uses naturally occurring cosmic-ray muons to image the internal density structure of large objects (buildings, volcanoes, cargo containers). Muons undergo multiple Coulomb scattering, with the scattering angle proportional to the areal density and atomic number of the traversed material. By measuring the incoming and outgoing muon trajectories, the density distribution inside the object can be tomographically reconstructed.
How to Build the System
Place tracking detectors (drift tubes, scintillator strips, resistive plate chambers, or GEM detectors) above and below (or around) the object to be imaged. Each detector station measures the position and angle of each cosmic-ray muon before and after it traverses the object. Typical cosmic-ray muon flux is ~10,000 muons/m²/min at sea level. Exposure times range from minutes (for dense nuclear materials) to months (for geological structures like volcanoes).
Common Reconstruction Algorithms
- Point of Closest Approach (POCA) voxel reconstruction
- Maximum Likelihood / Expectation Maximization (ML/EM) scattering tomography
- Angle Statistics Reconstruction (ASR) for material discrimination
- Binned scattering density reconstruction
- Deep-learning muon tomography for faster convergence with fewer muons
Common Mistakes
- Insufficient muon statistics for the desired spatial resolution (need long exposure)
- Detector alignment errors causing incorrect scattering angle measurements
- Not accounting for muon momentum spectrum (affects scattering angle distribution)
- Background tracks (electrons, low-momentum muons) contaminating the data
- POCA algorithm limitations in complex, non-point-like geometries
How to Avoid Mistakes
- Calculate required exposure time based on object size, density, and desired resolution
- Align detectors carefully using straight-through cosmic ray tracks as calibration
- Use momentum measurement (from curvature in a magnetic field) or momentum-dependent MCS model
- Apply track quality cuts (chi-squared, minimum number of detector hits) to reject background
- Use iterative reconstruction (ML/EM) rather than POCA for quantitative density imaging
Forward-Model Mismatch Cases
- The widefield fallback applies Gaussian blur, but muon tomography measures the scattering angle of cosmic-ray muons passing through the object — the scattering angle (Highland formula) encodes radiation length and density, not image blur
- Muon tomography uses natural cosmic-ray flux (~10,000 muons/m^2/min) with tracking detectors above and below the object — the widefield optical model has no connection to high-energy particle tracking or multiple Coulomb scattering physics
How to Correct the Mismatch
- Use the muon tomography operator that models multiple Coulomb scattering: incoming and outgoing muon tracks are measured, and the scattering angle distribution at each voxel encodes the local radiation length (related to material Z and density)
- Reconstruct using POCA (Point of Closest Approach) for quick imaging, or ML/EM iterative methods for quantitative density/Z mapping, using the correct scattering probability forward model
Experimental Setup — Signal Chain
Experimental Setup — Details
Key References
- Borozdin et al., 'Radiographic imaging with cosmic-ray muons', Nature 422, 277 (2003)
- Tanaka et al., 'Imaging the conduit size of the dome with cosmic-ray muons: The structure beneath Showa-Shinzan Lava Dome', Geophysical Research Letters 34, L22311 (2007)
Canonical Datasets
- Los Alamos muon tomography simulation benchmarks
- IAEA muon imaging reference data
Neutron Radiography / Tomography
Neutron imaging exploits the unique interaction of thermal neutrons with matter — neutrons are attenuated strongly by light elements (hydrogen, lithium, boron) while penetrating heavy elements (lead, iron) that are opaque to X-rays. The forward model follows Beer-Lambert: I = I_0 * exp(-integral(Sigma(s) ds)) where Sigma is the macroscopic cross-section. Tomographic reconstruction from multiple projection angles uses FBP or iterative methods. Neutron sources include research reactors and spallation sources. The lower flux compared to X-rays requires longer exposures (seconds) and results in lower spatial resolution (50-100 um).
Neutron Radiography / Tomography
Description
Neutron imaging exploits the unique interaction of thermal neutrons with matter — neutrons are attenuated strongly by light elements (hydrogen, lithium, boron) while penetrating heavy elements (lead, iron) that are opaque to X-rays. The forward model follows Beer-Lambert: I = I_0 * exp(-integral(Sigma(s) ds)) where Sigma is the macroscopic cross-section. Tomographic reconstruction from multiple projection angles uses FBP or iterative methods. Neutron sources include research reactors and spallation sources. The lower flux compared to X-rays requires longer exposures (seconds) and results in lower spatial resolution (50-100 um).
Principle
Neutron radiography and tomography image the transmission of a thermal or cold neutron beam through a sample. Neutrons interact with nuclei (not electrons), providing complementary contrast to X-rays: hydrogen-rich materials (water, polymers, organics) attenuate neutrons strongly, while metals like aluminum and lead are relatively transparent. Tomographic reconstruction from multiple projection angles yields 3-D maps of neutron attenuation.
How to Build the System
Access a research reactor or spallation neutron source with an imaging beamline (e.g., ICON at PSI, IMAT at ISIS, NIST BT-2). A collimated neutron beam (thermal or cold, 1-10 Å) passes through the sample, and a scintillator-camera system (⁶LiF/ZnS screen + sCMOS camera) records the transmitted intensity. Rotate the sample through 180° or 360° for tomography. Spatial resolution is typically 20-100 μm, limited by beam divergence and scintillator thickness.
Common Reconstruction Algorithms
- Filtered back-projection (FBP) adapted for neutron tomography
- Iterative reconstruction (SIRT, CGLS) for limited-angle or noisy data
- Beam hardening correction for polychromatic neutron spectra
- Scattering correction (point-scattered function approach)
- Neutron phase-contrast tomography (grating interferometry)
Common Mistakes
- Scattering from hydrogen-rich samples producing artifacts (halo around sample)
- Beam hardening (spectral hardening) not corrected for polychromatic beams
- Activation of sample materials, creating radiation safety issues post-experiment
- Gamma contamination in the beam degrading image quality
- Insufficient exposure time per projection, yielding noisy tomograms
How to Avoid Mistakes
- Apply scattering correction algorithms; use thin or diluted hydrogen-rich samples
- Correct beam hardening with polynomial methods or by using a velocity selector (monochromatic)
- Check sample activation potential before irradiation; use short-lived isotope-free materials
- Use gamma-blind detectors (⁶Li glass) or filters to reject gamma contamination
- Optimize exposure per projection for adequate SNR; total scan time often 2-8 hours
Forward-Model Mismatch Cases
- The widefield fallback applies optical Gaussian blur, but neutron tomography measures neutron transmission (I = I_0 * exp(-sigma_t * n * t)) — neutrons interact with nuclei, not electron clouds, giving completely different contrast (hydrogen-rich materials are opaque to neutrons but transparent to X-rays)
- Neutron attenuation depends on nuclear cross-sections that vary dramatically between isotopes (H, Li, B are strong absorbers) — the widefield model has no nuclear physics and cannot distinguish materials by their neutron interaction properties
How to Correct the Mismatch
- Use the neutron tomography operator implementing Beer-Lambert neutron transmission: y(theta,s) = I_0 * exp(-integral(Sigma_t(x,y) dl)) where Sigma_t is the macroscopic total cross-section
- Reconstruct using FBP or iterative methods (same algorithms as X-ray CT) but with neutron-specific attenuation coefficients — neutron imaging reveals hydrogen/water content, lithium batteries, and metallurgical features invisible to X-rays
Experimental Setup — Signal Chain
Experimental Setup — Details
Key References
- Kardjilov et al., 'Advances in neutron imaging', Materials Today 21, 652-672 (2018)
- IAEA, 'Neutron Imaging: A Non-Destructive Tool for Materials Testing', IAEA-TECDOC-1604 (2008)
Canonical Datasets
- PSI ICON neutron imaging benchmark data
- NIST neutron radiography reference images
Proton Radiography
Proton radiography/CT uses high-energy proton beams (100-250 MeV) to image the relative stopping power (RSP) of tissue, which is the quantity directly needed for proton therapy treatment planning. Unlike X-rays which measure attenuation, proton imaging measures the energy loss and scattering of individual protons as they traverse the object. Each proton's entry/exit position and angle are tracked, and the residual energy is measured. The RSP is reconstructed from many proton histories using iterative algorithms. Challenges include multiple Coulomb scattering (which blurs the spatial resolution to ~1 mm) and the need for single-proton tracking at high rates.
Proton Radiography
Description
Proton radiography/CT uses high-energy proton beams (100-250 MeV) to image the relative stopping power (RSP) of tissue, which is the quantity directly needed for proton therapy treatment planning. Unlike X-rays which measure attenuation, proton imaging measures the energy loss and scattering of individual protons as they traverse the object. Each proton's entry/exit position and angle are tracked, and the residual energy is measured. The RSP is reconstructed from many proton histories using iterative algorithms. Challenges include multiple Coulomb scattering (which blurs the spatial resolution to ~1 mm) and the need for single-proton tracking at high rates.
Principle
Proton radiography images the transmission and scattering of high-energy protons (50-800 MeV) through dense objects. Unlike X-rays, protons undergo significant multiple Coulomb scattering (MCS) in matter, which provides density and compositional contrast. Both transmission (energy loss) and scattering angle measurements contribute to image formation. Proton radiography can penetrate very dense materials (steel, depleted uranium) that are opaque to X-rays.
How to Build the System
Requires a high-energy proton accelerator facility (synchrotron or cyclotron delivering 200-800 MeV protons). The object is placed in the beam path between tracking detectors (silicon strip or GEM detectors) that measure each proton's position and angle before and after the object. A magnetic spectrometer (quadrupole lens system, e.g., at LANL pRad facility) focuses transmitted protons onto a scintillator + camera detector.
Common Reconstruction Algorithms
- Most Likely Path (MLP) estimation for proton CT reconstruction
- Filtered back-projection with scattering-angle weighting
- Algebraic reconstruction (ART) with MCS forward model
- Material discrimination from dual-parameter (transmission + scattering) analysis
- Deep-learning proton CT reconstruction for reduced view angles
Common Mistakes
- Ignoring multiple Coulomb scattering in the reconstruction model, causing blur
- Nuclear interaction losses (protons stopped or scattered out of detector acceptance)
- Insufficient proton statistics leading to noisy images
- Energy straggling not modeled, causing depth-of-field blur in radiography
- Detector alignment errors between upstream and downstream tracking systems
How to Avoid Mistakes
- Use MLP or cubic spline path estimation in iterative reconstruction algorithms
- Account for nuclear interaction losses in the forward model; filter outlier tracks
- Accumulate sufficient proton histories (>10⁶ for radiography, >10⁸ for proton CT)
- Include energy straggling in the forward model or use higher energy protons to reduce it
- Carefully align tracking detectors with survey or use track-based alignment algorithms
Forward-Model Mismatch Cases
- The widefield fallback applies Gaussian blur, but proton radiography measures energy loss and multiple Coulomb scattering (MCS) of high-energy protons traversing the object — the scattering angle distribution encodes areal density, not spatial blur
- Protons lose energy continuously (Bethe-Bloch formula: -dE/dx ~ Z/A * z^2/beta^2) and scatter via Coulomb interaction — the measurement combines transmission intensity, residual energy, and scattering angle, none of which are modeled by optical blur
How to Correct the Mismatch
- Use the proton radiography operator that models energy-dependent proton transport: energy loss via Bethe-Bloch stopping power and angular broadening via Highland MCS formula (theta_rms ~ 13.6 MeV/(p*v) * sqrt(t/X_0))
- Reconstruct water-equivalent path length (WEPL) maps from residual energy measurements, or use scattering radiography for material discrimination — essential for proton therapy treatment planning
Experimental Setup — Signal Chain
Experimental Setup — Details
Key References
- Schulte et al., 'Conceptual design of a proton computed tomography system for applications in proton radiation therapy', IEEE Trans. Nucl. Sci. 51, 866-872 (2004)
Canonical Datasets
- Simulated proton CT phantoms (Penfold et al.)