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.

Young massive star clusters spanning $\sim 10^4 - 10^8 M_\odot$ in mass have been observed to have similar surface brightness profiles. Recent hydrodynamical simulations of star cluster formation have also produced star clusters with this structure. We argue analytically that this type of mass distribution arises naturally in the relaxation from a hierarchically-clustered distribution of stars into a monolithic star cluster through hierarchical merging. We show that arbitrary initial profiles will tend to converge to a universal profile under hierarchical merging, owing to phase-space mixing obeying certain conservation constraints. We perform $N$-body simulations of a pairwise merger of model star clusters and find that mergers readily produce the shallow surface brightness profiles observed in young massive clusters. Finally, we simulate the relaxation of a hierarchically-clustered mass distribution constructed from an idealized fragmentation model. Assuming only power-law spatial and kinematic scaling relations, these numerical experiments are able to reproduce the surface density profiles of observed young massive star clusters. Thus we provide physical motivation for the structure of young massive clusters within the paradigm of hierarchical star formation. This has important implications for the structure of nascent globular clusters.

Intrinsic alignments (IA), correlations between the intrinsic shapes and orientations of galaxies on the sky, are both a significant systematic in weak lensing and a probe of the effect of large-scale structure on galactic structure and angular momentum. In the era of precision cosmology, it is thus especially important to model IA with high accuracy. Efforts to use cosmological perturbation theory to model the dependence of IA on the large-scale structure have thus far been relatively successful; however, extant models do not consistently account for time evolution. In particular, advection of galaxies due to peculiar velocities alters the impact of IA, because galaxy positions when observed are generally different from their positions at the epoch when IA is believed to be set. In this work, we evolve the galaxy IA from the time of galaxy formation to the time at which they are observed, including the effects of this advection, and show how this process naturally leads to a dependence of IA on the velocity shear. We calculate the galaxy-galaxy-IA bispectrum to tree level (in the linear matter density) in terms of the evolved IA coefficients. We then discuss the implications for weak lensing systematics as well as for studies of galaxy formation and evolution. We find that considering advection introduces nonlocality into the bispectrum, and that the degree of nonlocality represents the memory of a galaxy's path from the time of its formation to the time of observation. We discuss how this result can be used to constrain the redshift at which IA is determined and provide Fisher estimation for the relevant measurements using the example of SDSS-BOSS.

In a recent work based on 3200 stacked H$\alpha$ maps of galaxies at $z \sim 1$, Nelson et al.~find evidence for `coherent star formation': the stacked SFR profiles of galaxies above (below) the 'star formation main sequence' (MS) are above (below) that of galaxies on the MS at all radii. One might interpret this result as inconsistent with highly bursty star formation and evidence that galaxies evolve smoothly along the MS rather than crossing it many times. We analyze six simulated galaxies at $z\sim1$ from the Feedback in Realistic Environments (FIRE) project in a manner analogous to the observations to test whether the above interpretations are correct. The trends in stacked SFR profiles are qualitatively consistent with those observed. However, SFR profiles of individual galaxies are much more complex than the stacked profiles: the former can be flat or even peak at large radii because of the highly clustered nature of star formation in the simulations. Moreover, the SFR profiles of individual galaxies above (below) the MS are not systematically above (below) those of MS galaxies at all radii. We conclude that the time-averaged coherent star formation evident stacks of observed galaxies is consistent with highly bursty, clumpy star formation of individual galaxies and is not evidence that galaxies evolve smoothly along the MS.