sobereva 发表于 2022-5-3 11:50 明白了，谢谢老师 |
LuLu-Chem 发表于 2022-5-1 22:40 我不知道你的初始结构怎么构建的 如果只是按原子数比例，随便往里面堆积单原子，那肯定不行。比如有些C彼此离得太近，可能直接就变成C2分子了。这正如同优化单个分子，也不能随便就把原子按数目随便画上去就了事 |
sobereva 发表于 2022-4-30 07:27 老师，您说的对，我多跑了一会儿，换了不同的体系，发现还是有原子跑出来了 请教老师，像这种原子跑出来的情况一般考虑有哪些原因导致的呢？ |
LuLu-Chem 发表于 2022-4-29 11:07 这不能说明什么，很可能只是巧合 怀疑是截断能方面的问题的话，应当同时提升CUTOFF和REL_CUTOFF再试，或者做收敛性测试 |
sobereva 发表于 2022-4-29 09:33 嗯，明白了老师，我后面不再用这个关键词了。我的目的是实现快速降温，但又要保持体系的稳定性，所以后面考虑采用梯度降温的方式跑。 老师，关于原子跑出盒子边界的问题我调小了REL_CUTOFF值，就不再有原子跑出盒子了，不知道这样是否合理？ |
| TEMPERATURE_ANNEALING是在控温时之前的温度要乘的系数，不是你退火过程的目标温度，你对这个关键词的用法理解完全错误 |
吞木木 发表于 2022-4-28 18:44 感谢老师的回复，我去尝试看看。另外，老师发的帖子对我有很大启发，向您学习！ |
LuLu-Chem 发表于 2022-4-28 17:33 以下仅个人见解，不一定对，仅供参考。 假如你真的想实现退火的功能，可以考虑ANNEALING这个关键词，不过这个关键词仅仅适用于NVE体系，表示的是下一步与上一步温度的比值。这个关键词人为控制性过强，一般不考虑作为数据收集部分。 https://manual.cp2k.org/cp2k-8_2 ... N/MD.html#ANNEALING 而且个人认为温控这一块cp2k做的并不是很好，可能lammps或者vasp更适合你研究这种需要退火的体系 |
吞木木 发表于 2022-4-28 16:36 感谢老师的回复，因为研究的是固态体系，然后需要实现退火这一步，所以就尝试用了这个关键词，刚开始学习CP2K有很多不懂的地方，就想尝试下看看结果。想请教老师1273.15K条件下，初始的固体结构应该已经变为熔融态，在这一步应该用什么系综比较合适呢？ 还有一个问题，需要退火的这一步，研究体系应该是由熔融态转变为固态的过程，在这一步应该选用什么系综合适呢？退火温度应该通过什么方式设置呢？ |
|
https://manual.cp2k.org/cp2k-8_2-branch/CP2K_INPUT/MOTION/MD.html#TEMPERATURE_ANNEALING Specifies the rescaling factor for the external temperature.This scheme works only for the Langevin ensemble. 这个TEMPERATURE_ANNEALING关键词，是你这么用的？我很迷惑。 这里不是明确说明，这个关键词只适用LANGEVIN系综？你用的NVT啊 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-19 19:57 , Processed in 0.248009 second(s), 25 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites