用户工具

站点工具


adf:lihhyper

差别

这里会显示出您选择的修订版和当前版本之间的差别。

到此差别页面的链接

两侧同时换到之前的修订记录前一修订版
后一修订版
前一修订版
adf:lihhyper [2020/12/01 14:41] – [偶极矩] liu.junadf: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),选择该选项会同时计算α、β、γ),频率(ω<sub>1</sub>、ω<sub>2</sub>、ω<sub>3</sub>、liftime均为0,因为计算的是静态超极化率)。
  
 文中的计算结果为: 文中的计算结果为:
行 23: 行 23:
 {{ :adf:hyperpol002.png?650 }} {{ :adf:hyperpol002.png?650 }}
  
-由于计算的是静态超极化率,因此lifetime、$ω_1$$ω_2$$ω_3$均设置为0。如果计算动态极化、超极化,设置相应的频率、寿命即可。具体可以参考:[[adf:tpa2019]]+由于计算的是静态超极化率,因此lifetime、ω<sub>1</sub>、ω<sub>2</sub>、ω<sub>3</sub>均设置为0。如果计算动态极化、超极化,设置相应的频率、寿命即可。具体可以参考:[[adf:tpa2019]]
  
 如果是动态超极化率,激发态的寿命是用一个常见的现象学阻尼参数来近似的。这个值最好通过拟合分子的吸收数据来获得,不过这个值在相似的分子之间变化不大,因此估计值并不困难。文献"Journal of Chemical Physics 123, 174110 , 2005"中使用了0.004 Hatree。 如果是动态超极化率,激发态的寿命是用一个常见的现象学阻尼参数来近似的。这个值最好通过拟合分子的吸收数据来获得,不过这个值在相似的分子之间变化不大,因此估计值并不困难。文献"Journal of Chemical Physics 123, 174110 , 2005"中使用了0.004 Hatree。
行 63: 行 63:
  Isotropic        POLARIZABILITY        0.0000     atomic units  Isotropic        POLARIZABILITY        0.0000     atomic units
 </code> </code>
-这里分别列出了极化率的实部与虚部,各自都是3*3的矩阵,分别对应$α_{ij}$,其中i,j=x,y,z+这里分别列出了三次迭代的极化率的实部与虚部结果(看最后一个即可),各自都是3*3的矩阵,分别对应α<sub>ij</sub>,其中i,j=y,x,x
 可以看到 可以看到
-  * $α_{xx}$=20.5826 a.u. +  * α<sub>xx</sub>=20.5826 a.u. 
-  * $α_{zz}$=8.34178 a.u.+  * α<sub>zz</sub>=8.34178 a.u.
  
 与文献中20.59、8.21一致。 与文献中20.59、8.21一致。
行 73: 行 73:
 在out窗口,搜索 Hyperpolarizability tensor: 在out窗口,搜索 Hyperpolarizability tensor:
 <code> <code>
-Hyperpolarizability tensor:+ Hyperpolarizability tensor:
                                     Y              Z                                     Y              Z
         X         0.0000         0.0000      -345.0548         X         0.0000         0.0000      -345.0548
行 101: 行 101:
           Y         0.0000         0.0000         0.0000           Y         0.0000         0.0000         0.0000
           Z         0.0000         0.0000         0.0000           Z         0.0000         0.0000         0.0000
 + ------------------------------------------
 </code> </code>
-这里同样地列出了β张量的实部与虚部。β是三阶张量,有三个下标$β_{ijk}$,i,j,k=x,y,z。其中 +这里同样地列出了β张量的实部与虚部。β是三阶张量,有三个下标β<sub>ijk</sub>,i,j,k=x,y,z。其中 
-  * $β_{zxx}$ = -345.0757 a.u. +  * β<sub>zxx</sub> = -345.05 a.u. 
-  * $β_{zzz}$ = -360.9820 a.u.+  * β<sub>zzz</sub> = -360.9820 a.u.
  
 与文中结果大致一致,实际上与Dalton的结果更接近。不同版本的ADF计算方法也不同,尤其是2016到2019每个版本之间均有较大调整。 与文中结果大致一致,实际上与Dalton的结果更接近。不同版本的ADF计算方法也不同,尤其是2016到2019每个版本之间均有较大调整。
- 
 ====二阶超极化率γ==== ====二阶超极化率γ====
 在out中搜索Second hyperpolarizability tensor: 在out中搜索Second hyperpolarizability tensor:
行 189: 行 189:
                            0.0000         0.0000         0.0000                            0.0000         0.0000         0.0000
 </code> </code>
-同样列出了二阶超极化率的实部与虚部。γ为四阶张量,因此有四个下标$γ_{ijkl}$,i,j,k,l=x,y,z。其中: +同样列出了二阶超极化率的实部与虚部。γ为四阶张量,因此有四个下标γ<sub>ijkl</sub>,i,j,k,l=x,y,z。其中: 
-  * $γ_{xxxx}$ = -4238.9431 a.u. +  * γ<sub>xxxx</sub> = -4238.9431 a.u. 
-  * $γ_{xxyy}$ = -1426.0317 a.u. +  * γ<sub>xxyy</sub> = -1426.0317 a.u. 
-  * $γ_{xxzz}$ = 9904.3713 a.u. +  * γ<sub>xxzz</sub> = 9904.3713 a.u. 
-  * $γ_{zzzz}$ = 29981.6352 a.u.+  * γ<sub>zzzz</sub> = 29981.6352 a.u.
adf/lihhyper.1606804918.txt.gz · 最后更改: 2020/12/01 14:41 由 liu.jun

© 2014-2022 费米科技(京ICP备14023855号