Of the outstanding problems in astronomy, the nature of dark matter is certainly one of the most mysterious. Containing five times mover energy density than it luminous counterpart, dark matter has been shaping the large-scale structure of our Universe for billions of years. The expansion of accessible and accurate cosmological simulations has revolutionized how we visualize the imprint of dark matter in the structure of our Universe. In my Ph.D., I contributed to this revolution through the development and implementation of a new code, chplUltra: a parallel, portable, and efficient tool for HPC simulations of a promising dark matter candidate, Fuzzy or UltraLight Dark Matter (ULDM). ULDM is a well-motivated axion-like dark matter candidate whose incredibly small mass results in naturally cored profiles, thus ameliorating many of the small-scale problems of cold dark matter (CDM) while maintaining the same robust large-scale results. When unperturbed, the lowest energy solution of the ULDM system is a spherical “soliton” structure with a known mass density profile. A ULDM dark matter halo is formed through collisions of these solitons and has two characteristic parts: a central soliton core, and an “skirt” surrounding it. In order to investigate ULDM dynamics, I calculated the full spectrum of eigenstates for ULDM systems with approximately stationary potentials, thus allowing me to: 1) link qualitative behavior of soliton cores in ULDM simulations with superpositions of specific modes and 2) decompose chplUltra simulations of ULDM halos into individual eigenstates. Using this formalism, I investigated formation of halos through soliton collisions and the dependence of the final halo product on initial parameters. Crucially, this allowed me to explore how halo cores form and explain discrepancies in the literature surrounding the core-halo mass relation: a key prediction of ULDM. I was also able to comment on the composition of the halos’ skirts, including their qualitative behavior and eigenstate makeup, as a function initial binary parameters. Finally, I sketched out some of the exciting future directions for understanding ULDM through the language of its eigenstates.
Thesis Advisors: Nikhil Padmanabhan (nikhil.padmanabhan@yale.edu) and Richard Easther (University of Auckland)