Plasticity of Neutron Star Crusts

大幅に遅い歪み速度を用いた第一原理分子動力学法を用いることで、本研究は中性子星の地殻が破壊後に定常的な塑性流動の普遍的な領域を示すことを明らかにしており、この知見は、繰り返される地殻の破壊と再アニーリングがマグネターのバーストやフレアを駆動している可能性を示唆している。

要約

中性子星を、信じられないほど高密度で、猛烈な速さで回転する宇宙の巨人と想像してみてください。その表面のすぐ下には、原子核が非常に密に詰まり、山ほどの大きさの超硬質な角砂糖のような、固体結晶を形成している「地殻」が存在します。長い間、科学者たちは、もしこの地殻を強く押しすぎたり、ねじりすぎたりすれば、単にパキッと折れて壊れ、巨大なエネルギーの放出（星震のようなもの）を引き起こすと考えてきました。

この新しい研究は、強力なコンピュータ・シミュレーションを用いて、この宇宙の地殻をねじったときに正確に何が起こるのかを調べています。研究者たちはただ押しただけではありません。かつてないほど「ずっとゆっくり」と押し、そのプロセスを高精細に観察したのです。

以下に、日常的な例えを用いてその発見を説明します。

1. 「硬い」フェーズ（弾性）

輪ゴムを伸ばす場面を想像してみてください。最初は滑らかに伸び、手を離せば元に戻ります。中性子星の地殻も同様です。わずかなストレス（ねじれや圧縮）を加えると、それは完璧で硬いバネのように振る舞います。

• 発見: もし地殻が単一の完璧な結晶（欠陥のないダイヤモンドのようなもの）であれば、壊れるまでに約**11%まで伸びることができます。もしそれが多くの小さな結晶が集まったもの（多くの小石でできた花崗岩のような「多結晶」）であれば、はるかに早く、約5%**で屈服し始めます。

2. 「壊れる」地点

かつて、科学者たちは、地殻が限界に達すると、粉々に砕けて形を保てなくなると考えていました。

• 旧来の視点: 乾燥した小枝を想像してください。曲げていくと限界に達し、「パキッ」と音がして、バラバラに壊れます。
• 新しい発見: 研究者たちは、「多くの小石」バージョン（多結晶）の場合、単にパキッと折れて終わるのではないことを発見しました。一度5%の限界に達すると、それは壊れてバラバラになるのではなく、**「流動」**し始めるのです。

3. 「ハチミツ」フェーズ（塑性流動）

これが最も驚くべき部分です。地殻が屈服した後、それは崩壊しません。代わりに、濃いハチミツや**温かいタフィー（練り飴）**のように振る舞います。

• 例え: タフィーを引っ張る場面を想像してください。一度十分に引っ張って伸び始めると、それは切れることなく、どれほど強く引っ張っても滑らかに伸び続けます。地殻は「完全な塑性流動」の状態に入るのです。
• 結果: 地殻は、壊れたり硬くなったりすることなく、膨大な量（シミュレーションでは最大60%まで）のねじれや変形に耐えることができます。ただ着実に流れるのです。

4. なぜこのようなことが起こるのか？（自己修復する群衆）

なぜ地殻が「ハチミツ」に変わるのでしょうか？

• 比喩: 混み合ったダンスフロアを想像してください。もし、完璧に整列した群衆（完全な結晶）を無理に押し通ろうとすれば、あなたは行き詰まり、最終的に誰かが倒れるまで群衆が激しく押し返してきます（結晶が壊れる）。
• 新しい洞察: しかし、もしその群衆がすでに少し乱れていて（多くの小さな粒を持つ多結晶）、ゆっくりと押し進められた場合、人々（原子の欠陥）は自らを再配置します。彼らは、群衆がスムーズに動けるような、ちょうど良い「隙間」や「滑らかな経路」を作り出すのです。地殻は本質的に、圧力を処理するために自らを再編成しています。圧力を受け流しながら流れ続けるための、独自の「交通システム」を作り出しているのです。

5. スピードが重要

研究では、どのように「押すか」という速度が非常に重要であることが分かりました。

• 速い押し: もし押しが速すぎる場合（車の衝突事故のように）、地殻は再編成する時間がありません。その場合、地殻は脆いガラスのように振る舞い、粉々に砕けるか、あるいは乱れた無定形の泥状になります。これが、古い、より高速なシミュレーションが異なる結果を示していた理由です。
• 遅い押し: ゆっくりと押した場合（氷河が移動するように）、地殻は内部の「交通」を再編成する時間があるため、ハチミツのように滑らかに流れます。

6. これが星にとって何を意味するか

この論文は、中性子星の挙動は、その内部の地殻がどのような姿をしているかに依存することを示唆しています。

• もし地殻が巨大で完璧な単一結晶であるなら: 膨大なエネルギーを蓄積し、突発的かつ壊滅的に破壊される（星震やマグネターのフレアのようなもの）可能性があります。
• もし地殻が多くの小さな粒で作られているなら: エネルギーを緩やかに、時間をかけて放出しながら、ゆっくりと流れ、変形していく可能性があります。

著者らは、もし地殻が壊れた後に再び大きな結晶へと「修復」されるのであれば、このサイクルが繰り返される可能性があり、それが中性子星から観測されるさまざまな種類の爆発やフレアを説明できる可能性があると示唆しています。

要約すると: 中性子星の地殻は、単に砕け散る脆い岩石ではありません。もしそれが多くの小さな粒で構成されており、ゆっくりと押された場合、それは自己組織化された内部構造のおかげで、壊れることなく曲がったりねじれたりすることができる、超強力で流動的な液体のように振る舞うのです。

技術概要

技術要約：中性子星地殻の可塑性

問題提起
圧力イオン化された原子核からなる結晶性固体である中性子星地殻の弾性的性質は、星震（スタークエイク）、マグネターのバースト、および回転する星における「山」からの連続重力波放射を含む、重要な天体物理学的現象を制御している。完全結晶の線形弾性領域と破壊歪みについては理論的および分子動力学（MD）によって研究されてきたが、極めて大きな歪み下における地殻物質の挙動については、依然として理解が進んでいない。先行研究、特にHorowitzとKadau（HK09）は、高歪みにおける歪み軟化と「アモルファス化」を示唆したが、これらの研究は準静的な条件よりも数桁速い歪み速度を用いており、流体力学的な衝撃波効果を導入した可能性がある。さらに、地殻が粒（グレイン）として形成されるか、初期の破壊後に欠陥ネットワークを発達させる場合に、多結晶構造における弾性から塑性流動への遷移についても、収束した遅い歪み速度を用いて体系的に調査されてこなかった。

手法
著者らは、GPUスーパーコンピュータ上でLAMMPSコードを用いた第一原理分子動力学（MD）シミュレーションを用い、中性子星地殻の変形をモデル化している。

• 物理モデル: 地殻は、電荷 $eZ$と質量$A$ を持つ核の古典的なプラズマとしてモデル化され、相対論的電子縮退を考慮するためにユカワ遮蔽を伴う反発的なクーロンポテンシャルによって相互作用する。電子は明示的にはシミュレートされないが、トーマス・フェルミ近似を介して遮蔽を提供する。
• シミュレーションパラメータ: 周期境界条件を持つ立方体体積を使用する。シミュレーションは、$0.035 T_{melt}$ から $0.30 T_{melt}$（ここで $T_{melt}$ は融解温度）の温度範囲で行われる。
• 歪み速度: 本研究の重要な手法的な進展は、先行研究よりも数桁遅い歪み速度（$\dot{\epsilon}$）を使用している点であり、その範囲は $2 \times 10^{-7} \omega_p$ から $5.1 \times 10^{-4} \omega_p$ である（ここで $\omega_p$ はイオン・プラズマ振動数）。これにより、変形が準静的（$v_{box} \ll v_{el}$）となり、衝撃波を誘発することなく、弾性情報がシミュレーションボックス全体に伝播することを保証している。
• 構成: 完璧な体心立方（bcc）単結晶（$N \approx 128,000$）と、変動する数のグレイン（12個から96個のグレイン）を含む多結晶（$N$ は最大 $960,000$）を比較する。
• 解析: 応力テンソルは対相互作用力から計算される。欠陥の進化は、秩序だったバルク結晶と無秩序な欠陥／粒界を区別する局所結合秩序パラメータ $Q_6$ を用いて可視化される。

主な結果

1. 多結晶における収束した塑性流動: 十分に遅い歪み速度（$\dot{\epsilon} \le 2 \times 10^{-5} \omega_p$）でシミュレートされた多結晶では、材料は線形弾性から、せん断歪み $\epsilon \approx 0.05$ における完全塑性流動へと堅牢に遷移する。この降伏点を超えると、応力はさらなる歪みに関わらず一定（$\sigma_{yield} \approx 2 \times 10^{-3} n_i (eZ)^2/a_i$）に保たれ、非常に大きな歪み（$\epsilon = 0.60$）まで続く。この「弾性－完全塑性」挙動は、高速なシミュレーション（HK09）で報告された歪み軟化やアモルファス化とは対照的である。
2. 単結晶の挙動: 完璧な単結晶は、破壊歪み $\epsilon \approx 0.11$ まで線形弾性を示す。破壊時には、より単純な欠陥ネットワークのために確率性が高くなるものの、同様に塑性流動へと遷移する。破壊応力は歪み速度に敏感であり、より速い速度では、欠陥核生成の熱活性化のための時間が不足するため、応力のオーバーシュートが生じる。
3. 破壊後の塑性の普遍性: 本研究は、破壊後の塑性流動が初期の結晶構造（グレインサイズや数）に依存しないことを発見した。12個のグレインから始めても96個のグレインから始めても、シミュレーションは同じ降伏応力と塑性流動領域へと収束する。結晶は、課された歪み速度を補償するのに十分な密度の欠陥ネットワークを自己整合的に核生成し、初期の欠陥構成を事実上「忘れる」。
4. 歪み速度および温度依存性:
   • 歪み速度: 高速で収束していない速度では、応力のオーバーシュートが発生し、見かけのせん断弾性率は速度依存性を示す。収束した遅い速度では、降伏応力は歪み速度に依存しない。
   • 温度: 降歴応力は、シミュレートされた範囲において緩やかな温度依存性を示す（$\sigma_y \approx 1.5 - 2 \times 10^{-3} n_i (eZ)^2/a_i$）。有効せん断弾性率は、ほぼゼロ温度極限に収束している。
5. メカニズム: 塑性流動は転位運動と粒界滑りによって支配される。著者らは、この完全塑性を、欠陥密度がオーラン（Orowan）の式（$\dot{\epsilon} \propto b n v$）と、応力に比例する粘性抵抗とのバランスをとるように調整される定常状態であると解釈している。

意義および主張
本論文は、高速なシミュレーションの時間スケールの制限により未解決であった、中性子星地殻における収束した定常塑性流動の新しい領域を特定したと主張している。主な意義は、地殻物質が一度降伏すると、硬化することなく塑性的に流動するという発見にあり、この挙動は初期のグレイン構造に関わらず普遍的であることにある。

著者らは、これが天体物理学的モデリングに示唆を与えると考えている：

• マグネター活動: マグネターのバーストやフレアの多様性は、地殻内のグレインサイズの分布に関連している可能性がある。もし応力勾配がグレインサイズよりも小さい場合、地殻は突然降伏し（単結晶のように）、過渡現象を引き起こす可能性がある。もし勾配がより大きい場合、地殻は突発的な破滅的エネルギー放出を伴わずに塑性的に流動する。
• 循環的な破壊と修復: もし破壊された地殻が再アニールして大きな結晶を形成できるのであれば、地殻は弾性エネルギーを蓄積し、破壊し、塑性的に流動し、そして修復するというサイクルを繰り返すことができる。このサイクリングが、繰り返される星震やマグネターのバーストを駆動する可能性がある。

本研究は、現実的なレオロジーを備えた将来のメゾスケールモデリングのための基礎を提供し、純粋に弾性的な地殻、あるいはアモルファスな破壊という仮定を超え、欠陥ダイナミクス（転位）の詳細な特性評価の必要性を強調している。