Skip to main content
Sandbox Physics
Follow a radio pulse across the plasma

AC06 · L4

Read a Faraday spectrum without turning it into a distance map

Open Brentjens and de Bruyn’s “Faraday Rotation Measure Synthesis”, arXiv:astro-ph/0507349v2, dated 2 November 2005. Follow Sections 2, 4 and 6, then compare the finite-channel treatment in Schnitzeler and Lee, arXiv:1411.1412v1. The goal is to read a Faraday spectrum with its response function beside it.

Section 2 distinguishes depth from a fitted slope

Start with equations 3–9. Faraday depth is an electron-weighted magnetic integral from one emitting location to the observer. A single background source behind one screen has one depth, and its polarization angle has one RM slope. If emission and rotation occur together, radiation from different places accumulates different depths.

The function F assigns complex polarized emission to each Faraday depth. The observed complex polarization sums those contributions after each has rotated by an amount set by wavelength squared. Depth is a magnetic integral, not a geometric coordinate. Two separate places can have the same depth, and a reversing field can make depth vary non-monotonically along a path.

The telescope samples only a limited set of positive wavelength-squared values. Multiplying by that sampling window in wavelength space blurs the inferred depth distribution. Equations 36–38 give the discrete reconstruction and the rotation-measure transfer function, RMTF. Compute that response from the actual weights before interpreting a peak or a gap.

P(λ2)=∫F(ϕ)e2iϕλ2dϕ,R(ϕ)=∑jwje−2iϕ(λj2−λ02)∑jwjP(\lambda^2)=\int F(\phi)e^{2i\phi\lambda^2}d\phi,\qquad R(\phi)=\frac{\sum_jw_j e^{-2i\phi(\lambda_j^2-\lambda_0^2)}}{\sum_jw_j}
F carries polarized emission per Faraday-depth interval. The mean reference wavelength shifts complex phase without improving depth resolution. The recovered distribution is blurred by R.

Reading: Brentjens & de Bruyn (2005), Faraday Rotation Measure Synthesis, arXiv v2 · Schnitzeler & Lee (2014), Rotation measure synthesis revisited, arXiv v1

Figures 7 and 8 reveal aliases from the sampling

In Figure 7 compare two sets of eight sampling points. Their full wavelength-squared span is similar, yet their sidelobes differ strongly. Resolution depends on the broad span; ambiguity also depends on where the individual samples sit. Figure 8 changes those intervals and suppresses a resonance.

A second peak in a measured spectrum could be a second component or a sidelobe of the first. Overlay the response shifted to the strongest peak, inspect Q and U residuals, and test the candidate model through the original channel responses. Deconvolution can sharpen a representation, but its assumptions do not create new frequency coverage.

Reading: Brentjens & de Bruyn (2005), Faraday Rotation Measure Synthesis, arXiv v2 · Schnitzeler & Lee (2014), Rotation measure synthesis revisited, arXiv v1

Section 6 gives three different band limits

Figure 11 and equations 63–66 separate the full wavelength-squared span, the narrowest sampled wavelength, and each channel width. They control depth resolution, sensitivity to broad depth structure, and within-channel depolarization respectively. They are different scales and should appear as separate numbers in an observing plan.

For an ideal rectangular coverage from four hundred to eight hundred megahertz, the endpoints in wavelength squared are about 0.562 and 0.140 square meters. The paper’s approximate resolution is 8.22 radians per square meter, and the largest recoverable depth scale is about 22.4. Those estimates assume its rectangular weighting. The teaching detector has sparse frequency samples and its own transfer function; its exact peak width and aliases must be computed from those samples.

A broad emitting region can lose much of its polarized signal even if it looks unresolved by the main peak. Increasing the number of channels at fixed coverage helps averaging and depolarization, but does not substitute for access to shorter wavelengths when broad Faraday structure is the target.

δϕ≃23Δλ2,ϕlargest scale≃πλmin⁡2,∣ϕmax⁡∣≃3δλ2\delta\phi\simeq\frac{2\sqrt3}{\Delta\lambda^2},\qquad \phi_{\rm largest\ scale}\simeq\frac{\pi}{\lambda_{\min}^2},\qquad |\phi_{\max}|\simeq\frac{\sqrt3}{\delta\lambda^2}
These are approximate design estimates from the pinned paper, for its stated response. A measured RMTF supplies the actual resolution and sidelobes.

Reading: Brentjens & de Bruyn (2005), Faraday Rotation Measure Synthesis, arXiv v2 · Schnitzeler & Lee (2014), Rotation measure synthesis revisited, arXiv v1

Read the channel response in the 2014 paper

Turn to Section 2, equations 4–7 of Schnitzeler and Lee. Each channel averages a vector over a response whose integral is one. The normalized derotation-vector construction in that paper and the one-component least-squares amplitude fit in our chapter are related response-aware methods, but they use different amplitude normalization. Our numerical example is not an implementation of every estimator in the paper.

Compare the frequency top-hat with the wavelength-squared top-hat before borrowing a sinc correction. The discrepancy grows when a vector turns rapidly inside a channel. Fit measured Q and U through the full response of a severely averaged channel.

Reading: Schnitzeler & Lee (2014), Rotation measure synthesis revisited, arXiv v1 · Brentjens & de Bruyn (2005), Faraday Rotation Measure Synthesis, arXiv v2

A lens and a scattering tail ask different questions

Next open Cordes and colleagues, arXiv:1703.06580v1, 20 March 2017, Section 2 and Figures 1–4. Their Gaussian plasma lens is a specific model. Figure 1 shows that several ray paths can reach one observer. The gain maps predict frequency-dependent brightening and dimming as an observer or source moves across the lens pattern. These are model predictions, not an identified lens in our short CHIME cutout.

Random density fluctuations instead motivate a distribution of extra delays and a scattering kernel. Geyer and Karastergiou’s Section 2.1 explains the exponential case and the geometric assumptions behind it. Inferring a turbulence spectrum requires observations across frequency, a source-size and screen-geometry model, and attention to anisotropy or finite screens. A fitted inverse-fourth-power tail alone does not determine a turbulence spectrum.

Useful follow-up questions are concrete: do spectral features drift over repeated bursts, does a tail remain after instrumental smearing is included, and can one geometry explain timing and frequency structure together? These guide model comparison without choosing the radio-emission mechanism, which belongs to the separate source-physics route.

Reading: Cordes et al. (2017), Lensing of FRBs by plasma structures, arXiv v1 · Geyer & Karastergiou (2016), revisiting pulsar scattering

Move from the teaching fit to calibrated observations

Pandhi and colleagues’ arXiv:2401.17378v2, dated 2 May 2024, starts from CHIME/FRB baseband voltage data. Read Sections I.1 and II.1–II.2 for the required Stokes information and calibration. The first-catalog intensity picture used in our tutorial is a different data product.

For a new polarization result, note the coordinate and sign convention, channel response, calibration and foreground subtraction. Plot the sampling response beside the recovered depth spectrum, test alternative components in Q/U space, and keep any rest-frame conversion separate from the observed RM. This sequence is the reusable outcome of the reading guide.

Reading: Pandhi et al. (2024), CHIME/FRB polarimetry, arXiv v2 · Brentjens & de Bruyn (2005), Faraday Rotation Measure Synthesis, arXiv v2

Try the prediction

Find a sampling alias before naming a second cloud

  1. Export a default synthetic record. Compare secondary RM peaks with the known single-component injection.
  2. Choose coarse channels and inspect both Q/U and the fitted spectrum. Explain which lost information comes from channel averaging and which ambiguity comes from the gaps between samples.

Sources and calculation scope

  1. Pandhi et al. (2024), CHIME/FRB polarimetry, arXiv v2

    Version 2, 2 May 2024. Sections I.1, II.1–II.2: Stokes data, redshift weighting and calibration. The text of equation 2 lists kpc alongside coefficient 0.812; this course uses the dimensionally consistent pc normalization from NRAO.

  2. Cordes et al. (2017), Lensing of FRBs by plasma structures, arXiv v1

    Version 1, 20 March 2017; Section 2, equations 1–4 and Figures 1–4. Read as a Gaussian-lens prediction; no lens inversion or source-mechanism inference is claimed.

  3. Brentjens & de Bruyn (2005), Faraday Rotation Measure Synthesis, arXiv v2

    Version 2, 2 November 2005. Sections 2, 4 and 6; equations 3–9 and 63–66. Read for Faraday depth, the transfer function, angle aliases and observing-band limits; this is a reading guide and a one-component numerical experiment, not a reproduction of the paper’s sky observations.

  4. Schnitzeler & Lee (2014), Rotation measure synthesis revisited, arXiv v1

    Version 1, 5 November 2014. Section 2, equations 1–6, and Section 3: a finite frequency response differs from a rectangular response in wavelength squared. The simulation integrates rectangular frequency channels.

  5. Geyer & Karastergiou (2016), revisiting pulsar scattering

    Section 2.1, equations 1–3: the one-sided exponential assumes an infinite isotropic thin screen with a Gaussian angular distribution. Finite and anisotropic screens alter the response.

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