这里会显示出您选择的修订版和当前版本之间的差别。
两侧同时换到之前的修订记录前一修订版后一修订版 | 前一修订版 | ||
adf:lihhyper [2020/12/01 14:49] – [一阶超极化率β] 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、ω< |
如果是动态超极化率,激发态的寿命是用一个常见的现象学阻尼参数来近似的。这个值最好通过拟合分子的吸收数据来获得,不过这个值在相似的分子之间变化不大,因此估计值并不困难。文献" | 如果是动态超极化率,激发态的寿命是用一个常见的现象学阻尼参数来近似的。这个值最好通过拟合分子的吸收数据来获得,不过这个值在相似的分子之间变化不大,因此估计值并不困难。文献" | ||
行 46: | 行 46: | ||
极化率α是一个二阶张量,搜索 Polarizability tensor: | 极化率α是一个二阶张量,搜索 Polarizability tensor: | ||
< | < | ||
- | THE DIPOLE-DIPOLE POLARIZABILITY TENSOR: | + | |
+ | X | ||
+ | | ||
+ | 0.00000 | ||
+ | 0.00000 | ||
- | FREQUENCY CYCLE NR.: 1 AT 0.0000000000 HARTREES | + | Isotropic |
- | Y Z X | + | |
- | 20.582595 | + | |
- | 0.000000 | + | |
- | 0.000000 | + | X Y |
+ | 0.00000 0.00000 0.00000 | ||
+ | 0.00000 0.00000 0.00000 | ||
+ | 0.00000 0.00000 0.00000 | ||
- | + | Isotropic | |
- | FREQUENCY CYCLE NR.: 2 AT 0.0000000000 HARTREES | + | |
- | Y Z X | + | |
- | | + | |
- | 0.000000 | + | |
- | 0.000000 | + | |
- | + | ||
- | + | ||
- | FREQUENCY CYCLE NR.: 3 AT 0.0000000000 HARTREES | + | |
- | Y Z X | + | |
- | | + | |
- | 0.000000 | + | |
- | 0.000000 | + | |
</ | </ | ||
- | 这里分别列出了三次迭代的极化率的实部与虚部结果(看最后一个即可),各自都是3*3的矩阵,分别对应$α_{ij}$,其中i,j=y,x,x。 | + | 这里分别列出了三次迭代的极化率的实部与虚部结果(看最后一个即可),各自都是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一致。 | ||
行 78: | 行 73: | ||
在out窗口,搜索 Hyperpolarizability tensor: | 在out窗口,搜索 Hyperpolarizability tensor: | ||
< | < | ||
- | The STATIC hyperpolarizability | + | |
- | Non-zero components of beta: | + | X Y Z |
- | + | X X 0.0000 0.0000 | |
- | beta(z z z ) = | + | |
- | beta(y z z ) = -0.11166E-06 | + | Z |
- | beta(x z z ) = 0.11461E-06 | + | |
- | beta(z y z ) = -0.11166E-06 | + | |
- | beta(y y z ) = -345.05 -172.53 -0.29810E-29 | + | |
- | beta(z x z ) = 0.11461E-06 | + | Z 0.0000 |
- | beta(x x z ) = -345.05 | + | |
- | beta(z z y ) = | + | |
- | beta(y z y ) = -345.05 | + | |
- | beta(z y y ) = -345.05 -172.53 -0.29810E-29 | + | Z 0.0000 0.0000 |
- | beta(z z x ) = 0.11461E-06 | + | |
- | beta(x z x ) = -345.05 -172.53 -0.29810E-29 | + | IMAGINARY HYPERPOLARIZABILITY |
- | beta(z x x ) = -360.95 -172.53 -0.29810E-29 | + | ------------------------------------------ |
+ | | ||
+ | | ||
+ | | ||
+ | Z 0.0000 | ||
+ | |||
+ | | ||
+ | Y 0.0000 | ||
+ | | ||
+ | |||
+ | | ||
+ | | ||
+ | Z 0.0000 | ||
+ | | ||
</ | </ | ||
- | 这里同样地列出了β张量的实部与虚部。β是三阶张量,有三个下标$β_{ijk}$,i,j,k=x,y,z。其中 | + | 这里同样地列出了β张量的实部与虚部。β是三阶张量,有三个下标β<sub>ijk</ |
- | * $β_{zxx}$ = -345.05 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: | ||
行 182: | 行 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</ |