Gaussian和CP2K都完全可以
Gaussian的话通常需要你去猜迁移路径上各个过渡态（可能不止一个）结构然后优化过渡态。CP2K用NEB的话，可以自己给出迁移前后的结构并插点，如果中间有中间体、多个过渡态，在NEB能量变化变化曲线上会有所体现（之后再对每个基元过程用CI-NEB或dimer找轨道），因此适合不好估计中间有哪些过渡态/中间体的情况。也可以二者结合，比如CP2K做NEB确定过渡态粗略位置，作为Gaussian的初猜结构用opt=TS优化过渡态。