Modeling the neutral phases of the interstellar medium in galaxy simulations | Synapse