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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » vasp中有关自旋
51/1返回列表
查看: 14278  |  回复: 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的回帖

jpchou

专家顾问 (著名写手)
★ ★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
刘小伟110: 金币+1 2014-05-02 18:03:14
franch: 金币+5, 谢谢回帖交流。。。 2014-05-05 23:07:11
引用回帖:
16楼: Originally posted by 刘小伟110 at 2014-04-03 22:23:42
jpchou前辈好，好久不见，最近打篮球脑震荡了，就没思考问题。。。哈哈。。又来给你找麻烦了。。。
前面的问题都得到了解决，还想问您几个问题。。
1. 使用vasp tools的脚本文件，是不是可以不进行编译，直接复制 ...

>>就是在Vasp TST Tools上下载的脚本，比如关于过渡态插点的，可不可以直接用啊。。。还有就是如果想用CNEB搜素过渡态的话（现在只有asp自带的NEB），是不是要编译呢，不太会。。。

VTST 上面的 vtstscripts.tar.gz 是用 perl 寫的，跟 shell script 是不太一樣的東西~ 不過這兩者都是不需要編譯的


2. 现在用NEB搜寻过渡态，如果在反应物和产物之间线性插入4个image01 02 03 04，没有找到鞍点的话，继续插点，是在02 03之间的POSCAR之间插点呢，还是在02 03优化后的CONTCAR之间插点呢？

這....... 這要看你的結果是長什麼樣，看你的系統是長什麼樣
我是覺得那邊能量高~就插那一邊


3. 在这4个点中总有一个能量最高点，是不是要把每一次插点的能量最高点计算虚频呢？那岂不是要计算很多。。。判断是不是过渡态的标准有哪些？

判斷過渡態我也只知道要算 phonon
這可能要請更專業做過渡態研究的人?碇v講

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什么意思呢？

不是很肯定~所以我暫不回答
但幫你找到有人在討論這個
http://muchong.com/html/201006/2147544.html

5. 由于01 04比较接近与固定的两个端点，优化必定会向两端点进行，那岂不是鞍点只会出现在02 03中吗？存不存在插不出过渡态呢？
期待你的回复，再次感谢。。。。

所以只有四點是不保險的做法
七點或九點~這�?的奇數點會比較能看出趨勢
過渡態的找尋是件很辛苦的事
做 NEB 不見得一定能找到過渡態
我是覺得 potential energy surface (PES) 比較能找到過渡態
但 PES 計算很吃資源~~
赞一下(2人) 回复此楼

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

刘小伟110
17楼2014-04-04 00:46:26
 已阅   回复此楼   关注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 个回答
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭