Recently, I started to learning solving schrodinger equation numerically and I noticed that a lot of researchers implemented the [Gausssian expansion method](Gaussian expansion method for few-body systems - ScienceDirect) developed by Hiyama. The approach works efficiently in most cases but one problem occurs to me: how can I implement the boundary condition after solving the generalized eigen matrix problem(eq.11)? We know that for centre potential the reduced wavefunction(R(r)=u(r)/r) vanishes at the original point but the eigen vector from the generalized eigen matrix problem doesn't lead to the above result. Therefore, what should I do to apply the boundary condition to the eigen problem? Any help would be appreciated!!!
Just now I realized that the gaussian expansion method is applied to the 3-d schrodinger equation and the wavefunction at the original point is finite but zero(in the radial form). So the boundary condition doesn't apply to the generalized eigen matrix problem. The question is closed!
- Badges
- Science topic