[論文レビュー] ElATools: A tool for analyzing anisotropic elastic properties of the 2D and 3D materials

Shahram Yalameha, Zahra Nourbakhsh+1 more|arXiv (Cornell University)|May 15, 2021

Chemical and Physical Properties of Materials参考文献 3被引用数 5

ひとこと要約

ElATools は Fortran 90 で記述された使いやすい計算ツールであり、2次元および3次元材料の異方性弾性特性を、2次弾性剛性テンソルを用いて分析する。このツールは Voigt、Reuss、Hill の平均化法を用いて体積率、剪断率、ヤング率、ポisson比、異方性インデックスといった主要な機械的特性を計算し、統合されたプロットライブラリを介して音速および弾性異方性を可視化する。負のポisson比や負の線形圧縮率といった異常な挙動を同定でき、ナノスケールデバイス設計における予測的分析を可能にする。

ABSTRACT

We introduce a computational method and a user-friendly code with a terminal-based graphical user interface (GUI), named ElATools, developed to analyze mechanical and anisotropic elastic properties. ElATools enables facile analysis of the second-order elastic stiffness tensor of two-dimensional (2D) and three-dimensional (3D) crystal systems. It computes and displays the main mechanical properties including the bulk modulus, Young's modulus, shear modulus, hardness, p-wave modulus, universal anisotropy index, Chung-Buessem anisotropy index, log-Euclidean anisotropy parameter, Cauchy pressures, Poisson's ratio, and Pugh's ratio, using three averaging schemes of Voigt, Reuss, and Hill. It includes an online and offline database from the Materials Project with more than 13,000 elastic stiffness constants for 3D materials. The program supports output files of the well-known computational codes IRelast, IRelast2D, ElaStic, and AELAS. Four types of plotting and visualization tools are integrated to conveniently interface with GNUPLOT, XMGRACE, view3dscene, and plotly libraries, offering immediate post-processing of the results. It can efficiently identify anomalous mechanical properties, such as negative linear compressibility, negative Poisson's ratio, and highly-anisotropic elastic modulus in 2D and 3D materials, which are central properties to design and develop high-performance nanoscale electromechanical devices. Six case studies on selected material systems, namely, ZnAu$_2$ (CN)$_4$, CrB$_2$, $δ$-phosphorene, Pd$_2$O$_6$Se$_2$ monolayer, and GaAs, and a hypothetical set of systems with cubic symmetry are presented to demonstrate the descriptive and predictive capabilities of ElATools.

研究の動機と目的

• 結晶性材料における実験的弾性定数の不足を補うために、高スルーレットな計算分析を可能にする。
• 多様な 2次元および3次元材料にわたる機械的および異方性弾性特性の統合的プラットフォームを提供する。
• 先進的なナノスケールデバイスに不可欠な、負のポisson比や負の線形圧縮率といった異常な機械的挙動を同定する。
• 弾性テンソルの固有値解析を通じて、信頼性の高い機械的安定性評価を支援する。
• 主要な DFT コードおよび可視化ツールと連携し、即時の結果後処理を可能にする。

提案手法

• DFT を用いたコード（VASP、WIEN2k、Quantum Espresso および ElaStic や IRelast などの外部パッケージを含む）によって計算された 2次弾性剛性テンソルを入力として使用する。
• Voigt、Reuss、Hill の平均化法を適用し、体積率、剪断率、ヤング率などの有効な機械的特性を計算する。
• 普遍的異方性（A_U）、対数ユークリッド（A_L）、Chung-Buessem（A_CB）、Reuss（A_R）といった異方性インデックスを計算し、弾性異方性を定量化する。
• 行列の対角化によりクリストッフェル方程式を解き、音響位相速度および群速度（パワー流れの角度を含む）を計算する。
• G NUPLOT、X MGRACE、view3dscene、Plotly と連携し、弾性特性のリアルタイムな 2次元および3次元可視化を実現する。
• IRelast、IRelast2D、ElaStic、AELAS の出力形式からの入力をサポートし、Materials Project から得た 13,000 以上の弾性定数のオフラインデータベースにアクセスできる。

実験結果

リサーチクエスチョン

• RQ1 2次弾性剛性テンソルから出発して、2次元および3次元材料の異方性弾性特性を体系的に計算・可視化する方法は何か？
• RQ2 2次元および3次元材料の機械的安定性基準は、弾性テンソルの固有値に基づいてどのように定義されるか？
• RQ3 負のポisson比や負の線形圧縮率といった異常な機械的挙動を示す材料を最も信頼性高く同定するための異方性インデックスはどれか？
• RQ4 複雑な結晶対称性を示す材料において、クリストッフェル方程式を用いた音速異方性の予測精度はどの程度か？
• RQ5 多様な結晶対称性を有する材料を効率的に処理・可視化できる統合型計算ツールは、実際に実用的か？

主な発見

• ElATools は、ZnAu₂(CN)₄ が無限大の異方性（A_M = ∞）および負の線形圧縮率を示すことを特定し、異常な機械的応答を示していることを明らかにした。
• CrB₂ については、計算された A_M = 2.227 が中程度の異方性を示しており、体積率は最小 53.51 GPa から最大 119.22 GPa、剪断率は最小 125.04 GPa から最大 197.40 GPa の範囲に分布している。
• δ-リン膜では、ポisson比が -0.267（最小）に達し、ヤング率は最小 62.277 N/m から最大 142.877 N/m まで変動し、強い異方性が示された。
• Pd₂O₆Se₂ モノレイヤーでは、最大のポisson比 1.315 および大きな A_R = 4.257 を示し、強い弾性異方性を示した。
• GaAs では、SS モードにおける位相速度の異方性（A_M）が 1.35 に達し、群速度の異方性は最大 1.41 に達しており、音波伝播における顕著な方向依存性が確認された。
• このツールは、ZnAu₂(CN)₄ が -51.689 TPa⁻¹ の負の線形圧縮率を示すなど、負の線形圧縮率を示す材料、および δ-P が -0.267 の負のポisson比を示すなど、負のポisson比を示す材料を正しく同定し、異常な機械的挙動に対する予測能力を確認した。