采用Dimer算法精算过渡态结构时遇到CHGCAR路径报错问题

CedricLiang

各位老师好，我在计算NH3分子过渡态时先使用CI-NEB方法找到一个近似鞍点，但是此时有三个虚频，我想使用Dimer方法精算一下，使用VTST脚本neb2dim.pl脚本生成了MODECAR，提交作业后出现”forrtl: severe (29): file not found, unit 18, file /public1/home/scf2226/ltf/NH34/DIM/02/CHGCAR“的错误
MODECAR以及我的INCAR代码文件如下

1. SYSTEM=NH3
   
2. ISTART=0
   
3. ICHARG=2
   
4. 
   
5. LCLIMB=.TRUE.
   
6. IMAGES=3
   
7. SPRING=-5
   
8. 
   
9. NELM=60
   
10. EDIFF=1E-7
    
11. ENCUT=520
    
12. ALGO=Accurate
    
13. LDIAG=.FALSE.
    
14. 
    
15. IBRION=3
    
16. ISIF=2
    
17. POTIM=0
    
18. IOPT=2
    
19. NSW=100
    
20. EDIFFG=-0.05
    
21. ICHAIN=2
    
22. 
    
23. ISMEAR=0
    
24. SIGMA=0.05
    
25. 
    
26. PREC=Accurate
    
27. LREAL=.FALSE.
    
28. LWAVE=.FALSE.
    
29. LCHARG=.FALSE.
    
30. 
    
31. ISPIN=2
    
32. MAGMOM=4*0.6
    
33. 
    
34. DdR=5E-3
    
35. DRotMax=4
    
36. DFNMax=1.0
    
37. DFNMin=0.01

复制代码

脚本生成的MODECAR代码如下

1.     3.8957010937E-05     1.1247933573E-04     0.0000000000E+00
   
2.     7.1452132712E-05     1.0632939466E-04     0.0000000000E+00
   
3.    -9.7420986364E-06     1.3559805593E-04     0.0000000000E+00
   
4.    -4.6789626517E-04    -1.4331690046E-03    -9.9999883902E-01
   

复制代码