Most stars in our galaxy are M–dwarfs, much cooler and smaller than the sun. The ubiquitous nature of these stars is also paired with the formation of terrestrial exoplanets orbiting them. The strategic placement of M-dwarfs between main-sequence stars and brown dwarfs, their uniqueness as exoplanet analogs, and their dominating presence in the galactic stellar population make them priority targets for study. This work investigates outstanding questions, including the need to acquire constraints on their chemical compositions to decode formation processes, evolution, and interaction with companion objects. Chapter 1 lays out a broad background emphasizing the importance of studying the most populous star in the galaxy, their far-reaching implications, and primarily the numerous challenges in characterizing the atmospheres and environments of these stars. Chapter 2 investigates the influence of M-dwarf star spots propagating into spectra of transiting terrestrial planets, showing that inaccurate modeling of M-dwarf photospheres leads to significant bias when inferring atmospheric properties of companion exoplanets. These biases persist despite correcting M-dwarf spot signatures imprinted onto the exoplanetary spectra, even with high-fidelity JWST observations. This result emphasizes the need for improved stellar atmosphere models as the first step to improving our understanding of the companion planets. To address this, chapter 3 introduces SPHINX—a new stellar atmosphere model grid for M-dwarfs. SPHINX provides improved constraints on fundamental properties of benchmark M-dwarf systems (e.g., temperature, surface gravity, radius, and chemistry). The improvement is significant relative to the state-of-the-art stellar model grid available today. Chapter 4 expands this model, applying it to mid-to-late type M-dwarfs, and investigating chemical trends in their atmospheric properties. Using low-resolution observations, both archival data (from SpeX Prism Library Database) and from previous empirical studies; this chapter presents constraints on fundamental atmospheric properties of 71 low-mass, late-type M-dwarfs to understand spectroscopic degeneracies arising due to stellar activity, cloud/dust condensation and convection. With SPHINX models, the chemical properties of these stars are compared against main-sequence stars to acquire a more holistic understanding of M-dwarfs as a class—in the quest to ultimately characterize their companions.

The isotopic compositions of meteorites provide valuable insights into the earliest history of the Solar System and, in some cases, provide constraints on presolar components that contributed to the solar nebula. In the past decade or so, mass-independent isotope anomalies in titanium have become particularly important geochemical tracers to study the distinct isotopic reservoirs in the early Solar System. In particular, mass-independent anomalies in the most neutron-rich isotope of titanium (50Ti) have been used to distinguish between carbonaceous chondritic (CC) and non-carbonaceous chondritic (NC) materials. These two groupings likely represent distinct isotopic reservoirs in the inner (NC) and outer (CC) Solar System. However, while the titanium isotope compositions of CC and NC materials are distinct, each group's full range of compositional variability is poorly characterized. For example, only one CK carbonaceous chondrite group member has been analyzed thus far for its bulk Ti isotope composition. This work aims to characterize better the range of mass-independent Ti isotope compositions within and among the carbonaceous chondrites, which has implications for the degree and potential sources of Ti isotope heterogeneity in the early Solar System. Methods utilized in this study include column chromatography to purify Ti and high-precision multi-collector inductively coupled plasma mass spectrometry for measuring Ti isotope compositions. The Ti isotope compositions of bulk samples of nine carbonaceous chondrites are reported here. In addition, the bulk fractions of the meteorites used in this study were taken from homogenized powders of relatively large (~200 mg each) samples. This was done to assess whether variability in mass-independent Ti isotope compositions previously reported within some meteorites could be a sampling artifact. Results from this work show that the various CM2 chondrites and ungrouped carbonaceous chondrites have ε50Ti values that are similar, suggesting that the Ti in these samples was likely sourced from a common isotopic reservoir. On the other hand, the ε50Ti values reported for CI1 and CH/CBb bulk samples suggest that the parent bodies of these carbonaceous chondrite groups were formed in isotopic reservoir(s) distinct from that of the other CC groups in the early Solar System.

Direct imaging is a powerful tool in revealing the architectures of young planetary systems, clearly showing the structure of circumstellar disks. Circumstellar disks, similar to the asteroid belt, are critical elements of any planetary system, and the study of them is important to understanding planet formation. Disks around several main sequence stars have already been observed directly interacting with exoplanets in their respective systems. Imaging can help answer many of the key questions of how disks interact in their respective systems. The Gemini Planet Imager is a high contrast imaging instrument that has spatially resolved several circumstellar disks for the first time, many exhibiting tracers of ongoing planet formation or the presence of a perturbing exoplanet. With this new sample, population analyses of characteristics of disks can now be explored and compared to information at other wavelengths. Direct imaging is also a uniquely accessible tool in engaging students and the community in astronomy. In combination with a course-based undergraduate research experience, direct imaging has the ability to engage students in the process of doing research in a very accessible manner. In Chapter 1, I introduce the concepts related to circumstellar debris disks, further focusing on the sub-field of direct imaging and its value in understanding these systems and engaging students in astronomy. In Chapter 2, I present four images of newly-resolved debris disks in the Scorpius-Centaurus association, comparing their characteristics with many other spatially-resolved circumstellar disks within the moving group. In Chapter 3, I present a uniform analysis of debris disk structure using a consistent and empirically-informed modeling approach. In Chapter 4, I present my findings and experiences in developing and teaching a course-based undergraduate research experience for students in the country’s first online astronomy degree program centered on the direct imaging of brown dwarfs. In Chapter 5, I present my conclusions on the topics I have investigated and discuss future work within the field of direct imaging and its role in driving astronomy research and education forward.

White Dwarf stars are the stellar remnants of low mass stars which have completed their evolution. Nearly all stars will become white dwarfs. The interior of a white dwarf encapsulates its evolution history: unraveling a white dwarf’s internal structure constrains the physical events which occurred to construct its composition. Variable, or pulsating, white dwarfs emit pulsations which are sensitive to their internal stratification. Just as seismology reveals Earth’s interior, asteroseismology can reveal stellar interiors. The standard approach to construe an observed white dwarf’s chemical makeup is to match observed pulsation properties to theoretical stellar models. Observed white dwarf pulsation data has reached 6-7 significant digits of precision. As such, it is important for computational modeling to consider systematic offsets from initial conditions and theoretical uncertainties that are within the detectable threshold. By analyzing the magnitude of pulsation differences among various uncertainties from white dwarf models, one can place constraints on important theoretical uncertainties. In this thesis, I explore impacts on white dwarf pulsations that result from accounting for various uncertainties in computational models. I start by showing the importance of 22Ne, and its impact on the pulsations in Helium atmosphere white dwarfs. Next, I discuss how certain trapped modes of white dwarfs may yield a signal for the 12C(α,γ)16O reaction rate probability distribution function. This reaction occurs during the Helium core burning phase in stellar evolution, and chiefly determines the Carbon and Oxygen abundance of white dwarfs. Following this work, I show how overshooting impacts the pulsation signatures of the 12C(α, γ)16O reaction rate. I then touch on the analytical work I’ve done regarding educational research in the HabWorlds course offered at Arizona State University (ASU). I then summarize my conclusions from these efforts.

The balance between relative numbers, lifetime, and habitable zone (HZ) size of K stars (0.6 – 0.9 M⊙) in comparison with M (0.08 – 0.6 M⊙) and G (0.9 – 1.1 M⊙) stars makes them candidates to host “super-habitable” planets. Understanding the high- energy radiation environment of planets around these stars is crucial, since ultraviolet (UV) and X-ray radiation may cause severe photodissociation and ionization of the atmosphere, with the potential for complete erosion. In this thesis, I present the first broad study of the UV and X-ray evolution of K stars. I first focused on Galaxy Evolution Explorer (GALEX) and Ro ̈ntgen Satellit (ROSAT) photometric UV and X-ray evolutions of K stars and compared this with the age evolution of both early- (0.35 – 0.6 M⊙) and late-M (0.08 – 0.35 M⊙) stars. I found that the fractional UV and X-ray flux from M and K stars is similar; however, the wider and farther HZs of K stars mean that there is less incident UV radiation on HZ planets. Next, I led a spectroscopic study of 41 K stars using Hubble Space Telescope Cosmic Origins Spectrograph (HST/COS) data to show that the UV line and continua emission show no decrease in flux beyond 650 Myr whereas early-M star flux declines by 150 Myr; therefore, the K star intrinsic UV flux is greater than early-M stars after this time. I suggest that this phenomenon is related to K star rotational spin-down stalling. Lastly, I revisited the GALEX and ROSAT data with newly-available distances from the Gaia mission for both K and M stars. I find that the UV flux for K stars is an order of magnitude higher for M stars at all ages and the flux in their respective HZs is similar. However, K star X-ray flux is an order of magnitude less in the HZ than for M stars. The age of decline shows a dependency on wavelength, a phenomenon which is not seen in either the early- or late-M star data. These results suggest thatK stars may not exhibit quite the advantage as HZ planet host stars as the scientific community originally thought.

Brown dwarfs are a unique class of object which span the range between the lowest mass stars, and highest mass planets. New insights into the physics and chemistry of brown dwarfs comes from the comparison between spectroscopic observations, and theoretical atmospheric models. In this thesis, I present a uniform atmospheric retrieval analysis of the coolest Y, and late-T spectral type brown dwarfs using the CaltecH Inverse ModEling and Retrieval Algorithms (CHIMERA). In doing so, I develop a foundational dataset of retrieved atmospheric parameters including: molecular abundances, thermal structures, evolutionary parameters, and cloud properties for 61 different brown dwarfs. Comparisons to other modeling techniques and theoretical expectations from the James Webb Space Telescope (JWST) are made. Finally, I describe the techniques used to improve CHIMERA to run on Graphical Processing Units (GPUs), which directly enabled the creation of this large dataset.

The interactions that take place in the ionized halo of gas surrounding galaxies, known as the circumgalactic medium (CGM), dictates the host galaxy's evolution throughout cosmic time. These interactions are powered by inflows and outflows that enable the transfer of matter and energy, and are driven by feedback processes such as accretion, galactic winds, star formation and active galactic nuclei. Such feedback and the interactions that ensue leads to the formation of non-equilibrium chemistry in the CGM. This non-equilibrium chemistry is implied by observations that reveal the highly non-uniform distribution of lower ionization state species, such as Mg II and Si II, along with widespread higher ionization state material, such as O VI, that is difficult to match with equilibrium models. Given these observations, the CGM must be viewed as a dynamic, multiphase medium, such as occurs in the presence of turbulence. To better understand this ionized halo, I used the non-equilibrium chemistry package, MAIHEM, to perform hydrodynamic (HD) simulations. I carried out a suite of HD simulations with varying levels of artificially driven, homogeneous turbulence to learn how this influences the non-equilibrium chemistry that develops under certain conditions present in the CGM. I found that a level of turbulence consistent with velocities implied by observations replicated many observed features within the CGM, such as low and high ionization state material existing simultaneously. At higher levels of turbulence, however, simulations lead to a thermal runaway effect. To address this issue, and conduct more realistic simulations of this environment, I modeled a stratified medium in a Milky Way mass Navarro-Frenk-White (NFW) gravitational potential with turbulence that decreased radially. In this setup and with similar levels of turbulence, I alleviated the amount of thermal runaway that occurs, while also matching observed ionization states. I then performed magneto-hydrodynamic (MHD) simulations with the same model setup that additionally included rotation in the inner halo. Magnetic fields facilitate the development of an overall hotter CGM that forms dense structures within where magnetic pressure dominates. Ion ratios in these regions resemble detections and limits gathered from recent observations. Furthermore, magnetic fields allow for the diffusion of angular momentum throughout the extended disk and gas cooling onto the disk, allowing for the maintenance of the disk at late times.

The lives of high-mass stars end with core-collapse supernovae, which distribute energy and chemical elements into the interstellar medium. This process is integral to the Galactic ecosystem, since stars and planets will form from the enriched interstellar medium. Since most supernovae are detected at intergalactic distances, opportunities to examine them in detail are rare. Computer simulations and observations of supernova remnants are frequently employed to study these events and their influence on the universe. I explore the topic of supernovae using a multi-pronged approach, beginning with an examination of the core-collapse supernova engine. The radioisotopes 44Ti and 56Ni, produced in the innermost ejecta, provide a probe of this central engine. Using a three-dimensional supernova simulation with nucleosynthesis post-processing, I examine the production of these isotopes and their thermodynamic histories. Since production of 44Ti is especially sensitive to the explosion conditions, insights can be gained by comparing the model with 44Ti observations from supernova remnant Cassiopeia A. Next, I consider supernova remnants as potential sources of high-energy neutrinos within the Milky Way galaxy. The developing field of neutrino astronomy has yet to identify the origins of the diffuse neutrino flux first detected by the IceCube Neutrino Observatory in 2013. In principle, high-energy Galactic sources like supernova remnants could contribute measurably to this flux. I also consider Galactic open clusters, environments which are rich in supernovae and other energetic phenomena. Statistical analysis finds no evidence of causal association between these objects and the IceCube neutrino events. I conclude with a series of asymmetric three-dimensional supernova models, presented as a comparative analysis of how supernova morphology affects nucleosynthetic yields. Both real supernovae and simulations frequently exhibit aspherical morphologies, but the detailed thermodynamic consequences and the ultimate effects on yields are poorly understood. The simulations include symmetric and bipolar explosion geometries for both 15- and 20-solar-mass progenitor stars. Across the spectrum of models, I show how small changes in the peak temperatures and densities experienced by ejecta can influence the production of notable isotopes such as 44Ti.

Study of the early Universe is filled with many unknowns, one of which is the nature of the very first generation of stars, otherwise designated as "Population III stars". The early Universe was composed almost entirely of cold hydrogen and helium, with only trace amounts of any heavier elements. As such, these stars would have compositions very different from the stars we are able to observe today, which would in turn change how these stars functioned, as well as their lifespans. Population III stars are so old that the light they emitted has not yet reached us here on Earth. Yet we know they have to have existed, so how do we go about studying objects that we have not yet observed? And more importantly, is there a metallicity threshold at which stars begin to behave like the stars we observe today? These areas are where stellar modelling programs such as TYCHO8 and the Spanish Virtual Observatory's Theoretical Spectra Web Server (TSWS) come in. These programs allow astronomers to model the physics of Pop III stars. We can get a pretty good understanding of how these stars behaved, how long they lived, and the visual spectra they would have emitted. Such information is crucial to astronomers being able to search for remnants of these stars, and one day, the stars themselves.