This thesis contains four articles related to mathematical aspects of Density Functional Theory. In Paper A, the theoretical justification of density methods formulated with current densities is addressed. It is shown that the set of ground-states is determined by the ensemble-representable particle and paramagnetic current density. Furthermore, it is demonstrated that the Schrödinger equation with a magnetic field is not uniquely determined by its ground-state solution. Thus, a wavefunction may be the ground-state of two different Hamiltonians, where the Hamiltonians differ by more than a gauge transformation. This implies that the particle and paramagnetic current density do not determine the potentials of the system and, consequently, no Hohenberg-Kohn theorem exists for Current Density Functional Theory formulated with the paramagnetic current density. On the other hand, by instead using the particle density as data, we show that the scalar potential in the system's Hamiltonian is determined for a fixed magnetic field. This means that the Hohenberg-Kohn theorem continues to hold in the presence of a magnetic field, if the magnetic field has been fixed. Paper B deals with N-representable density functionals that also depend on the paramagnetic current density. Here the Levy-Lieb density functional is generalized to include the paramagnetic current density. It is shown that a wavefunction exists that minimizes the "free" Hamiltonian subject to the constraints that the particle and paramagnetic current density are held fixed. Furthermore, a convex and universal current density functional is introduced and shown to equal the convex envelope of the generalized Levy-Lieb density functional. Since this functional is convex, the problem of finding the particle and paramagnetic current density that minimize the energy is related to a set of Euler-Lagrange equations. In Paper C, an N-representable Kohn-Sham approach is developed that also include the paramagnetic current density. It is demonstrated that a wavefunction exists that minimizes the kinetic energy subject to the constraint that only determinant wavefunctions are considered, as well as that the particle and paramagnetic current density are held fixed. Using this result, it is then shown that the ground-state energy can be obtained by minimizing an energy functional over all determinant wavefunctions that have finite kinetic energy. Moreover, the minimum is achieved and this determinant wavefunction gives the ground-state particle and paramagnetic current density. Lastly, Paper D addresses the issue of a Hohenberg-Kohn variational principle for Current Density Functional Theory formulated with the total current density. Under the assumption that a Hohenberg-Kohn theorem exists formulated with the total current density, it is shown that the map from particle and total current density to the vector potential enters explicitly in the energy functional to be minimized. Thus, no variational principle as that of Hohenberg and Kohn exists for density methods formulated with the total current density. / <p>QC 20140523</p>