24小时热门版块排行榜    

查看: 16091  |  回复: 31
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

刘小伟110

木虫 (小有名气)


[交流] vasp中有关自旋

最近看来一些文献,也做了一些优化,但是遇到了一些问题,求各位大虾帮忙解答。。
1.文献中,分子或者原子在MgO,V2O5,Al2O3的吸附都加入了自旋,All calculations utilize spin polarization,为什么加入自旋呢,这些氧化物貌似没有磁性的?
2.如果加入自旋的话,是不是只在INCAR里设置ISPIN=2就行了呢,MAGMOM,和NUPDOWN,就不用设置,用默认值呢?
3.假如OH和H共吸附在Al2O3上,得到稳定吸附,是不是只考虑总能最低(total energy负值最小)就好呢,还需要考虑单重态和三重态问题吗?
4.如果三重态不如单重态稳定,那么是不是就可以去单重态的计算为最终结构,最终能量呢?
5.如果是单重态,在INCAR里可不可以不设置,选择默认值呢?
求大虾指教、。。。。
回复此楼

» 收录本帖的淘帖专辑推荐

VASP and MS 第一性原理相关文档 vasp 学习贴 第一性
vasp计算细节 Emuch 锂电池小白笔记 VASP

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

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

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yxlllzj1314

铁虫 (初入文坛)



小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
28楼: Originally posted by 刘小伟110 at 2015-05-26 22:53:28
jpchou前辈,好久好久不见,我又来麻烦您了呢。。。
现在想做关于LDA+U的计算,有几个问题想要请教您一下,还望您能指导!!
1.>所说的+U贯穿整个计算,是不是晶胞优化的时候,就开始+U计算了呢,+U与不+U的区 ...

我测试过,只要有效u值一样,得到的结果几乎没差别。一般j值设置为有效u值的10%左右。
30楼2016-01-04 21:34:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 32 个回答
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
fzx2008: 金币+5, 谢谢指导 2014-03-14 18:31:57
刘小伟110: 金币+10 2014-03-14 19:37:50
1.文献中,分子或者原子在MgO,V2O5,Al2O3的吸附都加入了自旋,All calculations utilize spin polarization,为什么加入自旋呢,这些氧化物貌似没有磁性的?

> 吸附物若是有自旋,吸附表面後 "可能" 自旋還是存在

2.如果加入自旋的话,是不是只在INCAR里设置ISPIN=2就行了呢,MAGMOM,和NUPDOWN,就不用设置,用默认值呢?

> "一般情況" ISPIN=2,vasp 就會自動幫你找到最佳的自旋態

3.假如OH和H共吸附在Al2O3上,得到稳定吸附,是不是只考虑总能最低(total energy负值最小)就好呢,还需要考虑单重态和三重态问题吗?

> 對於未知系統,為求研究完整性,會去確認自旋態
   無自旋,單重,三重,那一個總能最低~
   就能確認穩定吸附態

4.如果三重态不如单重态稳定,那么是不是就可以去单重态的计算为最终结构,最终能量呢?

> 是

5.如果是单重态,在INCAR里可不可以不设置,选择默认值呢?

> ? 默認值? 那一個參數的默認值?
   vasp 的參數 MAGMOM 的默認值是  "每顆原子初使 spin moment = 1"
   不代表整個系統算完就是單重態

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

2楼2014-03-14 17:23:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘小伟110

木虫 (小有名气)


引用回帖:
2楼: Originally posted by jpchou at 2014-03-14 17:23:00
1.文献中,分子或者原子在MgO,V2O5,Al2O3的吸附都加入了自旋,All calculations utilize spin polarization,为什么加入自旋呢,这些氧化物貌似没有磁性的?

> 吸附物若是有自旋,吸附表面後 "可能&quo ...

谢谢您的回复,有几个问题。。。
1.吸附物若是有自旋,吸附表面後 "可能" 自旋還是存在
> 那么也就是说吸附的分子之间存在自旋了?
2. "一般情況" ISPIN=2,vasp 就會自動幫你找到最佳的自旋態
>就是设置ISPIN=2,MAGMOM不用设置,为默认值,如果要设置的话,slab+absorbate共50个,MAGMOM=50*1对不对呢?还有就是默认值MAGMOM= NIONS*1.0 forISPIN= 2,NIONS是不是就是原子总数呢?
3. ? 默認值? 那一個參數的默認值?vasp 的參數 MAGMOM 的默認值是  "每顆原子初使 spin moment = 1"
   不代表整個系統算完就是單重態
>单重态的话,默认值是不是NUPDOWN呢,如果我设置ISPIN=2,NUPDOWN=1的话是不是单重态呢?=3的话是不是就是三重态呢?
求指教。。。
3楼2014-03-14 20:19:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
franch: 金币+5, 应助指数+1, 谢谢回帖交流。。。 2014-03-24 17:34:23
谢谢您的回复,有几个问题。。。
1.吸附物若是有自旋,吸附表面後 "可能" 自旋還是存在
> 那么也就是说吸附的分子之间存在自旋了?

嗯~ 因為根據你說的  "文献中,分子或者原子在MgO,V2O5,Al2O3的吸附都加入了自旋"
如果文獻是可信賴且正確的
那~會加入自旋計算表示系統一定有磁性
但~如同你說說,氧化物表面應該是沒有磁性
所以推論系統磁性是發生在吸附物上面
就我所知過渡金屬(transition metal)原子吸附的確是會發生這樣的事

2. "一般情況" ISPIN=2,vasp 就會自動幫你找到最佳的自旋態
>就是设置ISPIN=2,MAGMOM不用设置,为默认值,如果要设置的话,slab+absorbate共50个,MAGMOM=50*1对不对呢?还有就是默认值MAGMOM= NIONS*1.0 forISPIN= 2,NIONS是不是就是原子总数呢?

是,是
假設N顆原子~ 就用 MAGMOM=N*1
不過這樣做,做跟沒做是一樣的,因為默認值就是這樣
如果有兩種原子,第一種 12 顆,第二種 7 顆
那 MAGMOM 可寫成
MAGMOM = 12*1  7*0   或
MAGMOM =  2*1  10*2  7*0   或
MAGMOM = 4*2  3*5  7*1  5*6
就是要有 19 個數字給它就是了

3. ? 默認值? 那一個參數的默認值?vasp 的參數 MAGMOM 的默認值是  "每顆原子初使 spin moment = 1"
   不代表整個系統算完就是單重態
>单重态的话,默认值是不是NUPDOWN呢,如果我设置ISPIN=2,NUPDOWN=1的话是不是单重态呢?=3的话是不是就是三重态呢?

NUPDOWN 沒有默認值,你沒在 INCAR 內加入它, vasp 是不會去使用它的
NUPDOWN 就是 spin up 跟 spin down 之電子的差
NUPDOWN = 0 -> S=0,    2S+1 = 1 = SINGLET
NUPDOWN = 1 -> S=1/2, 2S+1 = 2 = DOUBLET
NUPDOWN = 2 -> S=1,    2S+1 = 3 = TRIPLET

[ Last edited by jpchou on 2014-3-14 at 21:35 ]

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

4楼2014-03-14 21:08:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 有没有仍没收到信息的 +6 德尚中行 2026-08-27 7/350 2026-08-30 12:53 by Bxndjjd
[基金申请] 国社科又开始会评了,不知道这次命运如何 +3 雨打竹帘 2026-08-30 4/200 2026-08-30 12:35 by 白大多亮
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +8 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:40 by l0VvVHGBGRLv
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 11:32 by l0VvVHGBGRLv
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:19 by l0VvVHGBGRLv
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 07:39 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 06:22 by ZPa0EcMwuECS
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 02:09 by ZPa0EcMwuECS
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:58 by ZPa0EcMwuECS
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-30 01:03 by ZPa0EcMwuECS
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +9 苏东坡二世 2026-08-23 9/450 2026-08-29 14:41 by hustersqt
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 基金不中,共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
信息提示
请填处理意见