The interaction between galaxies and the surrounding gas plays a key role in galaxy formation and evolution. Feedback processes driven by star formation and active galactic nuclei facilitate the exchange of mass and energy between the galaxy and the circumgalactic medium through inflowing and outflowing gas. These outflows have a significant impact on the star formation rate and metallicity of the galaxy. Observations of outflows have provided evidence that these outflows are multi-phase in nature, identifying both low energy ions such as Mg II and C III and high energy ions such as O VI. The underlying physics maintaining the two phases as well as the ionization mechanism for these phases remains unclear. In order to better understand galactic outflows, hydrodynamic simulations are used to study the evolution of wind-cloud interactions. In this work, I carried out a suite of magnetohydrodynamic simulations to characterize the influence of magnetic fields on the evolution and lifetime of cold clouds. I found magnetic fields either provided little improvement to cloud stability over other influences such as radiative cooling or accelerated cloud disruption by pushing cloud material in the direction orthogonal to the wind and magnetic fields. To investigate the ionization mechanism of the material within outflows I first considered estimating the column densities of various ions within wind-cloud simulations with the post-processing tool Trident. Under the assumption of ionization equilibrium, the simulations did not reproduce the observed absorption profiles demonstrating the need for a more detailed treatment of the ionization processes. I then performed a new set of simulations with the non-equilibrium chemistry solver, MAIHEM. The column densities produced in the non-equilibrium model alter the evolution of the cloud and highlight the increased ionization along the boundary of the cloud.

Stellar mass loss has a high impact on the overall evolution of a star. The amount<br/>of mass lost during a star’s lifetime dictates which remnant will be left behind and how<br/>the circumstellar environment will be affected. Several rates of mass loss have been<br/>proposed for use in stellar evolution codes, yielding discrepant results from codes using<br/>different rates. In this paper, I compare the effect of varying the mass loss rate in the<br/>stellar evolution code TYCHO on the initial-final mass relation. I computed four sets of<br/>models with varying mass loss rates and metallicities. Due to a large number of models<br/>reaching the luminous blue variable stage, only the two lower metallicity groups were<br/>considered. Their mass loss was analyzed using Python. Luminosity, temperature, and<br/>radius were also compared. The initial-final mass relation plots showed that in the 1/10<br/>solar metallicity case, reducing the mass loss rate tended to increase the dependence of final mass on initial mass. The limited nature of these results implies a need for further study into the effects of using different mass loss rates in the code TYCHO.

The spectra of brown dwarfs are key to exploring the chemistry and physics thattake place in their atmospheres. Late T dwarf (950 - 500 K) spectra are particularly
diagnostic due to their relatively cloud free atmospheres and deep molecular
bands. With the use of powerful atmospheric retrieval tools, these properties permit
constraints on molecular/atomic abundances and temperature profiles. Building
upon previous analyses on T and Y dwarfs (Line et al. 2017; Zalesky et al. 2019),
I present a uniform retrieval analysis of 50 T dwarfs via their low-resolution near infrared
spectra. This analysis more than doubles the sample of T dwarfs with retrieved
properties. I present updates on current compositional trends and thermal
profile constraints amongst the T dwarf population. My analysis shows that my collection
of objects form trends that are consistent with solar grid model expectations
for water, ammonia, methane, and potassium. I also establish a consistency between
the thermal structures of my objects with those of grid models. Moreover, I explore
the origin of gravity-metallicity discrepancies that are observed in some of my brown
dwarf candidates.

I examine the effects of metallicity on solar mass stellar evolution, trying to replicate a previous result in Windhorst et.al., 2018, in which a zer metallicity solar mass star did not reach the AGB, and thus may turn into a helium white dwarf. In trying to replicate this result, I used the M.E.S.A. stellar evolution code and was unable to reproduce this result. While M.E.S.A has undergone several updates since the previous result was obtained, more current evidence suggests that this may have been a one-time occurrence, as no helium white dwarfs were produced for low-metallicity models. Nonetheless, interesting results were obtained, including a lowest metallicity value for which CNO burning does not significantly contribute during the main sequence, 1 −10 Z , which produces noticeable effects on post main sequence evolution. All models are run with no rotation, one solar mass, and a series of MESA parameters kept constant, with the only exception being metallicity. Any metallicity value listed as Nd −10 is an absolute mass fraction, and Z is relative to solar metallicity, 2d*10 −2 .

Part I – I analyze a database of Smoothed Particle Hydrodynamics (SPH) simulations of collisions between planetary bodies and use the data to define semi-empirical models that reproduce remant masses. These models may be leveraged when detailed, time-dependent aspects of the collision are not paramount, but analytical intuition or a rapid solution is required, e.g. in ‘N-body simulations’. I find that the stratification of the planet is a non-negligible control on accretion efficiency. I also show that the absolute scale (total mass) of the collision may affect the accretion efficiency, with larger bodies more efficiently disrupting, as a function of gravitational binding energy. This is potentially due to impact velocities above the sound speed. The interplay of these dependencies implies that planet formation, depending on the dynamical environment, may be separated into stages marked by differentiation and the growth of planets more massive than the Moon.

Part II – I examine time-resolved neutron data from the Dynamic Albedo of Neutrons (DAN) instrument on the Mars Science Laboratory (MSL) Curiosity rover. I personally and independently developed a data analysis routine (described in the supplementary material in Chapter 2) that utilizes spectra from Monte Carlo N-Particle Transport models of the experiment and the Markov-chain Monte Carlo method to estimate bulk soil/rock properties. The method also identifies cross-correlation and degeneracies. I use data from two measurement campaigns that I targeted during remote operations at ASU. I find that alteration zones of a sandstone unit in Gale crater are markedly elevated in H content from the parent rock, consistent with the presence of amorphous silica. I posit that these deposits were formed by the most recent aqueous alteration events in the crater, since subsequent events would have produced matured forms of silica that were not observed. I also find that active dunes in Gale crater contain minimal water and I developed a Monte Carlo phase analysis routine to understand the amorphous materials in the dunes.

Learning how properties of galaxies such as star formation, galaxy interactions, chemical composition, and others evolve to produce the modern universe has long been a goal of extragalactic astronomy. In recent years, grism spectroscopy from the Hubble Space Telescope (HST) has provided a means to study these properties with spectroscopy while avoiding the limitations of ground-based observation. In this dissertation, I present several studies wherein I used HST G102 grism spectroscopy from the Faint Infrared Grism Survey (FIGS) to investigate these fundamental properties of galaxies and how they interact and evolve. In the first study, I combined the grism spectra with broadband photometry to produce a catalog of redshifts with improved accuracy, reducing the median redshift error from 3\% to 2\%. With this redshift catalog, I conducted a systematic search for galaxy overdensities in the FIGS fields, producing a list of 24 significant candidates. In the second study, I developed a method for identifying emission line galaxy (ELG) candidates from continuum-subtracted 1D spectra, and identified 71 ELGs in one FIGS field. In matching MUSE/VLT spectra, I measured the [OIII]$\lambda$4363 emission line for 14 FIGS ELGs, and used this to measure their $T_e$-based gas-phase metallicities. These ELGs show a low-metallicity offset on the Mass-Metallicity Relation, and I demonstrated that this offset can be explained by recent star formation. In the third study, I expanded the ELG search to all four FIGS fields, identifying 208 H$\alpha$, [OIII]$\lambda\lambda$4959,5007, and [OII]$\lambda\lambda$3727,3729 line emitters. I compiled a catalog of line fluxes, redshifts, and equivalent widths. I combined this catalog with the overdensity study to investigate a possible relationship between line luminosity, star formation, and an ELG's environment. In the fourth study, I usde 15 FIGS H$\alpha$ emitters and 49 ``green pea'' line emitters to compare H$\alpha$ and the far-UV continuum as tracers of star formation. I explored a correlation between the H$\alpha$-FUV ratio and the ratio of [OIII]$\lambda\lambda$4959,5007 to [OII]$\lambda\lambda$3727,3729 and its implications for star formation history.

The intention of this report is to use computer simulations to investigate the viability of two materials, water and polyethylene, as shielding against space radiation. First, this thesis discusses some of the challenges facing future and current manned space missions as a result of galactic cosmic radiation, or GCR. The project then uses MULASSIS, a Geant4 based radiation simulation tool, to analyze the effectiveness of water and polyethylene based radiation shields against proton radiation with an initial energy of 1 GeV. This specific spectrum of radiation is selected because it a component of GCR that has been shown by previous literature to pose a significant threat to humans on board spacecraft. The analysis of each material indicated that both would have to be several meters thick to adequately protect crew against the simulated radiation over a several year mission. Additionally, an analysis of the mass of a simple spacecraft model with different shield thicknesses showed that the mass would increase significantly with internal space. Thus, using either material as a shield would be expensive as a result of the cost of lifting a large amount of mass into orbit.

The composition of planets and their volatile contents are intimately connected to the structure and evolution of their parent protoplanetary disks. The transport of momentum and volatiles is often parameterized by a turbulent viscosity parameter $\alpha$, which is usually assumed to be spatially and temporally uniform across the disk. I show that variable $\alpha$(r,z) (where $r$ is radius, and $z$ is height from the midplane) attributable to angular momentum transport due to MRI can yield disks with significantly different structure, as mass piles up in the 1-10 AU region resulting in steep slopes of p $>$ 2 here (where p is the power law exponent in $\Sigma \propto r^{-p}$). I also show that the transition radius (where bulk mass flow switches from inward to outward) can move as close in as 3 AU; this effect (especially prominent in externally photoevaporated disks) may significantly influence the radial water content available during planet formation.

I then investigate the transport of water in disks with different variable α profiles. While radial temperature profile sets the location of the water snowline (i.e., inside of which water is present as vapor; outside of which, as ice on solids), it is the rates of diffusion and drift of small icy solids and diffusion of vapor across the snow line that determine the radial water distribution. All of these processes are highly sensitive to local $\alpha$. I calculate the effect of radially varying α on water transport, by tracking the abundance of vapor in the inner disk, and fraction of ice in particles and larger asteroids beyond the snow line. I find one α profile attributable to winds and hydrodynamical instabilities, and motivated by meteoritic constraints, to show considerable agreement with inferred water contents observed in solar system asteroids.

Finally, I calculate the timing of gap formation due to the formation of a planet in disks around different stars. Here, I assume that pebble accretion is the dominant mechanism for planetary growth and that the core of the first protoplanet forms at the water snow line. I discuss the dependence of gap timing to various stellar and disk properties.

The formation of the firsts stars some 100-300 Myr after the Big Bang marked the end of the cosmic darks ages and created the elemental building blocks of not only rocky planets but eventually us. Understanding their formation, lifetimes, and contributions to the evolution of our universe is one of the current frontiers in astronomy and astrophysics.

In this work I present an improved model for following the formation of Pop III stars, their effects on early galaxy evolution, and how we might search for them. I make use of a new subgrid model of turbulent mixing to accurately follow the time scales required to mix supernova (SN) ejecta -- enriched with heavy elements -- into the pristine gas. I implement this model within a large-scale cosmological simulation and follow the fraction of gas with metallicity below a critical value marking the boundary between Pop III and metal enriched Population II (Pop II) star formation. I demonstrate that accounting for subgrid mixing results in a Pop III stars formation rate that is 2-3 times higher than standard models with the same physical resolution.

I also implement and track a new "Primordial metals" (PM) scalar that tracks the metals generated by Pop III SNe. These metals are taken up by second generation stars and likely result in a subclass of carbon-enhanced, metal-poor (CEMP) stars. By tracking both regular metals and PM, I can model, in post-processing, the elemental abundances of simulation stars. I find good agreement between observations of CEMP-no Milky Way halo stars and second generation stars within the simulation when assuming the first stars had a typical mass of 60 M☉, providing clues as to the Pop III initial mass function.

Supernovae are vital to supplying necessary elements to forming bodies in our solar systems. This project studies the creation of a subset of these necessary elements, called short-lived radionuclides (SLRs). SLRs are isotopes with relatively short half-lives and can serve as heat sources for forming planetary bodies, and their traces can be used to date stellar events. Computational models of asymmetric supernovae provide opportunities to study the effect of explosion geometry on the SLR yields. We are most interested in the production of \iso{Al}{26}, \iso{Fe}{60}, and \iso{Ca}{41}, whose decayed products are found in our own solar system. To study the effect of explosion asymmetries in supernovae, we use TYCHO stellar evolution code, SNSHP smooth particle hydrodynamics code for 3D explosion simulations, Burn code for nucleosythesis post-processing, and Python code written to analyze the output of the post-processing code.