We present explicit formulas for arbitrary-order derivatives of the energy, grand potential, electron density, and higher-order response functions with respect to the number of electrons, and the chemical potential for any smooth and differentiable model of the energy versus the number of electrons. The resulting expressions for global reactivity descriptors (hyperhardnesses and hypersoftnesses), local reactivity descriptors (hyperFukui functions and local hypersoftnesses), and nonlocal response functions are easy to evaluate computationally. Specifically, the explicit formulas for global/local/nonlocal hypersoftnesses of arbitrary order are derived using Bell polynomials. Explicit expressions for global and local hypersoftness indicators up to fifth order are presented. (C) 2016 Elsevier B.V. All rights reserved.