過去・未来の星空 計算の方法

「過去未来の星空」コンテンツでは、かなり長い期間の恒星の運動を調べています。
どうやって計算しているのだろう、と興味のある方もいらっしゃるかと思いますので、ここにまとめてみました。

ここでの計算の大前提として、恒星は等速直線運動をしていることとします。実際には、星はすべて銀河系を公転していますから、直線ではなく楕円運動をしているでしょうし、恒星間の引力の働きもあると思います。ただ、現在から前後100万年という期間では、等速直線運動とみなしてよいのではないかと考えます。
銀河系の重力の影響を含めて計算している天文学者の方たちもいらっしゃいます。その結果と比べても大きく異なる結果は出ていないので、ここでの計算方法も、近似値として十分に使えると思います。

星の固有運動

遠い過去・未来の星空を調べるためには、星ひとつひとつがどのように動いているか、を、知る必要があります。これを「固有運動」と呼び、ひとつひとつの星の固有運動を計算できれば、過去・未来の星空を知ることができます。星の固有運動を計算する、といっても、特に難しいことはなく「時速1000キロメートルで南に向って飛んでいる飛行機は、2時間後には何キロメートル先にありますか?」という問題と本質は同じ。

算数の問題と違うところは、最初の「時速1000キロメートルで南に」という部分も観測によって調べなければならないところですが、現在では、過去未来の星空を計算するためにそこから始める必要はありません。天文台で観測した結果が、星のカタログ=星表にまとめられています。
すべての星表に掲載されているわけではありませんが、位置や明るさ、スペクトルなどのほかに、固有運動に関わるデータを掲載しているものがあります。

星表に掲載されている、過去未来の星の位置の計算に必要な情報は、

  • 星の位置         alpha  delta      赤径、赤緯
  • 1年間の固有運動量  mu alpha  mu delta   地球から見た星の動き
  • 視線速度         R0        近づいてきているか、離れていっているかの速度
  • 年周視差         pi         その星までの距離

この4つ。

図1 固有運動量

図1 固有運動量

これらを使えば、過去未来の星空を知ることができます。では、その方法を紹介していきましょう。

固有運動の簡単な求め方

上の4つの要素を見ると、そのものずばり「固有運動」量というものがあります。これは、まさに星の動きの量、地球から見て1年間に動いた量を示しています。これを使えば、星の動きを求めることができます。
星表に出ている固有運動量は一般的に、赤径(mu alpha)=横に動いた量、赤緯(mu delta)=縦に動いた量、の2つの要素に分けて掲載されています。

星の位置も、地球上の位置と同じで座標で表し、地上で言う経度が赤径、緯度が赤緯にあたります。現在の星の位置は主に、赤径・赤緯の二つの数値で示していますから、そのそれぞれにmu alpha(赤径分の移動量)、mu delta(赤緯分の移動量)を、必要な年数分、足してあげれば、その星の過去未来の位置が計算できます。

 過去未来の赤径 = 現在の赤径alpha + mu alpha(赤径) × 年数
 過去未来の赤緯 = 現在の赤緯delta + mu delta(赤緯) × 年数
 年数は、未来はプラス、過去ならマイナスです

この式でも、そう遠くない過去未来なら、星の位置を知ることができます。掛け算と足し算だけですから、本当に簡単に求めることができます。
けれども、何万年、何十万年という話になると、これではよろしくありません。

というのは、上の計算で使ったmu alpha(赤径)、mu delta(赤緯)は、地球から見た、現在の動く方向、そして速さだからです。星が空に張り付いている、という前提の式で、あくまでも近似式です。星は、私たちと同じ3次元を運動していますから、天球に映った星の動きを使ったのでは、正しい計算はできません。

正面から飛行機がこちらに向かってくるところを想像してみてください。遠くに見えているときは、位置はあまり変わらずにゆっくり動いているように見えます。近づくにつれてどんどん速くなり、頭上を通過するときには、あっという間に後へと位置が変わり、遠ざかるとともにまた、位置は変わらずゆっくり動いているように見えるでしょう。まっすぐ飛んでいれば、進む向きが変わって見えることはありませんが、同じ速さで飛んでいても、地上からでは速さはずいぶん変わって見えます。この飛行機の動きを「遠くにあるころの動き」「頭上を通過するときの動き」など、その瞬間に見えた速さだけを基にしては、正しく知ることができない、ということです。見えた動きから、飛行機の実際の速さを調べて、それを基にして見える位置を計算しなければ、間違った答えになってしまいます。星も、おなじこと。星の動きを長期間調べるには、星の実際の速さを知る必要があるのです。

図1 固有運動量(再掲)

図1 固有運動量(再掲)

全固有運動量と、固有運動の向き

もう一度、図1を見てください。図にあるpsi(プシー、プサイ)とmu(ミュー)は、地球から見た固有運動の向きと速さです。これからの計算のために、この2つを求めます。
固有運動に関する情報はmu alphamu delta。観測で知ることができた、赤径赤緯それぞれの固有運動量を「速さ」と「動いていく方向」に分けなおします。

まずは、固有運動の向きpsiを求めてみましょう。固有運動の向きpsiは、天の北極、北を0度として、東→南→西と360度まで測ります。直角三角形なので、三角比を使って、

 psi=tan( {mu alpha} / {mu delta} )     ( {mu delta} >= 0″ title=”{mu delta} >= 0″/>  で、I、IV象限、  <img src=  でII、III象限)

として求めます。
「星が等速直線運動をしている」という前提で計算していますから、この値はこのまま最後まで使います。

次に、1年間の地球から見た固有運動の速さ、全固有運動量muは、mu alphamu deltaをあわせたもの。

 mu=sqrt{mu alpha^2+mu delta^2}

として求めます。

図2 恒星運動の角距離

図2 恒星運動の角距離

宇宙空間での星の動き

つづいて図2。上で示したデータを使って星の動きの計算方法を調べていきましょう。

始め、恒星はAにいました。そして、時間が過ぎてt年後、Bまでやってきました、という図です。これは、地球から眺めたのではなく、宇宙空間から見ています。AからBに動いた、という動きは、宇宙空間でのもの。

宇宙空間でt年間かけて星はAからBへと動きましたが、この動きを地球から眺めていると、S(太陽)と書かれている点の角sigma(シグマ)分、動いたように見えるはずです。これを角距離sigmaと呼びます。
これがt年間ではなく、1年間ということであれば、先ほど求めた全固有運動量muと同じです。最初に紹介した簡単な固有運動の計算方法は、sigmaを求めるのに、このmuを単純にt倍していたことになりますが、それが問題だったわけです。正確にするには、宇宙空間での実際の星の運動を考えなければなりません。

この角距離sigmaが分かれば、恒星が進む方向はpsiで分かっていますから、t年後には恒星がどこにいるか、が分かります。
ただ、このsigmaは、最初に提示した4つのデータの中にはありませんから、ちょっと工夫をしないといけません。

図を見ると、ABを結んでいるVというのが、実際の宇宙空間でのt年間の星の動き。これを知りたいのですが、ここでVの代わりに、奥行きRと横向きv、ふたつの要素に分けて考えてみます(図では、t年間の動きなので、Rtvtとなっています) というのは、奥行き方向の速度R、これは視線速度と呼ばれ、観測で求められているからです。最初に提示した4つの情報のひとつ、星表に[km/s]単位で掲載されています。このデータから、視線方向の星の動きは計算できます。

もうひとつの横向き方向の速度v、これは接線速度と呼ばれ、こちらは残念ながらデータにありませんので、計算で求めないといけません。この計算には、全固有運動量muともうひとつ、年周視差piを使います。
ちょっと頭の体操を。

視線速度、接線速度

ではさっそく、接線速度の計算式を考えましょう。図3の三角形を使います。図2と同じですが、補助線の引き方が違っています。また、ABt年間ではなく1年間の動き、とします。t年間の動きの計算のために、1年間の量を求めるのです。

図3 視線速度 接線速度

図3 視線速度 接線速度

三角形ACSは直角三角形です。なので、三角比の定義から、接線方向の移動量ACは、地球と恒星の距離dに、sin muを掛けたもの。

AC = d sin mu     ①

muは角度ですから、星までの距離が決まらないと、移動量も距離の単位では求まらないということです。同じ角度の1度の動きでも、近ければゆっくりでしょうし、うんと遠ければうんと早いでしょう。
ACを求めるには、星までの距離dを決めてあげないといけませんが、ここで、

   sin pi = a / d     ②

という式を提示します。これは、年周視差piを求める式。aは 地球軌道半径、dは恒星までの距離。地球軌道半径aを、恒星までの距離dで割ったものが、地球の公転運動によって動いて見える恒星の角度 = 年周視差、ということです。
年周視差piから星までの距離dを求めるには、②式を変形して

  d = a/ {sin pi}       ③

①に代入すると、

  AC = a*({sin mu} / {sin pi})     ④

これで、1年間の接線方向の移動量ACを求めることができます。
これを[km/s]単位の接線速度にするには、地球軌道半径a[km]を、1年間の秒数nで割って、

  接線速度v =  (a/n)*({sin mu} / {sin pi})[km/s]

     地球軌道半径a = 149600000[km]
     1年間の秒数n = 86400*365.2422 = 31560000[s] より、

     (a/n) = 149600000 / 31560000 = 4.741[km/s]
     なので、
  接線速度v = (a/n)*({sin mu} / {sin pi}) = 4.741*({sin mu} / {sin pi}) [km/s]

と、求まりました。

でも、ちょっと待ってください。今、そもそも何を計算したいのか、というと、何千年、何万年という過去未来の星の位置です。それを[km/s]といった単位で計算しても、どのみち[光年/年]や、[パーセク/年]に変換する必要がでてくるでしょう。接線速度の単位は、今の(a/n)の単位で決まりますから、時間の単位は[年] 距離の単位は[光年]で計算することにしましょう。

式④の場合、時間は始めから1年間として計算していますからnは1[年]、地球軌道半径aを、1光年を[km]単位で表した9460000000000[km]で割って、

  a = {149600000 km}  / {9460000000000 km} = 0.000001581
  v[光年/年]=0.000001581 * {{sin mu} / { sin pi}}   ⑤

とします。

あわせて視線速度も[光年/年]に変換しておきましょう。元々が[km/s]ですから、1年間の秒数を[km]単位の1光年で割ったものを掛けてあげればよく、

  31560000/9460000000000 = 0.000003336
  R[光年/年]=0.000003336 * R

これで、接線速度、視線速度ともに[光年/年]単位で求めることができました。1

sigmaを求める

図2 恒星運動の角距離(再掲)

図2 恒星運動の角距離(再掲)

上で、[光年/年]単位の視線速度R、接線速度vが求まりました。これらを使ってt年間の角距離sigmaを求めます。
再び図2です。

 直角三角形BA prime Sで、三角比の定義より、

 d prime sin sigma = A prime B
 d prime cos sigma = (d + A A prime)

 A prime B は、t 年間の視線方向の移動距離 = v * t
 A A prime は、t 年間の視線方向の移動距離 = R * t

なので、

 d prime sin sigma = v * t         ⑥-1
 d prime cos sigma = (d + R * t)   ⑥-2

となります。このふたつの式から、t年間の角距離sigmaと、t年後の太陽系からの距離d primeを求めることができます。

tan sigma = {sin sigma} / {cos sigma}より、⑥-1式の値を⑥-2の値で割り、tan sigmaの値からsigmaを求めます。⑥-1、⑥-2式の符号で、どの象限かがわかるわけです。
求めたsigmaで、たとえばcos sigmaを計算して、⑥-2の値を割れば、d primeも求まります。単位は[光年]です。

ここまで長かったですが、t年間の恒星の移動角距離sigmaがようやく求まりました。sigmaと固有運動の向きpsiを使って、t年後の恒星の位置を計算しましょう。

固有運動で星を動かす

ここからは、球面三角法という、球面上の「正弦法則」「余弦法則」などを使って、天球上で星がどう動くのかを計算します。が、球面三角法は複雑で、説明が長くなりますので、答えだけをここに出しましょう。

t年後の星の位置を求めるために使うデータは、現在の星の位置alphadelta、固有運動方向psi、そして今求めたt年間の移動角距離sigmaです。t年後の星の位置は、alpha prime delta primeで表してあります。
⑦-1式や⑦-2式の(alpha prime - alpha)は、これでひとつの未知数、現在の赤径との差を表し、alpha prime - alphaという意味ではありません。ばらして別々に計算したりしないようにしてください。

 cos delta prime sin( alpha prime - alpha) =             sin sigma sin psi      ⑦-1
 cos delta prime cos( alpha prime - alpha) = cos delta cos sigma - sin delta sin sigma cos psi     ⑦-2
 sin delta prime           = sin delta cos sigma + cos delta sin sigma cos psi   ⑦-3

まず、⑦-1で求めた数値を⑦-2で求めた数値で割ってtan(alpha prime - alpha)にし、象限判断をして赤径の移動量(alpha prime - alpha)を求めます。
現在位置からの移動量なので、現在位置の赤径alphaに、求めた(alpha prime - alpha)を足して、t年後の赤径が計算できます。
t年後の赤緯delta primeは、90から-90の範囲なので、象限判断は不要で、⑦-3式からそのまま計算できます。

まとめ

これで、t年後の星の位置が計算できるようになりました。天文計算の中では簡単なものですが、恒星の数が多くなると、計算は面倒です。
このサイトで使っている星図は、山上企画にて開発した、上の計算と、その結果を星図として描画するプログラムを使用して作成しています。

最後に、この計算に必要な式をまとめておきましょう。

恒星の固有運動計算に必要なデータ:

  • 星の位置         alpha  delta      赤径、赤緯
  • 1年間の固有運動量  mu alpha  mu delta   地球から見た星の動き
  • 視線速度         R0        近づいてきているか、離れていっているかの速度
  • 年周視差         pi         その星までの距離

固有運動方向psi 全固有運動量mu

mu=sqrt{mu alpha^2+mu delta^2}
psi=tan( {mu alpha} / {mu delta} )     ( {mu delta} >= 0″ title=”{mu delta} >= 0″/>  で、I、IV象限、  <img src=  でII、III象限)

視線速度R、接線速度v 単位は[光年/年]:

R=0.0000003336 * R0
v = 0.000001581 * {{sin mu} / { sin pi}}

距離角sigma 移動後の距離d prime

計算期間tは、年を単位とする
現在の距離d = a / {sin pi}     a = {149600000 km}  / {9460000000000 km} として
d prime sin sigma = v * t
d prime cos sigma = (d + R * t)
tan sigma = {d prime sin sigma} / {d prime cos sigma}
d prime = {d prime cos sigma} / {cos sigma}

t年後の恒星の位置alpha prime delta prime

cos delta prime sin( alpha prime - alpha) =             sin sigma sin psi
cos delta prime cos( alpha prime - alpha) = cos delta cos sigma - sin delta sin sigma cos psi
sin delta prime           = sin delta cos sigma + cos delta sin sigma cos psi
( alpha prime - alpha)は、現在との差なので注意

[天文計算入門 P112-116  天体の位置計算 P43-46]

  1. こんな小さい数字ばかりで、アンダーフローしない? と気にされる方もいらっしゃるでしょう。パソコンであれば、浮動少数で有効数字は確保できてますし、乗算で使う数値なので、桁落ちは起きません。逆にkm/s単位では、オーバーフローは?という話にもなりますね