Rosa Dominguez-Tenreiro

Chemical evolution of galaxies – I. A composition-dependent SPH model for chemical evolution and cooling

We describe an smooth particle hydrodynamics (SPH) model for chemical enrichment and radiative cooling in cosmological simulations of structure formation. This model includes: (i) the delayed gas restitution from stars by means of a probabilistic approach designed to reduce the statistical noise and, hence, to allow for the study of the inner chemical structure of objects with moderately high numbers of particles; (ii) the full dependence of metal production on the detailed chemical composition of stellar particles by using, for the first time in SPH codes, the Q i j matrix formalism that relates each nucleosynthetic product to its sources and (iii) the full dependence of radiative cooling on the detailed chemical composition of gas particles, achieved through a fast algorithm using a new metallicity parameter ζ(T) that gives the weight of each element on the total cooling function. The resolution effects and the results obtained from this SPH chemical model have been tested by comparing its predictions in different problems with known theoretical solutions. We also present some preliminary results on the chemical properties of elliptical galaxies found in self-consistent cosmological simulations. Such simulations show that the above ζ-cooling method is important to prevent an overestimation of the metallicity-dependent cooling rate, whereas the Q ij formalism is important to prevent a significant underestimation of the [a/Fe] ratio in simulated galaxy-like objects.

Dark Matter Halo Structure in CDM Hydrodynamical Simulations

We have carried out a comparative analysis of the properties of dark matter haloes in N-body and hydrodynamical simulations. We analyse their density profiles, shapes and kinematical properties with the aim of assessing the effects that hydrodynamical processes might produce on the evolution of the dark matter component. The simulations performed allow us to reproduce dark matter haloes with high resolution, although the range of circular velocities is limited. We find that for haloes with circular velocities of [150–200] km s−1 at the virial radius, the presence of baryons affects the evolution of the dark matter component in the central region, modifying the density profiles, shapes and velocity dispersions. We also analyse the rotation velocity curves of disc-like structures and compare them with observational results.

A two-phase scenario for bulge assembly in ΛcDM cosmologies

We analyze and compare the bulges of a sample of L * spiral galaxies in hydrodynamical simulations in a cosmological context, using two different codes, P-DEVA and GASOLINE. The codes regulate star formation in very different ways, with P-DEVA simulations inputting low star formation efficiency under the assumption that feedback occurs on subgrid scales, while the GASOLINE simulations have feedback that drives large-scale outflows. In all cases, the marked knee shape in mass aggregation tracks, corresponding to the transition from an early phase of rapid mass assembly to a later slower one, separates the properties of two populations within the simulated bulges. The bulges analyzed show an important early starburst resulting from the collapse-like fast phase of mass assembly, followed by a second phase with lower star formation, driven by a variety of processes such as disk instabilities and/or mergers. Classifying bulge stellar particles identified at z = 0 into old and young according to these two phases, we found bulge stellar sub-populations with distinct kinematics, shapes, stellar ages, and metal contents. The young components are more oblate, generally smaller, more rotationally supported, with higher metallicity and less alpha-element enhanced than the old ones. These results are consistent with the current observational status of bulges, and provide an explanation for some apparently paradoxical observations, such as bulge rejuvenation and metal-content gradients observed. Our results suggest that bulges of L * galaxies will generically have two bulge populations that can be likened to classical and pseudo-bulges, with differences being in the relative proportions of the two, which may vary due to galaxy mass and specific mass accretion and merger histories.

DISK GALAXIES WITH BROKEN LUMINOSITY PROFILES FROM COSMOLOGICAL SIMULATIONS

We present SPH cosmological simulations of the formation of three disk galaxies with a detailed treatment of chemical evolution and cooling. The resulting galaxies have properties compatible with observations: relatively high disk-to-total ratios, thin stellar disks and good agreement with the Tully-Fisher and the luminosity-size relations. They present a break in the luminosity profile at 3:0 0:5 disk scale lengths, while showing an exponential mass profile without any apparent breaks, in line with recent observational results. Since the stellar mass profile is exponential, only differences in the stellar populations can be the cause of the luminosity break. Although we find a cutoff for the star formation rate imposed by a density threshold in our star formation model, it does not coincide with the luminosity break and is located at 4:3 0:4 disk scale lengths, with star formation going on between both radii. The color profiles and the age profiles are “U-shaped”, with the minimum for both profiles located approximately at the break radius. The SFR to stellar mass ratio increases until the break, explaining the coincidence of the break with the minimum of the age profile. Beyond the break we find a steep decline in the gas density and, consequently, a decline in the SFR and redder colors. We show that most stars (64-78%) in the outer disk originate in the inner disk and afterwards migrate there. Such stellar migrations are likely the main origin of the U-shaped age profile and, therefore, of the luminosity break. Subject headings: galaxies: formation — galaxies: evolution — galaxies: spiral — galaxies: stellar content — methods: N-body simulations

Galaxies going MAD: the Galaxy-Finder Comparison Project

With the ever-increasing size and complexity of fully self-consistent simulations of galaxy formation within the framework of the cosmic web, the demands upon object finders for these simulations have simultaneously grown. To this extent we initiated the Halo-Finder Comparison Project that gathered together all the experts in the field and has so far led to two comparison papers, one for dark matter field haloes, and one for dark matter subhaloes. However, as state-of-the-art simulation codes are perfectly capable of not only following the formation and evolution of dark matter but also accounting for baryonic physics, i.e. gas hydrodynamics, star formation, stellar feedback, etc., object finders should also be capable of taking these additional physical processes into consideration. Here we report - for the first time - on a comparison of codes as applied to the Constrained Local UniversE Simulation (CLUES) of the formation of the Local Group which incorporates much of the physics relevant for galaxy formation. We compare both the properties of the three main galaxies in the simulation (representing the Milky Way, Andromeda and M33) and their satellite populations for a variety of halo finders ranging from phase space to velocity space to spherical overdensity based codes, including also a mere baryonic object finder. We obtain agreement amongst codes comparable to (if not better than) our previous comparisons - at least for the total, dark and stellar components of the objects. However, the diffuse gas content of the haloes shows great disparity, especially for low-mass satellite galaxies. This is primarily due to differences in the treatment of the thermal energy during the unbinding procedure. We acknowledge that the handling of gas in halo finders is something that needs to be dealt with carefully, and the precise treatment may depend sensitively upon the scientific problem being studied.

The early phases of galaxy clusters formation in IR: coupling hydrodynamical simulations with GRASIL-3D

We compute and study the infrared and sub-mm properties of high redshift (

grasil-3d: an implementation of dust effects in the SEDs of simulated galaxies

z \gtrsim 1

Origin of the metallicity distribution in the thick disc

) simulated clusters and proto-clusters. The results of a large set of hydro-dynamical zoom-in simulations including active galactic nuclei (AGN) feedback, have been treated with the recently developed radiative transfer code GRASIL-3D, which accounts for the effect of dust reprocessing in an arbitrary geometry. Here, we have slightly generalized the code to adapt it to the present purpose. Then we have post-processed boxes of physical size 2 Mpc encompassing each of the 24 most massive clusters identified at z=0, at several redshifts between 0.5 and 3, producing IR and sub-mm mock images of these regions and SEDs of the radiation coming out from them. While this field is in its infancy from the observational point of view, rapid development is expected in the near future thanks to observations performed in the far IR and sub-mm bands. Notably, we find that in this spectral regime our prediction are little affected by the assumption required by this post-processing, and the emission is mostly powered by star formation rather than accretion onto super massive black hole (SMBH). The comparison with the little observational information currently available, highlights that the simulated cluster regions never attain the impressive star formation rates suggested by these observations. This problem becomes more intriguing taking into account that the brightest cluster galaxies (BCGs) in the same simulations turn out to be too massive. It seems that the interplay between the feedback schemes and the star formation model should be revised, possibly incorporating a positive feedback mode.

Primordial nucleosynthesis bounds on the Brans-Dicke theory

We introduce a new model for the spectral energy distribution of galaxies, GRASIL-3D, which includes a careful modelling of the dust component of the interstellar medium. GRASIL-3D is an entirely new model based on the formalism of an existing and widely applied spectrophotometric model, GRASIL, but specifically designed to be interfaced with galaxies with any arbitrarily given geometry, such as galaxies calculated by theoretical hydrodynamical galaxy formation codes. GRASIL-3D is designed to separately treat radiative transfer in molecular clouds and in the diffuse cirrus component. The code has a general applicability to the outputs of simulated galaxies, either from Lagrangian or Eulerian hydrodynamic codes. As an application, the new model has been interfaced to the P-DEVA and GASOLINE smoothed-particle hydrodynamic codes, and has been used to calculate the spectral energy distribution for a variety of simulated galaxies from UV to sub-millimeter wavelengths, whose comparison with observational data gives encouraging results. In addition, GRASIL-3D allows 2D images of such galaxies to be obtained, at several angles and in different bands.

NIHAO XIII: Clumpy discs or clumpy light in high-redshift galaxies?

Aims. Using a suite of cosmological chemodynamical disc galaxy simulations, we assess how (a) radial metallicity gradients evolve with scaleheight; (b) the vertical metallicity gradients change through the thick disc; and (c) the vertical gradient of the stellar rotation velocity varies through the disc. We compare with the Milky Way to search for analogous trends.Methods. We analyse five simulated spiral galaxies with masses comparable to the Milky Way. The simulations span a range of star formation and energy feedback strengths and prescriptions, particle- and grid-based hydrodynamical implementations, as well as initial conditions/assembly history. Disc stars are identified initially via kinematic decomposition, with a posteriori spatial cuts providing the final sample from which radial and vertical gradients are inferred.Results. Consistently, we find that the steeper, negative, radial metallicity gradients seen in the mid-plane flatten with increasing height away from the plane. In simulations with stronger (and/or more spatially-extended) feedback, the negative radial gradients invert, becoming positive for heights in excess of ~1 kpc. Such behaviour is consistent with that inferred from recent observations. Our measurements of the vertical metallicity gradients show no clear correlation with galactocentric radius, and are in good agreement with those observed in the Milky Way’s thick disc (locally). Each of the simulations presents a decline in rotational velocity with increasing height from the mid-plane, although the majority have shallower kinematic gradients than that of the Milky Way.Conclusions. Simulations employing stronger/more extended feedback prescriptions possess radial and vertical metallicity and kinematic gradients more in line with recent observations. The inverted, positive, radial metallicity gradients seen in the simulated thick stellar discs originate in a population of younger, more metal-rich, stars formed in situ, superimposed upon a background population of older migrators from the inner disc; the contrast provided by the former increases radially, due to the inside-out growth of the disc. A similar behaviour may be responsible for the same flattening as seen in the radial gradients with scaleheight in the Milky Way.