We investigate the stellar kinematics of a sample of galaxies extracted from the hydrodynamic cosmological Magneticum Pathfinder simulations out to $5$ half-mass radii. We construct differential radial stellar spin profiles quantified by the observationally widely used $\lambda_\mathrm{R}$ and the closely related $(V/\sigma)$ parameters. We find three characteristic profile shapes: profiles exhibiting a (i) peak within $2.5$ half-mass radii and a subsequent decrease (ii) continuous increase that plateaus at larger radii typically with a high amplitude (iii) completely flat behaviour typically with low amplitude, in agreement with observations. This shows that the kinematic state of the stellar component can vary significantly with radius, suggesting a distinct interplay between in-situ star formation and ex-situ accretion of stars. Following the evolution of our sample through time, we provide evidence that the accretion history of galaxies with decreasing profiles is dominated by the anisotropic accretion of low mass satellites that get disrupted beyond $ \sim 2.0$ half-mass radii, building up a stellar halo with non-ordered motion while maintaining the central rotation already present at $z=2$. In fact, at $z=2$ decreasing profiles are the predominant profile class. Hence, we can predict a distinct formation pathway for galaxies with a decreasing profile and show that the centre resembles an old embedded disk. Furthermore, we show that the radius of the kinematic transition provides a good estimation for the transition radius from in-situ stars in the centre to accreted stars in the halo.

The effect of galactic orbits on a galaxy's internal evolution within a galaxy cluster environment has been the focus of heated debate in recent years. To understand this connection, we use both the $(0.5 \,$Gpc)$^3$ and the Gpc$^3$ boxes from the cosmological hydrodynamical simulation set Magneticum Pathfinder. We investigate the velocity-anisotropy, phase space, and the orbital evolution of up to $\sim 5 \cdot 10^{5}$ resolved satellite galaxies within our sample of 6776 clusters with $M_{\mathrm{vir}} > 10^{14} \, \mathrm{M_{\odot}}$ at low redshift, which we also trace back in time. In agreement with observations, we find that star-forming satellite galaxies inside galaxy clusters are characterised by more radially dominated orbits, independent of cluster mass. Furthermore, the vast majority of star-forming satellite galaxies stop forming stars during their first passage. We find a strong dichotomy both in line-of-sight and radial phase space between star-forming and quiescent galaxies, in line with observations. The tracking of individual orbits shows that the star-formation of almost all satellite galaxies drops to zero within $1 \, \mathrm{Gyr}$ after in-fall. Satellite galaxies that are able to remain star-forming longer are characterised by tangential orbits and high stellar mass. All this indicates that in galaxy clusters the dominant quenching mechanism is ram-pressure stripping.

Cosmological hydrodynamical simulations are rich tools to understand the build-up of stellar mass and angular momentum in galaxies, but require some level of calibration to observations. We compare predictions at $z\sim0$ from the Eagle, Hydrangea, Horizon-AGN, and Magneticum simulations with integral field spectroscopic (IFS) data from the SAMI Galaxy Survey, ATLAS3D, CALIFA and MASSIVE surveys. The main goal of this work is to simultaneously compare structural, dynamical, and stellar population measurements in order to identify key areas of success and tension. We have taken great care to ensure that our simulated measurement methods match the observational methods as closely as possible. We find that the Eagle and Hydrangea simulations reproduce many galaxy relations but with some offsets at high stellar masses. There are moderate mismatches in $R_e$ (+), $\epsilon$ (-), $\sigma_e$ (-), and mean stellar age (+), where a plus sign indicates that quantities are too high on average, and minus sign too low. The Horizon-AGN simulations qualitatively reproduce several galaxy relations, but there are a number of properties where we find a quantitative offset to observations. Massive galaxies are better matched to observations than galaxies at low and intermediate masses. Overall, we find mismatches in $R_e$ (+), $\epsilon$ (-), $\sigma_e$ (-) and $(V/\sigma)_e$ (-). Magneticum matches observations well: this is the only simulation where we find ellipticities typical for disk galaxies, but there are moderate differences in $\sigma_e$ (-), $(V/\sigma)_e$ (-) and mean stellar age (+). Our comparison between simulations and observational data has highlighted several areas for improvement, such as the need for improved modelling resulting in a better vertical disk structure, yet our results demonstrate the vast improvement of cosmological simulations in recent years.

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.

State-of-the-art integral field surveys like $\mathrm{ATLAS^{3D}}$, SLUGGS, CALIFA, SAMI, and MaNGA provide large data sets of kinematical observations of early-type galaxies (ETGs), yielding constraints on the formation of ETGs. Using the cosmological hydrodynamical \textit{Magneticum Pathfinder} simulations, we investigate the paradigm of fast and slow rotating ETGs in a fully cosmological context. We show that the ETGs within the \textit{Magneticum} simulation are in remarkable agreement with the observations, revealing fast and slow rotators quantified by the angular momentum proxy $\lambda_{\mathrm{R}}$ and the flattening $\epsilon$ with the observed prevalence. Taking full advantage of the three-dimensional data, we demonstrate that the dichotomy between fast and slow rotating galaxies gets enhanced, showing an upper and lower population separated by an underpopulated region in the edge-on $\lambda_{\mathrm{R}}$-$\epsilon$ plane. Following the evolution of the $\lambda_{\mathrm{R}}$-$\epsilon$ plane through cosmic time, we find that, while the upper population is already in place at $z=2$, the lower population gets statistically significant below $z=1$ with a gradual increase. At least $50\%$ of the galaxies transition from fast to slow rotators on a short timescale, in most cases associated to a significant merger event. Furthermore, we connect the $M_{*}$-$j_{*}$ plane, quantified by the $b$-value, with the $\lambda_{\mathrm{R}}$-$\epsilon$ plane, revealing a strong correlation between the position of a galaxy in the $\lambda_{\mathrm{R}}$-$\epsilon$ plane and the $b$-value. Going one step further, we classify our sample based on features in their velocity map, finding all five common kinematic groups, also including the recently observed group of prolate rotators, populating distinct regions in the $\lambda_{\mathrm{R}}$-$b$ plane.

Which mechanism(s) are mainly driving nuclear activity in the centres of galaxies is a major unsettled question. In this study, we investigate the statistical relevance of galaxy mergers for fuelling gas onto the central few kpc of a galaxy, potentially resulting in an active galactic nucleus (AGN). To robustly address that, we employ large-scale cosmological hydrodynamic simulations from the Magneticum Pathfinder set, including models for BH accretion and AGN feedback. Our simulations predict that for luminous AGN ($L_{\rm AGN} > 10^{45} {\rm erg/s}$) at $z = 2$, more than 50 per cent of their host galaxies have experienced a merger in the last 0.5~Gyr. These high merger fractions, however, merely reflect the intrinsically high merger fractions of massive galaxies at $z=2$, in which luminous AGN preferentially occur. Apart from that, our simulations suggest that merger events are not the statistically dominant fuelling mechanism for nuclear activity over a redshift range $z=0-2$: irrespective of AGN luminosity, less than 20 per cent of AGN hosts have on average undergone a recent merger, in agreement with a number of observational studies. The central ISM conditions required for inducing AGN activity can be, but are not necessarily caused by a merger. Despite the statistically minor relevance of mergers, at a given AGN luminosity and stellar mass, the merger fractions of AGN hosts can be by up to three times higher than that of inactive galaxies. Such elevated merger fractions still point towards an intrinsic connection between AGN and mergers, consistent with our traditional expectation.

Massive early-type galaxies typically have two subpopulations of globular clusters (GCs) which often reveal radial colour (metallicity) gradients. Collating gradients from the literature, we show that the gradients in the metal-rich and metal-poor GC subpopulations are the same, within measurement uncertainties, in a given galaxy. Furthermore, these GC gradients are similar in strength to the {\it stellar} metallicity gradient of the host galaxy. At the very largest radii (e.g. greater than 8 galaxy effective radii) there is some evidence that the GC gradients become flat with near constant mean metallicity. Using stellar metallicity gradients as a proxy, we probe the assembly histories of massive early-type galaxies with hydrodynamical simulations from the Magneticum suite of models. In particular, we measure the stellar metallicity gradient for the in-situ and accreted components over a similar radial range as those observed for GC subpopulations. We find that the in-situ and accreted stellar metallicity gradients are similar but have a larger scatter than the metal-rich and metal-poor GC subpopulations gradients in a given galaxy. We conclude that although metal-rich GCs are predominately formed during the in-situ phase and metal-poor GCs during the accretion phase of massive galaxy formation, they do not have a strict one-to-one connection.

We apply the Jeans Anisotropic MGE (JAM) dynamical modelling method to SAGES Legacy Unifying Globulars and GalaxieS (SLUGGS) survey data of early-type galaxies in the stellar mass range $10^{10}<M_*/{\rm M}_{\odot}<10^{11.6}$ that cover a large radial range of $0.1-4.0$ effective radii. We combine SLUGGS and ATLAS$^{\rm 3D}$ datasets to model the total-mass profiles of a sample of 21 fast-rotator galaxies, utilising a hyperparameter method to combine the two independent datasets. The total-mass density profile slope values derived for these galaxies are consistent with those measured in the inner regions of galaxies by other studies. Furthermore, the total-mass density slopes ($\gamma_{\rm tot}$) appear to be universal over this broad stellar mass range, with an average value of $\gamma_{\rm tot}=-2.12\,\pm\,0.05$, i.e. slightly steeper than isothermal. We compare our results to model galaxies from the Magneticum and EAGLE cosmological hydrodynamic simulations, in order to probe the mechanisms that are responsible for varying total-mass density profile slopes. The simulated-galaxy slopes are shallower than the observed values by $\sim0.1-0.3$, indicating that the physical processes shaping the mass distributions of galaxies in cosmological simulations are still incomplete. For galaxies with $M_*>10^{10.7}{\rm M}_{\odot}$ in the Magneticum simulations, we identify a significant anticorrelation between total-mass density profile slopes and the fraction of stellar mass formed ex situ (i.e. accreted), whereas this anticorrelation is weaker for lower stellar masses, implying that the measured total mass density slopes for low-mass galaxies are less likely to be determined by merger activity.

Selecting disk galaxies from the cosmological, hydrodynamical simulation Magneticum Pathfinder we show that almost half of our poster child disk galaxies at $z=2$ show significantly declining rotation curves and low dark matter fractions, very similar to recently reported observations. These galaxies do not show any anomalous behavior, reside in standard dark matter halos and typically grow significantly in mass until $z = 0$, where they span all morphological classes, including disk galaxies matching present day rotation curves and observed dark matter fractions. Our findings demonstrate that declining rotation curves and low dark matter fractions in rotation dominated galaxies at $z=2$ appear naturally within the $\Lambda$CDM paradigm and reflect the complex baryonic physics, which plays a role at the peak epoch of star-formation. In addition, we find some dispersion dominated galaxies at $z=2$ which host a significant gas disk and exhibit similar shaped rotation curves as the disk galaxy population, rendering it difficult to differentiate between these two populations with currently available observation techniques.

Recent hydrodynamic cosmological simulations cover volumes up to Gpc^3 and resolve halos across a wide range of masses and environments, from massive galaxy clusters down to normal galaxies, while following a large variety of physical processes (star formation, chemical enrichment, AGN feedback) to allow a self-consistent comparison to observations at multiple wavelengths. Using the Magneticum simulations, we investigate the buildup of the diffuse stellar component (DSC) around massive galaxies within group and cluster environments. The DSC in our simulations reproduces the spatial distribution of the observed intracluster light (ICL) as well as its kinematic properties remarkably well. For galaxy clusters and groups we find that, although the DSC in almost all cases shows a clear separation from the brightest cluster galaxy (BCG) with regard to its dynamic state, the radial stellar density distribution in many halos is often characterized by a single Sersic profile, representing both the BCG component and the DSC, very much in agreement with current observational results. Interestingly, even in those halos that clearly show two components in both the dynamics and the spatial distribution of the stellar component, no correlation between them is evident.

We present evidence from cosmological hydrodynamical simulations for a co-evolution of the slope of the total (dark and stellar) mass density profile, gamma_tot, and the dark matter fraction within the half-mass radius, f_DM, in early-type galaxies. The relation can be described as gamma_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 gamma_tot at high redshifts and at lower masses for a given redshift; f_DM and gamma_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 gamma_tot and the stellar mass surface density Sigma_*. 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 gamma_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 gamma_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.

We present a case study of a early-type galaxy (ETG) hosting a kinematically distinct core (KDC) formed in a binary high resolution 1:1 spiral galaxy merger simulation. The runtime of the simulation is pushed up to 10Gyr to follow the complete evolution of various physical properties. To investigate the origin of the KDC, the stellar component residing within the KDC is dissected, revealing that the rotational signal is purely generated by stars that belong to the KDC for at least 0.5Gyr and are newly formed during the merging process. Following the orientation of the total stellar angular momentum of the KDC, we show that it performs a motion comparable to the precession of a gyroscope in a gravitational potential. We draw the conclusion that the motion of the KDC is a superposition of an intrinsic rotation and a global precession that gets gradually damped over cosmic time. Finally, the stability of the KDC over the complete runtime of the simulation is investigated by tracing the evolution of the widely used $\lambda_{R}$ parameter and the misalignment angle distribution. We find that the KDC is stable for about 3Gyr after the merger and subsequently disperses completely on a timescale of ~1.5Gyr.

Metals are ideal tracers of the baryonic cycle within halos. Their composition is a fossil record connecting the evolution of the various stellar components of galaxies to the interaction with the environment by in- and outflows. The Magneticum simulations allow to study halos across a large range of masses and environments, from massive galaxy clusters containing hundreds of galaxies down to isolated field galaxies. They include a detailed treatment of the chemo-energetic feedback from the stellar component and its evolution as well as feedback from the evolution of supermassive black holes. Following the detailed evolution of various metal species and their relative composition due to continuing enrichment of the IGM and ICM by SNIa, SNII and AGB winds of the evolving stellar population reveals the complex interplay of local star formation processes, mixing, global baryonic flows, secular galactic evolution and environmental processes. We present results from the Magneticum simulations on the chemical properties of simulated galaxies and galaxy clusters, carefully comparing them to observations. We show that the simulations already reach a very high level of realism within their complex descriptions of the chemo-energetic feedback, successfully reproducing a large number of observed properties and scaling relations. Our simulated galaxies clearly indicate that there are no strong secondary parameters (like star formation rates at fixed redshift) driving the scatter in these scaling relations. The remaining differences clearly point to detailed physical processes which have to be included into future simulations.

In the past years several authors studied the abundance of satellites around galaxies in order to better estimate the halo masses of host galaxies. To investigate this connection, we analyze galaxies with $M_\mathrm{star}\geq\,10^{10}\,M_{\odot}$ from the hydrodynamical cosmological simulation Magneticum. We find that the satellite fraction of centrals is independent of their morphology. With the exception of very massive galaxies at low redshift, our results do not support the assumption that the dark matter (DM) haloes of spheroidal galaxies are significantly more massive than those of disc galaxies at fixed $M_\mathrm{star}$. We show that the density-morphology-relation starts to build up at $z\sim2$ and is independent of the star-formation properties of central galaxies. We conclude that environmental quenching is more important for satellites than for centrals. Our simulations indicate that conformity is already in place at $z=2$, where formation redshift and current star-formation rate (SFR) of central and satellite galaxies correlate. Centrals with low SFRs have formed earlier (at fixed $M_\mathrm{star}$) while centrals with high SFR formed later, with typical formation redshifts well in agreement with observations. However, we confirm the recent observations that the apparent number of satellites of spheroidal galaxies is significantly larger than for disc galaxies. This difference completely originates from the inclusion of companion galaxies, i.e. galaxies that do not sit in the potential minimum of a DM halo. Thus, due to the density-morphological-relation the number of satellites is not a good tracer for the halo mass, unless samples are restricted to the central galaxies of DM haloes.

The outer stellar halos of galaxies contain vital information about the formation history of galaxies, since the relaxation timescales in the outskirts are long enough to keep the memory, while the information about individual formation events in the central parts has long been lost due to mixing, star formation and relaxation. To unveil some of the information encoded in these faint outer halo regions, we study the stellar outskirts of galaxies selected from a fully hydrodynamical high resolution cosmological simulation, called Magneticum. We find that the density profiles of the outer stellar halos of galaxies over a broad mass range can be well described by an Einasto profile. For a fixed total mass range, the free parameters of the Einasto fits are closely correlated. Galaxies which had more (dry) merger events tend to have lesser curved outer stellar halos, however, we find no indication that the amount of curvature is correlated with galaxy morphology. The Einasto-like shape of the outer stellar halo densities can also explain the observed differences between the Milky Way and Andromeda outer stellar halos.

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 self-consistently 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 \pm 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 non-active 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 an implementation of smoothed particle hydrodynamics (SPH) with improved accuracy for simulations of galaxies and the large-scale structure. In particular, we combine, implement, modify and test a vast majority of SPH improvement techniques in the latest instalment of the GADGET code. We use the Wendland kernel functions, a particle wake-up time-step limiting mechanism and a time-dependent scheme for artificial viscosity, which includes a high-order 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 largely 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 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, thinner and more extended and our results on galaxy clusters show 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 SPH, thus becoming the core of an efficient code for large cosmological simulations.

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 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 the Magneticum simulations. 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 percent smaller than the specific angular momentum of the total DM halo and shows a significant scatter.

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 now very 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 end - consistent 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.

The stellar velocity dispersion as a function of the galactocentric radius of an early-type galaxy can generally be well approximated by a power law $\sigma \propto r^{\beta}$. However, some observed dispersion profiles show a deviation from this fit at intermediate radii, usually between one and three $R_{\mathrm{eff}}$, where the velocity dispersion remains constant with radius, showing a bump-like behavior, which we term the "$\sigma$-"bump. To understand the origin of this $\sigma$-bump, we study a set of simulated early-type galaxies formed in major mergers. We find the $\sigma$-bump in all of our simulated early-type galaxies, with the size and position of the bump slightly varying from galaxy to galaxy, leading to the assumption that the bump is a characteristic of the major merger formation scenario. The feature can be seen both in the intrinsic and projected stellar velocity dispersions. In contrast to shells that form during the merger event but evolve with time and finally disappear, the $\sigma$-bump stays nearly constant with radius and is a permanent feature that is preserved until the end of the simulation. The $\sigma$-bump is not seen in the dark matter and gas components and we therefore conclude that it is a purely stellar feature of merger remnants.

Dynamical modeling and strong lensing data indicate that the total density profiles of early-type galaxies are close to isothermal, i.e., rho_tot ~ r^gamma with gamma approx -2. To understand the origin of this universal slope we study a set of simulated spheroids formed in isolated binary mergers as well as the formation within the cosmological framework. The total stellar plus dark matter density profiles can always be described by a power law with an index of gamma approx -2.1 with a tendency toward steeper slopes for more compact, lower-mass ellipticals. In the binary mergers the amount of gas involved in the merger determines the precise steepness of the slope. This agrees with results from the cosmological simulations where ellipticals with steeper slopes have a higher fraction of stars formed in situ. Each gas-poor merger event evolves the slope toward gamma ~ -2, once this slope is reached further merger events do not change it anymore. All our ellipticals have flat intrinsic combined stellar and dark matter velocity dispersion profiles. We conclude that flat velocity dispersion profiles and total density distributions with a slope of gamma ~ -2 for the combined system of stars and dark matter act as a natural attractor. The variety of complex formation histories as present in cosmological simulations, including major as well as minor merger events, is essential to generate the full range of observed density slopes seen for present-day elliptical galaxies.

We study the possible magnetization of cosmic voids by void galaxies. Recently, observations revealed isolated starforming galaxies within the voids. Furthermore, a major fraction of a voids volume is expected to be filled with magnetic fields of a minimum strength of about $10^{-15}$ G on Mpc scales. We estimate the transport of magnetic energy by cosmic rays (CR) from the void galaxies into the voids. We assume that CRs and winds are able to leave small isolated void galaxies shortly after they assembled, and then propagate within the voids. For a typical void, we estimate the magnetic field strength and volume filling factor depending on its void galaxy population and possible contributions of strong active galactic nuclei (AGN) which border the voids. We argue that the lower limit on the void magnetic field can be recovered, if a small fraction of the magnetic energy contained in the void galaxies or void bordering AGNs is distributed within the voids.