My work primarily focuses on the transformation Density Functional Theory (DFT) into a 1-Electron Reduced Density Matrix Functional Theory (1-RDMFT), which respectively utilize the electronic density and 1-Electron Reduced Density Matrix (1-RDM) instead of an exponentially complex wavefunction, for the calculation of electronic properties in strongly correlated molecules. The combination of these two approaches allows us to retain the favorable computational scaling of DFT with respect to system size while still being able to capture strong correlation effects through the 1-RDM. With this approach are able to capture multi-reference effects such as the Mott insulator transition of a linear chain of hydrogens as they are stretched apart or the orbital degeneracies of dissociated diatomics.