黑洞照片是怎样拍的?干涉成像技术简介

上一篇文章“黑洞都可以拍!现在的望远镜为什么那么牛?” 介绍了望远镜的最基本指标。

这里,我们详细介绍一下现代的射电望远镜拍摄技术:甚长基线干涉测量 Very-long-baseline interferometry (VLBI) 的基础。需要用到的数学工具是:矢量;部分积分和复数的概念。需要更多解释的请让我知道 :-)

图1:甚长基线干涉测量 [Very-long-baseline interferometry (VLBI)]。遥远的星体上某个点发出的电磁波信号以相同方向抵达相距一定距离的电磁波接收器(望远镜),对积累记录下的信号进行关联处理可以得到遥远星体的合成影像。

为了简单起见,我们只考虑相对位置是\vec{L}的两个射电望远镜,且它们只接收遥远的星体上发出来的波长为 \lambda 的电磁波信号。由于星体到地球的距离远远大于\vec{L}的大小,星体上某个点发出的电磁波信号会以相同的方向(表示为单位矢量 \hat{k})分别进入这两个望远镜①和②。

只要测量了电磁波的入射方向\hat{k},就能知道星体的(黑白)影像信息;想象一下,如果我们知道星星的指向(和\hat{k}平行),那是不是就可以自己画在纸上,然后告诉别人这就是你“拍”的星图?偷偷告诉你,射电望远镜基本上就是这么干的...

两个望远镜天线在某个时间t接收的电磁波信号分别可以写成
E_{1} = A \cos (\omega t + \delta) , E_{2} = A \cos (\omega t +\delta+ \alpha + \phi )。 其中,\omega是电磁波的角频率,和波长\lambda成反比关系;简单起见我们假设两个望远镜的接收强度(例如最大信号电压)A 都一样。\delta是由星体发出的电磁波的随机初始相位而决定的量。\alpha是可以用仪器人为地单独调整的相位(例如可以通过改变望远镜到信号存储点的距离而实现)。

利用上一篇讨论的,可以知道两个望远镜接收的电磁波有相位差 \phi = \vec{k}\cdot \vec{L} = \frac{2\pi}{\lambda} (\vec{L}\cdot \hat{k}) = 2\pi\frac{l} {\lambda} ,---公式(1)其中\vec{k} = \frac{2\pi}{\lambda} \hat{k}是电磁波的波矢量,波矢的大小|\vec{k}| = \frac{2\pi}{\lambda}表示每传播一个波长\lambda的距离,电磁波的相位就改变2\pi(360度)回到原来的状态(例如,就像水波的一个波峰移到下一个波峰)。而 l=\vec{L}\cdot \hat{k}是相对位置\vec{L}在电磁波传播方向的投影大小(见图1)。知道相位差\phi,就知道星体的方位信息。

要测量\phi,我们对两个电磁波信号进行关联干涉。一种类似于讨论水波衍射时的方法是,测量两个信号叠加后的强度:
信号相加的测量= \langle (E_1+ E_2)^2 \rangle = \frac{1}{T}\int_0^{T} (E_1+ E_2)^2 dt ,其中尖括号代表测量时对叠加的信号进行时间平均(就是说,测量某个很短时间内接收的电磁波能量)。因为是时间平均,随着时间t的快变振荡项可以被平均掉,我们得到:
信号相加的测量= \langle (E_1+ E_2)^2 \rangle \approx 2A^2 \cos(\phi/2+\alpha/2)^2 。--公式(2)\vec{L}的投影l改变\lambda/2这么多的时候,信号相加的测量结果就从最大值变成最小值。越大的距离\vec{L},就可以有越大的分辨率。

不过,虽然这个信号相加的方法很直观,但它有个实际应用上的问题:它对信号噪声敏感。因此,目前的射电望远镜用的是另外一种信号处理方法:信号相乘。这是
信号相乘的测量= \langle E_1 E_2 \rangle 。类似地,信号处理器过程中会丢掉时间t的快变振荡项,所以有:
信号相乘的测量项= \langle E_1E_2 \rangle \approx \frac{1}{2}A^2 \cos(\phi+\alpha) 。 --公式(3)如果信号带有噪声E_1 \rightarrow E_1 + (随机的噪声项) 。 只要噪声在两个望远镜上是独立无关联的、且平均值为零(这是噪声通常的特点),那么可以发现,噪声对我们的结果公式(3)没有影响。这种巧妙的方法其实在很多地方都有应用。

从公式(3)可以提取出相位信息\phi,进而利用已经知道的距离\vec{L}得到天空的指向\hat{k}的信息。通过测量不同的\vec{L}下的结果(例如选用不同地方的望远镜,等目标星体在天空不同位置的时候再测量),就能得出\hat{k}的方向(精度受限于前述的分辨率)。

明白了?那我们看复杂一点的实际情况。黑洞照片是个圈圈而不是一个点,不可能只有一个方向\hat{k}的信息就能给出一张照片图案。实际上,望远镜接收到的信号不是前面讲的那么简单,而是包括很多很多不同方向位置的电磁波信号:
E_{1} = \sum_{\hat{k}} A_\hat{k} \cos (\omega t + \delta_{\hat{k}}) , E_{2} = \sum_{\hat{k}} A_\hat{k} \cos (\omega t + \delta_{\hat{k}}+ \alpha+\phi_{\hat{k}})。其中,求和号\sum_{\hat{k}}表示后面的项在不同方向\hat{k} 的情况下进行相加,下标带方向\hat{k}的表明这些量的值和方向\hat{k}有关。 \delta_{\hat{k}}是不同方向的随机初始相位。对黑洞这类遥远的天体,不同方向\hat{k}的随机相位\delta_{\hat{k}}是没有关联的。\alpha就像前面讲的,是一个可以人为控制的参数。相位差\phi_\hat{k} = \vec{k}\cdot \vec{L}和之前的一样,只是加了个下标。我们看信号相乘方法的测量结果。平均掉与时间有关的快速振荡项和与随机相位\delta_{\hat{k}}有关的涨落项,可以得到,对于相对位置为\vec{L} 的关联信号是S(\vec{L})= \langle E_1E_2 \rangle \approx \sum_{\hat{k}} I(\hat{k}) \cos(\frac{2\pi}{\lambda}\hat{k}\cdot\vec{L}+\alpha) , --公式(4) 其中星体的发光强度随空间指向的分布I(\hat{k})=\frac{1}{2}A_{\hat{k}}^2就是我们需要“冲印”的照片。目标星体在不同时间出现在天空不同的位置,因为目标星体的平均指向是照片的参考点,所以公式(4)里面的相对位置\vec{L}的指向,会随着星体在天空不同的位置而相对地变化。另外,利用多个不同相对位置\vec{L}的望远镜,我们可以获得很多不一样的\vec{L}下的信号S(\vec{L})I(\hat{k})就可以从信号S(\vec{L})反推出来(具体是逆傅里叶变换,后面会详细解释)。

公式(4)可以用合奏音乐来理解:把\vec{L}想象成演奏的时间点或者曲谱,I(\hat{k}) \cos(\frac{2\pi}{\lambda}\hat{k}\cdot\vec{L}+\alpha)则看成具有不同音调\hat{k}的乐器,这样S(\vec{L})就是所有乐器合奏的音乐。通过听不同曲谱\vec{L}下的演奏,音乐家就能确定里面有什么样的乐器或者I(\hat{k})

下面看具体的计算。为了表达成与星体照片相关的结果,我们以星体的大概指向\hat{k}_c 为基准,把单位矢量\hat{k}写成
\hat{k} = \hat{k}_c- \vec{r} ,其中矢量\vec{r}是相对基准的偏移,并定义\tilde{I}(\vec{r})=I(\hat{k})。分别测量\alpha=0\alpha = \pi/2的结果,可以从S(\vec{L})定义它的复数形式 \tilde{S}(\vec{L})= e^{i \frac{2\pi}{\lambda} \hat{k}_c\cdot \vec{L} } (\langle E_1E_2 \rangle_{\alpha=0} + i \langle E_1E_2 \rangle_{\alpha=\pi/2} ) 。把求和写成傅里叶积分形式:\tilde{S}(\vec{L}) = \int \tilde{I}(\vec{r}) e^{i \frac{2\pi}{\lambda} \vec{r}\cdot \vec{L} } d\vec{r} 。 其逆变换就是图像信号:\tilde{I}(\vec{r}) \propto \int \tilde{S}(\vec{L}) e^{-i \frac{2\pi}{\lambda} \vec{r}\cdot \vec{L} } d\vec{L} 。 ---公式(5)要得到更全面的图像\tilde{I}(\vec{r}),如前面说的,需要用多个望远镜在不同位置不同相对角度测量。由于\vec{L}不能任意取值(例如它的大小目前没有大于地球尺寸),所以目前"拍到“的图像\tilde{I}(\vec{r})不可避免会有点模糊。

什么,没明白?没事,黑房里拿个甜甜圈用黄灯照着,把照相机故意调失焦,然后拍一下。。。是不是和人类第一组黑洞照片很像呢?

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 214,588评论 6 496
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 91,456评论 3 389
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 160,146评论 0 350
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 57,387评论 1 288
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 66,481评论 6 386
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 50,510评论 1 293
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,522评论 3 414
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,296评论 0 270
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,745评论 1 307
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,039评论 2 330
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,202评论 1 343
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,901评论 5 338
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,538评论 3 322
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,165评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,415评论 1 268
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,081评论 2 365
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,085评论 2 352

推荐阅读更多精彩内容