| 查看: 1577 | 回复: 7 | ||
[求助]
如何考虑自选污染造成的能量误差
|
| 我最近投了一篇计算反应机理的文章，其中涉及一个单重态的双自由基分子，用的UB3LYP方法，当时没考虑自旋污染，审稿意见回来说让提供S**2的值，发现高斯输出文件里面这个S**2=0.97，自旋淹没以后S**2=0.07.也就是说UB3LYP方法算这个单重态双自由基的时候会出现很大的自旋污染，查阅文献，有很多课题组开发了很多方法来消除自旋污染的影响，但是我们又没有这个程序，不知道大家有没有遇到过这个问题，都是怎么解决的，有没有什么简单的方法可以矫正。 |
回复此楼» 猜你喜欢
国内青年破解费曼难题，从分子层面阐释液体蒸发（气液相变）微观过程
已经有5人回复
-
研究发现一个新的临界系数，独立于临界压缩因子
已经有7人回复
-
物理化学论文润色/翻译怎么收费?
已经有185人回复
-
200年来人们一直无法破解的难题： 气液临界点状态方程
已经有7人回复
-
细胞培养，这22个细节一定要注意！（一）
已经有0人回复
-
关于引力延迟下双星系统角动量不守恒的疑惑---求高手解答！
已经有10人回复
-
CSC访学博后项目获批，外方学校暂停合作该怎么办？
已经有51人回复
-
七嗪类物质合成求助
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 液相色谱重复性差怎么解决
已经有23人回复
- 今年基础不发小册子了，怎么考啊
已经有23人回复
- 2011年注册化工工程师基础考试大纲
已经有25人回复
- 【求助】注册化工工程师的报考条件，报考时间以及考试科目都是什么？
已经有26人回复
- 【求助】XPS如何用碳校正
已经有21人回复
- 【求助】求2010年注册化工工程师基础考试的考试大纲
已经有4人回复
- 【专题】动态光散射基本原理及其在纳米科技中的应用——Zeta电位测量
已经有193人回复
- 【专家一席谈】电镜随谈之三：能谱（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