Predicting ionised gas emission in 3D with SKIRT. I. Framework and validation

Emission lines from ionised gas are key diagnostics of star formation, metallicity, and ionisation conditions in galaxies. Interpreting spatially resolved observations from integral-field surveys (e.g. MaNGA, MUSE, JWST/NIRSpec) and comparing them with hydrodynamical simulations requires 3D photoionisation models that handle realistic geometries, dust attenuation, and synthetic instrument output. We present a new photoionisation module for the Monte Carlo radiative transfer code SKIRT that predicts emission-line luminosities of ionised gas in 3D, combining pre-computed Cloudy tables for gas temperature and opacity with a direct calculation of ion fractions and line emissivities. The local ionising radiation field (1-6 Ryd) is characterised by log U and four spectral-shape ratios; Cloudy tables map these to temperature and opacity, converging through SKIRT's existing iteration cycle. An inline solver then determines ion fractions from the converged field and temperature and evaluates line emissivities. We validate against Cloudy on 60 spherical shell models and against COLT on a Milky Way-analogue galaxy. On the 1D grid, hydrogen recombination lines agree with Cloudy to within a few per cent (Halpha median ratio 0.97) and the forbidden lines to within 5%, except [S II] 6717 (1.23), whose offset traces a temperature overestimate near the ionisation front. In 3D, integrated luminosities agree with COLT to within 18% for the hydrogen lines and 2% for [N II], while [O III] and [S II] are elevated by 70 and 80%. Pixel-by-pixel correlation coefficients reach r >= 0.92, with luminosity-weighted scatter of 0.14-0.31 dex and broadly consistent BPT ratios. The module enables self-consistent synthetic observations in which ionised-gas emission lines, dust attenuation, and dust re-emission are computed in a single MCRT run, applicable to any hydrodynamical simulation.