计算化学公社

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

[Multiwfn使用咨询] 使用Multiwfn计算LEA时会出现两个完全不同的结果

[复制链接 Copy URL]

10

帖子

0

威望

67

eV
积分
77

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
老师您好,

最近我在用Multiwfn计算小分子的LEA值,但是在计算过程中发现一个问题,就是当我计算完一个分子的ESP后(默认条件下,即12-0),通过-1返回主菜单后,再直接计算LEA值(即12-2-4-0,这里用①代替),会与新开程序时计算得到的LEA值(这里用②代替)差很多。同样的,我在计算完ALIE后,通过-1返回主菜单再直接计算LEA的值(这里用③代替)也会与新开程序计算得到的LEA值差很多,不过值①与值③是完全相同的。在计算期间,输入的fch文件没有任何改动。

另外,在计算ALIE时则不会受到我上一步计算数值的影响。

不知道这是不是我计算步骤错误,我该采信哪个数值?

数值①和③即计算完ESP或ALIE后再计算得到的LEA:
Volume:   473.47773 Bohr^3  (  70.16218 Angstrom^3)
Estimated density according to mass and volume (M/V):    2.4633 g/cm^3
Minimal value:     -4.05953 eV,   Maximal value:     -2.69074 eV
Overall surface area:         327.54536 Bohr^2  (  91.72206 Angstrom^2)
Positive surface area:          0.00000 Bohr^2  (   0.00000 Angstrom^2)
Negative surface area:        327.54536 Bohr^2  (  91.72206 Angstrom^2)
Overall average value:  -0.1379001 a.u. (  -3.75245 eV,    -86.53368 kcal/mol)
Positive average value:        NaN a.u. (       NaN eV,          NaN kcal/mol)
Negative average value: -0.1379001 a.u. (  -3.75245 eV,    -86.53368 kcal/mol)
Overall variance:   0.000122 a.u.^2  (   0.09025 eV^2, 0.47994E+02 kcal/mol^2)
Positive variance:  0.000000 a.u.^2  (   0.00000 eV^2, 0.00000E+00 kcal/mol^2)
Negative variance:  0.000122 a.u.^2  (   0.09025 eV^2, 0.47994E+02 kcal/mol^2)


数值②即新开程序后计算得到的LEA:
olume:   473.47773 Bohr^3  (  70.16218 Angstrom^3)
Estimated density according to mass and volume (M/V):    2.4633 g/cm^3
Minimal value:    -52.31489 eV,   Maximal value:    -19.55196 eV
Overall surface area:         327.54536 Bohr^2  (  91.72206 Angstrom^2)
Positive surface area:          0.00000 Bohr^2  (   0.00000 Angstrom^2)
Negative surface area:        327.54536 Bohr^2  (  91.72206 Angstrom^2)
Overall average value:  -1.2132229 a.u. ( -33.01347 eV,   -761.30950 kcal/mol)
Positive average value:        NaN a.u. (       NaN eV,          NaN kcal/mol)
Negative average value: -1.2132229 a.u. ( -33.01347 eV,   -761.30950 kcal/mol)
Overall variance:   0.123165 a.u.^2  (  91.19884 eV^2, 0.48499E+05 kcal/mol^2)
Positive variance:  0.000000 a.u.^2  (   0.00000 eV^2, 0.00000E+00 kcal/mol^2)
Negative variance:  0.123165 a.u.^2  (  91.19884 eV^2, 0.48499E+05 kcal/mol^2)


因为得到的数值相差很大,暂时也找不到其它文献值做验证,所以向您进行请教。谢谢老师!

6万

帖子

99

威望

5万

eV
积分
120201

管理员

公社社长

2#
发表于 Post on 2020-7-29 23:32:50 | 只看该作者 Only view this author
注意这个提示


之前算其它格点数据时,高阶空轨道会被删除节约计算时间,自然会影响之后LEA的计算,因为LEA直接基于虚轨道算的。

但凡发现不同操作顺序下结果不同,一律以新开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

10

帖子

0

威望

67

eV
积分
77

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2020-7-30 11:18:18 | 只看该作者 Only view this author
本帖最后由 imew 于 2020-7-30 11:22 编辑

谢谢老师您的解答!

本版积分规则 Credits rule

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

GMT+8, 2025-8-18 05:19 , Processed in 0.165023 second(s), 22 queries , Gzip On.

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