计算晶胞优化时固定角度失效

qli0218 · 发表于 Post on 2024-7-1 16:32:38

本帖最后由 qli0218 于 2024-7-1 16:42 编辑

如题，在输入文件中：KEEP_ANGLES T 设置角度固定，但是发现计算过程中，角度却是发生变化的。
所以请问这个是怎么回事，为什么固定不了？

另外还有一个小问题，有没有办法让abc三个晶胞常数一起发生变化，想让盒子的abc一直相等。谢谢

sobereva · 发表于 Post on 2024-7-1 18:28:40

你的晶胞定义不满足习俗
习俗应当a轴平行于X轴。我利用gview重新对晶轴做了变换以满足这一点并产生了新的CP2K输入文件，现在没有任何问题，始终能保持都是60度

qli0218 · 发表于 Post on 2024-7-1 19:43:17

sobereva 发表于 2024-7-1 18:28
你的晶胞定义不满足习俗
习俗应当a轴平行于X轴。我利用gview重新对晶轴做了变换以满足这一点并产生了新的C ...

非常感谢

Uus/pMeC6H4-/キ · 发表于 Post on 6 hour ago

类似问题在https://github.com/cp2k/cp2k/issues/3384有解释了，在CP2K里变胞优化的算法设计就没有考虑兼容不符合约定的盒子。这里的约定“第一个晶格矢量A平行于X轴、第二个晶格矢量B在XY平面，三个晶格矢量A, B, C构成下三角形矩阵”在最新版本的手册里有所强调，且最新版本的CP2K在读到不符合约定的晶格矢量定义时会弹警告。

除了建模时就做晶格矢量变换获得符合约定的盒子外，也可以改成用&CELL字段下ABC和ALPHA_BETA_GAMMA直接定义晶胞参数的写法搭配&COORD字段下SCALED T用分数坐标定义原子位置的设置，此时CP2K会自动构建符合约定的盒子并得到正确的原子笛卡尔坐标。

若要在变胞优化时保持盒子的对称性【此处“对称性”仅指晶胞参数的关系，如RHOMBOHEDRAL（菱方）有三个晶格矢量长度相等、三个夹角相等且非直角的关系，与需要spglib结合原子位置识别的空间群不是一回事】，则需同时在&CELL/SYMMETRY设定对称性、&CELL_OPT/KEEP_SYMMETRY启用对称性的约束。

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

[CP2K] 计算晶胞优化时固定角度失效

浏览过的版块