Ferromagnetism in infinite neutron matter as well as beta equilibrated, charge neutral, dense, and infinite nuclear matter is investigated using a model of interacting baryons and mesons. The standard minimal couplings between the magnetic field and the particle charges as well as the baryon dipole moments are included in the Lagrangian density. Minimizing the energy density with respect to the magnetic field yields a self-consistent expression for the ferromagnetic field. We calculate the phase boundary at a given density by increasing the strength of the baryon dipole moments till the energy density of magnetized matter is lower than that of unmagnetized matter. We find that, depending on the density, it is crossed when the baryon dipole moments are increased by a factor of 35. It is also sensitive to the details of the nuclear matter parameterizations and crossing it induces a magnetic field of $\sim 10^{17}$ gauss.

The spin-polarized ferromagnetic state of a cold Fermi gas is investigated for interacting and non-interacting charge-neutral and $\beta$-equilibrated gases. The standard minimal couplings between the magnetic field and the fermions' charges and magnetic dipole moments define the fermions' interaction with the magnetic field. Assuming a variable coupling strength between the magnetic field and the fermion (baryon) dipole moments, it is shown that a ferromagnetized state can be achieved that corresponds to a lower energy spin-polarized state with a magnetic field entirely due to the gas's magnetic response. We find that, depending on the density, a very large increase in the baryon dipole moments is needed to achieve this ferromagnetized state. While the required increase seems unlikely, the induced magnetic field is of the order $\sim10^{17}$ gauss. Furthermore, while externally magnetized Fermi gases have an anisotropic pressure, the pressure of the ferromagnetized gas is completely isotropic and the thermodynamically preferred magnetized state.

Nuclear physics can be applied in various ways to the study of neutron stars. This thesis reports on one such application, where the relativistic mean-field approximation has been employed to calculate the equations of state of matter in the neutron star interior. In particular the equations of state of nuclear and neutron star matter of the NL3, PK1 and FSUGold parameter sets were derived. A survey of available literature on neutron stars is presented and we use the derived equations of state to reproduce the properties of saturated nuclear matter as well as the mass-radius relationship of a static, spherical symmetric neutron star. Results are compared to published values of the properties of saturated nuclear matter and to available observational data of the mass-radius relationship of neutron stars.