计算化学公社

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

[新手求助] 计算C,Ti原子的单点能不收敛

[复制链接 Copy URL]

84

帖子

0

威望

213

eV
积分
297

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
计算C,Ti原子的单点能不收敛,采用了解决scf 不收敛的方法依旧不收敛,麻烦各位老师帮我看看原因及解决方法

C.out

25.8 KB, 下载次数 Times of downloads: 7

TI.LOG

25.74 KB, 下载次数 Times of downloads: 2

1万

帖子

0

威望

7400

eV
积分
18164

Level 6 (一方通行)

2#
发表于 Post on 2023-11-18 17:23:56 | 只看该作者 Only view this author
自旋多重度设置错误,按Hund规则重新判断两者的自旋多重度,再算
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

84

帖子

0

威望

213

eV
积分
297

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2023-11-18 18:31:36 | 只看该作者 Only view this author
wzkchem5 发表于 2023-11-18 17:23
自旋多重度设置错误,按Hund规则重新判断两者的自旋多重度,再算

谢谢老师。我把C原子的换成3重度之后计算合理,那么,我存在一个新的问题,就是我想要计算TI4Ta4C12团簇的结合能,根据结合能的定义,是团簇的能量减去各原子的能量,在计算团簇时我计算了其单重度,三重度,和五重度,那么,计算结合能时我需要将原子的自旋多重度和团簇的自旋多重度相对应,然而C原子的其他多重度不合理,请问我该怎么处理呢?

1万

帖子

0

威望

7400

eV
积分
18164

Level 6 (一方通行)

4#
发表于 Post on 2023-11-18 19:39:36 | 只看该作者 Only view this author
abu123 发表于 2023-11-18 11:31
谢谢老师。我把C原子的换成3重度之后计算合理,那么,我存在一个新的问题,就是我想要计算TI4Ta4C12团簇 ...

谁告诉你自旋多重度要对应了?好好复习一下涉及共价键断裂时的结合能怎么算
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

84

帖子

0

威望

213

eV
积分
297

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2023-11-18 20:06:21 | 只看该作者 Only view this author
wzkchem5 发表于 2023-11-18 19:39
谁告诉你自旋多重度要对应了?好好复习一下涉及共价键断裂时的结合能怎么算

好的,谢谢老师,我再看看

84

帖子

0

威望

213

eV
积分
297

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2023-11-19 12:48:41 | 只看该作者 Only view this author
abu123 发表于 2023-11-18 20:06
好的,谢谢老师,我再看看

老师,请问涉及共价键断裂的结合能应该怎么算?我查了一些资料,但是很多都是直接给出结合能,没有查到计算团簇的这种情况下的结合能的方法。😱

3811

帖子

4

威望

8000

eV
积分
11891

Level 6 (一方通行)

MOKIT开发者

7#
发表于 Post on 2023-11-19 13:41:06 | 只看该作者 Only view this author
abu123 发表于 2023-11-19 12:48
老师,请问涉及共价键断裂的结合能应该怎么算?我查了一些资料,但是很多都是直接给出结合能,没有查到计 ...

没有“计算结合能时我需要将原子的自旋多重度和团簇的自旋多重度相对应”、“然而C原子的其他多重度不合理”这些说法。如果你要算TI4Ta4C12团簇全分解为一个个原子的话,那就是减去相应倍数的原子的能量之和。不过,PBE泛函是算不准的,只能说当你写对了原子的自旋多重度后,计算步骤没错。
自动做多参考态计算的程序MOKIT

1万

帖子

0

威望

7400

eV
积分
18164

Level 6 (一方通行)

8#
发表于 Post on 2023-11-19 17:26:09 | 只看该作者 Only view this author
abu123 发表于 2023-11-19 05:48
老师,请问涉及共价键断裂的结合能应该怎么算?我查了一些资料,但是很多都是直接给出结合能,没有查到计 ...

对于基态的结合能,断键前后的自旋多重度都应该选择让能量最低的自旋多重度(也有的认为应当选择让自由能最低的自旋多重度,大同小异)
激发态比较复杂,另说
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

84

帖子

0

威望

213

eV
积分
297

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2023-11-19 17:58:33 | 只看该作者 Only view this author
zjxitcc 发表于 2023-11-19 13:41
没有“计算结合能时我需要将原子的自旋多重度和团簇的自旋多重度相对应”、“然而C原子的其他多重度不合 ...

对,这里只是用pbe0泛函做个优化,然后大致进行判断与筛选。后续会做进一步精度的计算。这里就是计算出单个原子的能量,然后用相应倍数的原子求总能量,最后再作差

84

帖子

0

威望

213

eV
积分
297

Level 3 能力者

10#
 楼主 Author| 发表于 Post on 2023-11-19 18:00:08 | 只看该作者 Only view this author
zjxitcc 发表于 2023-11-19 13:41
没有“计算结合能时我需要将原子的自旋多重度和团簇的自旋多重度相对应”、“然而C原子的其他多重度不合 ...

说法方面我理解的还不够深刻,会慢慢学习正确表达,谢谢老师的改正

84

帖子

0

威望

213

eV
积分
297

Level 3 能力者

11#
 楼主 Author| 发表于 Post on 2023-11-19 18:29:17 | 只看该作者 Only view this author
wzkchem5 发表于 2023-11-19 17:26
对于基态的结合能,断键前后的自旋多重度都应该选择让能量最低的自旋多重度(也有的认为应当选择让自由能 ...

也就是说,假设Ti4Ta4C12团簇能量最低的自旋多重态为5重态,C原子为三重态,Ti原子为5重态,Ta原子为6重态,那么,要计算基态结合能时我需要分别用到这几个态的能量,我这么理解对吗?老师?

1万

帖子

0

威望

7400

eV
积分
18164

Level 6 (一方通行)

12#
发表于 Post on 2023-11-19 19:49:56 | 只看该作者 Only view this author
abu123 发表于 2023-11-19 11:29
也就是说,假设Ti4Ta4C12团簇能量最低的自旋多重态为5重态,C原子为三重态,Ti原子为5重态,Ta原子为6重 ...

对,因为基态结合能指的是基态原子形成基态分子的能量变化,所以原子和分子都要取各自的基态自旋多重度
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

84

帖子

0

威望

213

eV
积分
297

Level 3 能力者

13#
 楼主 Author| 发表于 Post on 2023-11-19 21:14:23 | 只看该作者 Only view this author
wzkchem5 发表于 2023-11-19 19:49
对,因为基态结合能指的是基态原子形成基态分子的能量变化,所以原子和分子都要取各自的基态自旋多重度

好的,我理解啦,谢谢老师

本版积分规则 Credits rule

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

GMT+8, 2024-11-26 07:29 , Processed in 0.192997 second(s), 23 queries , Gzip On.

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