24小时热门版块排行榜    

查看: 2240  |  回复: 10

danyinggian

新虫 (初入文坛)

[求助] 自旋污染

最近计算遇到这样子的:

体系自旋多重度为2
在计算的过程中总是出现这样的提示:
S**2 before annihilation A, after B

从而导致了计算终止。
用了限制性开壳层计算
还是同样的结果
不知道哪位大侠能帮忙解答一下
该如何做?

[ Last edited by gmy1990 on 2011-5-14 at 20:48 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wushidi

木虫 (小有名气)

【答案】应助回帖


小红豆(金币+1): 感谢交流 2011-05-16 12:28:06
danyinggian(金币+3): 2011-05-23 10:51:02
多重度是2的话,说明有单电子,那应该用非限制性的吧
2楼2011-05-16 11:05:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lihb734

铁杆木虫 (职业作家)

站在计算化学入门的门槛上

【答案】应助回帖

★ ★
小红豆(金币+2): 感谢交流 2011-05-16 12:28:22
S**2 before annihilation A, after B,
其实很多时候这个不是真正的出错问题。
楼主应该仔细分析一下这个问题的了
前途光明,出路难觅!
3楼2011-05-16 11:39:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

danyinggian

新虫 (初入文坛)

引用回帖:
Originally posted by wushidi at 2011-05-16 11:05:27:
多重度是2的话,说明有单电子,那应该用非限制性的吧

我一开始使用的的确是非限制性的。。
解决不了
出现了问题
然后我去看了郭勇的那本教程
里面有人建议用限制性的,
结果还是不行。。
4楼2011-05-16 21:48:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangmt

至尊木虫 (著名写手)

我叫MT

语焉不详,掩盖一切可以提供有用的信息,任谁也无法回答你的提问。
一群自以为正义凛然的年轻人将一切不能以科学解释的事情定性为封建迷信并大刀阔斧地进行消灭,其实这是修养不足学识浅薄的一种体现,也是可恶的偏执和愚蠢的自以
5楼2011-05-16 21:54:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mchen10

金虫 (正式写手)

S**2 before annihilation A, after B
是正常提示,S2经过projection得到S2A

出错应该是别的原因,请给出文件最后
“Error terminated ... l???.exe" 信息
否则大家只能玩猜的
更痛苦的是梦没醒路已经不能走了
6楼2011-05-17 03:36:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kaegi

金虫 (小有名气)

【答案】应助回帖

根据我的经验剩下的只能通过改变分子的构像(可能用处不大)和换DFT方法了。
7楼2011-05-18 21:05:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

danyinggian

新虫 (初入文坛)

引用回帖:
Originally posted by mchen10 at 2011-05-17 03:36:43:
S**2 before annihilation A, after B
是正常提示,S2经过projection得到S2A

出错应该是别的原因,请给出文件最后
“Error terminated ... l???.exe" 信息
否则大家只能玩猜的

Annihilation of the first spin contaminant:
S**2 before annihilation     0.7880,   after     0.7506
Convergence failure -- run terminated.
Error termination via Lnk1e in C:\G03W\l502.exe at Fri May 06 00:16:31 2011.
Job cpu time:  0 days 13 hours 49 minutes 43.0 seconds.
File lengths (MBytes):  RWF=    208 Int=      0 D2E=      0 Chk=      3 Scr=      1


这个是结尾。。。
麻烦帮我看下?谢谢~!
8楼2011-05-19 09:19:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yzwyxnb

银虫 (小有名气)

7楼说得对,换构象或者换DFT方法是可能可以解决问题的。
好好学习
9楼2011-05-19 10:16:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangmt

至尊木虫 (著名写手)

我叫MT

【答案】应助回帖

danyinggian(金币+15): 非常感谢,谢谢您的建议,我试下,不行的话得麻烦您再帮我看看 2011-05-19 19:16:51
引用回帖:
Originally posted by danyinggian at 2011-05-19 09:19:19:
Annihilation of the first spin contaminant:
S**2 before annihilation     0.7880,   after     0.7506
Convergence failure -- run terminated.
Error termination via Lnk1e in C:\G03W\l502.exe a ...

S**2=0.75是标准的二重态,所以,你这个计算跟自旋污染可以说没有任何关系。
Convergence failure -- run terminated以及Error termination via Lnk1e in C:\G03W\l502.exe 都说明了是优化过程没有收敛而已。
由于没有更多的信息,不好说应该怎么改进,看在你执着地问,我们就来玩玩猜猜看的游戏:
1、详细监控opt收敛的过程,如果能量有收敛的趋势,那么增加cycle的圈数,继续优化,可能会收敛;
2、如果收敛过程出现振荡,那么增加多少圈都是白费,要么改方法——比如b3lyp改bpw91或者干脆DFT改HF之类的,要么改基组——减小或加大基组试试看,如果有过渡金属就换个赝势试试看,阴离子的加个弥散等等,要么用QC的收敛方式;
3、其实最早该做的,是分析一下你的初始构象是不是合理
一群自以为正义凛然的年轻人将一切不能以科学解释的事情定性为封建迷信并大刀阔斧地进行消灭,其实这是修养不足学识浅薄的一种体现,也是可恶的偏执和愚蠢的自以
10楼2011-05-19 10:33:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 danyinggian 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 5/250 2026-08-08 21:54 by 13586093586
[基金申请] 好奇怪的filecode +3 布布和一二 2026-08-08 4/200 2026-08-08 17:46 by zhanghaozhu
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:22 by oEVWOejN9taj
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 17:02 by oEVWOejN9taj
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 16:47 by oEVWOejN9taj
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 16:47 by oEVWOejN9taj
[考博] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:27 by oEVWOejN9taj
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:19 by oEVWOejN9taj
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 15:59 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 WQyTGMbfH7Cx 2026-08-07 4/200 2026-08-08 15:02 by oEVWOejN9taj
[基金申请] 这样的filecode谁见过 +5 布布和一二 2026-08-08 8/400 2026-08-08 14:54 by archvillain
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:47 by oEVWOejN9taj
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:42 by oEVWOejN9taj
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:39 by oEVWOejN9taj
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
信息提示
请填处理意见