求助：VASP怎么算一个单分子的频率（无周期性边界条件）

wangy0822

大家好，请教大家一个问题，我现在想用VASP计算一个单分子的频率，没有周期性边界条件的那种，就和Gaussian里面的计算结果做一下对比，原因是审稿人建议的。请问大家怎么计算（泛函用B3LYP-D3(BJ)）,INCAR里应该怎么写B3LYP-D3(BJ)的泛函关键词。

kyuu

不可能 vasp没有静电平板

卡开发发

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
   
8. 
   
9. # dispersion corrections（DFT-D3(BJ)）
   
10. IVDW=12
    
11. 
    
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
    
20. 
    
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校正整体偶极矩，盒子应当使用足够大且为立方盒子，可以测试体系能量随着盒子扩大的收敛的情况来确定盒子的大小。