數值相對論

數值相對論(英語:numerical relativity)是廣義相對論的一個分支,旨在透過數值方法求解愛因斯坦場方程式,以模擬強重力場中的物理過程。相對論天文學中的物理系統,如重力塌縮中子星黑洞重力波等等,以及其他不能利用弱場低速情形中結論[a]近似的現象都可以利用數值相對論模擬。[1]

恆星塌縮為黑洞並釋放出重力波的過程中的一個時刻。圖中白色球體是黑洞的表觀視界,周圍的彩色波表示多極重力波(l=2,m=2)。右下角的圖是遠處的觀測者採集到的波形。

由於愛因斯坦方程式的複雜性與非線性,[b]這一領域的模擬需要特定的數值方法。[2]而計算量巨大的三維問題則需要藉助超級計算機解決。一些數學與天體物理學問題,比如密接聯星及其重力波的數值模擬,目前可以利用數值相對論求解。[1]

概述

 
首次被人類直接探測到的重力波,GW150914。圖中展示了實驗探測結果與利用數值相對論模擬的結果。

數值相對論的要旨在於,利用數值計算方法研究目前不能解析地描述的重力場。這些重力場可動可靜,既可能是真空場,也可能包含物質場。[c]靜止場問題主要研究場的穩定性。而動態重力場問題則分為兩個主要方面,初值問題與演化問題。兩種問題各自有不同的求解方法。[3]

數值相對論可用於研究宇宙學模型,比如重力塌縮的臨界現象及有黑洞或中子星參與的過程(特別是這兩種天體的併合與受到的擾動)。這些情形都涉及到對時空演化的追蹤。描述這些情形的愛因斯坦場方程式可以用不同形式表述。常用方法包含柯西問題方法、特徵線法[4]以及雷奇理論英語Regge theory法等。[5]這些方法都是從某個超曲面上重力場某一時刻的情形開始推導的,也就是說基於初始值按照時間的流向追蹤它向鄰近超曲面的演化。[6]

與其他數值分析問題一樣,數值相對論也需要注意數值解是否穩定並收斂,以及它們適用的初始與邊界條件。而度規與坐標條件的存在更使得數值相對論顯得尤其複雜。能否求出合適的數值解還會受到愛因斯坦方程式表述方式的影響。

由於所涉及的場有所不同,古典場論的許多方法在數值相對論中並不可用。不過,在應對大尺度問題時,數值相對論則與計算流體力學電動力學以及固體力學等計算科學有很多相似之處。數值相對論還涉及到數值分析並發計算偏微分方程式以及幾何學等數學領域。[7]

發展史

理論基礎

阿爾伯特·愛因斯坦在1915年系統地闡釋了廣義相對論[8]與狹義相對論中的情形類似,時間與空間在廣義相對論中統一為時空,而演化則服從愛因斯坦方程式。愛因斯坦方程式是一組耦合非線性偏微分方程組。在這個方程組誕生後的近一個世紀的時間裡,物理學家僅僅解出了少數幾個精確解。並且,其中許多的精確解是在高度對稱的條件下,也就是說在方程式被很大程度簡化時,推導出來的,比如均一且各向同性宇宙的傅里德曼解[9]

為了在更普遍的條件下研究愛因斯坦方程式,數值相對論應運而生。數值求解愛因斯坦方程式的必要前提是將四維時空還原為各自獨立的三維空間與一維時間,也就是時空的「3+1解構」。這個過程可以透過幾種或簡或繁的途徑實現。這個問題的首個突破是由理察·阿諾維特斯坦利·德塞爾以及查爾斯·W·米斯納於1950年代後期利用哈密頓力學沿著保羅·狄拉克給出的途徑實現的。他們對於愛因斯坦方程式的推演結果稱作ADM形式[10]出於技術原因,他們的結果並不適用於數值計算。它的雙曲性很弱,因而在實際計算中很難使用。不過數值相對論中許多用到3+1解構的求解方法還是與ADM形式緊密相關。這種解構會將愛因斯坦方程式整理為具有一定初始條件的柯西問題。這種形式已經是計算機數值求解慣常解決的問題。[11]

時空的坐標也不能唯一確定。即使固定了起始超曲面上的坐標,但在穿越到鄰近的超曲面時,時空坐標也會因不同的取法而不同。[d]這是數值相對論的一個特徵。這種度規上的自由度雖然並不會影響實際上的物理過程,但對過程的描述可能會因其改變。假設從起始的超曲面到鄰近的超曲面,時間與空間坐標上分別移動了  。位移描述方式上的不同,也就是說  的取法不同,會導致推導出的方程組形式發生變化。適當的取法會為數值方程式的數值求解帶來潛在的好處,但坐標和度規許多看來「自然」的取法會導致數值解不穩定,並導致模擬失敗。[12]

 
茨維·皮蘭

在ADM形式發表時,計算機的算力並不足以應對這一計算。1964年,哈恩和林德奎斯特首次嘗試數值求解愛因斯坦方程式。[13]隨後在1970年代,斯馬[14][15]與埃普利[16]繼續這方面的嘗試。這些早期的嘗試主要針對軸對稱(「2+1維」)時空的情況。茨維·皮蘭寫出了首個可以追蹤輻射重力波的柱對稱系統演化的程序。[17]皮蘭在此項工作中為數值相對論中許多概念奠基,比如自由演化與約束演化。它們是推導起始數據隨時間演化的兩種途徑。[18][19]對稱性的引入降低了內存與運算速度的門檻,使得科學家可以透過當時的超級計算機解決一些問題。[17]

早期結果

理察·史塔克與茨維·皮蘭在1980年代早期首次對於真實的天體物理學問題做了數值計算。[20]他們首次計算了旋轉黑洞在產生階段輻射出的重力波。在這個結果發表後的近二十年的時間裡,數值相對論少有結果發表。這可能還是因為計算機算力不夠。1990年代,美國「大挑戰英語Grand Challenges」中的雙黑洞項目組成功地模擬了兩個黑洞迎頭相撞的過程。他們利用了問題軸對稱的特點簡化了求解過程。項目組還在後期處理階段透過這個數值解計算了事件視界[21]

而在三維時空中數值求解愛因斯坦方程式方面最早的嘗試專注於不自轉的史瓦西黑洞。這種黑洞靜止且具有球對稱性。對於驗證數值相對論,這一問題十分理想:一是由於科學家非常了解這個問題的解析解,二是由於不自轉的黑洞的數值解必然對時間收斂,三是由於這個問題還包含了數值模擬中一種非常複雜的模擬對象——黑洞中心的重力奇異點。1995年,安尼諾斯等人發表了這一問題最早的結果之一。[22]在文章中,他們寫道:

目前的計算機缺少足夠的內存與算力以模擬解析的三維時空。這一點阻礙了三維數值相對論的發展。

理論的發展

在之後的幾年裡,隨著計算機功能更為強大,幾組研究人員發展了幾種可以提升運算效率的技術。拉撒路組使用了早前透過非線性ADM方程組對黑洞融合模擬的結果,基於單黑洞線性微擾理論,得到了更為穩定的程序。[23]而從黑洞模型入手,又有兩種可以繞過奇異點問題的技術:切除法與穿刺法。[24]2005年,普利托利烏斯綜合上述方法,採用了適當的坐標條件在模擬雙黑洞方面取得了重大突破。[25]這幾種新方法的穩定性在之後的幾年裡進一步提升,使得人們可以模擬兩個黑洞在併合前彼此旋轉幾十以至幾百轉的情形。除此之外,計算流體力學中的自適應計算網格細化方法也被引入到數值相對論中來。[26]

拉撒路計劃

繼「大挑戰」後,科學家又在拉撒路計劃(1998-2005)中利用天體物理學結果對雙黑洞融合做短期數值模擬。不過當時,在超級計算機上模擬這個系統演變的過程中,對於愛因斯坦方程式的數值求解無論如何都得不到穩定的結果。研究者前後結合了多種近似方法(比如後牛頓方法和單黑洞擾動)來對全過程做數值模擬。[23]

拉撒路計劃標誌著雙黑洞問題當時最為前沿的進展,給出了許多可以用於天體物理學研究的相當精準結果,比如重力波帶走的能量與角動量,[27][28]不同質量的黑洞融合時發出的脈衝[29]以及融合後的黑洞的最終質量、動量和角動量。[30]科學家還在這項計劃中計算了融合過程中釋出的重力波的具體形式。這項結果對於重力波探測器的研製非常重要。同時,他們還預言黑洞撞擊會瞬間產生宇宙中最為強大的能暴。在僅僅一秒內以重力波形式放出的能量足以大過銀河系自產生以來其中恆星放出的能量總和。而重力波輻射只占系統總質量減量的很小一部分。[31]

切除法

切除法在1990年代後期提出,[32]是指在模擬演化過程中將奇異點周圍的事件視界「切除」掉。理論上,由於因果律以及視界本身的性質,[e]這種方法不會影響到對於切除部分之外部分的演化的模擬。也就是說,即使沒有解出黑洞內的情況也不會妨礙到外部的精確求解。人們可以透過設置適當的邊界條件實現「切除」的過程。[33]

儘管切除法在實際應用中非常成功,但它也存在兩個小缺陷。首先,在使用切除法時,人們必須審慎地選用坐標條件。儘管物理效應可能不會傳播到視界外,但坐標效應可以。比如,在使用橢圓坐標條件時,黑洞內的網格變化會立即傳播到視界外。[34]這意味著,人們需要採用特徵速度小於光速的雙曲型坐標條件(比如諧和坐標條件)。[35]其次,在黑洞運動時,所切除的部分也必須隨黑洞一起運動。[33]

切除法在隨後幾年裡隨著校準條件的改進穩定性逐步提高。並且科學家也解決了切除部分隨運算區域運動的問題。[36][37][38][39][40][35]藉由切除法,雙黑洞運動軌跡及融合過程的首個長期穩定的計算結果於2005年發表。[25]

穿刺法

在「穿刺法」中,解分為包含黑洞奇異點的解析部分(也就是「刺」)[41]和奇異點之外數值構造的部分。這種方法是由布里爾-林德奎斯特算法(針對起始時靜止的黑洞)推廣而得。[42]而如果再做一般化處理則可得到鮑恩-約克算法(針對起始時旋轉且運動的黑洞)。[43]直到2005年,所有採用穿刺法的結果都需要刺在模擬過程中坐標位置不變。不過,彼此密聯的黑洞自然會受到重力的作用運動。這意味著坐標系本身或會「伸展」或會「扭曲」。這會導致模擬的部分階段出現數值不穩定。另外那種避開奇異點的方法在模擬物質塌縮形成黑洞的過程中,由於也會採用這種坐標條件,同樣也會造成類似效應。[44]

2005年,研究人員首次展示了刺在坐標系中運動的可行性,從而解決了這種方法早前存在的部分問題,使其可以用於精確追蹤黑洞的長期演化。[25][45][46]透過選擇合適的坐標條件並對奇異點周圍的物理場做粗略的估計,[f]人們可以得到繞著彼此旋轉的雙黑洞的數值解,並能精確算出它們釋放的重力輻射。[47]

自適應網格細化

自適應網格細化英語Adaptive mesh refinement是一種在數值相對論萌芽前即已發展的數值算法。它是在1980年代由喬普提克在研究純量場塌縮的臨界現象時引入的。在那種情形中,場參數處於黑洞形成與空間擴張的臨界點。[48][49]由於使用球對稱,這項工作針對的是一維情形。之後這種算法被推廣到二維情形中。[50]二維網格法也用於研究非均勻宇宙[51][52]和史瓦西黑洞。[53]目前,自適應網格已經成為數值相對論的標準工具,用於研究黑洞與其他稠密系統的融合,以及這種事件產生的重力波。[54][55]

近期進展

 
對於兩個中子星併合生成黑洞過程的可視化模擬(NASA,2014)

目前,數值相對論領域每年都會有數十以至數百篇的相關論文發表,內容涵蓋廣義相對論數學、重力波以及涉及旋轉黑洞的天體物理學問題。所使用的方法在對中子星、黑洞等星體組成的雙星系統和多黑洞系統的研究中也有運用。[56][57]從這些工作中,科學家預測兩個旋轉黑洞併合產生的黑洞的速度能達到4000以至10000 km/s,可以讓它穿越所有已知的星系。[58][59]而釋出的能量占到總體靜質量的8%。黑洞旋轉軸所可能發出的劇變則可以解釋從無線電星系中觀測到的相對論性噴流[60]這一領域另一個重要的研究方向是分類LIGOVIRGO不能探測到的黑洞併合產生的重力輻射。[61]

新近算法的精確程度足以用來驗證重力波的探測結果,比如人類首次直接探測到的重力波GW150914,實驗數據與數值模擬數據的誤差在4%以內。[62]

腳註

注釋

  1. ^ 比如基於愛因斯坦方程式精確解的後牛頓解構微擾理論
  2. ^ 比如在不同的方程式表述以及初始和邊界條件下,時間演化的雙曲性柯西問題適定性
  3. ^ 在廣義相對論中,重力場外的其他場都可稱為物質場。
  4. ^ 類似地,在狹義相對論中,時間流向與速度在不同的慣性系中也不同。
  5. ^ 視界內部的交互作用不會影響到外部的物理性質
  6. ^ 由於沒有物理效應可從黑洞逃逸,這種估計並不會影響到整體結果。

引注

  1. ^ 1.0 1.1 Baumgarte T. W., Shapiro S. L. Numerical Relativity 2010,Preface. Sec. What is numerical relativity?.
  2. ^ Choptuik M. W. Computational Methods in General Relativity: The Theory 2006.
  3. ^ Baumgarte T. W., Shapiro S. L. Numerical Relativity 2010,Sec. 2.6 The constraint and evolution equations and Introduction of Ch. 3 Constructing initial data.
  4. ^ Winicour, J., Characteristic Evolution and Matching, Living Reviews in Relativity, 2012, 15 (2), Bibcode:2012LRR....15....2W, doi:10.12942/lrr-2012-2, (原始內容存檔於2015年8月15日) 
  5. ^ Gentle, Adrian P., Regge Calculus: A Unique Tool for Numerical Relativity, General Relativity and Gravitation, 2002, 34 (10): 1701—1718 [2018-03-14], ISSN 0001-7701, doi:10.1023/A:1020128425143, (原始內容存檔於2018-03-14) 
  6. ^ Cook, G., Initial Data for Numerical Relativity, Living Reviews in Relativity, 2000, 3 (5), Bibcode:2000LRR.....3....5C, arXiv:gr-qc/0007085 , doi:10.12942/lrr-2000-5, (原始內容存檔於2006年7月12日) 
  7. ^ Baumgarte T. W., Shapiro S. L. Numerical Relativity 2010,Ch. 6 Numerical methods.
  8. ^ Einstein, Albert. Der Feldgleichungen der Gravitation. Sitzungsberiche der Deutschen Akademie der Wissenschaften zu Berlin, Klasse fur Mathematik, Physik, und Technik. 
  9. ^ Stephani, H.; Kramer, D.; MacCallum, M.; Hoenselaers, C.; Herlt, E., Sec. 1.1 What are exact solutions, and why study them?, Exact Solutions of Einstein's Field Equations, Cambridge Monographs on Mathematical Physics, Cambridge University Press, 2003 [2018-03-14], ISBN 9781139435024, (原始內容存檔於2016-05-29) 
  10. ^ Arnowitt, R.; Deser, S.; Misner, C. W., The dynamics of general relativity, L. Witten (編), Gravitation: An Introduction to Current Research, New York: Wiley: 227–265, 1962 
  11. ^ Baumgarte T. W., Shapiro S. L. Numerical Relativity 2010,Ch. 2 The 3+1 decompostion of Einstein’s equations.
  12. ^ Baumgarte T. W., Shapiro S. L. Numerical Relativity 2010,Ch. 4 Choosing coordinates: the lapse and shift.
  13. ^ Hahn, S. G.; Lindquist, R. W., The two-body problem in geometrodynamics, Annals of Physics, 1964, 29: 304–331, Bibcode:1964AnPhy..29..304H, doi:10.1016/0003-4916(64)90223-4 
  14. ^ Smarr, L. The Structure of General Relativity with a Numerical Example (Ph.D.論文). University of Texas, Austin. 1975. 
  15. ^ Smarr, L., Spacetimes generated by computers: Black holes with gravitational radiation, Annals of the New York Academy of Sciences, 1977, 302: 569, doi:10.1111/j.1749-6632.1977.tb37076.x 
  16. ^ Eppley, K. The numerical evolution of the collision of two black holes (Ph.D. Dissertation論文). Princeton University. 1975. 
  17. ^ 17.0 17.1 Piran, T., Cylindrical general relativistic collapse, Physical Review Letters, 1978, 41 (16): 1085, Bibcode:1978PhRvL..41.1085P, doi:10.1103/PhysRevLett.41.1085 
  18. ^ Alcubierre M. Introduction to 3+1 Numerical Relativity 2008,Sec. 2.6 Free versus constrained evolution.
  19. ^ Bona C., Palenzuela-Luque C., Bona-Casas C. Elements of Numerical Relativity and Relativistic Hydrodynamics 2009,Subsec. 2.3.3 Evolution strategies.
  20. ^ Stark, R. F.; Piran, T., Gravitational-wave emission from rotating gravitational collapse, Physical Review Letters, 1985, 55: 891–894, Bibcode:1985PhRvL..55..891S, doi:10.1103/PhysRevLett.55.891 
  21. ^ Matzner, R. A.; Seidel, H. E.; Shapiro, S. L.; Smarr, L.; Suen, W.-M.; Teukolsky, S. A.; Winicour, J., Geometry of a Black Hole Collision, Science, 1995, 270 (5238): 941–947, Bibcode:1995Sci...270..941M, doi:10.1126/science.270.5238.941 
  22. ^ Anninos, P.; Massó, J.; Seidel, E.; Suen, W.-M.; Towns, J., Three-dimensional numerical relativity: The evolution of black holes, Physical Review D, 1995, 52 (4): 2059–2082, Bibcode:1995PhRvD..52.2059A, arXiv:gr-qc/9503025 , doi:10.1103/PhysRevD.52.2059 
  23. ^ 23.0 23.1 Baker, J.; Campanelli, M.; Lousto, C. O., The Lazarus project: A pragmatic approach to binary black hole evolutions, Physical Review D, 2002, 65 (4): 044001, Bibcode:2002PhRvD..65d4001B, arXiv:gr-qc/0104063 , doi:10.1103/PhysRevD.65.044001 
  24. ^ Baumgarte T. W., Shapiro S. L. Numerical Relativity 2010,Ch. 13 Binary black hole evolution.
  25. ^ 25.0 25.1 25.2 Pretorius, F., Evolution of Binary Black-Hole Spacetimes, Physical Review Letters, 2005, 95 (12): 121101, Bibcode:2005PhRvL..95l1101P, arXiv:gr-qc/0507014 , doi:10.1103/PhysRevLett.95.121101 
  26. ^ Baumgarte T. W., Shapiro S. L. Numerical Relativity 2010,Subsec. 6.2.5 Mesh refinement.
  27. ^ Baker, J.; Brügmann, B.; Campanelli, M.; Lousto, C. O.; Takahashi, R., Plunge Waveforms from Inspiralling Binary Black Holes, Physical Review Letters, 2001, 87 (12): 121103, Bibcode:2001PhRvL..87l1103B, arXiv:gr-qc/0102037 , doi:10.1103/PhysRevLett.87.121103 
  28. ^ Baker, J.; Campanelli, M.; Lousto, C. O.; Takahashi, R., Modeling gravitational radiation from coalescing binary black holes, Physical Review D, 2002, 65 (12): 124012, Bibcode:2002PhRvD..65l4012B, arXiv:astro-ph/0202469 , doi:10.1103/PhysRevD.65.124012 
  29. ^ Campanelli, M., Understanding the fate of merging supermassive black holes, Classical and Quantum Gravity, 2005, 22: 387, Bibcode:2005CQGra..22S.387C, arXiv:astro-ph/0411744 , doi:10.1088/0264-9381/22/10/034 
  30. ^ Baker, J.; Campanelli, M.; Lousto, C. O.; Takahashi, R., Coalescence remnant of spinning binary black holes, Physical Review D, 2004, 69 (2): 027505, Bibcode:2004PhRvD..69b7505B, arXiv:astro-ph/0305287 , doi:10.1103/PhysRevD.69.027505 
  31. ^ Berti, E.; Cardoso, V.; Gonzalez, J. A.; Sperhake, U.; Hannam, M.; Husa, S.; Brügmann, B., Inspiral, merger, and ringdown of unequal mass black hole binaries: A multipolar analysis, Physical Review D, 2007, 76 (6): 064034 (1—40), Bibcode:2007PhRvD..76f4034B, arXiv:gr-qc/0703053 , doi:10.1103/PhysRevD.76.064034 
  32. ^ Alcubierre, M.; Brügmann, B., Simple excision of a black hole in 3+1 numerical relativity, Physical Review D, 2001, 63 (10): 104006, Bibcode:2001PhRvD..63j4006A, arXiv:gr-qc/0008067 , doi:10.1103/PhysRevD.63.104006 
  33. ^ 33.0 33.1 Baumgarte T. W., Shapiro S. L. Numerical Relativity 2010,Subsec. 13.1.2 Black hole excision.
  34. ^ Bona C., Palenzuela-Luque C., Bona-Casas C. Elements of Numerical Relativity and Relativistic Hydrodynamics 2009,Subsec. 6.1.3 Regular initial data.
  35. ^ 35.0 35.1 Shoemaker, D.; et al., Moving black holes via singularity excision, Classical and Quantum Gravity, 2003, 20 (16): 3729-3743, Bibcode:2003CQGra..20.3729S, arXiv:gr-qc/0301111 , doi:10.1088/0264-9381/20/16/313 
  36. ^ Bona, C.; Masso, J.; Seidel, E.; Stela, J., New formalism for numerical relativity, Physical Review Letters, 1995, 75 (4): 600–603, Bibcode:1995PhRvL..75..600B, arXiv:gr-qc/9412071 , doi:10.1103/PhysRevLett.75.600 
  37. ^ Cook, G. B.; et al., Boosted Three-Dimensional Black-Hole Evolutions with Singularity Excision, Physical Review Letters, 1998, 80: 2512–2516, Bibcode:1998PhRvL..80.2512C, arXiv:gr-qc/9711078 , doi:10.1103/PhysRevLett.80.2512 
  38. ^ Alcubierre, M., Hyperbolic slicings of spacetime: singularity avoidance and gauge shocks, Classical and Quantum Gravity, 2003, (20): 607–623, Bibcode:2003CQGra..20..607A, arXiv:gr-qc/0210050  
  39. ^ Alcubierre, M.; et al., Gauge conditions for long-term numerical black hole evolutions without excision, Physical Review D, 2003, 67 (8): 084023, Bibcode:2003PhRvD..67h4023A, arXiv:gr-qc/0206072 , doi:10.1103/PhysRevD.67.084023 
  40. ^ Brügmann, B.; Tichy, W.; Jansen, N., Numerical Simulation of Orbiting Black Holes, Physical Review Letters, 2004, 92 (21): 211101, Bibcode:2004PhRvL..92u1101B, arXiv:gr-qc/0312112 , doi:10.1103/PhysRevLett.92.211101 
  41. ^ Brandt, S.; Brügmann, B., A Simple Construction of Initial Data for Multiple Black Holes, Physical Review Letters, 1997, 78 (19): 3606-3609, Bibcode:1997PhRvL..78.3606B, arXiv:gr-qc/9703066 , doi:10.1103/PhysRevLett.78.3606 
  42. ^ Brill, D. R.; Lindquist, R. W., Interaction Energy in Geometrostatics, Physical Review, 1963, 131 (1): 471–476, Bibcode:1963PhRv..131..471B, doi:10.1103/PhysRev.131.471 
  43. ^ Bowen, J. M.; York Jr., J. W., Time-asymmetric initial data for black holes and black-hole collisions, Physical Review D, 1980, 21 (8): 2047–2056, Bibcode:1980PhRvD..21.2047B, doi:10.1103/PhysRevD.21.2047 
  44. ^ Baumgarte T. W., Shapiro S. L. Numerical Relativity 2010,Subsec. 13.1.1 Singularity avoiding coordinates.
  45. ^ Baker, J. G.; et al., Gravitational-Wave Extraction from an Inspiraling Configuration of Merging Black Holes, Physical Review Letters, 2006, 96 (11): 111102, Bibcode:2006PhRvL..96k1102B, arXiv:gr-qc/0511103 , doi:10.1103/PhysRevLett.96.111102 
  46. ^ Campanelli, M.; et al., Accurate Evolutions of Orbiting Black-Hole Binaries without Excision, Physical Review Letters, 2006, 96 (11): 111101, Bibcode:2006PhRvL..96k1101C, arXiv:gr-qc/0511048 , doi:10.1103/PhysRevLett.96.111101 
  47. ^ Baumgarte T. W., Shapiro S. L. Numerical Relativity 2010,Subsec. 13.1.3 The moving puncture method.
  48. ^ Choptuik, M. W. Experiences with an adaptive mesh refinement algorithm in numerical relativity. Evans, C.; Finn, L.; Hobill, D. (編). Frontiers in numerical relativity. Cambridge: Cambridge University Press. 1989. ISBN 0521366666. 
  49. ^ Choptuik, M. W., Universality and scaling in gravitational collapse of a massless scalar field, Physical Review Letters, 1993, 70 (1): 9–12, Bibcode:1993PhRvL..70....9C, doi:10.1103/PhysRevLett.70.9 
  50. ^ Choptuik, M. W.; et al., Critical collapse of the massless scalar field in axisymmetry, Physical Review D, 2003, 68 (4): 044007, Bibcode:2003PhRvD..68d4007C, arXiv:gr-qc/0305003 , doi:10.1103/PhysRevD.68.044007 
  51. ^ Hern, S. D. Numerical relativity and inhomogeneous cosmologies (Ph.D.論文). Cambridge University. 1999. 
  52. ^ Belanger, Z. B. Adaptive mesh refinement in the T2 symmetric spacetime (Master論文). Oakland University. 2001. 
  53. ^ Schnetter, E.; et al., Evolutions in 3D numerical relativity using fixed mesh refinement, Classical and Quantum Gravity, 2004, 21: 1465–1488, Bibcode:2004CQGra..21.1465S, arXiv:gr-qc/0310042 , doi:10.1088/0264-9381/21/6/014 
  54. ^ Imbiriba, B.; et al., Evolving a puncture black hole with fixed mesh refinement, Physical Review D, 2004, 70 (12): 124025, Bibcode:2004PhRvD..70l4025I, arXiv:gr-qc/0403048 , doi:10.1103/PhysRevD.70.124025 
  55. ^ Fiske, D. R.; et al., Wave zone extraction of gravitational radiation in three-dimensional numerical relativity, Physical Review D, 2005, 71 (10): 104036, Bibcode:2005PhRvD..71j4036F, arXiv:gr-qc/0503100 , doi:10.1103/PhysRevD.71.104036 
  56. ^ Etienne, Zachariah B.; Liu, Yuk Tung; Shapiro, Stuart L.; Baumgarte, Thomas W. Relativistic Simulations of Black Hole-Neutron Star Mergers: Effects of black-hole spin. Phys. Rev. D. 2009, 76 (4). Bibcode:2009PhRvD..79d4024E. arXiv:0812.2245 . doi:10.1103/PhysRevD.79.044024. 
  57. ^ Lousto, C. O.; Zlochower, Y., Foundations of multiple-black-hole evolutions, Physical Review D, 2008, 77 (2): 024034, Bibcode:2008PhRvD..77b4034L, arXiv:0711.1165 , doi:10.1103/PhysRevD.77.024034 
  58. ^ Campanelli M.; et al., Maximum Gravitational Recoil, Physical Review Letters, 2007, 98 (23): 231102, Bibcode:2007PhRvL..98w1102C, arXiv:gr-qc/0702133 , doi:10.1103/PhysRevLett.98.231102 
  59. ^ Healy, J.; et al., Superkicks in Hyperbolic Encounters of Binary Black Holes, Physical Review Letters, 2009, 102 (4): 041101, Bibcode:2009PhRvL.102d1101H, arXiv:0807.3292 , doi:10.1103/PhysRevLett.102.041101 
  60. ^ Campanelli, M.; et al., Spin flips and precession in black-hole-binary mergers, Physical Review D, 2007, 75 (6): 064030, Bibcode:2007PhRvD..75f4030C, arXiv:gr-qc/0612076 , doi:10.1103/PhysRevD.75.064030 
  61. ^ Hinder, I.; et al., Error-analysis and comparison to analytical models of numerical waveforms produced by the NRAR Collaboration, Classical and Quantum Gravity, 2013, 31 (2): 025012, Bibcode:2013CQGra..31b5012H, arXiv:1307.5307 , doi:10.1088/0264-9381/31/2/025012 
  62. ^ The LIGO Scientific Collaboration; The Virgo Collaboration, Tests of general relativity with GW150914 (PDF) ArXiv e-prints, 2016 [2018-03-15], Bibcode:2016arXiv160203841T, arXiv:1602.03841 , (原始內容 (PDF)存檔於2016-02-15) 

引用書目

延伸閱讀

外部連結