[論文レビュー] Elasticity in crystals with high density of local defects : insights from ultra-soft colloids

Saswati Ganguly, Gaurav P. Shrivastav+6 more|arXiv (Cornell University)|Aug 13, 2021

Material Dynamics and Properties参考文献 42被引用数 7

ひとこと要約

本稿では、高密度の点欠陥を有する結晶の弾性を記述する初等的理論枠組みを、超軟らかなコロイドをモデル系として開発する。粗粒化場を用いて欠陥拡散を連続体弾性に組み込むことで、理論は測定可能な弾性応答の変化を予測し、分子動力学シミュレーションによって検証される。その結果、欠陥は従来の弾性理論をはるかに超えて機械的挙動を顕著に変化させることを示している。

ABSTRACT

In complex crystals close to melting or at finite temperatures, different types of defects are ubiquitous and their role becomes relevant in the mechanical response of these solids. Conventional elasticity theory fails to provide a microscopic basis to include and account for the motion of point-defects in an otherwise ordered crystalline structure. We study the elastic properties of a point-defect rich crystal within a first-principles theoretical framework derived from microscopic equations of motion. This framework allows us to make specific predictions pertaining to the mechanical properties which we can validate through deformation experiments performed in Molecular Dynamics simulations.

研究の動機と目的

• 点欠陥の密度が非常に高い系において、従来の弾性理論が失敗する状況での結晶の微視的弾性理論の構築。
• 点欠陥（空孔および間隙原子）の動的挙動を、弾性応答の連続体記述に組み込むこと。
• 非理想結晶における微視的密度相関と巨視的弾性特性との間の関係を確立すること。
• 高欠陥濃度を有するモデル系において、原子的分子動力学シミュレーションと照合して理論的予測を検証すること。
• 欠陥拡散が流体力学的モードとして果たす役割と、弾性定数および波動伝播に与える影響を明確にすること。

提案手法

• 微視的運動方程式に基づく形式的枠組みを導出し、連続体弾性理論を欠陥密度を動的場として含める形に拡張する。
• 変位 u(r,t) および欠陥密度 cd(r,t) の粗粒化場を導入し、密度ゆらぎを δn(r,t) = −n₀∇·u(r,t) − δcd(r,t) で与える。
• 弾性応答係数を微視的2点相関関数（密度-密度相関関数および直接相関関数）の形で表現する。
• 平衡統計力学を用いて、弾性率を未変形結晶内の構造因子および対相関関数と関連付ける。
• 超軟らかなコロイドからなるクラスタ結晶で、高欠陥濃度を示すモデル系に対して分子動力学シミュレーションを実施する。
• シミュレーションでは制御された変形を加え、応力-ひずみ応答を測定することで、理論的予測の弾性定数および欠陥媒介による軟化の検証を行う。

実験結果

リサーチクエスチョン

• RQ1 高密度の点欠陥が、標準的弾性理論をはるかに超えて結晶固体の機械的応答にどのように影響を与えるか？
• RQ2 欠陥拡散が、非理想結晶の弾性応答において、どの程度流体力学的モードとして寄与するか？
• RQ3 経験的フィッティングを一切行わずに、微視的相関関数から欠陥豊富な結晶の弾性定数を予測できるか？
• RQ4 欠陥ゆらぎは、弾性波の伝播および長波長域の機械的挙動にどのように影響を与えるか？
• RQ5 欠陥濃度が、モデルのソフトコロイド結晶における有効せん断率および体積弾性率に定量的にどのように影響を与えるか？

主な発見

• 理論は、欠陥拡散と弾性ひずみ場の結合によって、欠陥豊富な結晶で弾性率が軟化することを成功裏に予測した。
• 欠陥密度ゆらぎが密度変化に顕著に寄与し、密度変化が単に変位場に起因するという従来の仮定が成立しないことを示した。
• 分子動力学シミュレーションにより、有効せん断率が欠陥濃度の増加に伴い減少することが確認され、理論的予測と定量的に一致した。
• この枠組みは、欠陥拡散が長波長域の弾性応答に影響を与える独立した流体力学的モードであることを特定した。これは、物性的流体力学と整合的である。
• 弾性応答が微視的直接相関関数に依存することを示した。この関数は結晶内の欠陥誘発構造変化を記述している。
• 理論により、硬球系（参考文献9）におけるような、硬さの極限でのキャンセル効果が、ソフトポテンシャル系では生じないことが明らかになった。このため、超軟らかなコロイドでは欠陥効果がより顕著に現れる。