Debye–Waller factor Created 2026-10-14 Updated 2026-11-05
The Debye–Waller factor suppresses coherent scattering because atoms fluctuate about their equilibrium positions. For Gaussian displacement , the amplitude factor is and the intensity factor is its squared modulus. It multiplies the lattice interference factor, so exact extrema of the full finite-lattice intensity need not coincide with the constructive-interference centers when this envelope varies appreciably.
Diffraction intensity 2026-10-25
For coherent elastic diffraction, the intensity associated with a given scattering amplitude is proportional to its squared modulus. Integrated powder diffraction intensities additionally include the powder diffraction multiplicity, Lorentz-polarization factor, Debye–Waller factor, absorption and instrument corrections. A peak's height is not its integrated intensity unless its profile width and shape are also specified.
Lamb-Mössbauer factor 2026-10-18
The Lamb-Mössbauer factor is the fraction in the coherent elastic line. In the harmonic-lattice approximation it is a Debye–Waller factor. Write the displacement of the emitting or absorbing atomic nucleus relative to its equilibrium site as , put , and define the angular wavenumber . Choose the equilibrium site so that . Its directional displacement variance is . A harmonic solid in thermal equilibrium has a Gaussian distribution for this displacement, whose characteristic function gives the elastic amplitude
Taking the squared magnitude gives the intensity fraction
The factor belongs to the amplitude exponent; it disappears when the amplitude is squared. For an isotropic solid, , so an equivalent expression is . Directional and three-dimensional mean squared displacements must not be substituted for one another without this factor.
An Einstein solid approximation replaces the vibrational spectrum by one positive angular frequency . For the directional quantum harmonic oscillator, with the same effective mass as in ,
Here is the Boltzmann constant and is temperature. This follows by taking the expected value of in each oscillator energy eigenstate, , and averaging with the Boltzmann distribution. Substitution gives
The zero-temperature limit retains zero-point energy and a nonzero displacement variance. Cooling therefore increases the recoil-free fraction towards a limit below one for finite . At high temperature, the exponent tends to . This is a single-frequency approximation to a harmonic lattice's coherent elastic weight, not a complete description of all zero-net-energy processes of an isolated oscillator. Real vibrational spectra, directional binding and anharmonic motion can change its predictions.
Extend the stipulated teaching system by choosing , so . Then , while at . If source and absorber are independent, the fraction of events eligible for recoil-free emission followed by recoil-free absorption contains the product ; this gives for two zero-temperature copies. It is not the total probability of absorption, which also depends on the number of resonant atomic nuclei, geometry and spectral line overlap.
Figure 1.
Free recoil separates the lines; lattice binding supplies an elastic fraction
. Original model using only the stipulated teaching parameters. Left: peak-normalized Lorentzian line shapes of the isolated atomic nucleus, separated by ; each physical width is , with an inset resolving one peak. Right: Lamb-Mössbauer factor in the Einstein solid approximation for three choices of . Every zero-temperature intercept is below one. No historical trace is reproduced. Details.
Initially approximate every tungsten atom by the same atomic scattering factor . Pairing the quarter-coordinate sites in the crystal basis structure factor gives
For cubic powder diffraction, increasing deflection corresponds to increasing , since . The first possibilities are
Family | | \ | | \ | | \ | | \ | | \ | |
Hence the first two strong rings have Miller indices
The powder diffraction multiplicity is six for , from three choices of nonzero axis and two signs. For it is , from six permutations and four sign choices. If the atomic scattering factors and angular corrections are taken as approximately equal for these nearby rings, their integrated diffraction intensity obey
This is a multiplicity-based estimate. An exact measured ratio also depends on the momentum-transfer dependence of the atomic scattering factor, the Lorentz-polarization factor, the Debye–Waller factor, absorption and instrument geometry. In particular, the source does not provide the beta-phase lattice constant or instrument corrections needed to calculate those angular factors. Peak heights additionally depend on peak widths; the estimate refers to integrated ring intensities.
The two body-centre/origin sites and the other six sites have different local environments. Their effective atomic scattering factors therefore need not be identical even in an elemental, chemically pure crystal. Write them as and , respectively. The crystal basis structure factor becomes
For the low-angle families this gives
Thus can produce a weak extra powder diffraction ring below . This is an accidental diffraction extinction from equal atomic scattering factors: its cancellation in part d is exact only in the identical-atom-scattering approximation. Small differences in the bonding electron density, or site-dependent Debye–Waller factors, can leave a small nonzero amplitude. For instance, taking and gives while ; the corresponding intensity suppression is quadratic in . No impurity is required. The and cancellations remain for arbitrary scalar in this two-site-class model.
There is a genuine periodicity mismatch in the source. Write . The intended periodic displacements in a harmonic chain satisfy , equivalently for unwrapped positions. The printed instead makes the potential , with equilibrium at equal positions, not . The following lattice derivation uses periodic displacements, as required by the requested position operator.
Put , , and Fourier transform the displacements and momenta with normalization . The spring quadratic form is diagonal, with factor . Thus each nonzero mode is a harmonic oscillator with
Creation and annihilation operators obey . Inverting the oscillator-coordinate expansion gives
The conjugate momentum is . Fourier orthogonality diagonalizes the Hamiltonian as . The normalized oscillator vacuum is annihilated by every ; applying creation operators gives number states, each added phonon contributing . The omitted zero mode is a free centre-of-mass coordinate, which must be fixed or treated separately rather than assigned an oscillator vacuum. Accordingly the reduced displacement commutator is .
For , separate creation and annihilation parts. Their commutator is a scalar, so the supplied Baker--Campbell--Hausdorff formula and vacuum expectations give
Summing the geometric lattice phases yields
The ratio is interpreted by continuity at its removable singularities. The interference factor has its maximum at , where all lattice phases coincide. The source's Fourier-sum hint should read , rather than .
The claim about the full matrix element also needs qualification: its Debye–Waller factor envelope varies with . The absolute global maximum is ; at nonzero reciprocal-lattice points the value is . Moreover the derivative of the envelope there is nonzero, whereas that of the interference factor is zero, so those points are not exact local maxima of the full finite-chain intensity. The reciprocal-lattice condition maximizes the interference factor, which is the usual constructive-scattering interpretation of the requested result.
To determine a crystal structure, scatter probes whose wavelength is comparable to atomic spacings, typically . X-ray diffraction, neutron diffraction and electron diffraction probe periodic order through coherent elastic scattering. With incident and outgoing wavevectors , define . Elastic scattering has and constructive interference occurs when , a reciprocal lattice vector. The equivalent Bragg's law is
where the detector scattering angle is , not . Peak positions reveal plane spacings and the unit cell; systematic absent reflections constrain centring and the atomic basis.
For X-ray diffraction, a laboratory tube or a synchrotron supplies a collimated beam. In a monochromatic single-crystal experiment a goniometer rotates the specimen and an area detector records reflection positions and integrated wave intensities. A powder experiment uses many randomly oriented crystallites; the reflections form cones recorded as rings or peaks versus , permitting indexing and determination of lattice constants without selecting one crystallite. For a cubic cell, for example, . The white-beam Laue method is especially useful for crystal orientation.
The crystal basis structure factor is
after geometric and thermal corrections, including the Debye–Waller factor. For X-rays, describes scattering from the atomic electron density. Thus positions provide lattice geometry while relative wave intensities constrain atomic positions and species. Wave intensities do not directly give the phases of : structure refinement or other constraints must address the phase problem.
Neutron diffraction uses a reactor or spallation source, wavelength selection and a detector bank or time-of-flight detection. Neutrons interact with nuclei through nuclear scattering lengths that vary between elements and isotopes rather than simply with atomic number. This makes them useful for light atoms, isotope contrast and distinguishing some neighbours in the periodic table; their magnetic interaction also probes magnetic order. Electron diffraction in a transmission electron microscope gives strong signals from very small thin crystallites, although strong multiple scattering often complicates quantitative wave intensity analysis. These techniques complement one another.
A phonon dispersion relation requires both momentum and energy transfer, so elastic diffraction alone is insufficient. In coherent inelastic neutron scattering, the measured transfers for phonon creation are
Here is the phonon wavevector reduced into a Brillouin zone, and labels an acoustic phonon branch or optical phonon branch. Negative energy transfer corresponds to annihilating an already occupied phonon. For neutrons , so wavelength and energy selection are linked. Typical thermal-neutron energies of a few to tens of millielectronvolts ( electronvolts) are well matched to many lattice vibrations.
A triple-axis spectrometer uses a monochromator crystal to select , the sample as the second scattering axis, and an analyser crystal to select before detection. Rotate these elements and the specimen to hold fixed while scanning ; inelastic peaks locate . Repeat at different along symmetry directions to map the phonon dispersion relation. The dynamic structure factor governs peak weights, including the selection of phonon polarization by the scattering geometry. A time-of-flight spectrometer instead infers outgoing neutron energy from travel time over a known flight path and records many directions simultaneously.
Inelastic X-ray scattering can also resolve phonons when very high energy resolution is available. Optical Raman scattering and Brillouin scattering measure vibrations near the zone centre because photon momentum transfers are small compared with a typical reciprocal lattice spacing; they complement, rather than generally map the entire zone like, momentum-resolved neutron measurements. Ultrasound measurements give the long-wavelength sound speed, the initial slope of an acoustic phonon branch.