请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 182|回复: 5

[Gaussian/gview] 萌新求解:关于BSSE的一个小问题

[复制链接]

12

帖子

0

威望

122

eV
积分
134

Level 2 能力者

发表于 2018-6-12 20:18:56 | 显示全部楼层 |阅读模式
想把DFTB+优化得到的结构用高斯B3LYP-D3BJ算BSSE,是直接用这个结构拿来算单点呢,还是需要把这个DFTB+得到的结构用高斯再优化一遍,再去算单点呢?

1万

帖子

25

威望

1万

eV
积分
31827

管理员

公社社长

发表于 2018-6-13 05:40:20 | 显示全部楼层
说明算的是什么体系、多大、目的何在

对于Gaussian做DFT直接算得动的体系,用DFTB+是多余的,就算你先用DFTB+优化过,为了更准确的结果,也建议在Gaussian里用像样的级别重新优化
仅对于计算相互作用能而言才需要考虑BSSE问题,用的结构肯定也是优化过的结构。若是Gaussian做得动counterpoise的情况,Gaussian也不会优化不动(至少降低级别后算得动),显然还是应当在Gaussian里再优化
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位。

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

12

帖子

0

威望

122

eV
积分
134

Level 2 能力者

 楼主| 发表于 2018-6-13 09:35:33 | 显示全部楼层
sobereva 发表于 2018-6-13 05:40
说明算的是什么体系、多大、目的何在

对于Gaussian做DFT直接算得动的体系,用DFTB+是多余的,就算你先用 ...

一共的原子数为164,DFTB+里有一个限制某些原子仅在某个方向运动的功能,我想计算一下在存在特定限制的条件下,一个小分子在一个碳管类的物体中穿越的过程。那么在运动的过程中会得到一个能量最低的一步。我想直接计算这一步结构在高斯里的相互作用能。我也试过再用高斯优化,但是它的结构产生了些变化。
所以我在犹豫,既然我想特定的计算DFTB+中产生结果的相互作用能,可不可以就用这个结构直接在高斯里算相互作用能?(虽然再用高斯优化会得到更准确的结果,但是这个结构却不是我计算的特定条件下的结构)

1万

帖子

25

威望

1万

eV
积分
31827

管理员

公社社长

发表于 2018-6-13 09:49:51 | 显示全部楼层
lira 发表于 2018-6-13 09:35
一共的原子数为164,DFTB+里有一个限制某些原子仅在某个方向运动的功能,我想计算一下在存在特定限制的条 ...

凑合可以接受,算是有一定正当理由
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位。

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

12

帖子

0

威望

122

eV
积分
134

Level 2 能力者

 楼主| 发表于 2018-6-13 14:45:29 | 显示全部楼层
sobereva 发表于 2018-6-13 09:49
凑合可以接受,算是有一定正当理由

好的,谢谢老师的指点。还有一个小问题,我想问一下,比如B3LYP泛函下使用的6-31g(d)和6-31+g(d),这个+号,是添加弥散作用的意思么?是不是单独优化一个小分子用不带+号的就可以,计算两个小分子相互作用的时候用带+号的?

1万

帖子

25

威望

1万

eV
积分
31827

管理员

公社社长

发表于 2018-6-13 19:12:37 | 显示全部楼层
lira 发表于 2018-6-13 14:45
好的,谢谢老师的指点。还有一个小问题,我想问一下,比如B3LYP泛函下使用的6-31g(d)和6-31+g(d),这个+ ...


仔细看
谈谈弥散函数和“月份”基组
http://sobereva.com/119
乱谈DFT-D
http://sobereva.com/83
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位。

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2018-6-25 10:42 , Processed in 0.181105 second(s), 23 queries .

快速回复 返回顶部 返回列表