Fourier Transform µSR

Jess H. Brewer

Canadian Inst. for Advanced Research
and Dept. of Physics & Astronomy,
Univ. of British Columbia, Vancouver, BC, Canada V6T 1Z1


This page is under construction at     /~jess/musr/
The Fourier transform (and its cousin, Maximum Entropy analysis) is the most tantalizing and frustrating technique in µSR, because it can never deliver what one wants from it - namely an undistorted picture of the local field distribution in the sample. We take our data (usually) in the time domain and often we are actually interested in the frequency domain; but you can't get there from here without a model, and FFT (the Fast Fourier Transform, a numerical method virtually always used to transform discrete time spectra into the Fourier-complement frequency spectra), like MaxEnt, is intrinsically model independent. More precisely (and this is the crux of the problem) it does employ a model in which every frequency component is presumed to be a completely independent parameter. Although this may sound like just what we want, in fact it is not. We do not really believe that two frequency bins separated by 0.1 MHz can be wildly different for a 10 µs time spectrum, but that is exactly what FFT assumes.

  1.   Apodization vs. Ringing:
  2.   The most annoying and irreducible example of the consequences of this defective model is the ``ringing'' due to different frequencies matching or mismatching the endpoint of a finite time spectrum. To combat this effect, one can use the technique of apodization, in which the time spectrum is multiplied by an ``envelope'' function (a sort of artifical relaxation function) which smoothly reduces the oscillatory signal to zero at the end of the time range. This artificially broadens the resulting frequency spectrum, of course,

    [From David Spencer's Ph.D. thesis, ~ 1980]

    Apodization is especially needed when, through a flaw in the counting electronics or just sheer bad luck in counting statistics, one time bin has an anomalously high or low value.

    [From Tanya Riseman's Ph.D. thesis, 1993]

    In this case the frequency spectrum will exhibit a huge oscillation over its entire range (think what a single ``spike'' in a frequency spectrum means about the time spectrum, and then reason backwards).

    [From Tanya Riseman's Ph.D. thesis, 1993]

  3.   Complex Amplitude:
  4.   Another aspect of the FFT that is often overlooked is the fact that its output is a mixture of a real part that can be quite sharp since it is the cosine function (which is an even function and therefore doesn't mind starting abruptly at t=0) and an imaginary (sine function) part which does mind, and which has an enormous extra broadening because it is the transform of the Heaviside function. If the two are simply added in quadrature (or squared to produce the power spectrum) this superfluous broadening wipes out any fine structure in the frequency spectrum. At high fields it is no simple matter to separate the real and imaginary parts; this can be done numerically, but only to limited precision, leaving a spectrum with unknown (and probably unknowable) distortions.

  5.   Qualitatively Useful Anyway!
  6.   Many other problems are encountered by the determined Fourier transformer seeking an aesthetically pleasing frequency spectrum; most of these can be overcome sufficiently to satisfy the observer, as illustrated here, but one should never forget that all model-independent transformations from the time domain into the frequency domain involve arbitrary and fundamentally unreliable quantitative decisions by the transformer, and should therefore be regarded as qualitative illustrations, never as reliable quantitative measurements.

  7.   Paramagnetic Systems:   where FT-µSR is indispensible!
  8.   Actually the first (and most unambiguously useful) applications of FT-µSR were to complex spin systems in which the muon polarization oscillates at many frequencies simultaneously. Of the many examples available, I chose the classic study of Mu* in semiconductors as an illustration. Whereas ESR (or EPR if you prefer) usually probes the spin Hamiltonian by sweeping the magnetic field through the resonances on the Breit-Rabi diagram at fixed RF frequency, FT-µSR reveals all the frequencies at once for a given magnetic field.
      In gallium arsenide (like many other semiconductors, fortuitously excluding Si, Ge and diamond) the nuclei have their own HF interaction with the muonium electron, broadening the lines so dramatically at low field that they cannot be observed in TF-µSR. Fortunately, at high field (in this case 1.2 T) the nuclear hyperfine (NHF) splittings are reduced to the point where the muonium HF structure can be seen clearly.
      These formulae show how the muon polarization evolves in strong transverse field (TF) and longitudinal field (LF) assuming (often incorrectly) a scalar (isotropic) Mu HF interaction.
      In most semiconductors and insulators muonium experiences an anisotropic spin hamiltonian.
      Another important use of FT-µSR is in the spectroscopy of radicals. Rather than show one of the many examples from muonium chemistry, I have chosen the rather exotic case of cis-polyacetylene, a solid composed of the cis isomer of polyacetylene (the trans isomer forms an organic semiconductor). Notice the partially resolved NHF structure in the solid; in liquids, these splittings tend to be averaged away by rapid rotation and tumbling of the molecules.

Jess H. Brewer
Last modified: Wed Sep 2 11:13:46 PDT 1998