これの続きです。
去る９月１７日に発見した新変光星、KaiV85を今日VSXに登録申請したら、すぐにOKが来ました。これは速い！わずか３日で登録完了。
https://www.aavso.org/vsx/index.php?view=detail.top&oid=624653
周期もEBで短いので全位相はすぐに撮れました。15.3-16.4等と少々くらいですがAPASSのデータが意外とたくさんあって、それもだいたいうまくライトカーブに載せられました。
KaiV85は記念すべきVSX登録５０個目の新変光星でした。
天文とも音楽とも関係のないテーマです。
一週間くらい前に3 SATテレビの番組で紹介されていて、大変に驚きました。
芸術、物理学、工学、解剖学などの多様な要素を組み合わせた素晴らしいものだと思います。ちょっとLeonardo da Vinciのような？あるいはどこか天空の城ラピュタに通じるような？
NASAから問い合わせがあったそうです。金星での探査機は太陽電池が使えそうもないし、どんな地面かわからないので、そこで動かせる機構を開発してもらいたいとのことです。
これの続きです。
このあたりにまだ変光星がありそうな気がして、少し西の方へ写野をずらして見ました。前回見つけたKaiV84は写野に入れたままで、二つの明るい星は写野の外にし、露光時間は2倍の60秒にして一晩中観測。
そうしたら案の定一つきれいなのが受かりました。
これも変動幅が大きいので、まあVSXに載っているだろうとあまり期待しないで星の座標を入れて、Searchボタンを押したらそこには既知の変光星はないと出ました！
これもASAS-SNは見つけていない！このところ続けざまに５つ見つかりましたが、そのうち３つはかなり変動幅が大きいものです。まだまだチャンスはありそうですね。頑張ろう。
15日に見つけたKaiV84ももちろん健在です。
倉敷の赤澤さんがずっと前に発見された新変光星を再調査なさるというので、どんな星か私も一度観測してみることにしました。
9月15日の観測ではその星の短期変動は結局受かりませんでしたが、その近くに一つ新変光星を見つけました。
KaiV84
結構大きな変光を示す上に13等台とまあ明るい方なので、これはもうASAS-SNの変光星カタログに載っているだろうと思いながらVSXでチェックしてみましたが、なんと記載されていません。嬉しいです！
いまやASAS-SN Sky Patrolのトロール網で、私が見つけられるような変光星はもう全部見つかってしまったと思っていましたが、まだまだその網に引っかかっていない星もあるようです。安心しました。これからも探してみよう。
の続きです。
この二晩の観測フレームから新変光星が見つかりました。
KaiV83 9月10日の観測
9月10日と11日を並べてみます。
二晩でかなり違うライトカーブが得られたので、比較星との色の違いなどではないことがはっきりします。ELLタイプでしょうか？
古い観測の記録をみたらこの星が変光星ではないか?と疑問符がついていました。
今vsoljのMLに報告しましたが、ブログにも記録しておきたいと思います。ご覧になったかたもいらっしゃるかもしれませんが、ご容赦ください。
もうBJDを使う時代かもしれませんが、まだまだHJDが普通に使われていますね。
日本変光星研究会編の「変光星観測」の１４１ページには日心補正の計算には補正をするソフトウェアを利用するのが一般的ですとありますが、JDをHJDに変換する時に何で皆さんはやっておられますか？
私がNgaさんに日心補正をどうやってと、尋ねましたら、このエクセルシートを使えばと送ってくださったのがBruton氏のサイト（下記参照）からダウンロードできるエクセルシートでした。
今回、このエクセルシートやインターネットで見つけられた他の手段を使って３つの場合について実際に計算して比較してみました。
下の表を御覧ください。
Date & Timeは仮想の観測日時と時刻です。最初の２つの観測日時は冬至と夏至の頃で、太陽と正反対の方向にある仮想の星を観測した時に補正値がどうなるか見ました。３番目のは黄道面から離れた星の場合を見ました。
JDはその時刻のJDです。
Objectは仮想の対象星です。
dt=HJD-JDは日心補正値です。
ConverterはJDをHJDに変換する手段で
Kaiは私が作成したR言語によるスクリプトで地球の位置を地球の楕円軌道のケプラー方程式をニュートン近似で解をもとめて補正値を計算しています。基本的にはBrutonのマクロの計算方法を手直しし、定数を改訂したものです。
Osmは大島さんが作成されたエクセルシートです。「Henden&Kaitchuck(1982) Astronomical Photometry のP.114の式を使っています。ケプラー方程式は直接には解いていませんが、（和書で水路部「こよみ便利帳」にあるような）多項式近似で動径を求めているようです。」とのことです。http://otobs.org/hiki/?heliocentricJD
AIJはAstroImageJに付属しているCoordinate Converterです。
BAAはこのサイトです。 http://britastro.org/computing/applets_dt.html
L.Cookはこのサイトで計算させることもできますし、エクセルをダウンロードすることもできます
http://www.lewcook.com/JDSUNTime.htm
Brutonはこのサイトで計算させることもできますし、エクセルをダウンロードすることもできます。 http://www.physics.sfasu.edu/astro/javascript/hjd.html
astrolibRはR言語の天文プログラムのパッケージでその中の変換プログラムです。
https://rdrr.io/cran/astrolibR/man/helio_jd.html
Bruton.revは彼のエクセルシートのマクロをホンのわずか修正した場合です。（下記参照）
結果を見てみました。
冬至の頃は地球は近日点の近くで夏至の頃は反対に遠日点の近くなので、その距離の差が補正値に現れるはずですが、
BrutonとastrolibRではどちらの場合も補正値が４９９秒と全く違わないので円軌道で計算していることがわかります。
他の５つでは距離の差が補正値に反映していて、しかも0.1秒台で一致しています。
astrolibRのにはびっくりしました。パッケージに入っているのが信じられません。他のプログラムでも問題があったことがあります。(JDに入れる値はしかもJD-2400000です。要注意。まあ使わないほうが良いプログラムです。）
Brutonのマクロでは地球の楕円軌道をケプラー方程式で解いて地球の位置を出しているのですが、得られた太陽と地球の距離reを最後の最後で入れ忘れていて実際には円軌道扱いのように見えるのです。
彼のマクロのこの行
correction = ausec * (earthx * objectx + earthy * objecty + earthz * objectz)
にreをこのように入れます。
correction = ausec * re * (earthx * objectx + earthy * objecty + earthz * objectz)
そうすればBruton.revにあるように他のconverterと極めてよく合うようになります。
もしBrutonのエクセルシートをお使いの場合はこのようにマクロを修正なさることをおすすめします。
なお、AstroImageJのConverterだとBJDも計算させられます。
Date & Time 2011.12.22. 0 UT 2011.06.22. 0 UT 2018.09.02. 12 UT
JD 2455917.5 245734.5 2458364.0
Object RA:06 00 00 Dec:23.5 RA:18 00 00 Dec:-23.5 RA:19 16 50 Dec:29 10 52
Converter dt=HJD-JD in sec dt=HJD-JD in sec dt=HJD-JD in sec
Kai 0.005681 490.8 0.005869 507.1 0.002706 233.8
Osm 0.005681 490.8 0.005869 507.1 0.002736 236.4
AIJ 0.005682 490.9 0.005869 507.1 0.002705 233.7
BAA 0.00568 490.9* 0.00587 507.1* 0.00271 233.8*
L.Cook 0.005681 490.9** 0.005869 507.1** 0.002706 233.8**
Bruton 0.005775 499.0 0.005776 499.0 0.002682 231.7
astrolibR 0.005774 498.9 0.005774 498.9 0.002611 225.6
Bruton.rev 0.005682 490.9 0.005870 507.2 0.002706 233.8
* from dt in minutes
**from dt in seconds
これの続きです。
２０２２年に見られるかもしれないというので、もうドキュメントの映画が作られていたようで、ここで見られます。TrailerとかVideoをクリックして見てください。