24小时热门版块排行榜

>论坛更新日志 (4067)
>文献求助 (431)
>导师招生 (375)
>虫友互识 (369)
>休闲灌水 (248)
>考博 (168)
>硕博家园 (125)
>招聘信息布告栏 (118)
>博后之家 (96)
>论文投稿 (90)
>基金申请 (71)
>考研 (66)
>公派出国 (65)
>论文道贺祈福 (55)
>找工作 (45)
>教师之家 (40)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

剑魄诗魂

木虫 (正式写手)

1ST强帖: 1
应助: 6 (幼儿园)
贵宾: 0.06
金币: 2685.1
散金: 590
红花: 7
帖子: 405
在线: 109.9小时
虫号: 1061894
注册: 2010-07-21
性别: GG
专业: 理论和计算化学

[求助] 求教表面吸附的频率计算问题

最近在尝试计算表面吸附结构的自由能变化，所以需要计算频率以得到ZPE。但是在实际的计算中发现，计算时间非常久（超过4天），而且CHG, CHGCAR, WAVECAR变得非常大（合计超过50G），所以不得不停止了计算。
在此想跟同行求教一下是否我的具体设置上存在什么问题所以导致不容易计算？我的体系是在金属团簇表面的吸附，具体策略是固定所有金属原子（坐标文件如下）而不固定吸附物。
POSCAR：
1.00000000000000
40.0000000000000000 0.0000000000000000 0.0000000000000000
   0.0000000000000000 40.0000000000000000 0.0000000000000000
   0.0000000000000000 0.0000000000000000 40.0000000000000000
M C O
79    1    2
Direct
  0.4993046399203109  0.4992548039135158  0.4992924124318137 F F F
  0.4990485543665031  0.5438863108801399  0.5435908768283824 F F F
  0.4995033566807336  0.4546647224233595  0.4550508152545660 F F F

  .. (为节约空间省略若干行）
  0.5877448409594663  0.4996622676310523  0.4992604420477491 F F F
  0.4996940843718747  0.4109563411793573  0.4996687239395428 F F F
  0.4989078996841698  0.5876795921551263  0.4988993667433306 F F F
  0.5944357731254115  0.6390016216759307  0.5953946665121848 T T T
  0.5823633152745984  0.6624175346966013  0.5822536135436280 T T T
  0.6066869688169945  0.6157805664098532  0.6086636088716872 T T T
INCAR如下：
ISTART = 1
ICHARG = 1
PREC = Normal
ALGO = Normal
ISYM = 0 !0-no symmetry, 2-use symmetry

NSW = 1    ! 0 ionic steps
ENCUT = 400  ! energy cutoff for planewaves, ENMAX*1.3 when high

ISPIN = 2 !1-nonspin, 2-spin
GGA = PE ! PAW-PW91

LVDW = True ! DFT-D
IVDW = 11    ! zero damping DFT-D3 method of Grimme

EDIFF  = 1E-03 ! Energy convergency
EDIFFG = -0.03
LREAL  = .True.
ISMEAR = 1 ! Metallic
SIGMA  = 0.2

! frequency
POTIM = 0.02
IBRION = 5  ! frequency calculation
NFREE = 2
还望有经验的同行不吝赐教，我这个体系如何设置能够顺利计算表面吸附分子的振动频率。

回复此楼