24小时热门版块排行榜    

查看: 16077  |  回复: 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的回帖

刘小伟110

木虫 (小有名气)


送红花一朵
引用回帖:
15楼: Originally posted by jpchou at 2014-03-24 02:26:24
1. 对于LREAL的设置,FALSE是在倒空间进行计算,;Auto 应该是在实空间进行吧(也有说是自动选择空间,不太清楚),但是如果用FALSE进行,vasp会提示用Auto比较好,但是好多文献都用平面波基组,而非积分球法啊?这 ...

jpchou前辈好,好久不见,最近打篮球脑震荡了,就没思考问题。。。哈哈。。又来给你找麻烦了。。。
前面的问题都得到了解决,还想问您几个问题。。
1. 使用vasp tools的脚本文件,是不是可以不进行编译,直接复制到vasp的根目录下进行呢?
> ? 看不懂?
   腳本~也就是 shell script 對吧
   shell script 是不用編輯
   只需 chmod 即可
>>就是在Vasp TST Tools上下载的脚本,比如关于过渡态插点的,可不可以直接用啊。。。还有就是如果想用CNEB搜素过渡态的话(现在只有asp自带的NEB),是不是要编译呢,不太会。。。
2. 现在用NEB搜寻过渡态,如果在反应物和产物之间线性插入4个image01 02 03 04,没有找到鞍点的话,继续插点,是在02 03之间的POSCAR之间插点呢,还是在02 03优化后的CONTCAR之间插点呢?
3. 在这4个点中总有一个能量最高点,是不是要把每一次插点的能量最高点计算虚频呢?那岂不是要计算很多。。。判断是不是过渡态的标准有哪些?
4.有人说看OUTCAR最后一步坐标下的
  tangential force (eV/A)        -0.044742
  left and right image   0.601145  0.557231 A
如果tangential force为0则为鞍点,这对不对呢?我有一个tangential force不为0,但是再插点,没有找到比这更高测度能量了。。。还有就是left and right image什么意思呢?
5. 由于01 04比较接近与固定的两个端点,优化必定会向两端点进行,那岂不是鞍点只会出现在02 03中吗?存不存在插不出过渡态呢?
期待你的回复,再次感谢。。。。
16楼2014-04-03 22:23:42
已阅   回复此楼   关注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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 7/350 2026-08-29 17:31 by 4FFAWE8HcgUD
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 17:30 by 4FFAWE8HcgUD
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 ASdOkHsho7FD 2026-08-28 8/400 2026-08-29 17:22 by 4FFAWE8HcgUD
[基金申请] 为什么到现在没收到通知? +4 tannykie 2026-08-29 4/200 2026-08-29 17:14 by lmz0216
[基金申请] 有没有仍没收到信息的 +5 德尚中行 2026-08-27 5/250 2026-08-29 16:16 by manplx
[论文投稿] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 13:51 by jCd0dEvKHShX
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 11:49 by jCd0dEvKHShX
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 11:29 by jCd0dEvKHShX
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +14 晴天加油 2026-08-26 15/750 2026-08-29 11:06 by biily
[基金申请] 投票:  有多少人是今天查系统知道结果的? +16 爱看书的可乐 2026-08-26 18/900 2026-08-29 10:18 by winsaint
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见