用Python繪制一個仿黑洞圖像

簡介

黑洞圖像大傢都知道,畢竟前幾年剛發佈的時候曾火遍全網,甚至都做成表情包瞭。

問題在於,憑什麼認為這就是黑洞的照片,而不是一個甜甜圈啥的給整模糊瞭得到的呢?有什麼理論依據嗎?

單位制

利用einsteinpy模塊中的Shadow類,可以實現黑洞吸積盤的發射強度,換言之,用理論模擬一下黑洞的照片應該是什麼樣子的。

22年5月份發佈的是人馬座A*的照片,這個"黑洞"的質量為8.26×1036kg,距離地球26673光年。

根據源碼推測,einstenipy中采用幾何單位制,光速和萬有引力常數設為1, 即c=G=1。現設幾何單位制下的時間、距離以及質量單位為S,M,KG,國際制單位為s,m,kg,若以時間單位秒為基本單位,即1s=1S,由於幾何單位制下光速Sc=1M/S,則可得到距離單位的換算關系

同理,根據萬有引力常數可得到質量單位的換算關系

所以質量和距離可以分別表示為

觀測繪圖

盡管我們在地球上觀測這個黑洞,但由於相機鏡頭等因素,實際上肯定不能用這麼大的尺度來觀測,所以稍微選一個差不多的距離就好

import astropy.units as u
from einsteinpy.rays import Shadow
mass = 1.98 * u.kg         
fov  = 100 * u.km
shadow = Shadow(mass=mass, fov=fov, n_rays=1000)

shadow就是其發射模型,其內部用於繪圖的主要成員有

  • fb1 x xx負半軸
  • fb2 x xx正半軸
  • intensity 發射強度

下面繪制一下不同距離處的發射強度

import matplotlib.pyplot as plt
plt.plot(shadow.fb1, shadow.intensity, "r")
plt.plot(shadow.fb2, shadow.intensity, "r")
plt.xlabel("Impact Paramter (b)")
plt.ylabel("Intensity (Emissivity)")
plt.title("Intensity Plot")
plt.grid()
plt.show()

結果如下圖所示,果然中間凹瞭下去

einsteinpy內嵌瞭繪圖函數,通過ShadowPlotter(shadow, is_line_plot=True)可以更加方便地繪制上述圖像,如果將is_line_plot設為False,則可得到強度分佈圖

from einsteinpy.plotting import ShadowPlotter
obj = ShadowPlotter(shadow=shadow, is_line_plot=False)
obj.plot()
obj.show()
plt.show()

這樣一看是不是覺得和前面的那個甜甜圈確有相似之處。

這個模型是根據Cosimo Bambi的文章編寫的,論文地址在這:Cosimo Bambi, 10.1103/PhysRevD.87.107501。本來合計著把這篇文章的公式啥的粘過來解讀一下,但考慮到貼在這也沒人看,想想還是算瞭。

到此這篇關於用Python繪制一個仿黑洞圖像的文章就介紹到這瞭,更多相關Python仿黑洞內容請搜索WalkonNet以前的文章或繼續瀏覽下面的相關文章希望大傢以後多多支持WalkonNet!

推薦閱讀: