ATG for Intra RH Add Endocyclic family #578
Open
+46
−1,283
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
This PR brings a few advancements to the
Intra_RH_Add_Endocyclic
reaction family. The discussion is motivated by this reactionC=CCC[CH]C=C>>[CH2]C1CCCC=C1
that had matched this templateCurrently, this family has no training reactions. Instead, a previously defined heuristic rule is used to estimate the rate for all reactions that match this family. The A-factor was arbitrarily set as 1e10 and the barrier was arbitrarily set to 5 kcal/mol. In my opinion, this rule at the root node is a gross overestimate as shown by the data below.
For this reaction, k(T=673K) = 4.757e8 1/s as seen on the RMG website
However, this reactant can also participate in an
Intra_R_Add_Exocyclic
reaction which has k(T=673K) = 8.48e4 1/s. The TST calculation was done in PR #570. This reaction actually has the radical attacking the double bound and since radicals are quite reactive, I would anticipate this reaction to be faster than theIntra_RH_Add_Endocyclic
family. In an attempt to slow down the reaction rate of this family, I propose changing the Ea from 5 kcal/mol to 20 kcal/mol which would make k(T=673K) = 6.40e3 1/s. 20 kcal/mol is also an arbitrary choice so I'm open to suggestions. For what it's worth, changing the rate here probably has only minimal impact to our users since this family is not listed as one of the default families that people would likely be using.Intra_RH_Add_Endo_rxn1.zip