宇宙の塵の塊の「跳ね返り」が衝突合体による微惑星形成を阻害する 大きくなるとくっつきにくくなる粉状体の衝突挙動を発見

国立研究開発法人海洋研究開発機構（以下「JAMSTEC」）付加価値情報創生部門 数理科学・先端技術研究開発センターの荒川 創太Young Research Fellow、東京工業大学 理学院 地球惑星科学系の奥住聡准教授、辰馬未沙子研究員、中本泰史教授らは、国立天文台の計算サーバ^[用語2]を用い、惑星の材料物質である固体微粒子の塊について、その衝突挙動を数値シミュレーションによって調べました。様々な大きさの塊について数値シミュレーションを実施した結果、塊が大きい場合に2つの塊が衝突合体する確率が低下することを明らかにしました。

「惑星の種」である微惑星は、原始惑星系円盤^[用語3]において固体微粒子同士が付着による衝突合体を繰り返すことで形成されたと考えられています。しかし、どのような条件の下で2つの微粒子の塊が衝突合体するのかという重要な問題が未解明のままでした。本研究の結果は、微粒子の塊が大きくなるにつれて衝突合体しにくくなるため、微粒子同士の衝突合体による成長のみで微惑星を形成することは困難であることを示唆しています。これは惑星形成プロセスを理解する上で重要な知見となります。

本成果は、「The Astrophysical Journal Letters」に7月6日付け（日本時間）で掲載されました。

惑星形成の第一歩は、原始惑星系円盤におけるマイクロメートルサイズの固体微粒子の成長です。その成長メカニズムのひとつとして衝突合体が考えられています。しかし、微粒子同士が衝突合体してできた塊がさらに合体して成長するための条件は長年謎でした。また、数値シミュレーションと室内実験の結果に大きな乖離があり、特に、数値シミュレーションと室内実験とでは、衝突時の跳ね返り現象を引き起こす塊の密度の条件が大きく異なっていました。例えば、これまでの数値シミュレーションでは塊の内部の空隙の割合が50%程度以下という高密度の場合にのみ跳ね返り現象が頻繁に見られた一方、室内実験では空隙の割合が90%程度という低密度の場合でも高い確率で跳ね返ることが報告されていました。

「惑星の種」である微惑星の形成モデルには、固体微粒子の衝突合体による成長という仮説と、微粒子の塊が原始惑星系円盤中で局所的に濃集し、その後、自己重力によって集積し形成されるという仮説の2つがあります（図1）。もし、衝突時に微粒子の塊が跳ね返らず付着するならば、衝突合体による成長によって微惑星が形成される可能性があります。一方、跳ね返りによって衝突合体が阻害される場合、衝突合体以外の成長メカニズムで微惑星を形成する必要があります。微粒子の塊の衝突挙動を理解することは、これらの2つの仮説のどちらが微惑星形成シナリオとして妥当か判断する鍵となります。

Image may be NSFW.
Clik here to view.

図1.

原始惑星系円盤中では、マイクロメートルサイズの固体微粒子が衝突合体によって巨視的な塊へと成長し、さらに微惑星と呼ばれる惑星の種の状態を経て惑星へと進化する。微粒子の塊から微惑星へと成長する段階では、衝突合体成長によってそのまま大きくなるという仮説と、原始惑星系円盤中での自己重力による微惑星形成という仮説の2つの可能性がある。（クレジット：JAMSTEC）

本研究では離散要素法^[用語4]を用いた衝突数値シミュレーションを実施し、構成粒子数が約1万から14万までの様々な大きさの塊について、衝突時の付着確率を調査しました。その結果、微粒子の塊は空隙の割合が同じならば半径が大きい方がより付着確率が低いことが明らかになりました。このシミュレーションでは、惑星が形成される原始惑星系円盤を想定し、塊は半径0.1 µmの氷の微粒子で構成されています。ここでは、空隙の割合が60%の塊同士をコンピュータ上で衝突させました。すると、塊の半径が微粒子半径の50倍以下の場合、塊同士が付着する確率は約90%にものぼりました。一方、それよりも塊が大きい場合は跳ね返りやすくなり、微粒子半径の70倍の大きさをもつ塊では付着確率は約50%まで下がったのです（図2）。

Image may be NSFW.
Clik here to view.

図2.

2つの微粒子の塊の衝突数値シミュレーションの様子と、付着確率と塊の半径の関係。半径0.1 µmの氷微粒子で塊が構成され、空隙の割合が60％の場合、塊の半径が微粒子半径の50倍よりも大きくなると付着確率が顕著に低下する（跳ね返りやすくなる）ことが示された。（クレジット：S. Arakawa et al. ）

今回の発見は、過去の数値シミュレーションと室内実験の結果の乖離についても定性的な説明を与えるものであり、大規模シミュレーションの重要性を示しています。過去の数値シミュレーションでは、構成粒子数が数万程度で、塊の半径が微粒子の半径の数十倍の比較的小さな塊について研究が行われてきました。一方、室内実験では、粒子数が数億以上で微粒子半径の数百倍から数千倍の半径をもつ塊を用いた研究が行われてきました。つまり、従来の数値シミュレーションと室内実験とでは用いられてきた塊の構成粒子数に大きな差があり、これが両者の衝突挙動の違いをもたらしていたと考えられます。

本研究の成果は、原始惑星系円盤での固体微粒子の衝突合体成長は、塊がある程度大きくなった段階で付着確率が低下し抑制されることを示唆しています。つまり、微粒子の衝突合体成長のみで微惑星を形成することは困難であり、微粒子の塊の原始惑星系円盤中での局所的な濃集など、微惑星の形成にはその他のプロセスの助けを借りる必要があるかもしれません。実際、アルマ望遠鏡^[用語5]による原始惑星系円盤の電波観測では固体微粒子の塊は100 µm程度の大きさで成長が止まっている可能性が指摘されており、本研究はこの観測結果に説明を与える可能性があります。

微粒子の塊の衝突挙動は、粉状体という基礎物理の対象として興味深いだけでなく、惑星形成という地球惑星科学の根幹をなす問題においても鍵となるテーマです。今回の研究では、塊の半径が大きいほど、つまり構成する微粒子の数が多いほど塊が“くっつきにくくなる”ことを示しました。これは衝突合体のみで微惑星が作られることは困難であるという、惑星形成の理解に重要な示唆を与える発見です。

しかし、本研究で用いた塊とは異なる密度の塊について、跳ね返りが起こる条件は未だ明らかではありません。例えば、本研究では塊内部の空隙の割合を60%で固定していましたが、現実の原始惑星系円盤中における塊内部の空隙の割合はもう少し高く低密度な可能性があります。そのような低密度の塊でも十分大きい場合には本研究の結果と同様に跳ね返りが生じると考えられますが、定量的な理解は得られていません。今後我々はJAMSTECが保有する地球シミュレータなどのスーパーコンピュータと大規模シミュレーションコードDEPTH^[用語6]を活用し、現時点では調査できていない室内実験と同規模の非常に大きな塊を用いた衝突数値シミュレーションを実施することで、跳ね返り現象を支配する物理の解明を目指します。

本研究で得られる成果は、惑星科学における新たな知見の創出のみならず、地球科学における流体や固体のシミュレーション技術の発展のほか、粉体などの粒子状物質のデジタルツインによる産業の発展においても大きく貢献することが期待されます。