Search

The distribution of chemical elements in the star-forming regions can store information on the chemical enrichment history of the galaxies. Negative metallicity gradients are expected in galaxies forming inside-out. However, observations show that the metallicity profiles can be broken. We aim to study the diversity of metallicity profiles that can arise in the current cosmological context and compare them with available observations. We also seek to identify the physical processes responsible for breaks in metallicity profiles by using two galaxies as case studies. We analyze central galaxies from the cosmological simulations of the CIELO project, within the stellar mass range [$10^{8.5}$, $10^{10.5}$] M$_\odot$ at $z=0$. A new algorithm, DB-A, was developed to fit multiple power laws to the metallicity profiles, enabling a flexible assessment of metallicity gradients in various galactic regions. The simulations include detailed modeling of gas, metal-dependent cooling, star formation, and supernova feedback. At $z=0$, we find diverse profile shapes, including inner and outer drops and rises, with some galaxies exhibiting double breaks. Gradient values align with observations. A temporal analysis of Local Group analogs shows inner and outer breaks occurring at all cosmic times, with outer breaks being more frequent. Metallicity gradients show high variability at high redshift, transitioning to mild evolution at lower redshift. Most inner breaks show central oxygen enhancement, linked to gas accretion and star formation. Inner drops result from disrupted gas due to feedback-driven outflows. Outer breaks with high metallicities arise from re-accreted material, extended star formation, and CGM-driven gas mixing. Outer drops are common at high redshift, linked to metal-poor gas accretion from cold flows. We highlight the complex interplay of these processes which often act together., Comment: 21 pages, 14 figures. Comments are welcome

The CIELO project introduces a novel set of chemo-dynamical zoom-in simulations designed to simultaneously resolve galaxies and their nearby environments. The initial conditions include a diverse range of cosmic structures, such as local groups, filaments, voids, and walls, allowing for a detailed exploration of galaxies within the broader context of the cosmic web. This study presents the initial conditions and characterizes the global properties of CIELO galaxies and their environments. It focuses on galaxies with stellar masses ranging from log [8,11] solar masses and examines key scaling relations, including the mass-size relation, the Tully-Fisher relation, and the mass-metallicity relation for both stars and star-forming gas. The DisPerSe algorithm was used to determine the positions of CIELO galaxies within the cosmic web, with a specific focus on the Pehuen haloes. The selection of local group volumes was guided by criteria based on the relative positions and velocities of the two primary galaxies. The Pehuen regions were chosen to map walls, filaments, and voids. Synthetic images in the SDSS i, r, and g bands were generated using the SKIRT radiative transfer code. Additionally, a dynamical decomposition was performed to classify galaxy morphologies into bulge, disc, and stellar halo components (abridged)., Comment: 13 pages, 12 figures, submitted. Comments are welcomed

Stellar halos around galaxies contain key information about their formation and assembly history. Using simulations, we can trace the origins of different stellar populations in these halos, contributing to our understanding of galaxy evolution. We aim to investigate the assembly of stellar halos and their chemical abundances in 28 galaxies from CIELO project with logMgal[9 and 11]Msun. Stellar halos were identified using the AM E method, focusing on the outer regions between the 1.5 optical radius and the virial radius. We divided the stellar populations based on their formation channel: exsitu, endodebris, and insitu, and analyzed their chemical abundances, ages, and spatial distributions. Additionally, we explored correlations between halo mass, metallicity, and alpha element enrichment. CIELO simulations reveal that stellar halos are predominantly composed of accreted material (exsitu and endodebris stars), in agreement with previous works. The mass fraction of these populations is independent of stellar halo mass, though their metallicities scale linearly with it. Exsitu stars tend to dominate the outskirts and be more alpha rich and older, while endodebris stars are more prevalent at lower radii and tend to be less alpha rich and slightly younger. Massive stellar halos require a median of five additional satellites to build 90 percent of their mass, compared to lower mass halos, which typically need fewer (median of 2.5) and lower-mass satellites and are assembled earlier. The diversity of accreted satellite histories results in well defined stellar halo mass metallicity and [alpha/Fe] [Fe/H] relations, offering a detailed view of the chemical evolution and assembly history of stellar halos. We find that the [alpha/Fe] [Fe/H] is more sensitive to the characteristics and star formation history of the contributing satellites than the stellar halo mass metallicity relationship, Comment: 16 pages, 14 figures, accepted for publication in A&A

Kinematically-persistent planes of satellites (KPPs) are fixed sets of satellites co-orbiting around their host galaxy, whose orbital poles are conserved and clustered across long cosmic time intervals. They play the role of 'skeletons', ensuring the long-term durability of positional planes. We explore the physical processes behind their formation in terms of the dynamics of the local Cosmic Web (CW), characterized via the so-called Lagrangian Volumes (LVs) built up around two zoom-in, cosmological hydro-simulations of MW-mass disk galaxy + satellites systems, where three KPPs have been identified. By analyzing the LVs deformations in terms of the reduced Tensor of Inertia (TOI), we find an outstanding alignment between the LV principal directions and KPP satellites' orbital poles. The most compressive local mass flows (along the $\hat{e}_3$ eigenvector) are strong at early times, feeding the so-called $\hat{e}_3$-structure, while the smallest TOI axis rapidly decreases. The $\hat{e}_3$-structure collapse marks the end of this regime and is the timescale for the establishment of satellite orbital pole clustering when the Universe is $\lesssim$ 4 Gyr old. KPP proto-satellites aligned with $\hat{e}_3$ are those whose orbital poles are either aligned from early times, or have been successfully bent at $\hat{e}_3$-structure collapse. KPP satellites associated to $\hat{e}_1$ tend to have early trajectories already parallel to $\hat{e}_3$. We show that KPPs can arise as a result of the $\Lambda$CDM-predicted large-scale dynamics acting on particular sets of proto-satellites, the same dynamics that shape the local CW environment., Comment: 27 pages, 12 figures, 2 tables. Accepted for publication in The Astrophysical Journal

In this work, we aim at investigating the morphology evolution of Milky Way mass-like dark matter haloes selected from the CIELO and IllustrisTNG Projects. The connection between halo shapes and their environment has been studied in previous works at z=0 but their connection remains yet to be fully understood. We focus on the evolution across cosmic time of the halo shapes and the relation with the infalling material, using hydrodynamical simulations. Our findings show that haloes tend to be more triaxial at earlier times as a consequence of stronger accretion in the direction of the filaments. As the haloes evolve towards a dominant isotropic accretion mode and relaxation, their shape at 20 percent of the virial mass becomes more spherical. In agreement with previous results, baryons have an important effect within the inner regions of the haloes, driving them from triaxial to rounder shapes. We also find a correlation between the strength of the quadrupole infalling mode and the degree of ellipticity of the haloes: as the filament strength decreases steadily with redshift, the haloes became more spherical and less elliptical.

We use two zoom-in $\Lambda$CDM hydrodynamical simulations of massive disk galaxies to study the possible existence of fixed satellite groups showing a kinematically-coherent behaviour across evolution (angular momentum conservation and clustering). We identify three such groups in the two simulations, defining kinematically-coherent, time-persistent planes (KPPs) that last at least from virialization to $z=0$ (more than 7 Gyrs). This proves that orbital pole clustering is not necessarily set in at low redshift, representing a long-lived property of galaxy systems. KPPs are thin and oblate, represent $\sim25-40\%$ of the total number of satellites in the system, and are roughly perpendicular to their corresponding central disk galaxies during certain periods, consistently with Milky Way $z=0$ data. KPP satellite members are statistically distinguishable from satellites outside KPPs: they show higher specific orbital angular momenta, orbit more perpendicularly to the central disk galaxy, and have larger pericentric distances, than the latter. We numerically prove, for the first time, that KPPs and the best-quality positional planes share the same space configuration across time, such that KPPs act as `skeletons' preventing the latter of being washed out in short timescales. In one of the satellite-host systems, we witness the late capture of a massive dwarf galaxy endowed with its own satellite system, also organized into a KPP configuration prior to its capture. We briefly explore the consequences this event has on the host's KPP, and on the possible enhancement of the asymmetry in the number of satellites rotating in one sense or the opposite within the KPP., Comment: 24 pages, 8 figures, 2 tables. Accepted for publication in The Astrophysical Journal

We search for detectable signatures of f(R) gravity and its chameleon screening mechanism in the baryonic and dark matter (DM) properties of simulated void galaxies. The enhancement of the gravitational acceleration can have a meaningful impact on the scaling relations as well as on the halo morphology. The galaxy rotational velocity field (calculated with the velocity of the gas disc and the acceleration fields) deviates from the typical values of the Tully-Fisher Relation (TFR) in GR. For a given stellar mass, f(R) gravity tends to produce greater maximum velocities. On the other hand, the mass in haloes in f(R) gravity is more concentrated than their counterparts in GR. This trend changes when the concentration is calculated with the dynamical density profile, which takes into account the unscreened outer regions of the halo. Stellar discs interact with the overall potential well in the central regions, modifying the morphology of the screening regions and reshaping them. We find a trend for galaxies with a more dominant stellar disc to deviate further from round screening regions. We find that small haloes are less triaxial and more round in f(R) than their GR counterparts. The difference between halo morphology becomes smaller in f(R) haloes whose inner regions are screened. These results suggest possible observables that could unveil modified gravity effects on galaxies in voids in future cosmological tests of gravity., Comment: 19 pages, 10 figures

The formation of galaxies can be understood in terms of the assembly patterns of each type of galactic component. To perform this kind of analysis, is necessary to define some criteria to separate those components. Decomposition methods based on dynamical properties are more physically motivated than photometry-based ones. We use the unsupervised Gaussian Mixture model of \texttt{galactic structure finder} to extract the components of a sub-sample of galaxies with Milky Way-like masses from the EAGLE simulations. A clustering in the space of first and second order dynamical moments of all identified substructures reveals five types of galaxy components: thin and thick disks, stellar halos, bulges and spheroids. We analyse the dynamical, morphological and stellar population properties of these five component types, exploring to what extent these properties correlate with each other, and how much they depend on the total galaxy stellar and dark matter halo masses. All galaxies contain a bulge, a stellar halo and a disk. 60% of objects host two disks (thin and thick), and 68% host also a spheroid. The dynamical disk-to-total ratio does not depend on stellar mass, but the median rotational velocities of the two disks do. Thin disks are well separated in stellar ages, [Fe/H] and $\alpha$-enhancement from the three dispersion-dominated components, while thick disks are in between. Except for thin disks, all components show correlations among their stellar population properties: older ages mean lower metallicities and larger $\alpha$-enhancement. Finally, we quantify the weak dependence of stellar population properties on each component's dynamics., Comment: 20 pages, 18 figures. Accepted for publication in MNRAS

We analyse the evolution of the oxygen abundance gradient of star-forming galaxies with stellar mass Mstar > 10^9 Mo in the EAGK simulation over the redshift range z=[0, 2.5]. We find that the median metallicity gradient of the simulated galaxies is close to zero at all z, whereas the scatter around the median increases with z. The metallicity gradients of individual galaxies can evolve from strong to weak and vice-versa, since mostly low-metallicity gas accretes onto the galaxy, resulting in enhanced star formation and ejection of metal enriched gas by energy feedback. Such episodes of enhanced accretion, mainly dominated by major mergers, are more common at higher z, and hence contribute to increasing the diversity of gradients. For galaxies with negative metallicity gradients, we find a redshift evolution of ~ -0.03 dex/kpc/\delta z$. A positive mass dependence is found at z< 0.5, which becomes slightly stronger for higher redshifts and, mainly, for Mstar < 10^9.5 Mo. Only galaxies with negative metallicity gradients define a correlation with galaxy size, consistent with an inside-out formation scenario. Our findings suggest that major mergers and/or significant gas accretion can drive strong negative or positive metallicity gradients. The first ones are preferentially associated with disc-dominated galaxies, and the second ones with dispersion-dominated systems. The comparison with forthcoming observations at high redshift will allow a better understanding of the potential role of metallicity gradients as a chemical probe of galaxy formation., Comment: 19 pages, 12 figures. Accepted MNRAS

We address the 'plane of satellites problem' by studying planar configurations around two disc galaxies with no late major mergers, formed in zoom-in hydro-simulations. Due to the current lack of good quality kinematic data for M31 satellites, we use only positional information. So far, positional analyses of simulations are unable to find planes as thin and populated as the observed ones. Moreover, they miss systematicity and detail in the plane-searching techniques, as well as in the study of the properties and quality of planes, both in simulations or real data. To fill this gap, i) we extend the 4-galaxy-normal density plot method (Pawlowski et al. 2013) in a way designed to efficiently identify the best quality planes (i.e., thin and populated) without imposing extra constraints on their properties, and ii), we apply it for the first time to simulations. Using zoom-in simulations allows us to mimic MW/M31-like systems regarding the number of satellites involved as well as the galactic disc mass and morphology, in view of possible disc effects. At all timesteps analyzed in both simulations we find satellite planar configurations that are compatible, along given time intervals, with all the spatial characteristics of observed planes identified using the same methodology. However, the fraction of co-orbiting satellites within them is in general low, suggesting time-varying satellite membership. We conclude that high-quality positional planes of satellites are not infrequent in LCDM-formed disc galaxies with a quiet assembly history. Detecting kinematically-coherent, time-persistent planes demands considering the full six-dimensional phase-space information of satellites., Comment: 26 pages, 13 figures, 2 tables. Accepted for publication in ApJ

We study the optical morphology of galaxies in a large-scale hydrodynamic cosmological simulation, the EAGLE simulation. Galaxy morphologies were characterized using non-parametric statistics (Gini, $M_{20}$, Concentration and Asymmetry) derived from mock images computed using a 3D radiative transfer technique and post-processed to approximate observational surveys. The resulting morphologies were contrasted to observational results from a sample of $\log_{10}(M_{*}/M_\odot) > 10$ galaxies at $z \sim 0.05$ in the GAMA survey. We find that the morphologies of EAGLE galaxies reproduce observations, except for asymmetry values which are larger in the simulated galaxies. Additionally, we study the effect of spatial resolution in the computation of non-parametric morphologies, finding that Gini and Asymmetry values are systematically reduced with decreasing spatial resolution. Gini values for lower mass galaxies are especially affected. Comparing against other large scale simulations, the non-parametric statistics of EAGLE galaxies largely agree with those found in IllustrisTNG. Additionally, EAGLE galaxies mostly reproduce observed trends between morphology and star formation rate and galaxy size. Finally, We also find a significant correlation between optical and kinematic estimators of morphologies, although galaxy classification based on an optical or a kinematic criteria results in different galaxy subsets. The correlation between optical and kinematic morphologies is stronger in central galaxies than in satellites, indicating differences in morphological evolution., Comment: 21 pages, 12 Figures. Accepted for publication in MNRAS

Current observational results show that both late-and-early-type galaxies follow tight mass-size planes, on which physical properties such as age, velocity dispersion and metallicities correlate with the scatter on the plane. We study the mass-size plane of galaxies in cosmological hydrodynamical simulations, as a function of velocity dispersion, age, chemical abundances, ellipticity and spin parameters with the aim at assessing to what extent the current cosmological paradigm can reproduce these observations and provide a physical interpretation of them. We select a sample of well-resolved galaxies from the (100 Mpc)^3 simulation of the EAGLE Project. This sample is composed by 508 spheroid-dominated galaxies and 1213 disc-dominated galaxies. The distributions of velocity dispersion, age, metallicity indicators and gradients and spin parameters across the mass-size plane are analysed. Furthermore, we study the relation between shape and kinematic parameters. The results are compared with observations. The mass-weighted ages of the EAGLE galaxies are found to vary along lines of constant velocity dispersion on the mass-size plane, except for galaxies with velocity dispersion larger than aprox 150 km s^(-1) . Negative age gradients tend to be found in extended disc galaxies in agreement with observations. However, the age distributions of early-type galaxies show a larger fraction with inverted radial profiles. The distribution of metallicity gradients does not show any clear dependence on this plane. Galaxies with similar spin parameters ({\lambda}) display larger sizes as their dynamical masses increase. Stellar-weighted ages are found to be good proxies for {\lambda} in galaxies with low ellipticity ({\epsilon}). Abridged, Comment: 10 pages, 10 figures, accepted for publication in A&A

Aims: It is well known that the presence of baryons affects the dark matter host haloes. Exploring the galaxy assembly history together with the dark matter haloes properties through time can provide a way to measure these effects. Methods: We study the properties of four Milky Way mass dark matter haloes from the Aquarius project during their assembly history, between $z = 0 - 4$. In this work, we use the SPH run from Scannapieco et al. (2009) and the dark matter only counterpart as case studies. To asses the robustness of our findings, we compare them with one of the haloes run using a moving-mesh technique and different sub-grid scheme. Results: Our results show that the cosmic evolution of the dark matter halo profiles depends on the assembly history of the baryons. We find that the dark matter profiles do not significantly change with time, hence they become stable, when the fraction of baryons accumulated in the central regions reaches 80 percent of its present mass within the virial radius. Furthermore, the mass accretion history shows that the haloes that assembled earlier are those that contain a larger amount of baryonic mass aforetime, which in turn allows the dark matter halo profiles to reach a stable configuration earlier. For the SPH haloes, we find that the specific angular momentum of the dark matter particles within the five percent of the virial radius at z = 0, remains approximately constant from the time at which 60 percent of the stellar mass is gathered. We explore different theoretical and empirical models for the contraction of the haloes through redshift. A model to better describe the contraction of the haloes through redshift evolution must depend on the stellar mass content in the inner regions., Comment: 17 pages, 14 figures. Accepted for publication in A&A. Comments are welcome

Galaxy formation is a very complex process in which many different physical mechanisms intervene. Within the LCDM paradigm processes such as gas inflows and outflows, mergers and interactions contribute to the redistribution of the angular momentum content of the structures. Recent observational results have brought new insights and also triggered several theoretical studies. Some of these new contributions will be analysed here., Comment: 6 pages, 4 figures, to be published in "Galactic Angular Momentum", IAU General Assembly, Focus Meeting 6, 2018

The metallicity gradients of the stellar populations in disc galaxies and their evolution store relevant information on the disc formation history and on those processes which could mix stars a posteriori, such as migration, bars and/or galaxy-galaxy interactions. We aim to investigate the evolution of the metallicity gradients of the whole stellar populations in disc components of simulated galaxies in a cosmological context. We analyse simulated disc galaxies selected from a cosmological hydrodynamical simulation that includes chemical evolution and a physically motivated Supernova feedback capable of driving mass-loaded galactic winds. We detect a mild evolution with redshift in the metallicity slopes of $-0.02 \pm 0.01$ dex~kpc$^{-1}$ from $z\sim 1$. If the metallicity profiles are normalised by the effective radius of the stellar disc, the slopes show no clear evolution for $z < 1$, with a median value of approximately $-0.23$ dex ~$r_{\rm reff}^{-1}$. As a function of stellar mass, we find that metallicity gradients steepen for stellar masses smaller than $\sim 10^{10.3} {\rm M_{\odot}}$ while the trend reverses for higher stellar masses, in the redshift range $z=[0,1]$. Galaxies with small stellar masses have discs with larger $r_{\rm reff}$ and flatter metallicity gradients than expected. We detect migration albeit weaker than in previous works. Our stellar discs show a mild evolution of the stellar metallicity slopes up to $z\sim 1,$ which is well-matched by the evolution calculated archeologically from the abundance distributions of mono-age stellar populations at $z\sim 0$. Overall, Supernova feedback could explain the trends but an impact of migration can not be totally discarded. Galaxy-galaxy interactions or small satellite accretions can also contribute to modify the metallicity profiles in the outer parts. [abridged], Comment: 8 pages, 6 figures, accepted for publication in A&A

We study present-day galaxy clustering in the EAGLE cosmological hydrodynamical simulation. EAGLE's galaxy formation parameters were calibrated to reproduce the redshift $z=0.1$ galaxy stellar mass function, and the simulation also reproduces galaxy colours well. The simulation volume is too small to correctly sample large-scale fluctuations and we therefore concentrate on scales smaller than a few mega parsecs. We find very good agreement with observed clustering measurements from the Galaxy And Mass Assembly (GAMA) survey, when galaxies are binned by stellar mass, color, or luminosity. However, low-mass red-galaxies are clustered too strongly, which is at least partly due to limited numerical resolution. Apart from this limitation, we conclude that EAGLE galaxies inhabit similar dark matter haloes as observed GAMA galaxies, and that the radial distribution of satellite galaxies as function of stellar mass and colour is similar to that observed as well., Comment: 18 pages, 15 figures. Accepted for publication in MNRAS

We study non-parametric morphologies of mergers events in a cosmological context, using the Illustris project. We produce mock g-band images comparable to observational surveys from the publicly available Illustris simulation idealized mock images at $z=0$. We then measure non parametric indicators: asymmetry, Gini, $M_{20}$, clumpiness and concentration for a set of galaxies with $M_* >10^{10}$ M$_\odot$. We correlate these automatic statistics with the recent merger history of galaxies and with the presence of close companions. Our main contribution is to assess in a cosmological framework, the empirically derived non-parametric demarcation line and average time-scales used to determine the merger rate observationally. We found that 98 per cent of galaxies above the demarcation line have a close companion or have experienced a recent merger event. On average, merger signatures obtained from the $G-M_{20}$ criteria anticorrelate clearly with the elapsing time to the last merger event. We also find that the asymmetry correlates with galaxy pair separation and relative velocity, exhibiting the larger enhancements for those systems with pair separations $d < 50$ h$^{-1}$ kpc and relative velocities $V < 350$ km s$^{-1}$. We find that the $G-M_{20}$ is most sensitive to recent mergers ($\sim0.14$ Gyr) and to ongoing mergers with stellar mass ratios greater than 0.1. For this indicator, we compute a merger average observability time-scale of $\sim0.2$ Gyr, in agreement with previous results and demonstrate that the morphologically derived merger rate recovers the intrinsic total merger rate of the simulation and the merger rate as a function of stellar mass., Comment: 19 pages, 20 figures. Published in MNRAS

The stellar metallicity gradients of disc galaxies provide information on the disc assembly, star formation processes and chemical evolution. They also might store information on dynamical processes which could affect the distribution of chemical elements in the gas-phase and the stellar components. We studied the stellar metallicity gradients of stellar discs in a cosmological simulation. We explored the dependence of the stellar metallicity gradients on stellar age and the size and mass of the stellar discs. We used galaxies selected from a cosmological hydrodynamical simulation performed including a physically-motivated Supernova feedback and chemical evolution. The metallicity profiles were estimated for stars with different ages. We confront our numerical findings with results from the CALIFA Survey. The simulated stellar discs are found to have metallicity profiles with slopes in global agreement with observations. Low stellar-mass galaxies tend to have a larger variety of metallicity slopes. When normalized by the half-mass radius, the stellar metallicity gradients do not show any dependence and the dispersion increases significantly, regardless of galaxy mass. Galaxies with stellar masses around $10^{10}$M$_{\odot}$ show steeper negative metallicity gradients. The stellar metallicity gradients correlate with the half-mass radius. However, the correlation signal is not present when they are normalized by the half-mass radius. Stellar discs with positive age gradients are detected to have negative and positive metallicity gradients, depending on the relative importance of the recent star formation activity in the central regions. The large dispersions in the metallicity gradients as a function of stellar mass could be ascribed to the effects of dynamical processes such as mergers/interactions and/or migration as well as those regulating the conversion of gas into stars. [abridged], Comment: 9 pages, 9 figures, accepted for publication in A&A

We study oxygen abundance profiles of the gaseous disc components in simulated galaxies in a hierarchical universe. We analyse the disc metallicity gradients in relation to the stellar masses and star formation rates of the simulated galaxies. We find a trend for galaxies with low stellar masses to have steeper metallicity gradients than galaxies with high stellar masses at z ~0. We also detect that the gas-phase metallicity slopes and the specific star formation rate (sSFR) of our simulated disc galaxies are consistent with recently reported observations at z ~0. Simulated galaxies with high stellar masses reproduce the observed relationship at all analysed redshifts and have an increasing contribution of discs with positive metallicity slopes with increasing redshift. Simulated galaxies with low stellar masses a have larger fraction of negative metallicity gradients with increasing redshift. Simulated galaxies with positive or very negative metallicity slopes exhibit disturbed morphologies and/or have a close neighbour. We analyse the evolution of the slope of the oxygen profile and sSFR for a gas-rich galaxy-galaxy encounter, finding that this kind of events could generate either positive and negative gas-phase oxygen profiles depending on their state of evolution. Our results support claims that the determination of reliable metallicity gradients as a function of redshift is a key piece of information to understand galaxy formation and set constrains on the subgrid physics., Comment: 12 pages, 8 figures, accepted MNRAS

Galaxy formation in the current cosmological paradigm is a very complex process in which inflows, outflows, interactions and mergers are common events. These processes can redistribute the angular momentum content of baryons. Recent observational results suggest that disc formed conserving angular momentum while elliptical galaxies, albeit losing angular momentum, determine a correlation between the specific angular momentum of the galaxy and the stellar mass. These observations provide stringent constraints for galaxy formation models in a hierarchical clustering scenario. We aim to analyse the specific angular momentum content of the disc and bulge components as a function of virial mass, stellar mass and redshift. We also estimate the size of the simulated galaxies and confront them with observations. We use cosmological hydrodynamical simulations that include an effective, physically-motivated Supernova feedback which is able to regulate the star formation in haloes of different masses. We analyse the morphology and formation history of a sample of galaxies in a cosmological simulation by performing a bulge-disc decomposition of the analysed systems and their progenitors. We estimate the angular momentum content of the stellar and gaseous discs, stellar bulges and total baryons. In agreement with recent observational findings, our simulated galaxies have disc and spheroid components whose specific angular momentum contents determine correlations with the stellar and dark matter masses with the same slope, although the spheroidal components are off-set by a fixed fraction. Abridged., Comment: 9 pages, 8 figures, A&A accepted

Using cosmological hydrodynamics simulations we study the angular momentum content of the simulated galaxies in relation with their morphological type. We found that not only the angular momentum of the disk component follow the expected theoretical relation, Mo, Mao & Whiye (1998), but also the spheroidal one, with a gap due to its lost of angular momentum, in agreement with Fall & Romanowsky (2013),. We also found that the galaxy size can plot in one general relation, despite the morphological type, as found by Kravtsov (2013)., Comment: 1 page, BAAA proceedings (accepted with minor changes), AAA Meeeting, september 2014

We analyse the dark, gas, and stellar mass assembly histories of low-mass halos (Mvir ~ 10^10.3 - 10^12.3 M_sun) identified at redshift z = 0 in cosmological numerical simulations. Our results indicate that for halos in a given present-day mass bin, the gas-to-baryon fraction inside the virial radius does not evolve significantly with time, ranging from ~0.8 for smaller halos to ~0.5 for the largest ones. Most of the baryons are located actually not in the galaxies but in the intrahalo gas; for the more massive halos, the intrahalo gas-to-galaxy mass ratio is approximately the same at all redshifts, z, but for the least massive halos, it strongly increases with z. The intrahalo gas in the former halos gets hotter with time, being dominant at z = 0, while in the latter halos, it is mostly cold at all epochs. The multiphase ISM and thermal feedback models in our simulations work in the direction of delaying the stellar mass growth of low-mass galaxies., Comment: Accepted for publication in BAAA, Vol. 56, 2013. 4 pages, 1 figure. More detailed results about this project can be found in the accompanying paper: arXiv 1308.2727

The mass assembly of a whole population of sub-Milky Way galaxies is studied by means of hydrodynamical simulations within the $\Lambda$-CDM cosmology. Our results show that while dark halos assemble hierarchically, in stellar mass this trend is inverted in the sense that the smaller the galaxy, the later is its stellar mass assembly on average. Our star formation and supernovae feedback implementation in a multi-phase interstellar medium seems to play a key role on this process. However, the obtained downsizing trend is not yet as strong as observations show., Comment: Accepted for publication in BAAA, Vol. 55, 2012. 4 pages, 1 figure, baaa.sty. More detailed results about this project can be found in the accompanying paper: arXiv 1308.2727

Cosmological hydrodynamical simulations are studied in order to analyse generic trends for the stellar, baryonic and halo mass assembly of low-mass galaxies (M_* < 3 x 10^10 M_sun) as a function of their present halo mass, in the context of the Lambda-CDM scenario and common subgrid physics schemes. We obtain that smaller galaxies exhibit higher specific star formation rates and higher gas fractions. Although these trends are in rough agreement with observations, the absolute values of these quantities tend to be lower than observed ones since z~2. The simulated galaxy stellar mass fraction increases with halo mass, consistently with semi-empirical inferences. However, the predicted correlation between them shows negligible variations up to high z, while these inferences seem to indicate some evolution. The hot gas mass in z=0 halos is higher than the central galaxy mass by a factor of ~1-1.5 and this factor increases up to ~5-7 at z~2 for the smallest galaxies. The stellar, baryonic and halo evolutionary tracks of simulated galaxies show that smaller galaxies tend to delay their baryonic and stellar mass assembly with respect to the halo one. The Supernova feedback treatment included in this model plays a key role on this behaviour albeit the trend is still weaker than the one inferred from observations. At z>2, the overall properties of simulated galaxies are not in large disagreement with those derived from observations., Comment: 19 pages, 12 figures. Accepted for publication in MNRAS: 6th August 2013. First submitted: 7th July 2012

We studied the evolution of the gas kinematics of galaxies by performing hydrodynamical simulations in a cosmological scenario. We paid special attention to the origin of the scatter of the Tully-Fisher relation and the features which could be associated with mergers and interactions. We extended the study by De Rossi et al. (2010) and analysed their whole simulated sample which includes both, gas disc-dominated and spheroid-dominated systems. We found that mergers and interactions can affect the rotation curves directly or indirectly inducing a scatter in the Tully-Fisher Relation larger than the simulated evolution since z=3. In agreement with previous works, kinematical indicators which combine the rotation velocity and dispersion velocity in their definitions lead to a tighter relation. In addition, when we estimated the rotation velocity at the maximum of the rotation curve, we obtained the best proxy for the potential well regardless of morphology., Comment: Accepted for publication in BAAA, Vol. 55, 2012. 4 pages, 1 figure, baaa.sty

Recent observational and theoretical works have suggested that the Tully-Fisher Relation might be generalised to include dispersion-dominated systems by combining the rotation and dispersion velocity in the definition of the kinematical indicator. Mergers and interactions have been pointed out as responsible of driving turbulent and disordered gas kinematics, which could generate Tully-Fisher Relation outliers. We intend to investigate the gas kinematics of galaxies by using a simulated sample which includes both, gas disc-dominated and spheroid-dominated systems. Cosmological hydrodynamical simulations which include a multiphase model and physically-motivated Supernova feedback were performed in order to follow the evolution of galaxies as they are assembled. Both the baryonic and stellar Tully-Fisher relations for gas disc-dominated systems are tight while, as more dispersion-dominated systems are included, the scatter increases. We found a clear correlation between $\sigma / V_{\rm rot}$ and morphology, with dispersion-dominated systems exhibiting the larger values ($> 0.7$). Mergers and interactions can affect the rotation curves directly or indirectly inducing a scatter in the Tully-Fisher Relation larger than the simulated evolution since $z \sim 3$. Kinematical indicators which combine rotation velocity and dispersion velocity can reduce the scatter in the baryonic and the stellar mass-velocity relations. Our findings also show that the lowest scatter in both relations is obtained if the velocity indicators are measured at the maximum of the rotation curve. Moreover, the rotation velocity estimated at the maximum of the gas rotation curve is found to be the best proxy for the potential well regardless of morphology., Comment: 16 pages, 10 figures, accepted for publication in A&A

In this work, we studied the stellar and baryonic Tully-Fisher relations by using hydrodynamical simulations in a cosmological framework. We found that supernova feedback plays an important role on shaping the stellar Tully-Fisher relation causing a steepening of its slope at the low-mass end, consistently with observations. The bend of the relation occurs at a characteristic velocity of approximately 100 km/s, in concordance with previous observational and theoretical findings. With respect to the baryonic Tully-Fisher relation, the model predicts a linear trend at z~0 with a weaker tendency for a bend at higher redshifts. In our simulations, this behaviour is a consequence of the more efficient action of supernova feedback at regulating the star formation process in smaller galaxies., Comment: Accepted for publication in BAAA, Vol. 53, 2010. 4 pages, baaa2010.sty

We have studied the origin of the Tully-Fisher relation by analysing hydrodynamical simulations in a Lambda-CDM universe. We found that smaller galaxies exhibit lower stellar masses than those predicted by the linear fit to high mass galaxies (fast rotators), consistently with observations. In this model, these trends are generated by the more efficient action of Supernova feedback in the regulation of the star formation in smaller galaxies. Without introducing scale-dependent parameters, the model predicts that the Tully-Fisher relation bends at a characteristic velocity of around 100 km/s, in agreement with previous observational and theoretical findings., Comment: 2 pages, to appear in the proceedings of the XIII Latin American Regional IAU Meeting, Morelia, Michoacan, Mexico, November 8-12, 2010

Recent observational results found a bend in the Tully-Fisher Relation in such a way that low mass systems lay below the linear relation described by more massive galaxies. We intend to investigate the origin of the observed features in the stellar and baryonic Tully-Fisher relations and analyse the role played by galactic outflows on their determination. Cosmological hydrodynamical simulations which include Supernova feedback were performed in order to follow the dynamical evolution of galaxies. We found that Supernova feedback is a fundamental process in order to reproduce the observed trends in the stellar Tully-Fisher relation. Simulated slow rotating systems tend to have lower stellar masses than those predicted by the linear fit to the massive end of the relation, consistently with observations. This feature is not present if Supernova feedback is turned off. In the case of the baryonic Tully-Fisher relation, we also detect a weaker tendency for smaller systems to lie below the linear relation described by larger ones. This behaviour arises as a result of the more efficient action of Supernovae in the regulation of the star formation process and in the triggering of powerful galactic outflows in shallower potential wells which may heat up and/or expel part of the gas reservoir., Comment: 10 pages, 9 figures, accepted for publication in A&A

We have resimulated the six galaxy-sized haloes of the Aquarius Project including metal-dependent cooling, star formation and supernova feedback. This allows us to study not only how dark matter haloes respond to galaxy formation, but also how this response is affected by details of halo assembly history. In agreement with previous work, we find baryon condensation to lead to increased dark matter concentration. Dark matter density profiles differ substantially in shape from halo to halo when baryons are included, but in all cases the velocity dispersion decreases monotonically with radius. Some haloes show an approximately constant dark matter velocity anisotropy with $ \beta \approx 0.1-02$, while others retain the anisotropy structure of their baryon-free versions. Most of our haloes become approximately oblate in their inner regions, although a few retain the shape of their dissipationless counterparts. Pseudo-phase-space densities are described by a power law in radius of altered slope when baryons are included. The shape and concentration of the dark matter density profiles are not well reproduced by published adiabatic contraction models. The significant spread we find in the density and kinematic structure of our haloes appears related to differences in their formation histories. Such differences already affect the final structure in baryon-free simulations, but they are reinforced by the inclusion of baryons, and new features are produced. The details of galaxy formation need to be better understood before the inner dark matter structure of galaxies can be used to constrain cosmological models or the nature of dark matter., Comment: 14 pages, 9 figures. Accepted MNRAS. Revised version includes discussion on resolution effects and minor changes.

We have studied the dark matter (DM) distribution in a approx 10^12 h^-1 M_sun mass halo extracted from a simulation consistent with the concordance cosmology, where the physics regulating the transformation of gas into stars was allowed to change producing galaxies with different morphologies. The presence of baryons produces the concentration of the DM halo with respect to its corresponding dissipationless run, but we found that this response does not only depend on the amount of baryons gathered in the central region but also on the way they have been assembled. DM and baryons affect each other in a complex way so the formation history of a galaxy plays an important role on its final total mass distribution. Supernova (SN) feedback regulates the star formation and triggers galactic outflows not only in the central galaxy but also in its satellites. Our results suggest that, as the effects of SN feedback get stronger, satellites get less massive and can even be more easily disrupted by dynamical friction, transferring less angular momentum. We found indications that this angular momentum could be acquired not only by the outer part of the DM halo but also by the inner ones and by the stellar component in the central galaxy. The latter effect produces stellar migration which contributes to change the inner potential well, probably working against further DM contraction. As a consequence of the action of these processes, when the halo hosts a galaxy with an important disc structure formed by smooth gas accretion, it is more concentrated than when it hosts a spheroidal system which experienced more massive mergers and interactions. (abridged), Comment: 16 pages, 16 figures, accepted to MNRAS

We analyse the dark matter (DM) distribution in a approx 10^12 M_sun halo extracted from a simulation consistent with the concordance cosmology, where the physics regulating the transformation of gas into stars was allowed to change producing galaxies with different morphologies. Although the DM profiles get more concentrated as baryons are collected at the centre of the haloes compared to a pure dynamical run, the total baryonic mass alone is not enough to fully predict the reaction of the DM profile. We also note that baryons affect the DM distribution even outside the central regions. Those systems where the transformation of gas into stars is regulated by Supernova (SN) feedback, so that significant disc structures are able to form, are found to have more concentrated dark matter profiles than a galaxy which has efficiently transformed most of its baryons into stars at early times. The accretion of satellites is found to be associated with an expansion of the dark matter profiles, triggered by angular momentum transfer from the incoming satellites. As the impact of SN feedback increases, the satellites get less massive and are even strongly disrupted before getting close to the main structure causing less angular momentum transfer. Our findings suggest that the response of the DM halo is driven by the history of assembly of baryons into a galaxy along their merger tree., Comment: 5 pages, 5 figures, accepted for publication in MNRAS

We use cosmological simulations to study the effects of supernova (SN) feedback on the dark matter distribution in galaxies. We simulate the formation of a Milky-Way type galaxy using a version of the SPH code GADGET2 which includes chemical enrichment and energy feedback by SN, a multiphase model for the gas component and metal-dependent cooling. We analyse the impact of the main three input SN feedback parameters on the amplitude and shape of the dark matter density profiles, focusing on the inner regions of the halo. In order to test the dependence of the results on the halo mass, we simulated a scale-down version of this system. First results of this ongoing work show that the dark matter distribution is affected by the feedback, through the redistribution of the baryons. Our findings suggest that the response of the dark matter halo could be the result of a combination of several physical parameters such as the amount of stellar mass formed at the centre, its shape, and probably the bursty characteristics of the star formation rate. As expected, we find that the dark matter haloes of small galaxies are more sensitive to SN feedback. Higher resolution simulations are being performed to test for numerical effects., Comment: 6 pages, 5 figures, Proceedings of IAU Symposium 254 "The Galaxy disk in a cosmological context", Copenhagen, June 2008