计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: Theor_Comp
打印 Print 上一主题 Last thread 下一主题 Next thread

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

[复制链接 Copy URL]

195

帖子

0

威望

699

eV
积分
894

Level 4 (黑子)

16#
 楼主 Author| 发表于 Post on 2017-12-20 13:37:25 | 只看该作者 Only view this author
sobereva 发表于 2017-12-20 13:34
你补的是氢,一个就够

sob老师,如果D和A之间连接的是个双键,打断这个双键,分别算D的电离势和A的电子亲和能,那此时应该分别在D和A上补几个氢呢?

195

帖子

0

威望

699

eV
积分
894

Level 4 (黑子)

17#
 楼主 Author| 发表于 Post on 2017-12-20 16:08:08 | 只看该作者 Only view this author
sobereva 发表于 2017-12-20 13:34
你补的是氢,一个就够

Sob老师,还有就是D-A分子优化后,D和A分别已有一定构型。那计算D的电离势和A的电子亲和能时,还需要分别重新单独优化D部分和单独优化A部分吗,还是说只要取D-A分子中D的构型(直接加上氢)和A的构型(直接加一个氢),直接算单点就行呢?

6万

帖子

99

威望

5万

eV
积分
120061

管理员

公社社长

18#
发表于 Post on 2017-12-20 16:14:15 | 只看该作者 Only view this author
Theor_Comp 发表于 2017-12-20 13:37
sob老师,如果D和A之间连接的是个双键,打断这个双键,分别算D的电离势和A的电子亲和能,那此时应该分别 ...

这种情况两个片段间耦合很厉害,其实打断这种做法已经意义不大了。
北京科音自然科学研究中心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

195

帖子

0

威望

699

eV
积分
894

Level 4 (黑子)

19#
 楼主 Author| 发表于 Post on 2017-12-20 16:29:43 | 只看该作者 Only view this author
sobereva 发表于 2017-12-20 16:14
这种情况两个片段间耦合很厉害,其实打断这种做法已经意义不大了。

Sob老师,还有就是D-A分子优化后,D和A分别已有一定构型。那计算D的电离势和A的电子亲和能时,还需要分别重新单独优化D部分和单独优化A部分吗,还是说只要取D-A分子中D的构型(直接加上氢)和A的构型(直接加一个氢),直接算单点就行呢?

6万

帖子

99

威望

5万

eV
积分
120061

管理员

公社社长

20#
发表于 Post on 2017-12-20 16:33:34 | 只看该作者 Only view this author
Theor_Comp 发表于 2017-12-20 16:29
Sob老师,还有就是D-A分子优化后,D和A分别已有一定构型。那计算D的电离势和A的电子亲和能时,还需要分别 ...

建议把氢的位置优化一下(毕竟gview加的氢的位置没有物理意义),其它部分保持冻结
北京科音自然科学研究中心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

195

帖子

0

威望

699

eV
积分
894

Level 4 (黑子)

21#
 楼主 Author| 发表于 Post on 2017-12-20 16:45:34 | 只看该作者 Only view this author
sobereva 发表于 2017-12-20 16:33
建议把氢的位置优化一下(毕竟gview加的氢的位置没有物理意义),其它部分保持冻结

其他部分如何冻结呢?
是把其他原子(除加的这个氢外)序号添加到坐标下面吗?opt=modredundant

X 1 F
X 2 F
X 3 F
等等

6万

帖子

99

威望

5万

eV
积分
120061

管理员

公社社长

22#
发表于 Post on 2017-12-20 18:11:17 | 只看该作者 Only view this author
Theor_Comp 发表于 2017-12-20 16:45
其他部分如何冻结呢?
是把其他原子(除加的这个氢外)序号添加到坐标下面吗?opt=modredundant

冻结的原子后头写-1,不冻结的写0

# B3LYP/6-31G** opt

Title Card Required

0 1
O      0          0.00000000    0.00000000   -0.11081188
H     -1          0.00000000    0.58397589    0.44324751
H     -1          0.00000000   -0.58397589    0.44324751

评分 Rate

参与人数
Participants 1
eV +4 收起 理由
Reason
cyx98 + 4 正解

查看全部评分 View all ratings

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

195

帖子

0

威望

699

eV
积分
894

Level 4 (黑子)

23#
 楼主 Author| 发表于 Post on 2017-12-20 18:15:52 | 只看该作者 Only view this author
sobereva 发表于 2017-12-20 18:11
冻结的原子后头写-1,不冻结的写0

# B3LYP/6-31G** opt

非常感谢您细致耐心的解答!

3

帖子

0

威望

27

eV
积分
30

Level 2 能力者

24#
发表于 Post on 2018-7-24 08:59:47 | 只看该作者 Only view this author
记得sob老师在博客里提到算电子亲和势要加弥散函数
那么在这里是只在算E(N+1)时在基组中加入弥散函数吗?

26

帖子

0

威望

697

eV
积分
723

Level 4 (黑子)

25#
发表于 Post on 2020-9-15 16:58:43 | 只看该作者 Only view this author
借楼问一下,在计算垂直电离能的时候,先优化基态的结构也就是电中性的结构,再取优化完成的结构改电荷量和自旋多重度,这样的计算过程正确吗?因为出现了l502的问题,参考解决scf不收敛的帖子也没搞定。。。

206

帖子

0

威望

4301

eV
积分
4507

科音成员

26#
发表于 Post on 2020-9-15 17:46:43 | 只看该作者 Only view this author
84015917 发表于 2020-9-15 16:58
借楼问一下,在计算垂直电离能的时候,先优化基态的结构也就是电中性的结构,再取优化完成的结构改电荷量和 ...

原则上没问题,就是根据定义去计算而已

评分 Rate

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

查看全部评分 View all ratings

1万

帖子

0

威望

8959

eV
积分
20702

Level 6 (一方通行)

27#
发表于 Post on 2020-9-15 18:22:16 | 只看该作者 Only view this author
84015917 发表于 2020-9-15 16:58
借楼问一下,在计算垂直电离能的时候,先优化基态的结构也就是电中性的结构,再取优化完成的结构改电荷量和 ...

阳离子很难收敛的话,建议考虑阳离子是否有显著的多参考态性质,如果有的话,需要用CASPT2等方法;如果中性分子基本是单参考态的,那么也可考虑EOM-IP-CCSD

评分 Rate

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

6万

帖子

99

威望

5万

eV
积分
120061

管理员

公社社长

28#
发表于 Post on 2020-9-28 13:52:59 | 只看该作者 Only view this author
gfcjyb 发表于 2018-7-24 08:59
记得sob老师在博客里提到算电子亲和势要加弥散函数
那么在这里是只在算E(N+1)时在基组中加入弥散函数吗?

要加都加,不同基组下能量求差没有意义

评分 Rate

参与人数
Participants 1
eV +4 收起 理由
Reason
一条君 + 4 谢谢

查看全部评分 View all ratings

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

22

帖子

0

威望

1151

eV
积分
1173

Level 4 (黑子)

29#
发表于 Post on 2021-11-7 20:29:22 | 只看该作者 Only view this author
借楼问一个小白问题,EA指中性分子与中性分子多了一个电子的物种之间的能量差,这个值EA=E(N)-E(N+1),这里计算时不用考虑电子的能量吗?

601

帖子

1

威望

4776

eV
积分
5397

Level 6 (一方通行)

30#
发表于 Post on 2021-11-7 20:43:55 | 只看该作者 Only view this author
yejlyejl 发表于 2021-11-7 20:29
借楼问一个小白问题,EA指中性分子与中性分子多了一个电子的物种之间的能量差,这个值EA=E(N)-E(N+1),这里 ...

算的不就是电子的能量嘛
现代化学以狄拉克的一句“一切化学问题业已解决”为嚆矢。滥觞于经验主义传统的期望正失去它们的借鉴意义。但面对看似不可达的通往天堂之阶梯,我想循伍德沃德“最好的模型是你底物的对映异构体”的信仰好过过早地振翮。
我们怀揣热忱的灵魂天然被赋予对第一性的追求,不屑于单一坐标的约束,钟情于势能面彼端的芬芳。但

本版积分规则 Credits rule

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

GMT+8, 2025-8-12 21:01 , Processed in 0.203669 second(s), 27 queries , Gzip On.

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