Bimodal ion energy distribution functions in a hydrogen magnetized plasma