vasp对不同结构的水滑石结构优化差别较大

风风光风光

各位老师好，我想问一下我分别使用水滑石结构的FeOOH和CoOOH构建了Co:Fe=1:1的水滑石结构，在结构优化后结构的差异比较大，想问问老师们这个的原因是什么，这两个结构做吸附能的计算结果会有差别吗？

可以看到不同的水滑石块优化出来的结构差别比较的大，Fe上加Co出来的结构歪七扭八的，而且Co上加Fe体系的能量要比Fe上加Co的能量低50eV，这是否说明Co上加Fe的结构才是对的？

不懂的实在太多，请老师们不吝赐教。

pal

第一张看起来原子好像乱了，2-4好像差别不大吧，3和4旋转个90度目测跟2差不多

风风光风光

pal 发表于 2023-5-23 11:42
第一张看起来原子好像乱了，2-4好像差别不大吧，3和4旋转个90度目测跟2差不多

图一是图二优化后的状态，图四是图三优化后的状态，然后二和三因为分别是用FeOOH和CoOOH做的，所以尺寸不一样。我觉得可能是块尺寸不同，所以导致优化后的1乱了4没乱

风风光风光

图1的能量我算出来是-293.6027702eV，磁矩是32.2144，图4的能量算出来是-341.33232ev，磁矩是33.1959

zmjsce

风风光风光 发表于 2023-5-31 16:16
图1的能量我算出来是-293.6027702eV，磁矩是32.2144，图4的能量算出来是-341.33232ev，磁矩是33.1959

理论上来说，水滑石应该有磁矩吗？
如果不应该有磁矩，那可能要考虑通过MAGMOM设置金属元素的磁性方向（通过正负制定）以将磁性收敛到正确的状态。
此外，还需要考虑DFT+U参数。

风风光风光

zmjsce 发表于 2023-5-31 23:13
理论上来说，水滑石应该有磁矩吗？
如果不应该有磁矩，那可能要考虑通过MAGMOM设置金属元素的磁性方向（ ...

我在MP的晶体数据库里看到Fe、Co的羟基氧化物上是有磁性的，所以在做二元水滑石的时候就延续之前的参数设置，加了ISPIN和MAGMOM的设置。
关于DFT+U的参数，我之前有看到说如果是对于能量的计算的话这个加和不加影响不大，我现在只算一个吸附能，所以没加，这个在磁性结构计算中会有影响吗？

看的帖子是这篇
vasp中用PBE算含有d和f电子的原子必须要加U吗？
http://bbs.keinsci.com/forum.php ... 6&fromuid=42950
(出处: 计算化学公社)