24小时热门版块排行榜    

查看: 1383  |  回复: 5

JAZZME

新虫 (初入文坛)

[求助] 怎么才能确保Mulliken spin densities是对的? 已有2人参与

如题,请问大家,在算单点能或者做优化时,怎么才能确保Mulliken spin densities是对的?急求,谢谢大家。

发自小木虫Android客户端
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Yan_Jordan

木虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
个人认为这个问题的说法是有问题的

如果说是整个计算的结果的正确性,那实际是你整个计算的问题,这关系到你分子(基团)的结构,基组,计算方法,初始条件(整体电荷,和自旋值等)定义。如果这些条件都没问题,那么得到的每个原子的Mulliken spin 就是你最后的计算结果,Mulliken 只是其中的一种(还算简单的)分析方法,也可以用其他的如Bader等。采取何种分析方法,取决于你要达到什么目的、

所以与其说是Mulliken spin densities 的正确性问题,更应该说是你本身计算方法、参数的问题
2楼2017-04-20 19:26:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

JAZZME

新虫 (初入文坛)

非常感谢您,是这样的,我现在算一个产物的单点能,可是它的能量总是高于过渡态,我检查了一下,就是自旋密度有问题,所以想求教怎么才能使这个结果正确?

发自小木虫Android客户端
3楼2017-04-20 21:13:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

JAZZME

新虫 (初入文坛)

之前通过改变构型,自旋密度对了,能量降下来了,可是换成现在的LACV3P++**基组,自旋密度又不对了,改变构型也不行,所以还有什么好的方法吗?

发自小木虫Android客户端
4楼2017-04-20 21:16:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Yan_Jordan

木虫 (小有名气)

【答案】应助回帖

引用回帖:
4楼: Originally posted by JAZZME at 2017-04-20 21:16:54
之前通过改变构型,自旋密度对了,能量降下来了,可是换成现在的LACV3P++**基组,自旋密度又不对了,改变构型也不行,所以还有什么好的方法吗?

在做开壳计算的时候,其实有时的确会出现很多问题!比如DFT本身在计算自旋就会有本身的问题,如果是分子计算一般会考虑用CASSCF,晶体计算会采用DFT+U。
还有就是本事DFT不同的functional得出的能量也是不一样的。

对于你的问题,你说的自旋密度是指Mulliken spin densities的值,还是用GViewer 做的(cube)然后用做等势面看的。在我看来后者要比Mulliken spin densities更具有说服力。

如果你是用DFT,我觉得你可以尝试用不同的functionnal算一下(你改变构型有做结构优化再计算能量吗?!)。或者CASSCF(这个其实有点麻烦,需要很多尝试,轨道分析等等)。
5楼2017-04-21 02:57:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xytk04

铁杆木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
你是指你的electronic configuration不是你想要的么?那样的话你可以调整orbital或者用fragment=guess来构建一个你想要的electronic configuration的initial guess然后再算,另外你的TS和这个Int的electronic configuration接近么?可不可以读TS的chk?
6楼2017-04-21 10:00:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 JAZZME 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 Z7qfCsdIiD1T 2026-11-10 4/200 2026-11-12 02:23 by 9fOTqj074tIg
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 Z7qfCsdIiD1T 2026-11-10 6/300 2026-11-12 02:23 by 9fOTqj074tIg
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 czAiolzR0IJC 2026-11-10 4/200 2026-11-12 02:12 by 9fOTqj074tIg
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 czAiolzR0IJC 2026-11-10 5/250 2026-11-12 02:12 by 9fOTqj074tIg
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 czAiolzR0IJC 2026-11-10 7/350 2026-11-12 02:03 by 9fOTqj074tIg
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 czAiolzR0IJC 2026-11-10 6/300 2026-11-12 02:01 by 9fOTqj074tIg
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 G2rhSvqvhNzl 2026-11-10 3/150 2026-11-12 01:50 by 9fOTqj074tIg
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 G2rhSvqvhNzl 2026-11-10 5/250 2026-11-12 01:50 by 9fOTqj074tIg
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 C6CRn66RYYXS 2026-11-10 4/200 2026-11-12 01:30 by 9fOTqj074tIg
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 C6CRn66RYYXS 2026-11-10 4/200 2026-11-12 01:28 by 9fOTqj074tIg
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 FChKX8meu1ds 2026-11-09 7/350 2026-11-12 01:06 by 9fOTqj074tIg
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 FChKX8meu1ds 2026-11-09 8/400 2026-11-12 00:55 by 9fOTqj074tIg
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 FChKX8meu1ds 2026-11-09 10/500 2026-11-12 00:46 by 9fOTqj074tIg
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 cqv1zKnzd7hr 2026-11-09 7/350 2026-11-12 00:44 by 9fOTqj074tIg
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 FChKX8meu1ds 2026-11-09 7/350 2026-11-11 17:31 by TgvkbZZl9tqB
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 FChKX8meu1ds 2026-11-09 6/300 2026-11-11 17:03 by TgvkbZZl9tqB
[考博] 27博士 +3 锦鲤送福至 2026-11-10 3/150 2026-11-11 15:53 by 可淡不可忘
[文学芳草园] 厉娜 +3 张丹枫1991 2026-11-09 5/250 2026-11-09 16:53 by 张丹枫1991
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 5sJFDKp2kwuI 2026-11-06 16/800 2026-11-09 11:47 by GwrQVVR6Oenx
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
信息提示
请填处理意见