计算化学公社

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

[Gaussian/gview] gaussian如何计算EA电子亲和势和IP电离能,求详细步骤

  [复制链接 Copy URL]

195

帖子

0

威望

699

eV
积分
894

Level 4 (黑子)

请问大家和sob老师,如果使用高斯计算EA和IP呀,求详细步骤。多谢大家。

评分 Rate

参与人数
Participants 1
eV +3 收起 理由
Reason
reddy + 3

查看全部评分 View all ratings

1万

帖子

0

威望

9857

eV
积分
22093

Level 6 (一方通行)

52#
发表于 Post on 2024-3-24 00:30:17 | 只看该作者 Only view this author
WZH1233 发表于 2024-3-23 16:07
sob老师,请教两种分子片段计结合后的VIP数值,一定处于两个分子单独的VIP数值中间吗?最近重复文献VIP(C ...

不一定。举例:一个HOMO比较高的中性分子,和一个HOMO比较低的阴离子结合,阴离子可能会因为库伦作用,把中性分子的HOMO推高,导致总体系的VIP比两个分子单独的VIP都要小。
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
本团队长期招收研究生,有意者可私信联系

20

帖子

0

威望

484

eV
积分
504

Level 4 (黑子)

51#
发表于 Post on 2024-3-23 23:07:26 | 只看该作者 Only view this author
本帖最后由 WZH1233 于 2024-3-23 23:12 编辑
sobereva 发表于 2017-12-20 13:34
你补的是氢,一个就够

sob老师,请教两种分子片段计结合后的VIP数值,一定处于两个分子单独的VIP数值中间吗?最近重复文献VIP(CNT)=4.05; VIP(PANI)=4.44;  VIP(CNT+PANI)=4.34,CNT和PANI的趋势以及结合后的数据趋势均和参考文献有差异,想请教下从哪方面查找原因呢?计算方参考的文献细节


202403232309005229..png (853 Bytes, 下载次数 Times of downloads: 32)

202403232309005229..png

33

帖子

0

威望

367

eV
积分
400

Level 3 能力者

50#
发表于 Post on 2024-3-4 20:02:30 | 只看该作者 Only view this author
wzkchem5 发表于 2024-2-27 22:23
这不叫EA,只有得到一个电子释放的能量才叫EA。得到两个电子释放的能量,可以拆成得到第一个电子释放的能 ...

好的 谢谢老师回复 概念理清楚啦

1万

帖子

0

威望

9857

eV
积分
22093

Level 6 (一方通行)

49#
发表于 Post on 2024-2-27 22:23:00 | 只看该作者 Only view this author
cabbage 发表于 2024-2-27 14:31
sob老师您好,如果我想要算EA,我的体系是中性的分子变成带两个负电荷的阴离子,是用中性的能量减去阴离 ...

这不叫EA,只有得到一个电子释放的能量才叫EA。得到两个电子释放的能量,可以拆成得到第一个电子释放的能量和得到第二个电子释放的能量,前者是中性分子的EA,后者是一价阴离子的EA,但两者的和不叫EA
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
本团队长期招收研究生,有意者可私信联系

1万

帖子

0

威望

9857

eV
积分
22093

Level 6 (一方通行)

48#
发表于 Post on 2024-2-27 22:21:04 | 只看该作者 Only view this author
syj0922 发表于 2023-3-21 01:53
sob老师,分析A和B分子(一个输入文件)的电子亲合能时,如果只让A得到电子,B分子保持不变,是不可以采用在 ...

不能。固定原子指的是让原子的位置不动,不是让原子上面的电子、轨道等也不动。
视研究目的,多数情况下通过调节初猜轨道占据数使得收敛到你想要的电子组态来解决,少数情况下用CDFT解决。什么时候用哪种方法见http://sobereva.com/271
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
本团队长期招收研究生,有意者可私信联系

33

帖子

0

威望

367

eV
积分
400

Level 3 能力者

47#
发表于 Post on 2024-2-27 21:31:51 | 只看该作者 Only view this author

sob老师您好,如果我想要算EA,我的体系是中性的分子变成带两个负电荷的阴离子,是用中性的能量减去阴离子的能量吗?感谢老师的解答

4

帖子

0

威望

313

eV
积分
317

Level 3 能力者

46#
发表于 Post on 2023-3-21 08:53:02 | 只看该作者 Only view this author
本帖最后由 syj0922 于 2023-3-21 08:54 编辑

sob老师,分析A和B分子(一个输入文件)的电子亲合能时,如果只让A得到电子,B分子保持不变,是不可以采用在输入文件中固定B分子的所有原子呢?

32

帖子

0

威望

375

eV
积分
407

Level 3 能力者

45#
发表于 Post on 2023-3-20 20:53:53 | 只看该作者 Only view this author
wzkchem5 发表于 2023-3-19 00:05
仔细阅读http://sobereva.com/543

感谢老师的指导,该问题已解决。

1万

帖子

0

威望

9857

eV
积分
22093

Level 6 (一方通行)

44#
发表于 Post on 2023-3-19 00:05:12 | 只看该作者 Only view this author
汪文焕 发表于 2023-3-18 13:43
非常感谢老师的解答,目前关于EA的计算方法,检索到两种。
一种是EA=E(N)-E(N-1),关于e的能量未加以讨 ...

仔细阅读http://sobereva.com/543
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
本团队长期招收研究生,有意者可私信联系

32

帖子

0

威望

375

eV
积分
407

Level 3 能力者

43#
发表于 Post on 2023-3-18 20:43:05 | 只看该作者 Only view this author
wzkchem5 发表于 2023-3-14 22:59
当体系的中性态和阴离子态的第一激发能都远大于kT的时候,按E(N)+E(e-)-E(N+1)算没问题。
但是假如中性 ...

非常感谢老师的解答,目前关于EA的计算方法,检索到两种。
一种是EA=E(N)-E(N-1),关于e的能量未加以讨论(可能影响较小都忽略了)。学生想问:N和N+1态间电子 能量相减,这一数据是Etot吗

另一种是EA=-ELUMO,学生不太理解这两种方法之间的区别与联系。

1万

帖子

0

威望

9857

eV
积分
22093

Level 6 (一方通行)

42#
发表于 Post on 2023-3-14 22:59:24 | 只看该作者 Only view this author
汪文焕 发表于 2023-3-14 15:00
老师您好,请问如果想要计算298K的EA,是不是EA=E(N)+E(e-)-E(N+1)
E(e-)有298K下的实验数据

当体系的中性态和阴离子态的第一激发能都远大于kT的时候,按E(N)+E(e-)-E(N+1)算没问题。
但是假如中性态和阴离子态之中至少有一个的第一激发能和kT可比的话,按E(N)+E(e-)-E(N+1)算可能会有一个问题。因为E(e-)之所以会非零,是因为电子占据了较高的平动能级。另一方面我们知道,一般计算分子在有限温度下的电子能,是假设电子不会被激发到激发态。那么问题来了,E(N)、E(N+1)是按电子永远处于基态来算的,E(e-)却考虑了电子被激发的贡献,可能会造成理论级别不统一。至于这个不统一是否会有实质性的影响,乃至于是否会导致E(N)+E(e-)-E(N+1)这个公式不可用(或者必须明确考虑了激发态对E(N)、E(N+1)的贡献后才可用),我就不太清楚了,还请其他人补充。
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
本团队长期招收研究生,有意者可私信联系

32

帖子

0

威望

375

eV
积分
407

Level 3 能力者

41#
发表于 Post on 2023-3-14 22:00:04 | 只看该作者 Only view this author
wzkchem5 发表于 2021-11-7 20:44
电子在0K下的能量为0,因为电子能本身就是以所有粒子相距无穷远且静止为能量零点,而0K下的电子恰好满足 ...

老师您好,请问如果想要计算298K的EA,是不是EA=E(N)+E(e-)-E(N+1)
E(e-)有298K下的实验数据

1万

帖子

0

威望

9857

eV
积分
22093

Level 6 (一方通行)

40#
发表于 Post on 2023-3-6 16:19:51 | 只看该作者 Only view this author
一条君 发表于 2023-3-6 03:56
(1)十分感谢详细解答,想再确认下您的意思
【[A+] +  —— [A] + 以及后续的所有接续的反应总体来说算 ...

(1)不是。
举例:假设你要算两步反应,[A+] + [B-] —— [A] + [B], [B] + [C] —— [BC]。则B-、B、BC里的B部分里带显著负电的原子需要加弥散。A+、A、C,以及BC里的C部分一般不需要
(2)对

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
一条君 + 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?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
本团队长期招收研究生,有意者可私信联系

191

帖子

0

威望

1098

eV
积分
1289

Level 4 (黑子)

39#
发表于 Post on 2023-3-6 10:56:24 | 只看该作者 Only view this author
本帖最后由 一条君 于 2023-3-6 11:35 编辑
wzkchem5 发表于 2023-3-5 17:03
(1)在考虑到我说的“中性分子、阳离子未必不需要加弥散”,以及阴离子也有个别的不需要加弥散,这两种 ...

(1)十分感谢详细解答,想再确认下您的意思
【[A+] + [B-] —— [A] + [B]以及后续的所有接续的反应总体来说算一条反应路径】,同一条反应路径里计算精度必须保持精度一致,即理论级别需统一。可否理解为您是建议[A+] 、[B-] 、 [A]、 [B]结构优化和单点都加弥散呢?
(2)如果要比较中性有机分子、阳离子+中性分子、阴离子的LUMO,应该两者同一级别比较,若一个加了弥散,两者结构优化和单点都加弥撒对吗
想把之前自己的误区纠正一下,再次感谢老师
科音学员

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

GMT+8, 2026-2-20 10:06 , Processed in 0.475721 second(s), 33 queries , Gzip On.

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