A linear‐scaling time‐dependent density functional theory is developed. The equation of motion is solved for the reduced single‐electron density matrix in the real time domain. Chebyshev expansion is used for integration in time domain. Filter diagonalization is implemented to determine the excited state energies. The locality of the reduced single‐electron density matrix is utilized to ensure computational time scales linearly with system size. We summarize these methods in this brief review.