For almost two decades the properties of 'dwarf' galaxies have challenged the cold dark matter (CDM) model of galaxy formation. Most observed dwarf galaxies consist of a rotating stellar disk embedded in a massive dark-matter halo with a near-constant-density core. Models based on the dominance of CDM, however, invariably form galaxies with dense spheroidal stellar bulges and steep central dark-matter profiles, because low-angular-momentum baryons and dark matter sink to the centres of galaxies through accretion and repeated mergers. Processes that decrease the central density of CDM halos 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. Here we report hydrodynamical simulations (in a framework 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.

We have carried out a hydrodynamical code comparison study of interacting multiphase fluids. The two commonly used techniques of grid and smoothed particle hydrodynamics (SPH) show striking differences in their ability to model processes that are fundamentally important across many areas of astrophysics. Whilst Eulerian grid based methods are able to resolve and treat important dynamical instabilities, such as Kelvin-Helmholtz or Rayleigh-Taylor, these processes are poorly or not at all resolved by existing SPH techniques. We show that the reason for this is that SPH, at least in its standard implementation, introduces spurious pressure forces on particles in regions where there are steep density gradients. This results in a boundary gap of the size of the SPH smoothing kernel over which information is not transferred.Comment: 15 pages, 13 figures, to be submitted to MNRAS. For high-resolution figures, please see http://www-theorie.physik.unizh.ch/~agertz

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

Simulations of the formation of late-type spiral galaxies in a cold dark matter (ΛCDM) universe have traditionally failed to yield realistic candidates. Here we report a new cosmological N-body/smooth particle hydrodynamic simulation of extreme dynamic range in which a close analog of a Milky Way disk galaxy arises naturally. Named "Eris," the simulation follows the assembly of a galaxy halo of mass M vir = 7.9 × 10 11 M with a total of N = 18.6 million particles (gas + dark matter + stars) within the final virial radius, and a force resolution of 120 pc. It includes radiative cooling, heating from a cosmic UV field and supernova explosions (blastwave feedback), a star formation recipe based on a high gas density threshold (n SF = 5 atoms cm −3 rather than the canonical n SF = 0.1 atoms cm −3 ), and neglects any feedback from an active galactic nucleus. Artificial images are generated to correctly compare simulations with observations. At the present epoch, the simulated galaxy has an extended rotationally supported disk with a radial scale length R d = 2.5 kpc, a gently falling rotation curve with circular velocity at 2.2 disk scale lengths of V 2.2 = 214 km s −1 , an i-band bulge-to-disk ratio B/D = 0.35, and a baryonic mass fraction within the virial radius that is 30% below the cosmic value. The disk is thin, has a typical H I-to-stellar mass ratio, is forming stars in the region of the Σ SFR -Σ H i plane occupied by spiral galaxies, and falls on the photometric Tully-Fisher and the stellar-mass-halo-virial-mass relations. Hot (T > 3 × 10 5 K) X-ray luminous halo gas makes up only 26% of the universal baryon fraction and follows a "flattened" density profile ∝ r −1.13 out to r = 100 kpc. Eris appears then to be the first cosmological hydrodynamic simulation in which the galaxy structural properties, the mass budget in the various components, and the scaling relations between mass and luminosity are all consistent with a host of observational constraints. A twin simulation with a low star formation density threshold results in a galaxy with a more massive bulge and a much steeper rotation curve, as in previously published work. A high star formation threshold appears therefore key in obtaining realistic late-type galaxies, as it enables the development of an inhomogeneous interstellar medium where star formation and heating by supernovae occur in a clustered fashion. The resulting outflows at high redshifts reduce the baryonic content of galaxies and preferentially remove low-angular-momentum gas, decreasing the mass of the bulge component. Simulations of even higher resolution that follow the assembly of galaxies with different merger histories shall be used to verify our results.