We present a series of high-resolution cosmological simulations 1 of galaxy formation to z = 0, spanning halo masses ∼ 10 8 − 10 13 M , and stellar masses ∼ 10 4 − 10 11 M . Our simulations include fully explicit treatment of the multi-phase ISM & stellar feedback. The stellar feedback inputs (energy, momentum, mass, and metal fluxes) are taken directly from stellar population models. These sources of feedback, with zero adjusted parameters, reproduce the observed relation between stellar and halo mass up to M halo ∼ 10 12 M . We predict weak redshift evolution in the M * − M halo relation, consistent with current constraints to z > 6. We find that the M * − M halo relation is insensitive to numerical details, but is sensitive to feedback physics. Simulations with only supernova feedback fail to reproduce observed stellar masses, particularly in dwarf and high-redshift galaxies: radiative feedback (photo-heating and radiation pressure) is necessary to destroy GMCs and enable efficient coupling of later supernovae to the gas. Star formation rates agree well with the observed Kennicutt relation at all redshifts. The galaxy-averaged Kennicutt relation is very different from the numerically imposed law for converting gas into stars, and is determined by self-regulation via stellar feedback. Feedback reduces star formation rates and produces reservoirs of gas that lead to rising late-time star formation histories, significantly different from halo accretion histories. Feedback also produces large short-timescale variability in galactic SFRs, especially in dwarfs. These properties are not captured by common "sub-grid" wind models.

The Feedback In Realistic Environments (FIRE) project explores feedback in cosmological galaxy formation simulations. Previous FIRE simulations used an identical source code ("FIRE-1") for consistency. Motivated by the development of more accurate numerics -including hydrodynamic solvers, gravitational softening, and supernova coupling algorithms -and exploration of new physics (e.g. magnetic fields), we introduce "FIRE-2", an updated numerical implementation of FIRE physics for the GIZMO code. We run a suite of simulations and compare against FIRE-1: overall, FIRE-2 improvements do not qualitatively change galaxy-scale properties. We pursue an extensive study of numerics versus physics. Details of the star-formation algorithm, cooling physics, and chemistry have weak effects, provided that we include metal-line cooling and star formation occurs at higher-than-mean densities. We present new resolution criteria for high-resolution galaxy simulations. Most galaxy-scale properties are robust to numerics we test, provided: (1) Toomre masses are resolved; (2) feedback coupling ensures conservation, and (3) individual supernovae are time-resolved. Stellar masses and profiles are most robust to resolution, followed by metal abundances and morphologies, followed by properties of winds and circum-galactic media (CGM). Central (∼kpc) mass concentrations in massive (> L * ) galaxies are sensitive to numerics (via trapping/recycling of winds in hot halos). Multiple feedback mechanisms play key roles: supernovae regulate stellar masses/winds; stellar mass-loss fuels late star formation; radiative feedback suppresses accretion onto dwarfs and instantaneous star formation in disks. We provide all initial conditions and numerical algorithms used.

We investigate the birth and evolution of Galactic isolated radio pulsars. We begin by estimating their birth space velocity distribution from proper-motion measurements of Brisken and coworkers. We find no evidence for multimodality of the distribution and favor one in which the absolute one-dimensional velocity components are exponentially distributed and with a three-dimensional mean velocity of 380 þ40 À60 km s À1. We then proceed with a Monte Carlobased population synthesis, modeling the birth properties of the pulsars, their time evolution, and their detection in the Parkes and Swinburne Multibeam surveys. We present a population model that appears generally consistent with the observations. Our results suggest that pulsars are born in the spiral arms, with a galactocentric radial distribution that is well described by the functional form proposed by Yusifov & Küçük, in which the pulsar surface density peaks at radius $3 kpc. The birth spin period distribution extends to several hundred milliseconds, with no evidence of multimodality. Models that assume the radio luminosities of pulsars to be independent of the spin periods and period derivatives are inadequate, as they lead to the detection of too many old simulated pulsars in our simulations. Dithered radio luminosities proportional to the square root of the spin-down luminosity accommodate the observations well and provide a natural mechanism for the pulsars to dim uniformly as they approach the death line, avoiding an observed pileup on the latter. There is no evidence for significant torque decay (due to magnetic field decay or otherwise) over the lifetime of the pulsars as radio sources ($100 Myr). Finally, we estimate the pulsar birthrate and total number of pulsars in the Galaxy.

We present an analysis of the galaxy-scale gaseous outflows from the FIRE (Feedback in Realistic Environments) simulations. This suite of hydrodynamic cosmological zoom simulations resolves formation of star-forming giant molecular clouds to z = 0, and features an explicit stellar feedback model on small scales. Our simulations reveal that high redshift galaxies undergo bursts of star formation followed by powerful gusts of galactic outflows that eject much of the ISM and temporarily suppress star formation. At low redshift, however, sufficiently massive galaxies corresponding to L*-progenitors develop stable disks and switch into a continuous and quiescent mode of star formation that does not drive outflows far into the halo. Mass-loading factors for winds in L*-progenitors are η ≈ 10 at high redshift, but decrease to η 1 at low redshift. Although lower values of η are expected as halos grow in mass over time, we show that the strong suppression of outflows with decreasing redshift cannot be explained by mass evolution alone. Circumgalactic outflow velocities are variable and broadly distributed, but typically range between one and three times the circular velocity of the halo. Much of the ejected material builds a reservoir of enriched gas within the circumgalactic medium, some of which could be later recycled to fuel further star formation. However, a fraction of the gas that leaves the virial radius through galactic winds is never regained, causing most halos with mass M h ≤ 10 12 M to be deficient in baryons compared to the cosmic mean by z = 0.