In this paper, we study the optimization algorithm to compute the smallest eigenvalue of the Schrödinger operator with volume constraint. A finite element discretization of this problem is established. We provide the error estimate for the numerical solution. The optimal solution can be approximated by a fixed point iteration scheme. Then a monotonic decreasing algorithm is presented to solve the eigenvalue optimization problem. Numerical simulations demonstrate the efficiency of the method.