-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathintroduction.tex
569 lines (454 loc) · 42.1 KB
/
introduction.tex
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
\chapter{绪论} \label{chapter:intro}
\epigraph{... Onde possiamo stimare che de stelle innumerabili sono altre tante lune, altre tanti globi terrestri, altre tanti mondi simili a questo.}{\textit{Giordano Bruno}}
\section{引言}
\subsection{地球与人类,无独有偶?}
由古至今,人类在宇宙中的所处的位置就一直魂牵梦绕着人们的思绪。这个难题既
包括了人类居住的地球在空间上究竟是否独特,另一方面也包含人类(或生命)在宇
宙中的唯一性与否。而正因人类拥有了实证与思考的能力,因而在这个问题上的观测
与探索也从未曾停歇。
首先,行星科学(Planetology)是一门古老的以观测为基础的科学。早在公元前,古
希腊等各大文明已有天文学家(如伊巴谷)统计了相当可观的夜空发光天体的运行与
方位\cite{Xuan1992, wikipedia_astronomer}。而在以托勒密为代表所提倡的
地心说主导了这个星球一千多年后,文艺复兴时期的哥白尼等人也通过观测提出了日
心说\cite{Copernicus1543}。自此地球在空间位置谜团便开始层层被揭开,
地球的面貌(图 \ref{fig:rwalone})也愈发完整地展现在了人类面前。
\afterpage{
\begin{savenotes}
\begin{figure}[ht]
\centering
\begin{subfigure}[b]{.45\textwidth}
\captionsetup{width=0.9\textwidth}
\centering
\includegraphics[width=0.95\textwidth]{figures/chapter1/fig1a_Apollo17_BlueMarble.jpg}
\caption[Apollo 17 号船员拍摄的南半地球,拍摄时间 1972 年 12 月 7 日。(图片版权:NASA)]{Apollo 17 号船员拍摄的南半地球\footnotemark[1],拍摄日期为 1972 年 12 月 7 日。(图片版权:NASA)}
\label{fig:apollo17earth}
\setcounter{footnote}{1}
\addtocounter{footnote}{+0}\footnotetext{\url{https://www.nasa.gov/content/blue-marble-image-of-the-earth-from-apollo-17}}
\end{subfigure}
\begin{subfigure}[b]{.45\textwidth}
\captionsetup{width=0.9\textwidth}
\centering
\includegraphics[width=0.95\textwidth]{figures/chapter1/fig1b_Voyager_PaleBlueDot.jpg}
\caption[Voyager 1 于 1996 年 9 月 12 日拍摄到的多色合成地球肖像照。(图片版权:NASA/JPL-Caltech)]{Voyager 1 于 1996 年 9 月 12 日拍摄到的多色合成地球肖像照\footnotemark[2]。(图片版权:NASA/JPL-Caltech)}
\label{fig:voyagerearth}
\addtocounter{footnote}{+1}\footnotetext{\url{http://photojournal.jpl.nasa.gov/catalog/PIA00452}}
\end{subfigure}
\caption{人类从太空看到的赖以居住和生存的母亲行星 --- 地球。\textbf{(a)} Apollo 17 号飞船航向月球的途中宇航员们回望并拍摄到的地球。这张图被称作「蓝宝石」,图中可以清晰的看到南半球的大气云层与被冰雪覆盖的南极点,另外还可看到非洲与亚洲大陆。\textbf{(b)} Voyager 1 航天器在距离太阳四十亿公里外拍摄到的地球。在最右侧的太阳散射光柱中,地球就只是不到一个像素的小白点,而在这颗蔚蓝的星球上曾居住并流逝着所有的人类足迹。 }
\label{fig:rwalone}
\end{figure}
\end{savenotes}
}
此外,人们往往会从哲学角度来思考地球与生命在宇宙中普遍与否。比如原子论的宇宙
描述早已开端于古希腊与古印度。随后,意大利哲学家 Bruno 亦于 1584 年提出无穷宇
宙论:太阳与地球在宇宙中并不特殊,不计其数的其他星体也和地球相差无几(文献\citen{Bruno1584})。
%\begin{quote}
%Onde possiamo stimare che de stelle innumerabili sono altre tante lune, altre tanti globi terrestri, altre tanti mondi simili a questo.
%\end{quote}
历史一脉相承至今,随着人们对太阳系内行星与地球的运动探求愈来愈充分,也许便不难
理解搜索系外行星与系外生命为何如此举足轻重。而且从技术层面上,人类也已经有足够
的能力将「第二颗地球」与「地外生命」放入可见的未来计划中\cite{WoolfAngel1998}。
\subsection{系外行星的定义}
二十一世纪以前,由于对冥王星的测量并不充分,科学界从传统层面上也似乎没有必要争
论如何去界定行星。尔后,随着更多的海外天体(Trans-Neptunian Objects,简称 TNOs)
被发现,冥王星是否还是太阳系第九大行星的争辩也愈演愈烈\cite{SternMitton2005}。
终于,国际天文学联合会(International Astronomical Union,简称 IAU)于 2006 年在
「B5 决议\footnotemark[3]」中给太阳系内的行星做出如下定义:
\setcounter{footnote}{3}
\footnotetext{\url{https://www.iau.org/static/resolutions/Resolution_GA26-5-6.pdf}}
\begin{enumerate}[leftmargin=1\parindent]
\item 必须拥有围绕太阳的轨道;
\item 质量应充分大,使得其自引力足够克服刚性力,以形成(近)球外形;
\item 另外,行星必须清空轨道临近的其它天体。
\end{enumerate}
根据此定义,冥王星不再属于行星范畴,而被称作一颗矮行星。那如今观测到的近 3500 颗
\footnote{来源 NASA Exoplanet Archive \url{http://exoplanetarchive.ipac.caltech.edu/index.html},
数据库更新至 2017 年 2 月 14 日}系外行星呢?从词源学上,行星(planet)一词取名自希腊语,
那么系外行星(exoplanet)自然也该遵循上述 2006-B5 决议。然而部分系外气态巨行星
(gas giant)质量已经大于十个木星质量($M_\tif{J}$),此类行星的质量上限在 IAU 决议中
也并未给出。实际上,早于冥王星争辩,IAU 系外行星工作小组就已规定真实质量小于 13
$M_\tif{J}$(或行星内部氘元素热核反应的质量下限\cite{Baraffe2003})的行星才被称作系外
行星\footnote{请查看链接 \url{http://home.dtm.ciw.edu/users/boss/IAU/div3/wgesp/}}。在本文
亦一律遵循以上两者来定义系外行星。
\subsection{系外行星扼要史}
人类从诞生之始就默默地注视着夜空的「神行者」--- 行星。而作为专注于研究系外行星系统
的科学(Exoplanetology),其研究对象则跳脱出传统可观测的太阳系各大行星,转向了银河
系内太阳周围临近的恒星。回顾历史,这门领域的兴起也并非一蹴即成,而是融入许多先人不
停探索、试错与更新的艰辛历程。
近代系外行星搜索开始于十九世纪中页,Jacob 在对 70 Ophiuchi 的天体测量数据中找到类似
行星的信号\cite{Jacob1855},然而 Moulton 后续观测立即证否了这颗行星\cite{Moulton1899}。
在随后的一个世纪内,对木星质量的系外行星探索逐渐变得活跃,如 van de Kamp 于 1963 年
将仪器系统信号误认为一颗围绕巴纳德恒星的行星(文献 \citen{vandeKamp1963})。
另外一面,Struve 在1952 年提出短周期类木系外行星可同时通过视向速度与测光探测的设想
\cite{Struve1952}。然而科学界当时理所应当都认为木星只能处在距离主星 5 AU 轨道上,因此
需要太长的时间跨度来确认探测\footnote{一般来说,由于引力所造成的周期信号被连续确认三
个周期以上便被认为是真实的天体信号。}。不过,恒星光谱学依然在其间飞速发展,如 HD 114762
被发现拥有一颗非常接近系外行星质量上限的褐矮伴星\cite{Latham1989}。
此时,来自美国的二人团队 Wolszczan 和 Frail 于 1992 年利用射电脉冲信号的周期变化率在
脉冲星 PSR1257+12 周围率先成功探测到包含两颗约 3 $M_\oplus$ 行星的行星系统
\cite{WolszczanFrail1992}。而此时 Walker 等人则在对主序恒星长达 12 年的视向速度监测中,
宣布没有迹象表明他们嗅探到系外木星信号\cite{Walkeretal1995}。
同一年,来自瑞士日内瓦天文台的 Mayor 和 Queloz 做出了系外行星领域里程碑式的发现(文
献\citen{MayorQueloz1995}): 他们利用坐落在法国 Haute-Provence 天文台的 ELODIE 光谱
仪从 150 颗恒星样本中成功探测到第一颗类太阳恒星周围的行星 --- 51 Peg b(见图
\ref{fig:51pegsi})。这颗轨道周期只有 4.23 天的类木行星于同一年被美国 Marcy 和 Butler 团
组证实,随后便引起了对此类行星(即热木星)形成理论机制广泛而又深远的讨论
(\S \ref{sec:pltfrmatntheory}),一直延续至今。
\begin{figure}[ht]
\centering
\includegraphics[width=0.95\textwidth]{figures/chapter1/fig2_51pegsi.jpg}
\caption[1995 年,日内瓦天文台天文学家 Mayor 和 Queloz 连续观测恒星 51 Pegsi 得到的视向速度曲线图,观测点已全数叠加至其伴星的周期 4.23 天。该系统主星是一颗类太阳恒星,因此视向速度曲线振幅所对应的最小质量为 0.47 $M_\tif{J}$ ,系人类首次发现到类太阳恒星周围的短周期类木行星。]{日内瓦天文台天文学家 Mayor 和 Queloz 连续观测恒星 51 Pegsi 得到的视向速度曲线图,观测点已全数叠加至其伴星的周期 4.23 天。该系统主星是一颗类太阳恒星,因此视向速度曲线振幅所对应的最小质量为 0.47 $M_\tif{J}$ ,系人类首次发现到类太阳恒星周围的短周期类木行星。图片取自他们于 1995 年发表在自然杂志的文献 \citen{MayorQueloz1995}。}
\label{fig:51pegsi}
\end{figure}
紧随其后的二十年内,系外行星领域各大新奇发现此起彼伏。图 \ref{fig:pltdiscyear} 为年度新
发现系外行星数目柱状图,该数目明显呈现出指数式增长,尤其是 2009 年 $Kepler$ 太空望远
镜\cite{Boruckietal2010} 的发射新带来了大批通过凌星法发现的系外行星。图\ref{fig:pltdiscyear}
中所展现的愈加多样化的样本也极大地丰富了我们对系外行星的认知。
\begin{figure}[hb]
\centering
\includegraphics[width=0.95\textwidth]{figures/chapter1/fig3_discoveryear.pdf}
\caption{系外行星发现数量趋势图,不同颜色代表不同的探测方法(参照 \S \ref{sec:detmeth}),可以看出该数目呈现指数式增长。(图片取自 Open Exoplanet Catalogue)}
\label{fig:pltdiscyear}
\end{figure}
\section{探测方法} \label{sec:detmeth}
天文学所研究的天体普遍离地球遥远,因而观测手段也主要集中在分析天体发射的光子上。
观测系外行星相较恒星或太阳系内行星,通常需要用到非常高分辨率、灵敏度和稳定性的仪器,
也因而会遇到诸多困难,甚至还得排除恒星自身活动的干扰。比如测量围绕类太阳的系外地球
需要精度为 1 m$\odot$s$^{-1}$的视向速度测量,相当于使用分辨率为 R = 100,000 的光谱仪器探测
主星光谱约 $10^{-6} \AA$ 的移动大小。参考书籍 \citen{Perryman2014exohb},本文大
致罗列目前主要探测系外行星的方法如下:
\subsection{视向速度}
假使恒星周围存在行星,那么它们便会同时绕着公共的质心(Center Of Mass,简称为 COM)做
开普勒运动。在观测中,这种三维运动可分解为视线平面内的二维运动与视线方向上的运动。视
向速度(Radial Velocity,简称 RV)法就是通过测量视线方向上恒星谱线的多普勒红(蓝)移来
探测周围行星的存在。
在简单的二体运动中(\S \ref{apdx:twobodyproblem}),可以得到一颗轨道半长径为 $a$,质量
为 $M_\tif{p}$ 的行星,可对质量为 $M_\tif{s}$ 的宿主恒星造成如下大小的半振幅视向速度 $K_1$
(单位:米每秒):
\begin{equation} \label{eq:rvk}
K_1 \simeq \frac{28.4 }{\sqrt{1-e^2}} \frac{M_\tif{p}\sin i}{M_\tif{J}} \, \bigg(\frac{M_\tif{p}+M_\tif{s}}{M_\odot}\bigg)^{-1/2} \big(\frac{a}{1 \tif{AU}}\big)^{-1/2} \ \ ,
\end{equation} \myequation{探测系外行星中视向速度法引起的主星速度振幅}
其中各物理量的意义参见附录\ref{apdx:twobodyproblem}。由上式可知,视向速度方法最重要的缺
陷即它只能测量(或拟合)出行星的最小质量(minimum mass,即$M_\tif{p}\sin i$),因为行星的轨
道倾角 $i$ 或视线方向夹角无法被测量。但与此同时,视向速度最明显优势是它可以精确测量行星
的轨道偏心率 $e$(图 \ref{fig:rvcurve} )。
\begin{figure}[h!]
\centering
\includegraphics[width=1.0\textwidth]{figures/chapter1/fig5_rvcurve.pdf}
\caption[行星对主星造成的视向速度曲线示意图,从左至右分别代表不同的行星偏心率和轨道近心点取向。此图版权归 Perryman M. 所有。]{行星对主星造成的视向速度曲线示意图,分别代表不同的行星偏心率和轨道近心点取向。从左至右分别示意系统 HD 73256\cite{Udry2003HD73256},HD 142022\cite{Eggenberger2006HD142022} 和 HD 4113\cite{Tamuz2008HD4113}。此图取自文献 \citen{Perryman2014}。}
\label{fig:rvcurve}
\end{figure}
另外由于视向速度需要通过尽可能多的主星谱线来确认谱线的位移,因而往往更容易探测到 FGKM
型主序星附近的大质量行星。在这里值得一提的是,由 Murphy 等人提出的激光频率梳法\cite{Murphyetal2007lasercomb}
与 Molecule Iodine Cell 或 ThAr Lamps 等传统的光谱仪器定标方法相比,可以更稳定重复的覆盖同
样的谱线区间,因此未来在视向速度方法上有相当可观的应用前景。
%优缺点,国际主流望远镜和团队?
\subsection{凌星法} \label{sec:transit}
从原理上来说,凌星法是通过测量行星遮挡主星时主星亮度(流量)的变化来探测行星的。如图
\ref{fig:transit} 所示,系外行星处于恒星与地球的连线时被称作凌星( transit)或主掩食(primary
eclipse),而当恒星处于系外行星与地球的连线时则被成为掩星(occultation)或次掩食
(secondary eclipse)\footnote{凌星被约定特指半径较小天体遮挡大天体,而掩星则相反。双星中
的掩食通常只用于大小相近天体的互相遮挡}。
\begin{figure}[ht]
\centering
\includegraphics[width=0.95\textwidth,angle=-90, trim={0 1cm 0 0}, scale=0.9]{figures/chapter1/fig6_transit.pdf}
\caption[凌星法探测行星示意图,上半部分表示不同的轨道时刻。下方则为整个系统观测到的流量变化。版权 Joshua Winn。]{凌星法探测行星示意图,凌星法探测行星示意图,上半部分表示不同的轨道时刻。下方则为整个系统观测到的流量变化。此图取自书籍 \citen{Seager2010exobook}。}
\label{fig:transit}
\end{figure}
事实上,掩食在太阳系内已属屡见,比如日、月食和金星凌日现象。虽然凌星法探测与视向速度
法在原理上被 Struve 于同一年提出\cite{Struve1952},但真正意义上的第一次观测到系外凌星却
比视向速度晚了五年\cite{Seager2010exobook}。当时视向速度已经发现了 64 颗系外行星,在对
其中 6 颗行星系统的凌星监测中,行星 HD 209458 b 的主掩食成功被观测到\cite{Henryetal2000, Charbonneauetal2000}。
这种滞后主要是因为凌星发生的几何空间概率较低($p = R_\tif{s} / a$),与此同时还要求观测的
相对测光精度至少达到 5 \textperthousand(对应于木星凌太阳)。凌星法也因此通常结合巡天计
划(尤其是大视场巡天)展开,因而提高测光精度、处理数据和去除假阳性信号也成了该方法的技
术难关\footnote{凌星法只能得到系外行星的候选体,一般还需额外使用其他方法联合认证该行星,
如凌星计时变化(Transit Timing Variation,简称 TTV)或者 RV。}。目前地面上比较成熟的巡天项
目,包括 The Trans-atlantic Exoplanet Survey(TrES)\cite{Alonsoetal2004TrES},XO
\cite{McCulloughetal2005XO},The Hungarian-made Automated Telescope Network(HAT)
\cite{Bakosetal2007HAT} 以及 Wide Angle Search for Planets(SuperWASP)
\cite{Pollaccoetal2006WASP},相应的空间巡天项目也有 Convection, Rotation and planetary
Transits(CoRoT)\cite{Bargeetal2008CoRoT} 与 $Kepler$\cite{Boruckietal2010}。
\subsection{天体测量}
\begin{figure}[h!]
\centering
\includegraphics[width=0.95\textwidth]{figures/chapter1/fig7_astrometry.pdf}
\caption[由于行星的存在,宿主恒星的天球坐标随着时间的变化示意图,版权所有:Perryman M。]{由于行星的存在,宿主恒星的天球坐标随着时间的变化示意图。此图取自文献 \citen{Perryman2014}。}
\label{fig:astrometry}
\end{figure}
和视向速度探测主星前后摆动不同,天体测量法(astrometry)着重于观测恒星在天球上的左右摆
动。行星对宿主恒星造成的来回振幅大小可用如下公式估算
\begin{equation} \label{eq:astrometry}
\alpha \simeq \bigg(\frac{M_\tif{p}}{M_\tif{s}}\bigg) \, \big(\frac{a}{1 \tif{AU}}\big) \, \bigg(\frac{d}{1 \tif{pc}}\bigg)^{-1} \tif{arcsec} \ \ .
\end{equation} \myequation{天体测量法探测系外行星}
如图 \ref{fig:astrometry} 所示,天体测量法只需要恒星在天球球面上的经度纬度信息就可探测到
行星,它也对主星的性质没有任何依赖,因而有其独特的优势,如可探测轨道倾角参数 $i$。
2013 年,Gaia 天体测量卫星成功发射,虽然目前来看科学数据尚未达到探测系外行星的精度
\cite{GaiaCo2016},期待今后几年会有振奋人心的成果。
\subsection{直接成像} \label{sec:drctimgmeth}
中国古话有道是「耳听为虚,眼见为实」,直接成像法便可把行星直观地展示出来。技术上这其实
并非易事:距太阳系几十个秒差距以外的行星与其主星的空间张角(angular separation)往往远
小于望远镜分辨率。即使空间上能分辨,恒星的黑体辐射也往往比行星高近十个量级(图
\ref{fig:contrast})。在实际的观测中,对行星直接成像往往需要星冕仪(coronagraph)或干涉仪
(interferometer)的辅助(如图 \ref{fig:hr8799})。
\begin{figure}[h]
\centering
\begin{subfigure}[b]{.48\textwidth}
\captionsetup{width=0.9\textwidth}
\centering
\includegraphics[width=0.95\textwidth]{figures/chapter1/fig8a_dicontrast.jpg}
\caption[仪器接受类太阳恒星在 10 pc 以外的黑体辐射(黄色实线),其余为太阳系行星与热木星在此恒星周围黑体发射与发射的叠加流量。在可见光波段行星与主星的对比度可差至十个量级,因此直接成像法一般选择在长波范围(如红外波段)探测系外行星。图片版权 Seager and Deming。]{可见光波段行星与主星的对比度可差至十个量级,因而直接成像法一般选择在长波范围(如红外波段)观测系外行星。图片摘自文献\citen{Seager2010}。}
\label{fig:contrast}
\end{subfigure}
\begin{subfigure}[b]{.48\textwidth}
\captionsetup{width=0.9\textwidth}
\centering
\includegraphics[width=0.95\textwidth]{figures/chapter1/fig8b_hr8799.jpg}
\caption[Keck 望远镜利用角向较差成像与自适应光学技术在近红外波段通过直接成像得到的 HR 8799 系统。中心恒星的大部分光已被星冕仪遮挡,留下小部分噪声,与远处四颗行星 。图片版权:Marois 等人。]{Keck 望远镜在近红外波段通过直接成像得到的 HR 8799 系统。中心恒星的大部分光已被星冕仪遮挡,留下小部分噪声,与远处四颗行星 。图片取自文献\citen{Marois2010}。}
\label{fig:hr8799}
\end{subfigure}
\caption{直接成像法探测系外行星。\textbf{(a)} 仪器接受类太阳恒星在 10 pc 以外的黑体辐射强度,其余为太阳系行星与热木星在此恒星周围黑体发射与发射的叠加流量。\textbf{(b)} Keck 望远镜利用角向较差与自适应光学技术直接成像观测到的 HR 8799 系统。}
\label{fig:directimage}
\end{figure}
直接成像方法往往偏向于探测年轻的系统,如 HR 8799\cite{Marois2008HR8799},
Fomalhaut\cite{Kalas2008} 和 $\beta$ Pictoris\cite{Lagrange2010},这是因为行星在其形成早期会
拥有较强的近红外辐射。目前用于此方法的望远镜有 Hubble Space Telescope(HST),Keck,
Very Large Telescope(VLT)和 Gemini Planet Imager(GPI)等
\cite{Marois2008HR8799,Chauvin2005,Macintosh2014}。
\subsection{其他方法}
\textbf{微引力透镜(microlensing):} 本方法原理可以追溯到 1936 年,Einstein 发表了一篇
计算前景星对视线方向上背景星引力放大率的文章\cite{Einstein1936}。而且在第一颗用微引力透镜
观测得的行星前,Mao 和 Paczynski 就已经从理论上提出了行星能在恒星引力透镜基础上造成额外
的微引力透镜效\cite{Mao1991}。如今已经有 42 个通过微引力透镜发现的系外行星系统,包括 2 个
双行星系统\footnote{数据截至 2017 年 2 月,参见 NASA Exoplanet Archive:\url{http://exoplanetarchive.ipac.caltech.edu/}}。
在役的仪器主要有 Optical Gravitational Lensing Experiment(OGLE)\cite{Udalskietal2002OGLE}
和 Microlensing Observations in Astrophysics(MOA)\cite{Bond2004MOA},它们一般会监测恒星
密度较大的区域如银河系核球(Galactic Bulge,简称 GB)。
\textbf{计时法(timing):} 传统计时法最典型的例子就是最早被发现的系统 --- 脉冲星
PSR1257+12 行星系统\cite{WolszczanFrail1992}。除了传统的 Pulsar Timing Variation 之外,可以
说任何理论上拥有稳定周期信号的恒星计时变化都可以用来探测潜在的系外行星:如凌星计时变化
TTV\cite{Ford2011TTV,Xie2013} 和星震时变\cite{Silvotti2007,Murphy2016}。此方法应用前景也非
常可观,特别是针对多行星系统动力学特征刻画\cite{HolmanMurray2005}。
除此以外,Perryman 列出了更为系统且详细的探测方法,对应的探测能力以及它们在将来可能的应
用\cite{Perryman2000},请参见下页图 \ref{fig:detmethod}。
\begin{sidewaysfigure}
\centering
\includegraphics[scale=0.95, angle=-90]{figures/chapter1/fig4_fmethods.pdf}
\caption[系外行星探测方法大全,图片版权 Michael Perryman。]{系外行星探测方法大全,关于图片说明请参见文献\citen{Perryman2000},此图摘录自该文献的第一幅插图。}
\label{fig:detmethod}
\end{sidewaysfigure}
\section{行星形成理论} \label{sec:pltfrmatntheory}
星际空间并不是空无一物的,英国天文学家 William Herchel 在 18 世纪末观测到恒星周围冷暗物质
吸收带,我们的原太阳就是在这般空间环境中孕育而成\cite{Spitzer1978}。随着原初分子云坍缩,
核心区域会演化成原恒星,而周围的物质会由于原初角动量守恒而沉降成盘状结构以及双极喷流
(bipolar jets),比如观测到的 Herbig-Haro 型天体\cite{Lequeux2005}。
在恒星形成的早期,现有的观测证据主要集中在 T-Tauri 型天体上(又叫 Young Stellar Objects,
简称 YSOs)。YSOs 通常会按照红外谱指数($\alpha_\tif{IR} = \d\log \lambda F_\lambda)/(\d\log \lambda$)
分为 Class 0,Class \RN{1}, Class \RN{2}, Class \RN{3} 四种类型\cite{Andreetal2000}(见图
\ref{fig:ysostage}),它们恰恰各自对应着不同的原恒星演化阶段。
\begin{figure}[t]
\centering
\includegraphics[width=0.95\textwidth]{figures/chapter1/fig9_ysostages.pdf}
\caption[YSOs 按照红外谱指数划分的四种主要类别,版权所有:Armitage P.。]{YSOs 按照红外谱指数划分的四种主要类别,图中背离黑体谱的红外突起部分又被称作红外超,它也是恒星是否存在盘的重要判据之一。此图取自书籍 \citen{Armitage2010}。}
\label{fig:ysostage}
\end{figure}
从某种程度上说,此刻的原恒星盘(protostellar disk)已可被称作原行星盘
\footnote{一般而言在恒星形成的早期,当恒星处于活跃吸积阶段时星周盘被称作原恒星盘,
而在随后的行星形成阶段则被称作原行星盘。}(ProtoPlanetary Disk,简称 PPD)。
根据文献\citen{Armitage2007},由于恒星对周围气体的吸积\cite{Pringle1981}、恒星的高能辐射
\cite{Johnstone1998,Hollenbach1994},以及气体盘自身的粘滞性
\cite{Shakura1973,Gammie1996},角动量会从从气体盘的内侧向外转移,原行星盘也会
逐渐耗散(图 \ref{fig:discevo})。2001 年 Haisch 等人则通过统计不同星团中的红外超恒星比
例\cite{Haisch2001},发现原行星盘的存活时标在不到百万年(约 10 Myr)的量级(见图 \ref{fig:pftimescale})。
这对类木行星行形有着非常重要的限制,因为如今普遍认为木星必须在气体盘消散前吸收足够的气体
以长成如今的质量大小\cite{WilliamsCieza2011}。关于星周盘的详细介绍请参考 \S \ref{chapter:form_evo}。
\begin{figure}[t]
\centering
\includegraphics[width=1.0\textwidth]{figures/chapter1/fig10_discevo.eps}
\caption[原行星盘演化中由 (a) 至 (d) 四个主要过程,分别为活跃吸积阶段、固体沉降阶段、盘消散阶段与残骸盘阶段。图片版权 Williams 和 Cieza。]{原行星盘演化中由 (a) 至 (d) 四个主要过程,分别为活跃吸积阶段、固体沉降阶段、盘消散阶段与残骸盘阶段。图片取自文献\citen{WilliamsCieza2011}。 }
\label{fig:discevo}
\end{figure}
\begin{figure}[t]
\centering
\includegraphics[width=0.9\textwidth]{figures/chapter1/fig11_pftimescale.eps}
\caption[不同星团红外超等效的原(恒)行星盘站所有星团成员的比例。此图可推断星团年龄老于10 Myr 后,恒星周围的气体盘已近乎完全消散。此图版权所有:Hern{\'a}ndez J. 等人。]{不同星团红外超等效的原(恒)行星盘站所有星团成员的比例。此图可推断星团年龄在 5 Myr 左右,恒星周围的气体盘已近乎消散。此图取自文献 \citen{Hernandez2007}。}
\label{fig:pftimescale}
\end{figure}
\subsection{基于太阳系的经典形成理论} \label{sec:clspftheory}
早在 18 世纪中页,关于行星形成假说就已纷繁多样。其中康德、拉普拉斯等人提出的「星云假说」
因可较好应用至太阳系行星系统中而发展壮大。表\ref{tbl:solarsystem} 给出了太阳系八大行星的物
理参数,由此表可粗略看出太阳系各大行星具有非常好的一致共面性,除了水星以外都处于近圆轨道,
且在火星和木星之间存在明显的物态性质差异。如果将每个行星的重元素平均到打散至行星之间的
空隙中,并且混入适量的氢、氦使得新的混合物拥有与太阳相等的金属丰度
\footnote{在天文学中,金属元素特指氢和氦以外的所有重元素。}([Fe/H]),这样得到的最小质量
原太阳行星盘被称作 Minimum Solar Mass Nebula\cite{Weidenschilling1977,Hayashi1981},
MMSN 气体盘的面密度拥有如下的幂率函数形式:
\begin{equation} \label{eq:mmsng}
\Sigma_g = 1.7 \,\times10^3 \, \bigg(\frac{r}{\tif{AU}}\bigg)^{-3/2} \tif{g} \,\tif{cm}^{-2} \ \ ;
\end{equation} \myequation{最小质量太阳星云气体面密度分布表达式}
相应的固体盘的面密度分布同样可描述为:
\begin{equation} \label{eq:mmsns}
\Sigma_s = 7.1 \, f_\tif{ice} \, \bigg(\frac{r}{\tif{AU}}\bigg)^{-3/2} \tif{g} \,\tif{cm}^{-2} \ \ ,
\end{equation} \myequation{最小质量太阳星云固体面密度分布表达式}
其中 $f_\tif{ice}$ 为雪线因子。在距离太阳 2.7 AU 以外,水分子等会凝结成冰雪等固体物质,因而
此因子会从 1 增加至约 4.2\cite{IdaLin2004}。从此最小太阳星云模型出发,假若在垂直方向上给它
一个标高使其变为三维盘,这就是经典太阳系行星形成模型的开端。根据文献
\citen{Dullemond2005},重元素物质(以微米到厘米级尘埃形式存在)会最先开始在垂直方向上沉
降到盘的中平面;紧接着固体物质(dust)开始碰撞结合
\cite{ppvibook2014,Weidenschilling1997,BlumWurm2008}或引力坍缩
\cite{Safronov1972,GoldreichWard1973,YoudinShu2002,ChiangYoudin2010}成星子(planetesimal)。
\begin{table}[t]
\centering
\caption{太阳系八大行星物理参数与轨道参数,轨道根数取值依照 J2000 平赤道参考系,数据源自 NASA/JPL 网站。}
\label{tbl:solarsystem}
\begin{tabular}{cccccccc}
\hline \hline
& $a$(AU) & $P$(days) & $M_\tif{p}$($M_\oplus$) & $e$ & $i$(deg) & $R_\tif{p}$($R_\oplus$) & $\rho$(g cm$^{-3}$) \\ \hline
水星 & 0.3871 & 88.0 & 0.0553 & 0.2056 & 7.00 & 0.383 & 5.427 \\
金星 & 0.7233 & 224.7 & 0.815 & 0.0068 & 3.39 & 0.949 & 5.243 \\
地球 & 1.000 & 365.2 & 1 & 0.0167 & 0.00 & 1 & 5.514 \\
火星 & 1.524 & 687.0 & 0.107 & 0.0934 & 1.85 & 0.532 & 3.933 \\
木星 & 5.203 & 4331 & 317.8 & 0.0484 & 1.30 & 11.2 & 1.326 \\
土星 & 9.537 & 10,747 & 95.2 & 0.0539 & 2.49 & 9.45 & 0.687 \\
天王星 & 19.19 & 30,589 & 14.5 & 0.0473 & 0.77 & 4.01 & 1.271 \\
海王星 & 30.07 & 59,800 & 17.1 & 0.0086 & 1.77 & 3.88 & 1.638 \\
\hline \hline
\end{tabular}
\end{table}
在 MMSN 模型下,星子的半径大概为公里量级。这样大小的星子和尘埃性质不同,它可以脱离气体
的阻尼力并通过自引力保持其结构。此时星子相互之间的引力作用变得显著,它们会经过雪崩增长与
寡头生长两个时期成长为原行星\cite{Greenberg1978,Kokubo1996,Rafikov2003,Lissauer1993}。原
行星(protoplanet,亦称做行星胚胎)会在接下来分别演化成气态巨行星和类地行星,其中形成气态巨行星一般被
认为原行星核通过吸积气体(core accretion)而形成\cite{Mizuno1980,BodenheimerPollack1986,Pollack1996},
而类地行星则需要经历更多的动力学作用过程才能形成\cite{Chambers1998}。
当然,行星形成也还存在很多问题与挑战,解释太阳系太行星还存在许多难关
\cite{Walsh2011,Morbidelli2005},数值模拟如 Nice 模型\cite{Gomes2005Nice}也并不能解释太阳系
细节构型(尤其是水的来源\cite{Raymond2009}),本文仅能作粗略介绍,图 \ref{fig:pfstage} 为上
述经典形成模型的核心各阶段过程的总结。
\begin{figure}[h!]
\centering
\includegraphics[width=1.0\textwidth]{figures/chapter1/fig12_pfsequence.pdf}
\caption[基于太阳系的传统行星形成理论模型在不同阶段的说明图,其中「米级障碍」属于至今为止的重大难题。图片版权 Michael Perryman。]{基于太阳系的传统行星形成理论模型在不同阶段的说明图,其中「米级障碍」属于至今为止的重大难题。此图取自文献 \citen{Perryman2014}。}
\label{fig:pfstage}
\end{figure}
\subsection{经典理论新挑战:系外行星} \label{sec: exopftheory}
由于观测极限,太阳系行星如今尚身处观测能力范围以外。如图\ref{fig:exomassper} 所示,
系外行星种类遍布多样,堪称百花齐放,且已确认数量仍在日趋增加。抛开观测选择
偏差(observational selection bias),太阳系与之相比似乎显得有些「格格不入」。 按
照图\ref{fig:exomassper},本文将系外行星按照其集中区域大致分为三类: 热木星族
(椭圆实线),冷木星族(方框实线)和超级地球族(椭圆虚线)。
\begin{figure}[t]
\centering
\includegraphics[width=1.0\textwidth]{figures/chapter1/fig13_exompplot.pdf}
\caption{现今探测到系外行星的周期 --- 质量散点图。其中不同颜色分别代表不同的探测方法:青色代表视向速度法,灰色是凌星法,黄色指代微引力透镜法,红色代表直接成像法,紫色则表示其他方法。黑色字母 Me,V,E,Ma,J,S,U,N 则分别表示太阳系内水星,金星,地球,火星,木星,土星,天王星和海王星。本文将系外行星进行人为类别划分包络在封闭曲线内,分别为热木星族(椭圆实线),冷木星族(方框实线)和超级地球族(椭圆虚线)。}
\label{fig:exomassper}
\end{figure}
\subsubsection{热木星族 Hot Jupiter Population}
当 Mayor 和 Queloz 与 1995 年发现第一颗围绕类太阳的系外行星 51 Peg b 时
\cite{MayorQueloz1995},整个行星学界都为之震惊,因为这颗类木星的轨道周期
只有 4.23 天。51 Peg b 单独用经典行星形成理论根本无法解释。由于此类行星非
常靠近其主星(轨道周期 $P \le10 $天),质量大于土星质量($0.3\,\tif{M}_\tif{J}$)故
而得名曰热木星。作为本册论文重要的研究对象,此类行星详细的介绍内容请
参见 \S \ref{chapter:data_stat}。
\subsubsection{冷木星族 Cold/Normal Jupiter Population}
与轨道距离主星较近的热木星相对应有一类行星被称作冷木星(或者常规类木星)。
此类行星在观测上拥有长于约一百天的轨道周期,有效温度也比热木星低得多,太阳
系木星就属于一颗典型冷木星。相比短周期行星,确认此区域的行星通常需要望远镜
连续观测几年甚至数十年,因而样本完备性也相对较差\cite{Cumming2008}。根据使用
Keck 望远镜进行的 Lick-Carnegie 行星搜索项目,Rowan 等人通过观测样本估算得到
此类行星的出现概率大概只有约 3\%(文献\citen{Rowan2016}),因而木星在人类想法
中先入为主、见惯非惯的概念也许并不周全。
另外值得提的一点是直接成像法(\S \ref{sec:drctimgmeth})可以观测到更长周期行星
,如 HR 8799 系统(文献\citen{Marois2008HR8799},详见图\ref{fig:hr8799})
等屈指可数的样本。传统的核吸积模型需要花很长时间来形成此类气巨星的胚胎核心。
相比之下,另辟蹊径的引力不稳定(Gravitational/Disk Instability,一般简称 GI )模型
\cite{Kuiper1951,Cameron1978,Boss1997}则可比较合理地解释这些巨行星是如何形成于
距离主星几十个天文单位的轨道上\cite{Durisen2007}。
\subsubsection{超级地球族 Super-Earth Population}
超级地球,有时又称迷你海王星(mini-Neptune),对于 RV 探测到的系外行星一般定义为
介于十个地球质量与海王星质量之间,或者对于凌星法探测到的行星,则其半径大约介于
二到四个地球半径\cite{Haghighipour2011}。在图\ref{fig:exomassper} 中,不难发现超级
地球十分常见与 1 天至 100 天之间的轨道,Batalha 等人利用 $Kepler$ 前 16 个月的数据
也推断超级地球为银河系最常见的行星类别\cite{Batalha2013}。
那么超级地球的形成历史究竟为何般?它们又是否可归至类地行星范畴呢?前面 \S
\ref{sec:clspftheory} 提到在核吸积模型下,气态巨行星首先会成长为约 10
$\tif{M}_\oplus$ 的固态核心后通过吸积气体而成长\cite{Guillot2005}。超级地球质量
也正好估算为此值附近\cite{Lissauer2011MRR},因而很自然的假设便是此类行星在
雪线之外形成并且轨道迁移至如今的位置\cite{Terquem2007,Kennedy2008,IdaLin2008}。
然而也有争论认为在 M 型矮星周围,超级地球完全可以于当地形成($in\:situ$ formation,
文献\citen{Laughlin2004,Kennedy2006,Hansen2013,Chiang2013,Boss2006})。
或许随着越来越多的超级地球内部结构被观测所限制后,此难题才可被解答\cite{Lissauer2014}。
另外,超级地球一般存在于多行星系统中\cite{Borucki2011},而且 Fabrycky 等人发现
这些行星倾向于聚集在平运动共振(Mean Motion Resonance,简称 MMR)的内边缘,
尤其是 3:2 和 2:1 平运动共振\cite{Fabrycky2014}。一时之间包括潮汐作用、共振结构
等解释众说纷纭\cite{LithwickWu2012,Lee2013,Batygin2013,Baruteau2013,Delisle2014,Chatterjee2015},
到如今也尚未弄清其中的动力学机制,以及此现象是否依赖于超级地球的形成过程。
\subsubsection{值得一提的行星系统} \label{sec:pexample}
%https://en.wikipedia.org/wiki/List_of_exoplanet_firsts
\begin{figure}[b]
\centering
\includegraphics[width=1.0\textwidth]{figures/chapter1/fig14_binaryplanet.pdf}
\caption{已知双星内卫星型(S-Type)与行星型(P-Type)行星系统中双星和行星轨道半长径分布图,数据来自 \url{http://openexoplanetcatalogue.com}。}
\label{fig:pibinary}
\end{figure}
除了以上三大类系外行星系统,本文额外汇总了一些有趣的行星系统(Planet of Interest,POI),
它们诡怪的轨道构型在某种程度上令人叹为观止,甚至挑战了现有的行星形成理论。
\textbf{共振系统} --- \textit{GJ 876}。 {} 2001 年,Marcy 等人发现 GJ 876 b 和 c 两颗类木行星周期比
接近 2:1\cite{Marcy2001}。几年后,Rivera 等人更是发现 GJ 876 的额外一颗行星 d 与前两颗
处于 Laplace MMR\cite{Rivera2010}。这和太阳系木星的内侧内伽利略卫星 Io,Europa 以及
Ganymede 的构型堪称如出一辙。理论上,此类轨道构型很有可能是行星在气体盘中迁移所致
\cite{KleyNelson2012,ZhangZhou2010},关于此原理的简单介绍请参考 \S \ref{sec:diskmig}。
\textbf{紧凑系统} --- \textit{Kepler-11}。 {} 它因其 6 颗行星的紧密轨道构型\cite{Lissauer2011}而被称作
太阳系的孪生子系统\cite{Zhou2012}。行星在如此紧凑的系统几乎处于稳定的边缘\cite{Mahajan2014}。
此紧密构型起源也是一直处于争辩中,关于 Kepler-11 行星系统形成理论详见文献\citen{Angelo2016}。
此类系统还有另外一些有代表性的例子,如 HD 10180\cite{Lovis2011} 与 TRAPPIST-1\cite{Gillon2017}。
\textbf{双星系统行星} --- \textit{$\gamma$ Cephei b}。 {} ,双星系统中的行星分为卫星型与行星型
\footnote{卫星型又称 Satellite/Circumprimary Type,构型为行星围绕双星中的一颗运转;行星型
或称 Planet/Circumbinary Type,构型为行星轨道围绕双星的质心。},其中$\gamma$ Cephei
\cite{Hatzes2003} 系前者,而 Kepler-16\cite{Doyle2011} 则属后者。如图 \ref{fig:pibinary} ,近距离双
星中的行星系统与单恒星周围相差迥异,这是因为存在伴星的引力干扰,导致行星系统形成过程大相近
庭。关于此类行星系统,请详见书籍\citen{Haghighipour2010}。
\textbf{高偏心率系统} --- \textit{HD 80606}。 {} 作为偏心率最大的几个系统之一,HD 80606 拥有 0.93
的轨道偏心率\cite{Naef2001},远心点居然是近心点距离的近三十倍。同属此类的系统还有 HD 4113
\cite{Tamuz2008},HD 80606 行星轨道法向和主星自转轴方向测量(spin-orbit measurement )
\cite{Pont2009}表明它的起源很可能与 Lidov-Kozai 机制\cite{Lidov1962,Kozai1962}有着密切关联
\cite{Wu2003},或许是热木星系统高偏心率迁移的证据(详见 \S \ref{sec:highemig})。
\begin{figure}[ht!]
\centering
\includegraphics[height=0.91\textheight]{figures/chapter1/fig15_hjmul.pdf}
\caption{所有拥有热木星的已确认多行星系统汇总图,点的相对大小正比于行星的质量。数据同样来自 \url{http://openexoplanetcatalogue.com}。}
\label{fig:hjwcomp}
\end{figure}
\textbf{极短周期行星} --- \textit{Kepler-78 b}。 {} 极短周期行星又称(Ultra Short Period Planets,
USPs),它们因周期通常在一天以内,和主星表面的距离非常近而得称。代表性行星有 Kepler-78
b(周期为 8.5 小时\cite{SanchisOjeda2013}),Kepler-70 b(周期仅 5.8 小时\cite{Charpinet2011})
等。由于它们受到主星强烈的辐射与引力作用,因而通常拥有非常大的密度,近年来也越发成为行星与
主星物理性质新实验基地\cite{Lopez2016,Moutou2016}。
\textbf{主序后恒星} --- \textit{HD 13189}\cite{Hatzes2005} 。 {} 虽然利用 RV 探测主序后恒星周围行星
的效率并不高\cite{Sato2005},但是此类行星依然对检验行星对恒星的质量、演化历史依赖度极为重要
\cite{Kennedy2008,Johnson2007b,Jones2014}。观测显示巨星周围的行星普遍距离主星较远
($a > 1\,\tif{AU}$),这也许正是主星演化吞噬近距轨道行星的证据\cite{Johnson2007a,Bowler2010}。
\textbf{其他诡怪系统} --- \textit{WASP-47}。 {} 该行星系统拥有一颗典型的热木星\cite{Hellier2012},
然而后续长期检测显示该系统有另外两颗小质量行星\cite{Becker2015,SanchisOjeda2015}以及一颗
长周期冷木行星\cite{NeveuVanMalle2016}。这和传统的认为热木星更倾向于成单的想法截然不同
\cite{Steffen2012},并且即使与其他包含热木星的多行星系统相比,WASP-47 也大有不同(如图
\ref{fig:hjwcomp}),现有的理论并不能完全解释该系统的构型,因而更完善的形成理论与观测限制
也更加迫在眉睫。
另外对疏散星团以及球状星团的巡天显示,Free-floating Planets(FFPs)也许是潜在的数量
最多的行星质量天体\cite{Lucas2000,Bihain2009,Sumi2011},星团环境作为系外(内)行星
的出生环境\cite{Adams2010,Liu2013} ,也必须考虑该环境对行星形成过程的反馈作用(图 \ref{fig:pfenv})。
\section{本文立意}
系外行星从无到有,再到现如今样本越来越丰富,正是因此前人不断去进化仪器、细化观测流程以及
大胆探索数据。因此本文亦从观测作为出发点,讲述怎么利用南极望远镜观测、如何去处理测光
数据处理、如何提高数据的精度(\S \ref{chapter:obs_red})从而服务于后续行星探测以及原恒星盘
搜索等科学目标(\S \ref{chapter:form_evo})。
此外,本文将于 \S \ref{chapter:data_stat} 章详细探讨系外热木星系统的统计性质,回顾解释这类行
星形成的困难之处。最后还包括如何利用轨道-自转不共面性(spin-orbit misalignment )来限制这些
系统的潮汐演化过程和参数,以及这些过程和参数蕴含或还原了哪些系统演化的物理效应。在本文的
最后一部分,将会总结与展望系外行星领域的未来计划。而如何从宏观层面研究行星系统的形成与演
化以及从微观层面去细致描述行星的性质,这些研究将有利于击破现有的系外行星难关。
\begin{figure}[t!]
\centering
\includegraphics[width=1.0\textwidth]{figures/chapter1/fig16_planetenv.pdf}
\caption[不同的恒星环境对系外行星系统形成结果的影响示意图,图片版权 Malmberg 2011。]{不同的恒星环境对系外行星系统形成结果的影响示意图,图片源自文献 \citen{Malmberg2011}。}
\label{fig:pfenv}
\end{figure}