Density-matrix based formulation of KS response theory with applications to linear and non-linear spectroscopies