Helicon plasma sources operating with hydrogen or deuterium might be attractive for fusion applications due to their higher power efficiency compared to inductive radiofrequency plasma sources. During recent years, the Resonant Antenna Ion Device (RAID) has been investigating the physics of helicon plasmas and the possibility of employing them to produce negative ions for Heating Neutral Beam injectors (HNBs). We present herein a fluid-Monte Carlo model describing plasma species transport of a typical helicon hydrogen plasma discharge. This work is motivated by the interest to better understand the basic physics of helicon plasma devices when operating in hydrogen and, in particular, the volume production of negative ions. This model is based on the synergy between two separately self-consistent approaches: a plasma fluid model calculating ion transport, and a Monte Carlo (MC) model, to determine neutral and rovibrational density profiles of H2. By introducing as model constraints the electron density and temperature profiles measured by Langmuir Probes, the densities of ion species (H+, H2+ , H3+ , H−) are computed in a 1.5D (dimensional) geometry. The estimate of the negative ion density profile represents a useful benchmark to be compared with dedicated diagnostics such as Cavity Ring-Down Spectroscopy and Langmuir probe laser photodetachment. Neutral gas particles (atoms and molecules) are calculated assuming a fixed plasma background. This gas-plasma decoupling is necessary due to the different timescales of plasma (microseconds) and gas kinetics (milliseconds).