计算化学公社

标题: 求助计算离子液体和金属阳离子结合能 [打印本页]

作者
Author: xiaobai1 时间: 2023-12-20 09:07
标题: 求助计算离子液体和金属阳离子结合能
我想计算我这个离子液体和Na+的结合能，主要是为了比较1位置的O和2位置的O哪个更容易与我的Na+结合。我的做法是不考虑PF6-，真空下Na+用def2tzvp，其他原子b3lyp，前后能量做差得到我的结合能，但是目前我计算的结果显示E大于0，想请教一下各位老师这个是什么问题，还有就是我这个需要考虑PF6-吗，因为我还有其他带O的离子液体，就是想比较一下各个离子液体与Na+的结合能，不知道加了PF6-会有多大的影响，索性就没加

作者
Author: sobereva 时间: 2023-12-21 04:12
"真空下Na+用def2tzvp，其他原子b3lyp" 根本不是并列关系。先分清楚基组和理论方法再说，看下文
计算化学中的一些常见不良写法和用词
http://sobereva.com/298（http://bbs.keinsci.com/thread-1358-1-1.html）

不考虑PF6-时，两个阳离子相互作用能为正极其正常

没具体输入文件，别人无从判断计算有什么可能的硬伤

作者
Author: xiaobai1 时间: 2023-12-21 08:49

sobereva 发表于 2023-12-21 04:12
"真空下Na+用def2tzvp，其他原子b3lyp" 根本不是并列关系。先分清楚基组和理论方法再说，看下文
计算化学 ...

谢谢sob老师，我应该是在b3lyp泛函下，Na+用def2tzvp，其他原子6-311+G（d，p）然后是在真空下做的加了em=gd3bj，Na+和我离子液体的结合的计算输入我附在了下面（因为篇幅太长我把原子坐标删除了，我是先用b3lyp优化，在M062X下算的单点能
E（A，离子液体）

# opt freq b3lyp/6-311+g(d,p) em=gd3bj

1 1

E（B，Na+）
# opt freq b3lyp/def2tzvp

Title Card Required

1 1
Na -11.12041566 -0.09400025 0.00000000

E(AB）
# opt freq b3lyp/genecp em=gd3bj

Title Card Required

2 1

C H O N 0
6-311+g(d,p)
****
Na 0
def2tzvp
****

Na 0
def2tzvp

作者
Author: sobereva 时间: 2023-12-22 01:38

xiaobai1 发表于 2023-12-21 08:49
谢谢sob老师，我应该是在b3lyp泛函下，Na+用def2tzvp，其他原子6-311+G（d，p）然后是在真空下做的加了em ...

Na单个原子根本没法优化

def2-TZVP本来就是全电子基组，用genecp并且试图定义赝势毫无意义

当前是对两个阳离子算的相互作用能，静电互斥那么强，而且你又没带溶剂模型而是算的真空情况，数值为正极其正常

作者
Author: xiaobai1 时间: 2023-12-22 08:45

sobereva 发表于 2023-12-22 01:38
Na单个原子根本没法优化

def2-TZVP本来就是全电子基组，用genecp并且试图定义赝势毫无意义

好的sob老师，关于def2-TZVP这个基组我再去查一查资料，sob关于您说的Na单个原子没法优化，那么我在算作用能E（Na+）的这部分时是不是只考虑一个单点能就好了呢

作者
Author: sobereva 时间: 2023-12-22 16:03

xiaobai1 发表于 2023-12-22 08:45
好的sob老师，关于def2-TZVP这个基组我再去查一查资料，sob关于您说的Na单个原子没法优化，那么我在算作 ...

单个原子可以做振动分析。单原子的平动运动对热力学量也有贡献，振动分析给出的热力学量会考虑到

欢迎光临计算化学公社 (http://bbs.keinsci.com/)