Search

Galactic winds are crucial to the cosmic cycle of matter, transporting material out of the dense regions of galaxies. Observations show the coexistence of different temperature phases in such winds, which is not easy to explain. We present a set of 3D shock-multicloud simulations that account for radiative heating and cooling at temperatures between $10^2\,\rm K$ and $10^7\,\rm K$. The interplay between shock heating, dynamical instabilities, turbulence, and radiative heating and cooling creates a complex multi-phase flow with a rain-like morphology. Cloud gas fragments and is continuously eroded, becoming efficiently mixed and mass loaded. The resulting warm mixed gas then cools down and precipitates into new dense cloudlets, which repeat the process. Thus, radiative cooling is able to sustain fast-moving dense gas by aiding condensation of gas from warm clouds and the hot wind. In the ensuing outflow, hot gas with temperatures $\gtrsim 10^6\,\rm K$ outruns the warm and cold phases, which reach thermal equilibrium near $\approx 10^4\,\rm K$ and $\approx 10^2\,\rm K$, respectively. Although the volume filling factor of hot gas is higher in the outflow, most of the mass is concentrated in dense gas cloudlets and filaments with these temperatures. More porous multicloud layers result in more vertically extended outflows, and dense gas is more efficiently produced in more compact layers. The cold phase is not accelerated by ram-pressure, but, instead, precipitates from warm and mixed gas out of thermal equilibrium. This cycle can explain the presence of high-velocity H\,{\sc i} gas with $N_{\rm H\,{\scriptstyle I}}=10^{19-21}\,\rm cm^{-2}$ and $\Delta v_{{\rm FWHM}}\lesssim37\,\rm km\,s^{-1}$ in the Galactic centre outflow., Comment: Accepted for publication in MNRAS, 20 pages, 13 figures. This version includes an additional clump-finding analysis in the Appendix

Although galactic winds play a critical role in regulating galaxy formation, hydrodynamic cosmological simulations do not resolve the scales that govern the interaction between winds and the ambient circumgalactic medium (CGM). We implement the Physically Evolved Wind (PhEW) model of Huang et al. (2020) in the GIZMO hydrodynamics code and perform test cosmological simulations with different choices of model parameters and numerical resolution. PhEW adopts an explicit subgrid model that treats each wind particle as a collection of clouds that exchange mass, metals, and momentum with their surroundings and evaporate by conduction and hydrodynamic instabilities as calibrated on much higher resolution cloud scale simulations. In contrast to a conventional wind algorithm, we find that PhEW results are robust to numerical resolution and implementation details because the small scale interactions are defined by the model itself. Compared to conventional wind simulations with the same resolution, our PhEW simulations produce similar galaxy stellar mass functions at $z\geq 1$ but are in better agreement with low-redshift observations at $M_* < 10^{11}M_\odot$ because PhEW particles shed mass to the CGM before escaping low mass halos. PhEW radically alters the CGM metal distribution because PhEW particles disperse metals to the ambient medium as their clouds dissipate, producing a CGM metallicity distribution that is skewed but unimodal and is similar between cold and hot gas. While the temperature distributions and radial profiles of gaseous halos are similar in simulations with PhEW and conventional winds, these changes in metal distribution will affect their predicted UV/X-ray properties in absorption and emission., 23 pages, 17 figures, MNRAS accepted

We report three-dimensional hydrodynamical simulations of shocks (${\cal M_{\rm shock}}\geq 4$) interacting with fractal multicloud layers. The evolution of shock-multicloud systems consists of four stages: a shock-splitting phase in which reflected and refracted shocks are generated, a compression phase in which the forward shock compresses cloud material, an expansion phase triggered by internal heating and shock re-acceleration, and a mixing phase in which shear instabilities generate turbulence. We compare multicloud layers with narrow ($\sigma_{\rho}=1.9\bar{\rho}$) and wide ($\sigma_{\rho}=5.9\bar{\rho}$) log-normal density distributions characteristic of Mach $\approx 5$ supersonic turbulence driven by solenoidal and compressive modes. Our simulations show that outflowing cloud material contains imprints of the density structure of their native environments. The dynamics and disruption of multicloud systems depend on the porosity and the number of cloudlets in the layers. `Solenoidal' layers mix less, generate less turbulence, accelerate faster, and form a more coherent mixed-gas shell than the more porous `compressive' layers. Similarly, multicloud systems with more cloudlets quench mixing via a shielding effect and enhance momentum transfer. Mass loading of diffuse mixed gas is efficient in all models, but direct dense gas entrainment is highly inefficient. Dense gas only survives in compressive clouds, but has low speeds. If normalised with respect to the shock-passage time, the evolution shows invariance for shock Mach numbers $\geq10$ and different cloud-generating seeds, and slightly weaker scaling for lower Mach numbers and thinner cloud layers. Multicloud systems also have better convergence properties than single-cloud systems, with a resolution of $8$ cells per cloud radius being sufficient to capture their overall dynamics., Comment: Accepted for publication in MNRAS; 19 pages plus an appendix; 14 figures

Motivated by observations of outflowing galaxies, we investigate the combined impact of magnetic fields and radiative cooling on the evolution of cold clouds embedded in a hot wind. We perform a collection of three-dimensional adaptive mesh refinement, magnetohydrodynamical simulations that span two resolutions, and include fields that are aligned and transverse to the oncoming, super-Alfv\'enic material. Aligned fields have little impact on the overall lifetime of the clouds over the non-magnetized case, although they do increase the mixing between the wind and cloud material by a factor of $\approx 3.$ Transverse fields lead to magnetic draping, which isolates the clouds, but they also squeeze material in the direction perpendicular to the field lines, which leads to rapid mass loss. A resolution study suggests that the magnetized simulations have somewhat better convergence properties than non-magnetized simulations, and that a resolution of 64 zones per cloud radius is sufficient to accurately describe these interactions. We conclude that the combined effects of radiative cooling and magnetic fields are dependent on field orientation, but are unlikely to enhance cloud lifetimes beyond the effect of radiative cooling alone., Comment: 15 pages, 14 figures, accepted to ApJ

We use combined South Pole Telescope (SPT)+Planck temperature maps to analyze the circumgalactic medium (CGM) encompassing 138,235 massive, quiescent 0.5 ≤ z ≤ 1.5 galaxies selected from data from the Dark Energy Survey (DES) and Wide-Field Infrared Survey Explorer (WISE). Images centered on these galaxies were cut from the 1.85 arcmin resolution maps with frequency bands at 95, 150, and 220 GHz. The images were stacked, filtered, and fit with a graybody dust model to isolate the thermal Sunyaev–Zel’dovich (tSZ) signal, which is proportional to the total energy contained in the CGM of the galaxies. We separated these M ⋆ = 1010.9 M ⊙–1012 M ⊙ galaxies into 0.1 dex stellar mass bins, detecting tSZ per bin up to 5.6σ and a total signal-to-noise ratio of 10.1σ. We also detect dust with an overall signal-to-noise ratio of 9.8σ, which overwhelms the tSZ at 150 GHz more than in other lower-redshift studies. We corrected for the 0.16 dex uncertainty in the stellar mass measurements by parameter fitting for an unconvolved power-law energy-mass relation, E therm = E therm , peak M ⋆ / M ⋆ , peak α , with the peak stellar mass distribution of our selected galaxies defined as M ⋆,peak = 2.3 × 1011 M ⊙. This yields an E therm , peak = 5.98 − 1.00 + 1.02 × 10 60 erg and α = 3.77 − 0.74 + 0.60 . These are consistent with z ≈ 0 observations and within the limits of moderate models of active galactic nucleus feedback. We also computed the radial profile of our full sample, which is similar to that recently measured at lower-redshift by Schaan et al.

Absorption line studies are essential to understanding the origin, nature, and impact of starburst-driven galactic outflows. Such studies have revealed a multiphase medium with a number of poorly-understood features leading to a need to study the ionization mechanism of this gas. To better interpret these observations, we make use of a suite of adaptive mesh refinement hydrodynamic simulations of cold, atomic clouds driven by hot, supersonic outflows, including the effect of radiative cooling, thermal conduction, and an ionizing background characteristic of a starbursting galaxy. Using a new analysis tool, Trident, we estimate the equilibrium column density distributions for ten different ions: H I, Mg II, C II, C III, C IV, Si III, Si IV, N V, O VI, and Ne VIII. These are fit to model profiles with two parameters describing the maximum column density and coverage, and for each ion we provide a table of these fit parameters, along with average velocities and line widths. Our results are most sensitive to Mach number and conduction efficiency, with higher Mach numbers and more efficient conduction leading to more compact, high column density clouds. We use our results to interpret down-the-barrel observations of outflows and find that the adopted ionization equilibrium model overpredicts column densities of ions such as Si IV and does not adequately capture the observed trends for N V and O VI, implying the presence of strong non equilibrium ionization effects., Comment: 16 pages, 6 figures, accepted for publication in ApJ

The Orion Star Forming Complex (OSFC) is a central target for the APOGEE-2 Young Cluster Survey. Existing membership catalogs span limited portions of the OSFC, reflecting the difficulty of selecting targets homogeneously across this extended, highly structured region. We have used data from wide field photometric surveys to produce a less biased parent sample of young stellar objects (YSOs) with infrared (IR) excesses indicative of warm circumstellar material or photometric variability at optical wavelengths across the full 420 square degrees extent of the OSFC. When restricted to YSO candidates with H < 12.4, to ensure S/N ~100 for a six visit source, this uniformly selected sample includes 1307 IR excess sources selected using criteria vetted by Koenig & Liesawitz and 990 optical variables identified in the Pan-STARRS1 3$��$ survey: 319 sources exhibit both optical variability and evidence of circumstellar disks through IR excess. Objects from this uniformly selected sample received the highest priority for targeting, but required fewer than half of the fibers on each APOGEE-2 plate. We fill the remaining fibers with previously confirmed and new color-magnitude selected candidate OSFC members. Radial velocity measurements from APOGEE-1 and new APOGEE-2 observations taken in the survey's first year indicate that ~90% of the uniformly selected targets have radial velocities consistent with Orion membership.The APOGEE-2 Orion survey will include >1100 bona fide YSOs whose uniform selection function will provide a robust sample for comparative analyses of the stellar populations and properties across all sub-regions of Orion., Accepted for publication in ApJS