“…Because the interaction |x| −d (d > 0) diverges when x → 0, the ground-state even wavefunction φ (+) of (31) The occupancies and their natural orbitals are determined by the following integral eigenequation ∞ −∞ ρ(x, y)v s (y)dy = λ s v s (x), which can be turned into an algebraic problem by discretizing the variables x and y with equal subintervals of length ∆y (see, for example, [14]). Thus, the eigenvalues of the matrix B = [∆yρ(m i , m j )] K×K , where m i = −c+ △yi, △y = 2c/(K − 1), i, j = 0, ..., K − 1, provide approximations to the K largest occupancies.…”