计算化学公社

标题: 在cp2k计算过程中的自旋多重度怎么设置呀？ [打印本页]

作者
Author: 栋栋家的小猪猪 时间: 2021-11-16 09:00
标题: 在cp2k计算过程中的自旋多重度怎么设置呀？
在当前的计算中设置了自旋，加了自旋多重度这个参数，然后每次设置分别为1/3/5... 发现随着自旋多重度的增加，体系能量一直在降低，要想得到最低的能量，是要一直测试下去吗？还是有其他比较好的方法能直接获得体系能量最低时的多重度的数值？

作者
Author: 丁越 时间: 2021-11-16 09:54
本帖最后由丁越于 2021-11-16 09:55 编辑

你的体系要是有机分子，一般情况下自旋多重度为1；除非涉及到过渡金属，你才需要设置不同的自旋多重度来判断哪个自旋多重度下能量最低；
当涉及到过渡金属含有磁性时，必须要指定磁矩。&MAGNETIZATION 结合UKS，磁矩大小可以根据MP上的数值设定，程序会自动收敛到合理值。此时自旋设不设定都没有必要了。
还有一种方法只适合于UKS，即不需要设定自旋多重度，让程序去自发寻找，该关键词在这有介绍https://manual.cp2k.org/cp2k-8_2 ... _RELAX_MULTIPLICITY

作者
Author: 栋栋家的小猪猪 时间: 2021-11-16 14:47

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

好的谢谢我试一试

作者
Author: 栋栋家的小猪猪 时间: 2021-11-16 16:04

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

您好，我想再问一下，就是如果打开UKS，那我先前使用的LSD是不是就可以注销啦？以及您所说的MP是不是自旋多重度的缩写呀？小白不太懂，盼回复，谢谢

作者
Author: 栋栋家的小猪猪 时间: 2021-11-16 16:40

栋栋家的小猪猪发表于 2021-11-16 16:04
您好，我想再问一下，就是如果打开UKS，那我先前使用的LSD是不是就可以注销啦？以及您所说的MP是不是自旋 ...

不用啦已经解决啦

作者
Author: 丁越 时间: 2021-11-16 19:46

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

UKS=LSD，说的是一个意思

作者
Author: 不会扣篮的后卫 时间: 2022-1-27 01:01

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

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

作者
Author: sobereva 时间: 2022-1-27 01:40

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

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

作者
Author: 不会扣篮的后卫 时间: 2022-1-27 06:05

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

好的，谢谢卢老师！

作者
Author: 丁越 时间: 2022-1-27 16:44

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

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

作者
Author: 不会扣篮的后卫 时间: 2022-1-27 23:47

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

谢谢丁老师！

作者
Author: TJU_RYC 时间: 2024-4-19 17:01

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

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

作者
Author: ljlljl 时间: 2024-7-6 11:02

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

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

作者
Author: 15689900879@163 时间: 2024-11-29 13:07

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

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

作者
Author: 丁越 时间: 2024-11-29 19:43

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

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

作者
Author: 15689900879@163 时间: 2024-11-29 22:14

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

好嘞，谢谢您~

作者
Author: ljh123 时间: 2024-12-20 15:28

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

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

作者
Author: zjxitcc 时间: 2024-12-20 15:36
本帖最后由 zjxitcc 于 2024-12-20 15:46 编辑

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

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

作者
Author: ljh123 时间: 2024-12-20 16:57

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

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