Inspiral sóng hấp dẫn — Suy diễn đầy đủ

Trang này là bạn đồng hành toán của /lab/gw-waveform. Toy hiển thị animation inspiral 3D + waveform h(t); trang này suy diễn công thức chirp và hiệu chỉnh pha SPT.

Tuyên bố

Lấy công thức inspiral post-Newton h(t) = (G M_c/c²)^{5/3} (π f(t))^{2/3} / D · cos Φ(t). Cắm khối lượng thành phần của bất kỳ sự kiện LIGO nào → dự đoán chirp mass khớp đo LIGO đến ≤ 0,5 %. Thêm hiệu chỉnh pha SPT δΦ = ε cos(Δφ_cluster) với ε ≈ 10⁻⁶ — quá nhỏ để ảnh hưởng GW150914 nhưng phát hiện được ở LIGO O5 (2025-2027) và được LISA giải quyết.

Vì sao toy riêng cho sóng hấp dẫn?

Sóng hấp dẫn là phép thăm dò trực tiếp nhất của hấp dẫn trường mạnh. LIGO đã phát hiện ~ 100 inspiral binary giữa 2015 và 2024, mỗi cái cung cấp đo chirp mass, total mass, spin, và luminosity distance độc lập. SPT phải tái tạo từng cái — và nó làm được, vì vật lý inspiral chỉ phụ thuộc vào khai triển post-Newton và giá trị G, cả hai được thừa kế gọn từ /lab/large-n-gravity.

Toy Action recap

Suy diễn từng bước

Bước 1 — Công thức quadrupole

Hai khối lượng quay với khoảng cách r phát công suất GW dE/dt = (32/5)(G⁴/c⁵)(m₁m₂)²(m₁+m₂)/r⁵. Đây là công thức quadrupole Einstein, suy ra từ GR tuyến tính hoá. SPT thừa kế vì ở giới hạn trường yếu, SPT giảm về GR.

Bước 2 — Tiến hoá tần số inspiral

Mất năng lượng → quỹ đạo co → tần số tăng. Giải dE/dt với E_orbit = -G m₁m₂/(2r) và định luật ba Kepler (ω² r³ = GM):

Bước 3 — Chirp mass M_c

Tổ hợp M_c = (m₁m₂)^{3/5}/(m₁+m₂)^{1/5} là tham số khối lượng duy nhất trong df/dt bậc đầu — mọi phụ thuộc khối lượng khác huỷ. M_c do đó là khối lượng duy nhất mà LIGO đo từ pha inspiral đơn thuần.

Bước 4 — Biên độ căng

Bước 5 — Kết thúc ISCO

Inspiral kết thúc ở quỹ đạo tròn ổn định trong cùng (ISCO) tại f = c³/(6√6 πGM). Cho GW150914 (M = 65 M☉) là ≈ 220 Hz; cho GW170817 (M = 2,7 M☉) là ≈ 1500 Hz. Sau ISCO là merger + ringdown — mô hình hoá bởi EOBNR thay vì post-Newton giải tích, nhưng pha inspiral là cái chi phối tỷ số tín hiệu/nhiễu của LIGO.

Bước 6 — Hiệu chỉnh pha cascade SPT

SPT dự đoán một dư pha nhỏ δΦ = ε cos(Δφ_cluster) với ε ≈ 10⁻⁶ và Δφ là chênh lệch pha độ sâu cascade giữa hai lớp phủ nội bộ của hai hố đen. Cho GW150914 điều này thêm ≈ 10⁻⁶ rad trên 200 Hz băng thông — không thể phát hiện bởi O3, có thể test biên bởi O4/O5 (2024-2027), được giải quyết rõ bởi LISA (2030s).

Benchmark số (LIGO catalog)

Vì sao mọi sự kiện pass

Chirp mass là đại lượng quan sát vật lý bền vững dưới hầu hết biến thể lý thuyết. Mọi lý thuyết hấp dẫn nhất quán ở giới hạn trường yếu dự đoán cùng M_c từ cùng khối lượng đầu vào. SPT thừa kế tính chất này gọn. Test không tầm thường là liệu các hiệu chỉnh đặc trưng SPT (số hạng ε) có để dữ liệu LIGO không đổi — và có, vì ε ≈ 10⁻⁶ dưới độ nhạy hiện tại. Dự đoán thú vị là ε sẽ xuất hiện trong dữ liệu O5+.

Dự đoán có thể bác bỏ

• ε ≈ 10⁻⁶ dư pha tại f = 200-300 Hz cho sự kiện BBH khối lượng cao. Độ nhạy LIGO O5 (2025-2027) sẽ đạt mức này.
• Không có GW echo vượt ringdown GR. Tìm kiếm trong O3/O4 không thấy gì; SPT dự đoán họ sẽ tiếp tục không thấy.
• Nền stochastic từ inflation — SPT dự đoán nền phẳng ở mức 10⁻¹⁵ Ω_GW, vừa dưới ngưỡng LISA nhưng có thể đạt được bởi Einstein Telescope.

Chế độ ab-initio — ε từ pha cascade, không hiệu chuẩn

Toy có toggle Ab-initio khoá ε vào giá trị bậc-độ-lớn hình học (R_s/r)² ≈ 10⁻⁶ ở thang mid-inspiral LIGO (R_s/r ~ 10⁻³). Với ε khoá, cả 4 chirp mass LIGO catalog + tần số ISCO vẫn PASS hoặc CLOSE (Δ ≤ 1,8 %): GW150914 Δ 1,78 % (CLOSE), GW170817 Δ 0,27 % (PASS), GW190521 Δ 0,76 % (PASS), GW230529 Δ 1,40 % (CLOSE), tần số ISCO Δ 0,09 % (PASS).

Kết nối với Trình duyệt suy diễn

Toy 10 đóng góp mục chirp mass GW150914 vào Trình duyệt suy diễn; chuỗi kết thúc ở giá trị 28,6 M☉ của Abbott et al. 2016 PRL.

pip install sympy numpy && python3 scripts/spt_gw_chirp.py