Monte Carlo (MC) simulations in ensembles with a fixed chemical potential or fugacity, for example the grand-canonical or the osmotic ensemble, are often used to compute phase equilibria. Chemical potentials can be computed either with an equation of state (EoS) or from molecular simulations. The accuracy of the computed chemical potentials depends on the quality of the (critical) parameters used in the EoS and the applied force field in the simulations. We investigated the consistency of both approaches for computing fugacities of the industrially relevant gases CO 2 , CH 4 , CO, H 2 , N 2 , and H 2 S. The critical temperature (T c), pressure (P c), and acentric factors (ω) of these gases are computed from MC simulations in the Gibbs ensemble. The effect of cutoff radius and tail corrections on the computed values of T c , P c , and ω is investigated. In addition, MC simulations in the Gibbs ensemble are used to compute the VLE of the 15 possible binary systems comprising the gases CO 2 , CH 4 , CO, H 2 , N 2 , and H 2 S, and the ternary systems CO 2 /CH 4 /H 2 S and CO 2 /CO/H 2. Binary interaction parameters (k ij) of these natural/synthesis gas mixtures are obtained by fitting the Peng-Robinson (PR) EoS to the binary VLE data from the MC simulations. The computed properties from the MC simulations are compared with the PR EoS, the GERG EoS, and experimental results. The MC results show that including tail corrections in the simulations is crucial to obtain accurate critical properties. The force fields used for the gases can reproduce the fugacities of the gases within 5% of