特任助教(国立天文台フェロー)らの研究チームはアプリケーション「Bonsai」を開発し、18600台のGPU1を用いた計算により、これまでで最大規模の天の川銀河進化の数値シミュレーションを成功させました。計算に用いた粒子数は約2400億であり、この計算により初めて、天の川銀河の星の観測データと直接比較可能なシミュレーションデータを得ることができるようになりました。(1秒間に2.477京回の計算)を達成しました。このシミュレーションの性能と科学的意義が評価され、2014年のゴードン・ベル賞5のファイナリストとして選出されました。本研究のシミュレーションで計算された銀河の姿。左半分の緑色の部分は、重力計算のための領域分割を可視化したもの。私たちの棲む太陽系は、天の川銀河の銀河円盤の中に存在しています。そのため銀河円盤を常に横から見ている状態にあり、天の川銀河がどのような姿をしているのかを直接見ることができません。銀河内の星の位置やその速度を正確に測ることで、棒状構造や渦巻き腕のような銀河の構造を知ることができます。JASMINE」(ジャスミン)7によって、まもなく天の川銀河内の星の位置と速度がこれまで以上に正確にわかってきます。この観測データを読み解くためには、星の運動が銀河のどのような構造を反映しているかを知っておく必要があります。これには銀河の構造と星の運動の関係を直接知ることができるシミュレーションが重要な役割を果たします。観測とシミュレーションと合わせることによって初めて、天の川銀河の構造とそのメカニズムが明らかになるのです。ハッブル宇宙望遠鏡で観測された渦巻銀河UGC12158。天の川銀河はこの銀河に似た形をしていると考えられている。中心付近に棒状の構造があり、その周りに渦巻き腕が取り巻いている。)実際の天の川銀河を構成する星は2000億~4000億個と言われています。シミュレーションと観測を比較するためには、本来、膨大な数の星や銀河を取り囲むダークマター8の間に働く重力を計算し、個々の星の運動を導きださなければなりません。これには莫大な計算能力を持ったスーパーコンピュータが必要になります。これまでのシミュレーションでは、1000個程度の星を1つの粒子に置き換えるなどの近似を用いた計算が行われて来ました。このようなシミュレーションでは、観測データとの直接比較はできません。来るべき位置天文衛星のデータを理解するためには、より現実的なシミュレーションが必要となるのです。Daint9とアメリカ・オークリッジ国立研究所のTitan10という、CPU11に加えてGPUを搭載したスーパーコンピュータを用いて、これまでで最大規模の天の川銀河の進化シミュレーションを行いました。シミュレーションに用いたGPUの最大の数は、Titanに搭載された18600台になります。これほど多くのGPUを使い、かつ効率よく計算を行うには、様々な工夫が必要となります。一般に、使用するGPUの数が多くなるほど、GPU間のデータ通信に時間がかかるようになり、計算全体の効率は落ちます。2倍の数のGPUを使っても、計算速度は2倍にはならず、それよりも遅くなってしまいます。今回、研究チームが開発したアプリケーション「Bonsai」では、シミュレーションの中でもっとも計算時間のかかる星と星の間に働く重力の計算を全てGPUで行い、通信とそれに関わる計算をCPUで行うことで、GPUとCPUを同時にフル活用できるように工夫しました。そのため、GPUは通信を待つことなく計算を続けることができ、多くのGPUを用いても高い計算効率を維持することができたのです。18600台のGPUを用いて行われた計算の実効性能は、単精度で24.77Pflops(1秒間に2.477京回計算)を達成しました。これは、現在公表されている中では最高の値です。今回の天の川銀河のシミュレーションで用いた星とダークマター粒子の数は合わせて最大約2400億個、星のみでは約200億個です。これは従来の計算より数百倍から数千倍多くの粒子を用いた、まさしく世界で最大規模の天の川銀河シミュレーションとなりました。天の川銀河の姿を知るためには、このシミュレーションで作られた天の川銀河を疑似観測して、実際の観測データと照らし合わせる事が必要になります。Gaiaなどで観測できる実際の天の川銀河の星の数が約10億天体ほどであることを考えると、今回の規模の計算で、初めてシミュレーションから得られた星の数が実際の天の川銀河の観測データの星ひとつひとつと比較できる数になりました。(上)シミュレーションによって再現された天の川銀河を可視化したもの。星を青や白で表している。ダークマター粒子はここでは表示していない。渦巻き腕などの銀河の構造が見えている。観測から得られる星の速度図(c)。横軸が銀河中心からの半径方向の速度、縦軸が銀河回転方向の速度を表す。本研究ではいくつかの条件で計算をしており、(b)では500億粒子(従来の計算の500倍の粒子数)で計算した結果を示している。以前までのシミュレーションの結果(a)に対して、本研究の結果(b)は多くの星の情報を得ることができている。(b)では銀河の構造を反映した速度のパターンを見て取ることができる。今回の研究で行われたシミュレーションは、そのパフォーマンスの良さと、多くのGPUを使って効率よく計算するための手法、そして科学的意義が評価され、2014年のゴードン・ベル賞ファイナリストとして選出されました。(国立天文台フェロー)の藤井通子氏はこの研究において、シミュレーションを行う銀河のモデル作りと計算結果の解析を担当しました。天の川銀河を再現するようなモデルにはどのようなパラメータ設定が最適か、これまでの藤井氏自身の研究成果も踏まえて選びました。現在は、計算結果からシミュレーションが正しく行われているかを評価し、実際の天の川銀河の観測と比較していく上で必要となるデータの解析をおこなっています。センターの石山智明研究員がアプリケーション開発で協力をしました。)、エフゲニー・ガブロフ(SURFsara)が協力しました。藤井氏は今回ゴードン・ベル賞にファイナリストに選出されたことについて、このように語ります。「私は昨年9月までオランダのライデン大学に所属しており、昨年から、今回のプロジェクトに参加しています。今回、2回目の挑戦でファイナリストに残ることができました。オランダのスーパーコンピュータは日本と比べると規模が小さく、いつでも大規模なスーパーコンピュータを使ってコードのテストをできるわけではありません。そこで、筑波大学のGPUスパコン「HA-PACS」12や、国立天文台のCPUスパコン「アテルイ」13も利用しながらコード開発を続けてきました。最終的に、最も大規模な計算を行ったTitanも、利用できる時間が非常に限られた中での挑戦でした。そのような不利な環境を乗り越えてのファイナリスト選出なので、非常に嬉しく思います。今後、データ解析という大変な仕事が待っていますが、天の川銀河の研究にとって大きな一歩となる結果を出せるよう頑張ります。今回のシミュレーションでは、これまでと比べ、非常に多くの数の星を扱っているため、観測とのより詳細な比較が可能となりました。今後、昨年打ちあげられたGaiaのデータとシミュレーションで得られた結果を比較していくことで、天の川銀河の棒状構造や渦巻き腕の位置や大きさ、またそれらの構造の進化過程が明らかになると期待されています。このシミュレーションによって作られる天の川銀河は、観測から直接見ることができない天の川銀河の構造の解明を可能にしてくれます。今回のシミュレーションで得られたデータの解析は現在進行中です。シミュレーションと観測が解き明かす天の川銀河の姿にご期待ください。Unit):コンピュータの中で、画像演算処理を行うための専用のプロセッサ。画像演算処理のためには、浮動小数演算を行う必要があり、そのためにGPUには多数の演算器が用意されている。おけるコンピュータの数値表現の形式の一つで,1つの数値を32ビットで表現するもの。表現できる値の範囲は-3.40282×1038~3.40282×1038で、精度は6桁である。よび、表現できる値の範囲は-1.79769×10308~1.79769×10308であり、精度は15桁である。実効性能:理論性能であるピーク性能に対して、あるアプリケーションを実行したときの計算性能。計算機の実質的な性能とされている。Second):flopsは1秒間に行った浮動小数点数演算の数で、コンピュータの性能を表す指標の一つ。P(ペタ)は1015(1千兆)倍。促進のため1987年に創設された賞。この賞は実用的な並列計算で高い性能を達したことに対して与えられる。達成した演算速度ではなく、科学技術研究に有用な計算の高速化を、並列化や最新のハードウエアの利用、計算法の改良といった総合的なアプローチによって実現したことが評価される。Gaia(ガイア):2013年12月に欧州宇宙機関(ESA)によって打ち上げられた観測機。約20等級までの明るさの星を観測し、約10億の恒星について位置と速度のデータを取得することを目標としている。観測が難しい3万光年先の天の川銀河の中心に近い、バルジと呼ばれる星が混んだ領域を赤外線で観測し、星の位置と速度を測ることで、天の川銀河の地図をつくること目的としている。ダークマター:重力相互作用だけが働く物質で、素粒子としての正体は解明されていない。宇宙初期に存在したダークマターの密度揺らぎが重力相互作用により成長して、至るところにさまざまなサイズの天体を形成し、その中で陽子や中性子といった通常物質を集め、現在の宇宙で観測される銀河のような構造を作ってきたという理論「低温ダークマターモデル」がある。これが、宇宙の構造形成を記述する標準的なモデルとして広く受け入れられている。XC30システムのスーパーコンピュータ。CPUとGPUを組み合わせた異種混合クラスタスパコンに分類される。Pflops(倍精度)。XK7システムのスーパーコンピュータ。CPUとGPUを組み合わせた異種混合クラスタスパコンに分類される。理論演算性能は27.1Pflops(倍精度)。500で2位となった。コンピュータを構成する部品のひとつで、プログラムに従って演算・計算を行う電子回路や演算ユニットのこと。コンピュータのまさしく頭脳にあたる部分である。スーパーコンピュータ。理論演算性能は1.166Pflops(倍精度)。天文シミュレーションプロジェクトが運用する天文学専用スーパーコンピュータ。XC30だが、GPUは搭載されずCPUのみで構成されている。今回の研究に用いられた当時の理論演算性能は502Tflops(倍精度,2014年8月まで)。