Modelling the cosmological Lyman-Werner background radiation field in the early Universe

The Lyman-Werner (LW) radiation field is a key ingredient in the chemothermal evolution of gas in the early Universe, as it dissociates H-2 molecules, the primary cooling channel in an environment devoid of metals and dust. Despite its important role, it is still not implemented in cosmological simulations on a regular basis, in contrast to the ionizing UV background. This is in part due to uncertainty in the source modelling, their spectra and abundance, as well as the detailed physics involved in the propagation of the photons and their interactions with the molecules. The goal of this work is to produce an accurate model of the LW radiation field at z = 6, by post-processing the physics-rich high-resolution FiBY simulation. Our novelties include updated cross-sections for H-2, H- and H-2(+) chemical species, IGM absorption by neutral Hydrogen and various spectral models for Population III and Population II stars. With our fiducial set of parameters, we show that the mean LW intensity steadily increases by three orders of magnitude from z similar to 23 to z similar to 6, while spatial inhomogeneities originate from massive star-forming galaxies that dominate the photon budget up to a distance of similar to 100 proper kpc. Our model can be easily applied to other simulations or semi-analytical models as an external radiation field that regulates the formation of stars and massive black hole seeds in high-z low-mass haloes.