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.

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.

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

Despite the insights gained in the last few years, our knowledge about the formation and evolution scenario for the spheroid-dominated galaxies is still incomplete. New and more powerful cosmological simulations have been developed that together with more precise observations open the possibility of more detailed study of the formation of early-type galaxies (ETGs). The aim of this work is to analyse the assembly histories of ETGs in a $\Lambda$-CDM cosmology, focussing on the archeological approach given by the mass-growth histories.We inspected a sample of dispersion-dominated galaxies selected from the largest volume simulation of the EAGLE project. This simulation includes a variety of physical processes such as radiative cooling, star formation (SF), metal enrichment, and stellar and active galactic nucleus (AGN) feedback. The selected sample comprised 508 spheroid-dominated galaxies classified according to their dynamical properties. Their surface brightness profile, the fundamental relations, kinematic properties, and stellar-mass growth histories are estimated and analysed. The findings are confronted with recent observations.The simulated ETGs are found to globally reproduce the fundamental relations of ellipticals. All of them have an inner disc component where residual younger stellar populations (SPs) are detected. A fraction of this inner-disc correlates with bulge-to-total ratio. We find a relation between kinematics and shape that implies that dispersion-dominated galaxies with low $V/\sigma_L$ (where $V$ is the average rotational velocity and $\sigma_L$ the one dimensional velocity dispersion) tend to have ellipticity smaller than $\sim 0.5$ and are dominated by old stars. Abridged

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.

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.

Understanding the formation and evolution of early-type, spheroid-dominated galaxies is an open question within the context of the hierarchical clustering scenario, particularly, in low-density environments. Our goal is to study the main structural, dynamical, and stellar population properties and assembly histories of field spheroid-dominated galaxies formed in a LCDM scenario to assess to what extend they are consistent with observations. We selected spheroid-dominated systems from a LCDM simulation that includes star formation, chemical evolution and Supernova feedback. A sample of 18 field systems with Mstar <= 6x10^10 Msun that are dominated by the spheroid component. For this sample we estimate the fundamental relations of ellipticals and then compared with current observations. The simulated spheroid galaxies have sizes in good agreement with observations. The bulges follow a Sersic law with Sersic indexes that correlate with the bulge-to-total mass ratios. The structural-dynamical properties of the simulated galaxies are consistent with observed Faber-Jackson, Fundamental Plane, and Tully-Fisher relations. However, the simulated galaxies are bluer and with higher star formation rates than observed isolated early-type galaxies. The archaeological mass growth histories show a slightly delayed formation and more prominent inside-out growth mode than observational inferences based on the fossil record method. The main structural and dynamical properties of the simulated spheroid-dominated galaxies are consistent with observations. This is remarkable since none of them has been tuned to be reproduced. However, the simulated galaxies are blue and star-forming, and with later stellar mass growth histories as compared to observational inferences. This is mainly due to the persistence of extended discs in the simulations. Abridged

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.

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]

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.

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]

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.

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.

High-mass X-ray binaries (HMXBs) might have contributed a non-negligible fraction of the energy feedback to the interstellar and intergalactic media at high redshift, becoming important sources for the heating and ionization history of the Universe. However, the importance of this contribution depends on the hypothesized increase in the number of HMXBs formed in low-metallicity galaxies and in their luminosities. In this work we test the aforementioned hypothesis, and quantify the metallicity dependence of HMXB population properties. We compile from the literature a large set of data on the sizes and X-ray luminosities of HMXB populations in nearby galaxies with known metallicities and star formation rates. We use Bayesian inference to fit simple Monte Carlo models that describe the metallicity dependence of the size and luminosity of the HMXB populations. We find that HMXBs are typically ten times more numerous per unit star formation rate in low-metallicity galaxies (12 + log(O/H) < 8, namely < 20% solar) than in solar-metallicity galaxies. The metallicity dependence of the luminosity of HMXBs is small compared to that of the population size. Our results support the hypothesis that HMXBs are more numerous in low-metallicity galaxies, implying the need to investigate the feedback in the form of X-rays and energetic mass outflows of these high-energy sources during cosmic dawn.

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).

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.

The standard disc formation scenario postulates that disc forms as the gas cools and flows into the centre of the dark matter halo, conserving the specific angular momentum. Major mergers have been shown to be able to destroy or highly perturb the disc components. More recently, the alignment of the material that is accreted to form the galaxy has been pointed out as a key ingredient to determine galaxy morphology. However, in a hierarchical scenario galaxy formation is a complex process that combines these processes and others in a non-linear way so that the origin of galaxy morphology remains to be fully understood. We aim at exploring the differences in the formation histories of galaxies with a variety of morphology, but quite recent merger histories, to identify which mechanisms are playing a major role. We analyse when minor mergers can be considered relevant to determine galaxy morphology. We also study the specific angular momentum content of the disc and central spheroidal components separately. We used cosmological hydrodynamical simulations that include an effective, physically motivated supernova feedback that is able to regulate the star formation in haloes of different masses. We analysed the morphology and formation history of a sample of 15 galaxies of a cosmological simulation. We performed a spheroid-disc decomposition of the selected galaxies and their progenitor systems. The angular momentum orientation of the merging systems as well as their relative masses were estimated to analyse the role played by orientation and by minor mergers in the determination of the morphology. We found the discs to be formed by conserving the specific angular momentum in accordance with the classical disc formation model. The specific angular momentum of the stellar central spheroid correlates with the dark matter halo angular momentum and determines a power law. Abridged

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.

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.

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.

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.

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.

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.

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.

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.

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)

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.

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.

Rotation curves of interacting galaxies often show that velocities are either rising or falling in the direction of the companion galaxy. We seek to reproduce and analyse these features in the rotation curves of simulated equal-mass galaxies suffering a one-to-one encounter, as possible indicators of close encounters. Using simulations of major mergers in 3D, we study the time evolution of these asymmetries in a pair of galaxies, during the first passage. Our main results are: (a) the rotation curve asymmetries appear right at pericentre of the first passage, (b) the significant disturbed rotation velocities occur within a small time interval, of ~ 0.5 Gyr h^-1, and therefore the presence of bifurcation in the velocity curve could be used as an indicator of the pericentre occurrence. These results are in qualitative agreement with previous findings for minor mergers and fly-byes.