24小时热门版块排行榜    

查看: 1385  |  回复: 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一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 XJXMuRjD7cK7 2026-11-13 3/150 2026-11-15 05:44 by xSrQNSSjcgoX
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 uKOSCocCb6Fm 2026-11-12 7/350 2026-11-15 02:48 by xSrQNSSjcgoX
[硕博家园] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 G2rhSvqvhNzl 2026-11-10 6/300 2026-11-15 00:27 by Equinoxhua
[考博] 27博士 +6 锦鲤送福至 2026-11-10 6/300 2026-11-14 21:23 by xupeii
[考博] 27博士申请 +5 白白白ly 2026-11-09 5/250 2026-11-13 05:01 by wlf1982
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 TgvkbZZl9tqB 2026-11-11 4/200 2026-11-12 12:28 by ECLTXnJnRSXf
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 Z7qfCsdIiD1T 2026-11-10 4/200 2026-11-12 11:22 by HajghKOMynmz
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 Z7qfCsdIiD1T 2026-11-10 8/400 2026-11-12 11:20 by HajghKOMynmz
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 czAiolzR0IJC 2026-11-10 8/400 2026-11-12 11:11 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 czAiolzR0IJC 2026-11-10 6/300 2026-11-12 11:03 by HajghKOMynmz
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 czAiolzR0IJC 2026-11-10 8/400 2026-11-12 10:49 by HajghKOMynmz
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 C6CRn66RYYXS 2026-11-10 8/400 2026-11-12 10:27 by HajghKOMynmz
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 YHiw1gou1Eda 2026-11-10 6/300 2026-11-12 10:08 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +7 FChKX8meu1ds 2026-11-09 12/600 2026-11-12 09:43 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 TgvkbZZl9tqB 2026-11-11 4/200 2026-11-12 08:15 by HajghKOMynmz
[基金申请] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 TgvkbZZl9tqB 2026-11-11 5/250 2026-11-12 08:15 by HajghKOMynmz
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 C6CRn66RYYXS 2026-11-10 4/200 2026-11-12 06:03 by 9fOTqj074tIg
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 FChKX8meu1ds 2026-11-09 9/450 2026-11-12 05:19 by 9fOTqj074tIg
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 FChKX8meu1ds 2026-11-09 7/350 2026-11-12 05:19 by 9fOTqj074tIg
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 zSbvVX9JFl9J 2026-11-09 8/400 2026-11-12 04:57 by 9fOTqj074tIg
信息提示
请填处理意见