最近在尝试计算表面吸附结构的自由能变化，所以需要计算频率以得到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
还望有经验的同行不吝赐教，我这个体系如何设置能够顺利计算表面吸附分子的振动频率。 返回小木虫查看更多