这里会显示出您选择的修订版和当前版本之间的差别。
两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版上一修订版两侧同时换到之后的修订记录 | ||
adf:tpa2019 [2020/03/11 11:05] – [结果查看] liu.jun | adf:tpa2019 [2022/10/31 19:16] – liu.jun | ||
---|---|---|---|
行 1: | 行 1: | ||
====== 计算动态超极化率、双光子吸收截面TPA Cross Section====== | ====== 计算动态超极化率、双光子吸收截面TPA Cross Section====== | ||
+ | |||
+ | 本例以MoS2为例。软件版本为AMS2019.3及之前的版本,AMS2020以后的版本请参考链接:[[adf: | ||
本文参考文献资料: | 本文参考文献资料: | ||
行 9: | 行 11: | ||
{{ : | {{ : | ||
- | * 精细结构常数α=e$^2$/4πε$_0$cħ,其中e电子电荷,ε$_0$真空介电常数,c真空光速,ħ约化普朗克常数。 | + | * 精细结构常数α=e<sup>2</ |
* 整数N取决于实验设置,很多文章中都取N=4 | * 整数N取决于实验设置,很多文章中都取N=4 | ||
* γ约化形式的二阶超极化率γ(-ω, | * γ约化形式的二阶超极化率γ(-ω, | ||
* ω频率 | * ω频率 | ||
* Im[]表示取虚部 | * Im[]表示取虚部 | ||
- | * σ$^{TPA}$单位为Göppert-Mayer (1 GM = 10$^{−50}$ cm$^4$ s photon$^{−1}$) | + | * σ<sup>TPA</ |
原子单位下: | 原子单位下: | ||
* α=1/ | * α=1/ | ||
行 21: | 行 23: | ||
* ADF得到的γ< | * ADF得到的γ< | ||
- | 如果我们得到原子单位的γ< | + | 如果我们得到原子单位的γ< |
我们仅仅选取其中1-NMe2为例计算TPA截面。 | 我们仅仅选取其中1-NMe2为例计算TPA截面。 | ||
行 35: | 行 37: | ||
SAOP是ADF独有的一种模型势,响应性质计算方面,一般优于其他LDA、GGA、metaGGA泛函,因此本文采用SAOP。 | SAOP是ADF独有的一种模型势,响应性质计算方面,一般优于其他LDA、GGA、metaGGA泛函,因此本文采用SAOP。 | ||
====关于基组:==== | ====关于基组:==== | ||
- | 原则上而言,基组越大越精确,但是根据文献中介绍,基组对结果影响不大,因此本文采用了很小的DZP基组。如果体系中包含较重的元素,还是建议为重元素设置TZP,甚至TZ2P、QZ4P基组。设置方法参考:[[adf: | + | 原则上而言,基组越大越精确,但是根据文献中介绍,基组对结果影响不大,因此本文采用了很小的DZP基组。如果体系中包含较重的元素,还是建议为重元素设置TZP,甚至TZ2P、QZ4P基组。设置方法参考:[[adf: |
====数值精度:==== | ====数值精度:==== | ||
这里为了快速计算得到结果,因此使用了默认值Normal,实际上建议选择Good。每提高一个等级,计算量都会数倍增长。 | 这里为了快速计算得到结果,因此使用了默认值Normal,实际上建议选择Good。每提高一个等级,计算量都会数倍增长。 | ||
行 96: | 行 98: | ||
上面数据中,Im[γ< | 上面数据中,Im[γ< | ||
- | 本例中,27个数据之和$γ^{TPA}$=824647744.9551 a.u.,ω = 1.56 eV = 1.56/ | + | 本例中,27个数据之和γ<sup>TPA</ |
- | 分别计算1.46~1.76 eV之间的结果({{: | + | 分别计算1.46~1.76 eV之间的结果,得到σ<sup>TPA</ |
{{ : | {{ : | ||
行 111: | 行 113: | ||
{{ : | {{ : | ||
- | |||
- | {{adf: | ||
=====致谢===== | =====致谢===== | ||
感谢该文献作者Abdou Boucekkine教授,帮助确定公式中Im[γ< | 感谢该文献作者Abdou Boucekkine教授,帮助确定公式中Im[γ< |