在cp2k计算过程中的自旋多重度怎么设置呀？

15689900879@163

丁越 发表于 2024-11-29 19:43
不用单独优化比较能量，直接当前结构做个单点比较就行。

好嘞，谢谢您~

ljh123

丁越 发表于 2024-11-29 19:43
不用单独优化比较能量，直接当前结构做个单点比较就行。

老师请问下，如果我想用AIMD考察一个金属络合物的化学变化，反应过程中自旋多重度会从5变成7。可是AIMD只能指定一个初始自旋多重度，并不能考虑势能面交叉的情况，可以请教下您该怎么跑AIMD考察这个化学反应吗？

zjxitcc

ljh123 发表于 2024-12-20 15:28
老师请问下，如果我想用AIMD考察一个金属络合物的化学变化，反应过程中自旋多重度会从5变成7。可是AIMD只 ...

如果能垒非常低（比如小于5 kcal/mol）或无能垒，你可以使用非绝热从头算分子动力学，程序在每条轨迹中的每一步会计算几个电子态。反之，AIMD在你能跑得动的轨迹时长内 难以跑出目标化学反应（考虑多个电子态更拖慢速度），但可采用常见研究反应机理的手段，即分别计算两个自旋多重度下的化学反应路径（过渡态、中间体、产物等），不涉及AIMD，然后寻找两个energy profile的MECP结构，最后分析反应如何发生。

ljh123

zjxitcc 发表于 2024-12-20 15:36
如果能垒非常低（比如小于5 kcal/mol）或无能垒，你可以使用非绝热从头算分子动力学，程序在每条轨迹中的 ...

老师我了解了下非绝热动力学，发现CP2K可以借助NewtonX计算surface hopping 动力学。但是目前有两个问题（1）surface hopping 这个关键词基本都出现在激发态计算里，我不清楚我这种仅仅是电子转移导致的自旋多重度变化适不适合用NewtonX计算。（2）我安装完NewtonX2.6后用自带的inp-testnx.pl进行测试，选择了test53(Dynamics, Non-adiabatic with cioverlap-od)。但是出现结果报错：sh: /share/home/tj85153/apps/newtonx-cs/bin/run_cioverlap_cp2k.pl: No such file or directory
run-cp2k.pl: Error in run_cioverlap_cp2k.pl
似乎应该是缺少run_cioverlap_cp2k.pl脚本，我在/apps/newtonx-cs/bin确实没找到这个脚本，但是找到了其他测试文件的相关脚本。请问这种情况是什么原因呢？非常感谢老师的回答