24小时热门版块排行榜    

查看: 7343  |  回复: 14

xiaowu787

木虫 (正式写手)


[交流] 【求助】VASP怎么考虑自旋极化?

想拟合一下Mo的各相的状态方程和单原子的能量,Mo是过渡元素,考虑自旋极化,是不是加上ISPIN=2就行了?其他还有特别需要注意的地方吗?请高手指点一下,谢谢
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★
zzy870720z(金币+2):谢谢指点 2010-11-19 07:59:26
xiaowu787(金币+2):谢谢 2010-11-19 14:01:28
xiaowu787(金币+3):谢谢 2010-11-19 18:43:32
首先,用ISPIN=2是必须的,
其次,可以先估算一下每个Mo原子的单电子数,然后指定一个magmom.
最后,如果你知道体系中有未成对电子,可以设置UNPDOWN(spinup--spindown)
2楼2010-11-19 03:32:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowu787

木虫 (正式写手)


引用回帖:
Originally posted by fegg7502 at 2010-11-19 03:32:40:
首先,用ISPIN=2是必须的,
其次,可以先估算一下每个Mo原子的单电子数,然后指定一个magmom.
最后,如果你知道体系中有未成对电子,可以设置UNPDOWN(spinup--spindown)

请问怎么估算一下每个Mo原子的单电子数,和最外层电子数有什么样的关系?不指定magmom,使用程序的默认值可以吗?谢谢您的回复.
3楼2010-11-19 18:47:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzy870720z(金币+1):谢谢指教 2010-11-19 21:37:01
引用回帖:
Originally posted by xiaowu787 at 2010-11-19 18:47:50:


请问怎么估算一下每个Mo原子的单电子数,和最外层电子数有什么样的关系?不指定magmom,使用程序的默认值可以吗?谢谢您的回复.

看Mo是几价的,,最外层有几个未成对电子。然后你设置MAGMOM就OK了
4楼2010-11-19 20:47:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘小伟110

木虫 (小有名气)



小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
2楼: Originally posted by fegg7502 at 2010-11-19 03:32:40
首先,用ISPIN=2是必须的,
其次,可以先估算一下每个Mo原子的单电子数,然后指定一个magmom.
最后,如果你知道体系中有未成对电子,可以设置UNPDOWN(spinup--spindown)

请问版主,加入自旋或者NUPDOWN的结果,都是让总能变得更低是不?
如果我加入NUPDOWN后,total energy升高啦,是不是说明我不该加呢
5楼2014-03-14 10:26:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
5楼: Originally posted by 刘小伟110 at 2014-03-14 10:26:48
请问版主,加入自旋或者NUPDOWN的结果,都是让总能变得更低是不?
如果我加入NUPDOWN后,total energy升高啦,是不是说明我不该加呢...

一切以total energy变的更负为标准,,NUPDOWN后total energy升高啦,说明体系变的不稳定了,不是说不该加,是说可能加错了,
6楼2014-03-14 10:29:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘小伟110

木虫 (小有名气)



小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
6楼: Originally posted by fegg7502 at 2014-03-14 10:29:50
一切以total energy变的更负为标准,,NUPDOWN后total energy升高啦,说明体系变的不稳定了,不是说不该加,是说可能加错了,...

谢谢版主的回复,我在用vasp计算Al2O3表面H2O吸附时,OH和H共吸附时共吸附能约为8eV,而OH单吸附为2eV,H单吸附能为1eV,显然1+2和8相差太大,得到了大家的质疑。。。我重新摆放OH和H的相对位置,排斥能也不大。。。我通过加入自旋和NUPDOWN=3来重新矫正,果然把能量给算上来了,但是把原先放热的,算成了吸热。。。我这样做对不对啊?是不是以后加入至加入ISPIN=2就好呢,MAGMOM要不要设置呢。。。把所有的计算都设置成NUPDOWN=1,行不。。。
7楼2014-03-14 11:46:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
7楼: Originally posted by 刘小伟110 at 2014-03-14 11:46:02
谢谢版主的回复,我在用vasp计算Al2O3表面H2O吸附时,OH和H共吸附时共吸附能约为8eV,而OH单吸附为2eV,H单吸附能为1eV,显然1+2和8相差太大,得到了大家的质疑。。。我重新摆放OH和H的相对位置,排斥能也不大。。。 ...

这个体系不用加这个NUPDOWN=1,,不用自旋也可以,,能量差别这么大合理的,
单纯羟基有个未配对电子、质子有个未配对电子,,两个配对,放出能量大

» 本帖已获得的红花(最新10朵)

8楼2014-03-14 12:09:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘小伟110

木虫 (小有名气)



小木虫: 金币+0.5, 给个红包,谢谢回帖
送红花一朵
引用回帖:
8楼: Originally posted by fegg7502 at 2014-03-14 12:09:12
这个体系不用加这个NUPDOWN=1,,不用自旋也可以,,能量差别这么大合理的,
单纯羟基有个未配对电子、质子有个未配对电子,,两个配对,放出能量大...

我看的文献中都加的自旋,我也没搞懂,心想没有磁性干嘛要加啊,而且计算的话也会收敛的很慢,所以没办法也要加。。。看vasp介绍,NUPDOWN的设置是固定磁矩,那和MAGMOM又有什么不同呢,对于非铁磁,是不是只设置ISPIN=2,MAGMOM=1就足够了呢。。
9楼2014-03-14 17:35:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘小伟110

木虫 (小有名气)


送红花一朵
引用回帖:
8楼: Originally posted by fegg7502 at 2014-03-14 12:09:12
这个体系不用加这个NUPDOWN=1,,不用自旋也可以,,能量差别这么大合理的,
单纯羟基有个未配对电子、质子有个未配对电子,,两个配对,放出能量大...

我看的文献中都加的自旋,我也没搞懂,心想没有磁性干嘛要加啊,而且计算的话也会收敛的很慢,所以没办法也要加。。。看vasp介绍,NUPDOWN的设置是固定磁矩,那和MAGMOM又有什么不同呢,对于非铁磁,是不是只设置ISPIN=2,MAGMOM=1就足够了呢。。
10楼2014-03-14 17:35:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
引用回帖:
10楼: Originally posted by 刘小伟110 at 2014-03-14 17:35:42
我看的文献中都加的自旋,我也没搞懂,心想没有磁性干嘛要加啊,而且计算的话也会收敛的很慢,所以没办法也要加。。。看vasp介绍,NUPDOWN的设置是固定磁矩,那和MAGMOM又有什么不同呢,对于非铁磁,是不是只设置I ...

只设置ISPIN=2
11楼2014-03-17 08:08:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘小伟110

木虫 (小有名气)



小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
11楼: Originally posted by fegg7502 at 2014-03-17 08:08:22
只设置ISPIN=2...

谢谢版主,再无耻的问您一个问题,对于我所计算的Al2O3表面H2O的吸附,是设置ISMEAR = 2,还是0,-5呢 ,这个有点搞不清,这影响能量吗, SIGMA是不是设置为0.05呢。。。求指教
12楼2014-03-18 15:32:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
12楼: Originally posted by 刘小伟110 at 2014-03-18 15:32:38
谢谢版主,再无耻的问您一个问题,对于我所计算的Al2O3表面H2O的吸附,是设置ISMEAR = 2,还是0,-5呢 ,这个有点搞不清,这影响能量吗, SIGMA是不是设置为0.05呢。。。求指教...

SIGMA设置为0.05,OK,,,
ISMEAR用fermi的那个就OK
13楼2014-03-19 08:17:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘小伟110

木虫 (小有名气)



小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
13楼: Originally posted by fegg7502 at 2014-03-19 08:17:10
SIGMA设置为0.05,OK,,,
ISMEAR用fermi的那个就OK...

那您的意思是优化用ISMEAR=-1了?
14楼2014-03-20 18:10:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

taoyanyaya

新虫 (小有名气)



小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
4楼: Originally posted by fegg7502 at 2010-11-19 20:47:11
看Mo是几价的,,最外层有几个未成对电子。然后你设置MAGMOM就OK了...

请问,不设置NUPDOWN,体系具有磁矩,设置了NUPDOWN=0后,体系磁矩是0,请问这样说明这个体系具不具有磁性呢?谢谢啦
15楼2015-01-20 17:32:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiaowu787 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见