Vật lý toán

Cơ học thiên thể: bài toán hai vật và ba vật

Cách định luật hấp dẫn Newton đưa chuyển động của hai vật quay quanh nhau về một quỹ đạo conic duy nhất, vì sao hành tinh quét diện tích bằng nhau trong thời gian bằng nhau, và vì sao thêm vật thứ ba phá vỡ bức tranh gọn gàng đó, mở ra hỗn loạn.

Trực giácHai vật rơi quanh nhau

Thả một quả bóng, nó rơi thẳng xuống. Ném nó theo phương ngang đủ nhanh, nó vẫn rơi — nhưng mặt đất cong đi bên dưới đúng bằng tốc độ đó, nên nó không bao giờ chạm đất. Đó chính xác là điều Mặt Trăng làm quanh Trái Đất, và Trái Đất làm quanh Mặt Trời: một cú rơi liên tục dưới tác dụng của hấp dẫn, kết hợp chuyển động ngang đủ nhanh để luôn hụt vật trung tâm. Cơ học thiên thể là toán học mô tả chính xác cách các vật rơi quanh nhau như vậy, từ một hành tinh quanh một ngôi sao cho tới một tàu vũ trụ len lỏi qua nhiều thế giới chỉ nhờ lực hấp dẫn.

Kéo

\theta

: cùng một góc cực đánh dấu một điểm trên vòng tròn cũng đánh dấu vị trí hành tinh trên quỹ đạo. Trên vòng tròn,

r

không đổi khi

\theta

quay; trên quỹ đạo thực,

r(\theta) = \dfrac{p}{1+e\cos\theta},\qquad p = \dfrac{h^{2}}{GM}

khiến

r

dao động giữa điểm gần nhất (

\theta=0

) và điểm xa nhất (

\theta=180^\circ

Phổ thôngĐiều Kepler đã quan sát được

Dựa trên hàng chục năm dữ liệu quan sát hành tinh bằng mắt thường của Tycho Brahe và không hề dùng giải tích, Johannes Kepler công bố ba quy luật thuần túy thực nghiệm khoảng năm 1609–1619: (1) hành tinh chuyển động trên elip với Mặt Trời ở một tiêu điểm, không phải ở tâm; (2) đoạn thẳng nối Mặt Trời và hành tinh quét những diện tích bằng nhau trong những khoảng thời gian bằng nhau, nên hành tinh chuyển động nhanh nhất khi gần Mặt Trời và chậm nhất khi xa; (3) bình phương chu kỳ quỹ đạo của hành tinh tỉ lệ với lập phương bán trục lớn của quỹ đạo. Kepler không có lời giải thích vì sao — điều đó phải chờ gần bảy mươi năm sau, tới định luật hấp dẫn và giải tích của Newton. Phần còn lại của bài này sẽ chứng minh cả ba định luật, và hơn thế, chỉ từ một định luật nghịch đảo bình phương duy nhất đó.

Đại họcĐưa hai vật về một vật

Định nghĩa: Khối lượng rút gọn và vị trí tương đối

Với hai chất điểm khối lượng $m_1$ , $m_2$ chỉ hút nhau bằng hấp dẫn, gọi $\mathbf r = \mathbf r_1 - \mathbf r_2$ là vectơ từ vật 2 tới vật 1, $M = m_1 + m_2$ là tổng khối lượng, và $\mu = \dfrac{m_1 m_2}{m_1+m_2}$ là khối lượng rút gọn. Trong hệ quy chiếu khối tâm, bài toán hai vật thật sự tương đương chính xác với một hạt hư cấu duy nhất khối lượng $\mu$ chuyển động trong trường nghịch đảo bình phương cố định sinh bởi khối lượng $M$ đặt tại gốc — vectơ tương đối $\mathbf r$ vẽ nên quỹ đạo của hạt hư cấu này, còn vị trí thật $\mathbf r_1,\mathbf r_2$ được khôi phục từ $\mathbf r$ bằng phép co giãn theo tỉ lệ khối lượng đơn giản.

\ddot{\mathbf r} = -\dfrac{GM}{r^{3}}\mathbf r

Ở đây $G$ là hằng số hấp dẫn Newton, $M=m_1+m_2$ , $\mathbf r$ là vectơ tách rời tương đối và $r=|\mathbf r|$ là khoảng cách giữa hai vật. Chú ý khối lượng rút gọn $\mu$ đã hoàn toàn triệt tiêu khỏi phương trình này: mọi vật rơi tự do, nhẹ hay nặng, đều theo cùng một quỹ đạo tương đối với cùng điều kiện đầu — hệ quả trực tiếp của quan sát của Galileo rằng mọi vật rơi với cùng tốc độ.

\mathbf h = \mathbf r \times \dot{\mathbf r}

Vectơ $\mathbf h = \mathbf r \times \dot{\mathbf r}$ là mômen động lượng riêng (mômen động lượng trên một đơn vị khối lượng rút gọn). Vì hấp dẫn là lực xuyên tâm — $\ddot{\mathbf r}$ luôn hướng dọc theo $\mathbf r$ — nên nó không sinh mômen lực quanh gốc tọa độ, khiến $\mathbf h$ trở thành ứng viên tự nhiên cho một đại lượng bảo toàn. Định lý 1 dưới đây sẽ chứng minh chặt chẽ điều đó, và chỉ riêng sự thật ấy buộc quỹ đạo phải nằm trong một mặt phẳng cố định và quét diện tích với tốc độ không đổi.

Nâng caoHình dạng quỹ đạo: phương trình Binet và các thiết diện conic

Để tìm hình dạng quỹ đạo, làm việc trong mặt phẳng quỹ đạo với tọa độ cực $(r,\theta)$ và đặt $u=1/r$ . Dùng $\dot\theta=hu^2$ (từ $h = r^{2}\dot\theta$ ), quy tắc dây chuyền cho $\dot r = \dfrac{dr}{d\theta}\dot\theta = \dfrac{d(1/u)}{d\theta}\,hu^{2} = -h\dfrac{du}{d\theta}$ , và đạo hàm thêm lần nữa, $\ddot r = -h\dfrac{d^{2}u}{d\theta^{2}}\dot\theta = -h^{2}u^{2}\dfrac{d^{2}u}{d\theta^{2}}$ .

\dfrac{d^{2}u}{d\theta^{2}} + u = \dfrac{GM}{h^{2}},\qquad u = \dfrac1r

Thay $\ddot r$ và $\dot\theta=hu^2$ vào phương trình chuyển động hướng kính $\ddot r - r\dot\theta^{2} = -GMu^{2}$ rồi chia cho $-h^{2}u^{2}$ cho đúng phương trình quỹ đạo Binet ở trên — một phương trình vi phân tuyến tính hệ số hằng cho $u(\theta)$ . Nghiệm tổng quát là $u(\theta) = \dfrac{GM}{h^{2}} + C\cos(\theta-\theta_0)$ ; chọn $\theta_0=0$ (hướng cận điểm) và viết hằng số tích phân là $C=\dfrac{GM}{h^{2}}e$ cho $r(\theta) = \dfrac{p}{1+e\cos\theta},\qquad p = \dfrac{h^{2}}{GM}$ — phương trình cực của một thiết diện conic với khối lượng hút nằm tại một tiêu điểm, không phải tâm. Tùy $e$ : $0\le e<1$ cho elip (quỹ đạo bị chặn, $e=0$ là đường tròn), $e=1$ cho parabol, và $e>1$ cho hyperbol.

Thiết diện conic theo độ lệch tâm
Độ lệch tâm	Loại quỹ đạo	Năng lượng riêng	Ví dụ
$0 \le e < 1$	Elip (đường tròn nếu $e=0$ )	Âm (bị chặn)	Trái Đất quanh Mặt Trời, $e\approx0.017$
$e = 1$	Parabol	Bằng không	Sao chổi chu kỳ dài cận biên
$e > 1$	Hyperbol	Dương (không bị chặn)	Vật thể liên sao ʻOumuamua

Phương trình quỹ đạo cho biết hình dạng, nhưng riêng năng lượng và mômen động lượng chưa rõ ràng cố định hướng của cận điểm. Điều đáng kinh ngạc là hấp dẫn Newton còn giấu thêm một vectơ bảo toàn nữa: vectơ Laplace–Runge–Lenz $\mathbf A = \dot{\mathbf r}\times\mathbf h - GM\hat{\mathbf r}$ , với $\hat{\mathbf r}=\mathbf r/r$ . Định lý 3 dưới đây chứng minh $\mathbf A$ không đổi và cho $|\mathbf A| = GMe$ , hướng từ tiêu điểm thẳng tới cận điểm — một đối xứng "ẩn" chỉ riêng có ở lực nghịch đảo bình phương (chỉ chia sẻ với dao động điều hòa trong các lực xuyên tâm), phản ánh một đối xứng quay bốn chiều tiềm ẩn của bài toán Kepler.

Đại họcChứng minh ba kết quả cổ điển

Định lý: Định luật hai Kepler (diện tích bằng nhau trong thời gian bằng nhau)

Dưới tác dụng của bất kỳ lực xuyên tâm nào (lực luôn hướng dọc theo $\mathbf r$ ), vectơ bán kính từ tâm lực tới vật chuyển động quét diện tích với tốc độ không đổi $\dfrac{dA}{dt} = \dfrac12 r^{2}\dot\theta = \dfrac{h}{2}$ .

Vì sao đúng?

Đây là định luật hai của Kepler, và chứng minh dưới đây cho thấy nó không liên quan gì riêng tới dạng nghịch đảo bình phương của hấp dẫn — nó chỉ đến từ việc lực là xuyên tâm (song song với vectơ vị trí), nên áp dụng như nhau cho bất kỳ lực xuyên tâm nào, dù có phải hấp dẫn hay không.

Chứng minh

Định nghĩa $\mathbf h = \mathbf r\times\dot{\mathbf r}$ . Đạo hàm, $\dot{\mathbf h} = \dot{\mathbf r}\times\dot{\mathbf r} + \mathbf r\times\ddot{\mathbf r} = \mathbf 0 + \mathbf r\times\ddot{\mathbf r}$ (số hạng đầu triệt tiêu vì một vectơ nhân có hướng với chính nó luôn bằng không). Với lực xuyên tâm, $\ddot{\mathbf r}$ song song với $\mathbf r$ — viết $\ddot{\mathbf r} = f(r)\hat{\mathbf r}$ với hàm vô hướng $f$ nào đó — nên $\mathbf r\times\ddot{\mathbf r} = f(r)\,\mathbf r\times\hat{\mathbf r} = \mathbf 0$ luôn. Vậy $\dot{\mathbf h} = \mathbf r\times\ddot{\mathbf r} = -\dfrac{GM}{r^{3}}(\mathbf r\times\mathbf r) = \mathbf 0$ : $\mathbf h$ là một vectơ cố định, không đổi cả về độ lớn lẫn hướng.

Vì $\mathbf h$ có hướng cố định và $\mathbf h = \mathbf r\times\dot{\mathbf r}$ luôn vuông góc với $\mathbf r$ , vectơ vị trí $\mathbf r$ mãi mãi bị giới hạn trong một mặt phẳng cố định duy nhất vuông góc với $\mathbf h$ — chuyển động là phẳng. Thiết lập tọa độ cực $(r,\theta)$ trong mặt phẳng đó. Viết $\mathbf r = r\hat{\mathbf r}$ và $\dot{\mathbf r} = \dot r\hat{\mathbf r} + r\dot\theta\hat{\boldsymbol\theta}$ , tích có hướng cho $\mathbf h = r\hat{\mathbf r}\times(\dot r\hat{\mathbf r}+r\dot\theta\hat{\boldsymbol\theta}) = r^{2}\dot\theta\,\hat{\mathbf z}$ , nên đại lượng vô hướng $h=r^{2}\dot\theta$ tự nó không đổi.

Trong thời gian $dt$ , vectơ bán kính quét một hình quạt mỏng gần như tam giác có diện tích $dA = \tfrac12 r\cdot(r\,d\theta) = \tfrac12 r^{2}\,d\theta$ (đáy $r\,d\theta$ , chiều cao $r$ , hệ số $\tfrac12$ của tam giác). Chia cho $dt$ cho đúng $\dfrac{dA}{dt} = \dfrac12 r^{2}\dot\theta = \dfrac{h}{2}$ . Vì $h$ không đổi, tốc độ quét diện tích $dA/dt$ không đổi theo thời gian — diện tích bằng nhau được quét trong thời gian bằng nhau, chứng minh xong định luật hai Kepler. $\blacksquare$

Định lý: Định luật ba Kepler

Với quỹ đạo bị chặn (elip) có bán trục lớn $a$ quanh tổng khối lượng $M$ , chu kỳ quỹ đạo $T$ thỏa mãn $T^{2} = \dfrac{4\pi^{2}}{GM}a^{3}$ .

Vì sao đúng?

Điều này liên hệ trực tiếp kích thước quỹ đạo với thời gian hoàn thành nó, không phụ thuộc vào độ lệch tâm — sự thật được dùng mỗi khi nhà thiên văn cân một ngôi sao hay hành tinh bằng cách đo thời gian quỹ đạo của một vật nhỏ hơn quanh nó.

Chứng minh

Theo Định lý 1, tốc độ quét diện tích $dA/dt=h/2$ không đổi, nên lấy tích phân trên một chu kỳ đầy đủ $T$ cho tổng diện tích bao quanh $A = \tfrac{h}{2}T$ . Với elip bán trục lớn $a$ và bán trục nhỏ $b$ , hình học cho $A=\pi ab$ . Cân bằng hai biểu thức: $A = \pi a b = \dfrac{h}{2}T$ .

Tiếp theo, liên hệ $h$ với hình dạng qua phương trình Binet: bán trực thông của quỹ đạo là $p=h^{2}/GM$ , và với elip có các hệ thức chuẩn $b=a\sqrt{1-e^{2}}$ và $p=a(1-e^{2})=b^{2}/a$ . Kết hợp $p=h^{2}/GM$ với $p=b^{2}/a$ cho $h^{2}=\dfrac{GMb^{2}}{a}$ , nên $h = \sqrt{\dfrac{GMb^{2}}{a}} = b\sqrt{\dfrac{GM}{a}}$ .

Thay biểu thức này của $h$ vào $\pi ab=\tfrac{h}{2}T$ và giải ra $T$ : $T=\dfrac{2\pi ab}{h}=\dfrac{2\pi ab}{b\sqrt{GM/a}}=2\pi a\sqrt{\dfrac{a}{GM}}=2\pi\sqrt{\dfrac{a^{3}}{GM}}$ . Bình phương hai vế cho $T^{2} = \dfrac{4\pi^{2}}{GM}a^{3}$ , chính là định luật ba Kepler — và đúng với mọi elip bất kể $e$ , vì $e$ đã triệt tiêu hoàn toàn. $\blacksquare$

Định lý: Bảo toàn vectơ Laplace–Runge–Lenz

Vectơ $\mathbf A = \dot{\mathbf r}\times\mathbf h - GM\hat{\mathbf r}$ không đổi theo thời gian với chuyển động dưới lực nghịch đảo bình phương $\ddot{\mathbf r} = -\dfrac{GM}{r^{3}}\mathbf r$ , có độ lớn $|\mathbf A| = GMe$ , và hướng từ tiêu điểm tới cận điểm.

Vì sao đúng?

Riêng năng lượng và mômen động lượng cố định kích thước và hình dạng của quỹ đạo Kepler nhưng không cố định hướng của nó trong mặt phẳng; vectơ Laplace–Runge–Lenz là đại lượng bảo toàn thêm giúp ghim vị trí cận điểm, và chính sự tồn tại của nó là đặc thù riêng của lực $1/r^2$ — hầu hết lực xuyên tâm khác không có vectơ như vậy, và đó chính là lý do vì sao hầu hết quỹ đạo dưới lực xuyên tâm (như trong bài toán ba vật hạn chế dưới đây) không nằm cố định trên một đường cong khép kín.

Chứng minh

Đạo hàm: $\dfrac{d\mathbf A}{dt} = \ddot{\mathbf r}\times\mathbf h + \dot{\mathbf r}\times\dot{\mathbf h} - GM\dfrac{d\hat{\mathbf r}}{dt}$ . Vì $\mathbf h$ không đổi (Định lý 1), số hạng giữa triệt tiêu. Dùng $\ddot{\mathbf r}=-\dfrac{GM}{r^{2}}\hat{\mathbf r}$ và $\mathbf h=\mathbf r\times\dot{\mathbf r}$ , đồng nhất thức tích ba vectơ $\hat{\mathbf r}\times(\mathbf r\times\dot{\mathbf r}) = \mathbf r(\hat{\mathbf r}\cdot\dot{\mathbf r}) - \dot{\mathbf r}(\hat{\mathbf r}\cdot\mathbf r)$ cho $\ddot{\mathbf r}\times\mathbf h = -\dfrac{GM}{r^{2}}\big[r\dot r\,\hat{\mathbf r} - r\dot{\mathbf r}\big] = GM\left(\dfrac{\dot{\mathbf r}}{r} - \dfrac{\dot r}{r}\hat{\mathbf r}\right)$ , với $\dot r=\hat{\mathbf r}\cdot\dot{\mathbf r}$ là tốc độ thay đổi vô hướng của khoảng cách $r$ (đã dùng $\hat{\mathbf r}\cdot\mathbf r=r$ ).

Mặt khác, đạo hàm trực tiếp $\hat{\mathbf r}=\mathbf r/r$ cho $\dfrac{d\hat{\mathbf r}}{dt} = \dfrac{\dot{\mathbf r}}{r} - \dfrac{\mathbf r\,\dot r}{r^{2}} = \dfrac{\dot{\mathbf r}}{r} - \dfrac{\dot r}{r}\hat{\mathbf r}$ — đúng bằng biểu thức trong ngoặc tìm được cho $\ddot{\mathbf r}\times\mathbf h/GM$ ở trên. Vậy $\ddot{\mathbf r}\times\mathbf h = GM\dfrac{d\hat{\mathbf r}}{dt}$ một cách đồng nhất.

Do đó $\dfrac{d\mathbf A}{dt} = GM\dfrac{d\hat{\mathbf r}}{dt} - GM\dfrac{d\hat{\mathbf r}}{dt} = \mathbf 0$ , chứng minh $\mathbf A$ không đổi. Tính $\mathbf A$ tại cận điểm, nơi $\dot{\mathbf r}$ hoàn toàn tiếp tuyến (vuông góc $\hat{\mathbf r}$ ) với tốc độ $v_p=h/r_p$ , cho $\dot{\mathbf r}\times\mathbf h$ hướng dọc $-\hat{\mathbf r}$ với độ lớn $h^{2}/r_p$ , nên $\mathbf A = \left(\dfrac{h^{2}}{r_p}-GM\right)\hat{\mathbf r}$ ; dùng $r_p=p/(1+e)=\dfrac{h^{2}}{GM(1+e)}$ từ phương trình quỹ đạo cho $\mathbf A = GMe\,\hat{\mathbf r}_p$ , xác nhận $|\mathbf A|=GMe$ hướng đúng về cận điểm. $\blacksquare$

Nâng caoVượt ra ngoài hai vật: bài toán ba vật hạn chế

Thêm một vật thứ ba khối lượng không đáng kể (tàu vũ trụ, tiểu hành tinh) chuyển động dưới hấp dẫn của hai khối lượng lớn $M_1,M_2$ quay quanh nhau theo quỹ đạo tròn. Trong hệ quy chiếu quay cùng hai khối lượng lớn với tốc độ góc quỹ đạo $\omega$ , chuyển động của vật nhỏ chịu chi phối bởi hấp dẫn cộng với lực ly tâm và lực Coriolis hư cấu của hệ quay, kết hợp thành thế hiệu dụng dưới đây.

\Omega(x,y) = -\dfrac{GM_1}{r_1} - \dfrac{GM_2}{r_2} - \dfrac12\omega^{2}(x^{2}+y^{2})

Định nghĩa: Năm điểm Lagrange

Năm điểm nơi lực hiệu dụng (gradient của $\Omega$ , kể cả số hạng Coriolis cho vật đứng yên trong hệ quay) triệt tiêu là các điểm Lagrange. $L_1$ nằm giữa hai khối lượng; $L_2$ nằm ngay ngoài khối lượng nhỏ hơn, phía xa $M_1$ ; $L_3$ nằm phía đối diện của $M_1$ so với $M_2$ . Với $m_2\ll M_1$ (ví dụ Trái Đất quanh Mặt Trời), cả $L_1$ và $L_2$ đều cách khối lượng nhỏ hơn xấp xỉ $r \approx a\left(\dfrac{m_2}{3M_1}\right)^{1/3}$ , với $a$ là khoảng cách quỹ đạo. $L_4$ và $L_5$ tạo thành tam giác đều với $M_1$ và $M_2$ , dẫn trước và theo sau khối lượng nhỏ $60^{\circ}$ .

$L_1$ , $L_2$ , $L_3$ luôn bất ổn định về động lực (điểm yên ngựa của $\Omega$ ) — một vật đặt ở đó sẽ trôi đi nếu không giữ vị trí. $L_4$ và $L_5$ , dù là điểm yên ngựa của chính $\Omega$ , trở nên thực sự ổn định khi tính cả lực Coriolis, miễn là tỉ lệ khối lượng thỏa tiêu chuẩn Routh $\mu = \dfrac{m_2}{m_1+m_2} < 0.03852$ . Hệ Mặt Trời–Sao Mộc thỏa mãn thoải mái điều này, đó là lý do vì sao hàng nghìn tiểu hành tinh Trojan bị khóa tại các điểm $L_4$ và $L_5$ của Sao Mộc, được phát hiện từ năm 1906 — hơn một thế kỷ sau khi Lagrange tiên đoán sự tồn tại của chúng bằng toán học.

Khác với bài toán hai vật, bài toán ba vật tổng quát không có nghiệm dạng đóng tương tự. Ngoài mười đại lượng bảo toàn cổ điển (năng lượng, ba thành phần mỗi loại của động lượng tuyến tính và mômen động lượng, và chuyển động đều của khối tâm), Henri Poincaré đã chứng minh trong bản ghi nhớ giành giải thưởng năm 1890 cho vua Oscar II của Thụy Điển rằng nói chung không tồn tại thêm tích phân giải tích đơn trị nào của chuyển động — những định luật bảo toàn "ẩn" thêm khiến bài toán hai vật giải được chính xác (như vectơ Laplace–Runge–Lenz ở trên) đơn giản là không còn khi thêm vật thứ ba. Khi đó quỹ đạo có thể phụ thuộc nhạy đến mức dự đoán dài hạn trở nên gần như bất khả thi trên thực tế, dù các phương trình nền tảng hoàn toàn tất định. Khám phá này, sinh ra từ nỗ lực của Poincaré nhằm giải quyết tính ổn định của hệ Mặt Trời, thường được xem là sự ra đời của lý thuyết hỗn loạn.

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

Các kỹ sư nhiệm vụ điều khiển tàu vũ trụ thật bằng cách ghép chính xác các khối hai vật đã chứng minh ở trên. Chuyển giao Hohmann dùng kết quả hình dạng quỹ đạo ở định lý 3 để chọn elip chuyển giao rẻ nhất giữa hai quỹ đạo tròn. Một kính viễn vọng đậu tại điểm Lagrange dùng lý thuyết trên để nằm yên (trong một vòng "halo" nhỏ) gần như không tốn nhiên liệu. Một cú bay ngang trợ lực hấp dẫn khai thác sự thật rằng, nhìn từ hệ quy chiếu của chính hành tinh bay ngang, cuộc chạm trán chỉ là một sự kiện tán xạ hai vật đàn hồi — chỉ có việc đổi hệ quy chiếu mới làm thay đổi năng lượng của tàu.

Ví dụ: Chuyển giao Hohmann tới Sao Hỏa

Quỹ đạo Trái Đất có $a_E=1\text{ AU}$ và Sao Hỏa có $a_M=1.524\text{ AU}$ (cả hai gần tròn và đồng phẳng). Dùng $T[\text{yr}]^{2} = a[\text{AU}]^{3}$ (định luật ba Kepler với khối lượng Mặt Trời gộp vào đơn vị), tìm thời gian chuyển giao một chiều của quỹ đạo chuyển giao Hohmann từ Trái Đất tới Sao Hỏa.

Lời giải

Elip chuyển giao phải tiếp xúc cả hai quỹ đạo tròn, nên bán trục lớn của nó là $a_t = \dfrac{a_E+a_M}{2} = \dfrac{1+1.524}{2} = 1.262\text{ AU}$ (quỹ đạo Trái Đất ở cận điểm của chuyển giao, quỹ đạo Sao Hỏa ở viễn điểm).

Áp dụng $T[\text{yr}]^{2} = a[\text{AU}]^{3}$ cho quỹ đạo chuyển giao: chu kỳ đầy đủ của nó là $T_t = a_t^{3/2}\text{ yr} = 1.262^{1.5}\text{ yr} \approx 1.418\text{ yr}$ .

Chuyến đi một chiều chỉ là nửa elip này, từ cận điểm tới viễn điểm, nên thời gian chuyển giao là $T_t/2 \approx 0.709\text{ yr} \approx 259$ ngày — gần đúng đáng kể với thời gian bay khoảng 7 tới 9 tháng mà các sứ mệnh Sao Hỏa thật sự dùng, vốn bay gần quỹ đạo Hohmann năng lượng tối thiểu này.

Ví dụ: Bay ngang trợ lực hấp dẫn qua Sao Mộc

Một tàu vũ trụ bay ngang qua Sao Mộc (tốc độ quỹ đạo nhật tâm $v_J=13.1\text{ km/s}$ ) trên quỹ đạo phía sau mà, trong giới hạn lý tưởng góc lệch rất lớn, được xấp xỉ tốt bởi một cú "nảy" đàn hồi khỏi hành tinh: $v_{\text{out}} = 2v_J - v_{\text{in}}$ . Nếu tàu đến với tốc độ nhật tâm $v_{\text{in}}=10.0\text{ km/s}$ cùng hướng chuyển động của Sao Mộc, tìm tốc độ rời đi $v_{\text{out}}$ và giải thích năng lượng thêm đến từ đâu.

Lời giải

Trong hệ quy chiếu đứng yên (gần quán tính) của chính Sao Mộc, cuộc chạm trán của tàu với hấp dẫn hành tinh là một sự kiện tán xạ lực xuyên tâm thuần túy: theo bảo toàn năng lượng trong hệ đó, tốc độ của tàu so với Sao Mộc là như nhau rất lâu trước và rất lâu sau cú bay ngang, chỉ hướng của nó thay đổi theo một góc lệch nào đó. Trong hình học lý tưởng phía sau với độ lệch tối đa, vận tốc tương đối đó đảo ngược hướng hoàn toàn.

Chuyển ngược lại hệ Mặt Trời bằng cách cộng vận tốc Sao Mộc $v_J$ biến sự đảo ngược hoàn toàn của vận tốc tương đối $u=v_{\text{in}}-v_J$ thành $v_{\text{out}} = -u+v_J = 2v_J-v_{\text{in}}$ , đúng bằng công thức đã cho. Thay số: $v_{\text{out}} = 2(13.1)-10.0 = 16.2\text{ km/s}$ , tăng $\Delta v = 6.2\text{ km/s}$ .

Năng lượng này không sinh ra từ hư không: trong hệ Mặt Trời, tốc độ quỹ đạo của chính Sao Mộc giảm đi một lượng nhỏ không thể nhận thấy (Sao Mộc nặng gấp khoảng $6\times10^{26}$ lần một tàu vũ trụ, nên bảo toàn động lượng chia đều sự trao đổi nhưng thay đổi vận tốc phía Sao Mộc là hoàn toàn không đáng kể) — tàu vũ trụ thực chất đã "mượn" một mẩu nhỏ động năng quỹ đạo khổng lồ của Sao Mộc. Cú bay ngang thật không bao giờ đạt độ đảo ngược $180^{\circ}$ lý tưởng đầy đủ, nên mức tăng thực tế chỉ là một phần của mức tối đa này, và một cú bay ngang phía trước (đến trước hành tinh thay vì sau) tạo hiệu ứng ngược lại chính xác, làm chậm tàu vũ trụ.

Lỗi hay gặp. Nghĩ rằng trợ lực hấp dẫn vi phạm bảo toàn năng lượng vì tàu tăng tốc 'miễn phí', không đốt nhiên liệu nào. Bảo toàn năng lượng phụ thuộc hệ quy chiếu đối với hệ nhiều vật: động năng của tàu được đo so với Mặt Trời, nhưng bản thân cuộc chạm trán bảo toàn năng lượng trong hệ Sao Mộc, nơi tốc độ tàu gần như không đổi chút nào. Nhảy thẳng tới tốc độ trong hệ Mặt Trời mà không qua hệ hành tinh trước sẽ che giấu năng lượng thực sự đến từ đâu. Luôn phân tích một cú bay ngang qua hai bước: trước tiên trong hệ đứng yên (gần quán tính) của chính hành tinh, nơi cuộc chạm trán là đàn hồi và tốc độ không đổi; sau đó chuyển ngược về hệ Mặt Trời bằng cách cộng vectơ vận tốc quỹ đạo hành tinh. 'Cú tăng tốc' đơn giản là sự phân bổ lại động năng giữa tàu vũ trụ và hành tinh nặng hơn rất nhiều, bị làm chậm một lượng không thể nhận thấy — không có gì được tạo ra từ hư không.

Hiện trạng nghiên cứu tính đến năm 2026

Kết quả không khả tích năm 1890 của Poincaré không kết thúc câu chuyện. Lý thuyết Kolmogorov–Arnold–Moser (KAM), phát triển những năm 1950-60, tinh chỉnh nó một cách ngoạn mục: nó cho thấy với hệ gần khả tích như hệ Mặt Trời (nơi hấp dẫn giữa các hành tinh chỉ là nhiễu loạn nhỏ so với lực hút chủ đạo của Mặt Trời), phần lớn điều kiện đầu — theo nghĩa lý thuyết độ đo chính xác — vẫn cho quỹ đạo tựa tuần hoàn, ổn định hiệu quả, trong khi một tập mỏng, fractal các quỹ đạo cộng hưởng thực sự trở nên hỗn loạn; điều này đồng thời giải thích tính ổn định dài hạn của hệ Mặt Trời và các đặc điểm như khe Kirkwood hỗn loạn trong vành đai tiểu hành tinh. Năm 2000, Alain Chenciner và Richard Montgomery đưa ra chứng minh biến phân chặt chẽ cho vũ điệu ba vật hình "số 8" — ba khối lượng bằng nhau mãi mãi đuổi theo nhau trên một đường cong số 8 cố định với tổng mômen động lượng bằng không — khơi lại hứng thú tìm kiếm thêm các nghiệm tuần hoàn kỳ lạ của bài toán $n$ vật, một cuộc tìm kiếm mang màu sắc tô pô vẫn tiếp diễn tới những năm 2020. Trong khi đó, các hệ ngoại hành tinh nhiều hành tinh quan sát bằng JWST và các thiết bị khác, nổi tiếng nhất là bảy thế giới đá của TRAPPIST-1 bị khóa trong một chuỗi cộng hưởng chuyển động trung bình chặt chẽ, đang được dùng tính đến năm 2026 làm phòng thí nghiệm sống cho đúng câu hỏi của Poincaré: một hệ hành tinh thật có thể nằm gần rìa hỗn loạn tới đâu mà vẫn tồn tại hàng tỉ năm, và điều đó ràng buộc gì về cách các hệ chặt như vậy hình thành và di cư.

Với elip bán trục lớn $a$ và độ lệch tâm $e$ , khối lượng hút đặt tại tiêu điểm, khoảng cách cận điểm (gần nhất) là bao nhiêu?

Một tiểu hành tinh quay quanh Mặt Trời với bán trục lớn $a=4\text{ AU}$ . Dùng $T[\text{yr}]^{2} = a[\text{AU}]^{3}$ , chu kỳ quỹ đạo của nó gần nhất với:

Cặp điểm Lagrange nào có thể ổn định tuyến tính (giữ được vật như tiểu hành tinh Trojan của Sao Mộc) khi tỉ lệ khối lượng thỏa tiêu chuẩn Routh $\mu = \dfrac{m_2}{m_1+m_2} < 0.03852$ ?

Vì sao Kính viễn vọng James Webb được đặt trong quỹ đạo halo quanh điểm $L_2$ Mặt Trời–Trái Đất thay vì quay quanh Trái Đất trực tiếp như Hubble?

Tài liệu tham khảo

Carl D. Murray, Stanley F. Dermott (1999). Solar System Dynamics
Alain Chenciner, Richard Montgomery (2000). A remarkable periodic solution of the three-body problem in the case of equal masses
NASA Science (2024). Webb's Orbit at Sun-Earth Lagrange Point 2 (L2)
NASA Science (2024). Basics of Spaceflight: A Gravity Assist Primer