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 µm 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 Al 2 O 3 grains to minimal destruction of Fe grains (only 20% destruction in the most extreme cases). C and SiO 2 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.

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.

Aims. Adaptive optics (AO) 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. Results. Several sub-arcsecond H 2 features and many [Fe ] "fingertips" on the projected outskirts of the flow show proper motions of ∼300 km s −1 . High-velocity, sub-arcsecond H 2 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 yr 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/Z = 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 (v shock ≤ 3000 km s −1 ), the amount of dust destruction is comparable across metallicities, or in some cases is decreased with increased metallicity. For higher shock velocities (v shock ≥ 5000 km s −1 ), an increase in metallicity from Z/Z = 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.

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 data sets 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 data set 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. TRIDENT was originally developed to aid in the interpretation of observations of the circumgalactic medium and intergalactic medium, but it remains a general tool applicable in other contexts.