research32 min studyLesson 986 in Physics

Inelastic neutron scattering and dynamical structure factors

Define energy transfer and correlation normalization before applying detailed balance or fluctuation–dissipation.

On this page

What you’ll understand

  • Define a structure factor per unit energy with positive sample energy gain.
  • Apply detailed balance with momentum reversal and a declared retarded-response sign.
  • Distinguish nuclear-density, magnetic-spin, and X-ray charge response from measured detector counts.

An inelastic scattering experiment transfers both momentum and energy to a material. Its intensity samples correlations of the operator to which the probe couples, filtered by kinematics, form factors, polarization, and instrument response. A bright feature is therefore not an unqualified density of excitations. Before comparing spectra or invoking a fluctuation–dissipation relation, one must declare which operator fluctuates, whether frequency or energy is the horizontal variable, and which sign means that the sample gains energy.

Specify the correlation and its Fourier convention

Use momentum transfer q=ki−kf\mathbf q=\mathbf k_i-\mathbf k_f and sample energy gain E=Ei−EfE=E_i-E_f. Positive E means neutron energy loss. Let Aq†=A−qA_{\mathbf q}^\dagger=A_{-\mathbf q} be an extensive Fourier component of the chosen sample operator, with N equivalent sites. Define the unsymmetrized structure factor per site and per energy by the correlation below. Its spectral representation has positive-energy transitions from initially occupied lower-energy states to higher-energy states, weighted by their thermal populations.

Integrating over all E gives the equal-time correlation divided by N. A spectrum defined per angular frequency instead differs by a factor of hbar, so normalization cannot be inferred from the symbol S alone. For a scalar Hermitian density the spectral weights are nonnegative; magnetic scattering involves a tensor of spin correlations and a transverse polarization projection. Elastic or disconnected contributions also need to be retained or explicitly subtracted when comparing integrated sum rules.

SA(q,E)=12πℏN∫−∞∞dt eiEt/ℏ⟨Aq(t)A−q(0)⟩,∫dE SA=⟨AqA−q⟩NS_A(\mathbf q,E)=\frac{1}{2\pi\hbar N}\int_{-\infty}^{\infty}dt\,e^{iEt/\hbar}\langle A_{\mathbf q}(t)A_{-\mathbf q}(0)\rangle,\qquad \int dE\,S_A=\frac{\langle A_{\mathbf q}A_{-\mathbf q}\rangle}{N}
S_A has units of A squared per energy per site, using the stated extensive-operator convention. This is an unsymmetrized equilibrium correlation, not a symmetrized noise spectrum.↔ Scroll wide equations horizontally.

Detailed balance and response need matching signs

Interchanging the initial and final states in a thermal spectral sum gives detailed balance with both momentum and energy reversed. At positive sample energy gain, the reverse process is suppressed because it requires a thermally populated excited state. Replacing minus q by plus q requires an additional symmetry such as inversion for the chosen channel; it must not be done silently in a nonreciprocal magnetic material. Dividing out a Bose factor does not manufacture information outside the measured kinematic range.

Define the retarded response to a perturbation minus field times the conjugate operator with the positive-i commutator convention below, and Fourier transform with the same positive exponential. Its dissipative imaginary part is positive at positive energy in a passive scalar channel. The spectral sum then produces the displayed fluctuation–dissipation factor. Other authors choose a negative-i Green function or integrate frequency rather than energy; their signs and hbar factors must be translated before their formulas are combined with these definitions.

χAR(q,t)=iθ(t)ℏN⟨[Aq(t),A−q(0)]⟩,SA(−q,−E)=e−E/(kBT)SA(q,E),χA′′(q,E)=π[1−e−E/(kBT)]SA(q,E)\chi_A^R(\mathbf q,t)=\frac{i\theta(t)}{\hbar N}\langle[A_{\mathbf q}(t),A_{-\mathbf q}(0)]\rangle,\quad S_A(-\mathbf q,-E)=e^{-E/(k_BT)}S_A(\mathbf q,E),\quad \chi_A''(\mathbf q,E)=\pi[1-e^{-E/(k_BT)}]S_A(\mathbf q,E)
chi is the Fourier transform of this retarded kernel with exp(iEt/hbar). The equations assume thermal equilibrium and consistent operator normalization. A different response-sign convention changes the stated imaginary-part sign.↔ Scroll wide equations horizontally.

The probe determines the operator and prefactors

For identical nuclei with real coherent scattering length b, neglecting isotope and spin incoherence, a Born single-scattering nuclear cross section is proportional to the density structure factor, with an outgoing-to-incoming wave-vector ratio. Magnetic neutrons instead probe spin components perpendicular to momentum transfer, multiplied by magnetic form factors and the appropriate magnetic coupling. Ordinary nonresonant X rays predominantly probe charge density with electronic form factors and polarization factors. Resonant X-ray channels can carry other effective operators, so none of these signals is universally interchangeable.

Detector efficiency, monitor normalization, sample mass, absorption, multiple scattering, backgrounds, and resolution must be handled before absolute spectral weights are compared. A finite-resolution elastic line can leak into low-energy inelastic intensity. Momentum-integrated sum rules require sufficient reciprocal-space and energy coverage, not merely one strong scan. Compare measured counts to a resolution-convolved cross section; a fitted dispersion alone does not establish the nature, uniqueness, or fractionalization of the excitation carrying that spectral weight.

d2σnucdΩ dEf=Nb2kfkiSn(q,Ei−Ef)\frac{d^2\sigma_{\rm nuc}}{d\Omega\,dE_f}=N b^2\frac{k_f}{k_i}S_n(\mathbf q,E_i-E_f)
Here A_q is the dimensionless nuclear number-density Fourier component and S_n uses the first equation. The differential cross section has units of area per steradian per energy. This simplified expression excludes magnetic, incoherent, absorption, and multiple-scattering contributions.↔ Scroll wide equations horizontally.

A positive-energy peak and its reverse process

In a hypothetical equilibrium scalar channel, E/(kBT)=ln⁡3E/(k_BT)=\ln3 and SA(q,E)=6 meV−1S_A(\mathbf q,E)=6\,\mathrm{meV}^{-1} for a dimensionless operator per site. Find the reversed structure factor and the dissipative response in the declared convention.

  1. The thermal factor is e−E/(kBT)=1/3e^{-E/(k_BT)}=1/3. Detailed balance compares the measured point with (−q,−E)(-\mathbf q,-E), not automatically with the same momentum.

  2. The reverse spectral weight is (1/3)6=2 meV−1(1/3)6=2\,\mathrm{meV}^{-1}. It is smaller because the reverse transition starts in a less populated excited state.

  3. The fluctuation–dissipation factor gives χA′′=π(1−1/3)6=4π meV−1\chi_A''=\pi(1-1/3)6=4\pi\,\mathrm{meV}^{-1}. The positive sign follows from the specified response convention and positive sample energy gain.

The result

The reversed weight is two inverse millielectronvolts, and the dissipative response is four pi inverse millielectronvolts in this normalization.

Explore the model

Interactive lab

Interactive lab · scattering

A prefactor is distinct from a material response

For each fresh paused run Reset and set coupling g and angle to (1,90)(1,90) and (1,0)(1,0). Predict schematic weights 1 and 2. This instrument uses only g to the fourth times one plus cosine squared; it is not a neutron or X-ray cross section, magnetic polarization projector, structure-factor calculation, or detector simulation.

A prefactor is distinct from a material responseA schematic term diagram with dimensionless weight g⁴(1 + cos²θ). Lines denote terms, not particle paths; this is not a process cross-section or loop calculation. Current readings: Coupling factor g⁴: 0.24; Angular factor 1 + cos²θ: 1.329; Illustrative weight: 0.319. All parameters are adjustable with the labeled controls below.Incoming termsOutgoing termsθSchematic interaction term · no particle trajectories
Incoming termsOutgoing termsAngle θ
Coupling factor g⁴
0.24
Angular factor 1 + cos²θ
1.329
Illustrative weight
0.319
t = 0 s

A schematic term diagram with dimensionless weight g⁴(1 + cos²θ). Lines denote terms, not particle paths; this is not a process cross-section or loop calculation. 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 measured intensity combines a response with probe-dependent prefactors.

Check your understanding

With positive E defined as sample energy gain, what is the general equilibrium detailed-balance relation?

Try it yourself

For a hypothetical purely coherent nuclear channel let Nb2=4 barnNb^2=4\,\mathrm{barn}, kf/ki=1/2k_f/k_i=1/2, and Sn=3 meV−1S_n=3\,\mathrm{meV}^{-1}. Compute the differential cross section. Would doubling background counts double this intrinsic S?

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

Give me a hint

Multiply the three physical factors before considering an additive instrumental background.

Compare with the solution

The cross section is 4(1/2)3=6 barn sr−1meV−14(1/2)3=6\,\mathrm{barn}\,\mathrm{sr}^{-1}\mathrm{meV}^{-1} under the stated normalization. Additive background changes measured counts and must be modeled or subtracted; it does not change the intrinsic structure factor of the stipulated sample channel.

Take it with you

  • Energy versus frequency normalization and the transfer sign must be stated explicitly.
  • Detailed balance and fluctuation–dissipation connect compatible equilibrium quantities.
  • Nuclear, magnetic, and charge probes couple to different operators and prefactors.

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.