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.

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.

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.

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.

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.

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.