Modeling Physical Processes at Galactic Scales and Above

What should these lectures be? The subject assigned to us is so broad that many books can be written about it. So, in planning these lectures I had several options. One would be to focus on a narrow subset of topics and to cover them in great detail. Such a subset necessarily would be highly personal and useful to a few read- ers at best. Another option would be to give a very shallow overview of the whole field, but then it wont be very much different from a highly compressed version of a university course (which anyone can take if they wish so). So, I decided to be selfish and to prepare these lectures as if I was teaching my own graduate student. Given my research interests, I selected what the student would need to know to be able to discuss science with me and to work on joint research projects. So, the story presented below is both personal and incomplete, but it does cover several subjects that are poorly represented in the existing textbooks (if at all). Some of topics I focus on below are closely connected, others are disjoint, some are just side detours on specific technical questions. There is an overlapping theme, however. Our goal is to follow the cosmic gas from large scales, low densities, (rel- atively) simple physics to progressively smaller scales, higher densities, closer rela- tion to galaxies, and more complex and uncertain physics. So, we (you - the reader, and me - the author) are going to follow a yellow brick road from the gas well be- yond any galaxy confines to the actual sites of star formation and stellar feedback. On the way we will stop at some places for a tour and run without looking back through some others. So, the road will be uneven, but I hope that some readers find it useful.

On the mass of the Local Group

We use recent proper motion measurements of the tangential velocity of M31, along with its radial velocity and distance, to derive the likelihood of the sum of halo masses of the Milky Way and M31. This is done using a sample halo pairs in the Bolshoi cosmological simulation of $Lambda$CDM cosmology selected to match properties and environment of the Local Group. The resulting likelihood gives estimate of the sum of masses of $M_{rm MW,200}+M_{rm M31,200}=$ $2.40_{-1.05}^{+1.95}times10^{12},M_{odot}$ ($90%$ confidence interval). This estimate is consistent with individual mass estimates for the Milky Way and M31 and is consistent, albeit somewhat on the low side, with the mass estimated using the timing argument. We show that although the timing argument is unbiased on average for all pairs, for pairs constrained to have radial and tangential velocities similar to that of the Local Group the argument overestimates the sum of masses by a factor of $1.6$. Using similar technique we estimate the total dark matter mass enclosed within $1$ Mpc from the Local Group barycenter to be $M_{rm LG}(r<1, {rm Mpc})=4.2_{-2.0}^{+3.4}times10^{12},M_{odot}$ ($90%$ confidence interval).

Satellites in MW-like hosts: Environment dependence and close pairs

Previous studies showed that an estimate of the likelihood distribution of the Milky Way halo mass can be derived using the properties of the satellites similar to the Large and Small Magellanic Clouds (LMC and SMC). However, it would be straightforward to interpret such an estimate only if the properties of the Magellanic Clouds (MCs) are fairly typical and are not biased by the environment. In this study we explore whether the environment of the Milky Way affects the properties of the SMC and LMC such as their velocities. To test for the effect of the environment, we compare velocity distributions for MC-sized subhalos around Milky Way hosts in a sample selected simply by mass and in the second sample of such halos selected with additional restrictions on the distance to the nearest cluster and the local galaxy density, designed to mimic the environment of the Local Group (LG). We find that satellites in halos in the LG-like environments do have somewhat larger velocities, as compared to the halos of similar mass in the sample without environmental constraints. We derive the host halo likelihood distribution for the samples in the LG-like envirionment and in the control sample and find that the environment does not significantly affect the derived likelihood. We use the updated properties of the SMC and LMC to derive the constraint on the MW halo mass $log{({rm M}_{200} /msol)}=12.06^{+0.31}_{-0.19}$ (90% confidence interval). We also explore the incidence of close pairs with relative velocities and separations similar to those of the LMC and SMC and find that such pairs are quite rare among $Lambda$CDM halos. Taking into account the close separation of the MCs in the Busha et al. 2011 method results in the shift of the MW halo mass estimate to smaller masses, with the peak shifting approximately by a factor of two.[Abridged]

Ultra-faint dwarf galaxies as a test of early enrichment and metallicity-dependent star formation

The tight relation of star formation with molecular gas indicated by observations and assumed in recent models implies that the efficiency with which galaxies convert their gas into stars depends on gas metallicity. This is because the abundance of molecular hydrogen is sensitive to the abundance of dust, which catalyzes the formation of H_2 and helps to shield it from dissociating radiation. In this study we point out that in the absence of significant pre-enrichment by Population III stars forming out of zero metallicity gas, such H_2-based star formation is expected to leave an imprint in the form of bi-modality in the metallicity distribution among dwarf galaxies and in the metallicity distribution of stars within individual galaxies. The bi-modality arises because when gas metallicity (and dust abundance) is low, formation of molecular gas is inefficient, the gas consumption time scale is long, and star formation and metal enrichment proceed slowly. When metallicity reaches a critical threshold value star formation and enrichment accelerate, which leads to rapid increase in both stellar mass and metallicity of galaxies. We demonstrate this process both using a simple analytical model and full cosmological simulations. In contrast, observed metallicity distributions of dwarf galaxies or stars within them are not bi-modal. We argue that this discrepancy points to substantial early stochastic pre-enrichment by population III stars to levels Z ~ 0.01 Z_sun in dense, star forming regions of early galaxies.

Resolving Gas Dynamics in the Circumnuclear Region of a Disk Galaxy in a Cosmological Simulation

Using a hydrodynamic adaptive mesh refinement code, we simulate the growth and evolution of a galaxy, which could potentially host a supermassive black hole, within a cosmological volume. Reaching a dynamical range in excess of 10 million, the simulation follows the evolution of the gas structure from super-galactic scales all the way down to the outer edge of the accretion disk. Here, we focus on global instabilities in the self-gravitating, cold, turbulence-supported, molecular gas disk at the center of the model galaxy, which provide a natural mechanism for angular momentum transport down to sub-pc scales. The gas density profile follows a power-law scaling as r^-8/3, consistent with an analytic description of turbulence in a quasi-stationary circumnuclear disk. We analyze the properties of the disk which contribute to the instabilities, and investigate the significance of instability for the galaxys evolution and the growth of a supermassive black hole at the center.