【答案】应助回帖

臾兮，臾兮，何时是归期！

31楼2011-05-26 13:53:28

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

kmw.8668

金虫 (正式写手)

归期守望者

应助: 9 (幼儿园)
金币: 1731.6
散金: 217
红花: 2
帖子: 668
在线: 235.9小时
虫号: 1008177
注册: 2010-04-29
专业: 凝聚态物性 II ：电子结构

引用回帖:

谢谢，可以把把你计算H2的INCAR 文件和POSCAR 文件发过来让我看看吗？为什麽我计算的和你的还有那麽大差距？我用的是GGA，你呢？

臾兮，臾兮，何时是归期！

32楼2011-05-26 13:56:12

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

胡同学

金虫 (正式写手)

应助: 4 (幼儿园)
金币: 1286.7
散金: 293
红花: 2
沙发: 1
帖子: 433
在线: 669.5小时
虫号: 1090636
注册: 2010-09-06
性别: GG
专业: 原子和分子物理

【答案】应助回帖

★ ★
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的回帖

胡同学

金虫 (正式写手)

应助: 4 (幼儿园)
金币: 1286.7
散金: 293
红花: 2
沙发: 1
帖子: 433
在线: 669.5小时
虫号: 1090636
注册: 2010-09-06
性别: GG
专业: 原子和分子物理

【答案】应助回帖

★
youzhizhe(金币+1): 谢谢交流。 2011-05-26 19:23:43

我用的是GGA(PBE)，我计算的结果和我看到的一篇论文得到的数据很接近，估计我计算正确的可能性比较大。

人生苦短，抓紧科研，抓紧找对象。

34楼2011-05-26 15:39:56

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

1ST强帖: 13
应助: 241 (大学生)
金币: 25701.7
散金: 602
红花: 166
帖子: 3782
在线: 873.8小时
虫号: 1007127
注册: 2010-04-27
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

★
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

至尊木虫 (职业作家)

1ST强帖: 13
应助: 241 (大学生)
金币: 25701.7
散金: 602
红花: 166
帖子: 3782
在线: 873.8小时
虫号: 1007127
注册: 2010-04-27
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

★
youzhizhe(金币+1): 谢谢交流。 2011-05-26 19:24:03
kmw.8668(金币+2): 2011-05-26 21:35:01

引用回帖:

H2的基态是非磁性的。即使你做的是自旋极化的计算，最后计算收敛到了非磁态。你检查你进行的自旋极化计算的结果，看一下总的磁矩是多少。

36楼2011-05-26 16:24:13

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

kmw.8668

金虫 (正式写手)

归期守望者

应助: 9 (幼儿园)
金币: 1731.6
散金: 217
红花: 2
帖子: 668
在线: 235.9小时
虫号: 1008177
注册: 2010-04-29
专业: 凝聚态物性 II ：电子结构

引用回帖:

Originally posted by valenhou001 at 2011-05-26 16:19:04:
通常是按键长的变化。

但是我在初始计算的时候，并不知道吸附的分子到底是以解离形式还是以分子的形式吸附，那我就认为它是以分子形式吸附吗？通常键长超过多大范围就说分子解离了呢？

臾兮，臾兮，何时是归期！

37楼2011-05-26 17:27:53

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

kmw.8668

金虫 (正式写手)

归期守望者

应助: 9 (幼儿园)
金币: 1731.6
散金: 217
红花: 2
帖子: 668
在线: 235.9小时
虫号: 1008177
注册: 2010-04-29
专业: 凝聚态物性 II ：电子结构

引用回帖:

Originally posted by valenhou001 at 2011-05-26 16:19:04:
通常是按键长的变化。

一个虫友说还可以按照吸附能，电子分布图来分析吸附的分子是否发生解离，请问你知道如何通过这两项进行分析吗？

臾兮，臾兮，何时是归期！

38楼2011-05-26 17:32:51

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

kmw.8668

金虫 (正式写手)

归期守望者

应助: 9 (幼儿园)
金币: 1731.6
散金: 217
红花: 2
帖子: 668
在线: 235.9小时
虫号: 1008177
注册: 2010-04-29
专业: 凝聚态物性 II ：电子结构

引用回帖:

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的回帖

胡同学

金虫 (正式写手)

应助: 4 (幼儿园)
金币: 1286.7
散金: 293
红花: 2
沙发: 1
帖子: 433
在线: 669.5小时
虫号: 1090636
注册: 2010-09-06
性别: GG
专业: 原子和分子物理

【答案】应助回帖

我其实也不大清楚，我用的是Medea的平台，设置好参数后，有些参数是默认的，它会自动得出那些结果。我也有些不是很清楚，但我知道那样设置没错