We produce radiation hydrodynamics models of an AGN 'torus' plus outflow on 1-100 pc scales. This large scale permits direct comparison with observations, provides justification for configurations used in radiation transfer models, and tests the sensitivity of results of smaller scale dynamical models. We find that anisotropic radiation from an AGN accretion disk can cause an outflow to evolve to become more polar, agreeing with the ubiquity of polar extended mid-IR emission, and the general geometry predicted by radiative transfer models. We also find the velocity maps can reproduce many features of observations, including apparent 'counter-rotation'.

Infrared interferometry of local AGN has revealed a warm (~300K-400K) polar dust structure that cannot be trivially explained by the putative dust torus of the unified model. This led to the development of the disk+wind scenario which comprises of a hot (~1000K) compact equatorial dust disk and a polar dust wind. This wind is assumed to be driven by radiation pressure and, therefore, we would expect that long term variation in radiation pressure would influence the dust distribution. In this paper we attempt to quantify if and how the dust distribution changes with radiation pressure. We analyse so far unpublished VLTI/MIDI data on 8 AGN and use previous results on 25 more to create a sample of 33 AGN. This sample comprises all AGN successfully observed with VLTI/MIDI. For each AGN, we calculate the Eddington ratio, using the intrinsic 2-10keV X-ray luminosity and black hole mass, and compare this to the resolved dust emission fraction as seen by MIDI. We tentatively conclude that there is more dust in the wind at higher Eddington ratios, at least in type 2 AGN where such an effect is expected to be more easily visible.

We have developed a new dynamical model of the torus region in active galactic nucleus (AGN), using a three-dimensional radiation hydrodynamics algorithm. These new simulations have the specific aim to explore the role of radiatively-driven outflows, which is hotly debated in current literature as a possible explanation for the observed infrared emission from the polar regions of AGN. In this first paper, we only consider radiative effects induced by the primary radiation from the AGN. The simulations generate a disk & outflow structure that qualitatively agrees with observations, although the outflow is radial rather than polar, likely due to the lack of radiation pressure from hot dust. We find cut-offs between the wind and disk at gas temperatures of 1000 K and dust temperatures of 100 K, producing kinematic signatures that can be used for interpretation of high resolution infrared observations. We also produce line emission maps to aid in the interpretation of recent ALMA observations and future JWST observations. We investigate a number of simulation parameters, and find that the anisotropy of the radiation field is equally important to the Eddington factor, despite the anisotropy often being assumed to have a single sometimes arbitrary form in many previous works. We also find that supernovae can have a small but significant impact, but only at extremely high star formation rates.

By implementing a dynamic wind-tunnel model in a smoothed-particle chemodynamic/hydrodynamic simulation suite, we have investigated the effects of ram pressure and tidal forces on dwarf galaxies similar to the Magellanic Clouds, within host galaxies with gas and dark matter halos that are varied, to compare the relative effects of tides and ram pressure. We concentrate on how the distributions of metals are affected by interactions. We find that while ram pressure and tidal forces have some effect on dwarf galaxy outflows, these effects do not produce large differences in the metal distributions of the dwarf disks other than truncation in the outer regions in some cases, and that confinement from the host galaxy gas halo appears to be more significant than ram pressure stripping. We find that stochastic variations in the star formation rate can explain the remaining variations in disk metal properties. This raises questions on the cause of low metallicities in dwarf galaxies.

To investigate the role of bars in the chemical evolution of isolated disc galaxies, we performed a series of 39 gas dynamical simulations of isolated barred and unbarred galaxies with various masses, initial gas fractions, and AGN feedback models. The presence of a bar drives a substantial amount of gas toward the central region of the galaxy. In the most massive galaxies, this results in a violent starburst, followed by a drop in star formation resulting from gas exhaustion. The time delay between Type Ia and Type II supernovae explosions means that barred galaxies experience a rapid increase in [O/H] in the central region, and a much more gradual increase in [Fe/H]. In unbarred galaxies, star formation proceeds at a slow and steady rate, and oxygen and iron are produced at steady rates which are similar except for a time offset. Comparing the abundance ratios in barred and unbarred galaxies with the same central stellar mass $M_*$, we find in barred galaxies an enhancement of 0.07 dex in [O/H], 0.05 dex in [Fe/H], and 0.05 dex in [O/Fe]. The [O/H] enhancement is in excellent agreement with observations from the SDSS. The initial gas fraction has very little effect on the abundance ratios in barred and unbarred galaxies, unless the galaxies experience a starburst. We considered AGN-host galaxies located near the bottom of the AGN regime, $M_*>3\times10^{10}M_\odot$, where AGN feedback dominates over supernovae feedback. We found that the impact of AGN feedback on the central abundances is marginal.

We present a numerical study of the impact of AGN accretion and feedback on the star formation history of barred disc galaxies. Our goal is to determine whether the effect of feedback is positive (enhanced star formation) or negative (quenched star formation), and to what extent. We performed a series of 12 hydrodynamical simulations of disc galaxies, 10 barred and 2 unbarred, with various initial gas fractions and AGN feedback prescriptions. In barred galaxies, gas is driven toward the centre of the galaxy and causes a starburst, followed by a slow decay, while in unbarred galaxies the SFR increases slowly and steadily. AGN feedback suppresses star formation near the central black hole. Gas is pushed away from the black hole, and collides head-on with inflowing gas, forming a dense ring at a finite radius where star formation is enhanced. We conclude that both negative and positive feedback are present, and these effects mostly cancel out. There is no net quenching or enhancement in star formation, but rather a displacement of the star formation sites to larger radii. In unbarred galaxies, where the density of the central gas is lower, quenching of star formation near the black hole is more efficient, and enhancement of star formation at larger radii is less efficient. As a result, negative feedback dominates. Lowering the gas fraction reduces the star formation rate at all radii, whether or not there is a bar or an AGN.

The mass-metallicity relation shows that the galaxies with the lowest mass have the lowest metallicities. As most dwarf galaxies are in group environments, interaction effects such as tides could contribute to this trend. We perform a series of smoothed particle hydrodynamics (SPH) simulations of dwarf galaxies in external tidal fields to examine the effects of tides on their metallicities and metallicity gradients. In our simulated galaxies, gravitational instabilities drive gas inwards and produce centralized star formation and a significant metallicity gradient. Strong tides can contribute to these instabilities, but their primary effect is to strip the outer low-metallicity gas, producing a truncated gas disk with a large metallicity. This suggests that the role of tides on the mass-metallicity relation is to move dwarf galaxies to higher metallicities.

We perform a series of smoothed particle hydrodynamics simulations of isolated dwarf galaxies to compare different metal mixing models. In particular, we examine the role of diffusion in the production of enriched outflows, and in determining the metallicity distributions of gas and stars. We investigate different diffusion strengths, by changing the pre-factor of the diffusion coefficient, by varying how the diffusion coefficient is calculated from the local velocity distribution, and by varying whether the speed of sound is included as a velocity term. Stronger diffusion produces a tighter [O/Fe]-[Fe/H] distribution in the gas, and cuts off the gas metallicity distribution function at lower metallicities. Diffusion suppresses the formation of low-metallicity stars, even with weak diffusion, and also strips metals from enriched outflows. This produces a remarkably tight correlation between "metal mass-loading" (mean metal outflow rate divided by mean metal production rate) and the strength of diffusion, even when the diffusion coefficient is calculated in different ways. The effectiveness of outflows at removing metals from dwarf galaxies and the metal distribution of the gas is thus dependent on the strength of diffusion. By contrast, we show that the metallicities of stars are not strongly dependent on the strength of diffusion, provided that some diffusion is present.

We examine simulations of isolated galaxies to analyse the effects of localised feedback on the formation and evolution of molecular clouds. Feedback contributes to turbulence and the destruction of clouds, leading to a population of clouds that is younger, less massive, and with more retrograde rotation. We investigate the evolution of clouds as they interact with each other and the diffuse ISM, and determine that the role of cloud interactions differs strongly with the presence of feedback: in models without feedback, scattering events dramatically increase the retrograde fraction, but in models with feedback, mergers between clouds may slightly increase the prograde fraction. We also produce an estimate of the viscous time-scale due to cloud-cloud collisions, which increases with increasing strength of feedback (~20 Gyr vs ~10 Gyr), but is still much smaller than previous estimates (~1000 Gyr); although collisions become more frequent with feedback, less energy is lost in each collision than in the models without feedback.

Observations have revealed cold gas with large velocity dispersions (~300 km/s) within the hot outflows of ultra-luminous infrared galaxies (ULIRGs). This gas may trace its origin to the Rayleigh-Taylor (RT) fragmentation of a super-bubble or may arise on smaller scales. We model a ULIRG outflow at two scales to recreate this gas in three-dimensional hydrodynamic simulations using FLASH. Although resolution is limited, these models successfully produce cold gas in outflows with large velocity dispersions. Our small-scale models produce this cold gas through RT fragmentation of the super-bubble wall, but the large-scale models produce the cold gas after hot bubbles fragment the disc's gas into cold clouds which are then accelerated by thermal pressure, and supplemented by cooling within the outflow. We produce simple mock spectra to compare these simulations to observed absorption spectra and find line-widths of ~250 km/s, agreeing with the lower end of observations.

Analytic estimates of the viscous time-scale due to cloud-cloud collisions have been as high as thousands of Gyr. Consequently, cloud collisions are widely ignored as a source of viscosity in galactic disks. However, capturing the hydrodynamics of discs in simple analytic models is a challenge, both because of the wide dynamic range and importance of 2D and 3D effects. To test the validity of analytic models we present estimates for the viscous time-scale that are measured from three dimensional SPH simulations of disc formation and evolution. We have deliberately removed uncertainties associated with star-formation and feedback thereby enabling us to place lower bounds on the time-scale for this process. We also contrast collapse simulations with results from simulations of initially stable discs and examine the impact of numerical parameters and assumptions on our work, to constrain possible systematics in our estimates. We find that cloud-collision viscous time-scales are in the range of 0.6-16 Gyr, considerably shorter than previously estimated. This large discrepency can be understood in terms of how the efficiency of collisions is included in the analytical estimates. We find that the viscous time-scale only depends weakly on the number of clouds formed, and so while the viscous time-scale will increase with increasing resolution, this effect is too weak to alter our conclusions.