计算化学公社

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

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

   关闭 [复制链接 Copy URL]

201

帖子

0

威望

3571

eV
积分
3772

Level 5 (御坂)

3001#
发表于 Post on 2017-2-23 13:10:00 | 只看该作者 Only view this author
sobereva 发表于 2017-2-23 08:30
原理上并没有错误,但化学意义不大

谢谢老师回复

175

帖子

0

威望

2136

eV
积分
2311

Level 5 (御坂)

3002#
发表于 Post on 2017-2-24 09:36:19 | 只看该作者 Only view this author
sob老师,我有一个关于genecp的问题想要请教您一下。网上说“Lanl2dz,这是双zeta赝势基组,只能跟6-31G系列的基组联用;如果要得到更精确的结果,可以用Stuggart赝势,这里有三zeta的,可以跟6-311++G**联用”。以Cu(CO)+为例的话,是不是若Cu采用了Lanl2dz基组,则C及O就不能采用6-311+G(d)基组了呢?即混合基组的输入是要求各基组的zeta数一致?
It's my life and it's now or never

6万

帖子

99

威望

6万

eV
积分
125197

管理员

公社社长

3003#
 楼主 Author| 发表于 Post on 2017-2-24 09:47:50 | 只看该作者 Only view this author
wrtgcn 发表于 2017-2-24 09:36
sob老师,我有一个关于genecp的问题想要请教您一下。网上说“Lanl2dz,这是双zeta赝势基组,只能跟6-31G系 ...


网上净瞎说
1 双zeta赝势基组没有必然性必须和双zeta全电子基组联用。如果是配合物,我建议给过渡金属用的赝势基组的zeta数>=给配体用的基组的zeta数。因为过渡金属在当前体系起最关键角色,基组用大一点也不会使计算耗时多增加多少。因此可以比如lanl2TZ(准3-zeta)结合6-31G*。PS:像是lanl08,是完全去收缩的,按网上那个歪理,简直没法跟任何常用的全电子基组一起用了。双zeta全电子基组绝非只有6-31G系列,def2-SVP、cc-pVDZ也是很常用的。
2 6-311++G**是很恶心的,多数情况根本不需要加弥散,尤其是氢的弥散更是几乎总是多余的,白费计算时间。详见
谈谈量子化学中基组的选择
http://sobereva.com/336
诸如Cu(CO)+,加弥散没有丝毫价值。

评分 Rate

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

查看全部评分 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

175

帖子

0

威望

2136

eV
积分
2311

Level 5 (御坂)

3004#
发表于 Post on 2017-2-24 09:57:49 | 只看该作者 Only view this author
sobereva 发表于 2017-2-24 09:47
网上净瞎说
1 双zeta赝势基组没有必然性必须和双zeta全电子基组联用。如果是配合物,我建议给过渡金属 ...

受教了,非常感谢您!
It's my life and it's now or never

108

帖子

0

威望

3260

eV
积分
3368

Level 5 (御坂)

3005#
发表于 Post on 2017-2-24 10:15:41 | 只看该作者 Only view this author
sob老师,我要计算50个原子左右两分子有机体系的相互作用能,希望用不同的方法基组做一下对比试验,分别使用B3LYP和M06-2X两种方法,基组选择为:6-311++G**、ma-def2-tzvp、aug-cc-pvtz。你看这样设计合理么?

6万

帖子

99

威望

6万

eV
积分
125197

管理员

公社社长

3006#
 楼主 Author| 发表于 Post on 2017-2-24 13:02:38 | 只看该作者 Only view this author
ShuangfeiZhu 发表于 2017-2-24 10:15
sob老师,我要计算50个原子左右两分子有机体系的相互作用能,希望用不同的方法基组做一下对比试验,分别使 ...

6-311++G**改成6-311+G**。氢的弥散总是不需要加,而且ma-TZVP里也同样没有氢的弥散。
如果你想检验基组依次增大对结果的影响,6-311+G**、ma-TZVP、aug-cc-pVTZ是合理的。不过中间再插一个jun-cc-pVTZ更好。见
谈谈弥散函数和“月份”基组
http://sobereva.com/119
北京科音自然科学研究中心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

175

帖子

0

威望

2136

eV
积分
2311

Level 5 (御坂)

3007#
发表于 Post on 2017-2-24 14:15:13 | 只看该作者 Only view this author
sob老师,我看了您写得关于赝势下波函数分析的博文,由于个人水平有限还是没能很好的理解。如果我想做第五周期的过渡元素,比如Mo的静电势(ESP)分析,对CHO等元素采用6-31G*,对Mo采用Lanl2dz基组,这样得到的ESP图是否可行?若将Mo的赝势基组换成小核的目的是为了更加准确吗?
It's my life and it's now or never

6万

帖子

99

威望

6万

eV
积分
125197

管理员

公社社长

3008#
 楼主 Author| 发表于 Post on 2017-2-24 14:37:52 | 只看该作者 Only view this author
wrtgcn 发表于 2017-2-24 14:15
sob老师,我看了您写得关于赝势下波函数分析的博文,由于个人水平有限还是没能很好的理解。如果我想做第五 ...

可行。lanl2dz对过渡金属本来就是小核。
北京科音自然科学研究中心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

175

帖子

0

威望

2136

eV
积分
2311

Level 5 (御坂)

3009#
发表于 Post on 2017-2-24 14:52:25 | 只看该作者 Only view this author
本帖最后由 wrtgcn 于 2017-2-24 15:27 编辑

您看我可以这样理解吗:(涉及赝势时)前四周期元素采用全电子基组不定义赝势;四周期后过渡金属采用lanl2dz,主族元素采用def2-SV(P)或cc-pVDZ并定义相应赝势,就可以做常见的波函数分析了;对主族元素也可以先用lanl2dz进行优化,再用cc-pVDZ-pp算一次单点能并进行波函数分析。
It's my life and it's now or never

6万

帖子

99

威望

6万

eV
积分
125197

管理员

公社社长

3010#
 楼主 Author| 发表于 Post on 2017-2-24 15:34:52 | 只看该作者 Only view this author
wrtgcn 发表于 2017-2-24 14:52
您看我可以这样理解吗:(涉及赝势时)前四周期元素采用全电子基组不定义赝势;四周期后过渡金属采用lanl2d ...

不对。def2到了第五周期开始就已经是赝势基组了,cc-pVDZ只对前四周期有定义。
主族特别是前四周期的不要用lanl2dz赝势基组
北京科音自然科学研究中心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

175

帖子

0

威望

2136

eV
积分
2311

Level 5 (御坂)

3011#
发表于 Post on 2017-2-24 15:42:13 | 只看该作者 Only view this author
sobereva 发表于 2017-2-24 15:34
不对。def2到了第五周期开始就已经是赝势基组了,cc-pVDZ只对前四周期有定义。
主族特别是前四周期的不 ...

谢谢您的回答。那对于I元素,我可以先用lanl2dz基组进行优化,再用cc-pVDZ-pp基组算单点能并做ESP图吗?
It's my life and it's now or never

6万

帖子

99

威望

6万

eV
积分
125197

管理员

公社社长

3012#
 楼主 Author| 发表于 Post on 2017-2-24 16:31:50 | 只看该作者 Only view this author
wrtgcn 发表于 2017-2-24 15:42
谢谢您的回答。那对于I元素,我可以先用lanl2dz基组进行优化,再用cc-pVDZ-pp基组算单点能并做ESP图吗?

优化也用cc-pVDZ-PP。lanl2DZ太烂

评分 Rate

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

查看全部评分 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

113

帖子

0

威望

271

eV
积分
384

Level 3 能力者

3013#
发表于 Post on 2017-2-24 19:15:26 | 只看该作者 Only view this author
本帖最后由 erya123 于 2017-2-24 21:41 编辑

内坐标的计算效率是不是就比直角坐标高?我查了查说大体系,对称性差的用cartesian coordinate。也不知道这个体系大有没有经验值,多少原子以上算大?

6万

帖子

99

威望

6万

eV
积分
125197

管理员

公社社长

3014#
 楼主 Author| 发表于 Post on 2017-2-25 01:09:45 | 只看该作者 Only view this author
erya123 发表于 2017-2-24 19:15
内坐标的计算效率是不是就比直角坐标高?我查了查说大体系,对称性差的用cartesian coordinate。也不知道这 ...

“大体系,对称性差的用cartesian coordinate”是胡说,不要听信网上菜鸟谣言
对于几何优化,冗余内坐标是效率最高的,优化所需步数一般最少。这里说过
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164
北京科音自然科学研究中心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

64

帖子

0

威望

737

eV
积分
801

Level 4 (黑子)

3015#
发表于 Post on 2017-2-25 08:17:04 | 只看该作者 Only view this author

过渡态计算出现的问题

老师,您好,我在计算过渡态的时候,结果是正常的,如图1所示,只有一个虚频,但是用优化好的过渡态去做IRC分析,得到如图2的结果,我的IRC命令行如图3所示,我不知道是什么原因导致的?

NOY$0)L7011Z7ED)%~(GIK5.png (34.58 KB, 下载次数 Times of downloads: 29)

NOY$0)L7011Z7ED)%~(GIK5.png

$O5FXH~A~RCRBS79S$X`88D.png (40.57 KB, 下载次数 Times of downloads: 32)

$O5FXH~A~RCRBS79S$X`88D.png

CARXVE_9`YM)IXULC2`8GYN.png (2.47 KB, 下载次数 Times of downloads: 28)

CARXVE_9`YM)IXULC2`8GYN.png

本版积分规则 Credits rule

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

GMT+8, 2026-2-26 18:35 , Processed in 0.226577 second(s), 26 queries , Gzip On.

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