[论文解读] Extended scheme for the projection of material tensors of arbitrary symmetry onto a higher symmetry tensor

A. Miguel|arXiv (Cornell University)|Aug 6, 2014

Elasticity and Material Modeling被引用 3

一句话总结

本文通过引入旋转优化，扩展了Moakher和Norris的张量投影方法，以找到与低对称性材料张量（如特殊准随机结构（SQS）超胞中的三斜弹性张量）最接近的对称张量（例如立方对称张量）。通过在投影前对张量进行旋转以最小化欧氏距离，该方法相比标准投影方法将误差降低了最多20%，从而更准确地提取出合金的宏观弹性常数。

ABSTRACT

I propose a straightforward generalization of the projection scheme for elastic tensors introduced by Moakher and Norris [J. Elasticity 85, 215 (2006)] that takes into account also rotations. The "closest" tensor of any desired symmetry to the original tensor of lower symmetry is "closer" in this generalized scheme. The method has an important application in the context of the special quasirandom structure (SQS) method for the computational modeling of alloys, whereby the supercell's symmetry, and therefore that of the tensors representing its properties, is reduced with respect to the material's underlying symmetry. The approach allows to extract the tensor components most representative of the macroscopic symmetry of the material. Although the approach is general, in the present case I apply it to the elastic tensor and give numerical examples. Simple approximate analytical expressions for cubic materials are also provided.

研究动机与目标

• 解决现有张量投影方法在将低对称性张量投影到高对称形式时忽略旋转自由度的局限性。
• 改进从具有降低对称性的超胞（如通过特殊准随机结构（SQS）方法生成的超胞）进行第一原理计算时提取宏观弹性性质的准确性。
• 为任何因无序或有限尺寸效应导致对称性轻微破缺的材料张量（例如弹性、压电张量）提供一个通用框架。
• 通过旋转引导的投影最小化无序超胞张量与其对称对应张量之间的距离，实现更可靠的本征材料性质提取。

提出的方法

• 该方法通过在标准投影算子应用前引入旋转步骤，推广了Moakher和Norris的投影方案。
• 将问题表述为在所有可能的旋转角度下，最小化原始张量的旋转版本与目标对称张量（例如立方对称）之间的欧氏距离。
• 使用由三个欧拉角（tx, ty, tz）参数化的旋转矩阵，对三斜弹性张量施加旋转，这些角度通过数值方法进行优化。
• 旋转后的张量以Voigt记法表示，并应用适当的归一化因子（例如C12为√2，C44为2），以保持张量范数不变。
• 通过Voigt记法中的三个基向量构建立方对称弹性张量的空间投影算子。
• 通过Mathematica中的NMinimize进行数值优化，以找到使旋转张量与其投影立方形式之间距离最小的旋转角度。

实验结果

研究问题

• RQ1 通过在投影过程中考虑旋转自由度，能否提升张量投影的准确性？
• RQ2 三斜弹性张量的最优旋转角度是什么，可使其与最近的立方张量之间的距离最小？
• RQ3 与标准投影相比，在无序超胞的弹性张量投影中引入旋转后，欧氏距离能减少多少？
• RQ4 旋转优化在从基于SQS的第一原理计算中提取宏观弹性常数方面，改善程度如何？
• RQ5 在新的旋转优化框架下，能否为立方材料推导出解析近似表达式？

主要发现

• 旋转优化投影将三斜弹性张量与其立方投影之间的欧氏距离从标准方法的91.0 GPa降低至69.997 GPa，降幅约为23%。
• 最优旋转角度约为tx = -0.0329 rad，ty = -0.0319 rad，tz = 0.1111 rad，表明微小但显著的旋转可改善对齐效果。
• 旋转并投影后，所得立方弹性常数为：C11 = C22 = C33 = 436.8 GPa，C12 = C13 = C23 = 228.3 GPa，C44 = C55 = C66 = 377.5 GPa。
• 通过Voigt记法中的适当归一化，该方法保持了张量范数，确保了距离度量的一致性。
• 该方法具有通用性，适用于任意对称性的材料张量，不限于弹性性质。
• 本文在新框架下为立方材料提供了简洁的解析表达式，可高效估算来自无序体系的对称张量分量。