计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: sobereva
打印 Print 上一主题 Last thread 下一主题 Next thread

[Multiwfn使用咨询] 简单Multiwfn使用问题和波函数分析问题答疑专贴

  [复制链接 Copy URL]

6万

帖子

99

威望

5万

eV
积分
120050

管理员

公社社长

91#
 楼主 Author| 发表于 Post on 2018-8-24 19:04:21 | 只看该作者 Only view this author
Theor_Comp 发表于 2018-8-24 17:20
Sob老师,能否用优化激发态过程中的最后一步的电子激发信息来直接算?也就是说,优化激发态最后一步得出 ...

原理上一样。但是Multiwfn不支持优化任务的输出文件做delta_r计算
北京科音自然科学研究中心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

195

帖子

0

威望

699

eV
积分
894

Level 4 (黑子)

92#
发表于 Post on 2018-8-24 19:20:31 | 只看该作者 Only view this author
sobereva 发表于 2018-8-24 19:04
原理上一样。但是Multiwfn不支持优化任务的输出文件做delta_r计算

Sob老师,我是这么想的:
1. 将激发态优化得到的输出文件中只保留最后一步的结果,其他部分删掉,那这样能用Multiwfn来做激发态结构下的delta_r了吗?

2. 如果用所得的结构做一次电子激发计算,那使用激发态的结构进行TD时,关键词是# TD b3lyp/6-31G* scrf=(iefpcm,solvent=water)还是# TD=eqsolv b3lyp/6-31G* scrf=(iefpcm,solvent=water),也就是说这种情况下要不要加eqsolv?

谢谢您!

6万

帖子

99

威望

5万

eV
积分
120050

管理员

公社社长

93#
 楼主 Author| 发表于 Post on 2018-8-24 19:37:02 | 只看该作者 Only view this author
Theor_Comp 发表于 2018-8-24 19:20
Sob老师,我是这么想的:
1. 将激发态优化得到的输出文件中只保留最后一步的结果,其他部分删掉,那这样 ...

1 如果删不好,Multiwfn照样识别不了。而且带着IOp(9/40=x)做优化,每一步都会输出一大坨信息,输出文件也会颇大。

2 写eqsolv对当前情况从原理上有益
北京科音自然科学研究中心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

146

帖子

0

威望

1329

eV
积分
1475

Level 4 (黑子)

94#
发表于 Post on 2018-8-31 08:51:12 | 只看该作者 Only view this author
本帖最后由 刘小英雄steven 于 2018-8-31 08:52 编辑

老师,这是手册里的EDD图的分析,我想知道O和H 之间形成了为啥有两块区域,一块是电子密度的富集代表作用力增强可以理解,另一块电子密度减少反应了O-H之间怎样的作用?是不是原子间只要有大块的富集区就代表作用力增强呢?谢谢老师

EDD图分析.jpg (117.77 KB, 下载次数 Times of downloads: 150)

EDD图分析.jpg

6万

帖子

99

威望

5万

eV
积分
120050

管理员

公社社长

95#
 楼主 Author| 发表于 Post on 2018-8-31 16:41:58 | 只看该作者 Only view this author
刘小英雄steven 发表于 2018-8-31 08:51
老师,这是手册里的EDD图的分析,我想知道O和H 之间形成了为啥有两块区域,一块是电子密度的富集代表作用力 ...

电子有増的地方必定伴随着电子减少的地方。减少的那块可以认为是氧上的电子云分布发生了极化、轨道的杂化,导致相应区域电子数减少所致

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
刘小英雄steven + 2 谢谢

查看全部评分 View all ratings

北京科音自然科学研究中心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

195

帖子

0

威望

699

eV
积分
894

Level 4 (黑子)

96#
发表于 Post on 2018-9-1 20:16:46 | 只看该作者 Only view this author
Sob老师,在使用Multiwfn绘制荧光光谱时,您的教程中“因此和绘制吸收光谱的关键不同点是要把S0至S1以上激发态的振子强度都设为0。做法就是先优化出S1的结构,然后在这个结构下做一次电子激发计算(默认的nstates=3实际上已经够了),然后将输出文件载入Multiwfn...“,
为什么优化S1后还需要再在该结构下做一次电子激发计算,优化S1过程中不是已经给出了吗?直接用优化S1的输出文件不可以吗?

490

帖子

2

威望

5122

eV
积分
5652

Level 6 (一方通行)

97#
发表于 Post on 2018-9-1 20:21:14 | 只看该作者 Only view this author
Theor_Comp 发表于 2018-9-1 20:16
Sob老师,在使用Multiwfn绘制荧光光谱时,您的教程中“因此和绘制吸收光谱的关键不同点是要把S0至S1以上激 ...

上面的回答不是已经说了吗。
Multiwfn不支持优化任务。我猜Multiwfn在读取输出文件中的组态系数的时候,是通过“->”和"<-"这两个符号来识别的,因此如果是激发态优化的任务,优化的每一步都会有这些组态系数需要读取,肯定就会出错了。
以上是个人的猜测,如有误,还请sob老师指出。

195

帖子

0

威望

699

eV
积分
894

Level 4 (黑子)

98#
发表于 Post on 2018-9-1 21:41:58 | 只看该作者 Only view this author
让你变成回忆 发表于 2018-9-1 20:21
上面的回答不是已经说了吗。
Multiwfn不支持优化任务。我猜Multiwfn在读取输出文件中的组态系数的时候, ...

但Multiwfn好像可以使用S1优化的输出文件来画荧光光谱,这个结果可不可以用?

490

帖子

2

威望

5122

eV
积分
5652

Level 6 (一方通行)

99#
发表于 Post on 2018-9-1 22:13:31 | 只看该作者 Only view this author
Theor_Comp 发表于 2018-9-1 21:41
但Multiwfn好像可以使用S1优化的输出文件来画荧光光谱,这个结果可不可以用?

我试了下好像是可以支持激发态优化的输出文件绘制光谱。而且激发能、阵子强度等信息的确也和优化最后一步的结果一致。我看了下sob老师上面的回答,应该是说一些分析激发态跃迁特征的功能,是不能支持优化的任务。

具体还是等着sob老师来回答。

6万

帖子

99

威望

5万

eV
积分
120050

管理员

公社社长

100#
 楼主 Author| 发表于 Post on 2018-9-2 02:25:39 | 只看该作者 Only view this author
Theor_Comp 发表于 2018-9-1 21:41
但Multiwfn好像可以使用S1优化的输出文件来画荧光光谱,这个结果可不可以用?

手册里有明示


不能用激发态优化任务的是主功能18(电子激发分析)中的那些需要组态系数的功能,比如空穴-电子分析,这和绘制光谱不是一码事
北京科音自然科学研究中心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

195

帖子

0

威望

699

eV
积分
894

Level 4 (黑子)

101#
发表于 Post on 2018-9-2 11:01:32 | 只看该作者 Only view this author

明白了,谢谢您!

210

帖子

0

威望

1853

eV
积分
2063

Level 5 (御坂)

102#
发表于 Post on 2018-9-10 13:58:00 | 只看该作者 Only view this author
请问sob老师,“做波函数分析时,一定要用小核赝势,而不能用大核赝势,否则价层密度、波函数不准确”的说法,有没有文献支持。谢谢老师。

6万

帖子

99

威望

5万

eV
积分
120050

管理员

公社社长

103#
 楼主 Author| 发表于 Post on 2018-9-10 18:41:11 | 只看该作者 Only view this author
alystone 发表于 2018-9-10 13:58
请问sob老师,“做波函数分析时,一定要用小核赝势,而不能用大核赝势,否则价层密度、波函数不准确”的说 ...

这个凭基本的逻辑和理论知识就能理解
非要参考文献可以看
J. Phys. Chem. A 2011, 115, 12879–12894
北京科音自然科学研究中心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

210

帖子

0

威望

1853

eV
积分
2063

Level 5 (御坂)

104#
发表于 Post on 2018-9-10 21:01:07 | 只看该作者 Only view this author
sobereva 发表于 2018-9-10 18:41
这个凭基本的逻辑和理论知识就能理解
非要参考文献可以看
J. Phys. Chem. A 2011, 115, 12879–12894

好的,谢谢sob老师。

62

帖子

0

威望

559

eV
积分
621

Level 4 (黑子)

105#
发表于 Post on 2018-9-11 10:27:27 | 只看该作者 Only view this author
几何优化后得到fchk文件,使用Multiwfn绘制分子表面时,出现了如下错误提示,请问是Multiwfn的设置有问题吗?多谢
Calculating electrostatic potential at surface vertices, please wait patiently
Since the input file type is .fch/fchk and "cubegenpath" parameter in settings.ini has been properly defined, now Multiwfn directly invokes cubegen to calculate electrostatic potential
Running: /home/ncvdse/gaussian/g16/cubegen   36 potential=SCF /home/ncvdse/chenglong/CS-LYS.fchk ESPresult.cub -5 h < cubegenpt.txt > nouseout
sh: 行 1: 63511 已放弃               (吐核)/home/ncvdse/gaussian/g16/cubegen 36 potential=SCF /home/ncvdse/chenglong/CS-LYS.fchk ESPresult.cub -5 h < cubegenpt.txt > nouseout
forrtl: severe (24): end-of-file during read, unit 10, file /home/ncvdse/chenglong/ESPresult.cub
Image              PC                Routine            Line        Source            
Multiwfn           00000000007D4BC7  Unknown               Unknown  Unknown
Multiwfn           00000000007FFA93  Unknown               Unknown  Unknown
Multiwfn           0000000000652909  Unknown               Unknown  Unknown
Multiwfn           00000000005457D1  Unknown               Unknown  Unknown
Multiwfn           00000000004091DE  Unknown               Unknown  Unknown
libc.so.6          00002B1290BE9445  Unknown               Unknown  Unknown
Multiwfn           00000000004090E9  Unknown               Unknown  Unknown


本版积分规则 Credits rule

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

GMT+8, 2025-8-12 12:11 , Processed in 1.131050 second(s), 26 queries , Gzip On.

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