| 查看: 1655 | 回复: 7 | ||
[求助]
如何考虑自选污染造成的能量误差
|
| 我最近投了一篇计算反应机理的文章，其中涉及一个单重态的双自由基分子，用的UB3LYP方法，当时没考虑自旋污染，审稿意见回来说让提供S**2的值，发现高斯输出文件里面这个S**2=0.97，自旋淹没以后S**2=0.07.也就是说UB3LYP方法算这个单重态双自由基的时候会出现很大的自旋污染，查阅文献，有很多课题组开发了很多方法来消除自旋污染的影响，但是我们又没有这个程序，不知道大家有没有遇到过这个问题，都是怎么解决的，有没有什么简单的方法可以矫正。 |
回复此楼» 猜你喜欢
七嗪类物质合成求助
已经有0人回复
-
需要一个圆偏振光诱导化学反应，请问使用圆偏振光作为光源，应该需要哪些设备搭建
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有106人回复
-
[电子教材]天津大学物理化学（第七版）上、下册
已经有3人回复
-
阴离子交换膜电解二氧化碳还原
已经有0人回复
-
九江学院2026年最新高层次人才招聘公告
已经有0人回复
-
深圳信息职业技术学院2026年教师招聘公告（最新）
已经有0人回复
-
求ENDNOTE破解软件安装包
已经有4人回复
-
研究生招生
已经有0人回复
-
河海大学 （211.双一流高校）钙钛矿光电实验室丁勇团队欢迎有志青年！
已经有3人回复
-
黑龙江省原子能研究院化学课题组招收2026级化学专业或化学相关专业硕士生
已经有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
