[论文解读] ElasTool: An automated toolkit for elastic constants calculation

Zhongli Liu|arXiv (Cornell University)|Feb 16, 2020

Machine Learning in Materials Science参考文献 29被引用 10

一句话总结

ElasTool 是一个用于通过 DFT 计算 2D 和 3D 晶体中二阶弹性常数（SOECs）的自动化 Python 工具包，采用三种应变矩阵集——OHESS、ASESS 和 ULICS——基于应力-应变法。它能够高效、准确地计算零温及高温下的弹性性质，包括弹性模量、各向异性、德拜温度和稳定性分析，显著提升了高通量材料设计与弹性研究的效率。

ABSTRACT

We here present, the ElasTool package, an automated toolkit for calculating the second-order elastic constants (SOECs) of any two- (2D) and three-dimensional (3D) crystal systems. It can utilize three kinds of strain-matrix sets, the high-efficiency strain-matrix sets (OHESS), the universal linear-independent coupling strains (ULICS) and the all-single-element strain-matrix sets (ASESS) to calculate the SOECs automatically. In an automatic manner, ElasTool can deal with both zero- and high-temperature elastic constants. The theoretical background and computational method of elastic constants, the package structure, the installation and run, the input/output files, the controlling parameters, and two representative examples of ElasTool are described detailedly. ElasTool is useful for either the exploration of materials' elastic properties or high-throughput new materials design. ElasTool is also available at our website: www.matdesign.cn.

研究动机与目标

• 开发一个自动化、高效的工具包，用于计算 2D 和 3D 晶体体系中的二阶弹性常数（SOECs）。
• 通过采用优化的应变矩阵集（OHESS、ASESS、ULICS），降低基于 DFT 的弹性常数计算的计算成本。
• 支持零温静态应力与高温动态应力（通过分子动力学）的计算。
• 提供全面的输出结果，包括弹性模量、声速、德拜温度、各向异性和机械稳定性分析。
• 通过简化弹性性质计算流程，实现高通量材料筛选与理论设计。

提出的方法

• 基于胡克定律的应力-应变法：σ_i = Σ C_ij ε_j，其中应力通过 DFT 计算获得，对应多种应变构型。
• 集成三种应变矩阵集：OHESS（优化高效率应变集）、ASESS（单元素应变矩阵）和 ULICS（通用线性无关耦合应变）。
• 自动化完整工作流程：输入解析、通过用户自定义命令（如 VASP）执行 DFT 应力计算、应力-应变拟合及弹性常数提取。
• 使用 DFT 计算的应力张量（来自 VASP 或其他代码），并通过最小二乘法拟合，从线性应力-应变响应中确定弹性常数。
• 支持常规晶胞与原胞输入，可选择静态或动态应力平均（例如来自 MD 模拟）。
• 输出结果包括 C_ij、体积累积模量（Voigt、Reuss、Voigt-Reuss-Hill 平均值）、剪切模量、杨氏模量、泊松比、声速、德拜温度、各向异性及 Born 弹性稳定性判据分析。

实验结果

研究问题

• RQ1 如何通过使用优化的应变矩阵集，提升基于 DFT 的弹性常数计算的计算效率？
• RQ2 能否开发一个自动化工作流，以支持在多种 2D 和 3D 晶体体系中同时处理零温与高温下的弹性性质？
• RQ3 与 ULICS 和 ASESS 相比，OHESS 应变矩阵集在最小化 DFT 计算时间的同时，其性能与精度如何？
• RQ4 ElasTool 能在多大程度上基于 DFT 应力数据准确预测机械稳定性、弹性各向异性和德拜温度？
• RQ5 ElasTool 在多大程度上通过减少弹性性质计算中的手动干预，支持高通量材料筛选？

主要发现

• ElasTool 能够成功使用三种应变矩阵集（OHESS、ASESS、ULICS）计算任意 2D 或 3D 晶体体系的二阶弹性常数（SOECs）。
• 对于金刚石（立方晶系），计算得到的 SOECs 为：C11 = 1055.04 GPa，C12 = 136.56 GPa，C44 = 567.76 GPa，平均压力为 -0.08 GPa。
• 体积累积模量的 Voigt 平均值（B_V）和 Reuss 平均值（B_R）均为 442.72 GPa，Voigt-Reuss-Hill 平均值（B_VRH）也为 442.72 GPa。
• 剪切模量值为：G_V = 524.36 GPa，G_R = 518.73 GPa，G_VRH = 521.54 GPa，杨氏模量 E = 1123.47 GPa，泊松比 ν = 0.0771。
• 声速为：V_S = 12.20 km/s，V_B = 11.24 km/s，V_P = 18.03 km/s，V_M = 13.31 km/s，德拜温度 T_D = 1761.03 K。
• 通过 Born 弹性稳定性准则确认该结构具有机械稳定性，弹性各向异性分别为 A_U = 0.0542 和 A_C = 0.0054。