Inelastic Neutron Scattering (INS)

• INS is a spectroscopy technique that measures momentum and energy transfers to directly probe the dynamical structure factor in materials.
• It is used to study magnetic excitations, phonon dynamics, and crystal-field transitions in systems like iridates, rare-earth compounds, and nanoclusters.
• Advanced implementations such as time-of-flight and polarized INS, along with specialized environments, enable precise extraction of microscopic parameters.

Inelastic neutron scattering (INS) is a momentum- and energy-resolved spectroscopy in which a neutron exchanges momentum transfer ${\bf Q}$ and energy transfer $\hbar\omega$ with a target. In the studies considered here, INS is used as a direct probe of the dynamical structure factor $S({\bf Q},\hbar\omega)$, of the scattering law $S(\alpha,\beta)$ for isotropic materials, and of related magnetic and vibrational response functions. Across contemporary applications, it resolves dispersive spin waves, nearly dispersionless crystal electric field (CEF) transitions, localized molecular spin levels, phonon polarization states, hydrogen vibrational overtones, and laser-driven nonequilibrium structure factors (Fujita et al., 2020, Cai et al., 2018, Anand et al., 2017, Wang et al., 27 Jul 2025, Shiga et al., 2024, Hua et al., 2023).

1. Fundamental observables and kinematics

In one common formulation, INS directly measures the dynamical structure factor

$S({\bf Q},\hbar\omega),$

as emphasized in the study of high-energy magnetic excitations in Sr$_3$Ir$_2$O$_7$ (Fujita et al., 2020). For isotropic materials, the inelastic response can also be written through the scattering kernel $S(\alpha,\beta)$, which governs the double differential cross section per atom,

$$\frac{\mathrm d^2 \sigma}{\mathrm dE'\,\mathrm d\Omega} = \frac{\sigma_b}{4\pi}\frac{k'}{k} S(Q,\omega),$$

with reduced variables

$\hbar\omega$0

The kinematically allowed region is

$\hbar\omega$1

where

$\hbar\omega$2

This representation is central in Monte Carlo transport treatments of INS (Cai et al., 2018).

For magnetic spectroscopy in the dipole approximation, one expression used for the INS cross section is

$\hbar\omega$3

with

$\hbar\omega$4

This form underlies the identification of magnetic CEF scattering in Nd$\hbar\omega$5Zn$\hbar\omega$6 (Yamamoto et al., 2023). For hydrogen in Pd, by contrast, the dominant quantity is the incoherent dynamic structure factor,

$\hbar\omega$7

which probes local vibrational motion and its quantum dynamics (Shiga et al., 2024).

Standard INS is treated in the coincidence-INS proposal as a first-order probe of single-spin dynamics, with the measured response proportional to a single-spin dynamical structure factor or susceptibility $\hbar\omega$8 (Su et al., 2020). This makes explicit a recurring feature of INS: the technique is defined not by a single class of excitations, but by a common kinematic framework that can access magnetic, orbital, and vibrational dynamics through different response functions.

2. Experimental realizations and instrumental modalities

The experimental implementations surveyed here span time-of-flight, triple-axis, polarized, high-pressure, and pump-probe configurations, using both powders and single crystals.

Modality	Example implementation	Primary use
Time-of-flight single-crystal INS	4SEASONS on Sr$\hbar\omega$9Ir$S({\bf Q},\hbar\omega)$0O$S({\bf Q},\hbar\omega)$1; Pelican on atacamite	Dispersive magnetic modes
Time-of-flight powder INS	MARI on Nd$S({\bf Q},\hbar\omega)$2Zn$S({\bf Q},\hbar\omega)$3 and CeNiC$S({\bf Q},\hbar\omega)$4; OSIRIS on $S({\bf Q},\hbar\omega)$5	CEF and localized levels
Polarized INS	Ba$S({\bf Q},\hbar\omega)$6K$S({\bf Q},\hbar\omega)$7Fe$S({\bf Q},\hbar\omega)$8As$S({\bf Q},\hbar\omega)$9	Spin-space anisotropy
Specialized environments	Zr-BMG high-pressure cell; laser-neutron HYSPEC	Pressure and nonequilibrium INS

Representative time-of-flight measurements include Sr$S(\alpha,\beta)$0Ir$S(\alpha,\beta)$1O$S(\alpha,\beta)$2 on 4SEASONS at J-PARC/MLF with $S(\alpha,\beta)$3 meV and $S(\alpha,\beta)$4 meV at $S(\alpha,\beta)$5 meV, using about 500 flux-grown single crystals with total mass $S(\alpha,\beta)$6 g mounted with the $S(\alpha,\beta)$7-axis parallel to the incident neutron beam to reduce absorption (Fujita et al., 2020). Powder INS on MARI at ISIS was used for Nd$S(\alpha,\beta)$8Zn$S(\alpha,\beta)$9 with incident energies $S({\bf Q},\hbar\omega),$0 and $S({\bf Q},\hbar\omega),$1 meV and for CeNiC$S({\bf Q},\hbar\omega),$2 with $S({\bf Q},\hbar\omega),$3, $S({\bf Q},\hbar\omega),$4, and $S({\bf Q},\hbar\omega),$5 meV (Yamamoto et al., 2023, Bhattacharyya et al., 2014). Very-high-resolution low-energy spectroscopy of the partially deuterated $S({\bf Q},\hbar\omega),$6 nanocluster employed OSIRIS with fixed final neutron energy $S({\bf Q},\hbar\omega),$7 meV and instrumental resolution about $S({\bf Q},\hbar\omega),$8–$S({\bf Q},\hbar\omega),$9 (Vaknin et al., 2014).

The measurements rely heavily on geometry control and background management. Low-$_3$0 integration is used to emphasize magnetic CEF scattering in Nd$_3$1Zn$_3$2, especially $_3$3 (Yamamoto et al., 2023). Nonmagnetic reference compounds—La$_3$4Hf$_3$5O$_3$6, LaNiC$_3$7, and YCo$_3$8Zn$_3$9—are measured under the same conditions to isolate phonon backgrounds or nuclear scattering (Anand et al., 2017, Bhattacharyya et al., 2014, Yamamoto et al., 2023). In the $_2$0 molecule, the Fe$_2$1 form factor changes by less than 5% over the measured $_2$2-range, enabling integration over $_2$3 to $_2$4 without losing meaningful magnetic information (Vaknin et al., 2014).

Polarization analysis extends INS from momentum-energy spectroscopy to spin-space discrimination. In Ba$_2$5K$_2$6Fe$_2$7As$_2$8, the spin-flip cross sections $_2$9, $_7$0, and $_7$1 are used to extract $_7$2 and $_7$3 (Zhang et al., 2013). Recent specialized environments further enlarge the experimental phase space: a Zr-based bulk metallic glass clamp cell provides higher neutron transparency and a broad, featureless background for high-pressure INS (Hayashida et al., 16 Feb 2026), while a synchronized laser-neutron capability at HYSPEC enables time-resolved INS of nonequilibrium magnons with a 515 nm, 5 ns pulsed laser (Hua et al., 2023).

3. Magnetic excitations: from magnons to continua

INS resolves both dispersive and localized magnetic excitations. In bilayer Sr$_7$4Ir$_7$5O$_7$6, the first INS study of high-energy magnetic excitations observed a dispersive magnon branch extending from about $_7$7 meV to $_7$8 meV. A broad peak appears near the antiferromagnetic zone center at 180 meV; on lowering energy to 130, 100, and 80 meV, the peak splits into two symmetrically displaced peaks in momentum; no visible magnetic signal below 70 meV was observed; and the gap energy at the antiferromagnetic zone boundary is $_7$9 meV at $S(\alpha,\beta)$0. The peak positions match the single-magnon dispersion previously obtained by RIXS, leading to the conclusion that INS and RIXS detect identical magnetic excitations in iridate oxides, as in La$S(\alpha,\beta)$1CuO$S(\alpha,\beta)$2 (Fujita et al., 2020).

Low-energy ordered-state spin dynamics are similarly accessible in frustrated magnets, but often require effective Hamiltonians. In atacamite, single-crystal INS at 1.5 K revealed two distinct dispersive spin-wave modes, an energy offset between them, and a spin gap of at least $S(\alpha,\beta)$3 meV; a cut at $S(\alpha,\beta)$4 r.l.u. shows peaks near 1.3 and 1.7 meV. These spectra are rationalized by a frustrated $S(\alpha,\beta)$5 zigzag-chain model for Cu(2) spins in an effective staggered mean field from Cu(1), with LSWT parameters approximately $S(\alpha,\beta)$6 meV, $S(\alpha,\beta)$7 meV, and $S(\alpha,\beta)$8 meV (Allen et al., 4 Aug 2025). In CsCrF$S(\alpha,\beta)$9, powder INS in the low-temperature phase showed magnetic excitations up to about 10 meV with nearly flat features around 0.5 meV and 1.5 meV, while the powder-averaged spectra at 0.8 K and 3.1 K were strikingly similar despite different magnetic propagation vectors. That unchanged dynamical spectrum across two ordered phases is presented as evidence for an order-by-disorder mechanism (Kikuchi et al., 2022).

Polarized INS can disentangle magnetic anisotropy that is hidden in unpolarized data. In optimally hole-doped Ba$$\frac{\mathrm d^2 \sigma}{\mathrm dE'\,\mathrm d\Omega} = \frac{\sigma_b}{4\pi}\frac{k'}{k} S(Q,\omega),$$0K$$\frac{\mathrm d^2 \sigma}{\mathrm dE'\,\mathrm d\Omega} = \frac{\sigma_b}{4\pi}\frac{k'}{k} S(Q,\omega),$$1Fe$$\frac{\mathrm d^2 \sigma}{\mathrm dE'\,\mathrm d\Omega} = \frac{\sigma_b}{4\pi}\frac{k'}{k} S(Q,\omega),$$2As$$\frac{\mathrm d^2 \sigma}{\mathrm dE'\,\mathrm d\Omega} = \frac{\sigma_b}{4\pi}\frac{k'}{k} S(Q,\omega),$$3, the normal state already shows $$\frac{\mathrm d^2 \sigma}{\mathrm dE'\,\mathrm d\Omega} = \frac{\sigma_b}{4\pi}\frac{k'}{k} S(Q,\omega),$$4 for energies below about 7 meV. In the superconducting state, the apparent difference between the previously reported spin gaps at $$\frac{\mathrm d^2 \sigma}{\mathrm dE'\,\mathrm d\Omega} = \frac{\sigma_b}{4\pi}\frac{k'}{k} S(Q,\omega),$$5 and $$\frac{\mathrm d^2 \sigma}{\mathrm dE'\,\mathrm d\Omega} = \frac{\sigma_b}{4\pi}\frac{k'}{k} S(Q,\omega),$$6 is traced to different spin gaps in the $$\frac{\mathrm d^2 \sigma}{\mathrm dE'\,\mathrm d\Omega} = \frac{\sigma_b}{4\pi}\frac{k'}{k} S(Q,\omega),$$7-axis and in-plane polarized channels, while the neutron spin resonance at $$\frac{\mathrm d^2 \sigma}{\mathrm dE'\,\mathrm d\Omega} = \frac{\sigma_b}{4\pi}\frac{k'}{k} S(Q,\omega),$$8 meV is essentially isotropic (Zhang et al., 2013).

Zero-dimensional spin systems furnish the complementary limit, in which INS resolves discrete, dispersionless transitions rather than bands. In the tridiminished-icosahedron $$\frac{\mathrm d^2 \sigma}{\mathrm dE'\,\mathrm d\Omega} = \frac{\sigma_b}{4\pi}\frac{k'}{k} S(Q,\omega),$$9 nano-cluster, a field-dependent low-energy excitation is assigned to a $\hbar\omega$00 fragment with an $\hbar\omega$01 ground state and powder-averaged $\hbar\omega$02, while two zero-field excitations near $\hbar\omega$03 meV and $\hbar\omega$04 meV, each split into three components, are assigned to $\hbar\omega$05 excitations of an $\hbar\omega$06 fragment with an $\hbar\omega$07 ground state (Vaknin et al., 2014). In the molecular dimer Ni$\hbar\omega$08, INS resolves three cold inter-multiplet transitions at $\hbar\omega$09, $\hbar\omega$10, and $\hbar\omega$11, from which the microscopic Hamiltonian parameters $\hbar\omega$12, $\hbar\omega$13, and $\hbar\omega$14 are extracted; Co$\hbar\omega$15 instead shows at least five peaks between about $\hbar\omega$16 and $\hbar\omega$17, consistent with a more complex high-spin $\hbar\omega$18 dimer (Ansbro et al., 2021).

INS also enters regimes where the elementary excitations are not conventional magnons. For the Ce pyrochlores Ce$\hbar\omega$19Sn$\hbar\omega$20O$\hbar\omega$21 and Ce$\hbar\omega$22Zr$\hbar\omega$23O$\hbar\omega$24, the proposed interpretation depends on whether the U(1) spin liquid is octupolar or dipolar. In an octupolar U(1) spin liquid, INS mainly selects the gapped spinon continuum; in a dipolar U(1) spin liquid, INS can also access the gapless gauge photon and a gapped electric-monopole continuum, and the low-energy spectrum is argued to be dominated by the electric monopole continuum because the $\hbar\omega$25 dual flux enhances the low-energy density of states whereas the photon spectral weight is linearly suppressed at low energies (Zhao et al., 2024). This indicates that INS line shapes can depend as strongly on emergent-gauge selection rules as on conventional exchange parameters.

4. Crystal-field, local vibrational, and thermodynamic spectroscopy

A major use of INS is the determination of local single-ion level schemes. In Nd$\hbar\omega$26Hf$\hbar\omega$27O$\hbar\omega$28, powder INS reveals nearly dispersionless CEF transitions at about 23.7, 34.7, and 106.5 meV, together with weaker modes at about 245.7, 265.9, and 311.9 meV. Their low-$\hbar\omega$29 intensity and suppression with increasing $\hbar\omega$30 identify them as magnetic. Fitting a Wybourne-tensor CEF Hamiltonian yields five Kramers doublets within the ground multiplet and a strongly Ising-like ground doublet with $\hbar\omega$31 and $\hbar\omega$32; after conversion to Stevens parameters, the negative dominant $\hbar\omega$33 meV is taken as the condition expected for a dipolar-octupolar doublet (Anand et al., 2017).

Closely related CEF spectroscopy appears in the caged compounds Nd$\hbar\omega$34Zn$\hbar\omega$35 $\hbar\omega$36 Co, Rh, Ir$\hbar\omega$37. At 5 K, powder INS shows two sharp, dispersionless excitations for each compound: 3.8 and 7.2 meV for Co, 3.1 and 5.8 meV for Rh, and 3.0 and 5.3 meV for Ir. Temperature-dependent spectral-weight transfer establishes a $\hbar\omega$38 ground doublet and two $\hbar\omega$39 excited quartets, specifically $\hbar\omega$40 for Co, $\hbar\omega$41 for Rh, and $\hbar\omega$42 for Ir. In NdCo$\hbar\omega$43Zn$\hbar\omega$44, a shoulder at 7.7 meV near the 7.2 meV CEF excitation is attributed to a bound state of the CEF and a low-lying optical phonon (Yamamoto et al., 2023).

In CeNiC$\hbar\omega$45, INS establishes two robust CEF excitations at about 8 and 30 meV, consistent with a localized Ce$\hbar\omega$46 $\hbar\omega$47 multiplet split into three Kramers doublets by an orthorhombic CEF. A low-energy feature near 4.5 meV appears only below $\hbar\omega$48 K and is assigned to spin waves. The same CEF analysis yields a calculated easy axis along the $\hbar\omega$49-axis and a quasielastic linewidth corresponding to an estimated Kondo temperature $\hbar\omega$50 K, so the reduced ordered moment is attributed to Kondo screening competing with antiferromagnetic order (Bhattacharyya et al., 2014).

INS is equally effective for local vibrational spectroscopy. For dilute hydrogen in fcc Pd, the target experimental spectrum at about 295 K contains a fundamental tone near 69 meV and a first overtone near 137 meV. Semiclassical Brownian chain molecular dynamics combined with a DFT-quality ANN potential reproduces the peak positions and intensities, including a shoulder near $\hbar\omega$51 meV attributed to coupling with Pd phonons. Classical MD fails to reproduce these features, whereas the semiclassical treatment shows that nuclear quantum effects act as a blue-shift for H in the octahedral site because the potential is strongly anharmonic with an even-symmetry quartic character (Shiga et al., 2024).

A different vibrational use of INS is the extraction of thermodynamic quantities from the phonon density of states. In diamond powders ranging from 5 nm detonation nanodiamond to macrodiamond, INS-derived PDOS was used to compute the constant-volume heat capacity

$\hbar\omega$52

and the isotopic beta-factor. The 5 nm and 40 nm samples show extra broad low-energy intensity below about 60 meV from adsorbed H-containing surface species, higher heat capacities than bulk diamond, and lower beta-factors; for those two nanodiamond sizes, about 50% of the heat-capacity difference from bulk is estimated to come from surface impurities rather than intrinsic phonon confinement (Shiryaev et al., 2020). This indicates that INS can connect microscopic vibrational spectra directly to equilibrium isotope thermodynamics.

5. Analysis frameworks, symmetry constraints, and inverse problems

Interpretation of INS commonly proceeds through explicit microscopic modeling. In atacamite, LSWT implemented with the Sunny package is fitted to the measured dispersions, relative intensities, spectral-weight asymmetry, and gap, whereas the full DFT-derived exchange network fails to reproduce the low-energy INS spectrum in LSWT (Allen et al., 4 Aug 2025). In CsCrF$\hbar\omega$53, powder-averaged LSWT spectra are computed with SpinW from a Hamiltonian including leg exchange $\hbar\omega$54, rung exchange $\hbar\omega$55, minute intertube coupling $\hbar\omega$56, Dzyaloshinskii-Moriya interaction $\hbar\omega$57, and single-ion anisotropy $\hbar\omega$58, yielding $\hbar\omega$59 meV, $\hbar\omega$60 meV, $\hbar\omega$61, $\hbar\omega$62, and $\hbar\omega$63 (Kikuchi et al., 2022). CEF studies similarly rely on Hamiltonian-based fitting, including Wybourne-tensor operators for Nd$\hbar\omega$64Hf$\hbar\omega$65O$\hbar\omega$66, cubic Stevens operators for Nd$\hbar\omega$67Zn$\hbar\omega$68, and orthorhombic Stevens operators for CeNiC$\hbar\omega$69 (Anand et al., 2017, Yamamoto et al., 2023, Bhattacharyya et al., 2014).

Beyond direct parameter extraction, point-group symmetry can fix the momentum-transfer dependence of magnetic INS intensities. For anisotropic spin clusters, a transition has a universal $\hbar\omega$70-dependence when the relevant operator irreducible representation appears exactly once in the decomposition of the site-spin operator space and the product $\hbar\omega$71 contains the totally symmetric irreducible representation. In such cases, the powder-averaged intensity can be written as

$\hbar\omega$72

with coefficients determined entirely by symmetry. This separation of dynamics from geometry is worked out explicitly for six- and eight-membered planar rings, the cube, and the icosahedron (Tabrizi, 2021).

INS data can also be inverted rather than merely fitted. In unconventional superconductors, the hour-glass spin-excitation spectrum is treated as an INS analogue of quasiparticle interference: the measured $\hbar\omega$73 positions of the magnetic quasiparticle-scattering peaks can be inverted to reconstruct the bulk Fermi surface and superconducting gap, using the relation

$\hbar\omega$74

for sign-changing superconducting gaps (Das et al., 2011). In neutron transport, the inverse problem takes a different form: rejection-based sampling of the kernel $\hbar\omega$75 produces continuous outgoing energy-angle distributions, agrees with analytical references, and gives acceptance above 96.9% when applied to conventional discrete double differential cross sections (Cai et al., 2018).

The same formalism extends to proposed higher-order and mode-selective variants of INS. Coincidence INS introduces two incident neutrons and two detectors in coincidence, with the coincidence probability governed by a two-spin Bethe-Salpeter wave function

$\hbar\omega$76

and is proposed as a direct probe of two-spin magnetic correlations beyond standard single-spin INS (Su et al., 2020). For chiral phonons, the calculated angle-resolved INS intensity distinguishes linear, elliptical, and circular phonons through characteristic polar patterns—“8-shaped” for linear modes, “open 8-shaped” for elliptical modes, and circular for chiral modes—while magnetic-field-induced mode splitting in CeF$\hbar\omega$77 is proposed as a direct probe of phonon magnetic moments and handedness (Wang et al., 27 Jul 2025). These developments suggest that INS analysis increasingly involves eigenvector-sensitive and correlation-sensitive observables, not only dispersion relations.

6. Constraints, specialized environments, and emerging regimes

INS remains technically constrained by neutron absorption, weak signal levels, and sample-environment backgrounds. The Sr$\hbar\omega$78Ir$\hbar\omega$79O$\hbar\omega$80 study stresses that iridates are difficult for INS because Ir has a large neutron absorption cross-section, which had made the observation of excitations above 100 meV rare (Fujita et al., 2020). High-pressure INS faces an analogous limitation from pressure-cell attenuation and structured background. The Zr-based bulk metallic glass hybrid clamp cell was designed specifically to address these problems: the Zr$\hbar\omega$81Al$\hbar\omega$82Ni$\hbar\omega$83Cu$\hbar\omega$84 material has reported neutron transmission of about 67% at 10 meV for 1 cm thickness, the empty-cell INS background is broad and featureless rather than sharply phononic, and a reference experiment on CsFeCl$\hbar\omega$85 showed transmission about 2.5 times larger than that of a conventional monobloc CuBe clamp cell at 2.96 meV (Hayashida et al., 16 Feb 2026).

Time resolution is a second frontier. A laser-neutron pump-probe capability at HYSPEC combines a nanosecond pulsed 515 nm laser, synchronized through custom ADCROC electronics, with INS on the quantum magnet Rb$\hbar\omega$86MnF$\hbar\omega$87. The method resolves a laser-induced nonequilibrium structure factor: on the magnon creation side the data match the equilibrium distribution at 3.6 K, but on the annihilation side there is extra intensity inconsistent with simple heating. This establishes that INS can follow nonequilibrium magnon populations rather than only equilibrium detailed balance (Hua et al., 2023).

Other studies extend the range of what may count as an INS excitation. A theoretical treatment of superconducting rings shows that a neutron can scatter inelastically from the magnetic field of a superconducting current, gaining kinetic energy in discrete steps while the number of trapped flux quanta in the ring decreases (Agafonov, 2010). In the Ce pyrochlores, the interpretation of low-energy continua remains an active issue: the reanalysis of Ce$\hbar\omega$88Zr$\hbar\omega$89O$\hbar\omega$90 argues that the low-energy INS spectrum is mostly the electric-monopole continuum rather than the gauge photon (Zhao et al., 2024). This suggests that, even when the experimental observable is fixed, the microscopic meaning of a continuum can remain model-dependent.

Taken together, these developments indicate that INS is no longer restricted to equilibrium single-magnon spectroscopy. Within the present corpus it serves as a direct probe of $\hbar\omega$91, a quantitative route to microscopic Hamiltonians and CEF wavefunctions, a method for extracting thermodynamic quantities from PDOS, a platform for symmetry-resolved and polarization-resolved spectroscopy, a route to nonequilibrium measurements, and a candidate probe of higher-order correlations and chiral phonon dynamics (Fujita et al., 2020, Anand et al., 2017, Shiryaev et al., 2020, Zhang et al., 2013, Hua et al., 2023, Su et al., 2020)