24小时热门版块排行榜    

查看: 1218  |  回复: 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 的主题更新
信息提示
请填处理意见