24小时热门版块排行榜    

查看: 1271  |  回复: 5

zhuzunwei

金虫 (小有名气)

[求助] 铁卟啉 自旋污染 DFT

用高斯,b3lyp方法,混合基组,C和H元素用6-31G,N, O和Fe用6-311+G(d),对铁的卟啉化合物进行优化。发现,当自旋多重度为1时,不会出现自旋污染,正常结束。当自旋多重度为3时,会出现自旋污染,不能正常结束。请问应该怎么处理?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ACSI

新虫 (初入文坛)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
zhuzunwei: 金币+10 2013-07-12 15:02:16
本帖内容被屏蔽

2楼2013-05-23 10:28:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhuzunwei

金虫 (小有名气)

谢谢你。我知道 金属应该用赝势基组。但,我使的这种方法是文献中用的。
3楼2013-05-23 10:49:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ylzhao

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
zhuzunwei: 金币+5 2013-05-24 08:17:34
zhuzunwei: 金币+10 2013-07-12 15:02:24
你的体系中Fe是+2还是+3啊,这个体系是挺难优化的。+2时,多重度为1,3,5, +3时为2,4,6,在关键词行加上iop(5/13=1)吧,这样自旋污染大程序也能运行,找到能量最低点的可能性就大一些。个人认为,有些多重度下,基态结构不存在,所以无论怎么算也不收敛(仅供参考)。
4楼2013-05-23 15:29:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dreamyeye

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
zhuzunwei: 金币+5, 有帮助 2013-05-24 15:16:51
zhuzunwei: 金币+5 2013-07-12 15:02:39
因为如果多重度为一,那么程序自动使用限制性dft方案,这是不可能出现自旋污染的,如果用非限制性,自选污染就可能会出现。要是自旋污染不是很大,应该也过得去吧,不收敛不一定是因为有自旋污染造成的。
5楼2013-05-24 12:59:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dreamyeye

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
zhuzunwei: 金币+5 2013-07-12 15:02:42
引用回帖:
4楼: Originally posted by ylzhao at 2013-05-23 15:29:31
你的体系中Fe是+2还是+3啊,这个体系是挺难优化的。+2时,多重度为1,3,5, +3时为2,4,6,在关键词行加上iop(5/13=1)吧,这样自旋污染大程序也能运行,找到能量最低点的可能性就大一些。个人认为,有些多重度下, ...

基态就是能量最低的多重态,因此,基态结构只可能在一个多重度出现。但是我认为你说的基态结构是指收敛的结构,这个每个多重度一定都有,只不过有可能很离谱而已。
6楼2013-05-24 13:01:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhuzunwei 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 jCd0dEvKHShX 2026-08-29 5/250 2026-08-30 19:48 by hZiFeudZyoGR
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jCd0dEvKHShX 2026-08-29 4/200 2026-08-30 19:35 by hZiFeudZyoGR
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 18:50 by hZiFeudZyoGR
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +7 gy1nBQXYQJqL 2026-08-29 8/400 2026-08-30 18:07 by hZiFeudZyoGR
[基金申请] 面上意见出来了 +9 黄鸟于飞Chao 2026-08-29 18/900 2026-08-30 16:47 by 黄鸟于飞Chao
[基金申请] 能否申诉? +4 echo8914667 2026-08-30 4/200 2026-08-30 16:03 by sundiv
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 jCd0dEvKHShX 2026-08-29 5/250 2026-08-30 14:03 by GkCVyODlx4mR
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 gy1nBQXYQJqL 2026-08-29 8/400 2026-08-30 12:03 by l0VvVHGBGRLv
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 11:51 by l0VvVHGBGRLv
[基金申请] 29号明天会评吗 +3 笨笨唐 2026-08-28 3/150 2026-08-30 09:29 by cww8181
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ASdOkHsho7FD 2026-08-28 11/550 2026-08-30 06:00 by ZPa0EcMwuECS
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 05:48 by ZPa0EcMwuECS
[基金申请] 国社科申报系统有变化 +3 kynobel 2026-08-23 3/150 2026-08-29 21:57 by 余韵清
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见