Hi How can I add additional term to Hamiltonian such (r/r0)^N (in which r_0 is constant and N is integer that I'll set it later)? In fact the final potential should be like this: V(r)=V_nuc(r) + V_Hartree(r) + V_xc(r) + (r/r0)^N and then make the necessary file (the files need to be changed) and run scf cycle.
Thanks -- *Yours sincerely* *Sajad Hamreh* *PhD Student of Condensed matter physics* *Tarbiat Modares University* *www.comphys.ir <http://www.comphys.ir>*
