计算化学公社

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

[Gaussian/gview] 求助计算H2O和Zn的结合能的问题

[复制链接 Copy URL]

63

帖子

0

威望

311

eV
积分
374

Level 3 能力者

请问各位老师,我在计算H2O和Zn的结合能时,先通过molclus确定了Zn-6H2O的初始构型后进行优化,我的想法是根据 EAB-EA-EB来确定结合能的大小,这里EAB我认为就是Zn-6H2O,然后算水的时候,我是重新生成了6个水的团簇构型并将优化后的作为初始构型,但是从优化的结果显示,6H2O的构型和Zn-6H2O中的水构型产生了很大的区别,不知道这样的优化结果是不是可信

H2O.png (23.54 KB, 下载次数 Times of downloads: 11)

H2O.png

Zn-H2O.png (38.11 KB, 下载次数 Times of downloads: 8)

Zn-H2O.png

1万

帖子

0

威望

9868

eV
积分
22108

Level 6 (一方通行)

2#
发表于 Post on 2022-11-29 17:16:47 | 只看该作者 Only view this author
本来就应该有很大区别,因为Zn和H2O的相互作用会改变H2O分子的取向
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?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

63

帖子

0

威望

311

eV
积分
374

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2022-11-29 19:38:55 | 只看该作者 Only view this author
wzkchem5 发表于 2022-11-29 17:16
本来就应该有很大区别,因为Zn和H2O的相互作用会改变H2O分子的取向

嗯嗯好的老师,我现在的做法是水溶液体系里优化并得到各部分的自由能校正项,在def2TZVP基组下(也在水溶液体系里)计算各部分的单点能,然后在相加作差,老师我得第一个问题是这里的Zn是不是得在真空环境下计算,然后去算一个溶解自由能呢?其次我在算Zn和其他有机分子的结合能的时候,是不是拿掉1个H2O换成我的有机分子用同样的方法去得到结合能呢?谢谢老师

1万

帖子

0

威望

9868

eV
积分
22108

Level 6 (一方通行)

4#
发表于 Post on 2022-11-29 20:30:11 | 只看该作者 Only view this author
xiaobai1 发表于 2022-11-29 12:38
嗯嗯好的老师,我现在的做法是水溶液体系里优化并得到各部分的自由能校正项,在def2TZVP基组下(也在水溶 ...

首先,你要算的是Zn2+还是Zn原子?
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?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

1016

帖子

0

威望

2126

eV
积分
3142

Level 5 (御坂)

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

5#
发表于 Post on 2022-11-29 21:08:31 | 只看该作者 Only view this author
wzkchem5 发表于 2022-11-29 17:16
本来就应该有很大区别,因为Zn和H2O的相互作用会改变H2O分子的取向

这个不太好处理啊,Zn2+的六水合物未必对称性很高,水分子之间说不定都有氢键,说不清是Zn2+和水分子的作用还是水分子之间的作用。
ORCA大法好!CP2K大法好!
嶺南粤音 泛粤典 https://jyutjam.org/

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

12

帖子

0

威望

1319

eV
积分
1331

Level 4 (黑子)

6#
发表于 Post on 2022-11-30 00:18:16 | 只看该作者 Only view this author
根据 EAB-EA-EB来确定结合能的大小的算法没错,但是这里的EAB,EA,EB的定义是整个团簇下优化后整体和某个片段的单点能,EA,EB不是各自独立时候优化后的单点能,结合能的计算可以参考社长的http://sobereva.com/542

63

帖子

0

威望

311

eV
积分
374

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2022-11-30 08:13:40 | 只看该作者 Only view this author
wzkchem5 发表于 2022-11-29 20:30
首先,你要算的是Zn2+还是Zn原子?

Zn2+

1万

帖子

0

威望

9868

eV
积分
22108

Level 6 (一方通行)

8#
发表于 Post on 2022-11-30 15:04:04 | 只看该作者 Only view this author

那一般只能算溶解自由能,或者气相里面的结合能。参见http://bbs.keinsci.com/forum.php ... =%BD%E1%BA%CF%C4%DC
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?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

63

帖子

0

威望

311

eV
积分
374

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2022-11-30 16:27:08 | 只看该作者 Only view this author
wzkchem5 发表于 2022-11-30 15:04
那一般只能算溶解自由能,或者气相里面的结合能。参见http://bbs.keinsci.com/forum.php?mod=viewthread& ...

老师您好,我大概看了下那篇帖子,我也是想比较我的几个有机小分子和Zn2+的结合能力强弱,帖子里提到了用真空下的结合能作为评判标准,这里的话我的单个Zn+还需要opt+freq来得到自由能校正项吗,还是直接在高精度下计算单点能作为E(Zn2+)就可以了

1万

帖子

0

威望

9868

eV
积分
22108

Level 6 (一方通行)

10#
发表于 Post on 2022-11-30 17:18:18 | 只看该作者 Only view this author
xiaobai1 发表于 2022-11-30 09:27
老师您好,我大概看了下那篇帖子,我也是想比较我的几个有机小分子和Zn2+的结合能力强弱,帖子里提到了用 ...

不需要opt。自由能校正项需要算,因为即使是单原子离子,自由能校正项也不是0,只不过转动、振动贡献为0而已。
当然如果只关心结合能的相对大小,那么Zn2+这边即使不加自由能校正项也不影响结论,但问题是算出来的量就不能叫结合自由能了。
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?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

63

帖子

0

威望

311

eV
积分
374

Level 3 能力者

11#
 楼主 Author| 发表于 Post on 2022-12-1 08:47:49 | 只看该作者 Only view this author
本帖最后由 xiaobai1 于 2022-12-1 08:51 编辑
wzkchem5 发表于 2022-11-30 17:18
不需要opt。自由能校正项需要算,因为即使是单原子离子,自由能校正项也不是0,只不过转动、振动贡献为0 ...

好的老师,那再计算高精度电子能的时候还需要用混合基组吗?我现在算高精度电子能的时候,对一般的C H O用的是M062X/6-311++g(d,p),Zn的话用了M062X/def2TZVP,不知道这样是否可行

1万

帖子

0

威望

9868

eV
积分
22108

Level 6 (一方通行)

12#
发表于 Post on 2022-12-1 15:09:03 | 只看该作者 Only view this author
xiaobai1 发表于 2022-12-1 01:47
好的老师,那再计算高精度电子能的时候还需要用混合基组吗?我现在算高精度电子能的时候,对一般的C H O ...

可以,但注意不要加不必要的弥散函数
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?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

63

帖子

0

威望

311

eV
积分
374

Level 3 能力者

13#
 楼主 Author| 发表于 Post on 2022-12-2 10:14:06 | 只看该作者 Only view this author
wzkchem5 发表于 2022-12-1 15:09
可以,但注意不要加不必要的弥散函数

好的 谢谢老师

本版积分规则 Credits rule

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

GMT+8, 2026-2-22 00:49 , Processed in 0.172265 second(s), 23 queries , Gzip On.

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