计算化学公社

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

[综合交流] NWChem计算碰撞截面问题

[复制链接 Copy URL]

6万

帖子

99

威望

5万

eV
积分
120189

管理员

公社社长

31#
发表于 Post on 2017-1-13 14:07:13 | 只看该作者 Only view this author
SAI 发表于 2017-1-13 09:53
sob老师,在做Mulliken分析这种波函数分析,您的博文和论文对基组的选择我好好地看过了。但是对于HF/DFT/ ...


DFT完全合适
如今不给出什么依据用MP2才显得莫名其妙,比DFT又耗时得多,精度多数情况还不如DFT
北京科音自然科学研究中心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

96

帖子

0

威望

4461

eV
积分
4557

Level 6 (一方通行)

32#
 楼主 Author| 发表于 Post on 2017-2-25 17:19:41 | 只看该作者 Only view this author
sobereva 发表于 2017-1-13 14:07
DFT完全合适
如今不给出什么依据用MP2才显得莫名其妙,比DFT又耗时得多,精度多数情况还不如DFT

我在自己算Mulliken population的时候,HF和MP2的计算结果有较好的一致性,但是B3LYP算出来的结果差异比较大,尤其是轨道能量数值明显较前两种算出来的结果要小,请问sob老师这是什么原因呢?

6万

帖子

99

威望

5万

eV
积分
120189

管理员

公社社长

33#
发表于 Post on 2017-2-26 02:04:40 | 只看该作者 Only view this author
SAI 发表于 2017-2-25 17:19
我在自己算Mulliken population的时候,HF和MP2的计算结果有较好的一致性,但是B3LYP算出来的结果差异比 ...


如果MP2的结果和HF一模一样,说明MP2下计算mulliken布居的时候其实用的是HF波函数。在MP2下做布居分析则需要产生MP2波函数,需要花很多额外时间。
MP2和B3LYP都考虑了电子相关,其结果应当接近,而都与HF有明显不同。即HF的结果是不合理的。
北京科音自然科学研究中心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

96

帖子

0

威望

4461

eV
积分
4557

Level 6 (一方通行)

34#
 楼主 Author| 发表于 Post on 2017-2-26 16:53:42 | 只看该作者 Only view this author
sobereva 发表于 2017-2-26 02:04
如果MP2的结果和HF一模一样,说明MP2下计算mulliken布居的时候其实用的是HF波函数。在MP2下做布居分析 ...

MP2和HF的结果只是比较一致,不是一模一样。相同基组下,DFT优化并算出来的轨道能量数值整体低于MP2优化算出来的结果。对于我要用的半经验方法来说,MP2或HF算出来的结果明显更加符合实验数据。那么我是不是只能说此处只是该半经验方法更适合MP2的布居分析结果?或者该如何解释呢

6万

帖子

99

威望

5万

eV
积分
120189

管理员

公社社长

35#
发表于 Post on 2017-2-26 17:10:14 | 只看该作者 Only view this author
SAI 发表于 2017-2-26 16:53
MP2和HF的结果只是比较一致,不是一模一样。相同基组下,DFT优化并算出来的轨道能量数值整体低于MP2优化 ...


MP2根本没有轨道能量这一说,你如果看到MP2轨道能量,那一定是HF轨道能量,如果和HF轨道能量不同,仔细对照输出
如果搞不清楚,就用DFT,甭用MP2,更甭用HF
北京科音自然科学研究中心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

96

帖子

0

威望

4461

eV
积分
4557

Level 6 (一方通行)

36#
 楼主 Author| 发表于 Post on 2017-2-26 17:59:10 | 只看该作者 Only view this author
sobereva 发表于 2017-2-26 17:10
MP2根本没有轨道能量这一说,你如果看到MP2轨道能量,那一定是HF轨道能量,如果和HF轨道能量不同,仔细 ...

那么MP2和HF的不同,是因为MP2和HF结构优化的不同吗?

6万

帖子

99

威望

5万

eV
积分
120189

管理员

公社社长

37#
发表于 Post on 2017-2-27 03:29:58 | 只看该作者 Only view this author
SAI 发表于 2017-2-26 17:59
那么MP2和HF的不同,是因为MP2和HF结构优化的不同吗?

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

96

帖子

0

威望

4461

eV
积分
4557

Level 6 (一方通行)

38#
 楼主 Author| 发表于 Post on 2017-2-28 09:55:03 | 只看该作者 Only view this author

sob老师,我下来验证了一下。我从输出中拷出了一个优化好的结构,然后在相同的基组下,MP2和HF方法pop=full看到的结果确实是相同的,那么结论正是如您所说的那样。可是还是同样的结构和基组,B3LYP下的结果依旧相差很远。我算的分子是丙酮,最外层的轨道HF和MP2算出来的能量是0.41 a.u. 左右,而B3LYP下算出来是0.24 a.u.左右,这个差异应该很不合理吧?

6万

帖子

99

威望

5万

eV
积分
120189

管理员

公社社长

39#
发表于 Post on 2017-2-28 10:09:59 | 只看该作者 Only view this author
SAI 发表于 2017-2-28 09:55
sob老师,我下来验证了一下。我从输出中拷出了一个优化好的结构,然后在相同的基组下,MP2和HF方法pop=fu ...

本来就应该有明显差异,HF没考虑电子相关,而B3LYP考虑了。
北京科音自然科学研究中心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

96

帖子

0

威望

4461

eV
积分
4557

Level 6 (一方通行)

40#
 楼主 Author| 发表于 Post on 2017-2-28 18:55:29 | 只看该作者 Only view this author
sobereva 发表于 2017-2-28 10:09
本来就应该有明显差异,HF没考虑电子相关,而B3LYP考虑了。

终于搞明白了,谢谢sob老师的讲解!

96

帖子

0

威望

4461

eV
积分
4557

Level 6 (一方通行)

41#
 楼主 Author| 发表于 Post on 2017-3-9 18:29:25 | 只看该作者 Only view this author
Mikasa 发表于 2016-11-10 09:55
1. 搜关键词 Quasi-Classical Trajectory (QCT),或者direct dynamics,或者看Prof. William L. Hase的文章 ...

您好,时隔很久终于渐渐的多了解了一些dynamics的相关知识,多谢您曾经在帖子下的指导。我以dynamics和cross section为关键词搜了一下文献,虽然大多看不太懂,但是也明白cross section可以由此计算。想再次请教您,如果我想入手这一部分的内容,自己编写程序还是利用现有的软件更好呢? 现有的软件以什么最为适合这一目标呢?非常感谢

220

帖子

0

威望

6004

eV
积分
6224

Level 6 (一方通行)

跳跳猪

42#
发表于 Post on 2017-3-10 09:17:08 | 只看该作者 Only view this author
SAI 发表于 2017-3-9 18:29
您好,时隔很久终于渐渐的多了解了一些dynamics的相关知识,多谢您曾经在帖子下的指导。我以dynamics和cr ...

自己编程可能工作量会比较大。
现有软件的话,我了解的不是太多,不过我推荐Truhlar教授课题组的ANT(https://comp.chem.umn.edu/ant/),可以跑直接动力学(超级耗时),也可以先搭建势能面之后再跑动力学(搭建PES略复杂)。
但是这种dynamics的计算量是特别大的,需要考虑一下计算资源能否接受。

评分 Rate

参与人数
Participants 1
eV +4 收起 理由
Reason
SAI + 4 谢谢!

查看全部评分 View all ratings

流年似水,浮生如梦。

96

帖子

0

威望

4461

eV
积分
4557

Level 6 (一方通行)

43#
 楼主 Author| 发表于 Post on 2017-3-10 10:00:20 | 只看该作者 Only view this author
本帖最后由 SAI 于 2017-3-10 11:04 编辑
Mikasa 发表于 2017-3-10 09:17
自己编程可能工作量会比较大。
现有软件的话,我了解的不是太多,不过我推荐Truhlar教授课题组的ANT(htt ...

您好,我去看了一下手册,好像ANT能算的cross section很多。我需要计算的主要是electron collision/interations with gas, 这一类截面是不是要简单一些?  顺便呼叫一下sob老师@sobereva

本版积分规则 Credits rule

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

GMT+8, 2025-8-18 00:08 , Processed in 0.226623 second(s), 21 queries , Gzip On.

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