[论文解读] 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 编写的用户友好型计算工具，利用二阶弹性刚度张量分析二维和三维材料的各向异性弹性特性。该工具通过 Voigt、Reuss 和 Hill 平均方法计算关键力学性能参数，如体变模量、剪切模量、杨氏模量、泊松比及各向异性指数，并利用集成的绘图库可视化声速和弹性各向异性。该工具可识别异常行为，如负泊松比和负线性压缩性，从而支持纳米器件设计的预测性分析。

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.

研究动机与目标

• 通过实现高通量计算分析，解决晶体材料中实验弹性常数稀缺的问题。
• 提供一个统一平台，用于分析多种二维和三维材料的力学及各向异性弹性特性。
• 检测对先进纳米器件至关重要的异常力学行为，如负泊松比和负线性压缩性。
• 通过弹性张量的特征值分析，支持可靠的力学稳定性评估。
• 与主流 DFT 代码及可视化工具集成，实现结果的即时后处理。

提出的方法

• 以通过 VASP、WIEN2k、Quantum Espresso 及外部工具如 ElaStic 和 IRelast 计算的二阶弹性刚度张量作为输入。
• 应用 Voigt、Reuss 和 Hill 平均方法，计算包括体变模量、剪切模量和杨氏模量在内的有效力学性能。
• 计算各向异性指数，如通用各向异性（A_U）、对数欧几里得（A_L）、Chung-Buessem（A_CB）和 Reuss（A_R），以量化弹性各向异性。
• 通过矩阵对角化求解 Christoffel 方程，计算声波相速度和群速度，包括能量传播方向角。
• 与 G NUPLOT、X MGRACE、view3dscene 和 Plotly 集成，实现实时二维和三维弹性特性可视化。
• 支持 IRelast、IRelax2D、ElaStic 和 AELAS 输出格式的输入，并可访问来自 Materials Project 的超过 13,000 个弹性常数的离线数据库。

实验结果

研究问题

• RQ1 如何系统地从二阶弹性刚度张量计算并可视化二维和三维材料的各向异性弹性特性？
• RQ2 基于弹性张量特征值的二维和三维材料的力学稳定性准则是什么？
• RQ3 哪些各向异性指数最可靠地识别出表现出异常力学行为（如负泊松比或负线性压缩性）的材料？
• RQ4 在具有复杂晶系对称性的材料中，使用 Christoffel 方程预测声速各向异性的准确性如何？
• RQ5 是否可以构建一个统一的计算工具，高效处理并可视化具有不同对称性的多种材料的弹性特性？

主要发现

• ElATools 成功识别出 ZnAu₂(CN)₄ 表现出无限各向异性（A_M = ∞）和负线性压缩性，表明其具有异常的力学响应。
• 对于 CrB₂，计算得到的 A_M = 2.227 表明中等程度的各向异性，体变模量范围为 53.51 GPa（最小值）至 119.22 GPa（最大值），剪切模量范围为 125.04 GPa（最小值）至 197.40 GPa（最大值）。
• 在 δ-黑磷中，泊松比最低达 -0.267，杨氏模量从 62.277 N/m（最小值）至 142.877 N/m（最大值），表明具有强烈的各向异性。
• Pd₂O₆Se₂ 单层表现出高达 1.315（最大值）的泊松比和较大的 A_R = 4.257，表明存在显著的弹性各向异性。
• 对于 GaAs，横声模（SS 模）的相速度各向异性（A_M）达到 1.35，群速度各向异性最高达 1.41，证实声波传播具有显著的方向依赖性。
• 该工具正确识别出具有负线性压缩性的材料（如 ZnAu₂(CN)₄，-51.689 TPa⁻¹）和负泊松比的材料（如 δ-P，-0.267），验证了其对异常力学行为的预测能力。