我是新来的，如果我缺乏适当的堆栈交换礼仪，请原谅我。所以，我想知道这里是否有人可以就我遇到的 Hartree-Fock 计算问题提供见解。我从头开始编写了一个代码，试图计算一维势能的 HF 能量和基础，该势势应该模拟核子-核子两体相互作用：

如果我为所谓的“无自旋”费米子运行代码，它会收敛并给出与文献中关于这种潜力的结果相似的结果。但是，如果我在问题中添加任何自旋自由度，SCF 过程不会收敛。相反，它到达一组两组或更多组特征向量，当用于计算新的 fock 矩阵时，它们最终将相互复制。我想相信这是因为我被困在 HF 泛函的局部最小值中。所以我的主要问题是：你们中有人遇到过这样的基本 HF 代码问题吗？考虑到在某些情况下我得到了很好的结果而在其他情况下没有结果，我觉得问题出在势的物理上，而不是在我的代码上。这在你眼里是合理的说法吗？