zhj45318403 发表于 2026-5-19 14:31
十分感谢老师。如果可以的话还想咨询一个与CI-NEB有关的问题。我目前计算200多步 始终没有收敛到我设定的 ...

1. 你可以自己在插点的文件夹里面看，也可以自己写脚本，vtst自带的有没有我没具体了解过。如果是把所有插点的能量一起罗列vtst有自带的脚本可以实现。能不能等dimer跑完要看你所有插点的收敛情况，比如如果是-0.1的收敛限，最高点到了-0.1，其他点0.2，0.3的样子也可以停掉，如果是-0.5的收敛限我建议等他跑完。如果跑飞了那就是没找到合适的路径，如果过渡态跑不出来大概率是初末态的机构有问题，或者是精度不够。
2. 参数我没看出来有什么明显的问题，EDIFF可以调到-7试试，不一定能跑出来。

nuonuo123456 发表于 2026-5-18 21:15
对的，我就是这个意思，而且可以clneb和dimer联用一般效率会更高。而且做neb的过程中要仔细观察结构，有 ...

十分感谢老师。如果可以的话还想咨询一个与CI-NEB有关的问题。我目前计算200多步 始终没有收敛到我设定的值，并且有越来越跑飞的趋势。我看到您说的和dimer联用，其它帖子也多次见到这个方法，我想额外的问题是：
1、是否可以罗列出所有轮次的力收敛结果，找到切线力接近于0 能量最高 最大力比较小的几个点，拿出来做dimer，而不用等全部ci-neb跑完
2、以下是我进行CI-NEB其中一个模型的INCAR ，如您有时间，能否帮我见到看看我一下关于CI-NEB的参数的设置是否合理？十分感激！


# --- Global Parameters ---
SYSTEM = Pd
ISTART = 0        
ISPIN  = 2        
PREC   = Normal  
ENCUT  = 450      
LWAVE  = .FALSE.   
LCHARG = .TRUE.
LAECHG = .FALSE.
ICHARG = 2   
LREAL  = Auto     
IVDW = 11
ISYM = 0   # 取消对称性

# --- Electronic Relaxation ---
ALGO   = Normal
ISMEAR = 0         
SIGMA  = 0.05
NELM   = 150      
EDIFF  = 1E-06     

# --- Ionic Relaxation ---
NSW    = 200      
IBRION = 3
POTIM = 0        
ISIF   = 2         
EDIFFG = -0.05     


# VTST Transitional state search
IOPT = 1
ICHAIN = 0
LCLIMB = .TRUE.
IMAGES = 5
SPRING = -5

# --- DFT+U Settings (The Soul of CeO2) ---
LDAU      = .TRUE.
LDAUTYPE  = 2
LDAUL     = -1  3 -1 -1 -1
LDAUU     = 0.0 5.0 0.0 0.0 0.0
LDAUJ     = 0.0 0.0 0.0 0.0 0.0
LMAXMIX   = 6  
LASPH  = .TRUE.   
MAGMOM = 94*0.0 27*0.0 1.2 4*0.0 1.2 8*0.0 1.2 2*0.0 1.2 2*0.0 0.0 0.0 4*0.0

# --- Parallel Performance ---
KPAR   = 4
NCORE  = 8   

zhj45318403 发表于 2026-5-18 20:57
谢谢老师回复！我理解您的意思是：其实结构优化后振动分析没有太大的必要，如果CI-NEB能跑通说明初、末态 ...

对的，我就是这个意思，而且可以clneb和dimer联用一般效率会更高。而且做neb的过程中要仔细观察结构，有时候它可能第一个点的能量比初态还低但是结构变化很小，或者是再次优化完能量或者结构差距很小这时候基本上可以默认是极小态，要灵活运用，多看结构。

nuonuo123456 发表于 2026-5-18 20:38
1. 个人认为频率分析即便是找极小态放开五层也太多了，你看文献几乎做甲烷吸附很少有找极小态的吧
2. 感觉 ...

谢谢老师回复！我理解您的意思是：其实结构优化后振动分析没有太大的必要，如果CI-NEB能跑通说明初、末态就是局域极小，如果跑不通就再优化就是了，没有单独去找极小态的必要。

仔细想来阅读文献确实也没有看到谁专门去找极小态的，谢谢您的指导！

1. 个人认为频率分析即便是找极小态放开五层也太多了，你看文献几乎做甲烷吸附很少有找极小态的吧
2. 感觉不太有必要
3. 手动选择毕竟好
另外，如果是甲烷吸附活化你肯定是以表面反应为主，我认为可以结构优化之后直接找过渡态，因为你过渡态CLNEB的过程中如果初末态不是极小态也可以再去重新优化。如果每一个吸附物种都这样做频率，付出和收获完全不成正比。