MathLabs

Vật lý toán

Hình học của lỗ đen

Nghiệm Schwarzschild và Kerr của phương trình trường Einstein: chân trời sự kiện, mặt cầu photon và quỹ đạo tròn ổn định trong cùng, độ lệch ánh sáng, định lý kỳ dị, và sóng hấp dẫn, gắn với kính viễn vọng chân trời sự kiện và LIGO.

Trực giácRơi vào một cái giếng trong không-thời gian

Hãy hình dung không-thời gian không phải là một sân khấu trống rỗng, mà như một tấm cao su co giãn: một khối lượng đặt lên đó làm tấm cao su lõm xuống, và mọi thứ ở gần — hành tinh, ánh sáng, thậm chí cả thời gian — đều phải di chuyển qua chỗ lõm đó. Dồn đủ khối lượng vào một vùng đủ nhỏ, chỗ lõm trở thành một cái giếng không đáy: ánh sáng rơi vào đó không thể leo trở ra được. Vùng không có đường quay lại đó là một lỗ đen, và biên của nó là chân trời sự kiện. Trang này tính toán chính xác cái giếng đó sâu bao nhiêu, cái gì rơi vào và cái gì chỉ bị bẻ cong quanh nó, và làm sao giờ đây ta có thể nhìn thấy bóng mà cái giếng đó đổ ra và nghe được gợn sóng khi hai cái giếng như vậy va vào nhau.

Một mặt hình phễu loe rộng ra và phẳng dần ở xa tâm, rồi cong dốc xuống và thu hẹp gần một cổ hẹp ở tâm, biểu diễn cách không gian giãn ra gần chân trời sự kiện của lỗ đen.
Mặt paraboloid Flamm z(r)=2rs(r−rs)z(r) = 2\sqrt{r_s(r-r_s)}: một cách nhúng hình học không gian quanh một khối lượng vào không gian 3 chiều thông thường, cho thấy khoảng cách bán kính giãn ra gần r=rsr=r_s thế nào. Cái phễu dốc dần chỉ là một hình ảnh minh họa, không phải hình học 4 chiều thật — nhưng cổ phễu tại r=rsr=r_s là thật: đó chính là chân trời sự kiện.

Phổ thôngCái giếng lớn cỡ nào? Bán kính Schwarzschild

Định nghĩa: Vận tốc thoát và bán kính Schwarzschild

Vận tốc thoát kiểu Newton khỏi một quả cầu khối lượng MM bán kính rr là vesc=2GM/rv_{\text{esc}}=\sqrt{2GM/r}: tốc độ tối thiểu để bay xa mãi mãi. Đặt vesc=cv_{\text{esc}}=c, tốc độ ánh sáng, và giải theo rr cho ta một bán kính đặc biệt, bán kính Schwarzschild rs=2GMc2r_s = \dfrac{2GM}{c^2} — bán kính mà tại đó ngay cả ánh sáng cũng không thoát được. Cách tính tắt kiểu Newton này cho ra đúng con số, dù lý do thực sự khiến ánh sáng không thoát được (độ cong không-thời gian, không phải một "lực hấp dẫn" mà ánh sáng phải chống lại) hoàn toàn mang tính tương đối tính; phép suy ra nghiêm túc ở dưới.

rs=2GMc2r_s = \dfrac{2GM}{c^2}

Ở đây GG là hằng số hấp dẫn Newton, MM là khối lượng, và cc là tốc độ ánh sáng. Với Mặt Trời, rs≈2.95 kmr_s\approx 2.95\text{ km} — Mặt Trời cần bị nén xuống cỡ một thị trấn nhỏ để trở thành lỗ đen; với Trái Đất, rs≈8.9 mmr_s\approx 8.9\text{ mm}, nhỏ hơn một quả nho.

Đại họcMêtric Schwarzschild

Định nghĩa: Mêtric Schwarzschild

Năm 1916, Karl Schwarzschild tìm ra nghiệm chân không đối xứng cầu duy nhất của phương trình trường Einstein. Trong hệ tọa độ (t,r,θ,ϕ)(t,r,\theta,\phi), khoảng cách không-thời gian là ds2=−(1−rsr)c2 dt2+(1−rsr)−1dr2+r2 dΩ2ds^2 = -\left(1-\dfrac{r_s}{r}\right)c^2\,dt^2 + \left(1-\dfrac{r_s}{r}\right)^{-1}dr^2 + r^2\,d\Omega^2, trong đó dΩ2=dθ2+sin⁡2θ dϕ2d\Omega^2 = d\theta^2+\sin^2\theta\,d\phi^2 là mêtric trên mặt cầu đơn vị và rs=2GM/c2r_s=2GM/c^2 như trên. Nó mô tả không-thời gian bên ngoài mọi khối lượng đối xứng cầu, không quay — một ngôi sao, một hành tinh, hay một lỗ đen — và trở về không-thời gian Minkowski phẳng khi r→∞r\to\infty.

ds2=−(1−rsr)c2 dt2+(1−rsr)−1dr2+r2 dΩ2ds^2 = -\left(1-\dfrac{r_s}{r}\right)c^2\,dt^2 + \left(1-\dfrac{r_s}{r}\right)^{-1}dr^2 + r^2\,d\Omega^2

Hai bán kính trông có vẻ kỳ dị: r=0r=0, nơi cả gttg_{tt} và grrg_{rr} đều bùng nổ, và r=rsr=r_s, nơi gtt→0g_{tt}\to0 và grr→∞g_{rr}\to\infty. Chỉ một trong hai điểm này là kỳ dị vật lý thật sự. Vô hướng Kretschmann K=RαβγδRαβγδ=48 G2M2c4r6K = R_{\alpha\beta\gamma\delta}R^{\alpha\beta\gamma\delta} = \dfrac{48\,G^2M^2}{c^4 r^6} là một độ đo độ cong không phụ thuộc tọa độ, được xây từ toàn bộ tensor Riemann; nó phân kỳ tại r=0r=0 (độ cong vô hạn thật sự, không thể tránh) nhưng vẫn hoàn toàn hữu hạn tại r=rsr=r_s. Vậy sự bùng nổ tại r=rsr=r_s là một kỳ dị tọa độ — sản phẩm phụ của hệ tọa độ Schwarzschild, không phải của bản thân không-thời gian — hệt như "kỳ dị" ở cực bắc của các đường kinh tuyến trên quả địa cầu chỉ là sản phẩm phụ của tọa độ vĩ độ–kinh độ, không phải khiếm khuyết của mặt cầu. Chuyển sang tọa độ Eddington–Finkelstein (thay tt bằng tọa độ null đi vào v=t+r∗v=t+r_*, r∗=r+rsln⁡∣r/rs−1∣r_*=r+r_s\ln|r/r_s-1|) hoặc tọa độ Kruskal–Szekeres mở rộng đầy đủ làm mọi thành phần mêtric trở nên chính quy xuyên suốt qua r=rsr=r_s, cho thấy đó là một mặt null hoàn toàn trơn: chân trời sự kiện.

Các bán kính chủ chốt quanh một lỗ đen Schwarzschild
Đại lượngCông thứcTheo đơn vị rsr_s
Chân trời sự kiệnrs=2GM/c2r_s=2GM/c^21 rs1\,r_s
Mặt cầu photon (quỹ đạo ánh sáng tròn không ổn định)rph=32rs=3GMc2r_{\text{ph}} = \dfrac{3}{2}r_s = \dfrac{3GM}{c^2}1.5 rs1.5\,r_s
Quỹ đạo tròn ổn định trong cùng (ISCO)rISCO=3rs=6GMc2r_{\text{ISCO}} = 3r_s = \dfrac{6GM}{c^2}3 rs3\,r_s
Bán kính bóng (như EHT chụp được)rshadow=332rs≈2.6 rsr_{\text{shadow}} = \dfrac{3\sqrt3}{2}r_s \approx 2.6\,r_s≈2.6 rs\approx2.6\,r_s

Đại họcQuỹ đạo quanh lỗ đen: mặt cầu photon và ISCO

Cả ánh sáng lẫn hạt có khối lượng chuyển động trong hình học Schwarzschild đều bảo toàn hai đại lượng dọc theo quỹ đạo trắc địa của chúng — một hằng số giống năng lượng và một hằng số giống mô men động lượng — hệt như trong bài toán Kepler, nhưng với một biến thể tương đối tính trong thế năng hiệu dụng không có đối tác Newton nào và tạo ra hai bán kính mang tính cột mốc.

Với một lỗ đen Schwarzschild, quỹ đạo tròn của photon chỉ tồn tại tại rph=32rs=3GMc2r_{\text{ph}} = \dfrac{3}{2}r_s = \dfrac{3GM}{c^2} luôn không ổn định; quỹ đạo tròn của hạt có khối lượng chỉ ổn định khi r≥rISCOr \ge r_{\text{ISCO}}, ổn định biên tại rISCO=3rs=6GMc2r_{\text{ISCO}} = 3r_s = \dfrac{6GM}{c^2} và không ổn định khi 3rs/2<r<3rs3r_s/2<r<3r_s.

Vì sao đúng?

Hãy nghĩ chuyển động xuyên tâm như một quả bóng lăn trong một giếng thế năng một chiều Veff(r)V_{\text{eff}}(r): quỹ đạo tròn nằm ở nơi giếng phẳng (dVeff/dr=0dV_{\text{eff}}/dr=0), và chúng chỉ ổn định ở nơi giếng cong lên (d2Veff/dr2>0d^2V_{\text{eff}}/dr^2>0) chứ không cong xuống. Vì thế năng tương đối tính có thêm một số hạng 1/r31/r^3 vốn không có trong hấp dẫn Newton, giếng phát triển một cực đại gần lỗ đen — bên trong bán kính đó không hề có quỹ đạo tròn ổn định nào, dù hạt quay nhanh đến đâu; vật chất chỉ đơn giản là rơi thẳng vào.

Chứng minh

Bước 1 (đại lượng bảo toàn và phương trình xuyên tâm). Dọc theo mọi trắc địa (giống thời gian κ=−c2\kappa=-c^2, hoặc null κ=0\kappa=0) nằm trong mặt phẳng xích đạo θ=π/2\theta=\pi/2, việc mêtric không phụ thuộc tt và ϕ\phi cho ta hai đại lượng bảo toàn E~=(1−rs/r)c2t˙\tilde E = (1-r_s/r)c^2\dot t và L~=r2ϕ˙\tilde L = r^2\dot\phi (chấm =d/dτ=d/d\tau). Thay vào điều kiện chuẩn hóa gμνx˙μx˙ν=κg_{\mu\nu}\dot x^\mu\dot x^\nu=\kappa và rút gọn, với chuyển động giống thời gian ta được (drdτ)2=E~2−(1−rsr)(c2+L~2r2)≡E~2−Veff2(r)\left(\frac{dr}{d\tau}\right)^2 = \tilde E^2 - \left(1-\frac{r_s}{r}\right)\left(c^2+\frac{\tilde L^2}{r^2}\right) \equiv \tilde E^2 - V^2_{\text{eff}}(r)

Bước 2 (điều kiện quỹ đạo tròn). Khai triển ra Veff2(r)=c2−c2rsr+L~2r2−L~2rsr3V^2_{\text{eff}}(r) = c^2 - \frac{c^2 r_s}{r} + \frac{\tilde L^2}{r^2} - \frac{\tilde L^2 r_s}{r^3} Một quỹ đạo tròn có rr không đổi, tức dr/dτ=0dr/d\tau=0 mọi lúc, đòi hỏi cả E~2=Veff2(r)\tilde E^2=V^2_{\text{eff}}(r) lẫn (để rr không trôi đi) dVeff2/dr=0dV^2_{\text{eff}}/dr=0: ddrVeff2=c2rsr2−2L~2r3+3L~2rsr4=0\frac{d}{dr}V^2_{\text{eff}} = \frac{c^2 r_s}{r^2} - \frac{2\tilde L^2}{r^3} + \frac{3\tilde L^2 r_s}{r^4} = 0

Bước 3 (giải ra mô men động lượng của quỹ đạo tròn). Nhân với r4r^4 rồi giải theo L~2\tilde L^2 cho ta mô men động lượng cần để duy trì quỹ đạo tròn tại bán kính rr: L~2(r)=c2rsr22r−3rs\tilde L^2(r) = \frac{c^2 r_s r^2}{2r-3r_s} Điều này đã mang thông tin về mặt cầu photon: khi r→(3/2)rsr\to (3/2)r_s từ phía trên, mẫu số 2r−3rs→0+2r-3r_s\to0^+ và L~2→+∞\tilde L^2\to+\infty — không có mô men động lượng hữu hạn nào duy trì được quỹ đạo tròn gần đến vậy, đúng là giới hạn khối lượng bằng không (photon) đạt được ở dưới.

Bước 4 (ổn định biên — ISCO). Khi rr giảm từ vô cực, L~2(r)\tilde L^2(r) trước tiên giảm, đạt cực tiểu, rồi bùng nổ tại r=3rs/2r=3r_s/2; một quỹ đạo tròn ổn định đúng nơi cần L~2\tilde L^2 tăng để rr co lại thêm (tức là dấu của dL~2/drd\tilde L^2/dr quyết định tính ổn định). Lấy đạo hàm kết quả Bước 3, dL~2dr=2c2rs r(r−3rs)(2r−3rs)2\frac{d\tilde L^2}{dr} = \frac{2c^2 r_s\,r(r-3r_s)}{(2r-3r_s)^2} triệt tiêu tại r=0r=0 (loại), r=3rsr=3r_s, và không nơi nào khác với r>3rs/2r>3r_s/2. Nghiệm trong duy nhất này là bán kính ổn định biên: rISCO=3rs=6GMc2r_{\text{ISCO}} = 3r_s = \dfrac{6GM}{c^2} Với r>3rsr>3r_s quỹ đạo ổn định (dL~2/dr>0d\tilde L^2/dr>0, rr lớn hơn cần mô men động lượng lớn hơn, như trong Kepler); với 3rs/2<r<3rs3r_s/2<r<3r_s nó không ổn định.

Bước 5 (mặt cầu photon từ trắc địa null). Với trắc địa null, cùng phép thay thế với κ=0\kappa=0 cho ta (drdλ)2=E2−L2r2(1−rsr)≡E2−Vγ2(r)\left(\frac{dr}{d\lambda}\right)^2 = E^2 - \frac{L^2}{r^2}\left(1-\frac{r_s}{r}\right) \equiv E^2 - V^2_{\gamma}(r) với Vγ2(r)=L2(1r2−rsr3)V^2_\gamma(r) = L^2\left(\frac{1}{r^2}-\frac{r_s}{r^3}\right) Cho đạo hàm của biểu thức trong ngoặc bằng không, ddr(1r2−rsr3)=−2r3+3rsr4=0  ⟹  r=3rs2\frac{d}{dr}\left(\frac{1}{r^2}-\frac{r_s}{r^3}\right) = -\frac{2}{r^3}+\frac{3r_s}{r^4}=0 \;\Longrightarrow\; r=\frac{3r_s}{2} xác định trực tiếp bán kính duy nhất của quỹ đạo tròn photon, rph=32rs=3GMc2r_{\text{ph}} = \dfrac{3}{2}r_s = \dfrac{3GM}{c^2} khớp với giới hạn tìm thấy ở Bước 3, và nó không ổn định vì biểu thức 1/r2−rs/r31/r^2-r_s/r^3 có một cực đại địa phương, không phải cực tiểu, tại đó — bất kỳ xê dịch vào trong nào cũng làm photon xoáy vào chân trời, bất kỳ xê dịch ra ngoài nào cũng làm nó thoát ra vô cực.

Một tia sáng đi qua khối lượng MM với thông số va chạm b≫rsb\gg r_s (điểm gần nhất nằm xa ngoài chân trời) bị lệch, ở bậc thấp nhất, một góc Δϕ=4GMc2b\Delta\phi = \dfrac{4GM}{c^2 b} — đúng bằng hai lần giá trị Newton 2GM/(c2b)2GM/(c^2b) mà một phép tính hạt ánh sáng ngây thơ sẽ dự đoán.

Vì sao đúng?

Hấp dẫn Newton, coi như một lực tác lên hạt chuyển động nhanh, dự đoán có độ lệch, nhưng nó chỉ tính đến độ cong của thời gian (đồng hồ chạy chậm gần khối lượng). Thuyết tương đối rộng thêm vào một đóng góp bằng nhau từ độ cong của không gian (thước co lại theo hướng xuyên tâm gần khối lượng), và hai hiệu ứng cộng lại, làm gấp đôi độ lệch — chính hệ số này đã giúp đoàn thám hiểm nhật thực năm 1919 phân biệt thuyết Einstein với thuyết Newton.

Chứng minh

Bước 1 (phương trình quỹ đạo cho ánh sáng). Với một trắc địa null nằm trong mặt phẳng xích đạo, cùng các đại lượng bảo toàn E=(1−rs/r)c2t˙E=(1-r_s/r)c^2\dot t và L=r2ϕ˙L=r^2\dot\phi dùng cho mặt cầu photon, thay vào gμνx˙μx˙ν=0g_{\mu\nu}\dot x^\mu\dot x^\nu=0 và chuyển từ đạo hàm theo τ\tau sang đạo hàm theo ϕ\phi của u≡1/ru\equiv1/r (dùng r˙=−L du/dϕ\dot r = -L\,du/d\phi), sau khi lấy đạo hàm thêm một lần cho ta phương trình quỹ đạo bẻ cong ánh sáng chuẩn d2udϕ2+u=3rs2u2,u≡1r\frac{d^2u}{d\phi^2}+u = \frac{3r_s}{2}u^2, \qquad u\equiv\frac1r

Bước 2 (bậc không: đường thẳng). Bỏ vế phải (đặt rs=0r_s=0, không-thời gian phẳng) cho u0′′+u0=0u_0''+u_0=0, nghiệm đi qua khoảng cách gần nhất bb tại ϕ=π/2\phi=\pi/2 là đường thẳng u0(ϕ)=sin⁡ϕbu_0(\phi) = \frac{\sin\phi}{b} — quỹ đạo không bị lệch dưới dạng tọa độ cực.

Bước 3 (nhiễu loạn bậc một). Viết u=u0+rsu1u=u_0+r_s u_1 và thay vào phương trình Bước 1, chỉ giữ số hạng tuyến tính theo rsr_s; phần u0′′+u0=0u_0''+u_0=0 triệt tiêu và còn lại là một dao động điều hòa cưỡng bức cho u1u_1: u1′′+u1=32b2sin⁡2ϕu_1''+u_1 = \frac{3}{2b^2}\sin^2\phi Một nghiệm riêng khớp với độ lệch đối xứng mong đợi ở hai bên điểm gần nhất là u1(ϕ)=12b2(1+cos⁡2ϕ)u_1(\phi) = \frac{1}{2b^2}\left(1+\cos^2\phi\right) nên quỹ đạo đầy đủ bậc một là u(ϕ)=sin⁡ϕb+rs2b2(1+cos⁡2ϕ)u(\phi) = \frac{\sin\phi}{b} + \frac{r_s}{2b^2}\left(1+\cos^2\phi\right)

Bước 4 (trích ra tổng góc lệch). Không gian phẳng có u→0u\to0 đúng tại ϕ=0\phi=0 và ϕ=π\phi=\pi; với số hạng sửa rsr_s, u→0u\to0 tại ϕ=−δ1\phi=-\delta_1 và ϕ=π+δ2\phi=\pi+\delta_2 với δ1,δ2\delta_1,\delta_2 nhỏ. Khai triển u(−δ1)=0u(-\delta_1)=0 đến bậc một theo δ1\delta_1 và rsr_s cho −δ1/b+rs/b2=0-\delta_1/b+r_s/b^2=0, nên δ1=rs/b\delta_1=r_s/b; tính tương tự ở đầu kia cho δ2=rs/b\delta_2=r_s/b. Tổng góc mà tiệm cận đi ra lệch khỏi song song ngược hướng hoàn toàn với tiệm cận đi vào là Δϕ=δ1+δ2=2rsb=4GMc2b\Delta\phi = \delta_1+\delta_2 = \frac{2r_s}{b} = \frac{4GM}{c^2 b} mà, khi thay lại rs=2GM/c2r_s=2GM/c^2, chính là công thức cần chứng minh.

Nâng caoKỳ dị là không thể tránh khỏi: định lý Penrose

Kỳ dị r=0r=0 bên trong một lỗ đen Schwarzschild có thể trông như một đặc thù của tính đối xứng cầu tuyệt đối được giả định khi giải phương trình Einstein — có lẽ một sự sụp đổ kém đối xứng hơn, thực tế hơn sẽ tránh được nó. Năm 1965, Roger Penrose chỉ ra điều ngược lại: một khi một mặt bị bẫy hình thành (một mặt kín mà cả hai chùm tia sáng đi ra lẫn đi vào đều hội tụ, không chỉ chùm đi vào như ở xa mọi lỗ đen), một kỳ dị là không thể tránh khỏi, với bất kỳ vật chất nào tuân theo các điều kiện năng lượng hợp lý, bất kể đối xứng. Đây là định lý kỳ dị hiện đại đầu tiên, và nó khởi động các phương pháp toàn cục, hình học chi phối thuyết tương đối toán học ngày nay.

Nếu không-thời gian chứa một mặt Cauchy không compact, tuân theo điều kiện năng lượng null, và chứa một mặt bị bẫy kín, thì không-thời gian là không đầy đủ trắc địa null về tương lai: ít nhất một tia sáng không thể kéo dài vô hạn về tương lai — một kỳ dị, theo nghĩa của Penrose.

Vì sao đúng?

Trong thuyết tương đối rộng, hấp dẫn luôn hội tụ ánh sáng chứ không làm nó phân kỳ (đó là điều điều kiện năng lượng mã hóa); một chùm tia sáng vốn đã hội tụ ở cả hai phía tại một mặt kín nào đó do đó sẽ tiếp tục hội tụ, và về mặt hình học các tia sáng hội tụ phải cắt nhau — nhưng hai tia sáng cắt nhau (một điểm liên hợp) bên trong một vùng mà một không-thời gian hyperbolic toàn cục, có tính chất tốt không thể lý giải được nếu không tạo ra một biên cho sự tồn tại của tia sáng đó. Định lý biến trực giác hội tụ này thành một bất đẳng thức chặt chẽ (phương trình Raychaudhuri) và một chướng ngại tô pô toàn cục.

Chứng minh

Bước 1 (mặt bị bẫy, bằng ký hiệu). Cho SS là một mặt 2 chiều kín, giống không gian, và cho θ\theta ký hiệu độ giãn nở của một chùm trắc địa null rời khỏi SS theo phương vuông góc — tốc độ tương đối mà diện tích tiết diện chùm tăng lên. Xa mọi khối lượng, chùm đi ra giãn nở (θout>0\theta_{\text{out}}>0) và chùm đi vào co lại (θin<0\theta_{\text{in}}<0). SS bị bẫy nếu thay vào đó θout<0,θin<0\theta_{\text{out}}<0, \qquad \theta_{\text{in}}<0 trên toàn bộ SS: ngay cả ánh sáng đi ra cũng đang bị kéo vào trong.

Bước 2 (phương trình Raychaudhuri). Lấy đạo hàm định nghĩa θ=∇aka\theta=\nabla_a k^a của độ giãn nở dọc theo chùm null có tiếp tuyến kak^a, và dùng định nghĩa của tensor Riemann để hoán đổi các đạo hàm, cho ta đẳng thức động học chính xác dθdλ=−θ22−σabσab−Rabkakb\frac{d\theta}{d\lambda} = -\frac{\theta^2}{2} - \sigma_{ab}\sigma^{ab} - R_{ab}k^ak^b trong đó σab\sigma_{ab} là tensor cắt của chùm. Phương trình này thuần túy hình học — một phát biểu về cách mọi họ tia sáng phải bẻ cong, suy ra từ định nghĩa độ cong, trước khi giả định bất kỳ vật lý nào.

Bước 3 (điều kiện năng lượng buộc hội tụ thêm). Vật chất hợp lý về mặt vật lý thỏa điều kiện năng lượng null Rabkakb≥0R_{ab}k^ak^b \ge 0 và số hạng cắt σabσab≥0\sigma_{ab}\sigma^{ab}\ge0 luôn đúng (nó là tổng bình phương); bỏ cả hai số hạng không âm ở vế phải đẳng thức Bước 2 cho ta bất đẳng thức dθdλ≤−θ22\frac{d\theta}{d\lambda} \le -\frac{\theta^2}{2}

Bước 4 (bất đẳng thức Riccati buộc bùng nổ trong thời gian hữu hạn). Trên một mặt bị bẫy, θ0<0\theta_0<0 tại λ=0\lambda=0. Miễn là θ\theta còn âm, chia bất đẳng thức cho θ2>0\theta^2>0 và sắp xếp lại cho ddλ(1θ)≥12  ⟹  1θ(λ)≥1θ0+λ2\frac{d}{d\lambda}\left(\frac{1}{\theta}\right) \ge \frac12 \;\Longrightarrow\; \frac{1}{\theta(\lambda)} \ge \frac{1}{\theta_0}+\frac{\lambda}{2} Vì vế phải tăng không giới hạn khi λ\lambda tăng trong khi 1/θ1/\theta (âm) phải thỏa cận dưới này, 1/θ1/\theta bị buộc tiến lên 0−0^-, tức θ→−∞\theta\to-\infty, tại một tham số afin hữu hạn nào đó λ∗≤2/∣θ0∣\lambda_*\le 2/|\theta_0|: chùm null vuông góc với SS phải phát triển một điểm liên hợp (một điểm caustic, nơi các tia sáng lân cận cắt nhau) trong thời gian afin hữu hạn.

Bước 5 (từ hội tụ cục bộ đến không đầy đủ toàn cục). Một kết quả chuẩn của lý thuyết cấu trúc nhân quả (được dùng nhưng không suy lại ở đây) cho thấy một trắc địa null trong không-thời gian hyperbolic toàn cục không thể duy trì tính achronal — không thể ở lại biên của tương lai của một tập — sau một điểm liên hợp. Kết hợp với mặt Cauchy không compact, điều này loại trừ khả năng trắc địa vòng trở lại hay kết thúc theo cách thông thường, và buộc ít nhất một trắc địa null vuông góc với SS phải không đầy đủ về tương lai: nó không thể kéo dài tới tham số afin lớn tùy ý. Theo định nghĩa, tính không đầy đủ này chính là điều Penrose gọi là kỳ dị — định lý cho thấy nó phải xảy ra mà không cần giả định đối xứng nào, bất cứ khi nào một mặt bị bẫy hình thành.

Đại họcỨng dụng thực tiễn và Ví dụ minh họa

Hình học lỗ đen thôi chỉ thuần túy lý thuyết từ những năm 2010. Kính viễn vọng chân trời sự kiện (EHT), một mạng lưới đĩa radio cỡ hành tinh hoạt động như một máy giao thoa duy nhất, chụp ảnh cái bóng mà một lỗ đen đổ lên vùng khí phát sáng quanh nó — một đĩa tối gần bán kính liên quan tới mặt cầu photon rshadow≈2.6 rsr_{\text{shadow}}\approx2.6\,r_s, không phải bản thân chân trời. Và LIGO–Virgo–KAGRA phát hiện sóng hấp dẫn phát ra khi hai lỗ đen xoáy vào nhau rồi sáp nhập, mã hóa khối lượng của chúng trong sự thay đổi tần số của sóng và biên độ căng cực nhỏ h∼10−21h\sim10^{-21}.

Ví dụ: Đo kích thước bóng của Sagittarius A*

Sagittarius A*, lỗ đen siêu khối lượng ở tâm thiên hà chúng ta, có khối lượng M≈4.0×106 M⊙M\approx4.0\times10^6\,M_\odot và ở khoảng cách D≈8 kpcD\approx8\text{ kpc}. Hãy ước lượng đường kính góc của bóng nó và so sánh với phép đo năm 2022 của kính viễn vọng chân trời sự kiện, khoảng 51.8 μas51.8\,\mu\text{as}.

Lời giải

Bước 1 (bán kính Schwarzschild). Với một khối lượng Mặt Trời, rs(M⊙)=2GM⊙/c2≈2.95 kmr_s(M_\odot)=2GM_\odot/c^2\approx2.95\text{ km}. Tỷ lệ tuyến tính theo khối lượng, rs≈4.0×106×2.95 km≈1.18×107 km=1.18×1010 mr_s\approx4.0\times10^6\times2.95\text{ km}\approx1.18\times10^7\text{ km}=1.18\times10^{10}\text{ m}.

Bước 2 (đường kính bóng). Bán kính bóng là rshadow=332rs≈2.6 rsr_{\text{shadow}} = \dfrac{3\sqrt3}{2}r_s \approx 2.6\,r_s nên đường kính là 2rshadow=2×2.6 rs≈5.2 rs≈6.14×1010 m2r_{\text{shadow}}=2\times2.6\,r_s\approx5.2\,r_s\approx6.14\times10^{10}\text{ m}.

Bước 3 (đổi sang góc). Với góc nhỏ, kích thước góc == kích thước vật lý // khoảng cách. Với D=8 kpc=8×3.086×1019 m=2.47×1020 mD=8\text{ kpc}=8\times3.086\times10^{19}\text{ m}=2.47\times10^{20}\text{ m}, góc là 6.14×1010/2.47×1020≈2.49×10−10 rad6.14\times10^{10}/2.47\times10^{20}\approx2.49\times10^{-10}\text{ rad}.

Bước 4 (đổi radian sang micro giây góc). Một radian là 2.063×1011 μas2.063\times10^{11}\,\mu\text{as}, nên đường kính bóng là 2.49×10−10×2.063×1011 μas≈51.3 μas2.49\times10^{-10}\times2.063\times10^{11}\,\mu\text{as}\approx51.3\,\mu\text{as} — khớp với giá trị EHT đo được 51.8 μas51.8\,\mu\text{as} trong phạm vi 1%1\%, một xác nhận ấn tượng rằng Sgr A* hoạt động đúng như hình học Schwarzschild (chính xác hơn, Kerr với spin gần bằng không) dự đoán.

Ví dụ: Khối lượng chirp và bậc độ lớn của căng cho GW150914

Phát hiện đầu tiên của LIGO, GW150914, đến từ hai lỗ đen khoảng m1≈36 M⊙m_1\approx36\,M_\odot và m2≈29 M⊙m_2\approx29\,M_\odot sáp nhập ở khoảng cách D≈410 MpcD\approx410\text{ Mpc}. Hãy tính khối lượng chirp M\mathcal M điều khiển sự quét tần số của sóng, và ước lượng bậc độ lớn của căng hh tại Trái Đất.

Lời giải

Bước 1 (khối lượng chirp). M=(m1m2)3/5(m1+m2)1/5\mathcal M = \dfrac{(m_1 m_2)^{3/5}}{(m_1+m_2)^{1/5}} Với m1m2≈1044 M⊙2m_1m_2\approx1044\,M_\odot^2 và m1+m2≈65 M⊙m_1+m_2\approx65\,M_\odot: (m1m2)3/5≈10440.6 M⊙1.2≈100 M⊙1.2(m_1m_2)^{3/5}\approx1044^{0.6}\,M_\odot^{1.2}\approx100\,M_\odot^{1.2} và (m1+m2)1/5≈650.2≈2.3(m_1+m_2)^{1/5}\approx65^{0.2}\approx2.3, cho M≈100/2.3 M⊙≈28 M⊙\mathcal M\approx100/2.3\,M_\odot\approx28\,M_\odot — rất gần với giá trị LIGO báo cáo khoảng 28 M⊙28\,M_\odot.

Bước 2 (thang độ dài của nguồn). Bán kính Schwarzschild của tổng khối lượng là rs,tot=2G(m1+m2)/c2≈65×2.95 km≈1.9×105 mr_{s,\text{tot}}=2G(m_1+m_2)/c^2\approx65\times2.95\text{ km}\approx1.9\times10^5\text{ m}: kích thước của hệ đang sáp nhập gần lúc hợp nhất.

Bước 3 (căng theo bậc độ lớn). Gần lúc sáp nhập tốc độ quỹ đạo tiến gần một phần đáng kể của cc; lấy v/c∼0.5v/c\sim0.5 và dùng tỷ lệ h∼rs,totD(vc)2h \sim \dfrac{r_{s,\text{tot}}}{D}\left(\dfrac{v}{c}\right)^2 với D=410 Mpc≈1.26×1025 mD=410\text{ Mpc}\approx1.26\times10^{25}\text{ m}: h∼(1.9×105/1.26×1025)×0.25≈3.8×10−21h\sim(1.9\times10^5/1.26\times10^{25})\times0.25\approx3.8\times10^{-21}.

Bước 4 (so sánh). Ước lượng bậc độ lớn này rơi trong khoảng vài lần so với căng LIGO thực đo được, h∼10−21h\sim10^{-21} — con số nhỏ bé này không phải dấu hiệu của một hiệu ứng yếu tại nguồn (gần lúc sáp nhập, không-thời gian ở đó cong mạnh gần như mức vật lý cho phép) mà chỉ đơn giản là sự pha loãng 1/D1/D của bức xạ hấp dẫn qua hàng trăm megaparsec.

Nâng caoLỗ đen quay: hình học Kerr và vùng ergosphere

Các lỗ đen thiên văn thực sự sinh ra từ các ngôi sao quay hoặc từ các vụ sáp nhập và mang mô men động lượng JJ, được tham số hóa bởi độ dài spin a=J/(Mc)a=J/(Mc). Roy Kerr tìm ra nghiệm chân không dừng, đối xứng trục duy nhất vào năm 1963. Sự quay tách mặt Schwarzschild duy nhất r=rsr=r_s thành hai mặt riêng biệt: một chân trời sự kiện ngoài tại r+=GM/c2+(GM/c2)2−a2r_+ = GM/c^2 + \sqrt{(GM/c^2)^2-a^2}, và bên ngoài nó là một mặt giới hạn dừng lớn hơn, dẹt ở hai cực rergo(θ)=GMc2+(GMc2)2−a2cos⁡2θr_{\text{ergo}}(\theta) = \dfrac{GM}{c^2} + \sqrt{\left(\dfrac{GM}{c^2}\right)^2 - a^2\cos^2\theta} Chạm vào chân trời ở hai cực (θ=0,π\theta=0,\pi) và phình ra tới r=rsr=r_s ở xích đạo (θ=π/2\theta=\pi/2), mặt này bao quanh ergosphere — một vùng bên ngoài chân trời nơi hiệu ứng kéo hệ quy chiếu mạnh đến mức không người quan sát nào có thể đứng yên tại (r,θ,ϕ)(r,\theta,\phi) cố định so với các ngôi sao ở xa mà không chuyển động nhanh hơn ánh sáng: bản thân không-thời gian đang bị cuốn xoay theo chiều spin của lỗ đen, và mọi thứ trong ergosphere buộc phải quay cùng chiều với nó. Vì ergosphere nằm ngoài chân trời, các hạt có thể đi vào đó rồi vẫn thoát ra vô cực, và Penrose (1969) chỉ ra rằng tách một hạt bên trong ergosphere cho phép mảnh thoát ra mang đi nhiều năng lượng hơn hạt ban đầu mang vào — khai thác năng lượng quay của lỗ đen.

Nâng caoGợn sóng của không-thời gian: sóng hấp dẫn tuyến tính hóa

Ở xa nguồn, nơi không-thời gian gần như phẳng, mêtric có thể viết như một nhiễu loạn nhỏ của mêtric Minkowski, gμν=ημν+hμνg_{\mu\nu}=\eta_{\mu\nu}+h_{\mu\nu} với ∣hμν∣≪1|h_{\mu\nu}|\ll1. Trong chuẩn điều hòa (Lorenz) ∂μhˉμν=0\partial^\mu\bar h_{\mu\nu}=0, với hˉμν=hμν−12ημνh\bar h_{\mu\nu}=h_{\mu\nu}-\tfrac12\eta_{\mu\nu}h là nhiễu loạn đảo vết, phương trình trường Einstein tuyến tính hóa thành một phương trình sóng chuẩn □hˉμν=−16πGc4Tμν\Box \bar h_{\mu\nu} = -\dfrac{16\pi G}{c^4}T_{\mu\nu} với □=−c−2∂t2+∇2\Box=-c^{-2}\partial_t^2+\nabla^2: các nhiễu loạn của mêtric lan truyền với tốc độ ánh sáng cc, có nguồn là tensor năng lượng–xung lượng TμνT_{\mu\nu} hệt như sóng điện từ có nguồn là dòng điện. Vì khối lượng và động lượng được bảo toàn, mô men đơn cực và lưỡng cực của một hệ cô lập không thể bức xạ; bức xạ bậc thấp nhất đến từ đạo hàm bậc hai theo thời gian — và trong công suất, đạo hàm bậc ba theo thời gian — của mô men tứ cực khối lượng QijQ_{ij}, cho ta công thức tứ cực của Einstein cho công suất bức xạ: P=G5c5⟨Q...ijQ...ij⟩P = \dfrac{G}{5c^5}\left\langle \dddot Q_{ij}\dddot Q^{ij}\right\rangle Hệ số đứng trước G/c5G/c^5 cực nhỏ (∼3×10−53 W−1\sim3\times10^{-53}\text{ W}^{-1}), đó là lý do chỉ những hệ tương đối tính, cỡ khối lượng Mặt Trời trở lên như các lỗ đen hay sao neutron đang sáp nhập mới bức xạ đủ mạnh để phát hiện được.

Nghiên cứuCâu hỏi mở và hướng nghiên cứu đang sống

Nếu khối lượng MM của một lỗ đen tăng gấp ba, bán kính Schwarzschild rsr_s, bán kính mặt cầu photon rphr_{\text{ph}}, và bán kính quỹ đạo tròn ổn định trong cùng rISCOr_{\text{ISCO}} của nó thay đổi thế nào?

Tại sao ta biết rằng r=rsr=r_s trong mêtric Schwarzschild là một kỳ dị tọa độ chứ không phải một kỳ dị độ cong vật lý?

Điều gì quyết định kích thước cái "bóng" tối mà Kính viễn vọng Chân trời Sự kiện chụp được quanh M87 và Sagittarius A?

Tại sao độ căng sóng hấp dẫn mà LIGO đo được trên Trái Đất lại nhỏ đến vậy (h∼10−21h\sim10^{-21}), dù một vụ sáp nhập lỗ đen đôi trong khoảnh khắc phát ra công suất lớn hơn toàn bộ các ngôi sao trong vũ trụ khả kiến cộng lại?

Tài liệu tham khảo

  1. Charles W. Misner, Kip S. Thorne, John Archibald Wheeler (1973). Gravitation
  2. Roger Penrose (1965). Gravitational Collapse and Space-Time Singularities · DOI:10.1103/PhysRevLett.14.57
  3. B. P. Abbott et al. (LIGO Scientific Collaboration and Virgo Collaboration) (2016). Observation of Gravitational Waves from a Binary Black Hole Merger · DOI:10.1103/PhysRevLett.116.061102
  4. Event Horizon Telescope Collaboration (2022). First Sagittarius A* Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole in the Center of the Milky Way · DOI:10.3847/2041-8213/ac6674
  5. Sergiu Klainerman, Jérémie Szeftel (2021). Kerr stability for small angular momentum · arXiv:2104.11857 [preprint, chưa bình duyệt]
  6. Geoffrey Penington (2019). Entanglement Wedge Reconstruction and the Information Paradox · arXiv:1905.08255 [preprint, chưa bình duyệt]