计算化学公社

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

[新手求助] UKS和TDDFT关于T1结构优化结果不一致

[复制链接 Copy URL]

9

帖子

0

威望

83

eV
积分
92

Level 2 能力者

结构优化,按照教程采用UKS和TDDFT分别进行T1结构优化,对比UKS和TDDFT结果发现优化的结构不一致,具体的表现为苯环结构采用UKS优化后还是平面结构(肉眼观察的),而采用TDDFT优化的结果是苯环结构发生明显的弯曲。
实验条件:光催化纯水,温度室温(25±2)。在这个条件下我的目标化合物是中性状态。分子包含CHONF五种元素。
计算步骤:
UKS关键词:# opt freq cam-b3lyp/6-31+g(d,p) scrf=(iefpcm,solvent=water)

M0004

0 3
TDDFT关键词:
# opt freq td=(triplet,nstates=7,root=1,) cam-b3lyp/6-31+g(d,p) scrf=(iefpcm,solvent=water)

M0004

0 1


计算过程:采用UKS计算时,跑了16次以后就收敛了。采用TDDFT计算时,第一次优化震荡后,取第一次震荡下的最后一次结构进行第二次优化的初始结构,一直重复这个操作,第三次算了106次后收敛。
计算结果:对比结构发现,TDDFT计算的结构中苯环部分发现明显的扭曲,而UKS方法苯环没有扭曲。
问题:
(1)我的计算关键词和选择的基组泛函是否合适。为什么会出现UKS和TD结果差别这么大。
(2)哪个结果更可信,为什么,是从哪些方面可以进行评价的。
(3)为什么TDDFT优化过程中总是频繁的震荡。

烦请各位老师解答指教

202607082034123060..png (622.98 KB, 下载次数 Times of downloads: 1)

两种方法结果对比

两种方法结果对比

202607082016029903..jpg (136.67 KB, 下载次数 Times of downloads: 1)

TDDFT优化过程

TDDFT优化过程

6万

帖子

99

威望

6万

eV
积分
128518

管理员

公社社长

2#
发表于 Post on 2026-7-8 20:41:36 | 只看该作者 Only view this author
T1结构优化一律以UKS为准,原理最简单、稳健,还能充分考虑轨道弛豫效应,速度还快得多得多
用6-31+G**远不如用6-311G**
北京科音自然科学研究中心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

9

帖子

0

威望

83

eV
积分
92

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2026-7-8 20:44:20 | 只看该作者 Only view this author
sobereva 发表于 2026-7-8 20:41
T1结构优化一律以UKS为准,原理最简单、稳健,还能充分考虑轨道弛豫效应,速度还快得多得多
用6-31+G**远 ...

谢谢sobereva老师。那如果我需要算单点能的话,其实S1(TD算法)和T1(UKS算法),算法不一致,这样的还能直接进行比较吗?

6万

帖子

99

威望

6万

eV
积分
128518

管理员

公社社长

4#
发表于 Post on 2026-7-8 21:19:56 | 只看该作者 Only view this author
给回忆点根烟 发表于 2026-7-8 20:44
谢谢sobereva老师。那如果我需要算单点能的话,其实S1(TD算法)和T1(UKS算法),算法不一致,这样的还 ...

UKS和TDDFT算的激发能会存在差异,不能直接对比。相差甚至比如0.2 eV都是正常的。原理上UKS算的更准
北京科音自然科学研究中心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

9

帖子

0

威望

83

eV
积分
92

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2026-7-9 18:34:35 | 只看该作者 Only view this author
sobereva 发表于 2026-7-8 21:19
UKS和TDDFT算的激发能会存在差异,不能直接对比。相差甚至比如0.2 eV都是正常的。原理上UKS算的更准

谢谢sobereva老师,我在看看咱们公社相关的帖子,看看是哪里出问题了。

9

帖子

0

威望

83

eV
积分
92

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 2026-7-9 20:31:17 | 只看该作者 Only view this author
我需要对比T1和S1的能垒来解释一下我的机理路径,请问有没有老师知道我目前的结构或者是算法如何调整能通过TD算出一个正常的T1结构,我重复了两三遍了,现在是真的麻了。

202607092029248494..png (71.27 KB, 下载次数 Times of downloads: 1)

202607092029248494..png

1万

帖子

0

威望

1万

eV
积分
22876

Level 6 (一方通行)

7#
发表于 Post on 2026-7-10 08:54:52 | 只看该作者 Only view this author
给回忆点根烟 发表于 2026-7-9 20:31
我需要对比T1和S1的能垒来解释一下我的机理路径,请问有没有老师知道我目前的结构或者是算法如何调整能通过 ...

什么叫正常?为什么认为你现在算的结构不正常?
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php,墙内下载ORCA:https://www.faccts.de/customer)开发团队成员
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
本团队长期招收研究生,有意者可私信联系

9

帖子

0

威望

83

eV
积分
92

Level 2 能力者

8#
 楼主 Author| 发表于 Post on 2026-7-10 08:57:01 | 只看该作者 Only view this author
wzkchem5 发表于 2026-7-10 08:54
什么叫正常?为什么认为你现在算的结构不正常?

老师,就是目前两种方法算出来结果不一致,而且我算了两遍。我不确定用哪一个才是正确的,不知道如何进行排除

本版积分规则 Credits rule

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

GMT+8, 2026-7-24 15:11 , Processed in 0.227478 second(s), 23 queries , Gzip On.

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