Supernovae (SN) generate hot gas in the interstellar medium (ISM), help setting the ISM structure and support the driving of outflows. It is important to resolve the hot gas generation for galaxy formation simulations at solar mass and sub-parsec resolution which realise individual supernova (SN) explosions with ambient densities varying by several orders of magnitude in a realistic multi-phase ISM. We test resolution requirements by simulating SN blast waves at three metallicities ($Z = 0.01, 0.1$ and $1 Z_{\odot}$), six densities and their respective equilibrium chemical compositions ($n=0.001$ cm$^{-3}$ - $100$ cm$^{-3}$), and four mass resolutions ($0.1$ - $100$ M$_{\odot}$), in three dimensions. We include non-equilibrium cooling and chemistry, a homogenous interstellar radiation field, and shielding with a modern pressure-energy smoothed particle hydrodynamics (SPH) method including isotropic thermal conduction and a meshless-finite-mass (MFM) solver. We find stronger resolution requirements for chemistry and hot phase generation than for momentum generation. While at $10$ M$_{\odot}$ the radial momenta at the end of the Sedov phase start converging, the hot phase generation and chemistry require higher resolutions to represent the neutral to ionised hydrogen fraction at the end of the Sedov phase correctly. Thermal conduction typically reduces the hot phase by $0.2$ dex and has little impact on the chemical composition. In general, our $1$, and $0.1$ M$_{\odot}$ results agree well with previous numerical and analytic estimates. We conclude that for the thermal energy injection SN model presented here resolutions higher than $10$ M$_{\odot}$ are required to model the chemistry, momentum and hot phase generation in a multi-phase ISM.

We describe a population of young star clusters (SCs) formed in a hydrodynamical simulation of a gas-rich dwarf galaxy merger resolved with individual massive stars at sub-parsec spatial resolution. The simulation is part of the GRIFFIN (Galaxy Realizations Including Feedback From INdividual massive stars) project. The star formation environment during the simulation spans seven orders of magnitude in gas surface density and thermal pressure, and the global star formation rate surface density ($\Sigma_\mathrm{SFR}$) varies by more than three orders of magnitude during the simulation. Young SCs more massive than $M_{\mathrm{*,cl}}\sim 10^{2.5}\,M_{\odot}$ form along a mass function with a power-law index $\alpha\sim-1.7$ ($\alpha\sim-2$ for $M_{\mathrm{*,cl}}\gtrsim10^{3}\,M_{\odot}$) at all merger phases, while the normalization and the highest SC masses (up to $\sim 10^6 M_{\odot}$) correlate with $\Sigma_\mathrm{SFR}$. The cluster formation efficiency varies from $\Gamma\sim20\%$ in early merger phases to $\Gamma\sim80\%$ at the peak of the starburst and is compared to observations and model predictions. The massive SCs ($\gtrsim10^4\,M_{\odot}$) have sizes and mean surface densities similar to observed young massive SCs. Simulated lower mass clusters appear slightly more concentrated than observed. All SCs form on timescales of a few Myr and lose their gas rapidly resulting in typical stellar age spreads between $\sigma\sim0.1-2$ Myr ($1\sigma$), consistent with observations. The age spreads increase with cluster mass, with the most massive cluster ($\sim10^6\, M_{\odot}$) reaching a spread of $5\, \mathrm{Myr}$ once its hierarchical formation finishes. Our study shows that it is now feasible to investigate the SC population of entire galaxies with novel high-resolution numerical simulations.

We present a novel and fast application to generate glass-like initial conditions for Lagrangian hydrodynamic schemes (e.g. Smoothed Particle Hydrodynamics (SPH)) following arbitrary density models based on weighted Voronoi tessellations and combine it with improved initial configurations and an additional particle reshuffling scheme. We show our application's ability to sample different kinds of density features and to converge properly towards the given model density as well as a glass-like particle configuration. We analyse convergence with iterations as well as with varying particle number. Additionally, we demonstrate the versatility of the implemented algorithms by providing an extensive test suite for standard (magneto-) hydrodynamic test cases as well as a few common astrophysical applications. We indicate the potential to bridge further between observational astronomy and simulations as well as applicability to other fields of science by advanced features such as describing a density model using gridded data for exampling from an image file instead of an analytic model.

We investigate the build-up of the galactic dynamo and subsequently the origin of a magnetic driven outflow. We use a setup of an isolated disc galaxy with a realistic circum-galactic medium (CGM). We find good agreement of the galactic dynamo with theoretical and observational predictions from the radial and toroidal components of the magnetic field as function of radius and disc scale height. We find several field reversals indicating dipole structure at early times and quadrupole structure at late times. Together with the magnetic pitch angle and the dynamo control parameters $R_{\alpha}$, $R_{\omega}$ and $D$ we present strong evidence for an $\alpha^2$-$\Omega$ dynamo. The formation of a bar in the centre leads to further amplification of the magnetic field via adiabatic compression which subsequently drives an outflow. Due to the Parker-Instability the magnetic field lines rise to the edge of the disc, break out and expand freely in the CGM driven by the magnetic pressure. Finally, we investigate the correlation between magnetic field and star formation rate. Globally, we find that the magnetic field is increasing as function of the star formation rate surface density with a slope between $0.3$ and $0.45$ in good agreement with predictions from theory and observations. Locally, we find that the magnetic field can decrease while star formation increases. We find that this effect is correlated with the diffusion of magnetic field from the spiral arms to the inter-arm regions which we explicitly include by solving the induction equation and accounting for non-linear terms.

Galaxy mergers are a fundamental part of galaxy evolution. To study the resulting mass distributions of different kinds of galaxy mergers, we present a simulation suite of 36 high-resolution isolated merger simulations, exploring a wide range of parameter space in terms of mass ratios (mu = 1:5, 1:10, 1:50, 1:100) and orbital parameters. We find that mini mergers deposit a higher fraction of their mass in the outer halo compared to minor mergers, while their contribution to the central mass distribution is highly dependent on the orbital impact parameter: for larger pericentric distances we find that the centre of the host galaxy is almost not contaminated by merger particles. We also find that the median of the resulting radial mass distribution for mini mergers differs significantly from the predictions of simple theoretical tidal-force models. Furthermore, we find that mini mergers can increase the size of the host disc significantly without changing the global shape of the galaxy, if the impact occurs in the disc plane, thus providing a possible explanation for extended low-surface brightness disks reported in observations. Finally, we find clear evidence that streams are a strong indication of nearly circular infall of a satellite (with large angular momentum), whereas the appearance of shells clearly points to (nearly) radial satellite infall.

We present simulations of isolated disc galaxies in a realistic environment performed with the Tree-SPMHD-Code Gadget-3. Our simulations include a spherical circum-galactic medium (CGM) surrounding the galactic disc, motivated by observations and the results of cosmological simulations. We present three galactic models with different halo masses between 10e10 Msol and 10e12 Msol, and for each we use two different approaches to seed the magnetic field, as well as a control simulation without a magnetic field. We find that the amplification of the magnetic field in the centre of the disc leads to a biconical magnetic outflow of gas that magnetizes the CGM. This biconical magnetic outflow reduces the star formation rate (SFR) of the galaxy by roughly 40 percent compared to the simulations without magnetic fields. As the key aspect of our simulations, we find that small scale turbulent motion of the gas in the disc leads to the amplification of the magnetic field up to tens of 10e-6 G, as long as the magnetic field strength is low. For stronger magnetic fields turbulent motion does not lead to significant amplification but is replaced by an alpha-omega dynamo. The occurance of a small scale turbulent dynamo becomes apparent through the magnetic power spectrum and analysis of the field lines' curvature. In accordance with recent observations we find an anti-correlation between the spiral structure in the gas density and in the magnetic field due to a diffusion term added to the induction equation.