跳到主要内容
Sandbox Physics
星际尘埃、红外光与偏振

星际尘埃 · L2

透射光与热辐射的偏振为什么相差九十度

取向一致的尘粒,对两个电场方向可以吸收得不一样。透过去的光偏向容易幸存的方向;沿同一视线,热辐射偏向容易发射的方向。先把这两个方向定义清楚,再谈偏振角,就能看出九十度从哪里来。

先说朝哪里看,再说偏振朝哪里

选定光的传播方向 n\boldsymbol n,电场就在垂直于它的平面里振动。给这个平面定两根互相垂直的轴 e1,e2\boldsymbol e_1,\boldsymbol e_2。Lab的红外光沿 +z+z 到达观察者,因此天空平面是 x,yx,y。主画面却是 x,zx,z 侧视投影;看侧面方便拖云,但不会改变实际观察者。

许多波的振动有随机相位,不能直接把电场振幅相加后当成强度。将复振幅与共轭相乘后取平均,得到两行两列的相干矩阵 JJ。按这里的归一化,矩阵迹就是强度或光束功率。两个对角元记录两种线偏振模式的功率;非对角元保留偏好方向斜着放时所需的信息。

Stokes II 是总量;QQ 为正,表示第一根轴上的功率较大,为负则第二根较大;UU 区分两个斜对角方向。这里采用实的吸收张量,没有产生圆偏振的机制,因此也不从它推断真实观测中圆偏振应当多大。

I=J11+J22,Q=J11−J22,U=2Re⁡J12I=J_{11}+J_{22},\qquad Q=J_{11}-J_{22},\qquad U=2\operatorname{Re}J_{12}
正负号依赖已声明的天空坐标。线偏振度由Q与U组成的向量长度除以I得到。

对应阅读:Hensley & Draine, arXiv:2208.12365v1 (2022) · Dust absorption, heat and polarization: standalone science

把取向差异写进材料

用单位向量 b\boldsymbol b 指定取向轴。教学材料的不透明度由标量 κ0(ν)\kappa_0(\nu) 乘一个正定张量 KK。对比参数 η\eta 从零到二分之一。沿取向轴的特征值为 1−η1-\eta,两个垂直方向都为 1+η/21+\eta/2。

一束光只用到横向电场,所以先把张量投影到该光束的天空平面。均匀尘埃柱质量为 Σ\Sigma 时,振幅传播矩阵是 A=exp⁡[−κ0ΣK⊥/2]A=\exp[-\kappa_0\Sigma K_\perp/2],输出为 Jout=AJinA†J_{\mathrm{out}}=AJ_{\mathrm{in}}A^\dagger。指数里的二分之一属于振幅;乘上共轭后已得到强度衰减,不能再把强度透射率平方一次。

这套张量规定了两个方向的吸收差异,没有求辐射力矩怎样让尘粒转起来,也没有计算取向随尺寸的变化或磁流体演化。场方向柄改变的是输入轴;柄的长度只是方向引导,不表示磁场强度。

K=(1+η2)1−3η2bbT,tr⁡K=3K=\left(1+\frac\eta2\right)\mathbf1-\frac{3\eta}{2}\boldsymbol b\boldsymbol b^{\mathsf T},\qquad\operatorname{tr}K=3
做角平均时,这个固定迹给出原来的标量冷却不透明度;沿某个特定方向的发射仍然可以各向异性。

对应阅读:Dust absorption, heat and polarization: standalone science · Hensley & Draine, arXiv:2208.12365v1 (2022)

算一块均匀薄板:幸存较多与发射较多

把取向轴放在第一根天空轴上,设 η=0.35\eta=0.35,两个模式的不透明度就是 0.65κ00.65\kappa_0 与 1.175κ01.175\kappa_0。用总强度为一的非偏振光照一块均匀板,每个模式初始分到一半。若参考光深 κ0Σ=1\kappa_0\Sigma=1,输出分别为 e−0.65/2e^{-0.65}/2 和 e−1.175/2e^{-1.175}/2。

相加得 I=0.415432I=0.415432,相减再除以总量得 Q/I=0.256632Q/I=0.256632。第一种模式吸收少,因此透射光偏向第一根轴。这是可独立验算的均匀板极限,不是对某颗恒星的拟合。

再沿同一方向看同样材料的光学薄热辐射。模式功率正比于各自的吸收能力,而不是透过去剩下的比例,于是 Q/I=(0.65−1.175)/(0.65+1.175)=−0.287671Q/I=(0.65-1.175)/(0.65+1.175)=-0.287671。负号表示热辐射的电场轴转向第二根轴,与透射偏振垂直;总强度仍然为正。

整团云需要把不同温度、不同权重的单元加起来,照明光线也未必指向观察者。均匀板例子解释局部符号,并给程序一个解析极限;它不能授权我们把不同光束平面里的角度直接拿来相减。

(QI)trans=tanh⁡ ⁣[τ2−τ12]\left(\frac QI\right)_{\mathrm{trans}}=\tanh\!\left[\frac{\tau_2-\tau_1}{2}\right]
均匀、对角吸收板,入射非偏振,第一种模式的光深较小。
透射偏向吸收较少的模式,薄热发射偏向垂直模式。
独立双模式解析板。热辐射横线指另一个满足光学薄条件的发射波长,不是按横轴透射光深求出的厚板热辐射解。

对应阅读:Dust absorption, heat and polarization: standalone science

为什么转九十度,两个Stokes量都反号

线偏振度为 p=Q2+U2/Ip=\sqrt{Q^2+U^2}/I。只要它非零,电场位置角可写成 χ=12atan2⁡(U,Q)\chi=\tfrac12\operatorname{atan2}(U,Q)。二分之一来自轴的对称性:一条线转过 180∘180^\circ,仍是同一条偏振轴。

固定天空坐标,把整套物理配置绕观察方向主动旋转 α\alpha,Stokes向量就转过 2α2\alpha。转 90∘90^\circ 后,QQ 与 UU 同时反号,图像相应旋转,而这套整体旋转的源—云问题中 II 不变。科学复算检验的是整个配置;只转固定云里的取向轴,可能改变吸光,不能要求总强度也不动。

把 b\boldsymbol b 反向,bbT\boldsymbol b\boldsymbol b^{\mathsf T} 不变。如果正沿取向轴看过去,横向两个特征值相等,偏振消失。偏振度为零时,位置角没有定义,图像中的方向棒应该消失,而不是显示一个看似精确的角度。

(Q′U′)=(cos⁡2α−sin⁡2αsin⁡2αcos⁡2α)(QU)\begin{pmatrix}Q'\\U'\end{pmatrix}=\begin{pmatrix}\cos2\alpha&-\sin2\alpha\\\sin2\alpha&\cos2\alpha\end{pmatrix}\begin{pmatrix}Q\\U\end{pmatrix}
这里是在固定天空轴中旋转物理配置;若改成旋转坐标轴,转动方向相反。

对应阅读:Dust absorption, heat and polarization: standalone science

一个偏振角,为什么还画不出三维磁场

偏振度小,可能是取向轴接近视线,也可能是排列较弱、同一波束内不同方向相互抵消,或者混进了非偏振辐射。教学云只用一个规定轴;真实云不一定这么简单。一个偏振度不能把这些原因全部区分开。

原生Astrodust模型包含随波长变化的材料响应、尺寸分布与拟合取向函数。其偏振截面带符号,也不必遵循这里的简单正对比张量。教学推导与完整材料预测,应该分别看。

在Lab中,可以先选沿视线的取向,再改成斜向。读出实际Stokes结果,而不要把拖动柄当成反演出来的场。真实反演还需要额外的几何与取向假设,并说明它们带来的不确定性。

对应阅读:Hensley & Draine, arXiv:2208.12365v1 (2022) · Dust absorption, heat and polarization: standalone science

练习与解释性答案

  1. 均匀板有 τ1=0.65\tau_1=0.65 与 τ2=1.175\tau_2=1.175。透射电场偏向哪根轴,为什么热辐射符号相反?

    展开解释

    透射偏向第一根轴,Q/I=+0.256632Q/I=+0.256632,因为该模式被移走得较少。沿同一方向的光学薄热发射按基尔霍夫定律偏向较容易吸收的模式,所以同一坐标中 Q/I=−0.287671Q/I=-0.287671。这是两条垂直的偏振轴,不是负的总功率。

  2. 红外偏振消失了,能否断定磁场为零?

    展开解释

    不能。即使在这个规定取向的简单模型中,沿取向轴观察也会没有偏振。真实数据还可能有弱排列、分辨不开的方向抵消和非偏振成分。这个模型根本不测磁场强度。

动手检验预测

沿着同一团尘埃追踪吸光、升温与偏振

  1. 选择斜向取向预设,在保存的热辐射Stokes数据里找到非零U,并联系规定的天空轴解释符号。
  2. 把倾角转到沿视线,确认偏振降至零,此时图像不应给出有意义的电场方向。
  3. 固定云,拖动场方向柄,同时比较热图和截获功率。说明这与把源、云和取向整体旋转有什么不同。

来源与计算范围

  1. Hensley & Draine, arXiv:2208.12365v1 (2022)

    固定首稿,第2.1至2.3、4至6节,式1至11/27与图6/11/12/18。材料截面、温度分布和规定的尺寸相关取向,与教学张量分别处理。2022年展望不作为当前状态。

  2. Dust absorption, heat and polarization: standalone science

    MIT自有科学代码与另行许可的CC0原生输入。离线实际重算十步,保留粗网格失败、独立体积/普朗克积分参照及可选首次联网真实练习。站点补充说明记录后续成功获取,原冻结包保持不变。

反馈这个页面有错误、看不明白,或者缺了什么?告诉我们
反馈类型

提交时会附带本页地址、Lab ID、网站语言和设备类型。请不要填写个人信息。反馈最多保存 12 个月。隐私政策