The evolution and distribution of the angular momentum of dark matter (DM) halos have been discussed in several studies over the past decades. In particular, the idea arose that angular momentum conservation should allow us to infer the total angular momentum of the entire DM halo from measuring the angular momentum of the baryonic component, which is populating the center of the halo, especially for disk galaxies. To test this idea and to understand the connection between the angular momentum of the DM halo and its galaxy, we use a state-of-the-art, hydrodynamical cosmological simulation taken from the set of Magneticum Pathfinder simulations. Thanks to the inclusion of the relevant physical processes, the improved underlying numerical methods, and high spatial resolution, we successfully produce populations of spheroidal and disk galaxies self-consistently. Thus, we are able to study the dependence of galactic properties on their morphology. We find that(1) the specific angular momentum of stars in disk and spheroidal galaxies as a function of their stellar mass compares well with observational results; (2) the specific angular momentum of the stars in disk galaxies is slightly smaller compared to the specific angular momentum of the cold gas, in good agreement with observations; (3) simulations including the baryonic component show a dichotomy in the specific stellar angular momentum distribution when splitting the galaxies according to their morphological type(this dichotomy can also be seen in the spin parameter, where disk galaxies populate halos with slightly larger spin compared to spheroidal galaxies); (4) disk galaxies preferentially populate halos in which the angular momentum vector of the DM component in the central part shows a better alignment to the angular momentum vector of the entire halo; and (5) the specific angular momentum of the cold gas in disk galaxies is approximately 40% smaller than the specific angular momentum of the total DM halo and shows a significant scatter.

We present an implementation of smoothed particle hydrodynamics (SPH) with improved accuracy for simulations of galaxies and the large-scale structure. In particular, we implement and test a vast majority of SPH improvement in the developer version of GADGET-3. We use the Wendland kernel functions, a particle wake-up time-step limiting mechanism and a time-dependent scheme for artificial viscosity including highorder gradient computation and shear flow limiter. Additionally, we include a novel prescription for time-dependent artificial conduction, which corrects for gravitationally induced pressure gradients and improves the SPH performance in capturing the development of gas-dynamical instabilities.We extensively test our new implementation in a wide range of hydrodynamical standard tests including weak and strong shocks as well as shear flows, turbulent spectra, gas mixing, hydrostatic equilibria and self-gravitating gas clouds. We jointly employ all modifications; however, when necessary we study the performance of individual code modules. We approximate hydrodynamical states more accurately and with significantly less noise than standard GADGET-SPH. Furthermore, the new implementation promotes the mixing of entropy between different fluid phases, also within cosmological simulations.Finally, we study the performance of the hydrodynamical solver in the context of radiative galaxy formation and non-radiative galaxy cluster formation. We find galactic disks to be colder and more extended and galaxy clusters showing entropy cores instead of steadily declining entropy profiles. In summary, we demonstrate that our improved SPH implementation overcomes most of the undesirable limitations of standard GADGET-SPH, thus becoming the core of an efficient code for large cosmological simulations.

In large scale cosmological hydrodynamic simulations simplified sub-grid models for gas accretion onto black holes and AGN feedback are commonly used. Such models typically depend on various free parameters, which are not well constrained. We present a new advanced model containing a more detailed description of AGN feedback, where those parameters reflect the results of recent observations. The model takes the dependency of these parameters on the black hole properties into account and describes a continuous transition between the feedback processes acting in the so-called radio-mode and quasar-mode. In addition, we implement a more detailed description of the accretion of gas onto black holes by distinguishing between hot and cold gas accretion. Our new implementations prevent black holes from gaining too much mass, particularly at low redshifts, so that our simulations are successful in reproducing the observed present-day black hole mass function. Our new model also suppresses star formation in massive galaxies slightly more efficiently than many state-of-the-art models. Therefore, the simulations that include our new implementations produce a more realistic population of quiescent and star-forming galaxies compared to recent observations, even if some discrepancies remain. In addition, the baryon conversion efficiencies in our simulation are -except for the high mass endconsistent with observations presented in literature over the mass range resolved by our simulations. Finally, we discuss the significant impact of the feedback model on the low-luminous end of the AGN luminosity function.

In the last few years, it became possible to observationally resolve galaxies with two distinct nuclei in their centre. For separations smaller than 10kpc, dual and offset active galactic nuclei (AGN) are distinguished: in dual AGN, both nuclei are active, whereas in offset AGN only one nucleus is active. To study the origin of such AGN pairs, we employ a cosmological, hydrodynamic simulation with a large volume of (182Mpc) 3 from the set of Magneticum Pathfinder Simulations. The simulation selfconsistently produces 35 resolved black hole (BH) pairs at redshift z = 2, with a comoving distance smaller than 10kpc. 14 of them are offset AGN and nine are dual AGN, resulting in a fraction of (1.2 ± 0.3)% AGN pairs with respect to the total number of AGN. In this paper, we discuss fundamental differences between the BH and galaxy properties of dual AGN, offset AGN and inactive BH pairs and investigate their different triggering mechanisms. We find that in dual AGN the BHs have similar masses and the corresponding BH from the less massive progenitor galaxy always accretes with a higher Eddington ratio. In contrast, in offset AGN the active BH is typically more massive than its non-active counterpart. Furthermore, dual AGN in general accrete more gas from the intergalactic medium than offset AGN and nonactive BH pairs. This highlights that merger events, particularly minor mergers, do not necessarily lead to strong gas inflows and thus, do not always drive strong nuclear activity.

We present evidence from cosmological hydrodynamical simulations for a co-evolution of the slope of the total (dark and stellar) mass density profile, γ tot , and the dark matter fraction within the half-mass radius, f DM , in early-type galaxies. The relation can be described as γ tot = A f DM + B for all systems at all redshifts. The trend is set by the decreasing importance of gas dissipation towards lower redshifts and for more massive systems. Early-type galaxies are smaller, more concentrated, have lower f DM and steeper γ tot at high redshifts and at lower masses for a given redshift; f DM and γ tot are good indicators for growth by "dry" merging. The values for A and B change distinctively for different feedback models, and this relation can be used as a test for such models. A similar correlation exists between γ tot and the stellar mass surface density Σ * . A model with weak stellar feedback and feedback from black holes is in best agreement with observations. All simulations, independent of the assumed feedback model, predict steeper γ tot and lower f DM at higher redshifts. While the latter is in agreement with the observed trends, the former is in conflict with lensing observations, which indicate constant or decreasing γ tot . This discrepancy is shown to be artificial: the observed trends can be reproduced from the simulations using observational methodology to calculate the total density slopes.