MathLabs

Hình học

Hình học symplectic

Hình học xây dựng từ một dạng đo diện tích có dấu, bối cảnh tự nhiên của cơ học cổ điển.

Trực giácTrực giác: diện tích có dấu trong không gian pha

Hãy hình dung vị trí và động lượng của một con lắc đang dao động được vẽ cùng nhau thành một điểm (q,p)(q,p) chuyển động trong mặt phẳng. Khi con lắc dao động, điểm này vạch ra một đường khép kín, và diện tích được bao bởi đường đó luôn giữ nguyên dù năng lượng phân bố thế nào giữa vị trí và động lượng. Hình học symplectic nghiên cứu các không gian được trang bị một dạng symplectic ω\omega, một quy tắc đo diện tích có dấu này, và các ánh xạ bảo toàn nó một cách chính xác — ngôn ngữ toán học tự nhiên của cơ học cổ điển.

Hình dựng 3D của một xuyến biểu diễn một tập mức bất biến của hệ Hamilton khả tích trong không gian pha.
Một xuyến bất biến trong không gian pha: với một hệ Hamilton khả tích, định lý Liouville–Arnold phân lá không gian pha thành các xuyến như thế này, mỗi xuyến mang tọa độ tác dụng-góc (I,θ)(I,\theta) mà trên đó dòng chảy là một phép quay theo đường thẳng.

Đại họcĐịnh nghĩa: đa tạp symplectic và trường vector Hamilton

Định nghĩa: Đa tạp symplectic

Một đa tạp symplectic là một cặp (M,ω)(M,\omega) trong đó MM là một đa tạp trơn có số chiều chẵn 2n2n và ω\omega là một dạng vi phân bậc 22 đóng và không suy biến: đóng nghĩa là dω=0d\omega=0, và không suy biến nghĩa là với mọi vector tiếp xúc khác không vv luôn tồn tại một vector tiếp xúc ww sao cho ω(v,w)≠0\omega(v,w)\neq 0.

dω=0d\omega = 0

Tính đóng (dω=0d\omega=0) nghĩa là ω\omega không có nguồn cục bộ nào, và là điều kiện cho phép khôi phục phương trình Hamilton một cách nhất quán từ bất kỳ hàm năng lượng trơn nào. Tính không suy biến nghĩa là ω\omega thiết lập một đẳng cấu tuyến tính giữa các vector tiếp xúc và các đối vector tại mọi điểm, biến vi phân dHdH của một hàm năng lượng HH thành một trường vector thực sự XHX_H thông qua phương trình dưới đây.

ιXHω=dH\iota_{X_H}\omega = dH

Trong tọa độ chuẩn ω=∑i=1ndqi∧dpi\omega=\sum_{i=1}^n dq_i\wedge dp_i, phương trình định nghĩa ιXHω=dH\iota_{X_H}\omega=dH khai triển thành các phương trình Hamilton quen thuộc q˙i=∂H/∂pi\dot q_i=\partial H/\partial p_i, p˙i=−∂H/∂qi\dot p_i=-\partial H/\partial q_i: hình thức luận symplectic trừu tượng và hình thức luận cơ học cổ điển thực chất là cùng một phát biểu, chỉ khác ký hiệu.

Hình học symplectic so với hình học Riemann
Tính chấtSymplectic (M,ω)(M,\omega)Riemann (M,g)(M,g)
Dạng song tuyến tínhPhản đối xứng: ω(v,w)=−ω(w,v)\omega(v,w)=-\omega(w,v)Đối xứng: g(v,w)=g(w,v)g(v,w)=g(w,v)
Mô hình địa phươngLuôn giống nhau: ω=∑i=1ndqi∧dpi\omega=\sum_{i=1}^n dq_i\wedge dp_i (Darboux)Độ cong có thể thay đổi theo từng điểm
Được bảo toàn bởi phép đẳng cấuDiện tích/thể tích ωn\omega^n, không phải độ dài hay gócĐộ dài, góc và khoảng cách trắc địa

Đại họcĐịnh lý then chốt: tính cứng và bảo toàn

Định lý: Định lý Darboux

Mọi đa tạp symplectic (M,ω)(M,\omega) có số chiều 2n2n đều đẳng cấu địa phương với mô hình chuẩn: quanh mỗi điểm luôn tồn tại các tọa độ (q1,…,qn,p1,…,pn)(q_1,\dots,q_n,p_1,\dots,p_n) mà trong đó ω=∑i=1ndqi∧dpi\omega=\sum_{i=1}^n dq_i\wedge dp_i. Nói riêng, đa tạp symplectic không mang bất kỳ bất biến địa phương nào tương tự độ cong.

Vì sao đúng?

Đa tạp Riemann có độ cong, một bất biến địa phương phân biệt mặt cầu với mặt phẳng ngay cả trong một lân cận rất nhỏ. Định lý Darboux nói rằng đa tạp symplectic không có dấu vân tay địa phương như vậy: nhìn gần, mọi đa tạp symplectic đều trông giống hệt không gian pha phẳng R2n\mathbb{R}^{2n} với dạng chuẩn của nó. Vì vậy, mọi nội dung thú vị của tô pô symplectic đều mang tính toàn cục.

Chứng minh

Bước 1 (thiết lập bài toán tuyến tính). Bằng một phép đổi tọa độ tuyến tính, trước hết ta sắp xếp để tại điểm được chọn p0p_0, dạng ω\omega bằng đúng dạng chuẩn ω0\omega_0 tại chính p0p_0, vì mọi dạng song tuyến tính phản đối xứng không suy biến trên R2n\mathbb{R}^{2n} đều có thể đưa về dạng chuẩn nhờ một cơ sở thích hợp.

Bước 2 (nội suy). Định nghĩa họ ωt=(1−t)ω0+tω\omega_t=(1-t)\omega_0+t\omega với t∈[0,1]t\in[0,1]. Vì ω0\omega_0 và ω\omega trùng nhau tại p0p_0 và cả hai đều đóng, không suy biến ở đó, nên ωt\omega_t cũng đóng và, sau khi thu nhỏ lân cận nếu cần, không suy biến với mọi tt.

Bước 3 (giải phương trình Moser). Vì ω\omega trừ ω0\omega_0 đóng và triệt tiêu tại p0p_0, bổ đề Poincaré cho ra một dạng bậc 1 σ\sigma với dσ=ω−ω0d\sigma=\omega-\omega_0. Thủ thuật Moser tìm một trường vector phụ thuộc thời gian XtX_t thỏa ιXtωt=−σ\iota_{X_t}\omega_t=-\sigma; vì ωt\omega_t không suy biến nên điều này xác định XtX_t duy nhất tại mỗi điểm.

Bước 4 (tích phân dòng chảy). Cho φt\varphi_t là dòng chảy sinh bởi XtX_t. Tính toán trực tiếp bằng công thức Cartan cho LXtωt=dιXtωt+ιXtdωt=−dσ=−(ω−ω0)\mathcal{L}_{X_t}\omega_t=d\iota_{X_t}\omega_t+\iota_{X_t}d\omega_t=-d\sigma=-(\omega-\omega_0), kết hợp với ω˙t=ω−ω0\dot\omega_t=\omega-\omega_0, nên ddt(φt∗ωt)=0\frac{d}{dt}(\varphi_t^*\omega_t)=0: phép kéo lùi φt∗ωt\varphi_t^*\omega_t không đổi theo tt và bằng ω0\omega_0 với mọi tt.

Bước 5 (kết luận). Đặt t=1t=1 ta được φ1∗ω=ω0\varphi_1^*\omega=\omega_0 trên lân cận đã thu nhỏ, nên các tọa độ kéo lùi theo φ1\varphi_1 chính là tọa độ Darboux cần tìm.

Cho XHX_H là trường vector Hamilton của HH trên một đa tạp symplectic (M,ω)(M,\omega) có số chiều 2n2n, được xác định bởi ιXHω=dH\iota_{X_H}\omega=dH, và cho ϕt\phi_t là dòng chảy của nó. Khi đó ϕt\phi_t bảo toàn dạng thể tích symplectic: ϕt∗(ωn)=ωn\phi_t^{*}(\omega^n)=\omega^n với mọi tt. Nói riêng, dòng chảy Hamilton bảo toàn thể tích không gian pha.

Vì sao đúng?

Trong cơ học cổ điển, định lý Liouville là lý do toán học khiến một khí gồm các hạt tuân theo phương trình Hamilton không thể tự nén vào một vùng nhỏ hơn của không gian pha: dòng chảy có thể kéo giãn và xoắn vùng đó thành một hình dạng kỳ dị, nhưng thể tích của nó không bao giờ co lại hay tăng lên. Điều này là nền tảng của cơ học thống kê và giải thích tại sao các bộ tích phân số ngây thơ không tôn trọng định luật bảo toàn này có thể tạo ra sự trôi năng lượng giả trong các mô phỏng dài.

Chứng minh

Bước 1 (công thức Cartan). Với mọi trường vector XX và dạng ω\omega, công thức Cartan cho LXω=d(ιXω)+ιX(dω)\mathcal{L}_X\omega=d(\iota_X\omega)+\iota_X(d\omega). Áp dụng với XX bằng XHX_H: vì ω\omega đóng nên dω=0d\omega=0, và vì ιXHω=dH\iota_{X_H}\omega=dH theo định nghĩa trường vector Hamilton, ta được LXHω=d(dH)+0\mathcal{L}_{X_H}\omega=d(dH)+0.

Bước 2 (bản thân dạng không đổi). Vi phân ngoài thỏa d∘d=0d\circ d=0 với mọi hàm, nên d(dH)=0d(dH)=0. Do đó LXHω=0\mathcal{L}_{X_H}\omega=0: dòng chảy Hamilton bảo toàn chính dạng symplectic, không chỉ thể tích của nó.

Bước 3 (chuyển sang lũy thừa cao nhất). Dạng thể tích là ωn=ω∧⋯∧ω\omega^n=\omega\wedge\cdots\wedge\omega (gồm nn thừa số). Quy tắc Leibniz cho đạo hàm Lie trên tích ngoài cho LXH(ωn)=n ωn−1∧LXHω\mathcal{L}_{X_H}(\omega^n)=n\,\omega^{n-1}\wedge\mathcal{L}_{X_H}\omega, và vì LXHω=0\mathcal{L}_{X_H}\omega=0 theo Bước 2, vế phải triệt tiêu: LXH(ωn)=0\mathcal{L}_{X_H}(\omega^n)=0.

Bước 4 (tích phân theo dòng chảy). Nếu ϕt\phi_t là dòng chảy của XHX_H, thì LXH(ωn)=0\mathcal{L}_{X_H}(\omega^n)=0 chính xác có nghĩa là ddtϕt∗(ωn)=0\frac{d}{dt}\phi_t^{*}(\omega^n)=0 với mọi tt. Vì ϕ0=id\phi_0=\mathrm{id} cho ϕ0∗(ωn)=ωn\phi_0^{*}(\omega^n)=\omega^n, tích phân cho thấy ϕt∗(ωn)=ωn\phi_t^{*}(\omega^n)=\omega^n với mọi tt, đó chính là định lý Liouville.

Nâng caoNâng cao: tính cứng symplectic và định lý không nén được của Gromov

Gromov chứng minh vào năm 1985 rằng một quả cầu symplectic B2n(r)={∣z∣≤r}B^{2n}(r)=\{|z|\le r\} có thể được nhúng symplectic vào hình trụ B2(R)×R2n−2B^2(R)\times\mathbb{R}^{2n-2} khi và chỉ khi r≤Rr\le R — chính xác cùng điều kiện như khi nhúng một đĩa 22 chiều bán kính rr vào một đĩa bán kính RR, dù nn lớn đến đâu. Chứng minh sử dụng các đường cong giả chỉnh hình và khai sinh ra các bất biến gọi là dung lượng symplectic, trong đó đơn giản nhất là bề rộng Gromov c(B2n(r))=πr2c(B^{2n}(r))=\pi r^2, diện tích của đĩa nhỏ nhất mà một quả cầu có thể bị nén qua.

c(B2n(r))=πr2c(B^{2n}(r)) = \pi r^2

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

Vì cơ học Hamilton được viết bằng ngôn ngữ symplectic, hình học symplectic xuất hiện ở bất cứ đâu một hệ vật lý bảo toàn năng lượng và cấu trúc không gian pha: cơ học thiên thể dùng nó để nghiên cứu độ ổn định lâu dài của quỹ đạo hành tinh, thiết kế máy gia tốc hạt dùng các ánh xạ symplectic để theo dõi chùm hạt qua hàng triệu vòng quay mà không bị tắt dần hay phình to giả tạo, và phần mềm động lực học phân tử dùng các bộ tích phân symplectic để các phân tử mô phỏng không được hay mất năng lượng chỉ vì sai số số học. Lý thuyết điều khiển tối ưu và lượng tử hóa hình học trong vật lý toán cũng xây dựng trực tiếp trên hình thức luận symplectic.

Ví dụ: Kiểm tra một phép quay là ánh xạ symplectic

Trên (R2,ω=dq∧dp)(\mathbb{R}^2,\omega=dq\wedge dp), xét phép quay ϕθ(q,p)=(qcos⁡θ−psin⁡θ, qsin⁡θ+pcos⁡θ)\phi_\theta(q,p)=(q\cos\theta-p\sin\theta,\ q\sin\theta+p\cos\theta) theo một góc cố định θ\theta. Liệu ϕθ\phi_\theta có phải là một phép symplectomorphism, nghĩa là có thỏa ϕθ∗ω=ω\phi_\theta^{*}\omega=\omega hay không?

Lời giải

Bước 1: viết ma trận Jacobi. Ánh xạ là tuyến tính với ma trận J=(cos⁡θ−sin⁡θsin⁡θcos⁡θ)J=\begin{pmatrix}\cos\theta&-\sin\theta\\ \sin\theta&\cos\theta\end{pmatrix} tác động lên (q,p)(q,p).

Bước 2: kéo lùi dạng. Với một ánh xạ tuyến tính trên R2\mathbb{R}^2, ta có ϕθ∗(dq∧dp)=det⁡(J) dq∧dp\phi_\theta^{*}(dq\wedge dp)=\det(J)\,dq\wedge dp, nên chỉ cần tính định thức của JJ.

Bước 3: tính định thức. det⁡(J)=cos⁡2θ+sin⁡2θ=1\det(J)=\cos^2\theta+\sin^2\theta=1 với mọi θ\theta, theo hệ thức Pythagore.

Bước 4: kết luận. Vì det⁡(J)=1\det(J)=1, nên ϕθ∗ω=ω\phi_\theta^{*}\omega=\omega với mọi θ\theta: các phép quay là những phép symplectomorphism, đúng như sự kiện hình học rằng phép quay bảo toàn diện tích.

Ví dụ: Euler symplectic so với Euler tường minh cho dao động điều hòa

Với dao động điều hòa H(q,p)=12(p2+q2)H(q,p)=\tfrac12(p^2+q^2), hãy so sánh một bước của phương pháp Euler symplectic qn+1=qn+hpn, pn+1=pn−hqn+1q_{n+1}=q_n+hp_n,\ p_{n+1}=p_n-hq_{n+1} với một bước của phương pháp Euler tường minh (ngây thơ) qn+1=qn+hpn, pn+1=pn−hqnq_{n+1}=q_n+hp_n,\ p_{n+1}=p_n-hq_n, với bước nhỏ h>0h>0. Phương pháp nào phù hợp với định lý Liouville?

Lời giải

Bước 1: viết ánh xạ Euler symplectic như một hàm của (qn,pn)(q_n,p_n). Thay qn+1=qn+hpnq_{n+1}=q_n+hp_n vào phương trình thứ hai ta được pn+1=pn−h(qn+hpn)=pn−hqn−h2pnp_{n+1}=p_n-h(q_n+hp_n)=p_n-hq_n-h^2p_n.

Bước 2: tính ma trận Jacobi của nó. ∂(qn+1,pn+1)∂(qn,pn)=(1h−h1−h2)\frac{\partial(q_{n+1},p_{n+1})}{\partial(q_n,p_n)}=\begin{pmatrix}1&h\\ -h&1-h^2\end{pmatrix}, với định thức 1⋅(1−h2)−h⋅(−h)=1−h2+h2=11\cdot(1-h^2)-h\cdot(-h)=1-h^2+h^2=1 với mọi hh.

Bước 3: tính ma trận Jacobi của Euler tường minh. Ở đó ∂(qn+1,pn+1)∂(qn,pn)=(1h−h1)\frac{\partial(q_{n+1},p_{n+1})}{\partial(q_n,p_n)}=\begin{pmatrix}1&h\\ -h&1\end{pmatrix}, với định thức 1+h21+h^2, luôn lớn hơn 11 với mọi h≠0h\neq0.

Bước 4: diễn giải. Theo định lý Liouville, dòng chảy chính xác có định thức Jacobi đúng bằng 11 ở mỗi bước. Euler symplectic khớp chính xác điều này bất kể hh, nên nó không làm phình thể tích không gian pha một cách giả tạo; Euler tường minh mở rộng diện tích theo hệ số 1+h21+h^2 ở mỗi bước, khiến năng lượng trôi dần lên sau nhiều bước — đây là lý do các mô phỏng dài hạn dùng các bộ tích phân symplectic.

Tại sao một đa tạp symplectic phải có số chiều chẵn 2n2n?

Với dao động điều hòa H(q,p)=12(p2+q2)H(q,p)=\tfrac12(p^2+q^2) và dạng chuẩn ω=dq∧dp\omega=dq\wedge dp, trường vector Hamilton XH=(q˙,p˙)X_H=(\dot q,\dot p) được xác định bởi ιXHω=dH\iota_{X_H}\omega=dH là gì?

Theo định lý Liouville, dòng chảy Hamilton trên một đa tạp symplectic bảo toàn điều gì?

Các mô phỏng dài hạn hệ mặt trời kéo dài hàng triệu năm dùng bộ tích phân symplectic thay vì các phương pháp đa dụng như Runge-Kutta chuẩn. Tại sao?

Tài liệu tham khảo

  1. Dusa McDuff, Dietmar Salamon (2017). Introduction to Symplectic Topology
  2. Mikhail Gromov (1985). Pseudo holomorphic curves in symplectic manifolds
  3. Pazit Haim-Kislev, Yaron Ostrover (2026). A Counterexample to Viterbo's Conjecture · arXiv:2405.16513