Abstract
In order to study structure of the crust in neutron stars, we develop a finite-temperature Skyrme-Hartree-Fock method in the full three-dimensional coordinate space using the Fermion operator expansion method. It provides us with a possible order-N approach to non-uniform neutron-star matters.