| 查看: 1997 | 回复: 7 | ||
[求助]
如何考虑自选污染造成的能量误差
|
| 我最近投了一篇计算反应机理的文章，其中涉及一个单重态的双自由基分子，用的UB3LYP方法，当时没考虑自旋污染，审稿意见回来说让提供S**2的值，发现高斯输出文件里面这个S**2=0.97，自旋淹没以后S**2=0.07.也就是说UB3LYP方法算这个单重态双自由基的时候会出现很大的自旋污染，查阅文献，有很多课题组开发了很多方法来消除自旋污染的影响，但是我们又没有这个程序，不知道大家有没有遇到过这个问题，都是怎么解决的，有没有什么简单的方法可以矫正。 |
回复此楼» 猜你喜欢
分子动力学 Plumed教程求助
已经有0人回复
-
研究生学术新锐奖｜宁波微纳材料研讨会参会通知
已经有2人回复
-
物理化学论文润色/翻译怎么收费?
已经有113人回复
-
申请2027年材料/化学类博士
已经有6人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有37人回复
-
求相关合成条件
已经有12人回复
-
2027春季申博自荐
已经有6人回复
-
青基的材料青椒们，可以把硕士联培生送我这里，要研二男生
已经有67人回复
-
大湾区大学刘天辉课题组招收2027级博士生（与中山大学物理学院联合培养）
已经有65人回复
-
申博
已经有2人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 液相色谱重复性差怎么解决
已经有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
