A Spitzer Space Telescope Exploration Science Program to Search for Y Dwarf Variability

技術要約：Y型矮星の変動性を探索するためのスピッツァー宇宙望遠鏡探索科学プログラム

問題提起
超低温矮星（スペクトル型M7以降）の大気は、原子による不透明度よりも、分子および凝縮物の不透明度によって支配されている。モデルは、L型およびT型矮星における光度変動が回転によって変調される不均質な雲構造に起因することを示唆しているが、Y型矮星は、その本質的な光度の暗さ（光学域および近赤外線波長において）のために研究が困難であった。Y型矮星に関するこれまでの中間赤外線研究は、サンプルサイズの小ささや光度測定の不確実性の大きさによって制限されてきた。さらに、L型およびT型矮星における変動性は雲の進化に関連付けられている一方で、内部熱流束が低く散逸時間スケールが長い、最も低温のY型矮星における変動を駆動するメカニズムについては、依然として制約が少ない。本研究の主な目的は、包括的なY型矮星のサンプルを用いて中間赤外線光度変動を探索・特性評価し、変動率、回転周期、および振幅の傾向を決定すること、そして、雲構造の変化がブラウン星のスペクトル系列全体にわたるこの変動の主要な駆動源であるという仮説を検証することである。

手法
著者らは、提案時（2012年）に分光的に確認されていた14個のY型矮星を対象とした、スピッツァー宇宙望遠鏡探索科学プログラム（プログラム#90015）を実施した。スピッツァーで観測されたすべてのY型矮星からなる完全なサンプルを作成するため、WD 0806−661BとWISE J085510.83−071442.5の2つの天体を追加し、合計サンプル数を16個とした。

• 観測: 対象天体は、赤外線観測装置（IRAC）を用い、[3.6] $\mu$mおよび[4.5] $\mu$mの両方の波長で「ステアリングモード」にて観測された。回転周期をカバーするため、各天体は24時間（各バンド12時間ずつ）連続してモニタリングされた。観測は、長期的な安定性とライトカーブの形態変化を調査するため、122日から464日の間隔を空けて第2エポックでも繰り返された。
• データ・リダクション: 基本較正データ（BCD）フレームを総電子数に変換した。重心決定には2次元ガウス関数を用い、アパーチャ測光は2ピクセル半径で行った。背景レベルは、汚染を避けるためにアニュラス（環状領域）ではなく、ソースフリーの領域から決定した。中央絶対偏差（MAD）の50倍を超える外れ値は除去された。
• 変動性解析: 変動源を特定するために、2つの異なる統計的手法が用いられた。
  1. ベイズモデル比較: マルコフ連鎖モンテカルロ（MCMC）サンプリングを用い、定常フラックス、単一正弦波、および二重正弦波の3つのモデルを比較した。ベイズ情報量基準（BIC）を使用して、最適なモデルを選択した。定常モデルと正弦波（または正弦波と二重正弦波）の間の$\Delta$BICが6を超えた場合に、そのソースを変動源とみなした。
  2. 周期図解析: 各ライトカーブに対してロンブ・スカーグ（Lomb-Scargle）周期図を生成した。パワーが5%の偽アラームレベル（95%信頼区間）を超えたソースを変動源として分類した。
• 完全性と信頼性: 頑健な変動率と信頼区間を導出するために、著者らはベイズ形式論（Lafrenière et al. 2007）を採用した。ランダムな位相、振幅、および周期を持つ10,000個の合成ライトカーブを生成し、解析パイプラインによって検出可能な割合を決定することで、各ターゲットのサーベイ完全性（$p_i$）を計算した。

主要な結果

• 変動率: 本研究により、かなりの割合のY型矮星が中間赤外線での変動を示すことが判明した。
  • エポック1: 変動率は[3.6]で$59\% \pm 15\%$、[4.5]で$64^{+10}_{-13}\%$であった（ベイズ解析ではそれぞれ19%および56%と低い値が得られ、手法への依存性が示された）。
  • エポック2: 変動率は[3.6]で$35^{+17}_{-11}\%$、[4.5]で$75^{+8}_{-15}\%$であった。
  • 変動率は[3.6]よりも[4.5]において一貫して高い。これは、Y型矮星がより明るい[4.5]において信号対雑音比（S/N比）が高いためと考えられる。
• ライトカーブの安定性: 多くのライトカーブは数ヶ月のタイムスケールにわたって安定していたが、3つの天体（WISE 0713−29、WISE 1405−55、WISE 1639−68）は、エポック間で振幅および／または形状の明らかな変化を示した。特にWISE 1405−55は、第1エポックの低レベルの変動から、第2エポックでの顕著な変動への移行、および[4.5]におけるライトカーブ形状の劇的な変化を示した。
• 回転周期: 複数のバンドおよびエポック間で一貫した周期から導出された、2.44から8.42時間の範囲にある5つの確実な回転周期を持つY型矮星が特定された。追加で3つの天体が暫定的な周期を有している。
• 振幅の傾向: 本研究は、L型およびT型矮星で観察される「スペクトル型が減少するにつれて最大変動振幅が増加する」という傾向が、初期型のY型矮星にも及んでいることを裏付けている。大振幅の変動（>3%）は、L型およびT型矮星よりもY型矮星においてより一般的である。
• 位相解析: 両方のバンドおよび両方のエポックで強い変動を示したWISE 0359−54については、[3.6]と[4.5]の間で位相差は検出されなかった。これは、これらのバンドが同様の低圧大気層を捉えていることと一致している。

意義および主張
本論文は、L、T、およびY型矮星の変動率（本研究とS. A. Metchev et al. 2015によるL/T型矮星のサーベイを組み合わせたもの）が、ブラウン星の変動は凝縮雲の水平および／または垂直構造の変化によって引き起こされるという標準的なパラダイムを弱く支持していると主張している。具体的には、L型での高い変動、雲の消失に伴うT型での低下、そして（KCl、Na$_2$S、水氷などの新しい凝縮物の形成に関連する）Y型での変動の回帰という観察されたパターンは、「雲が多い $\to$ 雲がない $\to$ 雲が多い」というスペクトル系列の進化と一致している。

著者らは、変動の検出が解析手法（ベイズ法 vs 周期図法）に大きく依存するものの、Y型矮星の高い変動率に関する全体的な結論は強固であることを強調している。これらの結果は、ジェイムズ・ウェッブ宇宙望遠鏡（JWST）を用いたY型矮星の変動に関する将来の光度および分光研究のベースラインとなるものである。本論文は、変動の物理的メカニズムを決定的に証明したと主張するものではなく、観察された傾向が雲駆動型のモデルと整合していることを述べるにとどめている。