24小时热门版块排行榜    

CyRhmU.jpeg
查看: 14268  |  回复: 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的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见