[論文レビュー] Neutron Star Crust and Molecular Dynamics Simulation

C. J. Horowitz, Joe Hughto+2 more|arXiv (Cornell University)|Sep 23, 2011

Pulsars and Gravitational Waves Research参考文献 10被引用数 4

ひとこと要約

この論文では、古典的分子動力学シミュレーションを用いて中性子星の外殻を研究し、急速な拡散と強いクーロン相互作用のおかげで、ほぼ完璧な結晶格子が形成されることを明らかにした。主な発見には、高い熱伝導性および電気伝導性、極めて高い破壊ひずみを示す極めて高い機械的強度、そして奇妙な核スパゲッティ相（pasta phases）の存在が含まれ、これらは中性子星の冷却、重力波放出、r-モードの安定性に影響を与える。

ABSTRACT

In this book chapter we review plasma crystals in the laboratory, in the interior of white dwarf stars, and in the crust of neutron stars. We describe a molecular dynamics formalism and show results for many neutron star crust properties including phase separation upon freezing, diffusion, breaking strain, shear viscosity and dynamics response of nuclear pasta. We end with a summary and discuss open questions and challenges for the future.

研究の動機と目的

• 極度の密度および高温下における中性子星外殻の凍結過程と構造的性質を理解すること。
• 分子動力学シミュレーションを用いて、外殻の機械的強度および熱輸送特性を特定すること。
• 核スパゲッティ相が中性子星の力学的挙動および観測可能な物理量に果たす役割を調査すること。
• 外殻の組成および構造が重力波放出およびr-モード振動に与える影響を調査すること。
• 中性子星における外殻の組成、相転移、散逸メカニズムに関する未解決の問題を解明すること。

提案手法

• 中性子星外殻をイオンと縮重電子のクーロンプラズマとしてモデル化するために、古典的分子動力学（MD）シミュレーションが用いられた。
• 相転移および構造的進化を捉えるために、長距離作用するクーロン相互作用と短距離作用する核力が組み込まれた。
• 集団励起（例：プラズマ振動）を研究するために、時間平均密度-密度相関関数から動的構造因子 $ S(q,\omega) $ が計算された。
• シミュレートされた結晶格子における機械的応力応答から、破壊ひずみおよびせん断率が導出された。
• イオンの分布と移動度の追跡により、結晶相および液体相における化学的分離および拡散速度が分析された。
• 外殻およびスパゲッティ相における空間的および時間的相関を特徴付けるために、静的および動的構造因子が計算された。

実験結果

リサーチクエスチョン

• RQ1 極度の密度および高温下で中性子星外殻はどのように凍結し、凍結過程における化学的分離を引き起こす要因は何か？
• RQ2 中性子星外殻の機械的強度はどの程度で、破壊ひずみおよびせん断率とどのように関係しているか？
• RQ3 奇妙な核スパゲッティ相（例：棒状、平板状、バブル構造）は、外殻の動的応答および輸送特性にどのように影響を与えるか？
• RQ4 非球形のスパゲッティ形状の存在をどのように観測可能な物理量で検出できるか？
• RQ5 熱伝導度や粘性といった外殻の性質は、中性子星の冷却および重力波放出にどのように影響を与えるか？

主な発見

• 急速な拡散のおかげで、欠际が修復され、非晶質固体の形成を防ぐため、中性子星外殻はほぼ完璧な結晶となる。
• 外殻は高い熱伝導性および電気伝導性を示し、降着が停止した後のX線観測で観測される急速な冷却と整合的である。
• 長距離作用するクーロン相互作用のおかげで、外殻の破壊ひずみは極めて高く、極めて高い機械的強度を示している。
• クーロン反発と核力の引力の競合によって生じる核スパゲッティ相（棒状、平板状、バブル構造など）が観測され、$ q = 0.05 $ fm⁻¹ において $ \omega \approx 0.8 $ MeV で明確なプラズマ振動が観察された。
• 動的構造因子 $ S(q,\omega) $ は、高波数 $ q $ でピークが広がっており、クーロンおよび核励起の間の結合が強まっていることを示している。
• 拡散および結晶化のダイナミクスが規則的で秩序ある格子構造を強く支持するため、非晶質外殻の存在はシミュレーションで除外された。