我正在尝试计算 ZnS 晶格的马德隆常数。方法如下: 晶格是面心立方,基为和。Madelung 常数是 其中的电荷(在这种情况下除以相邻块的数量), r_j到固定离子的距离。Evjen 的方法告诉我们,我们需要形成中性电荷块才能获得良好的收敛性。
我选择从两个立方体开始,原点放置在两个立方体共有的面的中心(所以我有一个大小为 1x1x2 的块)。然后,两个获得更高的精度,我在每个方向添加两个立方体(所以在第二次迭代中,我有一个大小为 3x3x4 的块)。
问题是这个方法收敛到,而它应该是。 这是正确的块扩展吗? 我已经尝试过其他块,我得到了相同的值。
我没有包含代码,因为我认为问题出在所选的块中。