计算化学公社

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

[Gaussian/gview] 判断一个物种是否可以形成并稳定存在的判据求助

[复制链接 Copy URL]

226

帖子

0

威望

878

eV
积分
1104

Level 4 (黑子)

各位老师好,如何判断一个物种可以在体系中存在。如:金属离子在水溶液中以[M(H2O)6]2+形式存在,当加入CL-时,随浓度增加,逐级取代水分子。只要每级反应的吉布斯自由能为负就能说明存在吗。

1万

帖子

0

威望

9002

eV
积分
20758

Level 6 (一方通行)

2#
发表于 Post on 2022-4-8 14:19:57 | 只看该作者 Only view this author
本帖最后由 wzkchem5 于 2022-4-8 07:22 编辑

对,但是注意计算出来的自由能变是标准自由能,需要用实验的各物种浓度换算成实际的自由能,再判断是否为负。
此外众所周知,没有严格不存在的物质,不管多么不稳定的物质都会有极少量存在,只不过可能测不出来而已。所以在判断的时候,还需要根据你的实验手段或者文献惯例,人为定义一个浓度基准,低于这个浓度就认为测不出来,这样才能判断一个物质是否能测出来存在。
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
本团队长期招收研究生,有意者可私信联系

226

帖子

0

威望

878

eV
积分
1104

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2022-4-9 13:13:26 | 只看该作者 Only view this author
本帖最后由 sun666 于 2022-4-9 20:52 编辑

游客,本帖隐藏的内容需要积分高于 500 才可浏览,您当前积分为 0

1万

帖子

0

威望

9002

eV
积分
20758

Level 6 (一方通行)

4#
发表于 Post on 2022-4-9 14:55:00 | 只看该作者 Only view this author
sun666 发表于 2022-4-9 06:13
谢谢老师指导。我目前遇到的就是您说的检测不到但计算结果自由能为负的问题。具体如下:CL-逐级取代[Ni(H ...

你理论计算的级别是什么?溶剂效应是怎么考虑的?每取代一个氯离子自由能具体负多少?
标准自由能换算成实验浓度下的自由能,一般物化教材都会教,加一个类似于RTlnc的项就行了
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
本团队长期招收研究生,有意者可私信联系

226

帖子

0

威望

878

eV
积分
1104

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2022-4-9 20:47:13 | 只看该作者 Only view this author
wzkchem5 发表于 2022-4-9 14:55
你理论计算的级别是什么?溶剂效应是怎么考虑的?每取代一个氯离子自由能具体负多少?
标准自由能换算成 ...

游客,本帖隐藏的内容需要积分高于 1500 才可浏览,您当前积分为 0

HL[OX9XTEED_OM$[}6YLPUP.png (22.16 KB, 下载次数 Times of downloads: 23)

HL[OX9XTEED_OM$[}6YLPUP.png

1万

帖子

0

威望

9002

eV
积分
20758

Level 6 (一方通行)

6#
发表于 Post on 2022-4-9 22:35:31 | 只看该作者 Only view this author
sun666 发表于 2022-4-9 13:47
**** 本内容被作者隐藏 ****

首先检查一下自旋多重度有没有设对。其次可以尝试对Cl-做显式溶剂化
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
本团队长期招收研究生,有意者可私信联系

226

帖子

0

威望

878

eV
积分
1104

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2022-4-12 10:18:38 | 只看该作者 Only view this author
本帖最后由 sun666 于 2022-4-12 10:30 编辑
wzkchem5 发表于 2022-4-9 22:35
首先检查一下自旋多重度有没有设对。其次可以尝试对Cl-做显式溶剂化

游客,本帖隐藏的内容需要积分高于 500 才可浏览,您当前积分为 0

IDR`3UI%VRM$I2L~}UM](QS.png (48.04 KB, 下载次数 Times of downloads: 32)

IDR`3UI%VRM$I2L~}UM](QS.png

ener.xlsx

11.63 KB, 下载次数 Times of downloads: 1

1万

帖子

0

威望

9002

eV
积分
20758

Level 6 (一方通行)

8#
发表于 Post on 2022-4-12 15:31:19 | 只看该作者 Only view this author
sun666 发表于 2022-4-12 03:18
**** 本内容被作者隐藏 ****

不是的,注意你这里计算的是 Cl-(用SMD描述溶剂) + nH2O(用SMD描述溶剂) -> Cl(H2O)n-(用SMD描述溶剂) 的自由能变化,所以你算的不是溶解自由能,而是SMD溶解自由能和显式溶剂化的溶解自由能之差。所以通俗但不太严格地说,这个结果不是证明Cl-不能发生溶剂化,而是证明SMD低估了Cl-(aq)的能量。(之所以这个说法不完全严格,是因为你的显式溶剂化计算没有做构象平均,做了构象平均以后刚才这句话才严格成立。)
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
本团队长期招收研究生,有意者可私信联系

1016

帖子

0

威望

2096

eV
积分
3112

Level 5 (御坂)

天上地下三山亖海五湖六合八荒九州之泛粤大典編輯委員會編外委員

9#
发表于 Post on 2022-4-12 15:34:52 | 只看该作者 Only view this author
wzkchem5 发表于 2022-4-12 15:31
不是的,注意你这里计算的是 Cl-(用SMD描述溶剂) + nH2O(用SMD描述溶剂) -> Cl(H2O)n-(用SMD描述溶剂) 的 ...

我不太明白这里,理论上这一步delta G 应该为0才对吧
ORCA大法好!CP2K大法好!
嶺南粤音 泛粤典 https://jyutjam.org/

山河自落蒼梧月,風雨猶驅草木兵。
san huɔɞ ki lɔk tʰɔŋ ŋ ᵑgut,fʊŋ ʝi ʝiu kʰui tʰɔu ᵐbʊk pɪŋ

1万

帖子

0

威望

9002

eV
积分
20758

Level 6 (一方通行)

10#
发表于 Post on 2022-4-12 15:38:33 | 只看该作者 Only view this author
chands 发表于 2022-4-12 08:34
我不太明白这里,理论上这一步delta G 应该为0才对吧

对(但是需要做构象平均才能严格为0)。所以这个deltaG不是0,说明要么是SMD不准,要么是显式溶剂化计算不准,一般前者概率大
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
本团队长期招收研究生,有意者可私信联系

1016

帖子

0

威望

2096

eV
积分
3112

Level 5 (御坂)

天上地下三山亖海五湖六合八荒九州之泛粤大典編輯委員會編外委員

11#
发表于 Post on 2022-4-12 17:23:14 | 只看该作者 Only view this author
wzkchem5 发表于 2022-4-12 15:38
对(但是需要做构象平均才能严格为0)。所以这个deltaG不是0,说明要么是SMD不准,要么是显式溶剂化计算 ...

谢谢!
ORCA大法好!CP2K大法好!
嶺南粤音 泛粤典 https://jyutjam.org/

山河自落蒼梧月,風雨猶驅草木兵。
san huɔɞ ki lɔk tʰɔŋ ŋ ᵑgut,fʊŋ ʝi ʝiu kʰui tʰɔu ᵐbʊk pɪŋ

226

帖子

0

威望

878

eV
积分
1104

Level 4 (黑子)

12#
 楼主 Author| 发表于 Post on 2022-4-12 18:18:47 | 只看该作者 Only view this author
本帖最后由 sun666 于 2022-4-12 18:20 编辑
wzkchem5 发表于 2022-4-12 15:31
不是的,注意你这里计算的是 Cl-(用SMD描述溶剂) + nH2O(用SMD描述溶剂) -> Cl(H2O)n-(用SMD描述溶剂) 的 ...

老师 我的理解是,对CL-进行显示溶剂化,是要对Cl-在水溶液中与H2O的配位个数进行探究。因此计算了CL-+nH2O=CL(H2O)n-的吉布斯自由能,弱为负则证明可以形成CL(H2O)n-,以此来判断CL-的配位水个数。excel中的▲G(smd)是通过热力学循环得到的CL-+nH2O=CL(H2O)n-反应自由能

1万

帖子

0

威望

9002

eV
积分
20758

Level 6 (一方通行)

13#
发表于 Post on 2022-4-12 18:51:27 | 只看该作者 Only view this author
sun666 发表于 2022-4-12 11:18
老师 我的理解是,对CL-进行显示溶剂化,是要对Cl-在水溶液中与H2O的配位个数进行探究。因此计算了CL-+nH ...

但是这里的问题在于Cl-+nH2O=Cl(H2O)n-里的Cl-必须是气相的Cl-。如果是溶剂里的Cl-,那么即使你不放水分子在周围,你算的也是和水有配位的Cl-。那么为了计算理论级别统一起见,H2O和Cl(H2O)n-也必须在气相里面算。这样固然忽略了第二配位层以及以外的水对配位数的贡献,但是这是不得已而为之的近似。
更加严格的方法是做一个水盒子,在里面放一个Cl-跑动力学,统计第一配位层的水分子数目。不过因为MD需要跑很长时间,水盒子也需要比较大,所以只能用低精度方法,有一定概率得不偿失。
不过对于你这个特定的体系,Cl(H2O)n-显然是很早就被研究透了的体系,所以你需要做的是去查文献,查一下水溶液里Cl-的第一配位层是什么样的,水分子的位置和取向是怎样的。而不要从头自己去探究这个问题。
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
本团队长期招收研究生,有意者可私信联系

226

帖子

0

威望

878

eV
积分
1104

Level 4 (黑子)

14#
 楼主 Author| 发表于 Post on 2022-4-12 19:03:43 | 只看该作者 Only view this author
本帖最后由 sun666 于 2022-4-12 19:44 编辑
wzkchem5 发表于 2022-4-12 18:51
但是这里的问题在于Cl-+nH2O=Cl(H2O)n-里的Cl-必须是气相的Cl-。如果是溶剂里的Cl-,那么即使你不放水分 ...

了解了,感谢老师。但是还是有疑问,那计算金属离子的配位水个数时也不能按着热力学循环计算吉布斯自由能来确定吗(因为即便计算的是M2+单个金属离子,在使用SMD时也会有配位水的作用,是这个意思吗)

1万

帖子

0

威望

9002

eV
积分
20758

Level 6 (一方通行)

15#
发表于 Post on 2022-4-12 20:16:36 | 只看该作者 Only view this author
sun666 发表于 2022-4-12 12:03
了解了,感谢老师。但是还是有疑问,那计算金属离子的配位水个数时也不能按着热力学循环计算吉布斯自由能 ...

不是说不能用热力学循环,而是不能直接计算溶液里的配位数,要么用气相团簇里的配位数近似代替溶剂里的配位数,要么加特别多的水分子跑MD确定
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
本团队长期招收研究生,有意者可私信联系

本版积分规则 Credits rule

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

GMT+8, 2025-8-17 07:11 , Processed in 0.195291 second(s), 23 queries , Gzip On.

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