Forming Disk Galaxies in Lambda CDM Simulations

We used fully cosmological, high resolution N-body + SPH simulations to follow the formation of disk galaxies with rotational velocities between 135 and 270 km/sec in a Lambda CDM universe. The simulations include gas cooling, star formation, the effects of a uniform UV background and a physically motivated description of feedback from supernovae. The host dark matter halos have a spin and last major merger redshift typical of galaxy sized halos as measured in recent large scale N--Body simulations. The simulated galaxies form rotationally supported disks with realistic exponential scale lengths and fall on both the I-band and baryonic Tully Fisher relations. An extended stellar disk forms inside the Milky Way sized halo immediately after the last major merger. The combination of UV background and SN feedback drastically reduces the number of visible satellites orbiting inside a Milky Way sized halo, bringing it in fair agreement with observations. Our simulations predict that the average age of a primary galaxy's stellar population decreases with mass, because feedback delays star formation in less massive galaxies. Galaxies have stellar masses and current star formation rates as a function of total mass that are in good agreement with observational data. We discuss how both high mass and force resolution and a realistic description of star formation and feedback are important ingredients to match the observed properties of galaxies.Comment: Revised version after the referee's comments. Conclusions unchanged. 2 new plots. MNRAS in press. 20 plots. 21 page

Forming disc galaxies in ΛCDM simulations

We used fully cosmological, high-resolution N-body + smooth particle hydrodynamic (SPH) simulations to follow the formation of disc galaxies with rotational velocities between 135 and 270 km s−1 in a Λ cold dark matter (CDM) universe. The simulations include gas cooling, star formation, the effects of a uniform ultraviolet (UV) background and a physically motivated description of feedback from supernovae (SNe). The host dark matter haloes have a spin and last major merger redshift typical of galaxy-sized haloes as measured in recent large-scale N-body simulations. The simulated galaxies form rotationally supported discs with realistic exponential scalelengths and fall on both the I band and baryonic Tully-Fisher relations. An extended stellar disc forms inside the Milky Way (MW)-sized halo immediately after the last major merger. The combination of UV background and SN feedback drastically reduces the number of visible satellites orbiting inside a MW-sized halo, bringing it in fair agreement with observations. Our simulations predict that the average age of a primary galaxy's stellar population decreases with mass, because feedback delays star formation in less massive galaxies. Galaxies have stellar masses and current star formation rates as a function of total mass that are in good agreement with observational data. We discuss how both high mass and force resolution and a realistic description of star formation and feedback are important ingredients to match the observed properties of galaxie

Bulgeless dwarf galaxies and dark matter cores from supernova-driven outflows

For almost two decades the properties of ‘dwarf’ galaxies have challenged the cold dark matter (CDM) model of galaxy formation^1. Most observed dwarf galaxies consist of a rotating stellar disk^2 embedded in a massive dark-matter halo with a near-constant-density core^3. Models based on the dominance of CDM, however, invariably form galaxies with dense spheroidal stellar bulges and steep central dark-matter profiles^(4,5,6,) because low-angular-momentum baryons and dark matter sink to the centres of galaxies through accretion and repeated mergers^7. Processes that decrease the central density of CDM halos^8 have been identified, but have not yet reconciled theory with observations of present-day dwarfs. This failure is potentially catastrophic for the CDM model, possibly requiring a different dark-matter particle candidate^9. Here we report hydrodynamical simulations (in a framework^(10) assuming the presence of CDM and a cosmological constant) in which the inhomogeneous interstellar medium is resolved. Strong outflows from supernovae remove low-angular-momentum gas, which inhibits the formation of bulges and decreases the dark-matter density to less than half of what it would otherwise be within the central kiloparsec. The analogues of dwarf galaxies—bulgeless and with shallow central dark-matter profiles—arise naturally in these simulations

Forming a large disc galaxy from a z < 1 major merger

Using high-resolution SPH simulations in a fully cosmological Λ cold dark matter context, we study the formation of a bright disc-dominated galaxy that originates from a ‘wet' major merger at z= 0.8. The progenitors of the disc galaxy are themselves disc galaxies that formed from early major mergers between galaxies with blue colours. A substantial thin stellar disc grows rapidly following the last major merger and the present-day properties of the final remnant are typical of early-type spiral galaxies, with an i-band bulge-to-disc ratio ∼0.65, a disc scalelength of 7.2 kpc, g−r= 0.5 mag, an H i linewidth (W20/2) of 238 km s−1 and total magnitude i=−22.4. The key ingredients for the formation of a dominant stellar disc component after a major merger are (i) substantial and rapid accretion of gas through cold flows followed at late times by cooling of gas from the hot phase, (ii) supernova feedback that is able to partially suppress star formation during mergers and (iii) relative fading of the spheroidal component. The gas fraction of the progenitors' discs does not exceed 25 per cent at z < 3, emphasizing that the continuous supply of gas from the local environment plays a major role in the regrowth of discs and in keeping the galaxies blue. The results of this simulation alleviate the problem posed for the existence of disc galaxies by the high likelihood of interactions and mergers for galaxy-sized haloes at relatively low

The Formation of a Realistic Disk Galaxy in Lambda Dominated Cosmologies

We simulate the formation of a realistic disk galaxy within the hierarchical scenario of structure formation and study its internal properties to the present epoch. We compare results from a LambdaCDM simulation with a LambdaWDM (2keV) simulation that forms significantly less small scale structure. We show how high mass and force resolution in both the gas and dark matter components play an important role in solving the angular momentum catastrophe claimed from previous simulations of galaxy formation within the hierarchical framework. The stellar material in the disk component has a final specific angular momentum equal to 40% and 90% of that of the dark halo in the LambdaCDM and LambdaWDM models respectively. The LambdaWDM galaxy has a drastically reduced satellite population and a negligible stellar spheroidal component. Encounters with satellites play only a minor role in disturbing the disk. Satellites possess a variety of star formation histories linked to mergers and pericentric passages along their orbit around the primary galaxy. In both cosmologies, the galactic halo retains most of the baryons accreted and builds up a hot gas phase with a substantial X-ray emission. Therefore, while we have been successful in creating a realistic stellar disk in a massive galaxy within the LambdaCDM scenario, energy injection emerges as necessary ingredient to reduce the baryon fraction in galactic halos, independent of the cosmology adopted. (abridged)Comment: ApJ in press. Images and movies at http://hpcc.astro.washington.edu/faculty/fabio/galform.html Significantly expanded revised version. (9 pages vs the original 4

Dark Matter signals from Draco and Willman 1: Prospects for MAGIC II and CTA

The next generation of ground-based Imaging Air Cherenkov Telescopes (IACTs) will play an important role in indirect dark matter searches. In this article, we consider two particularly promising candidate sources for dark matter annihilation signals, the nearby dwarf galaxies Draco and Willman 1, and study the prospects of detecting such a signal for the soon-operating MAGIC II telescope system as well as for the planned installation of CTA, taking special care of describing the experimental features that affect the detectional prospects. For the first time in such a study, we fully take into account the effect of internal bremsstrahlung, which has recently been shown to considerably enhance, in some cases, the gamma-ray flux at the high energies where Atmospheric Cherenkov Telescopes operate, thus leading to significantly harder annihilation spectra than traditionally considered. While the detection of the spectral features introduced by internal bremsstrahlung would constitute a smoking gun signature for dark matter annihilation, we find that for most models the overall flux still remains at a level that will be challenging to detect unless one adopts rather (though by no means overly) optimistic astrophysical assumptions about the distribution of dark matter in the dwarfs.Comment: 10 pages, 4 figures, minor changes, matches the published version (JCAP

The Via Lactea INCITE Simulation: Galactic Dark Matter Substructure at High Resolution

It is a clear unique prediction of the cold dark matter paradigm of cosmological structure formation that galaxies form hierarchically and are embedded in massive, extended dark halos teeming with self-bound substructure or "subhalos". The amount and spatial distribution of subhalos around their host provide unique information and clues on the galaxy assembly process and the nature of the dark matter. Here we present results from the Via Lactea INCITE simulation, a one billion particle, one million cpu-hour simulation of the formation and evolution of a Galactic dark matter halo and its substructure population.Comment: 10 pages, Proceedings of the SciDAC 2008 conference, (Seattle, July 13-17, 2008

The Observed and Predicted Spatial Distribution of Milky Way Satellite Galaxies

We review evidence that the census of Milky Way satellites similar to those known may be incomplete at low latitude due to obscuration and in the outer halo due to a decreasing sensitivity to dwarf satellites with distance. We evaluate the possible impact that incompleteness has on comparisons with substructure models by estimating corrections to the known number of dwarfs using empirical and theoretical models. If we assume that the true distribution of Milky Way satellites is uniform with latitude, then we estimate a 33% incompleteness in the total number of dwarfs due to obscuration at low latitude. Similarly, if we suppose that the radial distribution of Milky Way satellites matches that of M31, or that of the oldest sub-halos or the most massive sub-halos in a simulation, we estimate a total number of Milky Way dwarfs ranging from 1 -- 3 times the known population. Although the true level of incompleteness is quite uncertain, the fact that our extrapolations yield average total numbers of MW dwarfs that are realistically 1.5 -- 2 times the known population, shows that incompleteness needs to be taken seriously when comparing to models of dwarf galaxy formation. Interestingly, the radial distribution of the oldest sub-halos in a Lambda+CDM simulation of a Milky Way-like galaxy possess a close match to the observed distribution of M31's satellites, which suggests that reionization may be an important factor controlling the observability of sub-halos. We also assess the prospects for a new SDSS search for Milky Way satellites to constrain the possible incompleteness in the outer halo.Comment: 9 pages, 7 figures. Replaced with MNRAS accepted versio

The Origin and Evolution of the Mass-Metallicity Relationship for Galaxies: Results from Cosmological N-Body Simulations

We examine the origin and evolution of the mass-metallicity relationship (MZR, M-Z) for galaxies using high resolution cosmological SPH + N-Body simulations that include a physically motivated description of supernovae feedback and subsequent metal enrichment. We discriminate between two sources that may contribute to the origin of the MZR: 1) metal and baryon loss due to gas outflow, or 2) inefficient star formation at the lowest galaxy masses. Our simulated galaxies reproduce the observed MZR in shape and normalization both at z=0 and z=2. We find that baryon loss occurs due to UV heating before star formation turns on in galaxies with M_baryon < 10^8 M_sun, but that some gas loss due to supernovae induced winds is required to subsequently reproduce the low effective chemical yield observed in low mass galaxies. Despite this, we show that low star formation efficiencies, regulated by supernovae feedback, are primarily responsible for the lower metallicities of low mass galaxies and the overall M-Z trend. We find that the shape of the MZR is relatively constant with redshift, but that its normalization increases with time. Simulations with no energy feedback from supernovae overproduce metals at low galaxy masses by rapidly transforming a large fraction of their gas into stars. Despite the fact that our low mass galaxies have lost a majority of their baryons, they are still the most gas rich objects in our simulations due to their low star formation efficiencies.Comment: 3 figures. Accepted for ApJL, in pres

Observations of the GRB afterglow ATLAS17aeu and its possible association with GW170104

We report the discovery and multi-wavelength data analysis of the peculiar optical transient, ATLAS17aeu. This transient was identified in the skymap of the LIGO gravitational wave event GW170104 by our ATLAS and Pan-STARRS coverage. ATLAS17aeu was discovered 23.1hrs after GW170104 and rapidly faded over the next 3 nights, with a spectrum revealing a blue featureless continuum. The transient was also detected as a fading x-ray source by Swift and in the radio at 6 and 15 GHz. A gamma ray burst GRB170105A was detected by 3 satellites 19.04hrs after GW170104 and 4.10hrs before our first optical detection. We analyse the multi-wavelength fluxes in the context of the known GRB population and discuss the observed sky rates of GRBs and their afterglows. We find it statistically likely that ATLAS17aeu is an afterglow associated with GRB170105A, with a chance coincidence ruled out at the 99\% confidence or 2.6$\sigma$. A long, soft GRB within a redshift range of $1 \lesssim z \lesssim 2.9$ would be consistent with all the observed multi-wavelength data. The Poisson probability of a chance occurrence of GW170104 and ATLAS17aeu is $p=0.04$. This is the probability of a chance coincidence in 2D sky location and in time. These observations indicate that ATLAS17aeu is plausibly a normal GRB afterglow at significantly higher redshift than the distance constraint for GW170104 and therefore a chance coincidence. However if a redshift of the faint host were to place it within the GW170104 distance range, then physical association with GW170104 should be considered.Comment: 16 pages, 6 figures, accepted to Ap