sobereva 发表于 2024-4-18 07:17
X光衍射的晶体结构里的氢位置本来往往就不靠谱，仔细看
实验测定分子结构的方法以及将实验结构用于量子 ...

明白了，谢谢老师

默默 发表于 2024-4-17 16:02
sob老师好，打扰了。我用CP2K-2023.1对自己的体系(三氨基胍盐酸盐与2,5-二羟基对苯二甲醛反应)进行了计算 ...

X光衍射的晶体结构里的氢位置本来往往就不靠谱，仔细看
实验测定分子结构的方法以及将实验结构用于量子化学计算需要注意的问题
http://sobereva.com/569
另外，用合适的杂化泛函优化准确度更高，有些氢转移可能是纯泛函的SIE问题造成的虚假的现象。


看优化的目的

sobereva 发表于 2024-4-15 08:20
单独对一层正电骨架做优化还说得过去，堆叠结构还只考虑正电骨架部分，就完全不能接受了

可以拿一层正电 ...

sob老师好，打扰了。我用CP2K-2023.1对自己的体系(三氨基胍盐酸盐与2,5-二羟基对苯二甲醛反应)进行了计算。我先做了收敛性测试，然后在PBE-D3(BJ)/ DZVP-MOLOPT-SR-GTH下进行了结构优化，输入和输出文件我发在楼上了。我想向您请教两个问题：
（1）我试了下固定角度和不固定角度进行优化，结果两个结构的氢都向氯转移了，这和模型化合物的单晶不太一样，不知道是不是我计算设置的问题？
（2）另外固定晶胞角度和不固定晶胞角度进行优化哪种才是正确的呢？我看两种结构差别还挺大的，从优化最后一步的能量（还没算单点）来看固定晶胞角度得到的结构能量反而更低

sobereva 发表于 2024-4-15 08:20
单独对一层正电骨架做优化还说得过去，堆叠结构还只考虑正电骨架部分，就完全不能接受了

可以拿一层正电 ...

明白了，感谢sob老师，我尝试尝试

单独对一层正电骨架做优化还说得过去，堆叠结构还只考虑正电骨架部分，就完全不能接受了

可以拿一层正电骨架，计算原子电荷或者表面静电势，估计阴离子容易结合在哪（如果本身就容易判断，就不需要这一步），然后对堆叠结构中相应位置放入阴离子，之后跑动力学（可以模拟退火），之后再做结构优化，大概率能得到能量最低结构。

对正离子的结构进行优化是程序是？