[论文解读] Investigating elastic constants across diverse strain-matrix sets

Zhongli Liu, Wei, Ya-Dong+6 more|arXiv (Cornell University)|Jan 31, 2020

Microstructure and mechanical properties被引用 11

一句话总结

本文提出了优化高效应变-矩阵集（OHESS），一种基于密度泛函理论中应力-应变方法计算二阶弹性常数的新方法。通过在形变过程中保持晶体对称性，OHESS在高对称性三维晶体中实现了卓越的计算效率，相较于ULICS和ASESS在速度上表现更优，且不损失准确性。

ABSTRACT

Elastic constants and mechanical properties play a pivotal role across multiple disciplines and engineering applications. We introduced the optimized high-efficient strain-matrix set (OHESS) that determines the second-order elastic constants of materials using the stress-strain method. Herein, we systematically investigate the computational efficiency of OHESS across a broad range of crystal systems and compare it with other notable stress-strain approaches, such as the single-element strain-matrix sets and the universal linear-independent coupling strains. Notably, our data affirm the superior efficacy of OHESS among the strain-matrix sets under consideration. We believe OHESS will markedly improve computational efficiency in determining the elastic constants and mechanical properties, becoming an indispensable tool for material research, design, and high-throughput screening.

研究动机与目标

• 开发一种更高效的应变-矩阵集，用于通过应力-应变方法计算二阶弹性常数。
• 解决现有方法（如ULICS）中因对称性降低的应变集导致的高计算成本问题。
• 在多种三维和二维晶体体系中，评估并比较OHESS、ULICS与ASESS的精度与效率。
• 通过弹性常数实现更快的高通量材料筛选与稳定性评估。
• 优化应变集，使形变过程中保持对称性，从而降低计算开销。

提出的方法

• 提出OHESS作为应变-矩阵集，可在形变过程中保持晶体的固有对称性。
• 基于胡克定律的应力-应变关系，从施加应变下的应力张量计算弹性常数。
• 优化应变分量的数量与类型，以最小化计算成本，同时确保线性无关性与弹性常数的完整覆盖。
• 采用密度泛函理论（DFT）计算每个应变集下的应力张量，并通过足够的k点采样保证精度。
• 从计算时间与对称性保持两个方面，将OHESS与ULICS（通用线性无关耦合应变）和ASESS（替代应变集策略）进行比较。
• 在三维体系（立方至三斜）和二维体系（方形、六方、矩形晶格）中评估性能。

实验结果

研究问题

• RQ1 在各种晶体体系中，OHESS相较于ULICS和ASESS的计算效率如何？
• RQ2 在弹性常数计算中，应变应用时的对称性保持在多大程度上降低了计算成本？
• RQ3 在哪些晶体体系中，OHESS相较于现有应变-矩阵集展现出最大的效率优势？
• RQ4 OHESS中的应变集数量与ASESS和ULICS相比如何？这对性能有何影响？
• RQ5 OHESS是否在显著降低计算时间的同时，仍保持与ULICS和ASESS相当的精度？

主要发现

• 在大多数三维晶体体系中，OHESS的计算效率超过ULICS和ASESS的两倍以上，尤其在立方和六方等高对称性结构中表现更优。
• 在石墨烯和MoS2等二维体系中，OHESS的性能与ULICS和ASESS相当，但因应变集更少而展现出更优的效率。
• ULICS效率较低的原因在于其对称性降低的应变耦合，导致应力张量计算成本增加。
• ASESS效率较低的原因在于其在三维体系中固定使用六组应变、二维体系中固定使用三组应变，与晶体对称性无关。
• OHESS在显著降低计算时间的同时保持了高精度，使其成为高通量材料筛选与数据库构建的理想选择。
• 在低对称性体系（如单斜和三斜）中，由于形变后对称性差异减小，该方法的效率优势有所减弱。