计算化学公社

标题: 用CP2K做金刚石单轴(z)应变的仿真求助 [打印本页]

作者
Author: yb6231290 时间: 2021-11-12 21:02
标题: 用CP2K做金刚石单轴(z)应变的仿真求助
大家好，我想咨询下，用CP2K做金刚石单轴(z)应变的仿真，应该是调整晶格参数，如做z轴拉伸金刚石supercell3 3 3的10*10*10的晶胞，是将晶胞参数调整为10*10*10.1，然后再跑GEO_ OPT吗？还是做CELL_OPT，然后逐渐调整EXTERNAL_PRESSURE来实现啊？

作者
Author: sobereva 时间: 2021-11-13 00:38
看你是想考察晶胞z轴变化特定距离时原子坐标怎么变化，还是想考察特定外压下结构怎么变化，取决于你

作者
Author: yb6231290 时间: 2021-11-13 14:22

sobereva 发表于 2021-11-13 00:38
看你是想考察晶胞z轴变化特定距离时原子坐标怎么变化，还是想考察特定外压下结构怎么变化，取决于你

我是想计算晶体受到外压力或外拉力作用下，晶体的能量变化，现在用的external_pressure，发现能量没有什么变化，所以在考虑是不是需要做应变

作者
Author: sobereva 时间: 2021-11-14 03:27

yb6231290 发表于 2021-11-13 14:22
我是想计算晶体受到外压力或外拉力作用下，晶体的能量变化，现在用的external_pressure，发现能量没有什 ...

数值设更大再试

作者
Author: Aridea 时间: 2021-11-14 03:35
做之前最好做一下cell opt得到relaxed cell

作者
Author: yb6231290 时间: 2021-11-15 15:13

Aridea 发表于 2021-11-14 03:35
做之前最好做一下cell opt得到relaxed cell

谢谢回复，请问relaxed cell后再加geo，是可以直接上0.01 0.02 ... 0.05这样做吗？

作者
Author: yb6231290 时间: 2021-11-15 15:14
本帖最后由 yb6231290 于 2021-11-15 15:23 编辑

sobereva 发表于 2021-11-14 03:27
数值设更大再试

我做pressure已经加拉力到20GPa了，可以看到晶格常数的变化，但是里面不同位置，不同方向的缺陷(空位和掺杂)完全没有变化，不知道是什么原因啊，我用的PBE泛函，我看网上说这里应该把potential GTH改成ALL，您看是否需要呢？

作者
Author: sobereva 时间: 2021-11-16 01:11

yb6231290 发表于 2021-11-15 15:14
我做pressure已经加拉力到20GPa了，可以看到晶格常数的变化，但是里面不同位置，不同方向的缺陷(空位和掺 ...

跟potential完全没直接关系。
我不知道你说的没变化是怎么回事，晶格尺寸都变了，不可能原子坐标都一点不变，顶多是肉眼不好分辨的程度

作者
Author: yb6231290 时间: 2021-11-16 12:54
本帖最后由 yb6231290 于 2021-11-16 14:30 编辑

sobereva 发表于 2021-11-16 01:11
跟potential完全没直接关系。
我不知道你说的没变化是怎么回事，晶格尺寸都变了，不可能原子坐标都一点 ...

老师您好，是这样的，我想做111面的金刚石拉伸，同时观察金刚石内部不同方向NV的能量变化，这个在其他文献中的总结中可知：沿着111的NV方向和非111的NV方向随着z轴的拉力或压力，逐渐出现能量的分叉，文献结果如图所示。([img][/img]见附件图片)，文献中使用vasp计算，计算中采用了672个原子，0应变时，二者的能量差在0.0018eV，应变到了0.01能量差大概在0.1eV，未使用k点。
但是我自己做的沿111面拉伸，原子数也为672个，未使用k点，cutoff为400，泛函和基组分别为：GTH-PBE和DZVP-MOLOPT-SR-GTH，同时添加了dftd3色散，OT收敛精度1e-6，计算cellopt，external_pressure从0-20GPa，做了10组，其中10GPa时金刚石z轴应变约为0.01，但发现这两个nv方向的能量差始终保持在0.0019-0.0023eV，所以我想请问下到底是什么原因导致了cp2k和vasp差异如此的大？

作者
Author: xjc 时间: 2024-11-25 09:07
我想用cp2k做一个第一性原理分子动力学的拉伸，改变z轴方向的长度来进行拉伸，但是我不知道怎么写输入文件,能指教一下吗？

作者
Author: huangyh 时间: 2024-11-26 20:32

xjc 发表于 2024-11-25 09:07
我想用cp2k做一个第一性原理分子动力学的拉伸，改变z轴方向的长度来进行拉伸，但是我不知道怎么写输入文件, ...

您好，请问这个问题解决了嘛，我也想学习一下。

作者
Author: xjc 时间: 2024-12-12 14:10

huangyh 发表于 2024-11-26 20:32
您好，请问这个问题解决了嘛，我也想学习一下。

还没有呢，你有什么想法吗

作者
Author: amoszhou 时间: 2025-4-11 14:19

xjc 发表于 2024-12-12 14:10
还没有呢，你有什么想法吗

请问实现了吗

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