In this work, a time fractional Schrödinger equation with Caputo fractional derivative of order \(\alpha \in (0,1)\) is considered, whose solution exhibits a weak singularity at initial time. We divide the solution into two parts: imaginary part and real part, then the discrete system is constructed by the five-point difference scheme in space together with the L1 approximation of Caputo derivative on graded mesh in time. We provide stability and local error analysis for the discrete scheme, which shows that one can attain the optimal convergence order \(2-\alpha \) in positive time by selecting milder grading parameters. Numerical results are presented to verify the accuracy of the theoretical analysis.
