Search

We study the effect of mergers on the morphology of galaxies by means of the simulated merger tree approach first proposed by Moster et al. This method combines N-body cosmological simulations and semi-analytic techniques to extract realistic initial conditions for galaxy mergers. These are then evolved using high resolution hydrodynamical simulations, which include dark matter, stars, cold gas in the disc and hot gas in the halo. We show that the satellite mass accretion is not as effective as previously thought, as there is substantial stellar stripping before the final merger. The fraction of stellar disc mass transferred to the bulge is quite low, even in the case of a major merger, mainly due to the dispersion of part of the stellar disc mass into the halo. We confirm the findings of Hopkins et al., that a gas rich disc is able to survive major mergers more efficiently. The enhanced star formation associated with the merger is not localised to the bulge of galaxy, but a substantial fraction takes place in the disc too. The inclusion of the hot gas reservoir in the galaxy model contributes to reducing the efficiency of bulge formation. Overall, our findings suggest that mergers are not as efficient as previously thought in transforming discs into bulges. This possibly alleviates some of the tensions between observations of bulgeless galaxies and the hierarchical scenario for structure formation., Comment: MNRAS Accepted, 17 pages, 11 figures, 3 Tables

We use the NIHAO simulations to investigate the effects of baryonic physics on the time evolution of Dark Matter central density profiles. The sample is made of $\approx 70$ independent high resolution hydrodynamical simulations of galaxy formation and covers a wide mass range: 1e10< Mhalo <1e12, i.e., from dwarfs to L* . We confirm previous results on the dependence of the inner dark matter density slope, $\alpha$, on the ratio between stellar-to-halo mass. We show that this relation holds approximately at all redshifts (with an intrinsic scatter of ~0.18 in $\alpha$). This implies that in practically all haloes the shape of their inner density profile changes quite substantially over cosmic time, as they grow in stellar and total mass. Thus, depending on their final stellar-to-halo mass ratio, haloes can either form and keep a substantial density core (size~1 kpc), or form and then destroy the core and re-contract the halo, going back to a cuspy profile, which is even steeper than CDM predictions for massive galaxies (~1e12 Msun). We show that results from the NIHAO suite are in good agreement with recent observational measurements of $\alpha$ in dwarf galaxies. Overall our results suggest that the notion of a universal density profile for dark matter haloes is no longer valid in the presence of galaxy formation., Comment: 11 pages, 13 figures. Corrected typo in table 2 (middle row) with respect to the version published in MNRAS

We use semi-analytic models and cosmological merger trees to provide the initial conditions for multi-merger numerical hydrodynamic simulations, and exploit these simulations to explore the effect of galaxy interaction and merging on star formation (SF). We compute numerical realisations of twelve merger trees from z=1.5 to z=0. We include the effects of the large hot gaseous halo around all galaxies, following recent obervations and predictions of galaxy formation models. We find that including the hot gaseous halo has a number of important effects. Firstly, as expected, the star formation rate on long timescales is increased due to cooling of the hot halo and refuelling of the cold gas reservoir. Secondly, we find that interactions do not always increase the SF in the long term. This is partially due to the orbiting galaxies transferring gravitational energy to the hot gaseous haloes and raising their temperature. Finally we find that the relative size of the starburst, when including the hot halo, is much smaller than previous studies showed. Our simulations also show that the order and timing of interactions are important for the evolution of a galaxy. When multiple galaxies interact at the same time, the SF enhancement is less than when galaxies interact in series. All these effects show the importance of including hot gas and cosmologically motivated merger trees in galaxy evolution models., Comment: 19 pages, 15 figures, 6 tables. Accepted for publication in MNRAS

We show that the cool gas masses of galactic discs reach a steady state that lasts many Gyr after their last major merger in cosmological hydrodynamic simulations. The mass of disc gas, M$_{\rm gas}$, depends upon a galaxy halo's spin and virial mass, but not upon stellar feedback. Halos with low spin have high star formation efficiency and lower disc gas mass. Similarly, lower stellar feedback leads to more star formation so the gas mass ends up nearly the same irregardless of stellar feedback strength. Even considering spin, the M$_{\rm gas}$ relation with halo mass, M$_{200}$ only shows a factor of 3 scatter. The M$_{\rm gas}$--M$_{200}$ relation show a break at M$_{200}$=$2\times10^{11}$ M$_\odot$ that corresponds to an observed break in the M$_{\rm gas}$--M$_\star$ relation. The constant disc mass stems from a shared halo gas density profile in all the simulated galaxies. In their outer regions, the profiles are isothermal. Where the profile rises above $n=10^{-3}$ cm$^{-3}$, the gas readily cools and the profile steepens. Inside the disc, rotation supports gas with a flatter density profile except where supernova explosions disrupt the disc. Energy injection from stellar feedback also provides pressure support to the halo gas to prevent runaway cooling flows. The resulting constant gas mass makes simpler models for galaxy formation possible, either using a "bathtub" model for star formation rates or when modeling chemical evolution., Comment: To be submitted to MNRAS. Comments very welcome in the meantime

We use cosmological hydrodynamical zoom-in simulations with the SPH code gasoline of four haloes of mass M_{200} \sim 10^{13}\Msun to study the response of the dark matter to elliptical galaxy formation. Our simulations include metallicity dependent gas cooling, star formation, and feedback from massive stars and supernovae, but not active galactic nuclei (AGN). At z=2 the progenitor galaxies have stellar to halo mass ratios consistent with halo abundance matching, assuming a Salpeter initial mass function. However by z=0 the standard runs suffer from the well known overcooling problem, overpredicting the stellar masses by a factor of > 4. To mimic a suppressive halo quenching scenario, in our forced quenching (FQ) simulations, cooling and star formation are switched off at z=2. The resulting z=0 galaxies have stellar masses, sizes and circular velocities close to what is observed. Relative to the control simulations, the dark matter haloes in the FQ simulations have contracted, with central dark matter density slopes d\log\rho/d\log r \sim -1.5, showing that dry merging alone is unable to fully reverse the contraction that occurs at z>2. Simulations in the literature with AGN feedback however, have found expansion or no net change in the dark matter halo. Thus the response of the dark matter halo to galaxy formation may provide a new test to distinguish between ejective and suppressive quenching mechanisms., Comment: 21 pages, 18 figures, accepted to MNRAS

In a Universe where AGN feedback regulates star formation in massive galaxies, a strong correlation between these two quantities is expected. If the gas causing star formation is also responsible for feeding the central black hole, then a positive correlation is expected. If powerful AGNs are responsible for the star formation quenching, then a negative correlation is expected. Observations so far have mainly found a mild correlation or no correlation at all (i.e. a flat relation between star formation rate (SFR) and AGN luminosity), raising questions about the whole paradigm of "AGN feedback". In this paper, we report the predictions of the GALFORM semi-analytical model, which has a very strong coupling between AGN activity and quenching of star formation. The predicted SFR-AGN luminosity correlation appears negative in the low AGN luminosity regime, where AGN feedback acts, but becomes strongly positive in the regime of the brightest AGN. Our predictions reproduce reasonably well recent observations by Rosario et al., yet there is some discrepancy in the normalisation of the correlation at low luminosities and high redshifts. Though this regime could be strongly influenced by observational biases, we argue that the disagreement could be ascribed to the fact that GALFORM neglects AGN variability effects. Interestingly, the galaxies that dominate the regime where the observations imply a weak correlation are massive early-type galaxies that are subject to AGN feedback. Nevertheless, these galaxies retain high enough molecular hydrogen contents to maintain relatively high star formation rates and strong infrared emission., Comment: 9 pages, 6 figures, accepted for publication to MNRAS

We study how X-rays from stellar binary systems and the hot intracluster medium (ICM) affect the radiative cooling rates of gas in galaxies. Our study uses a novel implementation of gas cooling in the moving-mesh hydrodynamics code \textsc{arepo}. X-rays from stellar binaries do not affect cooling at all as their emission spectrum is too hard to effectively couple with galactic gas. In contrast, X-rays from the ICM couple well with gas in the temperature range $10^4 - 10^6$ K. Idealised simulations show that the hot halo radiation field has minimal impact on the dynamics of cooling flows in clusters because of the high virial temperature ($> 10^7$K), making the interaction between the gas and incident photons very ineffective. Satellite galaxies in cluster environments, on the other hand, experience a high radiation flux due to the emission from the host halo. Low mass satellites ($< 10^{12}\rm{M_\odot}$) in particular have virial temperatures that are exactly in the regime where the effect of the radiation field is maximal. Idealised simulations of satellite galaxies including only the effect of host halo radiation (no ram pressure stripping or tidal effects) fields show a drastic reduction in the amount of cool gas formed ($\sim 40\%$) on a short timescale of about $0.5$ Gyrs. A galaxy merger simulation including all the other environmental quenching mechanisms, shows about $20\%$ reduction in the stellar mass of the satellite and about $\sim 30\%$ reduction in star formation rate after $1$ Gyr due to the host hot halo radiation field. These results indicate that the hot halo radiation fields potentially play an important role in quenching galaxies in cluster environments., Comment: 14 pages, 12 figures, submitted to MNRAS

We present the first numerical simulations in coupled dark energy cosmologies with high enough resolution to investigate the effects of the coupling on galactic and sub-galactic scales. We choose two constant couplings and a time-varying coupling function and we run simulations of three Milky-Way-size halos ($\sim$10$^{12}$M$_{\odot}$), a lower mass halo (6$\times$10$^{11}$M$_{\odot}$) and a dwarf galaxy halo (5$\times$10$^{9}$M$_{\odot}$). We resolve each halo with several millions dark matter particles. On all scales the coupling causes lower halo concentrations and a reduced number of substructures with respect to LCDM. We show that the reduced concentrations are not due to different formation times, but they are related to the extra terms that appear in the equations describing the gravitational dynamics. On the scale of the Milky Way satellites, we show that the lower concentrations can help in reconciling observed and simulated rotation curves, but the coupling values necessary to have a significant difference from LCDM are outside the current observational constraints. On the other hand, if other modifications to the standard model allowing a higher coupling (e.g. massive neutrinos) are considered, coupled dark energy can become an interesting scenario to alleviate the small-scale issues of the LCDM model., Comment: 12 pages, 16 figures, submitted to MNRAS

Recently it has been shown that a large fraction of the dwarf satellite galaxies orbiting the Andromeda galaxy are surprisingly aligned in a thin, extended and kinematically coherent planar structure. The presence of such a structure seems to challenge the current Cold Dark Matter paradigm of structure formation, which predicts a more uniform distribution of satellites around central objects. We show that it is possible to obtain a thin, extended, rotating plane of satellites resembling the one in Andromeda in cosmological collisionless simulations based on the Cold Dark Matter model. Our new high resolution simulations show a correlation between the formation time of the dark matter halo and the thickness of the plane of satellites. Our simulations have a high incidence of satellite planes as thin, extended, and as rich as the one in Andromeda and with a very coherent kinematic structure when we select high concentration/early forming halos. By tracking the formation of the satellites in the plane we show that they have been mainly accreted onto the main object along thin dark matter filaments at high redshift. Our results show that the presence of a thin, extended, rotating plane of satellites is not a challenge for the Cold Dark Matter paradigm, but actually supports one of the predictions of this paradigm related to the presence of filaments of dark matter around galaxies at high redshift., Comment: 6 pages, 5 figures, accepted for publication in ApJ Letters, minor changes to match the accepted version

In this first paper we discuss the linear theory and the background evolution of a new class of models we dub SCDEW: Strongly Coupled DE, plus WDM. In these models, WDM dominates today's matter density; like baryons, WDM is uncoupled. Dark Energy is a scalar field $\Phi$; its coupling to ancillary CDM, whose today's density is $\ll 1\, \%$, is an essential model feature. Such coupling, in fact, allows the formation of cosmic structures, in spite of very low WDM particle masses ($\sim 100$ eV). SCDEW models yields Cosmic Microwave Background and linear Large Scale features substantially undistinguishable from $\Lambda$CDM, but thanks to the very low WDM masses they strongly alleviate $\Lambda$CDM issues on small scales, as confirmed via numerical simulations in the II associated paper. Moreover SCDEW cosmologies significantly ease the coincidence and fine tuning problems of $\Lambda$CDM and, by using a field theory approach, we also outline possible links with inflationary models. We also discuss a possible fading of the coupling at low redshifts which prevents non linearities on the CDM component to cause computational problems. The (possible) low-$z$ coupling suppression, its mechanism, and its consequences are however still open questions -not necessarily problems- for SCDEW models. The coupling intensity and the WDM particle mass, although being extra parameters in respect to $\Lambda$CDM, are found to be substantially constrained a priori so that, if SCDEW is the underlying cosmology, we expect most data to fit also $\Lambda$CDM predictions., Comment: 11 pages, 9 figures , accepted for publication on MNRAS; updated to match the published version; the companion paper can be found here http://arxiv.org/abs/1503.07867

In this second paper we present the first Nbody cosmological simulations of strongly coupled Dark Energy models (SCDEW), a class of models that alleviates theoretical issues related to the nature of dark energy. SCDEW models assume a strong coupling between Dark Energy (DE) and an ancillary Cold Dark Matter (CDM) component together with the presence of an uncoupled Warm Dark Matter component. The strong coupling between CDM and DE allows us to preserve small scale fluctuations even if the warm particle is quite light ($\approx 100$ eV). Our large scale simulations show that, for $10^{11}

We use the NIHAO (Numerical Investigation of Hundred Astrophysical Objects) cosmological simulations to study the effects of galaxy formation on key properties of dark matter (DM) haloes. NIHAO consists of $\simeq 90$ high-resolution SPH simulations that include (metal-line) cooling, star formation, and feedback from massive stars and SuperNovae, and cover a wide stellar and halo mass range: $10^6 < M_* / M_{\odot} < 10^{11}$ ( $10^{9.5} < M_{\rm halo} / M_{\odot} < 10^{12.5}$). When compared to DM-only simulations, the NIHAO haloes have similar shapes at the virial radius, R_{\rm vir}, but are substantially rounder inside $\simeq 0.1R_{\rm vir}$. In NIHAO simulations $c/a$ increases with halo mass and integrated star formation efficiency, reaching $\sim 0.8$ at the Milky Way mass (compared to 0.5 in DM-only), providing a plausible solution to the long-standing conflict between observations and DM-only simulations. The radial profile of the phase-space $Q$ parameter ($\rho/\sigma^3$) is best fit with a single power law in DM-only simulations, but shows a flattening within $\simeq 0.1R_{\rm vir}$ for NIHAO for total masses $M>10^{11} M_{\odot}$. Finally, the global velocity distribution of DM is similar in both DM-only and NIHAO simulations, but in the solar neighborhood, NIHAO galaxies deviate substantially from Maxwellian. The distribution is more symmetric, roughly Gaussian, with a peak that shifts to higher velocities for Milky Way mass haloes. We provide the distribution parameters which can be used for predictions for direct DM detection experiments. Our results underline the ability of the galaxy formation processes to modify the properties of dark matter haloes., Comment: 19 pages, 17 figures, analysis strongly improved, main conclusions unchanged, accepted for publication in MNRAS

We introduce project NIHAO (Numerical Investigation of a Hundred Astrophysical Objects), a set of 100 cosmological zoom-in hydrodynamical simulations performed using the gasoline code, with an improved implementation of the SPH algorithm. The haloes in our study range from dwarf to Milky Way masses, and represent an unbiased sampling of merger histories, concentrations and spin parameters. The particle masses and force softenings are chosen to resolve the mass profile to below 1% of the virial radius at all masses, ensuring that galaxy half-light radii are well resolved. Using the same treatment of star formation and stellar feedback for every object, the simulated galaxies reproduce the observed inefficiency of galaxy formation across cosmic time as expressed through the stellar mass vs halo mass relation, and the star formation rate vs stellar mass relation. We thus conclude that stellar feedback is the chief piece of physics required to limit the efficiency of star formation in galaxies less massive than the Milky Way., Comment: 14 pages, 8 figures, accepted for publication in MNRAS

We present the formation of a Kinematically Decoupled Core (KDC) in an elliptical galaxy, resulting from a major merger simulation of two disk galaxies. We show that although the two progenitor galaxies are initially following a prograde orbit, strong reactive forces during the merger can cause a short-lived change of their orbital spin; the two progenitors follow a retrograde orbit right before their final coalescence. This results in a central kinematic decoupling and the formation of a large-scale (~2 kpc radius) counter-rotating core (CRC) at the center of the final elliptical-like merger remnant (M*=1.3x10^11 Msun), while its outer parts keep the rotation direction of the initial orbital spin. The stellar velocity dispersion distribution of the merger remnant galaxy exhibits two symmetrical off-centered peaks, comparable to the observed "2-sigma galaxies". The KDC/CRC consists mainly of old, pre-merger population stars (older than 5 Gyr), remaining prominent in the center of the galaxy for more than 2 Gyr after the coalescence of its progenitors. Its properties are consistent with KDCs observed in massive elliptical galaxies. This new channel for the formation of KDCs from prograde mergers is in addition to previously known formation scenarios from retrograde mergers and can help towards explaining the substantial fraction of KDCs observed in early-type galaxies., Comment: 7 pages, 5 figures, accepted for publication in ApJ Letters

The radial density profiles of stellar galaxy discs can be well approximated as an exponential. Compared to this canonical form, however, the profiles in the majority of disc galaxies show downward or upward breaks at large radii. Currently, there is no coherent explanation in a galaxy formation context of the radial profile per se, along with the two types of profile breaks. Using a set of controlled hydrodynamic simulations of disc galaxy formation, we find a correlation between the host halo's initial angular momentum and the resulting radial profile of the stellar disc: galaxies that live in haloes with a low spin parameter {\lambda} <~ 0.03 show an up-bending break in their disc density profiles, while galaxies in haloes of higher angular momentum show a down-bending break. We find that the case of pure exponential profiles ({\lambda} ~ 0.035) coincides with the peak of the spin parameter distribution from cosmological simulations. Our simulations not only imply an explanation of the observed behaviours, but also suggest that the physical origin of this effect is related to the amount of radial redistribution of stellar mass, which is anti-correlated with {\lambda}., Comment: 6 pages, 3 figures, adapted to the published version in MNRAS-L

Models including an energy transfer from CDM to DE are widely considered in the literature, namely to allow DE a significant high-z density. Strongly Coupled cosmologies assume a much larger coupling between DE and CDM, together with the presence of an uncoupled warm DM component, as the role of CDM is mostly restricted to radiative eras. This allows us to preserve small scale fluctuations even if the warm particle, possibly a sterile neutrino, is quite light, O(100 eV). Linear theory and numerical simulations show that these cosmologies agree with LCDM on supergalactic scales; e.g., CMB spectra are substantially identical. Simultaneously, simulations show that they significantly ease problems related to the properties of MW satellites and cores in dwarfs. SC cosmologies also open new perspectives on early black hole formation, and possibly lead towards unificating DE and inflationary scalar fields., Comment: 4 pages, 6 figures, Proceeding of the "Neutrino Oscillation Workshop" NOW 2014, Conca Specchiulla, Otranto, Italy, 7-14 September 2014

Galaxies are not randomly distributed across the universe but showing different kinds of alignment on different scales. On small scales satellite galaxies have a tendency to distribute along the major axis of the central galaxy, with dependence on galaxy properties that both red satellites and centrals have stronger alignment than their blue counterparts. On large scales, it is found that the major axes of Luminous Red Galaxies (LRGs) have correlation up to 30Mpc/h. Using hydro-dynamical simulation with star formation, we investigate the origin of galaxy alignment on different scales. It is found that most red satellite galaxies stay in the inner region of dark matter halo inside which the shape of central galaxy is well aligned with the dark matter distribution. Red centrals have stronger alignment than blue ones as they live in massive haloes and the central galaxy-halo alignment increases with halo mass. On large scales, the alignment of LRGs is also from the galaxy-halo shape correlation, but with some extent of mis-alignment. The massive haloes have stronger alignment than haloes in filament which connect massive haloes. This is contrary to the naive expectation that cosmic filament is the cause of halo alignment., Comment: 4 pages, 3 figures, To appear in the proceedings of the IAU Symposium 308 "The Zeldovich Universe: Genesis and Growth of the Cosmic Web"

Observations have shown that the spatial distribution of satellite galaxies is not random, but rather is aligned with the major axes of central galaxies (CGs). The strength of the alignment is dependent on the properties of both the satellites and centrals. Theoretical studies using dissipationless N-body simulations are limited by their inability to directly predict the shape of CGs. Using hydrodynamical simulations including gas cooling, star formation, and feedback, we carry out a study of galaxy alignment and its dependence on the galaxy properties predicted directly from the simulations.We found that the observed alignment signal is well produced, as is the color dependence: red satellites and red centrals both show stronger alignments than their blue counterparts. The reason for the stronger alignment of red satellites is that most of them stay in the inner region of the dark matter halo where the shape of the CG better traces the dark matter distribution. The dependence of alignment on the color of CGs arises from the halo mass dependence, since the alignment between the shape of the central stellar component and the inner halo increases with halo mass. We also find that the alignment of satellites is most strongly dependent on their metallicity, suggesting that the metallicity of satellites, rather than color, is a better tracer of galaxy alignment on small scales. This could be tested in future observational studies., Comment: ApJ Letter, accepted. Four figures, no table. The resolution of Fig 1 was downgraded due to the limitation of file size. Updated to match the version in press

We introduce a mass dependent density profile to describe the distribution of dark matter within galaxies, which takes into account the stellar-to-halo mass dependence of the response of dark matter to baryonic processes. The study is based on the analysis of hydrodynamically simulated galaxies from dwarf to Milky Way mass, drawn from the MaGICC project, which have been shown to match a wide range of disk scaling relationships. We find that the best fit parameters of a generic double power-law density profile vary in a systematic manner that depends on the stellar-to-halo mass ratio of each galaxy. Thus, the quantity Mstar/Mhalo constrains the inner ($\gamma$) and outer ($\beta$) slopes of dark matter density, and the sharpness of transition between the slopes($\alpha$), reducing the number of free parameters of the model to two. Due to the tight relation between stellar mass and halo mass, either of these quantities is sufficient to describe the dark matter halo profile including the effects of baryons. The concentration of the haloes in the hydrodynamical simulations is consistent with N-body expectations up to Milky Way mass galaxies, at which mass the haloes become twice as concentrated as compared with pure dark matter runs. This mass dependent density profile can be directly applied to rotation curve data of observed galaxies and to semi analytic galaxy formation models as a significant improvement over the commonly used NFW profile., Comment: 10 pages, 4 figures. Accepted for publication in MNRAS

Observations have shown that the spatial distribution of satellite galaxies is not random, but aligned with the major axes of central galaxies. This alignment is dependent on galaxy properties, such that red satellites are more strongly aligned than blue satellites. Theoretical work done to interpret this phenomena has found that it is due to the non-spherical nature of dark matter halos. However, most studies over-predict the alignment signal under the assumption that the central galaxy shape follows the shape of the host halo. It is also not clear whether the color dependence of alignment is due to an assembly bias or an evolution effect. In this paper we study these problems using a cosmological N-body simulation. Subhalos are used to trace the positions of satellite galaxies. It is found that the shape of dark matter halos are mis-aligned at different radii. If the central galaxy shares the same shape as the inner host halo, then the alignment effect is weaker and agrees with observational data. However, it predicts almost no dependence of alignment on the color of satellite galaxies, though the late accreted subhalos show stronger alignment with the outer layer of the host halo than their early accreted counterparts. We find that this is due to the limitation of pure N-body simulations that satellites galaxies without associated subhalos ('orphan galaxies') are not resolved. These orphan (mostly red) satellites often reside in the inner region of host halos and should follow the shape of the host halo in the inner region., Comment: 12 pages, 11 figures, Published on ApJ

We present the evolution of the structure of relaxed cold dark matter haloes in the cosmology from the Planck satellite. Our simulations cover 5 decades in halo mass, from dwarf galaxies to galaxy clusters. Due to the increased matter density and power spectrum normalization the concentration mass relation in the Planck cosmology has a 20 percent higher normalization at redshift z=0 compared to WMAP cosmology. We confirm that CDM haloes are better described by the Einasto profile; for example, at scales near galaxy half-light radii CDM haloes have significantly steeper density profiles than implied by NFW fits. There is a scatter of 0.2 dex in the Einasto shape parameter at fixed halo mass, adding further to the diversity of CDM halo profiles. The evolution of the concentration mass relation in our simulations is not reproduced by any of the analytic models in the literature. We thus provide a simple fitting formula that accurately describes the evolution between redshifts z=5 to z=0 for both NFW and Einasto fits. Finally, the observed concentrations and halo masses of spiral galaxies, groups and clusters of galaxies at low redshifts are in good agreement with our simulations, suggesting only mild halo response to galaxy formation on these scales., Comment: 18 pages, 18 figures, accepted to MNRAS, minor changes to previous version

The annihilation of dark matter particles in the halo of galaxies may end up into gamma rays, which travel almost unperturbed till to their detection at Earth. This annihilation signal can exhibit an anisotropic behavior quantified by the angular power spectrum, whose properties strongly depend on the dark matter distribution and its clumpiness. We use high resolution pure dark matter N-body simulations to quantify the contribution of different components (main halo and satellites) to the global signal as a function of the analytical profile adopted to describe the numerical results. We find that the smooth main halo dominates the angular power spectrum of the gamma-ray signal up to quite large multipoles, where the sub-haloes anisotropy signal starts to emerge, but the transition multipole strongly depends on the assumed radial profile. The extrapolation down to radii not resolved by current numerical simulations can affect both the normalization and the shape of the gamma-ray angular power spectrum. For the sub-haloes described by an asymptotically cored dark matter distribution, the angular power spectrum shows an overall smaller normalization and a flattening at high multipoles. Our results show the criticality of the dark matter density profile shape in gamma-ray anisotropy searches, and evaluate quantitatively the intrinsic errors occurring when extrapolating the dark matter radial profiles down to spatial scales not yet explored by numerical simulations., Comment: 7 pages, 8 figures. It matches the version published in MNRAS

We present the Dark MaGICC project, which aims to investigate the effect of Dark Energy (DE) modeling on galaxy formation via hydrodynamical cosmological simulations. Dark MaGICC includes four dynamical Dark Energy scenarios with time varying equations of state, one with a self-interacting Ratra-Peebles model. In each scenario we simulate three galaxies with high resolution using smoothed particle hydrodynamics (SPH). The baryonic physics model is the same used in the Making Galaxies in a Cosmological Context (MaGICC) project, and we varied only the background cosmology. We find that the Dark Energy parameterization has a surprisingly important impact on galaxy evolution and on structural properties of galaxies at z=0, in striking contrast with predictions from pure Nbody simulations. The different background evolutions can (depending on the behavior of the DE equation of state) either enhance or quench star formation with respect to a LCDM model, at a level similar to the variation of the stellar feedback parameterization, with strong effects on the final galaxy rotation curves. While overall stellar feedback is still the driving force in shaping galaxies, we show that the effect of the Dark Energy parameterization plays a larger role than previously thought, especially at lower redshifts. For this reason, the influence of Dark Energy parametrization on galaxy formation must be taken into account, especially in the era of precision cosmology., Comment: 11 pages, 13 figures

We use high-resolution Eulerian simulations to study the stability of cold gas flows in a galaxy size dark matter halo (10^12 Msun) at redshift z=2. Our simulations show that a cold stream penetrating a hot gaseous halo is stable against thermal convection and Kelvin-Helmholtz instability. We then investigate the effect of a satellite orbiting the main halo in the plane of the stream. The satellite is able to perturb the stream and to inhibit cold gas accretion towards the center of the halo for 0.5 Gyr. However, if the supply of cold gas at large distances is kept constant, the cold stream is able to re-establish itself after 0.3 Gyr. We conclude that cold streams are very stable against a large variety of internal and external perturbations., Comment: 5 pages, 3 figures, accepted for publication in MNRAS Letters

We investigate the interplay between jets from Active Galactic Nuclei (AGNs) and the surrounding InterStellar Medium (ISM) through full 3D, high resolution, Adaptive Mesh Refinement simulations performed with the FLASH code. We follow the jet- ISM system for several Myr in its transition from an early, compact source to an extended one including a large cocoon. During the jet evolution, we identify three major evolutionary stages and we find that, contrary to the prediction of popular theoretical models, none of the simulations shows a self-similar behavior. We also follow the evolution of the energy budget, and find that the fraction of input power deposited into the ISM (the AGN coupling constant) is of order of a few percent during the first few Myr. This is in broad agreement with galaxy formation models employing AGN feedback. However, we find that in these early stages, this energy is deposited only in a small fraction (< 1%) of the total ISM volume. Finally we demonstrate the relevance of backflows arising within the extended cocoon generated by a relativistic AGN jet within the ISM of its host galaxy, previously proposed as a mechanism for self-regulating the gas accretion onto the central object. These backflows tend later to be destabilized by the 3D dynamics, rather than by hydrodynamic (Kelvin- Helmholtz) instabilities. Yet, in the first few hundred thousand years, backflows may create a central accretion region of significant extent, and convey there as much as a few millions of solar masses., Comment: Accepted in MNRAS - 16 pages, 12 figures - Multimedia available on the author's webpage: http://www.mpia.de/~cielo

We present a first study of the effect of local photoionising radiation on gas cooling in smoothed particle hydrodynamics simulations of galaxy formation. We explore the combined effect of ionising radiation from young and old stellar populations. The method computes the effect of multiple radiative sources using the same tree algorithm used for gravity, so it is computationally efficient and well resolved. The method foregoes calculating absorption and scattering in favour of a constant escape fraction for young stars to keep the calculation efficient enough to simulate the entire evolution of a galaxy in a cosmological context to the present day. This allows us to quantify the effect of the local photoionisation feedback through the whole history of a galaxy`s formation. The simulation of a Milky Way like galaxy using the local photoionisation model forms ~ 40 % less stars than a simulation that only includes a standard uniform background UV field. The local photoionisation model decreases star formation by increasing the cooling time of the gas in the halo and increasing the equilibrium temperature of dense gas in the disc. Coupling the local radiation field to gas cooling from the halo provides a preventive feedback mechanism which keeps the central disc light and produces slowly rising rotation curves without resorting to extreme feedback mechanisms. These preliminary results indicate that the effect of local photoionising sources is significant and should not be ignored in models of galaxy formation., Comment: Accepted for Publication in MNRAS, 13 pages, 13 figures

Over the last decade, warm dark matter (WDM) has been repeatedly proposed as an alternative scenario to the standard cold dark matter (CDM) one, potentially resolving several disagreements between the CDM model and observations on small scales. Here, we reconsider the most important CDM small-scale discrepancies in the light of recent observational constraints on WDM. As a result, we find that a conventional thermal (or thermal-like) WDM cosmology with a particle mass in agreement with Lyman-$\alpha$ is nearly indistinguishable from CDM on the relevant scales and therefore fails to alleviate any of the small-scale problems. The reason for this failure is that the power spectrum of conventional WDM falls off too rapidly. To maintain WDM as a significantly different alternative to CDM, more evolved production mechanisms leading to multiple dark matter components or a gradually decreasing small-scale power spectrum have to be considered., Comment: Typo corrected in Eq 1. Does not affect any results

We study the effect of warm dark matter (WDM) on hydrodynamic simulations of galaxy formation as part of the Making Galaxies in a Cosmological Context (MaGICC) project. We simulate three different galaxies using three WDM candidates of 1, 2 and 5 keV and compare results with pure cold dark matter simulations. WDM slightly reduces star formation and produces less centrally concentrated stellar profiles. These effects are most evident for the 1 keV candidate but almost disappear for $m_{\mathrm{WDM}}>2$ keV. All simulations form similar stellar discs independent of WDM particle mass. In particular, the disc scale length does not change when WDM is considered. The reduced amount of star formation in the case of 1 keV particles is due to the effects of WDM on merging satellites which are on average less concentrated and less gas rich. The altered satellites cause a reduced starburst during mergers because they trigger weaker disc instabilities in the main galaxy. Nevertheless we show that disc galaxy evolution is much more sensitive to stellar feedback than it is to WDM candidate mass. Overall we find that WDM, especially when restricted to current observational constraints ($m_{\mathrm{WDM}}>2$ keV), has a minor impact on disc galaxy formation., Comment: 13 pages, 9 figures, 2 tables; minor clarifications added in results section, conclusions unchanged; accepted for publication in MNRAS

Using cosmological galaxy formation simulations from the MaGICC project, spanning more than three magnitudes in stellar mass (~10^7-3x10^{10} Msun), we trace the baryonic cycle of infalling gas from the virial radius through to its participation in the star formation process. An emphasis is placed upon the temporal history of chemical enrichment during its passage through the corona and CGM. We derive the distributions of time between gas crossing the virial radius and being accreted to the star forming region (which allows mixing within the corona), as well as the time between gas being accreted to the star forming region and then forming stars (which allows mixing within the disc). Significant numbers of stars are formed from gas that cycles back through the hot halo after first accreting to the star forming region. Gas entering high mass galaxies is pre-enriched in low mass proto-galaxies prior to entering the virial radius of the central progenitor, with only small amounts of primordial gas accreted, even at high redshift (z~5). After entering the virial radius, significant further enrichment occurs prior to the accretion of the gas to the star forming region, with gas that is feeding the star forming region surpassing 0.1Z by z=0. Mixing with halo gas, itself enriched via galactic fountains, is thus crucial in determining the metallicity at which gas is accreted to the disc. The lowest mass simulation (Mvir~2x10^{10}Msun, with M*~10^7Msun), by contrast, accretes primordial gas through the virial radius and onto the disc at all times. Much like the classical analytical solutions to the `G-dwarf problem', overproduction of low-metallicity stars is ameliorated by the inefficiency of star formation. Finally, gas outflow/metal removal rates from star forming regions as a function of galactic mass are presented., Comment: published in MNRAS

We use 31 simulated galaxies from the MaGICC project to investigate the effects of baryonic feedback on the density profiles of dark matter (DM) haloes. The sample covers a wide mass range: 9.4e9

Quasars represent the brightest Active Galactic Nuclei (AGN) in the Universe and are thought to indicate the location of prodigiously growing Black Holes (BHs), with luminosities as high as 10^48 erg/sec. It is often expected though that such an extremely energetic process will take place in the most massive bound structures in the dark matter (DM) distribution. We show that in contrast to this expectation, in a galaxy formation model which includes AGN feedback, quasars are predicted to live in average DM halo environments with typical masses of a few times 10^12 Msun. This fundamental prediction arises from the fact that quasar activity (i.e., BH accretion with luminosity greater than 10^46 erg/sec) is inhibited in DM haloes where AGN feedback operates. The galaxy hosts of quasars in our simulations are identified with over massive (in gas and stars) spheroidal galaxies, in which BH accretion is triggered via a galaxy merger or secular processes. We further show that the z=0 descendants of high redshift (z~6) QSOs span a wide range of morphologies, galaxy and halo masses. The z~6 BHs typically grow only by a modest factor by the present day. Remarkably, high redshift QSOs never inhabit the largest DM haloes at that time and their descendants are very seldom found in the most massive haloes at z=0. We also show that observationally it is very likely to find an enhancement in the abundance of galaxies around quasars at z~5. However, these enhancements are considerably weaker compared to the overdensities expected at the extreme peaks of the DM distribution. Thus, it is very unlikely that a quasar detected in the $z\gtrsim5$ Universe pinpoints the location of the progenitors of superclusters in the local Universe., Comment: 11 pages, 6 figures, submitted to MNRAS (comments are welcome)

We present a cosmological hydrodynamical simulation of a representative volume of the Universe, as part of the Making Galaxies in a Cosmological Context (MaGICC) project. MaGICC uses a thermal implementation for supernova and early stellar feedback. This work tests the feedback model at lower resolution across a range of galaxy masses, morphologies and merger histories. The simulated sample compares well with observations of high redshift galaxies ($z \ge 2$) including the stellar mass - halo mass ($M_\star - M_h$ ) relation, the Galaxy Stellar Mass Function (GSMF) at low masses ($M_\star \lt 5 \times 10^{10} M_\odot$ ) and the number density evolution of low mass galaxies. The poor match of $M_\star - M_h$ and the GSMF at high masses ($M_\star \ge 5 \times 10^{10} M_\odot$ ) indicates supernova feedback is insufficient to limit star formation in these haloes. At $z = 0$, our model produces too many stars in massive galaxies and slightly underpredicts the stellar mass around $L_\star$ mass galaxy. Altogether our results suggest that early stellar feedback, in conjunction with supernovae feedback, plays a major role in regulating the properties of low mass galaxies at high redshift., Comment: Simulation extended to z=0, conclusions unchanged. Accepted for publication in MNRAS

We analyse the structure and chemical enrichment of a Milky Way-like galaxy with a stellar mass of 2 10^{10} M_sun, formed in a cosmological hydrodynamical simulation. It is disk-dominated with a flat rotation curve, and has a disk scale length similar to the Milky Way's, but a velocity dispersion that is ~50% higher. Examining stars in narrow [Fe/H] and [\alpha/Fe] abundance ranges, we find remarkable qualitative agreement between this simulation and observations: a) The old stars lie in a thickened distribution with a short scale length, while the young stars form a thinner disk, with scale lengths decreasing, as [Fe/H] increases. b) Consequently, there is a distinct outward metallicity gradient. c) Mono-abundance populations exist with a continuous distribution of scale heights (from thin to thick). However, the simulated galaxy has a distinct and substantive very thick disk (h_z~1.5 kpc), not seen in the Milky Way. The broad agreement between simulations and observations allows us to test the validity of observational proxies used in the literature: we find in the simulation that mono-abundance populations are good proxies for single age populations (<1 Gyr) for most abundances., Comment: 11 page, 10 figures, submitted. Comments welcome