The movies and images below show projections of hydrodynamical properties of gas in simulation runs investigating different AGN feedback. The galaxy-gas evolution movie can be seen by clicking on the image or the words "Density/Temperature". The "Temperature" (T) movie actually shows the specific internal energy (Uspec), where, T = Uspec * (gamma - 1) * MeanWeight / k_B. The image is the end or an intermediate state of the gas for the duration in each movie. The simulations use a Hubble parameter of h = 0.72. The length scale is depicted by the 2 white circles with the radius written beside, whose unit is (1/h) kpc. The times written at the top-left are in units of 1.36 Gyr.
The initial condition is a disk galaxy of total mass 1.09e+12 Msun, having the following components: dark matter halo of mass 1.04e+12 Msun represented by 300000 collisionless particles, gas disk of mass 3.68e+09 Msun represented by 50000 SPH particles, stellar disk of mass 3.31e+10 Msun represented by 25000 collisionless particles, and stellar bulge of mass 9.20e+09 Msun represented by 25000 collisionless particles. It is same as the MW isolated galaxy in Murante, G. et al. (2010, MNRAS, 405, 1491). A Black Hole of mass written in the run description below is created at time = 0, at the location of the gas particle with maximum density. The galaxy is then allowed to evolve, with the BH accreting gas and its feedback being coupled as - either thermal or kinetic energy.
All runs have metallicity dependent radiative cooling, star formation, stellar evolution, and chemical enrichment (which includes thermal feedback from SNII & SNIa). The subgrid model of BH mass accretion has a multiplication factor of alpha = 100 in the Bondi accretion rate. In the kinetic feedback implementation, AGN wind particles always coupled to hydrodynamical interactions.The images and movie frames have been generated using Gadget File Viewer.
Run01-NoFeedback-alpha100
Sanity Check
Case without any BH feedback, only with gas accretion onto BH according to the Bondi rate multiplied by factor alpha. BH grows in mass uninterruptedly, depleting gas from the whole galaxy. No well-defined outflow is created, few particles flow out (leave the disk) at later times.
Face-On View, [X - Y] plane
Edge-On View, [Y - Z] plane
Run04-Kin-Enrg-epsF_0.01
Sanity Check
Kinetic feedback from BH, Energy-driven Wind prescription, eps_f = 0.01. Wind particles always coupled to hydrodynamical interactions. Intermittent bipolar jet-like outflows created, when a large number of gas particles flow out. Small hole forms at galaxy center by depletion of gas.
Temperature :
Thermal feedback from BH, eps_f = 0.01. Some particles flow out, & some of them fall back creating Galactic Fountain. BH does not grow as much as kinetic case.
Run21-NoBH-ForComparison : No BH. Only radiative cooling, star-formation, stellar evolution and chemical enrichment. Note a difference in the plotted color scale of Log(Internal Energy) between this case and the other runs shown in this page. This non-BH simulation demonstrates unambiguously that a prominent gas outflow develops, starting from time = 1.74 plotted units, via thermal SNIa feedback in the stellar evolution model. This SNIa-driven outflow is present in all the runs below (on top of BH-driven outflows), making it harder to disentangle the cause in those cases.
Run01-KinEnrg-v5k-e0.01 : Kinetic feedback from BH, Energy-driven Wind prescription, v_w = 5,000 km/s, eps_f = 0.01.
Run02-KinEnrg-v5k-e0.05 : Kinetic feedback from BH, Energy-driven Wind prescription, v_w = 5,000 km/s, eps_f = 0.05.
Run03-KinEnrg-v5k-e0.25 : Kinetic feedback from BH, Energy-driven Wind prescription, v_w = 5,000 km/s, eps_f = 0.25.
Run04-KinEnrg-v10k-e0.01 : Kinetic feedback from BH, Energy-driven Wind prescription, v_w = 10,000 km/s, eps_f = 0.01.
Run05-KinEnrg-v10k-e0.05 : Kinetic feedback from BH, Energy-driven Wind prescription, v_w = 10,000 km/s, eps_f = 0.05.
Run06-KinEnrg-v10k-e0.25 : Kinetic feedback from BH, Energy-driven Wind prescription, v_w = 10,000 km/s, eps_f = 0.25.
Run07-Thermal-e0.01 : Thermal feedback from BH, eps_f = 0.01.
Run08-Thermal-e0.05 : Thermal feedback from BH, eps_f = 0.05.
Run09-Thermal-e0.002 : Thermal feedback from BH, eps_f = 0.002.
Run10-KinMom-v5k-e0.25 : Kinetic feedback from BH, Momentum-driven Wind prescription, v_w = 5,000 km/s, eps_f = 0.25.
Run11-KinMom-v5k-e1.0 : Kinetic feedback from BH, Momentum-driven Wind prescription, v_w = 5,000 km/s, eps_f = 1.0.
Run12-KinMom-v2.5k-e1.0 : Kinetic feedback from BH, Momentum-driven Wind prescription, v_w = 2,500 km/s, eps_f = 1.0.
Run13-IC_vel300-KinE-v5k-e0.05 : Kinetic feedback from BH, Energy-driven Wind prescription, v_w = 5,000 km/s, eps_f = 0.05.
Run14-IC_vel300-KinE-v10k-e0.25 : Kinetic feedback from BH, Energy-driven Wind prescription, v_w = 10,000 km/s, eps_f = 0.25.
Run15-IC_vel300-Thermal-e0.01 : Thermal feedback from BH, eps_f = 0.01.
Run16-IC_vel300-KinM-v2.5k-e1.0 : Kinetic feedback from BH, Momentum-driven Wind prescription, v_w = 2,500 km/s, eps_f = 1.0.
Run17-IC_vel75-KinE-v5k-e0.05 : Kinetic feedback from BH, Energy-driven Wind prescription, v_w = 5,000 km/s, eps_f = 0.05.
Run18-IC_vel75-KinE-v10k-e0.25 : Kinetic feedback from BH, Energy-driven Wind prescription, v_w = 10,000 km/s, eps_f = 0.25.
Run19-IC_vel75-Thermal-e0.01 : Thermal feedback from BH, eps_f = 0.01.
Run20-IC_vel75-KinM-v2.5k-e1.0 : Kinetic feedback from BH, Momentum-driven Wind prescription, v_w = 2,500 km/s, eps_f = 1.0.