No Arabic abstract
We perform controlled N-Body/SPH simulations of disk galaxy formation by cooling a rotating gaseous mass distribution inside equilibrium cuspy spherical and triaxial dark matter halos. We systematically study the angular momentum transport and the disk morphology as we increase the number of dark matter and gas particles from 10^4 to 10^6, and decrease the gravitational softening from 2 kpc to 50 parsecs. The angular momentum transport, disk morphology and radial profiles depend sensitively on force and mass resolution. At low resolution, similar to that used in most current cosmological simulations, the cold gas component has lost half of its initial angular momentum via different mechanisms. The angular momentum is transferred primarily to the hot halo component, by resolution-dependent hydrodynamical and gravitational torques, the latter arising from asymmetries in the mass distribution. In addition, disk-particles can lose angular momentum while they are still in the hot phase by artificial viscosity. In the central disk, particles can transfer away over 99% of their initial angular momentum due to spiral structure and/or the presence of a central bar. The strength of this transport also depends on force and mass resolution - large softening will suppress the bar instability, low mass resolution enhances the spiral structure. This complex interplay between resolution and angular momentum transfer highlights the complexity of simulations of galaxy formation even in isolated haloes. With 10^6 gas and dark matter particles, disk particles lose only 10-20% of their original angular momentum, yet we are unable to produce pure exponential profiles.
We present the McMaster Unbiased Galaxy Simulations (MUGS), the first 9 galaxies of an unbiased selection ranging in total mass from 5$times10^{11}$ M$_odot$ to 2$times10^{12}$ M$_odot$ simulated using n-body smoothed particle hydrodynamics (SPH) at high resolution. The simulations include a treatment of low temperature metal cooling, UV background radiation, star formation, and physically motivated stellar feedback. Mock images of the simulations show that the simulations lie within the observed range of relations such as that between color and magnitude and that between brightness and circular velocity (Tully-Fisher). The greatest discrepancy between the simulated galaxies and observed galaxies is the high concentration of material at the center of the galaxies as represented by the centrally peaked rotation curves and the high bulge-to-total ratios of the simulations determined both kinematically and photometrically. This central concentration represents the excess of low angular momentum material that long has plagued morphological studies of simulated galaxies and suggests that higher resolutions and a more accurate description of feedback will be required to simulate more realistic galaxies. Even with the excess central mass concentrations, the simulations suggest the important role merger history and halo spin play in the formation of disks.
We study the z=0 gas kinematics, morphology, and angular momentum content of isolated galaxies in a suite of cosmological zoom-in simulations from the FIRE project spanning $M_{star}=10^{6-11}M_{odot}$. Gas becomes increasingly rotationally supported with increasing galaxy mass. In the lowest-mass galaxies ($M_{star}<10^{8}M_{odot}$), gas fails to form a morphological disk and is primarily dispersion and pressure supported. At intermediate masses ($M_{star}=10^{8-10}M_{odot}$), galaxies display a wide range of gas kinematics and morphologies, from thin, rotating disks, to irregular spheroids with negligible net rotation. All the high-mass ($M_{star}=10^{10-11}M_{odot}$) galaxies form rotationally supported gas disks. Many of the halos whose galaxies fail to form disks harbor high angular momentum gas in their circumgalactic medium. The ratio of the specific angular momentum of gas in the central galaxy to that of the dark-matter halo increases significantly with galaxy mass, from $j_{rm gas}/j_{rm DM}sim0.1$ at $M_{star}=10^{6-7}M_{odot}$ to $j_{rm gas}/j_{rm DM}sim2$ at $M_{star}=10^{10-11}M_{odot}$. The reduced rotational support in the lowest-mass galaxies owes to (a) stellar feedback and the UV background suppressing the accretion of high-angular momentum gas at late times, and (b) stellar feedback driving large non-circular gas motions. We broadly reproduce the observed scaling relations between galaxy mass, gas rotation velocity, size, and angular momentum, but may somewhat underpredict the incidence of disky, high-angular momentum galaxies at the lowest observed masses ($M_{star}=(10^{6}-2times10^{7})M_{odot}$). In our simulations, stars are uniformly less rotationally supported than gas. The common assumption that stars follow the same rotation curve as gas thus substantially overestimates galaxies stellar angular momentum, particularly at low masses.
We study the spatially-resolved stellar specific angular momentum $j_*$ in a high-quality sample of 24 CALIFA galaxies covering a broad range of visual morphology, accounting for stellar velocity and velocity dispersion. The shape of the spaxel-wise probability density function of normalised $s=j_*/j_{*mean}$, PDF($s$), deviates significantly from the near-universal initial distribution expected of baryons in a dark matter halo and can be explained by the expected baryonic effects in galaxy formation that remove and redistribute angular momentum. Further we find that the observed shape of the PDF($s$) correlates significantly with photometric morphology, where late-type galaxies have PDF($s$) that is similar to a normal distribution, whereas early types have a strongly-skewed PDF($s$) resulting from an excess of low-angular momentum material. Galaxies that are known to host pseudobulges (bulge Sersic index $n_b <2.2$) tend to have less skewed bulge PDF($s$), with skewness $(b_{1rb})lesssim0.8$. The PDF($s$) encodes both kinematic and photometric information and appears to be a robust tracer of morphology. Its use is motivated by the desire to move away from traditional component-based classifications which are subject to observer bias, to classification on a galaxys fundamental (stellar mass, angular momentum) properties. In future, PDF($s$) may also be useful as a kinematic decomposition tool.
We perform N-Body/SPH simulations of disk galaxy formation inside equilibrium spherical and triaxial cuspy dark matter halos. We systematically study the disk properties and morphology as we increase the numbers of dark matter and gas particles from 10^4 to 10^6 and change the force resolution. The force resolution influences the morphological evolution of the disk quite dramatically. Unless the baryon fraction is significantly lower than the universal value, with high force resolution a gaseous bar always forms within a billion years after allowing cooling to begin. The bar interacts with the disk, transferring angular momentum and increasing its scale length. In none of the simulations does the final mass distribution of the baryons obey a single exponential profile. Indeed within a few hundred parsecs to a kiloparsec from the center the density rises much more steeply than in the rest of the disk, and this is true irrespective of the presence of the bar.
Feedback from energy liberated by gas accretion onto black holes (BHs) is an attractive mechanism to explain the exponential cut-off at the massive end of the galaxy stellar mass function (SMF). Semi-analytic models of galaxy formation in which this form of feedback is assumed to suppress cooling in haloes where the gas cooling time is large compared to the dynamical time do indeed achieve a good match to the observed SMF. Furthermore, hydrodynamic simulations of individual halos in which gas is assumed to accrete onto the central BH at the Bondi rate have shown that a self-regulating regime is established in which the BH grows just enough to liberate an amount of energy comparable to the thermal energy of the halo. However, this process is efficient at suppressing the growth not only of massive galaxies but also of galaxies like the Milky Way, leading to disagreement with the observed SMF. The Bondi accretion rate, however, is inappropriate when the accreting material has angular momentum. We present an improved accretion model that takes into account the circularisation and subsequent viscous transport of infalling material and include it as a subgrid model in hydrodynamic simulations of the evolution of halos with a wide range of masses. The resulting accretion rates are generally low in low mass ($lsim 10^{11.5} msun$) halos, but show outbursts of Eddington-limited accretion during galaxy mergers. During outbursts these objects strongly resemble quasars. In higher mass haloes, gas accretion occurs continuously, typically at $~10$ % of the Eddington rate, which is conducive to the formation of radio jets. The resulting dependence of the accretion behaviour on halo mass induces a break in the relation between galaxy stellar mass and halo mass in these simulations that matches observations.