Gaussian+Multiwfn算极化率相关值

**参考链接：**[使用Multiwfn分析Gaussian的极化率、超极化率的输出](http://sobereva.com/231 "使用Multiwfn分析Gaussian的极化率、超极化率的输出")
**文献举例：**（极化率各向异性、超极化率、HOMO-LUMO gap)
![](/media/202303/2023-03-04_195556_7989910.18794091057446694.png)
## 1.优化结构
>d 以下输入文件仅关键词供参考，请根据实际情况进行适当修改,最后空两行,对于自选多重度，是多次计算取能量最低的
```
%chk=SO4.chk
%nproc=12
%rwf=SO4.rwf

#p pbe1pbe/6-31g(d) opt pop=full

Title Card Required

-2 1
 S                  1.22767857   -0.17857143    0.00000000
 O                 -0.44232143   -0.17855085    0.00000000
 O                  1.78435488    0.60866687    1.36354954
 O                  1.78432614   -1.75306928    0.00000000
 O                  1.78435488    0.60866687   -1.36354954

```
## 2.计算极化性质
>d 注意#p输出，是pbe1不是pbel，可以在此处加pop=regular算轨道等
```
%chk=SO4_P.chk
%nprocshared=12
%rwf=SO4_P.rwf
#p pbe1pbe/6-31g(d) polar

Title Card Required

-2 1
 S                  0.00000000    0.00000000    0.00000000
 O                  0.87476200    0.87476200    0.87476200
 O                 -0.87476200   -0.87476200    0.87476200
 O                 -0.87476200    0.87476200   -0.87476200
 O                  0.87476200   -0.87476200   -0.87476200

```
##3.数据处理
下载上一部的的输出log文件,拖进Multiwfn
然后进入主功能200，选7，之后会看到下面的菜单：
>-3 Set unit in the output, current: a.u.
-2 Set output destination, current: Output to screen
-1 Toggle loading frequency-dependent result for options 1 and 7, current: No
0 Return
1 "Polar" + analytic 3-order deriv. (HF,DFT,semi-empirical)
2 "Polar" + analytic 2-order deriv. (MP2,TDDFT,CIS...)
3 "Polar=Cubic" + analytic 2-order deriv.
4 "Polar" + analytic 1-order deriv. (CISD,QCISD,CCSD,MP3,MP4(SDQ)...)
5 "Polar=DoubleNumer" or "Polar=EnOnly" + analytic 1-order deriv.
6 "Polar" + energy only (CCSD(T),QCISD(T),MP4(SDTQ),MP5...)
7 "Polar=gamma" + analytic 3-order deriv. (HF,DFT,semi-empirical)

默认情况下输出信息会显示到屏幕上，也可以选-2然后11来更改为输出到当前目录下的polar.txt当中。
当前例子用的是PBE1PBE，是DFT中常用的一种泛函，故对应的是前述第一种情况，因此当我们看到上面的菜单时应选1，然后马上会看到α和β的分析结果。
## 4.结果

![](/media/202303/2023-03-04_195137_3959810.1288609109914104.png)
>i P.S. 如果要看HOMO-LUMO gap,把fchk文件拖进Multiwfn,按0就会输出在文本框了：）

Gaussian+Multiwfn算极化率相关值

服务器已完成数据迁移，数据丢失请及时反馈站长，旧服务器数据保留至11月30日。

【温馨提示：本站文档可配置可见范围，如登录后可见、对特定群组可见等，看不到就是没权限】

注册码获取邮箱
work@liuyaoze.com

Gaussian+Multiwfn算极化率相关值

服务器已完成数据迁移，数据丢失请及时反馈站长，旧服务器数据保留至11月30日。

【温馨提示：本站文档可配置可见范围，如登录后可见、对特定群组可见等，看不到就是没权限】

注册码获取邮箱 work@liuyaoze.com

注册码获取邮箱
work@liuyaoze.com