计算化学公社

 找回密码 Forget password
 注册 Register
Views: 12824|回复 Reply: 16
打印 Print 上一主题 Last thread 下一主题 Next thread

[Multiwfn使用咨询] Multiwfn计算极化率闪退问题

[复制链接 Copy URL]

14

帖子

0

威望

1175

eV
积分
1189

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
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因为太大无法上传)。

cellobiose.out

1.32 MB, 下载次数 Times of downloads: 11

cellobiose.gjf

3.42 KB, 下载次数 Times of downloads: 25

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

2#
发表于 Post on 2017-11-10 17:48:08 | 只看该作者 Only view this author
你的这个任务都没做第二步的polar计算,优化之后直接就结束了。
用优化的结构,做一个标准的polar计算,别分两步。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

207

帖子

1

威望

3855

eV
积分
4082

Level 6 (一方通行)

3#
发表于 Post on 2017-11-10 19:41:41 | 只看该作者 Only view this author
原输入文件第二步错了,是link1不是link 1

14

帖子

0

威望

1175

eV
积分
1189

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2017-11-10 20:27:36 | 只看该作者 Only view this author
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还是闪退,原因应该是和现在一样,所以我才采用两步法进行的计算。
刚才我采用您的建议,再次进行了计算,结果还是一样。
具体输入和输出文件见附件。

cellobiose2.gjf

128 Bytes, 下载次数 Times of downloads: 18

cellobiose2.out

84.25 KB, 下载次数 Times of downloads: 18

14

帖子

0

威望

1175

eV
积分
1189

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2017-11-10 20:29:50 | 只看该作者 Only view this author
bomsaude 发表于 2017-11-10 19:41
原输入文件第二步错了,是link1不是link 1

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

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

6#
发表于 Post on 2017-11-10 21:44:16 | 只看该作者 Only view this author
用最新版本,解析得很正常。启动后载入文件,依次输入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
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

93

帖子

0

威望

1782

eV
积分
1875

Level 5 (御坂)

7#
发表于 Post on 2017-11-11 08:20:06 | 只看该作者 Only view this author
请问这个计算对二聚体适用吗?另外不同的结合方式所形成的dimer1和dimer2,在极化率上会体现差异吗?谢谢!

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

8#
发表于 Post on 2017-11-11 12:07:28 | 只看该作者 Only view this author
wei 发表于 2017-11-11 08:20
请问这个计算对二聚体适用吗?另外不同的结合方式所形成的dimer1和dimer2,在极化率上会体现差异吗?谢谢!

可以用。会
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

101

帖子

1

威望

645

eV
积分
766

Level 4 (黑子)

9#
发表于 Post on 2017-11-11 23:18:40 | 只看该作者 Only view this author

请问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的办法

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

10#
发表于 Post on 2017-11-12 01:09:15 | 只看该作者 Only view this author
五十八 发表于 2017-11-11 23:18
请问sobereva老师,我们现在想计算permittivity(static dielectric properties),使用的方法是 (ε-1 ...

这部分我不太清楚
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

14

帖子

0

威望

1175

eV
积分
1189

Level 4 (黑子)

11#
 楼主 Author| 发表于 Post on 2017-11-12 09:32:10 | 只看该作者 Only view this author
sobereva 发表于 2017-11-10 21:44
用最新版本,解析得很正常。启动后载入文件,依次输入200、7、1
Dipole moment:
X,Y,Z=    0.777182   -0 ...

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

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

12#
发表于 Post on 2017-11-12 09:57:09 | 只看该作者 Only view this author
flowerstudy 发表于 2017-11-12 09:32
sobereva老师,不知道您说的最新版指的是那个版本?
可惜我这里只有g9A和g9c版本,得不到所需结果。
不 ...


Multiwfn最新版本就是现在官网上能下载到的最新版本
你最后上传的输出文件无问题
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

14

帖子

0

威望

1175

eV
积分
1189

Level 4 (黑子)

13#
 楼主 Author| 发表于 Post on 2017-11-12 18:31:57 | 只看该作者 Only view this author
sobereva 发表于 2017-11-12 09:57
Multiwfn最新版本就是现在官网上能下载到的最新版本
你最后上传的输出文件无问题

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

14

帖子

0

威望

1175

eV
积分
1189

Level 4 (黑子)

14#
 楼主 Author| 发表于 Post on 2017-11-12 18:36:32 | 只看该作者 Only view this author
bomsaude 发表于 2017-11-10 19:41
原输入文件第二步错了,是link1不是link 1

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

167

帖子

0

威望

1039

eV
积分
1206

Level 4 (黑子)

15#
发表于 Post on 2019-1-16 14:04:52 | 只看该作者 Only view this author
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

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-15 18:43 , Processed in 0.447179 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list