计算化学公社

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

[CASTEP/Dmol3/MS] DMol3计算电催化过程自由能的问题

[复制链接 Copy URL]

183

帖子

2

威望

4981

eV
积分
5204

Level 6 (一方通行)

能科top50

本帖最后由 zorow 于 2019-12-27 11:12 编辑

各位老师好,本人目前在计算电催化过程的阶梯图。用DMol3完成了几何优化和频率计算,但计算自由能的时候有一点疑惑:
DMol3的输出文件中有不同温度下的自由能和熵(图1),按照DMol3的教程,自由能=E(电子总能)+G(输出文件中的自由能校正量);而按照公式G=E+ZPE-TS,我需要通过频率来计算电子总能、ZPE和TS(如图2)。请问这两种方法哪个更好?采用第二种方法的时候,通过tool-vibrational analysis得出的频率很多,一些为负的是不是就被排除了?另外一些强度很小的频率,是否也可以排除?谢谢大家!


微信截图_20191227105555.png (18.15 KB, 下载次数 Times of downloads: 76)

图1 DMol3输出文件中的G S

图1 DMol3输出文件中的G S

微信图片_20191226204435.png (18.14 KB, 下载次数 Times of downloads: 52)

图2 ZPE和TS的计算公式

图2 ZPE和TS的计算公式
In defeat, malice. In victory, revenge.

3809

帖子

3

威望

1万

eV
积分
20339

Level 6 (一方通行)

围观吃瓜群众

2#
发表于 Post on 2020-1-2 06:01:32 | 只看该作者 Only view this author
在vibrational analysis看没虚频的情况下,直接看G就行。需要注意的是,分子在不放到格子处理的时候,程序会处理平动和转动,自己加很容易忽略这些项。vibrational analysis应该并不是直接略去负的和强度小的模式。应该是根据特定的投影规则去除整体运动的贡献。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

183

帖子

2

威望

4981

eV
积分
5204

Level 6 (一方通行)

能科top50

3#
 楼主 Author| 发表于 Post on 2020-1-2 10:23:49 | 只看该作者 Only view this author
卡开发发 发表于 2020-1-2 06:01
在vibrational analysis看没虚频的情况下,直接看G就行。需要注意的是,分子在不放到格子处理的时候,程序 ...

感谢卡开发发老师的回复,新年快乐!“分子在不放到格子处理的时候”指的是吸附前气相分子的自由能吗?这个我其实准备查表处理的。
In defeat, malice. In victory, revenge.

3809

帖子

3

威望

1万

eV
积分
20339

Level 6 (一方通行)

围观吃瓜群众

4#
发表于 Post on 2020-1-2 13:04:48 | 只看该作者 Only view this author
zorow 发表于 2020-1-2 10:23
感谢卡开发发老师的回复,新年快乐!“分子在不放到格子处理的时候”指的是吸附前气相分子的自由能吗?这 ...

对,但建议即便这样最好别把分子放格子里面算。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

3809

帖子

3

威望

1万

eV
积分
20339

Level 6 (一方通行)

围观吃瓜群众

5#
发表于 Post on 2020-1-2 13:05:03 | 只看该作者 Only view this author
zorow 发表于 2020-1-2 10:23
感谢卡开发发老师的回复,新年快乐!“分子在不放到格子处理的时候”指的是吸附前气相分子的自由能吗?这 ...

对,但建议即便这样最好别把分子放格子里面算。

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
zorow + 5 谢谢

查看全部评分 View all ratings

日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

183

帖子

2

威望

4981

eV
积分
5204

Level 6 (一方通行)

能科top50

6#
 楼主 Author| 发表于 Post on 2020-1-2 14:56:44 | 只看该作者 Only view this author
卡开发发 发表于 2020-1-2 13:05
对,但建议即便这样最好别把分子放格子里面算。

好的老师,感谢您的指教!
In defeat, malice. In victory, revenge.

23

帖子

0

威望

147

eV
积分
170

Level 3 能力者

7#
发表于 Post on 2020-5-24 09:19:27 | 只看该作者 Only view this author
你好,你算自由能的时候,没有算这个过渡态吧

4

帖子

0

威望

17

eV
积分
21

Level 1 能力者

8#
发表于 Post on 2022-1-7 16:48:35 | 只看该作者 Only view this author
你好,我也遇到了相同的问题,也是电催化反应,用CASTEP模块计算了E,然后如何进行校正到G呢?我的理解是楼主你用的是dmol3计算出频率,固定住其他原子,除了吸附物,得到的要校正的G值对吗

17

帖子

0

威望

5938

eV
积分
5955

Level 6 (一方通行)

9#
发表于 Post on 2022-5-31 10:45:23 | 只看该作者 Only view this author
请问楼主,计算电催化的反应自由能变化的时候,液态水分子的自由能以及氢氧根等物质的自由能是怎么计算的?

本版积分规则 Credits rule

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

GMT+8, 2026-2-20 06:43 , Processed in 0.211032 second(s), 32 queries , Gzip On.

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