Traditional cosmological hydrodynamics simulations fail to spatially resolve the circumgalatic medium (CGM), the reservoir of tenuous gas surrounding a galaxy and extending to its virial radius. We introduce the technique of Enhanced Halo Resolution (EHR), enabling more realistic physical modeling of the simulated CGM by consistently forcing gas refinement to smaller scales throughout the virial halo of a simulated galaxy. We investigate the effects of EHR in the Tempest simulations, a suite of Enzo-based cosmological zoom simulations following the evolution of an L* galaxy, resolving spatial scales of 500 comoving pc out to 100 comoving kpc in galactocentric radius. Among its many effects, EHR (1) changes the thermal balance of the CGM, increasing its cool gas content and decreasing its warm/hot gas content; (2) preserves cool gas structures for longer periods; and (3) enables these cool clouds to exist at progressively smaller size scales. Observationally, this results in a boost in "low ions" like H I and a drop in "high ions" like O VI throughout the CGM. These effects of EHR do not converge in the Tempest simulations, but extrapolating these trends suggests that the CGM in reality is a mist consisting of ubiquitous, small, long-lived, cool clouds suspended in a hot medium at the virial temperature of the halo. Additionally, we explore the physical mechanisms to explain why EHR produces the above effects, proposing that it works both by (1) better sampling the distribution of CGM phases enabling runaway cooling in the denser, cooler tail of the phase distribution; and (2) preventing cool gas clouds from artificially mixing with the ambient hot halo and evaporating. Evidence is found for both EHR mechanisms occurring in the Tempest simulations.

We use a cosmological hydrodynamic simulation calculated with Enzo and the semi-analytic galaxy formation model (SAM) GAMMA to address the chemical evolution of dwarf galaxies in the early universe. The long-term goal of the project is to better understand the origin of metal-poor stars and the formation of dwarf galaxies and the Milky Way halo by cross-validating these theoretical approaches. We combine GAMMA with the merger tree of the most massive galaxy found in the hydrodynamic simulation and compare the star formation rate, the metallicity distribution function (MDF), and the age-metallicity relationship predicted by the two approaches. We found that the SAM can reproduce the global trends of the hydrodynamic simulation. However, there are degeneracies between the model parameters and more constraints (e.g., star formation efficiency, gas flows) need to be extracted from the simulation to isolate the correct semi-analytic solution. Stochastic processes such as bursty star formation histories and star formation triggered by supernova explosions cannot be reproduced by the current version of GAMMA. Non-uniform mixing in the galaxy's interstellar medium, coming primarily from self-enrichment by local supernovae, causes a broadening in the MDF that can be emulated in the SAM by convolving its predicted MDF with a Gaussian function having a standard deviation of ~0.2 dex. We found that the most massive galaxy in the simulation retains nearby 100% of its baryonic mass within its virial radius, which is in agreement with what is needed in GAMMA to reproduce the global trends of the simulation.

Hydrodynamical simulations are increasingly able to accurately model physical systems on stellar, galactic, and cosmological scales, however, the utility of these simulations is often limited by our ability to directly compare them with the datasets produced by observers: spectra, photometry, etc. To address this problem, we have created Trident}, a Python-based, open-source tool for post-processing hydrodynamical simulations to produce synthetic absorption spectra and related data. Trident} can (i) create absorption-line spectra for any trajectory through a simulated dataset mimicking both background quasar and down-the-barrel configurations, (ii) reproduce the spectral characteristics of common instruments like the Cosmic Origins Spectrograph, (iii) operate across the ultraviolet, optical and infrared using customizable absorption line lists, (iv) trace simulated physical structures directly to spectral features, (v) approximate the presence of ion species absent from the simulation outputs, (vi) generate column density maps for any ion, and (vii) provide support for all major astrophysical hydrodynamical codes. The focus of Trident's development is for using simulated datasets to better interpret observations of the circumgalactic medium (CGM) and intergalactic medium (IGM), but it remains a general tool applicable in other contexts.

Galactic chemical evolution is a topic that involves nuclear physics, stellar evolution, galaxy evolution, observation, and cosmology. Continuous communication and feedback between these fields is a key component in improving our understanding of how galaxies form and how elements are created and recycled in galaxies and intergalactic space. In this proceedings, we present the current state of the JINA-NuGrid chemical evolution pipeline. It is designed to probe the impact of nuclear astrophysics uncertainties on galactic chemical evolution, to improve our knowledges regarding the origin of the elements in a cosmological context, and to create the required interdisciplinary connections.

Aims. Adaptive optics images are used to test the hypothesis that the explosive BN/KL outflow from the Orion OMC1 cloud core was powered by the dynamical decay of a non-hierarchical system of massive stars. Methods. Narrow-band H2, [Fe II], and broad-band Ks obtained with the Gemini South multi-conjugate adaptive optics (AO) system GeMS and near-infrared imager GSAOI are presented. The images reach resolutions of 0.08 to 0.10", close to the 0.07" diffraction limit of the 8-meter telescope at 2.12 microns. Comparison with previous AO-assisted observations of sub-fields and other ground-based observations enable measurements of proper motions and the investigation of morphological changes in H2 and [Fe II] features with unprecedented precision. The images are compared with numerical simulations of compact, high-density clumps moving ~1000 times their own diameter through a lower density medium at Mach 1000. Results. Several sub-arcsecond H2 features and many [Fe II] 'fingertips' on the projected outskirts of the flow show proper motions of ~300 km/s. High-velocity, sub-arcsecond H2 knots ('bullets') are seen as far as 140" from their suspected ejection site. If these knots propagated through the dense Orion A cloud, their survival sets a lower bound on their densities of order 10^7 cm^-3, consistent with an origin within a few au of a massive star and accelerated by a final multi-body dynamic encounter that ejected the BN object and radio source I from OMC1 about 500 years ago. Conclusions. Over 120 high-velocity bow-shocks propagating in nearly all directions from the OMC1 cloud core provide evidence for an explosive origin for the BN/KL outflow triggered by the dynamic decay of a non-hierarchical system of massive stars. Such events may be linked to the origin of runaway, massive stars.

Following our previous work, we investigate through hydrodynamic simulations the destruction of newly-formed dust grains by sputtering in the reverse shocks of supernova remnants. Using an idealized setup of a planar shock impacting a dense, spherical clump, we implant a population of Lagrangian particles into the clump to represent a distribution of dust grains in size and composition. We vary the relative velocity between the reverse shock and ejecta clump to explore the effects of shock-heating and cloud compression. Because supernova ejecta will be metal-enriched, we consider gas metallicities from Z/Zsun = 1 to 100 and their influence on cooling properties of the cloud and the thermal sputtering rates of embedded dust grains. We post-process the simulation output to calculate grain sputtering for a variety of species and size distributions. In the metallicity regime considered in this paper, the balance between increased radiative cooling and increased grain erosion depends on the impact velocity of the reverse shock. For slow shocks (velocity less than or equal to 3000 km/s), the amount of dust destruction is comparable across metallicities, or in some cases is decreased with increased metallicity. For higher shock velocities (velocity greater than or equal to 5000 km/s), an increase in metallicity from Z/Zsun = 10 to 100 can lead to an additional 24% destruction of the initial dust mass. While the total dust destruction varies widely across grain species and simulation parameters, our most extreme cases result in complete destruction for some grain species and only 44% dust mass survival for the most robust species. These survival rates are important in understanding how early supernovae contribute to the observed dust masses in high-redshift galaxies.

Ejecta knot flickering, ablation tails, and fragmentation are expected signatures associated with the gradual dissolution of high-velocity supernova (SN) ejecta caused by their passage through an inhomogeneous circumstellar or interstellar medium. Such phenomena mark the initial stages of the gradual merger of SN ejecta with and the enrichment of the surrounding interstellar medium. Here we report on an investigation of this process through changes in the optical flux and morphology of several high-velocity ejecta knots located in the outskirts of the young core-collapse SN remnant Cassiopeia A using {\sl Hubble Space Telescope} images. Examination of WFPC2 F675W and combined ACS F625W + F775W images taken between June 1999 and December 2004 of several dozen debris fragments in the remnant's northeast ejecta stream and along the remnant's eastern limb reveal substantial emission variations ('flickering') over time scales as short as nine months. Such widespread and rapid variability indicates knot scale lengths ~ 10^15 cm and a highly inhomogeneous surrounding medium. We also identify a small percentage of ejecta knots located all around the remnant's outer periphery which show trailing emissions typically 0.2" to 0.7" in length aligned along the knot's direction of motion suggestive of knot ablation tails. We discuss the nature of these trailing emissions as they pertain to ablation cooling, knot disruption and fragmentation, and draw comparisons to the emission "strings" seen in eta Car. Finally, we identify several tight clusters of small ejecta knots which resemble models of shock induced fragmentation of larger SN ejecta knots caused by a high-velocity interaction with a lower density ambient medium.

We investigate through hydrodynamic simulations the destruction of newly-formed dust grains by sputtering in the reverse shocks of supernova remnants. Using an idealized setup of a planar shock impacting a dense, spherical clump, we implant a population of Lagrangian particles into the clump to represent a distribution of dust grains in size and composition. We then post-process the simulation output to calculate the grain sputtering for a variety of species and size distributions. We explore the parameter space appropriate for this problem by altering the over-density of the ejecta clumps and the speed of the reverse shocks. Since radiative cooling could lower the temperature of the medium in which the dust is embedded and potentially protect the dust by slowing or halting grain sputtering, we study the effects of different cooling methods over the time scale of the simulations. In general, our results indicate that grains with radii less than 0.1 microns are sputtered to much smaller radii and often destroyed completely, while larger grains survive their interaction with the reverse shock. We also find that, for high ejecta densities, the percentage of dust that survives is strongly dependent on the relative velocity between the clump and the reverse shock, causing up to 50% more destruction for the highest velocity shocks. The fraction of dust destroyed varies widely across grain species, ranging from total destruction of Al2O3 grains to minimal destruction of Fe grains (only 20% destruction in the most extreme cases). C and SiO2 grains show moderate to strong sputtering as well, with 38% and 80% mass loss. The survival rate of grains formed by early supernovae is crucial in determining whether or not they can act as the "dust factories" needed to explain high-redshift dust.

The periodic eclipses of the pre-main-sequence binary, KH 15D, have been explained by a circumbinary dust ring inclined to the orbital plane, which causes occultations of the stars as they pass behind the ring edge. We compute the extinction and forward scattering of light by the edge of the dust ring to explain (1) the gradual slope directly preceding total eclipse, (2) the gradual decline at the end of ingress, and (3) the slight rise in flux at mid-eclipse. The size of the forward scattering halo indicates that the dust grains have a radius of a ~ 6 (D/3 AU) microns, where D is the distance of the edge of the ring from the system barycenter. This dust size estimate agrees well with estimates of the dust grain size from polarimetry, adding to the evidence that the ring lies at several AU. Finally, the ratio of the fluxes inside and outside eclipse independently indicates that the ring lies at a few astronomical units.