No Arabic abstract
We review recent developments in the field of chemodynamical simulations of elliptical galaxies, highlighting (in an admittedly biased fashion) the work conducted with our cosmological N-body/SPH code GCD+. We have demonstrated previously the recovery of several primary integrated early-type system scaling relations (e.g. colour-magnitude relation, L_X-T_X-[Fe/H]_X) when employing a phenomenological AGN heating scheme in conjunction with a self-consistent treatment of star formation, supernovae feedback, radiative cooling, chemical enrichment, and stellar/X-ray population synthesis. Here we emphasise characteristics derived from the full spatial information contained within the simulated dataset, including stellar and coronal morphologies, metallicity distribution functions, and abundance gradients.
We investigate how the stellar and gas-phase He abundances evolve as functions of time within simulated star-forming disc galaxies with different star formation histories. We make use of a cosmological chemodynamical simulation for galaxy formation and evolution, which includes star formation, as well as energy and chemical enrichment feedback from asymptotic giant branch stars, core-collapse supernovae, and Type Ia supernovae. The predicted relations between the He mass fraction, $Y$, and the metallicity, $Z$, in the interstellar medium of our simulated disc galaxies depend on the past galaxy star formation history. In particular, $dY/dZ$ is not constant and evolves as a function of time, depending on the specific chemical element that we choose to trace $Z$; in particular, $dY/dX_{text{O}}$ and $dY/dX_{text{C}}$ increase as functions of time, whereas $dY/dX_{text{N}}$ decreases. In the gas-phase, we find negative radial gradients of $Y$, due to the inside-out growth of our simulated galaxy discs as a function of time; this gives rise to longer chemical enrichment time scales in the outer galaxy regions, where we find lower average values for $Y$ and $Z$. Finally, by means of chemical evolution models, in the galactic bulge and inner disc, we predict steeper $Y$ versus age relations at high $Z$ than in the outer galaxy regions. We conclude that, for calibrating the assumed $Y$-$Z$ relation in stellar models, C, N, and C+N are better proxies for the metallicity than O, because they show steeper and less scattered relations.
We present a new chemodynamical code - Ramses-CH - for use in simulating the self-consistent evolution of chemical and hydrodynamical properties of galaxies within a fully cosmological framework. We build upon the adaptive mesh refinement code Ramses, which includes a treatment of self-gravity, hydrodynamics, star formation, radiative cooling, and supernovae feedback, to trace the dominant isotopes of C, N, O, Ne, Mg, Si, and Fe. We include the contribution of Type Ia and II supernovae, in addition to low- and intermediate-mass asymptotic giant branch stars, relaxing the instantaneous recycling approximation. The new chemical evolution modules are highly flexible and portable, lending themselves to ready exploration of variations in the underpining stellar and nuclear physics. We apply Ramses-CH to the cosmological simulation of a typical L* galaxy, demonstrating the successful recovery of the basic empirical constraints regarding, [{alpha}/Fe]-[Fe/H] and Type Ia/II supernovae rates.
We analyse the kinematics and chemistry of the bulge stars of two simulated disc galaxies using our chemodynamical galaxy evolution code GCD+. First we compare stars that are born inside the galaxy with those that are born outside the galaxy and are accreted into the centre of the galaxy. Stars that originate outside of the bulge are accreted into it early in its formation within 3 Gyrs so that these stars have high [alpha/Fe] as well as having a high total energy reflecting their accretion to the centre of the galaxy. Therefore, higher total energy is a good indicator for finding accreted stars. The bulges of the simulated galaxies formed through multiple mergers separated by about a Gyr. Since [alpha/Fe] is sensitive to the first few Gyrs of star formation history, stars that formed during mergers at different epochs show different [alpha/Fe]. We show that the [Mg/Fe] against star formation time relation can be very useful to identify a multiple merger bulge formation scenario, provided there is sufficiently good age information available. Our simulations also show that stars formed during one of the merger events retain a systematically prograde rotation at the final time. This demonstrates that the orbit of the ancient merger that helped to form the bulge could still remain in the kinematics of bulge stars.
(Abridged) We present an investigation of kinematical imprints of AGN feedback on the Warm Ionized gas Medium (WIM) of massive early-type galaxies (ETGs). To this end, we take a two-fold approach that involves a comparative analysis of Halpha velocity fields in 123 local ETGs from the CALIFA integral field spectroscopy survey with 20 simulated galaxies from high-resolution hydrodynamic cosmological SPHgal simulations. The latter were re-simulated for two modeling setups, one with and another without AGN feedback. In order to quantify the effects of AGN feedback on gas kinematics we measure three parameters that probe deviations from simple regular rotation using the kinemetry package. These indicators trace the possible presence of distinct kinematic components in Fourier space (k3,5/k1), variations in the radial profile of the kinematic major axis (sigma_PA), and offsets between the stellar and gas velocity fields (Delta Phi). These quantities are monitored in the simulations from a redshift 3 to 0.2 to assess the connection between black hole accretion history, stellar mass growth and kinematical perturbation of the WIM. Observed local massive galaxies show a broad range of irregularities, indicating disturbed warm gas motions, irrespective of being classified via diagnostic lines as AGN or not. Simulations of massive galaxies with AGN feedback generally exhibit higher irregularity parameters than without AGN feedback, more consistent with observations. Besides AGN feedback, other processes like major merger events or infalling gas clouds can lead to elevated irregularity parameters, but they are typically of shorter duration. More specifically, k3,5/k1 is most sensitive to AGN feedback, whereas Delta Phi is most strongly affected by gas infall.
The kinematics of stars and planetary nebulae in early type galaxies provide vital clues to the enigmatic physics of their dark matter halos. We fit published data for fourteen such galaxies using a spherical, self-gravitating model with two components: (1) a Sersic stellar profile fixed according to photometric parameters, and (2) a polytropic dark matter halo that conforms consistently to the shared gravitational potential. The polytropic equation of state can describe extended theories of dark matter involving self-interaction, non-extensive thermostatistics, or boson condensation (in a classical limit). In such models, the flat-cored mass profiles widely observed in disc galaxies are due to innate dark physics, regardless of any baryonic agitation. One of the natural parameters of this scenario is the number of effective thermal degrees of freedom of dark matter (F_d) which is proportional to the dark heat capacity. By default we assume a cosmic ratio of baryonic and dark mass. Non-Sersic kinematic ideosyncrasies and possible non-sphericity thwart fitting in some cases. In all fourteen galaxies the fit with a polytropic dark halo improves or at least gives similar fits to the velocity dispersion profile, compared to a stars-only model. The good halo fits usually prefer F_d values from six to eight. This range complements the recently inferred limit of 7<F_d<10 (Saxton & Wu), derived from constraints on galaxy cluster core radii and black hole masses. However a degeneracy remains: radial orbital anisotropy or a depleted dark mass fraction could shift our models preference towards lower F_d; whereas a loss of baryons would favour higher F_d.