The analysis presented in this study demonstrates the capability of the HhPF method with multi-Gaussian filters to accurately model liquid-gas coexistence across a wide range of thermodynamic conditions, capturing both macroscopic phase behavior and microscopic structural details. We have successfully integrated multi-Gaussian filters into the Hamiltonian hybrid particle-field (HhPF) methodology, enhancing its capability to address phase coexistence phenomena. This advancement is analogous to the evolution from GCM to MGCM and has been implemented within the HyMD simulation package. By utilizing a linear combination of two Gaussian filters, we demonstrate that the HhPF approach can produce potentials with attractive components similar to Lennard-Jones potentials, which are essential for modeling phase coexistence. This development extends the applicability of the hybrid particle-field methodology to phase coexistence studies without the need for modifications to the equation of state or the introduction of additional parameters.