Uus/pMeC6H4-/キ 发表于 2025-9-26 16:37
如果了解LDH类晶体结构中非羟基氧与羟基氧的键长特征，那是可以辅助判断，但是我不熟。键角与二面角的问 ...

好的，非常感谢，现在我遇到的问题是SCF收敛不了，CP2K试了，QE也试了，不知道下一步咋办

renzhogn424 发表于 2025-9-26 15:22
非常感谢！！！前阵子出差一下忘了回复了，不好意思！
根据您的建议，我调整了下，看起来感觉也挺合理 ...

如果了解LDH类晶体结构中非羟基氧与羟基氧的键长特征，那是可以辅助判断，但是我不熟。键角与二面角的问题结合VSEPR指出的四面体四配位、八面体六配位等结构化学基础知识相对容易快速判断一点。

从化学式看现实的晶体结构里氢的数量以及硫酸根的数量本来就不确定（编辑：这个老帖子也有提到），更别说晶体学无序的位置了，我怀疑就算把表征方法从测不准氢原子位置的X光衍射换成测得准的中子衍射之类（http://sobereva.com/569）也只能得到一堆分数占据的氢原子。

这个cif文件似乎是来自AMCSD的，和COD-9015059基本一致。晶胞内有16个铁原子和32个氧原子，可以大致区分为下图所示红绿蓝高亮的三组：

如果参考layered double hydroxide里氢氧根的成键情况，图中绿色的16个氧偏离与之最靠近（或者说可视化软件通过一定几何依据生成的辅助视觉效果的“成键”，见http://sobereva.com/414）的三个铁原子构成的平面较远，更有可能是四面体四配位氧的位置，刚好补氢变成羟基之后之后组成也符合+3价铁对应的FeOOH。如果往蓝色的另外16个氧上补氢，就容易产生较小的Fe-O-H键角而不太合适。

补氢可以用GaussView的Builder - Add Valence功能，一个个点击上述标绿色的氧，然后再调一调键长键角二面角使之指向空隙。