[論文レビュー] Molecular Dynamics Simulation of Shear Moduli for Coulomb Crystals

C. J. Horowitz, J. Hughto|ArXiv.org|Dec 15, 2008

Advanced ceramic materials synthesis参考文献 1被引用数 43

ひとこと要約

本稿では、中性子星の地殻におけるクーロン結晶の分子動力学（MD）シミュレーションを提示し、電子スクリーニングを組み込んで、以前のモンテカルロ研究よりも正確にせん断率を計算する。主な結果として、オガタらの結果と比較して、電子スクリーニングの影響により有効せん断率が約10%低下することが判明し、Γ ≥ 175 の範囲で有効せん断率 μ_eff ≈ (0.1106 − 28.7/Γ^1.3)(nZ²e²/a) のフィット式が有効である。

ABSTRACT

Torsional (shear) oscillations of neutron stars may have been observed in quasiperiodic oscillations of Magnetar Giant Flares. The frequencies of these modes depend on the shear modulus of neutron star crust. We calculate the shear modulus of Coulomb crystals from molecular dynamics simulations. We find that electron screening reduces the shear modulus by about 10% compared to previous Ogata et al. results. Our MD simulations can be extended to calculate the effects of impurities and or polycrystalline structures on the shear modulus.

研究の動機と目的

• 電子スクリーニング効果を組み込むことで、中性子星の地殻におけるせん断率の計算精度を向上させること。
• 統計的精度がより高い分子動力学を用いることで、以前のモンテカルロシミュレーションの限界を克服すること。
• 今後の破壊ひずみのシミュレーションの基盤を提供すること。破壊ひずみは、中性子星のマウントの安定性および重力波放射に支配的である。
• トーション振動に関連するクーロン結晶における弾性特性に及ぼす電子スクリーニングの影響を定量化すること。
• 今後の研究において、MDシミュレーションを用いて中性子星の地殻における不純物および多結晶構造の研究を可能にすること。

提案手法

• イオンがスクリーニングされたクーロン力 v(r) = Z²e²/r · e^(-r/λ_e) で相互作用する1成分プラズマ（OCP）の分子動力学シミュレーションを実行する。
• 縮重電子ガスにおける電子媒介スクリーニングをモデル化するため、電子スクリーニング長 λ_e = π^(1/2)/(2e(3π²n_e)^(1/3)) を用いる。
• 体積を一定に保つ6つの変形（D1–D6）を用い、ひずみ u_ik を誘導し、全エネルギーの二階微分 d²V_tot/dε² を通じて弾性定数 b_11 と c_44 を計算する。
• 有効せん断率 μ_eff = (2b_11 + 3c_44)/5 を、角度平均された弾性応答として計算する。
• ポテンシャルに有限のカットオフ距離 R_cut を導入し、R_cut → ∞ への外挿を実施することで、長距離相互作用の影響を評価する。
• 大規模なシミュレーション（N = 3456 および N = 9826 個のイオン）と有限サイズスケーリングを用い、系サイズに起因する系統的誤差を最小限に抑える。

実験結果

リサーチクエスチョン

• RQ1 電子スクリーニングは、中性子星の地殻におけるクーロン結晶のせん断率に、以前のスクリーニングなしの計算と比較してどのように影響するか？
• RQ2 有限サイズの系が、MDシミュレーションにおける計算されたせん断率に及ぼす定量的影響は何か？
• RQ3 b_11 と c_44 は、カットオフ距離 R_cut を大きくするにつれてどのように収束するか？
• RQ4 Coulombパラメータ Γ ≥ 175 の範囲で、有効せん断率 μ_eff を正確に記述できる単純な解析的フィット式が存在するか？
• RQ5 電子スクリーニングは、中性子星の地殻におけるせん断波の位相速度にどのような影響を及ぼすか？

主な発見

• 電子スクリーニングにより、有効せん断率 μ_eff は、オガチらのスクリーニングなしの結果と比較して約10%低下する。
• MDシミュレーションにより、μ_eff ≈ (0.1106 − 28.7/Γ^1.3)(nZ²e²/a) が得られ、Γ ≥ 175 の範囲でフィット誤差は ≤ 2% である。
• N = 9826 個のイオンの場合、有限サイズ効果は約1%と推定され、表2の μ_eff は表1の値と比較して約1%小さい。
• b_11 と c_44 の収束には、R_cut > 12λ_e が必要であり、長距離スクリーニングの重要性が強調される。
• 角度平均せん断率 μ_eff はスクリーニングの影響で約5%低下し、これはトーション振動の周波数をわずかに低減させる。
• 結果は、ゼロ温度への外挿計算と整合しており、電子スクリーニングが地殻を著しく柔ららかくするという事実を確認している。