计算化学公社

标题: 晶体不同方向加压优化结构相同求助 [打印本页]

作者
Author: 破晓时的曙光 时间: 2022-3-4 14:21
标题: 晶体不同方向加压优化结构相同求助
我想计算一个晶体xyz三个方向分别加压力，晶体结构如何变化，采用cp2K优化晶体结构的时候，设置external pressure这个三乘三的矩阵，能够实现不同方向加压吗？例如改变11，22，33矩阵元实现在x，y，z三个方向加不同的压力，为啥我x加0.1 Gpa，y加 0.1 Gpa，z加0.1 Gpa这三种情况优化的结构都是一样的呢？我的PRESSURE_TOLERANCE设置的是1 bar，比我的压力小多了。难道说这个external pressure不能实现不同方向加压吗？

作者
Author: sobereva 时间: 2022-3-4 18:33
加大压力数量级，让不同方向压力差更大再试。本身固体就极难被压缩

作者
Author: 破晓时的曙光 时间: 2022-3-7 15:53

sobereva 发表于 2022-3-4 18:33
加大压力数量级，让不同方向压力差更大再试。本身固体就极难被压缩

我加大了压力，从0.1 Gpa加大到0.5 Gpa，xyz不同方向加压优化结果才会不同，那这个结果能说明压力较小的时候，不同方向加压导致的晶体结构相同吗？

作者
Author: sobereva 时间: 2022-3-7 18:10

破晓时的曙光发表于 2022-3-7 15:53
我加大了压力，从0.1 Gpa加大到0.5 Gpa，xyz不同方向加压优化结果才会不同，那这个结果能说明压力较小的 ...

说明压力不够大的时候，压力差导致的形变可忽略不计。这跟具体体系特征有关，固体普遍难以压缩

作者
Author: 破晓时的曙光 时间: 2022-3-8 15:31
本帖最后由破晓时的曙光于 2022-3-8 15:38 编辑

sobereva 发表于 2022-3-7 18:10
说明压力不够大的时候，压力差导致的形变可忽略不计。这跟具体体系特征有关，固体普遍难以压缩

我用CP2K,x加0.1 Gpa，y加 0.1 Gpa，z加0.1 Gpa这三种情况优化的结构都是一样的,但是换成CASTEP，x加0.1 Gpa，y加 0.1 Gpa，z加0.1 Gpa这三种情况优化的结构又都不一样，比较符合直觉，如x方向加压，a轴变短了，y轴加压，b轴变短了。这好奇怪。是不是CP2K优化的算法在加小的压力时不太好，我就算改小了优化的步长也不行，而CASTEP的结果更合理。

作者
Author: sobereva 时间: 2022-3-9 04:56

破晓时的曙光发表于 2022-3-8 15:31
我用CP2K,x加0.1 Gpa，y加 0.1 Gpa，z加0.1 Gpa这三种情况优化的结构都是一样的,但是换成CASTEP，x加0.1 ...

改改压力收敛限再试

另外，CP2K和CASTEP用的计算方法没有可比性，基函数、赝势方面完全不同，也可能是理论方法层面而非数值设定层面的原因，也未必CASTEP就是合理的
原理上不可能存在“算法在加小的压力时不太好”这种问题

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