N-representability and stationarity in time-dependent density-functional theory

To construct an N-representable time-dependent density-functional theory, a generalization to the time domain of the Levy-Lieb (LL) constrained-search algorithm is required. That the action is only stationary in the Dirac-Frenkel variational principle eliminates the possibility of basing the search on the action itself. Instead, we use the norm of the partial functional derivative of the action in the Hilbert space of the wave functions in place of the energy of the LL search. The electron densities entering the formalism are N-representable, and the resulting universal action functional has a unique stationary point in the density at that corresponding to the solution of the Schrödinger equation. The original Runge-Gross (RG) formulation is subsumed within the current formalism. Concerns in the literature about the meaning of the functional derivatives and the internal consistency of the RG formulation are allayed by clarifying the nature of the functional derivatives entering the formalism.