|
发布时间: 2023-12-20 09:07
正文摘要:我想计算我这个离子液体和Na+的结合能，主要是为了比较1位置的O和2位置的O哪个更容易与我的Na+结合。我的做法是不考虑PF6-，真空下Na+用def2tzvp，其他原子b3lyp，前后能量做差得到我的结合能，但是目前我计算的结果 ... |
xiaobai1 发表于 2023-12-22 08:45 单个原子可以做振动分析。单原子的平动运动对热力学量也有贡献，振动分析给出的热力学量会考虑到 |
sobereva 发表于 2023-12-22 01:38 好的sob老师，关于def2-TZVP这个基组我再去查一查资料，sob关于您说的Na单个原子没法优化，那么我在算作用能E（Na+）的这部分时是不是只考虑一个单点能就好了呢 |
xiaobai1 发表于 2023-12-21 08:49 Na单个原子根本没法优化 def2-TZVP本来就是全电子基组，用genecp并且试图定义赝势毫无意义 当前是对两个阳离子算的相互作用能，静电互斥那么强，而且你又没带溶剂模型而是算的真空情况，数值为正极其正常 |
sobereva 发表于 2023-12-21 04:12 谢谢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 |
|
"真空下Na+用def2tzvp，其他原子b3lyp" 根本不是并列关系。先分清楚基组和理论方法再说，看下文 计算化学中的一些常见不良写法和用词 http://sobereva.com/298（http://bbs.keinsci.com/thread-1358-1-1.html） 不考虑PF6-时，两个阳离子相互作用能为正极其正常 没具体输入文件，别人无从判断计算有什么可能的硬伤 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2024-11-25 19:29 , Processed in 2.209397 second(s), 26 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites