Intermediate and massive stars drive fast and powerful isotropic winds that interact with the winds of nearby stars in star clusters and the surrounding interstellar medium (ISM). Wind-ISM collisions generate astrospheres around these stars that contain hot $T\sim 10^7$ K gas that adiabatically expands. As individual bubbles expand and collide they become unstable, potentially driving turbulence in star clusters. In this paper we use hydrodynamic simulations to model a densely populated young star cluster within a homogeneous cloud to study stellar wind collisions with the surrounding ISM. We model a mass-segregated cluster of 20 B-type young main sequence stars with masses ranging from 3--17 $M_{\odot}$. We evolve the winds for $\sim$11 kyrs and show that wind-ISM collisions and over-lapping wind-blown bubbles around B-stars mixes the hot gas and ISM material generating Kolmogorov-like turbulence on small scales early in its evolution. We discuss how turbulence driven by stellar winds may impact the subsequent generation of star formation in the cluster

Scientific visualization tools are currently not optimized to create cinematic, production-quality representations of numerical data for the purpose of science communication. In our pipeline \texttt{Estra}, we outline a step-by-step process from a raw simulation into a finished render as a way to teach non-experts in the field of visualization how to achieve production-quality outputs on their own. We demonstrate feasibility of using the visual effects software Houdini for cinematic astrophysical data visualization, informed by machine learning clustering algorithms. To demonstrate the capabilities of this pipeline, we used a post-impact, thermally-equilibrated Moon-forming synestia from \cite{Lock18}. Our approach aims to identify "physically interpretable" clusters, where clusters identified in an appropriate phase space (e.g. here we use a temperature-entropy phase-space) correspond to physically meaningful structures within the simulation data. Clustering results can then be used to highlight these structures by informing the color-mapping process in a simplified Houdini software shading network, where dissimilar phase-space clusters are mapped to different color values for easier visual identification. Cluster information can also be used in 3D position space, via Houdini's Scene View, to aid in physical cluster finding, simulation prototyping, and data exploration. Our clustering-based renders are compared to those created by the Advanced Visualization Lab (AVL) team for the full dome show "Imagine the Moon" as proof of concept. With \texttt{Estra}, scientists have a tool to create their own production-quality, data-driven visualizations.

Supernova (SN) feedback plays a vital role in the evolution of galaxies. While modern cosmological simulations capture the leading structures within galaxies, they struggle to provide sufficient resolution to study small-scale stellar feedback, such as the detailed evolution of SN remnants. It is thus common practice to assume subgrid models that are rarely extended to low metallicities, and which routinely use the standard solar abundance pattern. With the aid of 1-d hydrodynamical simulations, we extend these models to consider low metallicities and non-solar abundance patterns as derived from spectra of Milky Way stars. For that purpose, a simple, yet effective framework has been developed to generate non-solar abundance pattern cooling functions. We find that previous treatments markedly over-predict SN feedback at low metallicities and show that non-negligible changes in the evolution of SN remnants of up to $\approx 50\%$ in $cooling\; mass$ and $\approx 27\%$ in $momentum\; injection\; from\; SN\; remnants$ arise from non-solar abundance patterns. We use our simulations to quantify these results as a function of metallicity and abundance pattern variations and present analytic formulae to accurately describe the trends. These formulae have been designed to serve as subgrid models for SN feedback in cosmological hydrodynamical simulations.

We employ hydrodynamical simulations to investigate the underlying mechanism responsible for the low levels of gas and dust in globular clusters. Our models examine the competing effects of energy and mass supply from the various components of the evolved stellar population for globular clusters 47 Tucanae, M15, NGC 6440, and NGC 6752. Ignoring all other gas evacuation processes, we find that the energy output from the stars that have recently turned off the main sequence are capable of effectively clearing the evolved stellar ejecta and producing intracluster gas densities consistent with current observational constraints. This result distinguishes a viable gas and dust evacuation mechanism that is ubiquitous in globular clusters. In addition, we extend our analysis to probe the efficiency of pulsar wind feedback in globular clusters. The detection of intracluster ionized gas in 47~Tucanae allows us to place particularly strict limits on pulsar wind thermalization efficiency, which must be extremely low in the cluster's core in order to be in accordance with the observed density constraints.

We analyze the IllustrisTNG simulations to study the mass, volume fraction and phase distribution of gaseous baryons embedded in the knots, filaments, sheets and voids of the Cosmic Web from redshift $z=8$ to redshift $z=0$. We find that filaments host more star-forming gas than knots, and that filaments also have a higher relative mass fraction of gas in this phase than knots. We also show that the cool, diffuse Intergalactic Medium (IGM; $T<10^5 \, {\rm K}$, $ n_{\rm H}<10^{-4}(1+z) \, {\rm cm^{-3}}$) and the Warm-Hot Intergalactic Medium (WHIM; $ 10^5 \, {\rm K} <T<10^7 \, {\rm K}$, $ n_{\rm H} <10^{-4}(1+z)\, {\rm cm^{-3}}$) constitute $\sim 39\%$ and $\sim 46\%$ of the baryons at redshift $z=0$, respectively. Our results indicate that the WHIM may constitute the largest reservoir of {\it missing} baryons at redshift $z=0$. Using our Cosmic Web classification, we predict the WHIM to be the dominant baryon mass contribution in filaments and knots at redshift $z=0$, but not in sheets and voids where the cool, diffuse IGM dominates. We also characterise the evolution of WHIM and IGM from redshift $z=4$ to redshift $z=0$, and find that the mass fraction of WHIM in filaments and knots evolves only by a factor $\sim 2$ from redshift $z=0$ to $z=1$, but declines faster at higher redshift. The WHIM only occupies $4-11\%$ of the volume at redshift $0\leq z \leq 1$. We predict the existence of a significant number of currently undetected OVII and NeIX absorption systems in cosmic filaments which could be detected by future X-ray telescopes like Athena.

The thermodynamic structure of hot gas in galaxy clusters is sensitive to astrophysical processes and typically difficult to model with galaxy formation simulations. We explore the fraction of cool-core (CC) clusters in a large sample of $370$ clusters from IllustrisTNG, examining six common CC definitions. IllustrisTNG produces continuous CC criteria distributions, the extremes of which are classified as CC and non-cool-core (NCC), and the criteria are increasingly correlated for more massive clusters. At $z=0$, the CC fractions for $2$ criteria are in reasonable agreement with the observed fractions but the other $4$ CC fractions are lower than observed. This result is partly driven by systematic differences between the simulated and observed gas fraction profiles. The simulated CC fractions with redshift show tentative agreement with the observed fractions, but linear fits demonstrate that the simulated evolution is steeper than observed. The conversion of CCs to NCCs appears to begin later and act more rapidly in the simulations. Examining the fraction of CCs and NCCs defined as relaxed we find no evidence that CCs are more relaxed, suggesting that mergers are not solely responsible for disrupting CCs. A comparison of the median thermodynamic profiles defined by different CC criteria shows that the extent to which they evolve in the cluster core is dependent on the CC criteria. We conclude that the thermodynamic structure of galaxy clusters in IllustrisTNG shares many similarities with observations, but achieving better agreement most likely requires modifications of the underlying galaxy formation model.

We introduce the IllustrisTNG project, a new suite of cosmological magnetohydrodynamical simulations performed with the moving-mesh code AREPO employing an updated Illustris galaxy formation model. Here we focus on the general properties of magnetic fields and the diffuse radio emission in galaxy clusters. Magnetic fields are prevalent in galaxies, and their build-up is closely linked to structure formation. We find that structure formation amplifies the initial seed fields ($10^{-14}$ comoving Gauss) to the values observed in low-redshift galaxies ($1-10\,\mu{\rm G}$). The magnetic field topology is closely connected to galaxy morphology such that irregular fields are hosted by early-type galaxies, while large-scale, ordered fields are present in disc galaxies. Using two simple models for the energy distribution of relativistic electrons we predict the diffuse radio emission of $280$ clusters with a baryonic mass resolution of $1.1\times 10^{7}\,{\rm M_{\odot}}$, and generate mock observations for VLA, LOFAR, ASKAP and SKA. Our simulated clusters show extended radio emission, whose detectability correlates with their virial mass. We reproduce the observed scaling relations between total radio power and X-ray emission, $M_{500}$, and the Sunyaev-Zel'dovich $Y_{\rm 500}$ parameter. The radio emission surface brightness profiles of our most massive clusters are in reasonable agreement with VLA measurements of Coma and Perseus. Finally, we discuss the fraction of detected extended radio haloes as a function of virial mass and source count functions for different instruments. Overall our results agree encouragingly well with observations, but a refined analysis requires a more sophisticated treatment of relativistic particles in large-scale galaxy formation simulations.

Mass and energy injection throughout the lifetime of a star cluster contributes to the gas reservoir available for subsequent episodes of star formation and the feedback energy budget responsible for ejecting material from the cluster. In addition, mass processed in stellar interiors and ejected as winds has the potential to augment the abundance ratios of currently forming stars, or stars which form at a later time from a retained gas reservoir. Here we present hydrodynamical simulations that explore a wide range of cluster masses, compactnesses, metallicities and stellar population age combinations in order to determine the range of parameter space conducive to stellar wind retention or wind powered gas expulsion in star clusters. We discuss the effects of the stellar wind prescription on retention and expulsion effectiveness, using MESA stellar evolutionary models as a test bed for exploring how the amounts of wind retention/expulsion depend upon the amount of mixing between the winds from stars of different masses and ages. We conclude by summarizing some implications for gas retention and expulsion in a variety of compact ($\sigma_v \gtrsim 20 \, {\rm km s^{-1}}$) star clusters including young massive star clusters ($10^5 \lesssim M/M_\odot \lesssim 10^7$, $age \lesssim 500$~Myrs), intermediate age clusters ($10^5 \lesssim M/M_\odot \lesssim 10^7$, $age \approx 1-4$~Gyrs), and globular clusters ($10^5 \lesssim M/M_\odot \lesssim 10^7$, $age \gtrsim 10$~Gyrs).

The distribution of elements in galaxies provides a wealth of information about their production sites and their subsequent mixing into the interstellar medium. Here we investigate the distribution of elements within stars in the IllustrisTNG simulations. In particular, we analyze the abundance ratios of magnesium and europium in Milky Way-like galaxies from the TNG100 simulation (stellar masses ${\log} (M_\star / {\rm M}_\odot) \sim 9.7 - 11.2$). As abundances of magnesium and europium for individual stars in the Milky Way are observed across a variety of spatial locations and metallicities, comparison with the stellar abundances in our more than $850$ Milky Way-like galaxies provides stringent constraints on our chemical evolutionary methods. To this end we use the magnesium to iron ratio as a proxy for the effects of our SNII and SNIa metal return prescription, and a means to compare our simulated abundances to a wide variety of galactic observations. The europium to iron ratio tracks the rare ejecta from neutron star -- neutron star mergers, the assumed primary site of europium production in our models, which in turn is a sensitive probe of the effects of metal diffusion within the gas in our simulations. We find that europium abundances in Milky Way-like galaxies show no correlation with assembly history, present day galactic properties, and average galactic stellar population age. In general, we reproduce the europium to iron spread at low metallicities observed in the Milky Way, with the level of enhancement being sensitive to gas properties during redshifts $z \approx 2-4$. We show that while the overall normalization of [Eu/Fe] is susceptible to resolution and post-processing assumptions, the relatively large spread of [Eu/Fe] at low [Fe/H] when compared to that at high [Fe/H] is very robust.

We introduce the first two simulations of the IllustrisTNG project, a next generation of cosmological magnetohydrodynamical simulations, focusing on the optical colors of galaxies. We explore TNG100, a rerun of the original Illustris box, and TNG300, which includes 2x2500^3 resolution elements in a volume twenty times larger. Here we present first results on the galaxy color bimodality at low redshift. Accounting for the attenuation of stellar light by dust, we compare the simulated (g-r) colors of 10^9 < M*/Msun < 10^12.5 galaxies to the observed distribution from the Sloan Digital Sky Survey (SDSS). We find a striking improvement with respect to the original Illustris simulation, as well as excellent quantitative agreement in comparison to the observations, with a sharp transition in median color from blue to red at a characteristic M* ~ 10^10.5 Msun. Investigating the build-up of the color-mass plane and the formation of the red sequence, we demonstrate that the primary driver of galaxy color transition in the TNG model is supermassive blackhole feedback in its low-accretion state. Across the entire population we measure a median color transition timescale dt_green of ~1.6 Gyr, a value which drops for increasingly massive galaxies. We find signatures of the physical process of quenching: at fixed stellar mass, redder galaxies have lower SFRs, gas fractions, and gas metallicities; their stellar populations are also older and their large-scale interstellar magnetic fields weaker than in bluer galaxies. Finally, we measure the amount of stellar mass growth on the red sequence. Galaxies with M* > 10^11 Msun which redden at z<1 accumulate on average ~25% of their final z=0 mass post-reddening; at the same time, ~18% of such massive galaxies acquire half or more of their final stellar mass while on the red sequence.

The IllustrisTNG project is a new suite of cosmological magneto-hydrodynamical simulations of galaxy formation performed with the Arepo code and updated models for feedback physics. Here we introduce the first two simulations of the series, TNG100 and TNG300, and quantify the stellar mass content of about 4000 massive galaxy groups and clusters ($10^{13} \leq M_{\rm 200c}/M_{\rm sun} \leq 10^{15}$) at recent times ($z \leq 1$). The richest clusters have half of their total stellar mass bound to satellite galaxies, with the other half being associated with the central galaxy and the diffuse intra-cluster light. The exact ICL fraction depends sensitively on the definition of a central galaxy's mass and varies in our most massive clusters between 20 to 40% of the total stellar mass. Haloes of $5\times 10^{14}M_{\rm sun}$ and above have more diffuse stellar mass outside 100 kpc than within 100 kpc, with power-law slopes of the radial mass density distribution as shallow as the dark matter's ( $-3.5 < \alpha_{\rm 3D} < -3$). Total halo mass is a very good predictor of stellar mass, and vice versa: at $z=0$, the 3D stellar mass measured within 30 kpc scales as $\propto (M_{\rm 500c})^{0.49}$ with a $\sim 0.12$ dex scatter. This is possibly too steep in comparison to the available observational constraints, even though the abundance of TNG less massive galaxies ($< 10^{11}M_{\rm sun}$ in stars) is in good agreement with the measured galaxy stellar mass functions at recent epochs. The 3D sizes of massive galaxies fall too on a tight ($\sim$0.16 dex scatter) power-law relation with halo mass, with $r^{\rm stars}_{\rm 0.5} \propto (M_{\rm 500c})^{0.53}$. Even more fundamentally, halo mass alone is a good predictor for the whole stellar mass profiles beyond the inner few kpc, and we show how on average these can be precisely recovered given a single mass measurement of the galaxy or its halo.

Hydrodynamical simulations of galaxy formation have now reached sufficient volume to make precision predictions for clustering on cosmologically relevant scales. Here we use our new IllustrisTNG simulations to study the non-linear correlation functions and power spectra of baryons, dark matter, galaxies and haloes over an exceptionally large range of scales. We find that baryonic effects increase the clustering of dark matter on small scales and damp the total matter power spectrum on scales up to k ~ 10 h/Mpc by 20%. The non-linear two-point correlation function of the stellar mass is close to a power-law over a wide range of scales and approximately invariant in time from very high redshift to the present. The two-point correlation function of the simulated galaxies agrees well with SDSS at its mean redshift z ~ 0.1, both as a function of stellar mass and when split according to galaxy colour, apart from a mild excess in the clustering of red galaxies in the stellar mass range 10^9-10^10 Msun/h^2. Given this agreement, the TNG simulations can make valuable theoretical predictions for the clustering bias of different galaxy samples. We find that the clustering length of the galaxy auto-correlation function depends strongly on stellar mass and redshift. Its power-law slope gamma is nearly invariant with stellar mass, but declines from gamma ~ 1.8 at redshift z=0 to gamma ~ 1.6 at redshift z ~ 1, beyond which the slope steepens again. We detect significant scale-dependencies in the bias of different observational tracers of large-scale structure, extending well into the range of the baryonic acoustic oscillations and causing nominal (yet fortunately correctable) shifts of the acoustic peaks of around ~5%.

The fundamental metallicity relation (FMR) is a postulated correlation between galaxy stellar mass, star formation rate (SFR), and gas-phase metallicity. At its core, this relation posits that offsets from the mass-metallicity relation (MZR) at a fixed stellar mass are correlated with galactic SFR. In this Letter, we quantify the timescale with which galactic SFRs and metallicities evolve using hydrodynamical simulations. We find that Illustris and IllustrisTNG predict that galaxy offsets from the star formation main sequence and MZR evolve over similar timescales, are often anti-correlated in their evolution, evolve with the halo dynamical time, and produce a pronounced FMR. In fact, for a FMR to exist, the metallicity and SFR must evolve in an anti-correlated sense which requires that they evolve with similar time variability. In contrast to Illustris and IllustrisTNG, we speculate that the SFR and metallicity evolution tracks may become decoupled in galaxy formation models dominated by globally-bursty SFR histories, which could weaken the FMR residual correlation strength. This opens the possibility of discriminating between bursty and non-bursty feedback models based on the strength and persistence of the FMR -- especially at high redshift.

We analyze scaling relations and evolution histories of galaxy sizes in TNG100, part of the IllustrisTNG simulation suite. Observational qualitative trends of size with stellar mass, star-formation rate and redshift are reproduced, and a quantitative comparison of projected r-band sizes at 0~<z<~2 shows agreement to much better than 0.25dex. We follow populations of z=0 galaxies with a range of masses backwards in time along their main progenitor branches, distinguishing between main-sequence and quenched galaxies. Our main findings are as follows. (i) At M_{*,z=0}>~10^{9.5}Msun, the evolution of the median main progenitor differs, with quenched galaxies hardly growing in median size before quenching, whereas main-sequence galaxies grow their median size continuously, thus opening a gap from the progenitors of quenched galaxies. This is partly because the main-sequence high-redshift progenitors of quenched z=0 galaxies are drawn from the lower end of the size distribution of the overall population of main-sequence high-redshift galaxies. (ii) Quenched galaxies with M_{*,z=0}>~10^{9.5}Msun experience a steep size growth on the size-mass plane after their quenching time, but with the exception of galaxies with M_{*,z=0}>~10^{11}Msun, the size growth after quenching is small in absolute terms, such that most of the size (and mass) growth of quenched galaxies (and its variation among them) occurs while they are still on the main-sequence. After they become quenched, the size growth rate of quenched galaxies as a function of time, as opposed to versus mass, is similar to that of main-sequence galaxies. Hence, the size gap is retained down to z=0.

The coevolution of galaxies and their metal content serves as an important test for galaxy feedback models. We analyze the distribution and evolution of metals within the IllustrisTNG simulation suite with a focus on the gas-phase mass-metallicity relation (MZR). We find that the IllustrisTNG model broadly reproduces the slope and normalization evolution of the MZR across the redshift range $0<z<2$ and mass range $10^9 < M_*/\mathrm{M}_\odot < 10^{10.5}$. We make predictions for the high redshift ($2<z<10$) metal content of galaxies which is described by a gradual decline in the normalization of the metallicity with an average high redshift ($z>2$) evolution fit by $\mathrm{d\;log(Z)}/\mathrm{dz} \approx - 0.064$. Our simulations indicate that the metal retention efficiency of the interstellar medium (ISM) is low: a majority of gas-phase metals ($\sim$ 85 per cent at $z=0$) live outside of the ISM, either in an extended gas disk, the circumgalactic medium, or outside the halo. Nevertheless, the redshift evolution in the simulated MZR normalization is driven by the higher gas fractions of high redshift galaxies, not by changes to the metal retention efficiency. The scatter in the simulated MZR contains a clear correlation with the gas-mass or star formation rate of the system, in agreement with the observed fundamental metallicity relation. The scatter in the MZR is driven by a competition between periods of enrichment- and accretion-dominated metallicity evolution. We expect that while the normalization of the MZR declines with redshift, the strength of the correlation between metallicity and gas-mass at fixed stellar mass is not a strong function of redshift. Our results indicate that the "regulator" style models are best suited for simultaneously explaining the shape, redshift evolution, and existence of correlated scatter with gas fraction about the MZR.

The distribution of metals in the intra-cluster medium encodes important information about the enrichment history and formation of galaxy clusters. Here we explore the metal content of clusters in IllustrisTNG - a new suite of galaxy formation simulations building on the Illustris project. Our cluster sample contains 20 objects in TNG100 - a ~(100 Mpc)^3 volume simulation with 2x1820^3 resolution elements, and 370 objects in TNG300 - a ~(300 Mpc)^3 volume simulation with 2x2500^3 resolution elements. The z=0 metallicity profiles agree with observations, and the enrichment history is consistent with observational data going beyond z~1, showing nearly no metallicity evolution. The abundance profiles vary only minimally within the cluster samples, especially in the outskirts with a relative scatter of ~15%. The average metallicity profile flattens towards the center, where we find a logarithmic slope of -0.1 compared to -0.5 in the outskirts. Cool core clusters have more centrally peaked metallicity profiles (~0.8 solar) compared to non-cool core systems (~0.5 solar), similar to observational trends. Si/Fe and O/Fe radial profiles follow positive gradients. The outer abundance profiles do not evolve below z~2, whereas the inner profiles flatten towards z=0. More than ~80% of the metals in the intra-cluster medium have been accreted from the proto-cluster environment, which has been enriched to ~0.1 solar already at z~2. We conclude that the intra-cluster metal distribution is uniform among our cluster sample, nearly time-invariant in the outskirts for more than 10 Gyr, and forms through a universal enrichment history.

We introduce an updated physical model to simulate the formation and evolution of galaxies in cosmological, large-scale gravity+magnetohydrodynamical simulations with the moving mesh code AREPO. The overall framework builds upon the successes of the Illustris galaxy formation model, and includes prescriptions for star formation, stellar evolution, chemical enrichment, primordial and metal-line cooling of the gas, stellar feedback with galactic outflows, and black hole formation, growth and multi-mode feedback. In this paper we give a comprehensive description of the physical and numerical advances which form the core of the IllustrisTNG (The Next Generation) framework. We focus on the revised implementation of the galactic winds, of which we modify the directionality, velocity, thermal content, and energy scalings, and explore its effects on the galaxy population. As described in earlier works, the model also includes a new black hole driven kinetic feedback at low accretion rates, magnetohydrodynamics, and improvements to the numerical scheme. Using a suite of (25 Mpc $h^{-1}$)$^3$ cosmological boxes we assess the outcome of the new model at our fiducial resolution. The presence of a self-consistently amplified magnetic field is shown to have an important impact on the stellar content of $10^{12} M_{\rm sun}$ haloes and above. Finally, we demonstrate that the new galactic winds promise to solve key problems identified in Illustris in matching observational constraints and affecting the stellar content and sizes of the low mass end of the galaxy population.

We study the population of supermassive black holes (SMBHs) and their effects on massive central galaxies in the IllustrisTNG cosmological hydrodynamical simulations of galaxy formation. The employed model for SMBH growth and feedback assumes a two-mode scenario in which the feedback from active galactic nuclei occurs through a kinetic, comparatively efficient mode at low accretion rates relative to the Eddington limit, and in the form of a thermal, less efficient mode at high accretion rates. We show that the quenching of massive central galaxies happens coincidently with kinetic-mode feedback, consistent with the notion that active supermassive black cause the low specific star formation rates observed in massive galaxies. However, major galaxy mergers are not responsible for initiating most of the quenching events in our model. Up to black hole masses of about $10^{8.5}\,{\rm M}_\odot$, the dominant growth channel for SMBHs is in the thermal mode. Higher mass black holes stay mainly in the kinetic mode and gas accretion is self-regulated via their feedback, which causes their Eddington ratios to drop, with SMBH mergers becoming the main channel for residual mass growth. As a consequence, the quasar luminosity function is dominated by rapidly accreting, moderately massive black holes in the thermal mode. We show that the associated growth history of SMBHs produces a low-redshift quasar luminosity function and a redshift zero black hole mass-stellar bulge mass relation in good agreement with observations, whereas the simulation tends to over-predict the high-redshift quasar luminosity function.

The inefficiency of star formation in massive elliptical galaxies is widely believed to be caused by the interactions of an active galactic nucleus (AGN) with the surrounding gas. Achieving a sufficiently rapid reddening of moderately massive galaxies without expelling too many baryons has however proven difficult for hydrodynamical simulations of galaxy formation, prompting us to explore a new model for the accretion and feedback effects of supermassive black holes. For high accretion rates relative to the Eddington limit, we assume that a fraction of the accreted rest mass energy heats the surrounding gas thermally, similar to the `quasar mode' in previous work. For low accretion rates, we invoke a new, pure kinetic feedback model which imparts momentum into the surrounding gas in a stochastic manner. These two modes of feedback are motivated both by theoretical conjectures for the existence of different types of accretion flows as well as recent observational evidence for the importance of kinetic AGN winds in quenching galaxies. We find that a large fraction of the injected kinetic energy in this mode thermalises via shocks in the surrounding gas, thereby providing a distributed heating channel. In cosmological simulations, the resulting model produces red, non star-forming massive elliptical galaxies, and achieves realistic gas fractions, black hole growth histories and thermodynamic profiles in large haloes.

The rapid growth in scale and complexity of both computational and observational astrophysics over the past decade necessitates efficient and intuitive methods for examining and visualizing large datasets. Here we discuss some newly developed tools to import and manipulate astrophysical data into the three dimensional visual effects software, {\it Houdini}. This software is widely used by visual effects artists, but a recently implemented Python API now allows astronomers to more easily use Houdini as a visualization tool. This paper includes a description of features, work flow, and various example visualizations. The project website, www.ytini.com, contains Houdini tutorials and links to the Python script Bitbucket repository aimed at a scientific audience to simplify the process of importing and rendering astrophysical data.

The rapid growth in scale and complexity of both computational and observational astrophysics over the past decade necessitates efficient and intuitive methods for examining and visualizing large datasets. Here, I present {\it AstroBlend}, an open-source Python library for use within the three dimensional modeling software, {\it Blender}. While {\it Blender} has been a popular open-source software among animators and visual effects artists, in recent years it has also become a tool for visualizing astrophysical datasets. {\it AstroBlend} combines the three dimensional capabilities of {\it Blender} with the analysis tools of the widely used astrophysical toolset, {\it yt}, to afford both computational and observational astrophysicists the ability to simultaneously analyze their data and create informative and appealing visualizations. The introduction of this package includes a description of features, work flow, and various example visualizations. A website - www.astroblend.com - has been developed which includes tutorials, and a gallery of example images and movies, along with links to downloadable data, three dimensional artistic models, and various other resources.

The r-process nuclei are robustly synthesized in the material ejected during a neutron star binary merger (NSBM), as tidal torques transport angular momentum and energy through the outer Lagrange point in the form of a vast tidal tail. If NSBM are indeed solely responsible for the solar system r- process abundances, a galaxy like our own would require to host a few NSBM per million years, with each event ejecting, on average, about 5x10^{-2} M_sun of r-process material. Because the ejecta velocities in the tidal tail are significantly larger than in ordinary supernovae, NSBM deposit a comparable amount of energy into the interstellar medium (ISM). In contrast to extensive efforts studying spherical models for supernova remnant evolution, calculations quantifying the impact of NSBM ejecta in the ISM have been lacking. To better understand their evolution in a cosmological context, we perform a suite of three-dimensional hydrodynamic simulations with optically-thin radiative cooling of isolated NSBM ejecta expanding in environments with conditions adopted from Milky Way-like galaxy simulations. Although the remnant morphology is highly complex at early times, the subsequent radiative evolution that results from thermal instability in atomic gas is remarkably similar to that of a standard supernova blast wave. This implies that sub-resolution supernova feedback models can be used in galaxy-scale simulations that are unable to resolve the key evolutionary phases of NSBM blast waves. Among other quantities, we examine the radius, time, mass and kinetic energy content of the NSBM remnant at shell formation as well as the momentum injected to the ISM. We find that the shell formation epoch is attained when the swept-up mass is about 10^3 M_sun, at this point the mass fraction of r-process material is drastically enhanced up to two orders of magnitude in relation to a solar metallicity ISM.

During galaxy mergers the gas falls to the center, triggers star formation, and feeds the rapid growth of supermassive black holes (SMBHs). SMBHs respond to this fueling by supplying energy back to the ambient gas. Numerical studies suggest that this feedback is necessary to explain why the properties of SMBHs and the formation of bulges are closely related. This intimate link between the SMBH's mass and the large scale dynamics and luminosity of the host has proven to be a difficult issue to tackle with simulations due to the inability to resolve all the relevant length scales simultaneously. In this paper we simulate SMBH growth at high-resolution with {\it FLASH}, accounting for the gravitational focusing effects of nuclear star clusters (NSCs), which appear to be ubiquitous in galactic nuclei. In the simulations, the NSC core is resolved by a minimum cell size of about 0.001 pc or approximately $10^{-3}$ of the cluster's radius. We discuss the conditions required for effective gas funneling to occur, which are mainly dominated by a relationship between NSC velocity dispersion and the local sound speed, and provide a sub-grid prescription for the augmentation of central SMBH accretion rates in the presence of NSCs. For the conditions expected to persist in the centers of merging galaxies, the resultant large central gas densities in NSCs should produce drastically enhanced embedded SMBH accretion rates - up to an order of magnitude increase can be achieved for gas properties resembling those in large-scale galaxy merger simulations. This will naturally result in faster black hole growth rates and higher luminosities than predicted by the commonly used Bondi-Hoyle-Lyttleton accretion formalism.

We examine the consequences of a model in which relativistic jets can be triggered in quiescent massive black holes when a geometrically thick and hot accretion disk forms as a result of the tidal disruption of a star. To estimate the power, thrust and lifetime of the jet, we use the mass accretion history onto the black hole as calculated by detailed hydrodynamic simulations of the tidal disruption of stars. We go on to determine the states of the interstellar medium in various types of quiescent galactic nuclei, and describe how this external matter can affect jets propagating through it. We use this information, together with a two-dimensional hydrodynamic model of the structure of the relativistic flow, to study the dynamics of the jet, the propagation of which is regulated by the density stratification of the environment and by its injection history. The breaking of symmetry involved in transitioning from one to two dimensions is crucial and leads to qualitatively new phenomena. Many of the observed properties of the Swift 1644+57/GRB 110328A event can be understood as resulting from accretion onto and jets driven by a $10^6 M_\odot$ central mass black hole following the disruption of sun-like star. With the inclusion of a stochastic contribution to the luminosity due to variations in the feeding rate driven by instabilities near the tidal radius, we find that our model can explain the X-ray light curve without invoking a rarely-occurring deep encounter. In conjunction with the number density of black holes in the local universe, we hypothesize that the conditions required to produce the Swift event are not anomalous, but are in fact representative of the jet-driven flare population arising from tidal disruptions. [abridged]

Accretion studies have been focused on the flow around bodies with point mass gravitational potentials, but few general results are available for non-point mass distributions. Here, we study the accretion flow onto non-divergent, core potentials moving through a background medium. We use Plummer and Hernquist potentials as examples to study gas accretion onto star clusters, dwarf and large galaxy halos and galaxy clusters in a variety of astrophysical environments. The general conditions required for a core potential to collectively accrete large quantities of gas from the external medium are derived using both simulations and analytic results. The consequences of large mass accumulation in galaxy nuclei, dwarf galaxies and star clusters are twofold. First, if the gas cools effectively star formation can be triggered, generating new stellar members in the system. Second, if the collective potential of the system is able to alter the ambient gas properties before the gas is accreted onto the individual core members, the augmented mass supply rates could significantly alter the state of the various accreting stellar populations and result in an enhanced central black hole accretion luminosity.

Here we show that the overabundance of ultra-luminous, compact X-ray sources (ULXs) associated with moderately young clusters in interacting galaxies such as the Antennae and Cartwheel can be given an alternative explanation that does not involve the presence of intermediate mass black holes (IMBHs). We argue that gas density within these systems is enhanced by the collective potential of the cluster prior to being accreted onto the individual cluster members and, as a result, the aggregate X-ray luminosity arising from the neutron star cluster members can exceed $>10^{39} {\rm erg s^{-1}}$. Various observational tests to distinguish between IMBHs and accreting neutron star cusps are discussed.

We report new precision measurements of the properties of our Galaxy's supermassive black hole. Based on astrometric (1995-2007) and radial velocity (2000-2007) measurements from the W. M. Keck 10-meter telescopes, a fully unconstrained Keplerian orbit for the short period star S0-2 provides values for Ro of 8.0+-0.6 kpc, M_bh of 4.1+-0.6x10^6 Mo, and the black hole's radial velocity, which is consistent with zero with 30 km/s uncertainty. If the black hole is assumed to be at rest with respect to the Galaxy, we can further constrain the fit and obtain Ro = 8.4+-0.4 kpc and M_bh = 4.5+-0.4x10^6 Mo. More complex models constrain the extended dark mass distribution to be less than 3-4x10^5 Mo within 0.01 pc, ~100x higher than predictions from stellar and stellar remnant models. For all models, we identify transient astrometric shifts from source confusion and the assumptions regarding the black hole's radial motion as previously unrecognized limitations on orbital accuracy and the usefulness of fainter stars. Future astrometric and RV observations will remedy these effects. Our estimates of Ro and the Galaxy's local rotation speed, which it is derived from combining Ro with the apparent proper motion of Sgr A*, (theta0 = 229+-18 km/s), are compatible with measurements made using other methods. The increased black hole mass found in this study, compared to that determined using projected mass estimators, implies a longer period for the innermost stable orbit, longer resonant relaxation timescales for stars in the vicinity of the black hole and a better agreement with the M_bh-sigma relation.