Aims: Mapping the interstellar medium in 3D provides a wealth of insights into its inner working. The Milky Way is the only galaxy, for which a detailed 3D mapping can be achieved in principle. In this paper, we reconstruct the dust density in and around the local super-bubble. Methods: The combined data from surveys such as Gaia, 2MASS, PANSTARRS, and ALLWISE provide the necessary information to make detailed maps of the interstellar medium in our surrounding. To this end, we use variational inference and Gaussian processes to model the dust extinction density, exploiting its intrinsic correlations. Results: We reconstruct a highly resolved dust map, showing the nearest dust clouds at a distance of up to 400\,pc with a resolution of 1\,pc. Conclusions: Our reconstruction provides insights into the structure of the interstellar medium. We compute summary statistics of the spectral index and the 1-point function of the logarithmic dust extinction density, which may constrain simulations of the interstellar medium that achieve similar resolution.

The recent discovery of high redshift dusty galaxies implies a rapid dust enrichment of their interstellar medium (ISM). To interpret these observations, we run a cosmological simulation in a 30$h^{-1}$ cMpc/size volume down to $z \approx 4$. We use the hydrodynamical code dustyGadget, which accounts for the production of dust by stellar populations and its evolution in the ISM. We find that the cosmic dust density parameter ($\Omega_{\rm d}$) is mainly driven by stellar dust at $z \gtrsim 10$, so that mass- and metallicity-dependent yields are required to assess the dust content in the first galaxies. At $z \lesssim 9$ the growth of grains in the ISM of evolved systems (Log$(M_{\star}/M_{\odot})>8.5$) significantly increases their dust mass, in agreement with observations in the redshift range $4 \lesssim z < 8$. Our simulation shows that the variety of high redshift galaxies observed with ALMA can naturally be accounted for by modeling the grain-growth timescale as a function of the physical conditions in the gas cold phase. In addition, the trends of dust-to-metal (DTM) and dust-to-gas (${\cal D}$) ratios are compatible with the available data. A qualitative investigation of the inhomogeneous dust distribution in a representative massive halo at $z \approx 4$ shows that dust is found from the central galaxy up to the closest satellites along polluted filaments with $\rm Log({\cal D}) \leq -2.4$, but sharply declines at distances $d \gtrsim 30$ kpc along many lines of sight, where $\rm Log({\cal D}) \lesssim -4.0$.

We present a new dust extension to the Monte Carlo radiative transfer code crash, which enables it to simulate the propagation of ionizing radiation through mixtures of gas and dust. The new code is applied to study the impact of dust absorption on idealized galactic H II regions and on small scale reionization. We find that H II regions are reduced in size by the presence of dust, while their inner temperature and ionization structure remain largely unaffected. In the small scale reionization simulation, dust hardens ionization fronts and delays the overlap of ionized bubbles. This effect is found to depend only weakly on the assumed abundance of dust in underdense regions.

We investigate the impact of quasars on the ionization of the surrounding intergalactic medium (IGM) with the radiative transfer code \texttt{CRASH4}, now accounting for X-rays and secondary electrons. After comparing with analytic solutions, we post-process a cosmic volume ($ \approx 1.5\times 10^4$ Mpc$^3 h^{-3}$) containing a ULAS J1120+0641-like quasar (QSO) hosted by a $5 \times 10^{11} {\rm M}_\odot h^{-1}$ dark matter (DM) halo. We find that: (i) the average HII region ($R\sim3.2$~pMpc in a lifetime $t_f = 10^7$~yrs) is mainly set by UV flux, in agreement with semi-analytic scaling relations; (ii) a largely neutral ($x_{\textrm{HII}} < 0.001$), warm ($T\sim 10^3$~K) tail extends up to few Mpc beyond the ionization front, as a result of the X-ray flux; (iii) LyC-opaque inhomogeneities induce a line of sight (LOS) scatter in $R$ as high as few physical Mpc, consistent with the DLA scenario proposed to explain the anomalous size of the ULAS J1120+0641 ionized region. On the other hand, with an ionization rate $\dot{N}_{\gamma,0} \sim 10^{57}$~s$^{-1}$, the assumed DLA clustering and gas opacity, only one LOS shows an HII region compatible with the observed one. We deduce that either the ionization rate of the QSO is at least one order of magnitude lower or the ULAS J1120+0641 bright phase is shorter than $10^7$~yrs.