In this article, density functionals for Coulomb systems subjected to electric and magnetic fields are developed. The density functionals depend on the particle density ρ and paramagnetic current density jp. This approach is motivated by an adapted version of the Vignale and Rasolt formulation of current density functional theory, which establishes a one‐to‐one correspondence between the nondegenerate ground‐state and the particle and paramagnetic current density. Definition of N‐representable density pairs (ρ,jp) is given and it is proven that the set of v‐representable densities constitutes a proper subset of the set of N‐representable densities. For a Levy–Lieb‐type functional Q(ρ,jp), it is demonstrated that (i) it is a proper extension of the universal Hohenberg–Kohn functional
FHK(ρ,jp) to N‐representable densities, (ii) there exists a wavefunction ψ0 such that
Q(ρ,jp)=(ψ0,H0ψ0)L2, where H0 is the Hamiltonian without external potential terms, and (iii) it is not convex. Furthermore, a convex and universal functional F(ρ,jp) is studied and proven to be equal the convex envelope of Q(ρ,jp). For both Q and F, we give upper and lower bounds. © 2014 Wiley Periodicals, Inc.