using LinearAlgebra
b2_position = b_position ⊗ b_position
x1 = embed(b2_position, 1, x)
x2 = embed(b2_position, 2, x);
r = x1 - x2
d = abs.(diag(r.data)).^-1
d[d.==Inf] .= d[2]
data = sparse(diagm(0 => d))
H_coulomb = SparseOperator(b2_position, b2_position, data
运行显示倒数第二行出现标题的错误，是怎么回事？
内存不够了，你不认识英文吗？另外，你的例子无法运行。
前面还有程序，错误显示是说out of memory error在倒数第二行，是程序运行到倒数第二行时，内存不够了吗?
是吧，我猜是，你运行的时候自己监测下内存占用