| 查看: 1687 | 回复: 7 | ||
[求助]
如何考虑自选污染造成的能量误差
|
| 我最近投了一篇计算反应机理的文章，其中涉及一个单重态的双自由基分子，用的UB3LYP方法，当时没考虑自旋污染，审稿意见回来说让提供S**2的值，发现高斯输出文件里面这个S**2=0.97，自旋淹没以后S**2=0.07.也就是说UB3LYP方法算这个单重态双自由基的时候会出现很大的自旋污染，查阅文献，有很多课题组开发了很多方法来消除自旋污染的影响，但是我们又没有这个程序，不知道大家有没有遇到过这个问题，都是怎么解决的，有没有什么简单的方法可以矫正。 |
回复此楼» 猜你喜欢
298分 070300求调剂
已经有15人回复
-
硕士研究生调剂
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有60人回复
-
湖南大学张世国团队招26年第二批博士—离子液体、高分子与凝胶材料
已经有0人回复
-
湖南农业大学能源材料创新团队招收2025化学专业研究生（学硕）调剂
已经有0人回复
-
2026年西华大学材料学院-材料近净成形与表面技术研究团队-招收调剂研究生
已经有9人回复
-
调剂信息
已经有0人回复
-
河北农大生物与医药专业招收调剂，名额充足
已经有0人回复
-
齐鲁工业大学 国家重点实验室-新能源、锂离子电池、材料、化学方向
已经有0人回复
-
杰青/长江团队招收调剂
已经有0人回复
-
中南林业科技大学罗勇锋教授团队2026年招收材料类调剂研究生
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 液相色谱重复性差怎么解决
已经有23人回复
- 今年基础不发小册子了，怎么考啊
已经有23人回复
- 2011年注册化工工程师基础考试大纲
已经有25人回复
- 【求助】注册化工工程师的报考条件，报考时间以及考试科目都是什么？
已经有26人回复
- 【求助】XPS如何用碳校正
已经有21人回复
- 【求助】求2010年注册化工工程师基础考试的考试大纲
已经有4人回复
- 【专题】动态光散射基本原理及其在纳米科技中的应用——Zeta电位测量
已经有194人回复
- 【专家一席谈】电镜随谈之三：能谱（EDS）的一些问题
已经有232人回复
virtualzx
木虫 (著名写手)
- QC强帖: 2
- 应助: 263 (大学生)
- 金币: 7161.3
- 红花: 54
- 帖子: 1605
- 在线: 317.6小时
- 虫号: 2069080
- 注册: 2012-10-18
- 性别: GG
- 专业: 理论和计算化学
2楼2013-11-08 02:52:04
枪下游魂
木虫 (著名写手)
- 应助: 718 (博后)
- 金币: 8336.9
- 散金: 175
- 红花: 80
- 帖子: 1498
- 在线: 542.6小时
- 虫号: 2620258
- 注册: 2013-08-28
- 专业: 理论和计算化学
3楼2013-11-08 08:21:56
4楼2013-11-08 08:38:29
5楼2013-11-08 11:07:14
sobereva
至尊木虫 (著名写手)
本人已永久离开小木虫
- QC强帖: 17
- 应助: 226 (大学生)
- 贵宾: 5.4
- 金币: 2180.5
- 红花: 217
- 帖子: 1241
- 在线: 366.8小时
- 虫号: 281271
- 注册: 2006-09-23
- 专业: 宇宙学
6楼2013-11-08 16:20:36
7楼2013-11-08 18:52:34
virtualzx
木虫 (著名写手)
- QC强帖: 2
- 应助: 263 (大学生)
- 金币: 7161.3
- 红花: 54
- 帖子: 1605
- 在线: 317.6小时
- 虫号: 2069080
- 注册: 2012-10-18
- 性别: GG
- 专业: 理论和计算化学
【答案】应助回帖
★ ★ ★
gmy1990: 金币+3 2013-11-09 07:19:39
gmy1990: 金币+3 2013-11-09 07:19:39
|
我觉得用自旋算符计算DFT的自旋本征值本来就有点问题。理论上来说，KS波函数只能给出和正确波函数一样的电子密度，而计算自旋本征值需要一阶密度矩阵。可以把KS波函数做HF波函数来进行近似计算，不过结果是什么意义就说不清楚了，只能说近似等于自旋本征值吧。 在HF水平上，用ROHF可以得到正确的自旋态，但是必须要正确指定Roothaan系数a和b才行。而且一般需要手工选择初始轨道，比Unrestricted的计算要复杂一些。不过不像ROHF，因为用的是KS轨道，RODFT不能保证没有自旋污染，只能说是近似吧。我甚至想不出能够严格地计算DFT自旋污染程度的方法。 |
8楼2013-11-09 02:22:46
