Phase Diagrams of Diatomic Molecules Using the Gibbs Ensemble Monte Carlo Method