计算化学公社

 找回密码 Forget password
 注册 Register
Views: 762|回复 Reply: 12
打印 Print 上一主题 Last thread 下一主题 Next thread

[新手求助] 请教关于有机分子与锂离子的结合能与吉布斯自由能计算

[复制链接 Copy URL]

21

帖子

0

威望

67

eV
积分
88

Level 2 能力者

各位大佬们好,我是刚接触gaussian计算的小白,有两个问题:
1.最近根据下面的文献使用b3lyp-D3/6-311g+(d,p)的基组进行opt与freq的计算,我再对比其他文献的时候发现,有机分子在嵌入锂离子后,homo与lumo能级均增大,而我的计算结果均是减小,但是能级差确实正常的;
2.其次,我使用上面基组计算freq的结果中的吉布斯根据下面公式计算氧化还原电位的时候,发现结果更是天差地别,如果这种有机分子体系在计算吉布斯自由能时应该选择怎么样的基组?
本体结构的log文件无法上传file:///C:\Users\LENOVO\Documents\Tencent Files\992628774\Image\C2C\`VJ(~CNR$DD(QBESW]XKPXV.png

1li.log

497.63 KB, 下载次数 Times of downloads: 0

Advanced Science - 2023 - Shu - Rational Molecular Design of Redox‐\Active Carbo.pdf

3.42 MB, 下载次数 Times of downloads: 4

112

帖子

0

威望

471

eV
积分
583

Level 4 (黑子)

2#
发表于 Post on 2025-5-23 00:35:44 | 只看该作者 Only view this author
其一,弥散函数不要乱用,参考http://sobereva.com/119
其二,D3校正可以考虑换成D3(BJ)
其三,在B3LYP-D3(BJ)优化的结构基础上,用m062x/def2tzvp做高精度单点能计算
路漫漫其修远兮,吾将上下而求索

6万

帖子

99

威望

5万

eV
积分
120060

管理员

公社社长

3#
发表于 Post on 2025-5-23 07:23:01 | 只看该作者 Only view this author
“本体结构的log文件无法上传”  认真看置顶的社员必读贴,该压缩的必须压缩

6-311g+(d,p)要好好写成6-311+G(d,p)

仔细把此文看了
谈谈重复不出来计算化学文献里的数据的可能原因
http://sobereva.com/678http://bbs.keinsci.com/thread-38911-1-1.html

北京科音自然科学研究中心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

1万

帖子

0

威望

8949

eV
积分
20692

Level 6 (一方通行)

4#
发表于 Post on 2025-5-23 09:19:45 | 只看该作者 Only view this author
你画的地方不是Li的最优配位点吧
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

21

帖子

0

威望

67

eV
积分
88

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2025-5-23 11:21:19 | 只看该作者 Only view this author
gezq 发表于 2025-5-23 00:35
其一,弥散函数不要乱用,参考http://sobereva.com/119
其二,D3校正可以考虑换成D3(BJ)
其三,在B3LYP-D ...

感谢感谢,我把计算参数修改了再试试

21

帖子

0

威望

67

eV
积分
88

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 2025-5-23 11:22:36 | 只看该作者 Only view this author
wzkchem5 发表于 2025-5-23 09:19
你画的地方不是Li的最优配位点吧

我是计算过静电势,按照结果找到的点位,只有羰基附近电势是负的,应该没错吧

21

帖子

0

威望

67

eV
积分
88

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 2025-5-23 11:23:01 | 只看该作者 Only view this author
sobereva 发表于 2025-5-23 07:23
“本体结构的log文件无法上传”  认真看置顶的社员必读贴,该压缩的必须压缩

6-311g+(d,p)要好好写成 ...

好的好的,谢谢sob老师提醒,下次不会再犯了

1万

帖子

0

威望

8949

eV
积分
20692

Level 6 (一方通行)

8#
发表于 Post on 2025-5-24 12:27:34 | 只看该作者 Only view this author
bean233 发表于 2025-5-23 11:22
我是计算过静电势,按照结果找到的点位,只有羰基附近电势是负的,应该没错吧

同等条件下,螯合比单齿配位更稳定,静电势不能判断配位到哪里可以形成螯合作用。再者你的分子也可能先脱质子再结合Li+。
该做构型搜索就做构型搜索,不要以为靠静电势判断就一定能得到正确的结果。

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
longsay + 1 赞!

查看全部评分 View all ratings

Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

26

帖子

0

威望

161

eV
积分
187

Level 3 能力者

9#
发表于 Post on 2025-7-11 23:12:48 | 只看该作者 Only view this author
wzkchem5 发表于 2025-5-24 12:27
同等条件下,螯合比单齿配位更稳定,静电势不能判断配位到哪里可以形成螯合作用。再者你的分子也可能先脱 ...

老师,怎么保证这种又有柔性分子又是团簇的体系找到的是全局最小值点呢

1万

帖子

0

威望

8949

eV
积分
20692

Level 6 (一方通行)

10#
发表于 Post on 2025-7-12 18:46:52 | 只看该作者 Only view this author
QJW 发表于 2025-7-11 23:12
老师,怎么保证这种又有柔性分子又是团簇的体系找到的是全局最小值点呢

不能严格保证,只能让找到全局最小值点的概率无限接近100%。
用ORCA的GOAT和DOCKER模块结合,可以解决这种既需要搜索分子内构象又需要搜索分子间构型的情形。参见https://pubs.acs.org/doi/10.1021/jacs.5c05709?ref=PDF
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

26

帖子

0

威望

161

eV
积分
187

Level 3 能力者

11#
发表于 Post on 2025-7-14 16:06:03 | 只看该作者 Only view this author
wzkchem5 发表于 2025-7-12 18:46
不能严格保证,只能让找到全局最小值点的概率无限接近100%。
用ORCA的GOAT和DOCKER模块结合,可以解决这 ...

老师,请问论文里的GFN2-xTB一定要调参数吗,还有对于溶剂中的反应,怎么判断是用显式还是隐式溶剂化模型呢,辛苦老师!

1万

帖子

0

威望

8949

eV
积分
20692

Level 6 (一方通行)

12#
发表于 Post on 2025-7-14 17:23:50 | 只看该作者 Only view this author
QJW 发表于 2025-7-14 16:06
老师,请问论文里的GFN2-xTB一定要调参数吗,还有对于溶剂中的反应,怎么判断是用显式还是隐式溶剂化模 ...

如果不调参数的话,要做好被问以下两个问题的准备:(1)GFN2-xTB描述溶剂里的锂离子,准确度如何?(2)调参所需时间是否远小于实际计算(包括后续DFT opt freq等)所需时间,如果是,为什么不花少量计算时间调参来改善精度?
如果这两个问题可以回答,就可以不调参。
和溶质的作用远强于和其他溶剂分子作用的溶剂,显式描述,包括但不限于(1)直接和金属配位的溶剂,及(2)和溶质形成强氢键的溶剂(比较弱的氢键不算数)。其他溶剂隐式描述

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
QJW + 5 牛!

查看全部评分 View all ratings

Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

26

帖子

0

威望

161

eV
积分
187

Level 3 能力者

13#
发表于 Post on 2025-7-14 20:49:02 | 只看该作者 Only view this author
wzkchem5 发表于 2025-7-14 17:23
如果不调参数的话,要做好被问以下两个问题的准备:(1)GFN2-xTB描述溶剂里的锂离子,准确度如何?(2) ...

辛苦老师!

本版积分规则 Credits rule

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

GMT+8, 2025-8-12 17:42 , Processed in 1.340999 second(s), 24 queries , Gzip On.

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