计算化学公社

 找回密码 Forget password
 注册 Register

求助醋酸阴离子和铅离子在DMF溶剂环境中的结合能问题,计算数值比文献偏大

查看数: 140 | 评论数: 13 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2025-8-18 15:01

正文摘要:

我想计算一个醋酸阴离子(-1)和一个Pb离子(+2)的结合能,溶剂是DMF。 首先,我采用# opt freq m062x/def2svp em=gd3 scrf=(smd,solvent=n,n-dimethylformamide)来分别优化了醋酸阴离子和铅离子的结构,没有虚频 ...

回复 Reply

Solitude198 发表于 Post on 3 day ago
Big-Yellow-Dog 发表于 2025-8-18 21:53
大佬!我好像知道了,得在后面定义赝势

genecp是用来使用混合基组的
对于只有Pb离子的情况就可以直接

#p m062x/def2tzvp em=dg3 scrf=(smd,solvent=n,n-DiMethylFormamide)

到醋酸铅的时候再用genecp

看这个 http://sobereva.com/60
Big-Yellow-Dog 发表于 Post on 4 day ago
wzkchem5 发表于 2025-8-18 20:36
不能,因为即使优化结构加了弥散函数,也可以算得动,而且加上肯定会让结果变好,所以没有理由不加。
甚 ...

好的好的,感谢
wzkchem5 发表于 Post on 4 day ago
Big-Yellow-Dog 发表于 2025-8-18 17:21
大佬,还有一个小问题,优化结构时,可以不加弥散函数吗?在后面计算单点能的时候再加上,这样可以吗?:h ...

不能,因为即使优化结构加了弥散函数,也可以算得动,而且加上肯定会让结果变好,所以没有理由不加。
甚至可以说,如果你只算这么小的体系的话,完全应该用sf-X2C处理相对论效应,用赝势都可能被质疑
wzkchem5 发表于 Post on 4 day ago
本帖最后由 wzkchem5 于 2025-8-18 20:36 编辑
Big-Yellow-Dog 发表于 2025-8-18 18:12
大佬们!我现在发现了新的问题,我对CHO用def2tzvpd,函数是从basissets.turbomole网站上粘贴的,对Pb用def ...

其中一个计算没有正确地读入赝势信息
Big-Yellow-Dog 发表于 Post on 4 day ago
wzkchem5 发表于 2025-8-18 17:11
所谓的基组保持一致,指的是:如果反应物里某个原子A用了基组X,那么产物里对应的原子也要用基组X。不是 ...

大佬,还有一个小问题,优化结构时,可以不加弥散函数吗?在后面计算单点能的时候再加上,这样可以吗?
Big-Yellow-Dog 发表于 Post on 4 day ago
wzkchem5 发表于 2025-8-18 17:11
所谓的基组保持一致,指的是:如果反应物里某个原子A用了基组X,那么产物里对应的原子也要用基组X。不是 ...

好的好的,感谢
wzkchem5 发表于 Post on 4 day ago
Big-Yellow-Dog 发表于 2025-8-18 17:04
非常感谢回复,文献上是采用B3LYP 来进行结构优化. 考虑到铅是重原子,采用 LANL2DZ basis set 优化结构, ...

所谓的基组保持一致,指的是:如果反应物里某个原子A用了基组X,那么产物里对应的原子也要用基组X。不是说反应物里所有原子都要用一个基组。
以计算Pb2+ + 2(CH3COO-) -> (CH3COO)2Pb为例,基组保持一致指的是所有计算的Pb必须用同一个基组,所有计算的O必须用同一个基组,但Pb和O可以用不同的基组,也可以一个加弥散而另一个不加
Big-Yellow-Dog 发表于 Post on 4 day ago
pal 发表于 2025-8-18 16:51
计算能量保持基组一致,要跟文献里的结果对比,选择对应的方法和基组

非常感谢回复,文献上是采用B3LYP 来进行结构优化. 考虑到铅是重原子,采用 LANL2DZ basis set 优化结构, 这个 def2-TZVP basis set 来进行单点能计算。文献里没有考虑溶剂环境,算的是气相,可是实际情况就是在溶剂里面的,我感觉还是得考虑溶剂环境的。


我用# m062x/def2tzvp pop=(nbo) empiricaldispersion=gd3 scrf=(smd,solvent=n,n-dimethylformamide) 重新算了一下,没加弥散函数,阴离子,阳离子,复合物都没加弥散函数,能量计算保持了方法一致,最后算出来结果是-2.17eV,这个感觉还靠点谱,溶剂影响下,结合能就是会减弱。但是计算阴离子单点能不加弥散函数能行吗?


现在问题在阳离子上,Pb2+要是采用弥散函数,能量偏差很大,要是不加弥散函数,阴离子计算单点能又需要加弥散函数,该怎么统一呢?
pal 发表于 Post on 4 day ago
计算能量保持基组一致,要跟文献里的结果对比,选择对应的方法和基组

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

GMT+8, 2025-8-22 13:22 , Processed in 0.181140 second(s), 25 queries , Gzip On.

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