应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » vasp中有关自旋
51/1返回列表
查看: 14267  |  回复: 31
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

刘小伟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

» 猜你喜欢

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

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

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

查看全部散金贴
1楼2014-03-14 10:41:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘小伟110

木虫 (小有名气)
送红花一朵
引用回帖:
17楼: Originally posted by jpchou at 2014-04-04 00:46:26
>>就是在Vasp TST Tools上下载的脚本，比如关于过渡态插点的，可不可以直接用啊。。。还有就是如果想用CNEB搜素过渡态的话（现在只有asp自带的NEB），是不是要编译呢，不太会。。。

VTST 上面的 vtstscri ...

jpchou前辈，好久不见,你还好吧？最近忙的很，也不知在忙什么，难道我为去虫友互识蓄势。。。。唉，瞎忙。。。每次找你都得添点麻烦，过意不去。。。我想问您几个问题，麻烦您给解释一下，谢谢了。。。
1. 先前问过您，我在VTST上面下载的vtstscripts.tar.gz ，我该怎么用呢，就是要怎么安装呢，因为我没安装过vasp，所以也不知道把这个压缩文件放在哪里，主要是想用CINEB，都说NEB算出来的能垒偏小。。。这个问题蛮幼稚的，您就骂我小白吧。。
2. 我在搜寻过渡态，能量最高值进行频率计算，总能找到唯一虚频，得到的有点容易，我在想这个f/i的值有要求吗？是不是要大于某个值才算呢？是不是我弛豫原子数目有点少的原因（只弛豫吸附分子）
3. 我在计算中间体的时候，假如为H2O,中间体为HO,H。找到H和OH单一的最稳定吸附位，在各自稳定位计算最稳定的共吸附，但是优化完之后，OH和H又结合为H2O了，这是怎么回事呢，是不是我要改变一下共吸附时OH和H的相对高度呢？
4. 在什么情况下，需要考虑远程作用范德华力的影响呢，这个范德华力是不是就是周期性条件吸附分子之间的相互作用？
5. 最近遇到了一个灾难，绝对是一个灾难。。。要计算分子吸附过程的速率常数。吸附过程没有能垒也没有过渡态，怎么用TST理论计算呢，文献中把计算反应的能垒Ea改为ΔG(计算过态度的G与气态分子之间的G之差)，这该怎么理解呢，怎么计算呢
6. 还有就是，前辈，我用小木虫问您，这是否是合理且有礼貌的方式呢。。。比如邮箱什么呢，我尽量不扰您，要不，我不思考，也会显得太幼稚。。。哈哈
期待您的解答，再次感谢。。。
vasp中有关自旋
306EVTWPH8R@GUS8`_S8YG2.jpg
赞一下 回复此楼
18楼2014-05-02 19:17:13
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 32 个回答

jpchou

专家顾问 (著名写手)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫: 金币+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"
   不代表整個系統算完就是單重態
赞一下(3人) 回复此楼

» 本帖已获得的红花（最新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的回帖

jpchou

专家顾问 (著名写手)
★ ★ ★ ★ ★ ★
小木虫: 金币+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 ]
赞一下(6人) 回复此楼

» 本帖已获得的红花（最新10朵）

刘小伟110
4楼2014-03-14 21:08:35
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 32 个回答
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭