Quasi-particle spectrum of nano-scale conventional and unconventional superconductors under magnetic field

We have developed a numerical method to solve the Bogoliubov-de Gennes equation for nano-scaled isotropic and d-wave superconductors. It is based on the finite element method, and therefore applicable to arbitrary geometries. We argue the difference of the local density of states between isotropic and a d-wave superconducting square plate. For d-wave case, it appears as intrinsic surface bound states.