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

栋栋家的小猪猪 · 发表于 Post on 2021-11-16 09:00:29

在当前的计算中设置了自旋，加了自旋多重度这个参数，然后每次设置分别为1/3/5... 发现随着自旋多重度的增加，体系能量一直在降低，要想得到最低的能量，是要一直测试下去吗？还是有其他比较好的方法能直接获得体系能量最低时的多重度的数值？

ljh123 · 发表于 Post on 2024-12-20 16:57:39

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确实没找到这个脚本，但是找到了其他测试文件的相关脚本。请问这种情况是什么原因呢？非常感谢老师的回答

zjxitcc · 发表于 Post on 2024-12-20 15:36:06

本帖最后由 zjxitcc 于 2024-12-20 15:46 编辑

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

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

ljh123 · 发表于 Post on 2024-12-20 15:28:08

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

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

15689900879@163 · 发表于 Post on 2024-11-29 22:14:42

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

好嘞，谢谢您~

丁越 · 发表于 Post on 2024-11-29 19:43:27

15689900879@163 发表于 2024-11-29 13:07
老师您好，我也遇到自旋多重度选择的问题了，体系含有一个三线态氧气分子，和一个OH自由基。我在计算化学 ...

不用单独优化比较能量，直接当前结构做个单点比较就行。

15689900879@163 · 发表于 Post on 2024-11-29 13:07:59

丁越发表于 2021-11-16 09:54
你的体系要是有机分子，一般情况下自旋多重度为1；除非涉及到过渡金属，你才需要设置不同的自旋多重度来判 ...

老师您好，我也遇到自旋多重度选择的问题了，体系含有一个三线态氧气分子，和一个OH自由基。我在计算化学群里，有位老师建议我分别设置自旋多重度为2或4。请问我是在CP2K设置GEO_OPT任务，看同一结构进行优化的能量，哪个更低，选择哪个自旋多重度吗？还是跑AIMD任务呢，比较ener文件的势能大小呢？

ljlljl · 发表于 Post on 2024-7-6 11:02:48

丁越发表于 2021-11-16 09:54
你的体系要是有机分子，一般情况下自旋多重度为1；除非涉及到过渡金属，你才需要设置不同的自旋多重度来判 ...

老师你好，第一种方法是指只需要指定磁矩并结合UKS，自选多重度随便设定一个值吗？

TJU_RYC · 发表于 Post on 2024-4-19 17:01:26

丁越发表于 2021-11-16 09:54
你的体系要是有机分子，一般情况下自旋多重度为1；除非涉及到过渡金属，你才需要设置不同的自旋多重度来判 ...

您好，想请教一下我有一个体系含有磁性元素掺杂MoS2（17个Mo,1个Ni,36个S)，需要设置磁矩吗？我设置的自旋多重度除了1其他都不收敛

不会扣篮的后卫 · 发表于 Post on 2022-1-27 23:47:00

丁越发表于 2022-1-27 16:44
Materials project （https://materialsproject.org/）

谢谢丁老师！

丁越 · 发表于 Post on 2022-1-27 16:44:46

不会扣篮的后卫发表于 2022-1-27 01:01
丁老师，请问一下，MP是代表什么？

Materials project （https://materialsproject.org/）

不会扣篮的后卫 · 发表于 Post on 2022-1-27 06:05:22

sobereva 发表于 2022-1-27 01:40
大抵是指Material Project（https://www.materialsproject.org/）

好的，谢谢卢老师！

sobereva · 发表于 Post on 2022-1-27 01:40:49

不会扣篮的后卫发表于 2022-1-27 01:01
丁老师，请问一下，MP是代表什么？

大抵是指Material Project（https://www.materialsproject.org/）

不会扣篮的后卫 · 发表于 Post on 2022-1-27 01:01:08

丁越发表于 2021-11-16 09:54
你的体系要是有机分子，一般情况下自旋多重度为1；除非涉及到过渡金属，你才需要设置不同的自旋多重度来判 ...

丁老师，请问一下，MP是代表什么？

丁越 · 发表于 Post on 2021-11-16 19:46:42

栋栋家的小猪猪发表于 2021-11-16 16:40
不用啦已经解决啦

UKS=LSD，说的是一个意思

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

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

浏览过的版块