We consider a one-dimensional model where the confining potential is parametrized/obtained from finite element calculations. The bare Coulomb interaction is divergent in 1D (REF) and it is customary to use a smoothed Coulomb interaction given by
u(x1,x2)=α√(x1−x2)2+a2,where α adjusts the strength of the interaction and a removes the singularity at x1=x2.
The single-particle functions are chosen as the eigenfunctions of the single-particle Hamiltonian
(−d2dx2+v(x))ψp(x)=ϵpψp(x).