肖志峰,林文斌

(南華大學 數(shù)理學院,湖南 衡陽 421001)

0 引 言

廣義相對論最重要的應(yīng)用之一是理解時空曲率對天體運動的影響[1]。廣義相對論描述了牛頓力學框架下所不能描述的宇宙大范圍的行為，開啟了天體物理全新的研究領(lǐng)域。在愛因斯坦的廣義相對論中[2]，時空是一個四維的偽黎曼流形，物質(zhì)與時空的相互作用是非線性的，因此一般而言，廣義相對論的計算是非常復(fù)雜的，這也使得對廣義相對論的理解更加困難。

基于上述原因，愛因斯坦等人在二十世紀三十年代發(fā)展了一種攝動方法，能夠描述廣義相對論對太陽系中的天體演化所引起的細微變化，A.Einstein、L.Infeld、B.Hoffmann[3]在1938年得到的結(jié)果就是一階后牛頓(1PN)的結(jié)果，這也標志著“相對論天體力學”這一新領(lǐng)域的誕生。此后，相對論天體力學的研究一直方興未艾，對后牛頓展開進行了大量不同的分析，同時提出了許多新的更有力的工具來研究相對論天體力學[4-5]。盡管相對論天體力學取得了長足的進步，但得到后牛頓展開的解析仍然是一項艱巨的任務(wù)。特別是，到目前為止，軌道運動和旋轉(zhuǎn)運動之間的相互作用還沒有得到足夠的研究。

本文著重研究了帶自旋雙星的運動問題的二階后牛頓(2PN)展開，目標是基于李級數(shù)攝動理論[6-7]得到一個運動粒子的軌道角動量的z分量關(guān)于時間t演化的解析解，這個方法可以得到系統(tǒng)的長期演化行為，并且可以拓展到更高階的后牛頓近似展開。

1 哈密頓量

本文的推導(dǎo)從自旋雙星問題的2PN哈密頓量開始。在折合質(zhì)量與折合質(zhì)量坐標系中，自旋雙星問題的2PN哈密頓量為[8]

H=HN+εH1+ε2H2

(1)

這里ε=1/c2是選擇的一個小參數(shù)，其中c表示光速。在攝動理論中

(2)

(3)

H1的表達式為

H1=HPN+HSO+HSS,

(4)

這里

(5)

(6)

是軌道-自旋相互作用項，且

(7)

是自旋耦合作用項。

H2的表達式為

(8)

(9)

(10)

(11)

(12)

(13)

(14)

2 作用角變量

為了使用李級數(shù)攝動法來研究哈密頓量式(9)～式(13)，我們需要將哈密頓量表示在作用角變量坐標中的攝動形式。對于軌道運動變量，一般選擇Delaunay元素集來表示。而對于旋轉(zhuǎn)運動變量，一般選擇Serret-Aadoyer(SA)變量。這兩組變量都是正則變量。

對于Delaunay元素集，可以通過以下標準關(guān)系來引入[10]

H=Gcosih=Ω

(15)

這里a是軌道半長軸，e是軌道偏心率，i的軌道傾角，ω是近心點幅角，M是平近點角，Ω是升交點黃經(jīng)。L是半主軸的歸一化平方根，G是軌道角動量的歸一化，H是軌道角動量的歸一化的z分量。在方程中，(L,G,H)表示廣義動量，(l,g,h)表示廣義坐標。

SA變量[11]是描述一個剛體關(guān)于旋轉(zhuǎn)角和旋轉(zhuǎn)角動量的相互關(guān)系。在討論中，物體是一個理想的球體，所以SA變量的形式可以表示為

(16)

(17)

(18)

(19)

(20)

r=r(L,G,l),

(21)

Gxy=Gxy(G,H),

(22)

(23)

f=f(L,G,l)。

(24)

現(xiàn)在可以使用李級數(shù)[6-7]來分析哈密頓量。李級數(shù)擾動方法的目的是通過坐標的擬正則變換來簡化問題的哈密頓量，這個變換取決于生成函數(shù)χ。李級數(shù)變換如下

(25)

(26)

因此需要解同態(tài)方程

{HN,χ}+H1=K。

(27)

(28)

(29)

在方程中，其本質(zhì)表為平近點角的運動，因此，通過李級數(shù)變換且在方程中確定了χ后，哈密頓量變?yōu)槿缦滦问?/p>

(30)

(31)

(32)

通過式(31)、式(32)的變換后，最終形式的哈密頓量為

(33)

這里F0、F1、F2僅僅是動量的函數(shù)

(34)

(35)

(36)

(37)

現(xiàn)在可以來分析最終形式的哈密頓量了。首先將哈密頓量代入正則方程得到的運動方程為

(38)

和

(39)

將方程和方程中關(guān)于H的方程

(40)

結(jié)合起來，可以得到

(41)

這里H′是哈密頓常數(shù)。通過計算，可以證明，方程(41)中右邊關(guān)于H5和H6的項都抵消了。因此方程(41)的右邊是一個關(guān)于H的四次多項式[12]，即

f4(H)=ε2F12-(H′-HN-εF0-ε2F2)2=a0H4+a1H3+a2H2+a3H+a4,

(42)

其中a0,a1,a2,a3,a4是代入各項后整理得到的系數(shù)。重寫方程為

(43)

這里H0是H的初始值，t0是初始時間。除符號外，方程的左邊是標準的橢圓積分。通過Weierstrass公式[12]，可以寫為

(44)

這里a是一個任意常數(shù)，于是有

(45)

這里ξ(z)≡ξ(z;g2,g3)定義為關(guān)于不變量g2、g3的Weierstrass橢圓函數(shù)

(46)

其中f4的導(dǎo)數(shù)表示關(guān)于變量H的導(dǎo)數(shù)。H關(guān)于時間t的演化方程的解為

(47)

3 結(jié) 論

在本文中，文章基于李級數(shù)攝動理論分析了帶自旋雙星問題的演化。基于李級數(shù)攝動理論，得到了橢圓函數(shù)形式的解析解，從而可以分析系統(tǒng)的長期擾動效應(yīng)，這是用數(shù)值方法所不能做到的。這個方法也可以推廣到更高階的后牛頓展開，以得到更精確的結(jié)果。