李林森

(東北師范大學(xué)物理學(xué)院，長春 130024)

自從狄拉克(Dirac)(1937)[1-2]提出了引力常數(shù)隨時(shí)間變化后，這種理論在天文學(xué)和地球物理學(xué)領(lǐng)域顯示很大影響。首先，Kholchenikov和Fracassini(1968)[3]根據(jù)狄拉克的大數(shù)假說研究了具有引力常數(shù)變化的二體問題。他們將理論應(yīng)用于行星經(jīng)度變化的估計(jì)。Hut和Verhulst(1976)[4]研究了引力常數(shù)隨時(shí)間減少的二體問題，并推出了平均運(yùn)動和平均經(jīng)度的長期變化的方程式，給出了真經(jīng)度的性質(zhì)[4]。同年，Bishop和Landsberg(1976)[5]從能量守恒定律推出了具有引力常數(shù)隨時(shí)間變化的牛頓引力定律的運(yùn)動方程公式。隨后Maite(1978)[6]利用此理論研究在圓形軌道中引力常數(shù)隨時(shí)間變化和地球表面的有效溫度之間的關(guān)系，但他沒有研究引力常數(shù)變化和圓形軌道演變之間的關(guān)系。本文在這方面做了研究。一般來說，食雙星軌道并非都為圓形軌道或軌道偏心率為零，但有些密近食雙星的軌道為圓形軌道，所以本文的理論結(jié)果適用于密近食雙星軌道的演變。

1 引力常數(shù)隨時(shí)間變化產(chǎn)生的軌道根數(shù)變化(圓形軌道情形)

Maiti(1978)根據(jù)Bishop和Landsberg(1976)給出的公式推出了具有引力常數(shù)隨時(shí)間變化的矢量運(yùn)動方程[5]：

(1)

(2)

其中r為圓軌道半徑；n為平均運(yùn)動(軌道平均角速度)。 將(1)式寫成標(biāo)量方程：

(3)

(4)

將G(t)展開為Taylor級數(shù)：

(5)

其中G(t0)=G0為牛頓引力常數(shù)(不隨時(shí)間變化)。

(6)

(7)

將方程(6)～(7)同牛頓運(yùn)動方程相比較后可得后牛頓項(xiàng)：

(8)

(9)

因此后牛頓攝動加速度的徑向和切向分量是：

(10)

(11)

(δFr)PN和(δFθ)PN可以認(rèn)為具有引力常數(shù)G0的牛頓運(yùn)動方程的改正項(xiàng)或?yàn)樽円Τ?shù)G(t)的攝動項(xiàng)。因此，S和T就是具有G0的高斯攝動方程的攝動加速度。 利用以時(shí)間為自變量的高斯方程[7]：

其中a為軌道半長軸；e為軌道偏心率；E為偏近點(diǎn)角；p=a(1-e2)。

本文研究圓形軌道情形。對于圓形軌道：a=r(圓軌道半徑)；e=0；E=θ=nt；p=r。將這些代入軌道半長徑和軌道偏心率隨時(shí)間變化的高斯攝動方程[7]，得：

(12)

(13)

將(10)～(11)的S和T代入(12)～(13)式，得：

(14)

(15)

(16)

再推出軌道周期P的變化。將n=2π/P代入Kepler第三定律并對時(shí)間微分，得：

將(14)代入上方程，積分上述方程，可得：

(17)

下面推出軌道偏心率的變化。積分方程(15)后，可得到軌道偏心率e的變化：

(18)

結(jié)果表示：軌道半徑和周期以線性形式隨時(shí)間而縮減，圓形軌道的形狀以周期和混合長期項(xiàng)隨時(shí)間而變化。

2 對5個(gè)密近食雙星軌道半徑和軌道周期的長期演變的上下限估計(jì)

表1 5個(gè)密近食雙星軌道半徑和軌道周期長期演變的上下限估計(jì)

3 討論和結(jié)論

(1) 由(16)式和(17)式可知，軌道半徑和軌道周期以線性形式隨時(shí)間縮小。

(4) 因二體問題在共面上運(yùn)動，垂直軌道平面的攝動加速度分量W=0，將其代入攝動方程后軌道傾角的i和升變點(diǎn)經(jīng)度Ω不變。

(5) 同其他4種效應(yīng)相比較。如果將e=0(圓軌道情形)代入作者在文[10-13]中橢圓軌道情形中的雙星多方模型攝動效應(yīng)、后牛頓效應(yīng)、二體自轉(zhuǎn)效應(yīng)和引力輻射阻尼效應(yīng)的式子，可得引力常數(shù)變化對圓形軌道半徑和軌道偏心率的效應(yīng)同另外4種效應(yīng)的對比如表2。

表2 在圓形軌道情形中引力常數(shù)變化對軌道半徑和軌道偏心率產(chǎn)生的效應(yīng)和其他4種效應(yīng)的對比

由表2可以看出，只有引力輻射阻尼對軌道半徑產(chǎn)生的長期效應(yīng)同引力常數(shù)變化產(chǎn)生的長期效應(yīng)相一致。而對軌道偏心率除前者有混合長期項(xiàng)外兩者皆有周期項(xiàng)。兩者在長期效應(yīng)可在量級上比較。

[1] Dirac P A M.The cosmological constant[J].Nature,1937,139:323.

[2] Dirac P A M.A New basis for cosmology[J].Proc R Soc London A,1938,A165:199-208.

[3] Kholchevnikov C,Fracassins M.La Probleme dex edeux corps avec G variable selon l’hypothese de Dirac[J].Conferenze Dell’Osservatorio Astronomico di Milano Serie 1，1968, 9:1-50.

[4] Hut P,Verhulst F.The two-body problem with a decreasing gravitational constant[J]. M N R A S,1976,177:545-549.

[5] Bishop N T,Landsberg P T.Time-Varing Newtonian gravity and universal motion[J]. Nature,1976,264(25):346-347.

[6] Maiti S R.Time-varying gravitation in Newtonian theory[J].MNRAS,1978,185:293-295.

[7] Brouwer D,Clemence C M.Method of Celestial Mechanics[M].Academic press,New York and London,1961.

[8] http://binaries.boulder.swri.edu/atlas/.

[9] AI-Rawaf A S.Limits on the changes in the gravitational constant from nucleosynthesis [J].Astrophys &Space Sci,2007,310(1-2):173-175.

[10] 李林森.雙星多方模型的形狀及其對同步子星軌道要素的攝動影響[J].云南天文臺刊，1997,(4):9-167.

Li Linsen.Configuration of the Polytropic Model of Binary Star and Its Perturbation Influence on the Orbital Elements of the Synchronous Components in Close Binary Sysetm[J].Publications of The Yunnan Observatory,1997,(4):9-17.

[11] 李林森.天體軌道要素變化的后牛頓效應(yīng)[J].中國科學(xué)，A輯(數(shù)學(xué)物理學(xué)天文學(xué)技術(shù)科學(xué))1988，31(5):523-529.

[12] 李林森.雙星系中兩子星的自轉(zhuǎn)對雙星軌道變化的后牛頓效應(yīng)[J].天文學(xué)報(bào)，2001，42(4):428-435.

Li Linsen.The Post-newtonian Effects Due to Rotation of Stars on the Variation of the Orbital Elements of the Component in Binary System[J].Acta Astronomica Sinica,2001，42(4):428-435.

[13] 李林森.引力輻射阻尼對雙星軌道要素變化的影響[J].物理學(xué)報(bào)，1989，38(11):1877-1881.

Li Linsen.The effect of gravitational radiation damping on the Variation of the orbital elements of binary systems[J].Acta Phrysica Sinica,1989,38(11):1877-1881.

[14] Khajavi M,Edalate M T,Jassur D M Z.BVRI photometry and light curve analysis of V W Cep[J].Astrophys & Space Sci,2002,282(4):645-653.