Far-Ultraviolet (FUV) radiation greatly exceeds ultraviolet, supernovae and winds in the energy budget of young star clusters but is poorly modelled in galaxy simulations. We present results of the first full isolated galaxy disk simulations to include FUV radiation self-consistently. 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 is crucial for producing warm neutral gas. 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, reproducing the tight Kennicutt-Schmidt relation observed for an individual galaxy, without producing too much hot ionized medium and with less disruption to the ISM compared to supernovae alone.

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.