I want to try to model the electornic structure of atoms in a confining potential, eg. add a term to the Hamiltonian either of a constant value $V$ after certain distance from the origin or just the form $r^n$. I've came across this example which adds an electric field to the molecule, and it is ...