Skip to main content
Sandbox Physics

FR07-B · Different wells behind the same pattern

Nuclear Optical Potential

Can one scattering pattern determine a nuclear potential? Several nuclear potentials can fit the same finite set of scattering counts. Their predictions often separate at angles that were left out of the fit, so a wider angular range can distinguish the competing explanations.

Send a neutron wave toward a model calcium nucleus. Compare potentials that fit the same forward-angle counts, then use backward angles to separate them.

Interactive modelNuclear Optical Potential
Elastic cross section2.32 b2.32\,\mathrm b
Loss from elastic channel1.10 b1.10\,\mathrm b
Poisson deviance4.954.95
At the detector39.4 fm2 sr−139.4\,\mathrm{fm}^{2}\,\mathrm{sr}^{-1}

Saved synthetic Poisson record · seed 3707

Similar forward, different backward

Five forward angles enter the fit.

Gold shows your potential; violet, blue, coral and green show candidates. Pale points are fitted, gray points held out; bars approximate counting errors.

dσdΩ [fm2 sr−1],1 b=100 fm2\frac{d\sigma}{d\Omega}\ [\mathrm{fm}^{2}\,\mathrm{sr}^{-1}],\quad 1\,\mathrm b=100\,\mathrm{fm}^{2}

Upper bars use ten angles; lower bars use five. Shorter bars mean smaller deviance.

Inspect the depth profile

At each depth, radius and absorption are searched on the same finite grid. Each curve subtracts its own lowest grid value. Gold uses five angles; blue uses all ten. Multiple valleys show the remaining ambiguity. The grid is coarse and these curves do not define calibrated confidence intervals.

Scattering normalization and sources

Vertical counting-error bars use the square root of count plus one divided by exposure. They are an approximate visual scale; the fit uses the Poisson count likelihood directly.

Δ ⁣(dσdΩ)=n+1E\Delta\!\left(\frac{d\sigma}{d\Omega}\right)=\frac{\sqrt{n+1}}{\mathcal E}

The residual panel automatically rescales symmetrically to include every unmodified value, and its ticks show the current range. Each count residual is divided by the square root of count plus one.

Candidate numbers follow the exported array. These local fits illustrate ambiguity; their spread does not constitute a confidence interval.

U(r)={−V−iWr<R0r≥R,f(θ)=12ik∑ℓ=012(2ℓ+1)(Sℓ−1)Pℓ(cos⁡θ)U(r)=\begin{cases}-V-iW&r<R\\0&r\geq R\end{cases},\quad f(\theta)=\frac1{2ik}\sum_{\ell=0}^{12}(2\ell+1)(S_\ell-1)P_\ell(\cos\theta)
σel=πk2∑ℓ(2ℓ+1)∣Sℓ−1∣2,σreac=πk2∑ℓ(2ℓ+1)(1−∣Sℓ∣2)\sigma_{\rm el}=\frac\pi{k^2}\sum_\ell(2\ell+1)|S_\ell-1|^2,\quad\sigma_{\rm reac}=\frac\pi{k^2}\sum_\ell(2\ell+1)(1-|S_\ell|^2)

The entrance wave travels to the right; detector angles are measured from that forward direction in the centre-of-mass frame. The reduced mass uses the neutron mass and an explicitly approximate calcium mass of forty atomic mass units. Twelve orbital partial waves are retained beyond the s wave. The record stays at one energy when a potential changes.

The grid search and local coordinate refinement use the Poisson count likelihood. A second search fits all ten angles. The omitted channels are represented only by lost elastic flux; absorption does not specify the probabilities of particular reactions.

Study nuclear structure, scattering and stellar rates

Physics tutorial

Nuclear Optical Potential

BackgroundA neutron wave changes phase as it crosses an attractive nuclear potential. A complex part of that potential also removes flux from the elastic channel.

Why it mattersWe measure scattering counts, then work backward to a potential. Finite angular coverage can leave very different interiors consistent with the same observations.

Start with the essentials

Focus question
Can one scattering pattern determine a nuclear potential?
One-sentence intuition
Good agreement at the measured angles can leave competing potentials that disagree sharply at unmeasured angles.

Core mathematical model

A potential with attraction and absorption

U(r)={−V−iWr<R0r≥R,V,W≥0U(r)=\begin{cases}-V-iW&r<R\\0&r\geq R\end{cases},\quad V,W\geq0

The real depth changes the elastic wave phase. The negative imaginary part reduces the outgoing elastic flux. Radius is measured in femtometres; the energy and both depths are in MeV.

From partial waves to an angular distribution

f(θ)=12ik∑ℓ=012(2ℓ+1)(Sℓ−1)Pℓ(cos⁡θ)dσdΩ=∣f(θ)∣2\begin{aligned}f(\theta)&=\frac1{2ik}\sum_{\ell=0}^{12}(2\ell+1)(S_\ell-1)P_\ell(\cos\theta)\\\frac{d\sigma}{d\Omega}&=|f(\theta)|^2\end{aligned}

Each orbital angular momentum contributes a complex scattering amplitude. Their interference makes the angle dependence. The same matrices determine the exterior wave shown on the stage.

Compare predictions with counts

D=2∑i[λi−ni+niln⁡niλi],λi=E ∣f(θi)∣2D=2\sum_i\left[\lambda_i-n_i+n_i\ln\frac{n_i}{\lambda_i}\right],\quad\lambda_i=\mathcal E\,|f(\theta_i)|^2

The saved count is n; its predicted Poisson mean is lambda. Exposure combines beam, target, accepted solid angle and efficiency. A zero observed count contributes twice the predicted mean. A smaller deviance indicates a better fit to the selected angles.

Common difficulties

Absorption refers to one channel

Typical misconceptionAn absorbed neutron has vanished from the full physical system.

Better mental modelThe imaginary potential represents probability flowing into omitted reaction and excitation channels. This model reports their combined loss from the elastic channel.

A fitted curve can have several explanations

Typical misconceptionA small fit statistic proves the potential parameters are uniquely known.

Better mental modelThe forward counts constrain combinations of phases. Several potentials can reproduce those combinations. Additional angles or energies test the alternatives.

Run the experiment

  1. 01

    Compare the forward-angle twins

    Predict whether a much shallower well must spoil the fit. Compare Forward-angle twins with A shallower look-alike.

    What to observe: The potential profiles change substantially. The forward points and their fitted curves stay close, and the forward deviance hardly changes.
  2. 02

    Measure where predictions separate

    Drag the detector toward a backward angle. Then turn on Use all ten angles in the comparison.

    What to observe: The held-out counts penalize the shallow candidate much more strongly. These new angles carry information missing from the forward pattern.
  3. 03

    Fit the wider angular range

    Choose Backward angles decide, then export the record and fits. Compare the included angles and the unchanged saved counts.

    What to observe: A separate search fits all ten saved angles. Candidate spread describes alternatives found by this finite search; it does not supply a statistical confidence interval.

Related experiments and articles

Nuclear physics · Reading path

From nuclear binding to stellar reaction rates

Articles on this path that use this experiment, followed by the other experiments they suggest.

  1. ArticleWhere the wave goes after a neutron meets a nucleus
  2. ArticleWhy different nuclear potentials can fit the same scattering
FeedbackSomething broken, confusing, or missing on this page?Tell us
What kind of feedback?

Sent with this page’s address, its Lab ID, the site language and your device type. Please leave out personal details. Feedback is kept for up to 12 months.Privacy policy