The method described in this paper opens a new way to exploration of the vital cDFT derivatives: the linear response function 𝜔(𝐫, 𝐫 ′ ) and the softness kernel 𝑠(𝐫, 𝐫 ′ ). Applications 𝜔(𝐫, 𝐫 ′ ) in chemistry have been subject of many efforts, hampered by the nonlocal character of this derivative. The Fukui function f(r) and the dual descriptor f (r), another energy derivative of the third order, have been found more promising for chemical use, due to their local character. This present work has demonstrated that the non-local response functions may be replaced (to some extent) by exploring the local derivatives of higher order. This resulted not only in calculations of the local descriptors, but also the method for their phenomenological use has been proposed: calculation of the disturbances of the state functions Δ𝐸 𝑁 , Δ𝜇 𝑁 , Δ𝑁 𝜇 . While the first two describe merely the polarization of a system due to the external potential, Δ𝑁 𝜇 describes the truly chemical effect: the charge flow. The nonlocal kernels 𝜔(𝐫, 𝐫 ′ ) and 𝑠(𝐫, 𝐫 ′ ) have not been neglected but replaced by the higher order local derivative: [𝜕𝑠(𝐫)/𝜕𝑁] 𝑣 . Calculation thereof is then necessary for the completeness of the method.