Recent development in particle-particle random phase approximation (pp-RPA) broadens the perspective on ground state correlation energies (van Aggelen et al. Phys. Rev. A, 88, 030501 (2013)) and N ± 2 excitation energies (Yang, et al. J.Chem. Phys.). So far Hartree-Fock and approximated density-functional orbitals have been utilized to evaluate the pp-RPA equation. In this paper, to further explore the fundamentals and the potential use of pairing matrix dependent functionals, we present the linear-response time-dependent density-functional theory with pairing elds with both adiabatic and frequency-dependent kernels. This theory is related to the density-functional theory and time-dependent density-functional theory for superconductors, but is applied to normal non-superconducting systems for our purpose. Due to the lack of the proof of the one-to-one mapping between the pairing matrix and the pairing eld for time-dependent systems, the linear-response theory is established based on the representability assumption of the pairing matrix. The linear response theory justi es the use of approximated density-functionals in the pp-RPA equation. This work sets the fundamentals for future density-functional development to enhance the description of ground state correlation energies and N ± 2 excitation energies.