先日ブラックホールの撮影に成功したというニュースがありましたね。そのブラックホールの大きさをプログラムで求めてみました。解説動画はこちらまずはブラックホールって何なのって話ですがウィキペディアによると極めて高密度かつ大質量で強い重力のために物質だけでなく光さえ脱出することができない天体である。元は星だったりとかそんなやーつですね。そのあとの説明では周囲は非常に強い重力によって時空が著しくゆがめられある半径より内側では脱出速度が光速を超えてしまう。この半径をシュヴァルツシルト半径この半径を持つ球面を事象の地平面(シュヴァルツシルト面)と呼ぶ。この事象の地平面がブラックホールの大きさになるわけですがその大きさを調べるにはシュヴァルツシルト半径が分かればその直径も分かりますね。シュヴァルツシルト半径のwikiの記事も見ると1916年カール・シュヴァルツシルトはアインシュタインの重力場方程式の解を求め非常に小さく重い星があったとするとその星の中心からのある半径の球面内では光も脱出できなくなるほど曲がった時空領域が出現することに気づいた。その半径をシュワルツシルト半径と呼びシュワルツシルト半径よりも小さいサイズに収縮した天体はブラックホールと呼ばれる。それを求める公式があるそうです。これが関係する公式です。計算していきましょう。地球を縮めたらどれくらいの大きさでブラックホールになるのかです。重力定数と光速度はscipyのconstantsというライブラリの中に定数が有るのでそれを用いると計算がスムーズです。constantsお次は半径を求める式をプログラムで当てはめていきます。100)0.8870107529844293さて結果は・・・0.88センチ地球を1センチ程度に圧縮したらブラックホールになり得るのだそう。(C**2)print(R)2954.0291803731太陽は約3キロメートルほどでブラックホールになるんだそう。逆に半径から質量を求めてみましょう。G)print(M)5.9724e+24はいこれで半径から質量も出すことができましたね。最後にこの前撮影されたブラックホールの半径を求めます。(C**2)print(R/1000)19201189672.42515約192億キロです。デカすぎワロタwwwwwよくそんなの撮れましたね。最後の賭けでブラックホールに入ってこの世界の謎を解き明かしてラストは・・・うーーん奥が深かった。理論物理学者のキップ・ソーンが科学コンサルタント兼製作総指揮を務め理論武装も完璧であろうと思われる。物理学の話でもあるけど親子の絆を描く作品でもある不朽の名作ぜひお試しあれ。CodeWhispererで爆速コーディングコーディングに迷った際に役に立ちWEB検索する時間を無くして作業時間を効率化できます。