We compare the results of various cosmological gas-dynamical codes used to simulate the formation of a galaxy in the Λ cold dark matter structure formation paradigm. The various runs (13 in total) differ in their numerical hydrodynamical treatment [smoothed particle hydrodynamics (SPH), moving mesh and adaptive mesh refinement] but share the same initial conditions and adopt in each case their latest published model of gas cooling, star formation and feedback. Despite the common halo assembly history, we find large code-to-code variations in the stellar mass, size, morphology and gas content of the galaxy at z= 0, due mainly to the different implementations of star formation and feedback. Compared with observation, most codes tend to produce an overly massive galaxy, smaller and less gas rich than typical spirals, with a massive bulge and a declining rotation curve. A stellar disc is discernible in most simulations, although its prominence varies widely from code to code. There is a well-defined trend between the effects of feedback and the severity of the disagreement with observed spirals. In general, models that are more effective at limiting the baryonic mass of the galaxy come closer to matching observed galaxy scaling laws, but often to the detriment of the disc component. Although numerical convergence is not particularly good for any of the codes, our conclusions hold at two different numerical resolutions. Some differences can also be traced to the different numerical techniques; for example, more gas seems able to cool and become available for star formation in grid-based codes than in SPH. However, this effect is small compared to the variations induced by different feedback prescriptions. We conclude that state-of-the-art simulations cannot yet uniquely predict the properties of the baryonic component of a galaxy, even when the assembly history of its host halo is fully specified. Developing feedback algorithms that can effectively regulate the mass of a galaxy without hindering the formation of high angular momentum stellar discs remains a challenge

We present a first study of the effect of local photoionising radiation on gas cooling in smoothed particle hydrodynamics (SPH) simulations of galaxy formation. We explore the combined effect of ionising radiation from young and old stellar populations. The method computes the effect of multiple radiative sources using the same tree algorithm used for gravity, so it is computationally efficient and well resolved. The method foregoes calculating absorption and scattering in favour of a constant escape fraction for young stars to keep the calculation efficient enough to simulate the entire evolution of a galaxy in a cosmological context to the present day. This allows us to quantify the effect of the local photoionisation feedback through the whole history of a galaxy's formation. The simulation of a Milky Way-like galaxy using the local photoionisation model forms ∼ 40 % less stars than a simulation that only includes a standard uniform background UV field. The local photoionisation model decreases star formation by increasing the cooling time of the gas in the halo and increasing the equilibrium temperature of dense gas in the disc. Coupling the local radiation field to gas cooling from the halo provides a preventive feedback mechanism which keeps the central disc light and produces slowly rising rotation curves without resorting to extreme feedback mechanisms. These preliminary results indicate that the effect of local photoionising sources is significant and should not be ignored in models of galaxy formation.

We examine gas accretion and subsequent star formation in representative galaxies from the McMaster Unbiased Galaxy Simulations ). Accreted gas is bimodal with a natural temperature division at 10 5 K, near the peak of the cooling curve. Cold-mode accretion dominates inflows at early times, creating a peak in total accretion at redshift z=2-4 and declining exponentially below z∼2. Hot-mode accretion peaks near z=1-2 and declines gradually. Hot-mode exceeds cold-mode accretion at z∼1.8 for all four galaxies rather than when the galaxy reaches a characteristic mass. Cold-mode accretion can fuel immediate star formation, while hot-mode accretion preferentially builds a large, hot gas reservoir in the halo. Late-time star formation relies on reservoir gas accreted 2-8 Gyr prior. Thus, the reservoir allows the star formation rate to surpass the current overall gas accretion rate. Stellar feedback cycles gas from the interstellar medium back into the hot reservoir. Stronger feedback results in more gas cycling, gas removal in a galactic outflow and less star formation overall, enabling simulations to match the observed star formation history. For lower mass galaxies in particular, strong feedback can delay the star formation peak to z=1-2 from the accretion peak at z=2-4.

We present TREVR (Tree-based REVerse Ray Tracing), a general algorithm for computing the radiation field, including absorption, in astrophysical simulations. TREVR is designed to handle large numbers of sources and absorbers; it is based on a tree data structure and is thus suited to codes that use trees for their gravity or hydrodynamics solvers (e.g. Adaptive Mesh Refinement). It achieves computational speed while maintaining a specified accuracy via controlled lowering of the resolution of both sources and rays from each source. TREVR computes the radiation field in order N log N source time without absorption and order N log N source log N time with absorption. These scalings arise from merging sources of radiation according to an opening angle criterion and walking the tree structure to trace a ray to a depth that gives the chosen accuracy for absorption. The absorption-depth refinement criterion is unique to TREVR. We provide a suite of tests demonstrating the algorithm's ability to accurately compute fluxes, ionization fronts and shadows.

Far-Ultraviolet (FUV) radiation greatly exceeds ultraviolet, supernovae and winds in the energy budget of young star clusters but is poorly modeled in galaxy simulations. We present results of the first isolated galaxy disk simulations to include photo-electric heating of gas via dust grains from FUV radiation self-consistently, using a ray-tracing approach that calculates optical depths along the source-receiver sight-line. This is the first science application of the TREVR radiative transfer algorithm. We find that FUV radiation alone cannot regulate star formation. However, FUV radiation produces warm neutral gas and is able to produce regulated galaxies with realistic scale heights. FUV is also a long-range feedback and is more important in the outer disks of galaxies. We also use the super-bubble feedback model, which depends only on the supernova energy per stellar mass, is more physically realistic than common, parameter-driven alternatives and thus better constrains supernova feedback impacts. FUV and supernovae together can regulate star formation without producing too much hot ionized medium and with less disruption to the ISM compared to supernovae alone.