计算化学公社

标题: 求助VASP频率计算不输出频率结果 [打印本页]

作者
Author: HydroKagami 时间: 2024-6-3 17:45
标题: 求助VASP频率计算不输出频率结果
本帖最后由 HydroKagami 于 2024-6-3 17:46 编辑

如题，初学第一性原理计算，本人用MS建模的slab，用VASP进行结构优化之后，把CONTCAR里的原子全部固定之后去跑频率计算，但是OUTCAR里没有频率结果，只是提示一个报错，本人才疏学浅不知如何解决。报错的行具体如下：
(, 下载次数 Times of downloads: 25)
INCAR如下:(使用VASPKIT生成)
Global Parameters
ISTART =  1          (Read existing wavefunction, if there)
ISPIN  =  1          (Non-Spin polarised DFT)
# ICHARG =  11       (Non-self-consistent: GGA/LDA band structures)
LREAL  = .FALSE.    (Projection operators: automatic)
# ENCUT  =  400       (Cut-off energy for plane wave basis set, in eV)
# PREC =  Accurate (Precision level: Normal or Accurate, set Accurate when perform structure lattice relaxation calculation)
LWAVE  = .TRUE.       (Write WAVECAR or not)
LCHARG = .TRUE.       (Write CHGCAR or not)
ADDGRID= .TRUE.       (Increase grid, helps GGA convergence)
LASPH  = .TRUE.       (Give more accurate total energies and band structure calculations)
PREC = Accurate    (Accurate strictly avoids any aliasing or wrap around errors)
# LVTOT  = .TRUE.    (Write total electrostatic potential into LOCPOT or not)
# LVHAR  = .TRUE.    (Write ionic + Hartree electrostatic potential into LOCPOT or not)
# NELECT =          (No. of electrons: charged cells, be careful)
# LPLANE = .TRUE.    (Real space distribution, supercells)
# NWRITE = 2          (Medium-level output)
# KPAR = 2          (Divides k-grid into separate groups)
# NGXF = 300       (FFT grid mesh density for nice charge/potential plots)
# NGYF = 300       (FFT grid mesh density for nice charge/potential plots)
# NGZF = 300       (FFT grid mesh density for nice charge/potential plots)

Frequence Calculations
NSW =  1          (number of ionic steps. Make it odd.)
ISMEAR =  0          (gaussian smearing method)
SIGMA  =  0.05       (please check the width of the smearing)
IBRION =  5          (frequence calculation)
POTIM  =  0.02       (displacement step)
NFREE  =  2          (displacement freedom)

POSCAR（VASPKIT固定）：
Generated by VASPKIT code
1.000000
10.3800001143999996 0.0000000000000000 0.0000000000000000
-5.1900000571999998 8.9893437903999995 0.0000000000000000
0.0000000000000000 0.0000000000000000 25.0000000000000000
C N V
18    3    1
Selective Dynamics
Direct
0.3039590519589040 0.2071975652400279 0.2967582088793958 F  F  F    C001
0.3047049371180406 0.0731165633495648 0.2962870468882421 F  F  F    C002
0.2265830277899737 0.2816734379680033 0.2966698900877863 F  F  F    C003
0.2232236174417393 0.4174747370204481 0.2960252811022558 F  F  F    C004
0.0831091356790135 0.2826594704593303 0.2966814477426620 F  F  F    C005
0.7967285725152284 0.0747058435044500 0.2969968525072611 F  F  F    C006
0.9311667383602088 0.2091603740785780 0.2970289078392365 F  F  F    C007
0.2994917012742503 0.7063788050961396 0.2957182146761603 F  F  F    C008
0.2994457283867668 0.5688070224426257 0.2960434225015924 F  F  F    C009
0.4370812377568701 0.7064491800596492 0.2960136298626033 F  F  F    C010
0.2238344146280367 0.7819777042950078 0.2950311989857927 F  F  F    C011
0.2277453825435529 0.9213540325947079 0.2957398069492900 F  F  F    C012
0.0844095463493495 0.7779763704612455 0.2957443898058305 F  F  F    C013
0.7985969213100964 0.7017608576559132 0.2967203881135836 F  F  F    C014
0.9326320942957227 0.7009588390702401 0.2962792452006333 F  F  F    C015
0.7242056138586624 0.7792369615200402 0.2965928912695702 F  F  F    C016
0.7232525682763451 0.9227621095302445 0.2965995650960039 F  F  F    C017
0.5884187651355139 0.7826492627554023 0.2959376698822165 F  F  F    C018
0.4502132915769860 0.2128551215805163 0.2950607412273050 F  F  F    N001
0.7908597525325524 0.2150078256468580 0.2957625472005763 F  F  F    N002
0.7928757125249969 0.5554492053600404 0.2950252694514350 F  F  F    N003
0.6775122496871936 0.3282287013109806 0.3018034407305694 F  F  F    V001

OSZICAR如下：
   N    E                   dE          d eps    ncg    rms       rms(c)
DAV: 1    0.810402468710E+03 0.81040E+03 -0.24967E+04 752 0.775E+02
DAV: 2    0.780869303921E+02 -0.73232E+03 -0.68898E+03 936 0.208E+02
DAV: 3 -0.178195384206E+03 -0.25628E+03 -0.24886E+03 888 0.122E+02
DAV: 4 -0.202893909942E+03 -0.24699E+02 -0.24494E+02  1008 0.578E+01
DAV: 5 -0.203478977321E+03 -0.58507E+00 -0.57816E+00 896 0.714E+00 0.248E+01
DAV: 6 -0.224962861125E+03 -0.21484E+02 -0.57788E+01  1000 0.401E+01 0.483E+01
DAV: 7 -0.182193177584E+03 0.42770E+02 -0.10641E+02  1000 0.135E+02 0.152E+01
DAV: 8 -0.178720129909E+03 0.34730E+01 -0.36694E+01 992 0.163E+01 0.917E+00
DAV: 9 -0.177916054479E+03 0.80408E+00 -0.35281E+00 936 0.719E+00 0.482E+00
DAV:  10 -0.178016025112E+03 -0.99971E-01 -0.10182E+00 936 0.515E+00 0.683E+00
DAV:  11 -0.177875555548E+03 0.14047E+00 -0.64304E+00 976 0.482E+00 0.575E+00
DAV:  12 -0.177768380693E+03 0.10717E+00 -0.34094E-01 896 0.156E+00 0.535E+00
DAV:  13 -0.177586412912E+03 0.18197E+00 -0.23599E-01  1024 0.187E+00 0.425E+00
DAV:  14 -0.177521298066E+03 0.65115E-01 -0.11695E+00 992 0.186E+00 0.223E+00
DAV:  15 -0.177487359253E+03 0.33939E-01 -0.27112E-01 920 0.831E-01 0.872E-01
DAV:  16 -0.177483034376E+03 0.43249E-02 -0.78727E-02  1000 0.339E-01 0.733E-01
DAV:  17 -0.177513992821E+03 -0.30958E-01 -0.31346E-01 976 0.474E-01 0.276E+00
DAV:  18 -0.177482818043E+03 0.31175E-01 -0.21105E-01 936 0.394E-01 0.190E-01
DAV:  19 -0.177483117102E+03 -0.29906E-03 -0.11392E-02 920 0.210E-01 0.253E-01
DAV:  20 -0.177484448803E+03 -0.13317E-02 -0.14991E-03 984 0.126E-01 0.404E-01
DAV:  21 -0.177487326044E+03 -0.28772E-02 -0.12514E-03 808 0.925E-02 0.314E-01
DAV:  22 -0.177488547501E+03 -0.12215E-02 -0.70546E-04 920 0.778E-02 0.299E-01
DAV:  23 -0.177488809592E+03 -0.26209E-03 -0.50188E-04  1000 0.395E-02 0.134E-01
DAV:  24 -0.177489904882E+03 -0.10953E-02 -0.75097E-04 936 0.539E-02 0.119E-01
DAV:  25 -0.177491174388E+03 -0.12695E-02 -0.15335E-04 776 0.304E-02 0.120E-01
DAV:  26 -0.177491635130E+03 -0.46074E-03 -0.42728E-05 816 0.195E-02 0.107E-01
DAV:  27 -0.177491833462E+03 -0.19833E-03 -0.82076E-05 936 0.173E-02 0.336E-02
DAV:  28 -0.177492084338E+03 -0.25088E-03 -0.17304E-05 752 0.134E-02 0.852E-03
DAV:  29 -0.177492306559E+03 -0.22222E-03 -0.11054E-05 704 0.821E-03 0.851E-03
DAV:  30 -0.177492383472E+03 -0.76914E-04 -0.85308E-06 704 0.464E-03 0.134E-02
DAV:  31 -0.177492412566E+03 -0.29094E-04 -0.16131E-06 584 0.371E-03 0.701E-03
DAV:  32 -0.177492430961E+03 -0.18396E-04 -0.12967E-06 568 0.233E-03 0.211E-03
DAV:  33 -0.177492437144E+03 -0.61830E-05 -0.16959E-07 552 0.153E-03 0.334E-03
DAV:  34 -0.177492440250E+03 -0.31061E-05 -0.15612E-07 520 0.114E-03 0.130E-03
DAV:  35 -0.177492442009E+03 -0.17589E-05 -0.50929E-08 520 0.870E-04 0.119E-03
DAV:  36 -0.177492442535E+03 -0.52575E-06 -0.23728E-08 552 0.531E-04
1 F= -.17749244E+03 E0= -.17747755E+03  d E =-.297929E-01

作者
Author: nullzy 时间: 2024-6-3 17:50
原子全固定了，频率当然算不了了。

作者
Author: HydroKagami 时间: 2024-6-3 17:55

nullzy 发表于 2024-6-3 17:50
原子全固定了，频率当然算不了了。

谢谢老师回复，我以前是用MS算第一性原理，之前同样的体系的自由能是用老师给的脚本，那里面要求用ase固定，也能得出结果，我就照葫芦画瓢固定了这些原子。那请问老师这是完全不需要固定还是要把xy平面给固定了，只留下高度进行频率计算？

作者
Author: nullzy 时间: 2024-6-3 20:08

HydroKagami 发表于 2024-6-3 17:55
谢谢老师回复，我以前是用MS算第一性原理，之前同样的体系的自由能是用老师给的脚本，那里面要求用ase固 ...

不知道你是不是做理论催化的，领域不一样可能计算方式有区别，做催化的时候一般只算吸附中间体跟气态分子的自由能，slab的ZPE跟TS我们一般视为0。

作者
Author: HydroKagami 时间: 2024-6-3 20:42

nullzy 发表于 2024-6-3 20:08
不知道你是不是做理论催化的，领域不一样可能计算方式有区别，做催化的时候一般只算吸附中间体跟气态分子 ...

我是做电催化的，大一学生，刚入门计算，用dmol3算了一段时间开始深入了解，在学习vasp。之前用MS计算自由能是完全交给老师写的脚本当黑箱操作的，但是明确的知道，算台阶的时候第一步单独载体不吸附中间体的自由能也是要固定后交给脚本计算得出的

作者
Author: vasp冲冲 时间: 2024-6-3 21:45

HydroKagami 发表于 2024-6-3 20:42
我是做电催化的，大一学生，刚入门计算，用dmol3算了一段时间开始深入了解，在学习vasp。之前用MS计算自 ...

俺楼上说的来吧，只算吸附后固定slab的频率就行了

作者
Author: 乐平 时间: 2024-6-3 22:08
回忆一下中学物理，频率是不是在振动里学的？

如果你把所有原子都固定了，那怎么振动呢？

作者
Author: 努力学计算 时间: 2024-7-24 08:50

nullzy 发表于 2024-6-3 17:50
原子全固定了，频率当然算不了了。

请问老师，我是固定了底层原子来搜的过渡态，成功搜出来了。然后用的这个固定底层原子的模型结构计算频率，在outmol文件里也显示finished successfully，但频率分析的时候，强度那一列都是0，这是因为固定了底层，所以强度都为0嘛，但outmol文件最后提示成功完成，所以我不太明白是不是与底层固定有关，计算频率是不是需要都不固定呀？就是只固定底层原子也不行

作者
Author: 努力学计算 时间: 2024-7-24 08:50

乐平发表于 2024-6-3 22:08
回忆一下中学物理，频率是不是在振动里学的？

如果你把所有原子都固定了，那怎么振动呢？

欢迎光临计算化学公社 (http://bbs.keinsci.com/)