Computation of derivatives (gradient and Hessian) of a fidelity function is one of the most crucial steps in many optimization algorithms. Having access to accurate methods to calculate these derivatives is even more desired where the optimization process requires propagation of these calculations over many steps, which is in particular important in optimal control of spin systems. Here we propose a novel numerical approach, ESCALADE (Efficient Spin Control using Analytical Lie Algebraic Derivatives) that offers the exact first and second derivatives of the fidelity function by taking advantage of the properties of the Lie group of 2 × 2 Hermitian matrices, SU(2), and its Lie algebra, the Lie algebra of skew-Hermitian matrices, su(2). A full mathematical treatment of the proposed method along with some numerical examples are presented.