---
title: "Faraday 深度谱该怎样读"
description: "沿固定版本论文读 RM synthesis，计算频带限制，再把等离子体透镜与散射模型接到可检验的观测。"
canonical_url: "https://sandboxphysics.com/zh/stories/interstellar-radio/faraday-depth/"
language: "zh-CN"
page_type: "story"
alternate_languages:
  en: "https://sandboxphysics.com/stories/interstellar-radio/faraday-depth/"
last_modified: "2026-10-28"
---

[追踪穿过等离子体的射电脉冲](https://sandboxphysics.com/zh/stories/interstellar-radio/)

AC06 · L4

# Faraday 深度谱该怎样读

Sandbox Physics · 2026-10-28

打开 Brentjens 与 de Bruyn 的《Faraday Rotation Measure Synthesis》，arXiv:astro-ph/0507349v2，2005 年 11 月 2 日第二版。依次读第 2、4、6 节，再对照 Schnitzeler 与 Lee 的 arXiv:1411.1412v1。读深度谱时，始终把采样响应放在旁边。

## 开始前需要什么

完成偏振矢量一章。这里的傅里叶变换，是把许多试探深度下的旋转矢量相加；卷积，则表示真实分布被已知响应模糊。

## 第 2 节先分清深度与斜率

先看式 3–9。Faraday 深度，是从某个发射位置到观测者之间，电子加权的磁场积分。一个背景源隔着一块旋转屏时，可以对应单个深度，偏振角也有单个 RM 斜率。若一路都在发射、也都在旋转，不同位置发出的波就会积累不同深度。

函数 F 把复偏振发射分配给不同 Faraday 深度。观测到的复偏振，是这些分量各自转过与波长平方相关的角度以后，再相加。深度是磁积分，不是几何距离。两个相隔很远的位置可以有相同深度；磁场反转，还会让深度沿路不再单调。

望远镜只测有限的一些正波长平方。在波长空间乘上这个采样窗口，反推的深度分布就会被模糊。式 36–38 给出离散重建与旋转量传递函数 RMTF。解释峰或缺口之前，先按真实权重计算这条响应。

$$
P(\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 表示每个 Faraday 深度区间的偏振发射。平均参考波长改变复相位，不会提高深度分辨率。恢复的分布被 R 模糊。

对应阅读：Brentjens & de Bruyn (2005), Faraday Rotation Measure Synthesis, arXiv v2 · Schnitzeler & Lee (2014), Rotation measure synthesis revisited, arXiv v1

## 图 7 与图 8 展示采样带来的多解

图 7 比较两组各八个采样点。波长平方的总跨度相近，旁瓣却明显不同。总跨度决定大致分辨率，采样点具体落在哪些位置，还会决定混叠。图 8 调整了这些间隔，压低一个共振。

实测谱里的第二个峰，可能是第二个分量，也可能是第一个峰的旁瓣。把响应平移到主峰上叠加比较，再检查 Q/U 残差，把原来的通道响应纳入候选模型。反卷积可以让表示更尖锐，采用的假设却不会增加实际频率覆盖。

对应阅读：Brentjens & de Bruyn (2005), Faraday Rotation Measure Synthesis, arXiv v2 · Schnitzeler & Lee (2014), Rotation measure synthesis revisited, arXiv v1

## 第 6 节的三个频带尺度

图 11 与式 63–66 分开总波长平方跨度、最短波长，以及每个通道的宽度。它们分别影响深度分辨率、对宽深度结构的灵敏度、通道内退偏振。观测设计里应把三种尺度分开计算。

理想矩形覆盖从四百到八百兆赫时，波长平方两端约为 0.562 与 0.140 平方米。按论文的近似式，深度分辨率为 8.22 弧度每平方米，最大可恢复深度尺度约 22.4。这些估计采用论文的矩形权重。教学探测器采样稀疏，有自己的传递函数，确切峰宽与混叠要按这些采样点另算。

一个发射区即使看起来没有被主峰分辨开，也可能已经丢掉很多偏振信号。固定频带时增加通道数，有助于采样和通道退偏振；若要看很宽的 Faraday 结构，仍需获得更短波长的观测。

$$
\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}
$$

这些是固定论文在所述响应下的设计近似值，实际分辨率和旁瓣应由测得的 RMTF 给出。

对应阅读：Brentjens & de Bruyn (2005), Faraday Rotation Measure Synthesis, arXiv v2 · Schnitzeler & Lee (2014), Rotation measure synthesis revisited, arXiv v1

## 再读 2014 年论文的通道响应

转到 Schnitzeler 与 Lee 第 2 节式 4–7。一个通道平均一段偏振矢量，响应的积分为一。论文的归一化反旋转矢量，与前章单分量最小二乘振幅拟合，都考虑了有限通道，但振幅归一方法不同。课程数值例子采用前章推导的单分量最小二乘方法。

借用 sinc 校正前，先比较频率矩形和波长平方矩形。矢量在通道内转得很快时，两者差别会变明显。拟合应让 Q、U 通过实际响应；严重平均过的通道，需要整个频段的响应来描述。

对应阅读：Schnitzeler & Lee (2014), Rotation measure synthesis revisited, arXiv v1 · Brentjens & de Bruyn (2005), Faraday Rotation Measure Synthesis, arXiv v2

## 透镜与散射尾各问什么

接着读 Cordes 等人的 arXiv:1703.06580v1，2017 年 3 月 20 日第一版，第 2 节与图 1–4。高斯等离子体透镜是一个具体模型。图 1 显示多条射线可以到达同一观测者；增益图预测，源或观测者穿过透镜图样时，不同频率会明暗起伏。这些是模型预测，没有把课程的短 CHIME 切片认作透镜事件。

随机密度起伏则引出额外时延分布与散射核。Geyer 与 Karastergiou 的 2.1 节解释了指数尾及其几何前提。若要推湍流谱，需要跨频率观测、源尺寸和屏几何，还要考虑各向异性与有限屏。只拟合出一条负四次方长尾，确定不了湍流谱。

继续观察时，可以问得很具体：反复爆发的谱特征会不会漂移？计入仪器展宽以后，长尾还在不在？一个几何能不能同时解释计时和频谱结构？这些问题帮助比较传播模型；射电源的发射机制留给独立的源物理路线。

对应阅读：Cordes et al. (2017), Lensing of FRBs by plasma structures, arXiv v1 · Geyer & Karastergiou (2016), revisiting pulsar scattering

## 从教学拟合走向标定观测

Pandhi 等人的 arXiv:2401.17378v2，2024 年 5 月 2 日第二版，从 CHIME/FRB 基带电压出发。读 I.1、II.1–II.2 节，了解需要哪些 Stokes 信息，以及怎样校准。教程中的首目录强度图，是另一种数据产品。

读下一份偏振结果时，先记坐标和正负约定、通道响应、校准方式与前景扣除。把采样响应放在恢复的深度谱旁边，在 Q/U 空间比较候选分量，最后再单独处理观测 RM 与本地参考系的换算。

对应阅读：Pandhi et al. (2024), CHIME/FRB polarimetry, arXiv v2 · Brentjens & de Bruyn (2005), Faraday Rotation Measure Synthesis, arXiv v2

## 动手检验预测

### [先找采样混叠，再解释第二块云](https://sandboxphysics.com/zh/labs/interstellar-radio-propagation/)

1. 导出默认合成样本。比较旋转量次峰与已知只有一个分量的注入模型。
2. 选粗通道，观察 Q/U 与拟合谱。解释哪部分损失来自通道平均，哪部分多解来自采样间的空隙。

[上一篇：离线重跑一次真实 CHIME 射电闪光](https://sandboxphysics.com/zh/stories/interstellar-radio/chime-intensity/)

## 来源与计算范围

1. [Pandhi et al. (2024), CHIME/FRB polarimetry, arXiv v2](https://arxiv.org/html/2401.17378v2)

   第二版，2024 年 5 月 2 日。I.1、II.1–II.2 节：Stokes 数据、红移权重与校准。式 2 的说明把 kpc 与系数 0.812 并列，课程采用 NRAO 量纲一致的 pc 归一。

2. [Cordes et al. (2017), Lensing of FRBs by plasma structures, arXiv v1](https://arxiv.org/pdf/1703.06580v1)

   第一版，2017 年 3 月 20 日。第 2 节、式 1–4 与图 1–4，按高斯透镜模型预测阅读，没有实施透镜反演或源机制推断。

3. [Brentjens & de Bruyn (2005), Faraday Rotation Measure Synthesis, arXiv v2](https://arxiv.org/html/astro-ph/0507349v2)

   第二版，2005 年 11 月 2 日。第 2、4、6 节，式 3–9、63–66。核读 Faraday 深度、传递函数、角度混叠与频带限制。课程做单分量数值实验，并未复现论文的天空观测。

4. [Schnitzeler & Lee (2014), Rotation measure synthesis revisited, arXiv v1](https://arxiv.org/html/1411.1412v1)

   第一版，2014 年 11 月 5 日。第 2 节式 1–6 及第 3 节：有限频率通道与波长平方上的矩形响应不同。实验直接积分矩形频率通道。

5. [Geyer & Karastergiou (2016), revisiting pulsar scattering](https://academic.oup.com/mnras/article/462/3/2587/2589386)

   2.1 节，式 1–3：单边指数响应采用无限大、各向同性薄屏与高斯角分布。有限屏和各向异性会改变响应。