Simulating polymorphic phase behavior using reaction ensemble Monte Carlo