24小时热门版块排行榜    

查看: 1632  |  回复: 47
本帖产生 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的网格,消除简并度
人生苦短,抓紧科研,抓紧找对象。
33楼2011-05-26 15:37:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

胡同学

金虫 (正式写手)

【答案】应助回帖


youzhizhe(金币+1): 谢谢交流。 2011-05-26 19:23:43
我用的是GGA(PBE),我计算的结果和我看到的一篇论文得到的数据很接近,估计我计算正确的可能性比较大。
人生苦短,抓紧科研,抓紧找对象。
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:
看文献上经常说以解离的形式吸附,我想问一下您,怎样去判断解离呢?除了键长发生变化外,还有没有其它的方法?

通常是按键长的变化。
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的基态是非磁性的。即使你做的是自旋极化的计算,最后计算收敛到了非磁态。你检查你进行的自旋极化计算的结果,看一下总的磁矩是多少。
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 的主题更新
信息提示
请填处理意见