这里会显示出您选择的修订版和当前版本之间的差别。
两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版 | ||
adf:lihhyper [2020/12/01 14:28] – [模型] liu.jun | adf:lihhyper [2022/01/20 21:18] (当前版本) – [二阶超极化率γ] liu.jun | ||
---|---|---|---|
行 7: | 行 7: | ||
本文旨在比较各种计算方法。作为一个最简单的入手案例,我们选择LiH分子。按照作者选取的泛函(LDA)、基组(SZ)、计算方法(Damped | 本文旨在比较各种计算方法。作为一个最简单的入手案例,我们选择LiH分子。按照作者选取的泛函(LDA)、基组(SZ)、计算方法(Damped | ||
- | Cubic Response Theory,在GUI中对应Properties - Hyperolarizability - Calculate - γ(damped 2nd hyperpol),选择该选项会同时计算α、β、γ),频率($ω_1$、$ω_2$、$ω_3$、liftime均为0,因为计算的是静态超极化率)。 | + | Cubic Response Theory,在GUI中对应Properties - Hyperolarizability - Calculate - γ(damped 2nd hyperpol),选择该选项会同时计算α、β、γ),频率(ω< |
文中的计算结果为: | 文中的计算结果为: | ||
行 23: | 行 23: | ||
{{ : | {{ : | ||
- | 由于计算的是静态超极化率,因此lifetime、$ω_1$、$ω_2$、$ω_3$均设置为0。如果计算动态极化、超极化,设置相应的频率、寿命即可。具体可以参考: | + | 由于计算的是静态超极化率,因此lifetime、ω< |
如果是动态超极化率,激发态的寿命是用一个常见的现象学阻尼参数来近似的。这个值最好通过拟合分子的吸收数据来获得,不过这个值在相似的分子之间变化不大,因此估计值并不困难。文献" | 如果是动态超极化率,激发态的寿命是用一个常见的现象学阻尼参数来近似的。这个值最好通过拟合分子的吸收数据来获得,不过这个值在相似的分子之间变化不大,因此估计值并不困难。文献" | ||
行 35: | 行 35: | ||
SCM - Output,窗口底部输入dipole moment回车搜索得到: | SCM - Output,窗口底部输入dipole moment回车搜索得到: | ||
< | < | ||
- | |||
| | ||
行 64: | 行 63: | ||
| | ||
</ | </ | ||
- | 这里分别列出了极化率的实部与虚部,各自都是3*3的矩阵,分别对应$α_{ij}$,其中i,j=x,y,z。 | + | 这里分别列出了三次迭代的极化率的实部与虚部结果(看最后一个即可),各自都是3*3的矩阵,分别对应α<sub>ij</ |
可以看到 | 可以看到 | ||
- | * $α_{xx}$=20.5826 a.u. | + | * α<sub>xx</ |
- | * $α_{zz}$=8.34178 a.u. | + | * α<sub>zz</ |
与文献中20.59、8.21一致。 | 与文献中20.59、8.21一致。 | ||
行 74: | 行 73: | ||
在out窗口,搜索 Hyperpolarizability tensor: | 在out窗口,搜索 Hyperpolarizability tensor: | ||
< | < | ||
- | Hyperpolarizability tensor: | + | Hyperpolarizability tensor: |
| | ||
| | ||
行 102: | 行 101: | ||
Y | Y | ||
Z | Z | ||
+ | | ||
</ | </ | ||
- | 这里同样地列出了β张量的实部与虚部。β是三阶张量,有三个下标$β_{ijk}$,i,j,k=x,y,z。其中 | + | 这里同样地列出了β张量的实部与虚部。β是三阶张量,有三个下标β<sub>ijk</ |
- | * $β_{zxx}$ = -345.0757 a.u. | + | * β<sub>zxx</ |
- | * $β_{zzz}$ = -360.9820 a.u. | + | * β<sub>zzz</ |
与文中结果大致一致,实际上与Dalton的结果更接近。不同版本的ADF计算方法也不同,尤其是2016到2019每个版本之间均有较大调整。 | 与文中结果大致一致,实际上与Dalton的结果更接近。不同版本的ADF计算方法也不同,尤其是2016到2019每个版本之间均有较大调整。 | ||
- | |||
====二阶超极化率γ==== | ====二阶超极化率γ==== | ||
在out中搜索Second hyperpolarizability tensor: | 在out中搜索Second hyperpolarizability tensor: | ||
行 190: | 行 189: | ||
| | ||
</ | </ | ||
- | 同样列出了二阶超极化率的实部与虚部。γ为四阶张量,因此有四个下标$γ_{ijkl}$,i,j,k,l=x,y,z。其中: | + | 同样列出了二阶超极化率的实部与虚部。γ为四阶张量,因此有四个下标γ<sub>ijkl</ |
- | * $γ_{xxxx}$ = -4238.9431 a.u. | + | * γ<sub>xxxx</ |
- | * $γ_{xxyy}$ = -1426.0317 a.u. | + | * γ<sub>xxyy</ |
- | * $γ_{xxzz}$ = 9904.3713 a.u. | + | * γ<sub>xxzz</ |
- | * $γ_{zzzz}$ = 29981.6352 a.u. | + | * γ<sub>zzzz</ |