Rename ketoenol to 1,3 sigmatropic rearrangement (Part 3)#532
Conversation
|
Although this branch was based on main locally, I must've accidentally pushed to master so this PR includes commits that were already merged in by @xiaoruiDong in #527. Sorry about that. I've now updated my settings to push to main by default. |
d49d1d2 to
9d3d5d8
Compare
amarkpayne
left a comment
There was a problem hiding this comment.
The first 5 commits here look good to me once @kspieks updates the typo in the file name he found. @xiaoruiDong wants the PR to include a commit that forbids a structure. Once @kspieks includes that I think just that 1 commit can be reviewed and then this can be merged in
Previous template: R1=R2-O3-R4 => R4-R1-R2=O3 or R1=R2-S3-R4 => R4-R1-R2=S3 New template: R1=R2-R3-R4 => R4-R1-R2=R3
Thanks to Xiaorui for recommending this paper: https://pubs.acs.org/doi/full/10.1021/jp0212545 Values from 3-parameter Arrhenius fits are taken from the high-pressure limit rate constants in Table 4
Copies the .gitignore from RMG-Py, ARC, and T3
9d3d5d8 to
5661573
Compare
|
@amarkpayne thanks for reviewing @xiaoruiDong and I discussed offline that it was better to try finding an example reaction rather than forbid the group. Thanks to Xiaorui for finding this paper. I added a training reaction for the group we discussed and as expected, it has a very large barrier (which is in agreement with the Wikipedia description that 1,3 sigmatropic reactions with a C=C will be nearly impossible). This PR should be ready for review :) Here is the new notebook for fitting the rate trees with the additional training reactions: |
xiaoruiDong
left a comment
There was a problem hiding this comment.
@kspieks Thanks, Kevin! I was hoping to slightly improve the rate estimation for [1,3]-sigmatropic H migration. However, I just found the new addition from the paper only improves the entry for reactions involving a ring. Root_N-1R!H-inRing still looks terrible though. I guess that's all we can do for now.
For future reference, we don't forbid R!H=R!H-C-H is because we want this reaction to appear in the mechanism (i.e., rate corresponding to the antarafacial TS may be significant). If we simply forbid this reaction, it will never show up, never get into people's attention, and never get refined.
As described in #524 and #527, this is the third PR in a three-part series for adding reactions to and updating the ketoenol family.
This PR updates the reaction template to be more generalized.
Previous template: R1=R2-O3-R4 => R4-R1-R2=O3 or R1=R2-S3-R4 => R4-R1-R2=S3
New template: R1=R2-R3-R4 => R4-R1-R2=R3
This PR also adds 7 new training reactions that match the new template and refit the rate tree.
The zip file below contains the 3-parameter Arrhenius fits from Arkane, the notebook to create the training reactions, and the notebook to refit the rate tree.
Note that some of the Arrhenius fits have fairly large uncertainty in the parameters, with the worst case being
dA = *|/ 656.605, dn = +|- 0.860771, dEa = +|- 4.43858 kJ/molfor rxn006539. Despite the uncertainty in the paramters, the average percentage error across all 50 TST data points used for fitting is 130% i.e. only about a factor of 2 error in the rate constant, which is often quite acceptable. At the other end of the spectrum, rxn009176 had the least uncertainty in the parameters from fitting and also the lowest average percentage error in the rate constant (just 11% error across all 50 points).20211010_1,3_sigmatropic_rearrangement.zip
Edit: Do NOT use this^ zip file. Instead, use the updated one below!