Past experiments have revealed several unusual properties about interstitial hydrogen atoms in niobium. Absorption isotherms showed that niobium absorbs a large amount of hydrogen without changing its crystal structure. These isotherms also revealed that the interactions between hydrogen atoms in niobium are a combination of long-range attraction and short-range repulsion and exhibit many-body characteristics. Other experiments reported the facile thermal diffusion of hydrogen and deuterium in niobium. Contrary to the classical theory of diffusion, these experiments revealed a break in the activation energy of hydrogen diffusion at low temperatures, but no such break was reported for deuterium. Finally, experiments report a phenomenon called electromigration, where hydrogen atoms inside niobium respond to weak electric fields as if they had a positive effective charge. These experimental results date back to when tools like density functional theory (DFT) and modern high-performance computing abilities did not exist. Therefore, the current understanding of these properties is primarily based on inferences from experimental results. Understanding these properties at a deeper level, besides being scientifically important, can profoundly affect various applications involving hydrogen separation and transport. The high-level goal of this work is to use first-principles methods to explain the discussed properties of interstitial hydrogen in niobium. DFT calculations were used to study hydrogen atoms' site preference in niobium and its effect on the cell shape and volume of the host cell. The nature and origin of the interactions between hydrogen atoms were studied through interaction energy, structural, partial charge, and electronic densities of state analysis. A phenomenological model with fewer parameters than traditional models was developed and fit to the experimental absorption data. Thermodynamic quantities such as the enthalpy and entropy of hydrogen dissolution in niobium were derived from this model. The enthalpy of hydrogen dissolution in niobium was also calculated using DFT by sampling different geometric configurations and performing an ensemble-based averaging. Further work is required to explain the observed isotope effects for hydrogen diffusion in niobium and the electromigration phenomena. Applications of the niobium-hydrogen system require studying hydrogen's behavior on niobium's surface.