计算化学公社

标题: 求助：VASP怎么算一个单分子的频率（无周期性边界条件） [打印本页]
作者
Author: wangy0822    时间: 2023-1-9 20:55
标题: 求助：VASP怎么算一个单分子的频率（无周期性边界条件）
大家好，请教大家一个问题，我现在想用VASP计算一个单分子的频率，没有周期性边界条件的那种，就和Gaussian里面的计算结果做一下对比，原因是审稿人建议的。请问大家怎么计算（泛函用B3LYP-D3(BJ)）,INCAR里应该怎么写B3LYP-D3(BJ)的泛函关键词。
作者
Author: 卡开发发    时间: 2023-1-10 02:44
本帖最后由 卡开发发 于 2023-1-10 03:53 编辑
  1. PREC=N
  2. ENCUT=...
  3. ...
  4. # multipole corrections
  5. # Ref: https://www.vasp.at/wiki/index.php/Monopole_Dipole_and_Quadrupole_corrections
  6. LDIPOL=T
  7. IDIPOL=4

  9. # dispersion corrections（DFT-D3(BJ)）
  10. IVDW=12

  12. # B3LYP(VWN3)
  13. # Ref: https://www.vasp.at/wiki/index.php/List_of_hybrid_functionals#B3LYP
  14. LHFCALC=T
  15. GGA=B3
  16. AEXX=0.2
  17. AGGAX=0.72
  18. AGGAC=0.81
  19. ALDAC=0.19

  21. # vibrational analysis (phonon)
  22. # Ref: https://www.vasp.at/wiki/index.php/IBRION#IBRION.3D5_and_6:_second_derivatives.2C_Hessian_matrix_and_phonon_frequencies_.28finite_differences.29.
  23. IBRION=5
复制代码


注：
1、因为Gaussian中用的B3LYP是基于VWN3的所以使用GGA=B3，如果是VWN5对照考虑使用GGA=B5
2、IVDW=12我检查了源码，GGA=B3的情况指向LEXCH=11，理论上说程序有内置的DFTD3参数，手册上没提到包含B3LYP的参数，这点我建议慎重测试，你可以在OUTCAR看到打印的参数。
3、我不确定后续版本的考虑IVDW的时候是否可以使用线性响应方法计算声子，也许也不能，这里采用IBRION=5通过有限差分计算声子。
4、分子情况使用IDIPOL=4校正整体偶极矩，盒子应当使用足够大且为立方盒子，可以测试体系能量随着盒子扩大的收敛的情况来确定盒子的大小。
作者
Author: kyuu    时间: 2023-1-13 19:41
不可能 vasp没有静电平板



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3