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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 表面吸附
484/5返回列表上一页12345下一页
查看: 1632  |  回复: 47
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 3 个 1ST强帖 ，点击这里进行查看

kmw.8668

金虫 (正式写手)

归期守望者
引用回帖:
Originally posted by 胡同学 at 2011-05-26 10:54:01:
我计算一个H2分子的能量（也是放到10x10x10的晶胞中）不加ISPIN＝2这一项，得出的能量值为-6.699907ev,加spin则为-6.685360ev，H原子加spin为-1.103034ev,不加spin为-0.062981ev。我认为加spin的时候，结果中包含 ...

一个老师刚刚跟我说H是一种特殊的元素，计算的时候需要加ISPIN＝2这一项的
赞一下 回复此楼
臾兮，臾兮，何时是归期！
31楼2011-05-26 13:53:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kmw.8668

金虫 (正式写手)

归期守望者
引用回帖:
Originally posted by 胡同学 at 2011-05-26 10:54:01:
我计算一个H2分子的能量（也是放到10x10x10的晶胞中）不加ISPIN＝2这一项，得出的能量值为-6.699907ev,加spin则为-6.685360ev，H原子加spin为-1.103034ev,不加spin为-0.062981ev。我认为加spin的时候，结果中包含 ...

谢谢，可以把把你计算H2的INCAR 文件和POSCAR 文件发过来让我看看吗？为什麽我计算的和你的还有那麽大差距？我用的是GGA，你呢？
赞一下 回复此楼
臾兮，臾兮，何时是归期！
32楼2011-05-26 13:56:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

胡同学

金虫 (正式写手)

【答案】应助回帖

★ ★
kmw.8668(金币+5): 谢谢 2011-05-26 16:54:45
youzhizhe(金币+2): 谢谢交流。 2011-05-26 19:23:29
# SCF input for VASP
# Note that VASP uses the FIRST occurence of a keyword
SYSTEM = H  (P1) ~ 73 minimized: Untitled (VASP)

    PREC = Normal
   ENCUT = 250.000
  IBRION = -1
     NSW = 0
  NELMIN = 2
   EDIFF = 1.0e-05
  EDIFFG = -0.02
VOSKOWN = 1
  NBLOCK = 1
    NELM = 60
    ALGO = Normal (blocked Davidson)
   ISPIN = 2
MAGMOM = 0.106
  INIWAV = 1
  ISTART = 0
  ICHARG = 2
   LWAVE = .TRUE.
  LCHARG = .TRUE.
ADDGRID = .FALSE.
   LVTOT = .TRUE.
    LELF = .TRUE.
  ISMEAR = 1
   SIGMA = 0.2
   LREAL = .FALSE.
RWIGS = 0.32

H  (P1) ~ 73 minimized: Untitled (VASP)
    1.0
            7.14484571 0.00000000 0.00000000
            0.00000013 5.70609702 0.00000000
            0.00000017 0.00000017 7.85990351
  1
Direct
        0.50000000 0.50000000 0.50000000

# SCF input for VASP
# Note that VASP uses the FIRST occurence of a keyword
SYSTEM = (H)2  (P1) ~ 142 minimized: Untitled (VASP)

    PREC = Normal
   ENCUT = 250.000
  IBRION = -1
     NSW = 0
  NELMIN = 2
   EDIFF = 1.0e-05
  EDIFFG = -0.02
VOSKOWN = 1
  NBLOCK = 1
    NELM = 60
    ALGO = Normal (blocked Davidson)
   ISPIN = 2
MAGMOM = 2*0.005
  INIWAV = 1
  ISTART = 0
  ICHARG = 2
   LWAVE = .TRUE.
  LCHARG = .TRUE.
ADDGRID = .FALSE.
   LVTOT = .TRUE.
    LELF = .TRUE.
  ISMEAR = 1
   SIGMA = 0.2
   LREAL = .FALSE.
RWIGS = 0.32

(H)2  (P1) ~ 142 minimized: Untitled (VASP)
    1.0
            6.98452495 0.00000000 0.00000000
            0.00000011 5.90303441 0.00000000
            0.00000011 0.00000009 5.90303441
  2
Direct
        0.52170012 0.50000000 0.50000000
        0.41429988 0.50000000 0.50000000

H原子计算时我用的是10*8*11的网格，消除简并度
赞一下(1人) 回复此楼
人生苦短，抓紧科研，抓紧找对象。
33楼2011-05-26 15:37:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

胡同学

金虫 (正式写手)

【答案】应助回帖

★
youzhizhe(金币+1): 谢谢交流。 2011-05-26 19:23:43
我用的是GGA(PBE)，我计算的结果和我看到的一篇论文得到的数据很接近，估计我计算正确的可能性比较大。
赞一下(1人) 回复此楼
人生苦短，抓紧科研，抓紧找对象。
34楼2011-05-26 15:39:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

【答案】应助回帖

★
kmw.8668(金币+2): 谢谢 2011-05-26 16:50:58
youzhizhe(金币+1): 谢谢交流。 2011-05-26 19:23:52
引用回帖:
Originally posted by kmw.8668 at 2011-05-26 10:37:50:
看文献上经常说以解离的形式吸附，我想问一下您，怎样去判断解离呢？除了键长发生变化外，还有没有其它的方法？

通常是按键长的变化。
赞一下(1人) 回复此楼
35楼2011-05-26 16:19:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

【答案】应助回帖

★
youzhizhe(金币+1): 谢谢交流。 2011-05-26 19:24:03
kmw.8668(金币+2): 2011-05-26 21:35:01
引用回帖:
Originally posted by 胡同学 at 2011-05-26 10:54:01:
我计算一个H2分子的能量（也是放到10x10x10的晶胞中）不加ISPIN＝2这一项，得出的能量值为-6.699907ev,加spin则为-6.685360ev，H原子加spin为-1.103034ev,不加spin为-0.062981ev。我认为加spin的时候，结果中包含 ...

H2的基态是非磁性的。即使你做的是自旋极化的计算，最后计算收敛到了非磁态。你检查你进行的自旋极化计算的结果，看一下总的磁矩是多少。
赞一下(1人) 回复此楼
36楼2011-05-26 16:24:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kmw.8668

金虫 (正式写手)

归期守望者
引用回帖:
Originally posted by valenhou001 at 2011-05-26 16:19:04:
通常是按键长的变化。

但是我在初始计算的时候，并不知道吸附的分子到底是以解离形式还是以分子的形式吸附，那我就认为它是以分子形式吸附吗？通常键长超过多大范围就说分子解离了呢？
赞一下 回复此楼
臾兮，臾兮，何时是归期！
37楼2011-05-26 17:27:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kmw.8668

金虫 (正式写手)

归期守望者
引用回帖:
Originally posted by valenhou001 at 2011-05-26 16:19:04:
通常是按键长的变化。

一个虫友说还可以按照吸附能，电子分布图来分析吸附的分子是否发生解离，请问你知道如何通过这两项进行分析吗？
赞一下 回复此楼
臾兮，臾兮，何时是归期！
38楼2011-05-26 17:32:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kmw.8668

金虫 (正式写手)

归期守望者
引用回帖:
Originally posted by 胡同学 at 2011-05-26 15:37:11:
# SCF input for VASP
# Note that VASP uses the FIRST occurence of a keyword
SYSTEM = H  (P1) ~ 73 minimized: Untitled (VASP)

    PREC = Normal
   ENCUT = 250.000
  IBRION = -1
     NSW = 0
...

请问MAGMOM＝2*0.005和RWIGS＝0.32这两个参数你是怎末确定它们的值的？
赞一下 回复此楼
臾兮，臾兮，何时是归期！
39楼2011-05-26 18:54:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

胡同学

金虫 (正式写手)

【答案】应助回帖

我其实也不大清楚，我用的是Medea的平台，设置好参数后，有些参数是默认的，它会自动得出那些结果。我也有些不是很清楚，但我知道那样设置没错
赞一下 回复此楼
人生苦短，抓紧科研，抓紧找对象。
40楼2011-05-26 20:16:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kmw.8668 的主题更新
484/5返回列表上一页12345下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭