| 查看: 16096 | 回复: 31 | |||||||||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||||||||
[交流]
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 |
» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有19人回复
-
国社科又开始会评了，不知道这次命运如何
已经有9人回复
-
能否申诉？
已经有5人回复
-
为什么到现在没收到通知？
已经有5人回复
-
有没有仍没收到信息的
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- VASP自旋轨道耦合
已经有26人回复
- vasp SOC能带计算时无法分自旋
已经有5人回复
- vasp自旋极化计算时，初始设置自旋能否不被优化
已经有6人回复
- VASP计算-考虑自旋极化
已经有3人回复
- 产品中有水，老板让先旋蒸再那个精馏，为什么不能直接精馏，通过精馏除去水分呢
已经有11人回复
- 利用VASP计算材料的力学性质，用不用考虑 GGA+U 和 自旋 ？？？
已经有3人回复
- 求助vasp中的自旋轨道耦合
已经有3人回复
- 【求助】菜鸟问有关VASP作业提交的问题！
已经有21人回复
- 【求助】VASP怎么考虑自旋极化？
已经有14人回复
- 【交流】大家说说反铁磁的种类和特征吧
已经有19人回复
- 【求助】vasp检验Ti的赝势kpoint如何设置
已经有13人回复
- 【求助】vasp 离子驰豫与电子自旋的关系
已经有11人回复
- 【求助成功】VASP中，如何固定磁矩做结构优化
已经有14人回复
- 【求助】VASP考虑自旋轨道耦合的话，如何编译
已经有17人回复
- 【求助】doml3寻找磁矩
已经有10人回复
送红花一朵 |
jpchou前辈，谢谢您回复，我仔细想了一下，仍有几个问题想要请教您，望您能给与一定的解答。 1. 其实我想用Henkelman's group开发的小程式，在过渡态计算时，自动插点；用 Henkelman's group开发的 CNEB + vasp去寻找过渡态。后者要自己重新编译vasp吗？ 2. 在进行热力学计算时，转动熵或者转动配分函数中，都要涉及到转动惯量的计算，多原子分子的转动惯量该怎么计算呢，比如H2O的各个方向转动惯量Ix，Iy，Iz，是不是有个表呢 3. 前几天我计算了一个稳定气态分子的吸附，并找了一下频率，按理说，应该没有虚频的，难道不稳定，可是已经达到收敛标准了啊，怎么会有虚频呢，能不能消除呢？我是不是过分执着啦 4. 范德华力校正是不是就是传说中的偶极校正呢，如果一个非极性分子吸附在金属表面上，是分子与面的作用，还是分子周期性的作用呢？但是我原来看过，偶极校正就是最顶层版和最底层板是否对称。我该怎么理解呢 5. 要是算Bader电荷分析的话，我的版本有点低，算不了，能不能自己升级呢，还是要重新自己买vasp高版本呢 期待你的解答，再次膜拜，万分感谢。。。 |
20楼2014-05-14 21:54:09
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
fzx2008: 金币+5, 谢谢指导 2014-03-14 18:31:57
刘小伟110: 金币+10 2014-03-14 19:37:50
小木虫: 金币+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
|
谢谢您的回复，有几个问题。。。 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
★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
franch: 金币+5, 应助指数+1, 谢谢回帖交流。。。 2014-03-24 17:34:23
小木虫: 金币+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
10