  Linear-response time-dependent density functional theory (LR-TDDFT)
simulations of disordered extended systems require averaging over different
snapshots of ion configurations to minimize finite size effects due to the
snapshot--dependence of the electronic density response function and related
properties. We present a consistent scheme for the computation of the
macroscopic Kohn-Sham (KS) density response function connecting an average over
snapshot values of charge density perturbations to the averaged values of KS
potential variations. This allows us to formulate the LR-TDDFT within the
adiabatic (static) approximation for the exchange-correlation (XC) kernel for
disordered systems, where the static XC kernel is computed using the direct
perturbation method [Moldabekov et al., J. Chem. Theory Comput. 19, 1286
(2023)]. The presented approach allows one to compute the macroscopic dynamic
density response function as well as the dielectric function with a static XC
kernel generated for any available XC functional. The application of the
developed workflow is demonstrated for the example of warm dense hydrogen. The
presented approach is applicable for various types of extended disordered
systems such as warm dense matter, liquid metals, and dense plasmas.
