A family of new iteration methods is presented for designing quantum optimal controls to manipulate the transition probability. Theoretical analysis shows that these new methods exhibit quadratic and monotonic convergence. Numerical calculations verify that for these new methods, within very few steps, the optimized objective functional comes close to its convergent limit.