Remarkable progress has been made over the last few decades in furthering our understanding of the growth of cosmic structure. Nonetheless, there remains a great deal of uncertainty regarding the precise details of the complex baryonic physics that regulate galaxy formation. Any theory of star formation in galaxies must encompass the radiative cooling of gas into dark matter haloes, the formation of a turbulent, multiphase interstellar medium (ISM), the efficiency with which molecular gas is able to collapse into cores and ultimately stars, and the subsequent interaction of those stars with the gas through ionizing radiation, winds and supernova (SN) explosions. Given the highly non-linear nature of the problem, numerical simulations provide an invaluable tool with which to study galaxy formation. Yet, even with contemporary computational resources, the inherently large dynamical range of spatial scales that must be tackled makes the development of such models extremely challenging, inevitably leading to the adoption of `subgrid' approximations at some scale. In this thesis, I explore new methods of incorporating the physics of star formation and stellar feedback into high resolution hydrodynamic simulations of galaxies. I first describe a new implementation of star formation and SN feedback that I have developed for the state-of-the-art moving mesh code Arepo. I carry out a detailed study into various classes of subgrid SN feedback schemes commonly adopted in the literature, including injections of thermal and/or kinetic energy, two parametrizations of delayed cooling feedback and a 'mechanical' feedback scheme that injects the appropriate amount of momentum depending on the relevant scale of the SN remnant (SNR) resolved. All schemes make use of individually time-resolved SN events. Adopting isolated disk galaxy setups at different resolutions, with the highest resolution runs reasonably resolving the Sedov-Taylor phase of the SNR, I demonstrate that the mechanical scheme is the only physically well-posed method of those examined, is efficient at suppressing star formation, agrees well with the Kennicutt-Schmidt relation and leads to converged star formation rates and galaxy morphologies with increasing resolution without fine tuning any parameters. However, I find that it is difficult to produce outflows with high enough mass loading factors at all but the highest resolution. I discuss the various possible solutions to this effect, including improved modelling of star formation. Moving on to a more self-consistent setup, I carry out a suite of cosmological zoom-in simulations of low mass haloes at very high resolution, performed to z = 4, to investigate the ability of SN feedback models to produce realistic galaxies. The haloes are selected in a variety of environments, ranging from voids to crowded locations. In the majority of cases, SN feedback alone has little impact at early times even in low mass haloes ($\sim10^{10}\,\mathrm{M_\odot}$ at z = 0). This appears to be due largely to the build up of very dense gas prior to SN events, suggesting that other mechanisms (such as other stellar feedback processes) are required to regulate ISM properties before SNe occur. The effectiveness of the feedback also appears to be strongly dependent on the merger history of the halo. Finally, I describe a new scheme to drive turbulence in isolated galaxy setups. The turbulent structure of the ISM very likely regulates star formation efficiencies on small scales, as well as affecting the clustering of SNe. The large range of potential drivers of ISM turbulence are not fully understood and are, in any case, unlikely to arise ab initio in a whole galaxy simulation. I therefore neglect these details and adopt a highly idealised approach, artificially driving turbulence to produce an ISM structure of my choice. This enables me to study the effects of a given level of ISM turbulence on global galaxy properties, such as the fragmentation scale of the disk and the impact on SN feedback efficiencies. I demonstrate this technique in the context of simulations of isolated dwarfs, finding that moderate levels of turbulent driving in combination with SN feedback can produce a steady-state of star formation rates and global galaxy properties, rather than the extremely violent SN feedback that is produced by a rapidly fragmenting disk.

Aiming at understanding the role of stellar feedback in galaxy evolution, I present a study of the hot interstellar medium in several representative galaxies, based primarily on X-ray observations as well as theoretical modelling. I find that, in the massive disk galaxies NGC2613 and M104, the observed amount of hot gas is much less than that predicted by current galaxy formation models. Such a discrepancy suggests a lack of appropriate treatments of stellar/AGN feedback in these models. I also find that stellar feedback, primarily in the form of mass loss from evolved stars and energy released from supernovae, and presumably consumed by the hot gas, is largely absent from the inner regions of M104, a galaxy of a substantial content of evolved stars but little current star formation. A natural understanding of this phenomenon is that the hot gas is in the form of a galactic-scale outflow, by which the bulk of the stellar feedback is transported to the outer regions and perhaps into the intergalactic space. A comparison between the observed sub-galactic gas structures and model predictions indicate that this outflow is probably subsonic rather than being a classical supersonic galactic wind. Such outflows are likely prevalent in most early-type galaxies of intermediate masses in the present-day universe and thus play a crucial role in the evolution of such galaxies. For the first time I identify the presence of diffuse hot gas in and around the bulge of the Andromeda Galaxy (M31), our well-known neighbor. Both the morphology and energetics of the hot gas suggest that it is also in the form of a large-scale outflow. Assisted with multiwavelength observations toward the circumnuclear regions of M31, I further reveal the relation between the hot gas and other cooler phases of the interstellar medium. I suggest that thermal evaporation, mostly likely energized by Type Ia supernovae, acts to continuously turn cold gas into hot, a process that naturally leads to the inactivity of the central supermassive blackhole as well as the launch of the hot gas outflow. Such a mechanism plays an important role in regulating the multi-phase interstellar medium in the circumnuclear environment and transporting stellar feedback to the outer galactic regions.

<p>The effect of stellar feedback on the interstellar medium is investigated using numerical simulation. In particular, the roles of supernova feedback and ionization feedback on the star formation rate and structure of the interstellar medium are compared. We use Enzo, an adaptive mesh code, and employ the MUSCL-Hancock hydrodynamics scheme to run simulations of a section of a stratified galactic disk. A turbulent velocity field is imposed in the central region of the disk and self-gravity is applied. Star clusters are formed when density and temperature conditions are met, which, in turn, provide ionization and supernova feedback into the interstellar medium. Simulations were run with and without supernova and ionization feedback and the runs are compared. Ionization feedback is found to dominate over supernova feedback in regulating star formation rates. With no feedback, all the gas is converted to stars by 200 Myr. When supernova feedback is added, 98% of the gas is used to create stars by 300 Myr. With ionization feedback instead, at 1 Gyr into the run, only 30% of the gas is in stars. Even with supernova feedback added to ionization feedback, the gas converted to stars is just 29% at 1 Gyr. Very strong supernovae take this fraction down to 25%. The star formation rates in the runs with supernova feedback are consistent with the low end of the Kennicutt-Schmidt relation, while the runs without ionization feedback have star formation rates that are an order of magnitude larger. Gas phase masses and volumes produced in the ionization runs are broadly consistent with observations.</p> / Master of Science (MSc)

In this thesis we test and expand our current knowledge of Star Formation Laws (SF laws) in the extreme environment of dwarf irregular galaxies. We focus on the SF characteristics of our 18 galaxies sample, extending current investigations of the Schmidt-Kennicutt law to the low luminosity, low metallicity regime. The Hi data used in this project have been observed, calibrated and imaged according to the LITTLE THINGS Survey prescription to which I brought my own contribution as a member of the team. Apart from high resolution, VLA data in B, C and D array configurations, this project makes use of an extensive set of multi- wavelength data (H , FUV, 24 m, 3.6 m, V-band and K-band). Molecular gas in dwarfs is very difficult to observe, mainly because due to the low metallicity environment, we lose our only molecular tracer, the CO which becomes under luminous. Therefore the gas distribution is represented by Hi gas only. We create our Star Formation Rate (SFR) maps mainly based on FUV maps because our analysis shows that FUV is the SF tracer that allows us the most extensive sampling of the SFR surface density (SFRD) and Hi surface density relation. The main results of our study are: Whereas in spiral galaxies Bigiel et al. (2008) have found a one to one relation between star formation rate and molecular gas and no relation between the SFR and the neutral gas, in a small sample of dwarfs as well as in the outskirts of spiral galaxies Bigiel et al. (2010b) has found that SFRD does correlate with Hi surface density. We confirm the existence of the SFRD vs. Hi surface density relation in dwarf irregular galaxies and a linear fitting through all our data (all 18 galaxies combined) yields a power law relation ΣSFR ∝ Σ1.87±0.3/HI . We find that the interiors of Hi shells, at 400 pc scales, become resolved and show up in SFRD versus Hi surface density plots although within the shell interior we have SFRD values but no Hi surface density related to them. Thus, the points originating from those regions contribute significantly to the increase of the scatter in the plot. We show that by excluding those points the correlation between SFRD and Hi surface density improves between 10% and 20%. Eight of the 18 galaxies in our sample have Hi maxima higher than the 10M pc-2 value found by Bigiel et al. (2008) for spiral galaxies. Krumholz et al. (2011) predicted that the 10M pc-2 threshold is metallicity dependent in galaxies with sub-solar metallicity, however the theoretically predicted values for our galaxies only match the observed Hi maxima in one case (DDO168). We find that metallicity cannot be the only factor setting the Hi to H2 transition. In fact, we find evidence that the higher the interstellar radiation field (ISRF), the higher the Hi maximum is, hence we suggest that the ISRF should also be taken into consideration in predicting the Hi to H2 transition threshold. We find that even tighter than the SFRD vs. Hi surface density relation is the SFRD vs. V-band surface density relation. Unlike the SFRD vs. Hi surface density relation the SFRD vs. V-band surface density relation follows a power law and can be written as follows: ΣSFR ∝ (10^μv)^-0.43±0.03. The SFRD vs. V-band surface density relation suggests that the existing stars also play a role in the formation of the next generation of stars. Within our sample of dwarf galaxies the average pressure per resolution element and the SFRD are in a 1:1 linear relation: ΣSFR ∝ P_h^1.02±0.05. A similar relation has been found by Blitz & Rosolowsky (2006) for the low-pressure regimes of spiral galaxies. In conclusion we find that in the extreme environments of dwarf galaxies the metal deficiency and the lack of the classic SF stimulators (spiral arms, shear motions) do not impede the star forming process. In these galaxies, dust-shielding becomes predominantly self-shielding and there is plenty of Hi available to achieve this additional task. Existing stars assume the role of pressure enhancers, which in turn will stimulate SF without the need of spiral arms or shear motion.

No description available.