\lambda_D=\left(\lambda_{D e}^{-2}+\lambda_{D i}^{-2}\right)^{-1 / 2}, \lambda_{D j}=\left(\varepsilon_0 k T_j / e^2 n_{j o}\right)^{1 / 2}, \quadλD=(λDe−2+λDi−2)−1/2,λDj=(ε0kTj/e2njo)1/2, and ωpd=(e2Zd2ndo/ε0md)1/2. For the case of cold dust (Td→0) and in the long wavelength limit (k2λD2≪1), the dispersion relation reduces to ω/k=cDA, where cDA is the dust acoustic speed given by cDA2=(mdkTiniondoZd2)