计算化学公社

标题: 如何用multiwfn制作CP2K搜索过渡态（包含初始、过渡态两个结构）的输入文件？ [打印本页]

作者
Author: lao7 时间: 2024-8-13 22:06
标题: 如何用multiwfn制作CP2K搜索过渡态（包含初始、过渡态两个结构）的输入文件？
各位老师好！
想用multiwfn制作一个cp2k的过渡态计算的文件。但是，我读了之前的有些帖子，说是需要制作一个初始结构、一个过渡态结构。但是，multiwfn制作dimer搜索过渡态inp输入文件时，似乎只能包含一个结构。问题是：
1、如果multiwfn支持dimer计算过渡态，那么两个结构怎么通过multiwfn弄到一个文件里呢？
2、如果把初始结构稍微摆弄一下，使其结构倾向于具有过渡态结构特征，那么是否这么一个文件能否用来搜索过渡态？我的意思都是用multiwfn制作输入文件（自己不会写输入文件，太多文字内容）
谢谢各位！

作者
Author: sobereva 时间: 2024-8-14 05:14
CP2K支持dimer方法基于单个初猜结构找过渡态（用法类似于Gaussian中的opt=TS），也支持NEB类型的方法产生反应路径，其中CI-NEB可以同时得到过渡态。NEB类型的方法最起码需要提供反应物和产物的结构（但并不需要非得精确），也可以再额外提供反应路径中间一个或多个点的结构。方法原理参考下文
过渡态、反应路径的计算方法及相关问题
http://sobereva.com/44
用dimer时只能提供一个结构作为过渡态的初猜，越接近实际要找的过渡态越好。但对于初猜结构要求得比起opt=TS更为宽松，哪怕是在与过渡态相连的极小点附近也不是不可以，但显然优化到想要的过渡态的成功几率远低于基于较好的过渡态初猜、总耗时也高得多。

PS：北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）里有充分、非常详细的用dimer和NEB类方法找过渡态、产生反应路径的讲解和实例，包括用Multiwfn产生输入文件的全面细致的说明。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)