计算化学公社

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

[综合交流] ISC速率常数套文献公式与数据算出来与文献不一致?

[复制链接 Copy URL]

273

帖子

0

威望

4430

eV
积分
4703

Level 6 (一方通行)

本帖最后由 函数与激情 于 2021-1-22 22:26 编辑

大家好,在2019年的一篇JCC文章中基于经典Marcus理论计算了一些分子的ISC与RISC速率常数,公式中所用的参数都在文章中给出(见附件),我在重复他们的数据时算出来的ISC速率常数差的挺多(两个数量级),以下是我的计算过程,如有错误之处,希望这方面比较熟悉的老师不吝赐教,谢谢。

h_bar = 6.582119514e-16 eV·s (此处将普朗克常数替换为了4.1356676969×10-15 eV·s, 数据来源:https://baike.baidu.com/item/%E6 ... %B8%E6%95%B0/812256)
V = 3.471557E-5 eV (Tab 3. 分子1的旋轨耦合矩阵元是0.28 cm-1,这里替换为eV单位)
T = 300 K (计算细节提到)
KB*T = 0.025852 eV (数据来源:https://baike.baidu.com/item/%E7 ... C%E5%B8%B8%E6%95%B0)
delta_E = 0.132 eV (见Tab 3. 第一行)
lambda_T = 0.1244 (见Tab 3. 第一行)

#!/usr/bin/python

import math

pi = 3.1415926535
V = 3.471557E-5
delta_E = 0.132
lambda_T = 0.1244
Boltzmann_constant2 = 0.025852

k1 = (2*pi*(V**2))/h_bar
k2 = (1/(4*pi*Boltzmann_constant2*lambda_T))**0.5
k3 = math.exp(-((delta_E+lambda_T)**2)/(4*Boltzmann_constant2*lambda_T))
k = k1*k2*k3
print(k)



sp210122_222417.png (22.3 KB, 下载次数 Times of downloads: 57)

sp210122_222417.png

sp210122_222405.png (70.93 KB, 下载次数 Times of downloads: 46)

sp210122_222405.png

jcc.25808.pdf

1.1 MB, 下载次数 Times of downloads: 95

2425

帖子

1

威望

6197

eV
积分
8642

Level 6 (一方通行)

2#
发表于 Post on 2021-1-23 13:09:47 | 只看该作者 Only view this author
lambda_T or lambda_S
check out the difference.
High-Performance Computing for You
为您专属定制的高性能计算解决方案

更多讯息,请访问:
https://labitc.top
http://tophpc.top:8080
电邮: ask@hpc4you.top

2425

帖子

1

威望

6197

eV
积分
8642

Level 6 (一方通行)

3#
发表于 Post on 2021-1-23 13:11:21 | 只看该作者 Only view this author
By the way,
import math

then you can use math.pi,

you can import hbar either,
maybe from other module.
High-Performance Computing for You
为您专属定制的高性能计算解决方案

更多讯息,请访问:
https://labitc.top
http://tophpc.top:8080
电邮: ask@hpc4you.top

273

帖子

0

威望

4430

eV
积分
4703

Level 6 (一方通行)

4#
 楼主 Author| 发表于 Post on 2021-1-24 10:38:16 | 只看该作者 Only view this author
abin 发表于 2021-1-23 13:11
By the way,
import math

感谢您的回复,问题已经得到了解决。

正向ISC与反向ISC过程的delta_E绝对值是相同的,但是正负号不同,因此产生了错误。

2425

帖子

1

威望

6197

eV
积分
8642

Level 6 (一方通行)

5#
发表于 Post on 2021-1-24 14:38:53 | 只看该作者 Only view this author
函数与激情 发表于 2021-1-24 10:38
感谢您的回复,问题已经得到了解决。

正向ISC与反向ISC过程的delta_E绝对值是相同的,但是正负号不同 ...

物理化学中肯定有明确讲过,
delta_G 之类的怎么计算。
High-Performance Computing for You
为您专属定制的高性能计算解决方案

更多讯息,请访问:
https://labitc.top
http://tophpc.top:8080
电邮: ask@hpc4you.top

7

帖子

0

威望

45

eV
积分
52

Level 2 能力者

6#
发表于 Post on 2021-7-28 16:57:12 | 只看该作者 Only view this author
假如T1的势能曲线如虚线所示,lambda_T数值为负,请问如何处理?

202107281649079930..png (51.22 KB, 下载次数 Times of downloads: 55)

势能曲线

势能曲线

202107281653033817..png (50.82 KB, 下载次数 Times of downloads: 61)

势能曲线

势能曲线

2425

帖子

1

威望

6197

eV
积分
8642

Level 6 (一方通行)

7#
发表于 Post on 2021-7-28 17:41:41 | 只看该作者 Only view this author
zhanghaoxu 发表于 2021-7-28 16:57
假如T1的势能曲线如虚线所示,lambda_T数值为负,请问如何处理?

注意适用条件。
High-Performance Computing for You
为您专属定制的高性能计算解决方案

更多讯息,请访问:
https://labitc.top
http://tophpc.top:8080
电邮: ask@hpc4you.top

7

帖子

0

威望

45

eV
积分
52

Level 2 能力者

8#
发表于 Post on 2021-7-28 19:38:38 | 只看该作者 Only view this author
请教。下图取自Acc. Chem. Res. 2000, 33, 139-145。问题1:CP点与LS、HS这三个能级与翻转速率之间的关系是什么?
问题2:从JCC的公式看,翻转速率取决于LS、HS两能级之差和两个势能面的另外两个点,即与两个势能面的弯曲度有关,所以,翻转平衡(两态分布)是不是本应考虑势能曲面上的一系列结构、而有另外的公式?
问题3:有没有文献,具体计算并比较了某决速步的SSR速率和TSR速率,或者,估计了两者的量级关系?


问题4:对于低压下的团簇反应,团簇实际处于热状态,对SSR和TSR的反应速率之比有何影响?

202107281720086909..png (73.97 KB, 下载次数 Times of downloads: 64)

202107281720086909..png

7

帖子

0

威望

45

eV
积分
52

Level 2 能力者

9#
发表于 Post on 2021-7-28 19:39:26 | 只看该作者 Only view this author
abin 发表于 2021-7-28 17:41
注意适用条件。

谢谢!

2425

帖子

1

威望

6197

eV
积分
8642

Level 6 (一方通行)

10#
发表于 Post on 2021-7-28 20:59:59 | 只看该作者 Only view this author

仅供参考

DOI: 10.7503/cjcu20210175
DOI: 10.1039/C8DT04882G
DOI: 10.1021/acs.inorgchem.8b00705

以及这个视频提到部分内容:
https://www.cailiaoren.com/vinfo ... id=193&vid=1753

祝好.
High-Performance Computing for You
为您专属定制的高性能计算解决方案

更多讯息,请访问:
https://labitc.top
http://tophpc.top:8080
电邮: ask@hpc4you.top

本版积分规则 Credits rule

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

GMT+8, 2026-2-27 03:39 , Processed in 0.399552 second(s), 23 queries , Gzip On.

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