计算化学公社

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

[Gaussian/gview] 求助!UV-vis光谱的计算值和实验数据相差较大怎么解决??

[复制链接 Copy URL]

7

帖子

0

威望

211

eV
积分
218

Level 3 能力者

求助!UV-vis光谱的计算值和实验数据相差较大怎么解决??

202204140954006688..png (138.42 KB, 下载次数 Times of downloads: 54)

计算值

计算值

202204140953101686..png (44.23 KB, 下载次数 Times of downloads: 24)

实验值

实验值

RClSMDoUV.gjf

1.72 KB, 下载次数 Times of downloads: 36

输入文件

RClSMDoUV.out

45.02 KB, 下载次数 Times of downloads: 6

输出文件

7

帖子

0

威望

211

eV
积分
218

Level 3 能力者

12#
 楼主 Author| 发表于 Post on 2022-4-14 15:48:02 | 只看该作者 Only view this author
exity 发表于 2022-4-14 10:34
如果你这个结构已经优化过了,就改成这个关键词再看看看

# PBE1PBE def2tzvp td=(nstates=100,root=1) s ...

谢谢您!

368

帖子

0

威望

4677

eV
积分
5045

Level 6 (一方通行)

11#
发表于 Post on 2022-4-14 15:16:47 | 只看该作者 Only view this author
xianz 发表于 2022-4-14 15:05
谢谢您!计算后峰强度与实验值不匹配有什么解决方法吗?

峰强度的话再多算一些态再看看,200附近算的态有点少了

368

帖子

0

威望

4677

eV
积分
5045

Level 6 (一方通行)

10#
发表于 Post on 2022-4-14 15:13:05 | 只看该作者 Only view this author
xianz 发表于 2022-4-14 14:56
谢谢!按照您的建议重新计算后峰强度和实验值不匹配,这样有什么解决方法吗?

这个结果和实验符合的相当不错了,实验上300附近那个峰是有振动分辨的,你要算出那两个峰的话需要算振动分辨的吸收光谱。200之前的没有是你算的态的数目比较少。

记住,你所用的理论方法里面采用了很多近似,而且计算的垂直激发能和吸收光谱的最大峰值不是严格对应的关系,所以实验和理论之间必然会存在差异,只要计算结果和实验结果在一定的误差范围以内,那么理论计算的结果就可以认为是可靠的,可以在此基础上进行分析。

对于吸收谱和荧光谱,如果实验值和理论值的差别在0.2 eV以内,那么就可以认为符合得很好了,在0.3 eV以内也是可以接受的,0.3~0.4 eV稍微就有点大,大于0.4 eV就要认真分析原因,看看是什么原因造成的。如果你看到文献里实验和理论符合的很好,比如只差了几个nm那种,基本上是误差抵消的结果,不是说理论本身就有那么高的精度。

7

帖子

0

威望

211

eV
积分
218

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2022-4-14 15:05:19 | 只看该作者 Only view this author

谢谢您!计算后峰强度与实验值不匹配有什么解决方法吗?

7

帖子

0

威望

211

eV
积分
218

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2022-4-14 14:56:27 | 只看该作者 Only view this author
本帖最后由 xianz 于 2022-4-14 14:59 编辑
pwzhou 发表于 2022-4-14 10:07
你这是没指定泛函和基组,所以是用HF/STO-3G计算的,这是没法用的,实际这样的体系至少先用B3LYP/TZVP计算 ...

谢谢!按照您的建议重新计算后峰强度和实验值不匹配,这样有什么解决方法吗?

202204141459128882..png (106.69 KB, 下载次数 Times of downloads: 20)

202204141459128882..png

393

帖子

1

威望

5093

eV
积分
5506

Level 6 (一方通行)

7#
发表于 Post on 2022-4-14 11:44:25 | 只看该作者 Only view this author
QuantumicGuy 发表于 2022-4-14 11:18
PBE0很好用,一般20个态就够了,他的体系也很小。

我算个甲基橙50个态和100个态的深紫外吸收区都有显著不同,它这个分子设100个态看看比较保险,如果算力够的话,20个态目测是不够的。

120

帖子

0

威望

2361

eV
积分
2481

Level 5 (御坂)

#未来可寄

6#
发表于 Post on 2022-4-14 11:20:51 | 只看该作者 Only view this author
参考

202204141120245923..png (22.7 KB, 下载次数 Times of downloads: 23)

202204141120245923..png

120

帖子

0

威望

2361

eV
积分
2481

Level 5 (御坂)

#未来可寄

5#
发表于 Post on 2022-4-14 11:18:50 | 只看该作者 Only view this author
exity 发表于 2022-4-14 10:34
如果你这个结构已经优化过了,就改成这个关键词再看看看

# PBE1PBE def2tzvp td=(nstates=100,root=1) s ...

PBE0很好用,一般20个态就够了,他的体系也很小。

368

帖子

0

威望

4677

eV
积分
5045

Level 6 (一方通行)

4#
发表于 Post on 2022-4-14 10:44:25 | 只看该作者 Only view this author
exity 发表于 2022-4-14 10:34
如果你这个结构已经优化过了,就改成这个关键词再看看看

# PBE1PBE def2tzvp td=(nstates=100,root=1) s ...

计算吸收一般而言并没有什么特殊的理由不用Gaussian默认的IEFPCM而使用SMD

393

帖子

1

威望

5093

eV
积分
5506

Level 6 (一方通行)

3#
发表于 Post on 2022-4-14 10:34:49 | 只看该作者 Only view this author
如果你这个结构已经优化过了,就改成这个关键词再看看看

# PBE1PBE def2tzvp td=(nstates=100,root=1) scrf=(smd,solvent=methanol)

368

帖子

0

威望

4677

eV
积分
5045

Level 6 (一方通行)

2#
发表于 Post on 2022-4-14 10:07:24 | 只看该作者 Only view this author
你这是没指定泛函和基组,所以是用HF/STO-3G计算的,这是没法用的,实际这样的体系至少先用B3LYP/TZVP计算一下看看。

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

GMT+8, 2026-1-24 20:13 , Processed in 2.025984 second(s), 24 queries , Gzip On.

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