. The whole body of energy derivatives in the representation of atomic coordinates, both in the closed (N=const.) and open (µ=const.) systems has been formulated, their renormalization has also been demonstrated. The vital result of that effort focused directly on atoms identified by the H-F forces was the discovery of vibrational softening of molecules with perspectives for applications to explosive reactions. The linear response function has been essential for bridging the energy derivatives in the basic CDFT formalism and their counterparts built on