We present an update to the multiphase SPH galaxy formation code by Scannapieco et al. We include a more elaborate treatment of the production of metals, cooling rates based on individual element abundances, and a scheme for the turbulent diffusion of metals. Our SN feedback model now transfers energy to the ISM in kinetic and thermal form, and we include a prescription for the effects of radiation pressure from massive young stars on the ISM. We calibrate our new code on the well studied Aquarius haloes and then use it to simulate a sample of 16 galaxies with halo masses between 1 × 10 11 and 3 × 10 12 M ⊙ . In general, the stellar masses of the sample agree well with the stellar mass to halo mass relation inferred from abundance matching techniques for redshifts z = 0 − 4. There is however a tendency to overproduce stars at z > 4 and to underproduce them at z < 0.5 in the least massive haloes. Overly high SFRs at z < 1 for the most massive haloes are likely connected to the lack of AGN feedback in our model. The simulated sample also shows reasonable agreement with observed star formation rates, sizes, gas fractions and gas-phase metallicities at z = 0 − 3. Remaining discrepancies can be connected to deviations from predictions for star formation histories from abundance matching. At z = 0, the model galaxies show realistic morphologies, stellar surface density profiles, circular velocity curves and stellar metallicities, but overly flat metallicity gradients. 15 out of 16 of our galaxies contain disk components with kinematic disk fraction ranging between 15 and 65 %. The disk fraction depends on the time of the last destructive merger or misaligned infall event. Considering the remaining shortcomings of our simulations we conclude that even higher kinematic disk fractions may be possible for ΛCDM haloes with quiet merger histories, such as the Aquarius haloes.

We use proper motions and parallaxes from the new reduction of Hipparcos data and Geneva-Copenhagen radial velocities for a complete sample of ~15000 main-sequence and subgiant stars, and new Padova isochrones to constrain the kinematics and star-formation history of the solar neighbourhood. We rederive the solar motion and the structure of the local velocity ellipsoids. When the principal velocity dispersions are assumed to increase with time as t^\beta, the index \beta is larger for \sigma_W (\beta_W~0.45) than for \sigma_U (\beta_U~0.31). For the three-dimensional velocity dispersion we obtain \beta=0.35. We exclude saturation of disc heating after ~3 Gyr as proposed by Quillen & Garnett(2000). Saturation after >~4 Gyr combined with an abrupt increase in velocity dispersion for the oldest stars cannot be excluded. For all our models the star-formation rate is declining, being a factor 2-7 lower now than it was at the beginning. Models in which the SFR declines exponentially favour very high disc ages between 11.5 and 13 Gyr and exclude ages below ~10.5 Gyr as they yield worse fits to the number density and velocity dispersion of red stars. Models in which the SFR is the sum of two declining exponentials representing the thin and thick discs favour ages between 10.5 and 12 Gyr with a lower limit of ~10.0 Gyr. Although in our models the star-formation rate peaked surprisingly early, the mean formation time of solar-neighbourhood stars is later than in ab-initio models of galaxy formation, probably on account of weaknesses in such models.Comment: 18 pages, 18 figures, accepted for publication in MNRA

We analyze the radial distribution of HI gas for 23 disk galaxies with unusually high HI content from the Bluedisk sample, along with a similar-sized sample of "normal" galaxies. We propose an empirical model to fit the radial profile of the HI surface density, an exponential function with a depression near the center. The radial HI surface density profiles are very homogeneous in the outer regions of the galaxy; the exponentially declining part of the profile has a scale-length of ∼ 0.18 R1, where R1 is the radius where the column density of the HI is 1 M ⊙ pc −2 . This holds for all galaxies, independent of their stellar or HI mass. The homogenous outer profiles, combined with the limited range in HI surface density in the nonexponential inner disk, results in the well-known tight relation between HI size and HI mass. By comparing the radial profiles of the HI-rich galaxies with those of the control systems, we deduce that in about half the galaxies, most of the excess gas lies outside the stellar disk, in the exponentially declining outer regions of the HI disk. In the other half, the excess is more centrally peaked. We compare our results with existing smoothed-particle hydrodynamical simulations and semi-analytic models of disk galaxy formation in a Λ Cold Dark Matter universe. Both the hydro simulations and the semi-analytic models reproduce the HI surface density profiles and the HI size-mass relation without further tuning of the simulation and model inputs. In the semi-analytic models, the universal shape of the outer HI radial profiles is a consequence of the assumption that infalling gas is always distributed exponentially. The conversion of atomic gas to molecular form explains the limited range of HI surface densities in the inner disk. These two factors produce the tight HI mass-size relation.