这里会显示出您选择的修订版和当前版本之间的差别。
两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版 | ||
atk:弹性常数 [2018/05/16 14:54] – [运用DFT计算弹性常数] fermi | atk:弹性常数 [2018/08/22 21:39] (当前版本) – [弹性常数] fermi | ||
---|---|---|---|
行 1: | 行 1: | ||
====== 弹性常数 ====== | ====== 弹性常数 ====== | ||
- | Virtual NanoLab 和 Atomistix ToolKit 可以提供非常简单的方法计算**任意块体的弹性常数**。该方法具有普遍适用性,极易与密度泛函理论(ATK-DFT)、半经验方法(ATK-SE)、经典势(ATK-ForceField)结合使用。在本教程中,您将学习如何借助 | + | QuantumATK |
<WRAP center important 100%> | <WRAP center important 100%> | ||
行 14: | 行 14: | ||
应力和应变张量总是组成对称的 3×3 矩阵,因此它们可以通过 Voigt notation 表示法用 6 个矢量简洁地表达。 | 应力和应变张量总是组成对称的 3×3 矩阵,因此它们可以通过 Voigt notation 表示法用 6 个矢量简洁地表达。 | ||
- | $$\boldsymbol{\sigma} = ( \sigma_{xx}, | + | $$ |
+ | \boldsymbol{\sigma} = ( \sigma_{xx}, | ||
+ | $$ | ||
和 | 和 | ||
- | $$\boldsymbol{\varepsilon} = ( \varepsilon_{xx}, | + | $$ |
+ | \boldsymbol{\varepsilon} = ( \varepsilon_{xx}, | ||
+ | $$ | ||
应力张量和给定应变矢量的线性关系可以表达为: | 应力张量和给定应变矢量的线性关系可以表达为: | ||
- | $$\boldsymbol{\sigma} = \boldsymbol{C} \cdot \boldsymbol{\varepsilon} ,$$ | + | $$ |
+ | \boldsymbol{\sigma} = \boldsymbol{C} \cdot \boldsymbol{\varepsilon} , | ||
+ | $$ | ||
这里的 $\boldsymbol{C}$ 是包含**弹性常数**的 6×6 对称矩阵。 | 这里的 $\boldsymbol{C}$ 是包含**弹性常数**的 6×6 对称矩阵。 | ||
行 28: | 行 34: | ||
随着晶体对称性的增高,矩阵 $\boldsymbol{C}$ 中独立项数目越少。例如,在立方晶体中只有 $C_{11}$、$C_{12}$ 和 $C_{44}$ 这 3 个独立项: | 随着晶体对称性的增高,矩阵 $\boldsymbol{C}$ 中独立项数目越少。例如,在立方晶体中只有 $C_{11}$、$C_{12}$ 和 $C_{44}$ 这 3 个独立项: | ||
- | $$\begin{split}\begin{pmatrix} | + | $$ |
+ | \begin{split}\begin{pmatrix} | ||
\sigma_1 \\ | \sigma_1 \\ | ||
\sigma_2 \\ | \sigma_2 \\ | ||
行 53: | 行 60: | ||
\varepsilon_5 \\ | \varepsilon_5 \\ | ||
\varepsilon_6 \\ | \varepsilon_6 \\ | ||
- | \end{pmatrix}\end{split}$$ | + | \end{pmatrix}\end{split} |
+ | $$ | ||
为了获得弹性常数,需要在模拟的晶胞上沿选定的应变矢量上施加一个小的形变 $\eta$,计算产生的应力张量。通过拟合每个 Vogit 应力和对应的应变矢量组成的 $\sigma_i(\eta)$ 曲线就可得到线性应力分布情况。然后将弹性常数作为计算线性方程组的最小二乘解,同时考虑晶体的对称性。 | 为了获得弹性常数,需要在模拟的晶胞上沿选定的应变矢量上施加一个小的形变 $\eta$,计算产生的应力张量。通过拟合每个 Vogit 应力和对应的应变矢量组成的 $\sigma_i(\eta)$ 曲线就可得到线性应力分布情况。然后将弹性常数作为计算线性方程组的最小二乘解,同时考虑晶体的对称性。 | ||
行 147: | 行 155: | ||
日志文件底部显示有弹性常数矩阵,且与使用了 Stillinger-Weber 势的文献中的结果高度吻合 <color # | 日志文件底部显示有弹性常数矩阵,且与使用了 Stillinger-Weber 势的文献中的结果高度吻合 <color # | ||
+ | $$ | ||
| | ||
| | ||
行 153: | 行 161: | ||
| | ||
\end{align}\end{split} | \end{align}\end{split} | ||
+ | $$ | ||