check kinetics for pdep reactions before adding to core

rmgpy/rmg/model.py

@@ -1929,6 +1929,13 @@ def update_unimolecular_reaction_networks(self):
         while index < core_reaction_count:
             reaction = self.core.reactions[index]
             if isinstance(reaction, PDepReaction):
+                if np.isnan(reaction.kinetics.get_rate_coefficient(1000.0,1e5)):
+                    # this reaction kinetics are undefined, dont add to the core
+                    logging.warning(f"pdep reaction {reaction} has invalid kinetics and will be removed from the core: {reaction.kinetics}")
+                    self.core.reactions.remove(reaction)
+                    core_reaction_count -= 1
+                    index += 1
+                    continue
                 for reaction2 in self.core.reactions[index + 1:]:
                     if isinstance(reaction2,
                                   PDepReaction) and reaction.reactants == reaction2.products and reaction.products == reaction2.reactants: