ASAS113804-0026.5の多重周期がようやく見つかった!
ようやく今年の観測からフーリエ解析でこのDSCTの多重周期が見つかりました。
結論から行くと、2015年の観測データをPDM解析して見つかっていた周期、0.06100日と、昨年並びに今年の観測データをPDM解析して見つかっていた周期、0.05805日両方が、今年の観測データに入っているのです。
今年の観測データをdetrendしないでも、この二つの周期で位相図を作るときれいにできます。
まず今年の観測のPDM解析で見つかる周期で。
もう一つの周期で。
2015年の観測データから2番めの周期で位相図を作ると
というわけで、2015年には2番めの周期が強かったのが、2016年2017年には1番目の周期が強かったようです。
DN Ophのデータ改訂をVSXに登録
1928年発見のこの変光星。ほぼ90年ぶりにようやくタイプや周期が判明してカタログに記載できました。
https://www.aavso.org/vsx/index.php?view=detail.top&oid=20622
これもデジカメおじさんのおかげです。ありがとうございました!
ASASデータのMAG_0やMAG_1などの扱い
セバスチャンが返事をくれてMAG_Nについて説明してくれました。
規則
星の光度が12等より暗い場合は MAG_0
11等から12等の間では MAG_1
10等から11等の間では MAG_2
9等から10等の間では MAG_3
9等より明るい場合は MAG_4 の列を使うこと
だそうです。
えーっ、今までこんなことを考えておりませんでした。
ただ、調査してきた星はだいたいくらいのが多かったのでどちらにしてもMAG_0の列を使っていたと思います。
そこで今も研究中のASAS J113804-0026.5の場合は10等台なのでMAG_2の列を使わなければならないと、データを見直してみることにしました。
そうしたら、なんとMAG_2の列が先頭に来ています!
ASASデータの方でちゃんと正しい列を先頭にして他のを後ろにしてくれているのです。
ASASデータをHJD_magの表に変換する自作のRスクリプトでは読み込んだデータの第2列を光度として使うようにしているので、知らず知らず自動的に正しいデータを使っていました。良かった。ほ~っ。
DN Oph
デジカメおじさんのブログでDN Ophが3月に明るかったという報告を読みました。
デジカメおじさんは、VSXには変光星のタイプも周期も不明となっているが、DN OphのASASデータでははっきりミラ型の変光星であることがわかることも指摘されています。
そこで、ASASデータをダウンロードして周期をPDMで求めました。周期は216.8日。
ただ、極大光度が変化しているのでどの値を極大光度にするかとか、極小光度はASASではわからないのでどうするべきかなど、セバスチャンにメールを書きました。
そうしたら極大光度は一番明るくなったときの極大光度を採用、極小光度はVizieRのデータにあるSPM4.0に載っている17.4Vを使おうとのことでした。
彼が指摘してくれたことで、前から気になっていたのですが、ASASデータのMAG_0、MAG_1、MAG_2などの列の扱いです。
これはDN OphのASASデータの一部です。
1つのHJDにつき、MAG_0からMAG_4までとMER_0からMER_4まであります。
MAGは各アパーチャーに対応してるとありますので、アパーチャーの半径を変えて測光したものと理解しています。(本当にそうかは知りません。)
そしてMERはそのアパーチャーに対応したエラーの値と理解しています。
以前から気がついているのですが、MERを見るとMER_2の値がほぼいつも他のMERの値より小さいのです。そこで僕はASASのHJD mag errの表に直すスクリプトでは基本的にはMAG_2とMER_2を使うようにしていますが、場合によってはMAG_0とMER_0を使えるようにもしています。
今回はMAG_0の方で位相図を作って元期もそのデータ中一番明るい日を選びました。
下の表はASASデータで一番明るい日の付近の抜粋で、MAG_0で11.544等の日、4572.82です。
HJD MAG_0 MAG_1
そしたら、セバスチャンが光度が12等より明るい場合はMAG_1を使うようにと言ってきました。MAG_1ではHJD 4574.79の11.594の方が明るいのです。
どうして?とメールを今彼に書いたところです。
変光範囲が小さいEWとELLの見分け方
ちょうど1年前、IR Gemがアウトバーストして、それを観測した時に近くの星が変光していることに気が付き、KaiV36と番号をふりました。これが当時のブログ。
先日、病気になる前に何晩か観測することができ、そのデータと昨年のデータなどを組み合わせて正確な周期を出し、位相図も一応満足の行くものができたのでVSXに登録申請しました。前もってセバスチャンにも意見を聞いておいたので問題なくすぐにOKが来ました。これが位相図です。
セバスチャンに最初に意見を求めた時に、私はタイプを変光範囲が小さいEWと思って書いたのですが、彼からはっきりとこれはELLと来たので、「2ヶ月くらい前にKaiV60というのを登録したがあれはEWと言われたけれど、これもそれと似てるんでは?と質問しました。そうしたら、
「EWは極大がshallow浅く、極小はsharp尖い、ELLはその点で対称だ。
EWの周期は0.3日あたりが典型的。」
などと説明してくれました。下は変光範囲が小さいEW タイプのKaiV60の位相図です。
なるほど、こう並べてみると違いははっきりしていて言われたとおりです。これからは一人で違いを見分けられるだろうか?
初めてOV Booを観測
ほぼ一週間ぶりのブログと観測です。風邪を引いてしまって休んでいました。その間天気はひどかったので、「観測できるのにー!」というストレスは幸いありませんでした。昨晩、体調も天気も回復したのでOV Booが上がってくるまでOQ Gemを観測して時間を潰しました。(しかし食外の観測もたいせつです。今晩また第一極小があります。)
VSOLJの皆さんの観測のライトカーブを見て知ってはいましたが、すごい変化に富んだライトカーブにびっくりしています。
昨晩、例年のように冬時間から夏時間への切り替わりがありました。観測のコンピューターがどうなっていたか、ヨーロッパ時間午前2時前後のフレームの時刻を見てみましたが無事にUTのまま継続してFITSに入っていました。(ヨーロッパ時間は2時から3時に飛ぶ)
昨10月の夏時間から冬時間への切り替わりの時に日付を間違えて報告したりして、Katさんにご迷惑をおかけしたのを思い出します。
雲間からOQ Gemの復光を観測できた
久しぶりの極小がある日ですが、天気予報は曇り。でもところどころ雲間から薄雲があるものの、星が見えます。なんとか極小から増光しているところをとらえられました。
これで周期はまず、7.2160日が確定。
風邪気味なので切り上げて寝ることにしよう。