计算化学公社

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

[综合交流] 简单量子化学问题答疑专帖(旧版,已关帖)

   关闭 [复制链接 Copy URL]

212

帖子

0

威望

948

eV
积分
1160

Level 4 (黑子)

13006#
发表于 Post on 2021-4-12 16:24:22 | 只看该作者 Only view this author
sobereva 发表于 2021-4-12 08:43
一般如此
有的时候也可以把电子能量变化、焓的变化和自由能的变化都给出,用来更深入分析不同成分对自由 ...

多谢老师!请问老师,gaussian16直接写关键词freq得到的“Sum of electronic and thermal Free Energies=”吉布斯自由能是0K还是273.15K,还是298.15K?

57

帖子

0

威望

479

eV
积分
536

Level 4 (黑子)

13007#
发表于 Post on 2021-4-12 16:37:00 | 只看该作者 Only view this author
黄舒伟 发表于 2021-4-12 15:57
请问我用G16的 # pm7 opt scrf=solvent=water  这个参数算出来的RESP电荷很离谱,我是应该用G16生成的log文 ...

懂了,这里应该继续用Log文件去算。刚才用一个小分子去验证了,发现:通过比较AAA_b3lyp_opt.chg和 AAA_pm7_opt.chg和AAA1_pm7_opt_b3lyp.chg
发现,半经验的RESP对于这种程度的体系是不可以接受的,但是可以用来opt体系,然后用b3lyp直接计算单点得到的resp和耗时长很多的opt b3lyp差不大多.(AAA_pm7_opt.chg和AAA1_pm7_opt_b3lyp.chg相差不大)

108

帖子

0

威望

2879

eV
积分
2987

Level 5 (御坂)

13008#
发表于 Post on 2021-4-12 22:16:29 | 只看该作者 Only view this author

不同理论级别下计算的HOMO-LUMO gap不同

尊敬的Sob,各位老师,大家晚上好。
我用Gaussian在不同理论计算级别下得到不同.fch文件,用Multiwfn主功能0计算得到的HOMO-LUMO gap数值不同,请问这个时候怎么选择哪个级别下的HOMO-LUMO gap?
我在计算HOMO-LUMO gap时,有时还得到α-HOMO-LUMO gap和β-HOMO-LUMO gap,数值也略有不同,讨论分子化学稳定性的时候,该选取哪个?
在此先谢过各位老师啦~

1665

帖子

5

威望

4788

eV
积分
6553

Level 6 (一方通行)

喵星人

13009#
发表于 Post on 2021-4-12 22:22:31 | 只看该作者 Only view this author
当然不同了,而且讨论稳定性这么复杂的问题怎么可能靠一个homo-lumo gap

6万

帖子

99

威望

6万

eV
积分
125179

管理员

公社社长

13010#
 楼主 Author| 发表于 Post on 2021-4-13 09:03:29 | 只看该作者 Only view this author
sxh3231796 发表于 2021-4-12 15:12
请问下sob老师,用同样的基组和方法优化出来的结果有细微的差距,不是很大,能量相差在0.006个HF能量,是可 ...

那叫Hartree,不叫HF能量
0.006 Hartree的差异一点也不小,要么计算设定不同,要么优化到了不同结构,要么没有都收敛到稳定波函数
北京科音自然科学研究中心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

6万

帖子

99

威望

6万

eV
积分
125179

管理员

公社社长

13011#
 楼主 Author| 发表于 Post on 2021-4-13 09:06:13 | 只看该作者 Only view this author
独舞的印迹 发表于 2021-4-12 10:07
卢老师,请问《为什么优化和振动分析不需要用大基组》一文的内容有没有相关的文献可以引用。以及下面这个体 ...

没有什么可引用的,这是专业搞计算的人都知道的基本常识,普遍见于各种主流期刊上的计算文章

基本可以。没必要加弥散。单点算得动的话用def2-TZVP
北京科音自然科学研究中心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

6万

帖子

99

威望

6万

eV
积分
125179

管理员

公社社长

13012#
 楼主 Author| 发表于 Post on 2021-4-13 09:06:50 | 只看该作者 Only view this author
dali 发表于 2021-4-12 16:24
多谢老师!请问老师,gaussian16直接写关键词freq得到的“Sum of electronic and thermal Free Energies= ...

当前温度下的
北京科音自然科学研究中心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

6万

帖子

99

威望

6万

eV
积分
125179

管理员

公社社长

13013#
 楼主 Author| 发表于 Post on 2021-4-13 09:08:32 | 只看该作者 Only view this author
黄舒伟 发表于 2021-4-12 16:37
懂了,这里应该继续用Log文件去算。刚才用一个小分子去验证了,发现:通过比较AAA_b3lyp_opt.chg和 AAA_p ...

本来Multiwfn就不支持半经验计算的波函数
非要用半经验层面的波函数的话,目前Multiwfn只能结合xtb做GFN-xTB计算产生的波函数,参看比如
巨大体系的范德华表面静电势图的快速绘制方法
http://sobereva.com/481http://bbs.keinsci.com/thread-13140-1-1.html
不过GFN-xTB的波函数质量比DFT的有很大差距,不建议基于GFN-xTB波函数算RESP电荷

半经验用于结构的预优化是可以的
北京科音自然科学研究中心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

45

帖子

2

威望

1824

eV
积分
1909

Level 5 (御坂)

13014#
发表于 Post on 2021-4-13 10:30:35 | 只看该作者 Only view this author
本帖最后由 独舞的印迹 于 2021-4-13 10:33 编辑
sobereva 发表于 2021-4-13 09:06
没有什么可引用的,这是专业搞计算的人都知道的基本常识,普遍见于各种主流期刊上的计算文章

基本可以 ...

谢谢老师,
请问算单点能用def2-TZVP是指CHON用def2-TZVP Cu还用sdd吗?(def2-TZVP前四周期不都是全电子基组吗?)
既笨又菜,轻喷。
http://doc.liuyaoze.com

6万

帖子

99

威望

6万

eV
积分
125179

管理员

公社社长

13015#
 楼主 Author| 发表于 Post on 2021-4-13 10:33:35 | 只看该作者 Only view this author
独舞的印迹 发表于 2021-4-13 10:30
谢谢老师,
请问算单点能用def2-TZVP是指给CHON吗?(def2-TZVP前四周期不都是全电子基组吗?)

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

45

帖子

2

威望

1824

eV
积分
1909

Level 5 (御坂)

13016#
发表于 Post on 2021-4-13 10:43:44 | 只看该作者 Only view this author

谢谢老师,
看了您的博文说第四周期用不用赝势基组都可以,那Cu用def2-TZVP是因为比用赝势基组算的好吗?
既笨又菜,轻喷。
http://doc.liuyaoze.com

6万

帖子

99

威望

6万

eV
积分
125179

管理员

公社社长

13017#
 楼主 Author| 发表于 Post on 2021-4-13 11:29:53 | 只看该作者 Only view this author
独舞的印迹 发表于 2021-4-13 10:43
谢谢老师,
看了您的博文说第四周期用不用赝势基组都可以,那Cu用def2-TZVP是因为比用赝势基组算的好吗 ...

不要把赝势基组一概而论,必须明确说具体是什么赝势基组
有什么理由不用def2-TZVP?本来就是质量很好的基组
北京科音自然科学研究中心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

6万

帖子

99

威望

6万

eV
积分
125179

管理员

公社社长

13018#
 楼主 Author| 发表于 Post on 2021-4-13 11:31:08 | 只看该作者 Only view this author
yang2020 发表于 2021-4-12 22:16
尊敬的Sob,各位老师,大家晚上好。
我用Gaussian在不同理论计算级别下得到不同.fch文件,用Multiwfn主功能 ...

计算级别不同结果怎么可能相同

什么级别适合计算这个看
正确地认识分子的能隙(gap)、HOMO和LUMO
http://sobereva.com/543http://bbs.keinsci.com/thread-16758-1-1.html

开壳层体系讨论稳定性时建议用概念密度泛函理论里的软度代替HOMO-LUMO gap说事。否则用alpha还是beta的HOMO-LUMO gap讨论都不合适。
北京科音自然科学研究中心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

212

帖子

0

威望

948

eV
积分
1160

Level 4 (黑子)

13019#
发表于 Post on 2021-4-13 13:05:26 | 只看该作者 Only view this author

谢谢老师,请问如果命令行没写温度关键词,gaussian默认输出的是0K、273.15K、298.15K哪个温度下的热力学数据?

47

帖子

0

威望

206

eV
积分
253

Level 3 能力者

13020#
发表于 Post on 2021-4-13 15:13:47 | 只看该作者 Only view this author
想请教一下

什么是高阶鞍点?

本版积分规则 Credits rule

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

GMT+8, 2026-2-25 16:36 , Processed in 0.198759 second(s), 20 queries , Gzip On.

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