How X-rays heat the IGM in different 21-cm simulation codes: a comparison between Licorice and Beorn

宇宙を、主に水素でできた、あらゆる方向に広がる巨大で目に見えないガスの海だと想像してみてください。ビッグバンの後、最初の数百百万年間、この海は暗く、冷たく、静まり返っていました。その後、最初の星や銀河が、霧の中の灯台のように点火しました。それらの強烈な放射線はガスを加熱し始め、原子から電子を引き剥がし始めました。これは「再電離」と呼ばれるプロセスです。天文学者たちは、この時代がどのように展開したかを知る手がかりとなるため、この時代を観察しようと切望しています。これを行うために、彼らは「スクエア・キロメートル・アレイ（SKA）」と呼ばれる巨大な電波望遠鏡を建設しています。これは、微かな囁きである「21cm信号」を聴き取るための巨大な耳として機能します。この信号は、中性水素から放出される特定の電波であり、その強さはガスの温度によって変化します。ガスが冷たい場合と温かい場合では、信号の見え方が異なります。この信号を測定することで、科学者たちは宇宙の加熱の歴史をマッピングしたいと考えています。

しかし、問題があります。宇宙はあまりに広大で複雑であるため、リアルタイムで直接観察することはできません。私たちは、数十億年前の出来事をシミュレートするためのコンピュータモデルを構築しなければなりません。これらのシミュレーションは、いわばデジタルなタイムマシンです。しかし、どんなタイムマシンも、それが従うルール次第でその精度が決まります。すべてのガスの粒子や光の線を追跡する非常に詳細なシミュレーションもありますが、それにはスーパーコンピュータで数週間かかります。一方で、より高速で、ショートカット（簡略化）を用いるものもあり、使いやすい反面、精度が低い可能性があります。大きな疑問は、これらのショートカットが、結果を信頼できるほど十分に機能するのかどうかということです。もし2つの異なるシミュレーションコードが、ガスがどれほど熱くなったかについて異なる物語を語ったとしたら、どちらを信じるべきでしょうか？ロマン・メリオット氏とそのチームによる新しい研究が解決しようとしているのは、このパズルです。

チームは、X線による宇宙ガスの加熱をどのように扱うかを検証するために、2つの非常に異なる「タイムマシン」を直接対決させることにしました。一方の側には、高精細な3DシミュレーションであるLicoriceがありました。これはフルモーションのビデオゲームのように機能し、3次元でガスがどのように動き、光がどのように伝わるかを計算します。もう一方の側には、より高速で単純な、「1次元的」なアプローチを用いるBeornというコードがありました。Licoriceを、個々の材料を一つずつ混ぜ合わせながら複雑なシ Stew（シチュー）を作る熟練のシェフとするならば、Beornは、平均的なレシピに基づいて食事全体を加熱する電子レンジのようなものです。研究者たちは、電子レンジが、シェフのシチューと同じくらい美味しい食事を作れるかどうかを確認したかったのです。具体的には、最初の星々からのX線が水素ガスをどのように温めたかという点においてです。

公平な戦いにするために、科学者たちは単にそれぞれのコードにランダムな設定で走らせたわけではありません。代わりに、詳細なLiciceの「材料」（星々とその光の出力）を取り出し、それをBeornに投入しました。これにより、両方のコードが、異なる手法を用いながらも、全く同じ宇宙を作ろうとしていることが保証されました。そして、ガスの温度、21cm信号の全体的な強さ、そしてその信号が空のいたるところでどのように変動するかという統計的なパターンを比較しました。

結果は、朗報と警告が混ざり合ったものでした。全体像、つまりガスの平均温度や信号の全体的な強さを見たとき、2つのコードは非常によく一致していました。それは、電子レンジもシェフも、シチューを適切な平均温度に仕上げることに成功したかのようでした。しかし、細部に目を向けると、違いが現れ始めました。熱の分布が同じではなかったのです。詳細な3Dシミュレーションでは、熱は斑状（パッチ状）であり、非常に熱い場所もあれば冷たい場所もあるという複雑な質感を持っていました。一方、より高速な1Dシミュレーションでは、熱はより滑らかで均一でした。

この熱の「質感」における小さな違いが、驚くほど重要であることが判明しました。チームが「パワースペクトル」（異なるスケールで信号がどの程度変動するかを測定する高度な方法）を計算したところ、2つのコードの間には約**30%**の相違が見つかりました。これを例えるなら、一口食べただけでケーキのレシピを推測しようとしている場合、味に30%の差があれば、材料の推測を完全に間違えてしまうということです。

研究チームは、なぜコードが一致しないのかを突き止めるためにさらに深く掘り下げました。彼らは、高速なコードで使用されているショートカットこそが原因であることを発見しました。具体的には、1Dコードはガスの密度がどこでも一定である（例えるなら、海が完全に平坦であると仮定すること）と想定し、X線が実際にどれだけの熱を注入するかについても簡略化された公式を使用していました。実際には、ガスは塊（クランプ）状になっており、加熱はガスの正確な位置に依存します。チームが詳細な3Dコードに対して、これらと同じ簡略化されたショートカットを強制的に適用させたところ、結果は1Dコードに近づきました。これにより、これらの近似がエラーの源であったことが裏付けられました。

最後に、チームは究極の問いを投げかけました。「もし私たちが将来のSKA望遠鏡のデータを分析するために、これらの異なるシミュレーションを使用した場合、間違った答えを得ることになるのだろうか？」彼らは、望遠鏡が100時間観測するというシナリオをシミュレートしました。その結果、コード間の違いが、最初の星々の特性の推定に重大なバイアスをもたらすことが分かりました。統計学的に言えば、その誤差は1シグマ（信頼度の指標）よりも大きく、パラメータによってはそれ以上でした。これは、もし天文学者が、現実のデータを解釈するために、より高速で単純なコードを使用していたとしたら、最初の星々が実際とは異なっていると自信を持って結論付けてしまう可能性があることを意味しています。

この研究は、高速なコードは有用であり、「平均的な」物語は正しく捉えているものの、SKAの時代においてはまだ精密さが足りない、と結論付けています。パワースペクトルの30%の差と、それに伴うデータのバイアスは、極めて高い精度で宇宙を測定しようとしているとき、単にショートカットを信頼することはできないことを示唆しています。著者らは、現在のより感度の低い望遠鏡にはこれらのコードは十分通用するものの、SKAが聞き取りを開始する前に、高速なコードを洗練させるか、あるいはそのエラーを修正する新しい方法を開発するために、多くの作業が必要であると示唆しています。これは、宇宙の夜明けを理解しようとする探求において、細部にこそ真実が宿るということを思い出させてくれます。