In this paper, an adaptive numerical method is proposed for solving a 2D Schrödinger equation with an imaginary time propagation approach. The differential equation is first transferred via a Wick rotation to a real time-dependent equation, whose solution corresponds to the ground state of a given system when time approaches infinity. The temporal equation is then discretized spatially via a finite element method, and temporally utilizing a Crank–Nicolson scheme. A moving mesh strategy based on harmonic maps is considered to eliminate possible singular behaviour of the solution. Several linear and nonlinear examples are tested by using our method. The experiments demonstrate clearly that our method provides an effective way to locate the ground state of the equations through underlying eigenvalue problems.