计算化学公社

标题: Multiwfn计算极化率闪退问题 [打印本页]

作者
Author: flowerstudy 时间: 2017-11-10 11:38
标题: Multiwfn计算极化率闪退问题
Sobereva老师您好！
这两天我在用Multiwfn计算物质的极化率（使用的方法是您的帖子：http://sobereva.com/231）。整个量化计算过程中没有出现任何问题，只是在用到Multiwfn分析时就会出现闪退的问题。

为了排除原因，我用Gaussian 09 Revision C.01和Gaussian 09 Revision A.02两个版本都进行了计算，分析时我也用了Multiwfn3.3.7和Multiwfn3.4.1分别进行了分析，但都出现了闪退。

最后，我打开了源代码进行了分析，原来程序中需要在.out文件中查找到“SCF Polarizability”，“Isotropic polarizability”等才可进行接下来的计算，但是我也在.out中查找过，压根就没有出现“SCF Polarizability”，“Isotropic polarizability”等语句，我想这也是程序闪退的原因吧。

现在想问Sobereva老师该采用怎样的解决方案？附件是我的量化计算的语句及计算输出结果，方便您分析问题之所在（*.chk因为太大无法上传）。

作者
Author: sobereva 时间: 2017-11-10 17:48
你的这个任务都没做第二步的polar计算，优化之后直接就结束了。
用优化的结构，做一个标准的polar计算，别分两步。

作者
Author: bomsaude 时间: 2017-11-10 19:41
原输入文件第二步错了，是link1不是link 1

作者
Author: flowerstudy 时间: 2017-11-10 20:27

sobereva 发表于 2017-11-10 17:48
你的这个任务都没做第二步的polar计算，优化之后直接就结束了。
用优化的结构，做一个标准的polar计算，别 ...

谢谢，sobereva老师！
其实，之前我的输入是这样写的：
%nprocshared=4
%mem=1GB
%chk=D:\g09\cellobiose.chk
#p opt b3lyp/6-31g(d) polar maxdisk=10GB geom=connectivity
分子信息
但是结果就是和现在的结果一样，Multiwfn还是闪退，原因应该是和现在一样，所以我才采用两步法进行的计算。
刚才我采用您的建议，再次进行了计算，结果还是一样。
具体输入和输出文件见附件。

作者
Author: flowerstudy 时间: 2017-11-10 20:29

bomsaude 发表于 2017-11-10 19:41
原输入文件第二步错了，是link1不是link 1

谢谢你，指出我的错误！但刚才我分两步进行了计算，过程中并没有使用link1，但结果还是和之前一样！

作者
Author: sobereva 时间: 2017-11-10 21:44
用最新版本，解析得很正常。启动后载入文件，依次输入200、7、1
Dipole moment:
X,Y,Z= 0.777182 -0.311715 -0.166146 Norm= 0.853688

Static polarizability:
XX=    193.121000
XY=    -2.467410
YY=    172.012000
XZ=    -1.212630
YZ=    -1.134180
ZZ=    119.313000
Isotropic average polarizability:    161.482000
Polarizability anisotropy (definition 1):    66.043414
Polarizability anisotropy (definition 2):    65.842218
Eigenvalues of polarizability tensor: 119.26689 171.75794 193.42118
Polarizability anisotropy (definition 3):    47.908766

Note: It is well known that the sign of hyperpolarizability of Gaussian 09 sh
ld be multiplied by -1, the outputs below have already been corrected

Static first hyperpolarizability:
XXX=    -134.302000
XXY=    -10.618300
XYY=    24.248000
YYY=    21.799400
XXZ=    -11.513500
XYZ=    -4.273800
YYZ=       1.642120
XZZ=    -4.240150
YZZ=       1.164080
ZZZ=    22.728100

Beta_X=    -114.29415  Beta_Y=    12.34518  Beta_Z=    12.85672
Magnitude of first hyperpolarizability:    115.675630
Projection of beta on dipole moment:    -111.061274
Beta ||    :    -66.636764
Beta ||(z)  :       7.714032
Beta _|_(z) :       2.571344

作者
Author: wei 时间: 2017-11-11 08:20
请问这个计算对二聚体适用吗？另外不同的结合方式所形成的dimer1和dimer2，在极化率上会体现差异吗？谢谢！

作者
Author: sobereva 时间: 2017-11-11 12:07

wei 发表于 2017-11-11 08:20
请问这个计算对二聚体适用吗？另外不同的结合方式所形成的dimer1和dimer2，在极化率上会体现差异吗？谢谢！

可以用。会

作者
Author: 五十八 时间: 2017-11-11 23:18

sobereva 发表于 2017-11-11 12:07
可以用。会

请问sobereva老师，我们现在想计算permittivity（static dielectric properties），使用的方法是（ε-1）/（ε+2） = （ρ/M）*Pm
其中pm = （Na/（3*∑0））*（α+(μ^2)/(3kT）)
而α=1/3Tr[axx ayy azz]（使用au）
但是，VASP给出的结果得到了很奇怪的PM ，使用的是IDIPOL=（1，2，3，4）（4， EXTRA_FIELD=0==》GS）估计是不对
有没有VASP 不用LCACLPOL或者RPA（GW）或者DFPT得到PM的办法

作者
Author: sobereva 时间: 2017-11-12 01:09

五十八发表于 2017-11-11 23:18
请问sobereva老师，我们现在想计算permittivity（static dielectric properties），使用的方法是（ε-1 ...

这部分我不太清楚

作者
Author: flowerstudy 时间: 2017-11-12 09:32

sobereva 发表于 2017-11-10 21:44
用最新版本，解析得很正常。启动后载入文件，依次输入200、7、1
Dipole moment:
X,Y,Z= 0.777182 -0 ...

sobereva老师，不知道您说的最新版指的是那个版本？
可惜我这里只有g9A和g9c版本，得不到所需结果。
不知道是我的输入问题呢，还是这两个版本就不能计算polar呢？

作者
Author: sobereva 时间: 2017-11-12 09:57

flowerstudy 发表于 2017-11-12 09:32
sobereva老师，不知道您说的最新版指的是那个版本？
可惜我这里只有g9A和g9c版本，得不到所需结果。
不 ...

Multiwfn最新版本就是现在官网上能下载到的最新版本
你最后上传的输出文件无问题

作者
Author: flowerstudy 时间: 2017-11-12 18:31

sobereva 发表于 2017-11-12 09:57
Multiwfn最新版本就是现在官网上能下载到的最新版本
你最后上传的输出文件无问题

好的，谢谢sobereva老师，是我大意了。问题已解决！

作者
Author: flowerstudy 时间: 2017-11-12 18:36

bomsaude 发表于 2017-11-10 19:41
原输入文件第二步错了，是link1不是link 1

好的，谢谢bomsaude！是我大意了！
其实输入的核心错误就是您所提到的那个link1之间加了空格，导致第二层计算无法进行了，到第一层算完后就结束了！

作者
Author: qianbao 时间: 2019-1-16 14:04

sobereva 发表于 2017-11-10 21:44
用最新版本，解析得很正常。启动后载入文件，依次输入200、7、1
Dipole moment:
X,Y,Z= 0.777182 -0 ...

借该贴想请问一下sob老师，这里的定义一和定义二需要看哪一个呢，请问指得都是极化率的各向异性吗？
Polarizability anisotropy (definition 1): 66.043414
Polarizability anisotropy (definition 2): 65.842218

作者
Author: sobereva 时间: 2019-1-16 23:25

qianbao 发表于 2019-1-16 14:04
借该贴想请问一下sob老师，这里的定义一和定义二需要看哪一个呢，请问指得都是极化率的各向异性吗？
Pol ...

都可以，手册里都写了公式，你写文章的时候也写明怎么定义的就行。
是

作者
Author: qianbao 时间: 2019-1-17 08:17

sobereva 发表于 2019-1-16 23:25
都可以，手册里都写了公式，你写文章的时候也写明怎么定义的就行。
是

收到，谢谢老师！

欢迎光临计算化学公社 (http://bbs.keinsci.com/)