Abstract The present work develops a full particle-based model that couples the particle-in-cell plus Monte Carlo collision (PIC-MCC) simulation for plasma dynamics and the direct simulation Monte Carlo (DSMC) method for neutral dynamics in a synergistic iterative manner. This new model overcomes the slow convergence issue in the conventional direct coupling approach caused by the disparity of the time scales between the plasma and neutral dynamics. This model is applied to simulate the behavior of xenon (Xe) and its potential alternatives, krypton (Kr) and argon (Ar), in the discharge chamber of a miniature direct current (DC) ion thruster. The results show that a stable discharge is difficult to achieve for Kr and Ar under the operating conditions optimal for Xe. While increasing the discharge voltage can effectively improve the stability of discharge for Kr and Ar, other common strategies such as changing the magnetic field strength, propellant flow rate, and cathode current are not successful. The propellant utilization efficiency and discharge efficiency are affected by both discharge voltage and propellant flow rate. A maximum utilization efficiency and an optimal discharge efficiency are observed for all three propellants, with the values decreasing in the order of Xe, Kr, and Ar. Moreover, the discharge voltage corresponding to the optimal efficiency is inversely proportional to the square root of the propellant mass, indicating that the ion diffusional loss to the wall, rather than the ionization energy, is the dominant factor affecting the discharge performance for alternative propellants in a miniature DC thruster.