Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

June 2019 Training reaction fixing hackathon #340

Open
wants to merge 36 commits into
base: main
Choose a base branch
from
Open
Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
36 commits
Select commit Hold shift + click to select a range
4b6b19d
H_abstraction tree preparation
mjohnson541 Aug 15, 2018
466aec2
Database Label Fix: H_Abstraction, NH3 and NH2
oscarwumit Apr 25, 2019
fb102c6
Database Label Fix: H_Abstraction, CH4 and CH3
oscarwumit Apr 25, 2019
02d175b
Database Label Fix: H_Abstraction, H2O and OH
oscarwumit Apr 25, 2019
069c2ae
Database Label Fix: Disproportionation, CH3
oscarwumit Apr 25, 2019
14622d3
nitrogen dictionary fixes
mjohnson541 May 23, 2019
bccb8d3
reaction fixes H_Abstraction
mjohnson541 May 23, 2019
fdc87ac
fix entries for H_abs for tree generation
mjohnson541 May 27, 2019
703cd36
fix issue with use of CH3 vs CH3_r1 in species dictionary for disprop…
mjohnson541 May 27, 2019
7951f23
Script to identify a list of reactions including indices that produce…
ajocher Jun 12, 2019
b980c13
fixup! Script to identify a list of reactions including indices that …
ajocher Jun 13, 2019
bd528a3
add collision limit analysis and display bimolecular reactions
mjohnson541 Jun 13, 2019
78edd84
print collision limit rate
mjohnson541 Jun 13, 2019
6cc0c6f
Fixes to kineticsTrainingAnalysis
mliu49 Jun 13, 2019
40aa6ad
add notebook for finding duplicates
mjohnson541 Jun 13, 2019
1ff7658
Another fix to notebook
mliu49 Jun 13, 2019
af15267
Another fix
mliu49 Jun 13, 2019
bfc312e
Remove Training Reactions that Violate Collision Limit in Substitution_O
amarkpayne Jun 13, 2019
905f7fc
Remove Collision Limit Violators in intra_substitutionS_cyclization
amarkpayne Jun 13, 2019
328dcf2
deleted reaction 824 from Intra_R_Add_Endocyclic
yunsiechung Jun 13, 2019
e7c5578
Remove bad estimates from Intra_R_Add_Exocyclic training reactions
mliu49 Jun 13, 2019
00f0734
Fix species for Intra_R_Add_Exocyclic training rxn from rate rule
mliu49 Jun 13, 2019
b3d14d9
Remove Duplicate Reaction for R_Addition_COm
amarkpayne Jun 13, 2019
07b3324
Removed duplicate training reaction for R_Addition_CSm.
ajocher Jun 14, 2019
6642630
Updated reference and A factor for reaction 7 in 1+2_Cycloaddition.
ajocher Jun 14, 2019
ef9fcef
Remove Collision Limit Violator in intra_H_migration.
ajocher Jun 14, 2019
b8cf8f9
Commented out training reactions including species with thermo data t…
ajocher Jun 14, 2019
a1cf72f
Commented out training reactions including species with thermo data t…
ajocher Jun 14, 2019
014c467
Commented out training reactions including species with thermo data t…
ajocher Jun 14, 2019
eca12ec
Removed one more training reaction that violate collision limit in Su…
ajocher Jun 14, 2019
b375770
Remove duplicate reactions for R_Addition_MultipleBond, mainly with r…
ajocher Jun 14, 2019
612c575
Remove duplicate reactions for intra_H_migration, mainly with rank 11…
ajocher Jun 14, 2019
4c3a990
H_Abstraction maximum rate testing and duplicate removal
mjohnson541 Jun 16, 2019
ea555ff
fix i-butanol training reactions
mjohnson541 Jun 16, 2019
11ec581
remove GAV based training reactions
mjohnson541 Jun 16, 2019
41fa5ee
Remove Duplicates and Rank 0 Small GA entries from Intra_R_Add_Endocyc.
amarkpayne Jun 17, 2019
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Original file line number Diff line number Diff line change
Expand Up @@ -149,18 +149,19 @@
label = "CH2 + C2H2 <=> C3H4",
degeneracy = 1.0,
kinetics = Arrhenius(
A = (1.77e+15, 'cm^3/(mol*s)'),
A = (1.06e+15, 'cm^3/(mol*s)'),
n = -0.662,
Ea = (0.157737, 'kJ/mol'),
T0 = (1, 'K'),
Tmin = (200, 'K'),
Tmax = (2000, 'K'),
),
rank = 5,
shortDesc = u"""Polino [carbene,acetylene]""",
shortDesc = u"""Polino et al. J. Phys. Chem. A 2013, 117, 12677−12692.""",
longDesc =
u"""
Converted to training reaction from rate rule: carbene;mb_tb_unsub
Daniela Polino, Stephen J. Klippenstein, Lawrence B. Harding, and Yuri Georgievskii,Predictive Theory for the Addition and Insertion Kinetics of 1CH2
Reacting with Unsaturated Hydrocarbons, J. Phys. Chem. A 2013, 117, 12677−12692.
""",
)

Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -80,6 +80,13 @@ multiplicity 2
3 H u0 p0 c0 {1,S}
4 *4 H u0 p0 c0 {1,S}

CH3_r1
multiplicity 2
1 *1 C u1 p0 c0 {2,S} {3,S} {4,S}
2 H u0 p0 c0 {1,S}
3 H u0 p0 c0 {1,S}
4 H u0 p0 c0 {1,S}

H
multiplicity 2
1 *1 H u1 p0 c0
Expand Down
28 changes: 14 additions & 14 deletions input/kinetics/families/Disproportionation/training/reactions.py
Original file line number Diff line number Diff line change
Expand Up @@ -159,7 +159,7 @@

entry(
index = 5,
label = "CH3 + C2H5 <=> CH4 + C2H4",
label = "CH3_r1 + C2H5 <=> CH4 + C2H4",
degeneracy = 3.0,
kinetics = Arrhenius(
A = (6.57e+14, 'cm^3/(mol*s)', '*|/', 1.1),
Expand Down Expand Up @@ -636,7 +636,7 @@

entry(
index = 20,
label = "CH3 + C3H7-2 <=> CH4 + C3H6-2",
label = "CH3_r1 + C3H7-2 <=> CH4 + C3H6-2",
degeneracy = 2.0,
kinetics = Arrhenius(
A = (2.3e+13, 'cm^3/(mol*s)', '*|/', 1.7),
Expand Down Expand Up @@ -1151,7 +1151,7 @@

entry(
index = 37,
label = "CH3 + C4H9-2 <=> CH4 + C4H8-2",
label = "CH3_r1 + C4H9-2 <=> CH4 + C4H8-2",
degeneracy = 1.0,
kinetics = Arrhenius(
A = (6.02e+12, 'cm^3/(mol*s)', '*|/', 2),
Expand Down Expand Up @@ -1457,7 +1457,7 @@

entry(
index = 46,
label = "CH3 + C3H5-2 <=> CH4 + C3H4",
label = "CH3_r1 + C3H5-2 <=> CH4 + C3H4",
degeneracy = 1.0,
kinetics = Arrhenius(
A = (3.01e+12, 'cm^3/(mol*s)', '*|/', 3),
Expand Down Expand Up @@ -1814,7 +1814,7 @@

entry(
index = 57,
label = "CH3O + CH3 <=> CH4 + CH2O",
label = "CH3O + CH3_r1 <=> CH4 + CH2O",
degeneracy = 1.0,
kinetics = Arrhenius(
A = (8.49e+13, 'cm^3/(mol*s)'),
Expand Down Expand Up @@ -2606,7 +2606,7 @@

entry(
index = 88,
label = "H2N2 + CH3 <=> CH4 + HN2",
label = "H2N2 + CH3_r1 <=> CH4 + HN2",
degeneracy = 2.0,
kinetics = Arrhenius(
A = (3.2e+06, 'cm^3/(mol*s)'),
Expand Down Expand Up @@ -2738,7 +2738,7 @@

entry(
index = 94,
label = "H3N2 + CH3 <=> CH4 + H2N2-2",
label = "H3N2 + CH3_r1 <=> CH4 + H2N2-2",
degeneracy = 2.0,
kinetics = Arrhenius(
A = (1.64e+06, 'cm^3/(mol*s)'),
Expand Down Expand Up @@ -2894,7 +2894,7 @@

entry(
index = 101,
label = "H2NO + CH3 <=> CH4 + HNO-2",
label = "H2NO + CH3_r1 <=> CH4 + HNO-2",
degeneracy = 2.0,
kinetics = Arrhenius(
A = (3.2e+06, 'cm^3/(mol*s)'),
Expand Down Expand Up @@ -3050,7 +3050,7 @@

entry(
index = 108,
label = "H2NO-3 + CH3 <=> CH4 + HNO-3",
label = "H2NO-3 + CH3_r1 <=> CH4 + HNO-3",
degeneracy = 1.0,
kinetics = Arrhenius(
A = (1.6e+06, 'cm^3/(mol*s)'),
Expand Down Expand Up @@ -3186,7 +3186,7 @@

entry(
index = 114,
label = "CH3 + CH2N <=> CH4 + CHN",
label = "CH3_r1 + CH2N <=> CH4 + CHN",
degeneracy = 2.0,
kinetics = Arrhenius(
A = (1.62e+06, 'cm^3/(mol*s)'),
Expand Down Expand Up @@ -3362,7 +3362,7 @@

entry(
index = 122,
label = "CH2N-3 + CH3 <=> CH4 + CHN-2",
label = "CH2N-3 + CH3_r1 <=> CH4 + CHN-2",
degeneracy = 1.0,
kinetics = Arrhenius(
A = (820000, 'cm^3/(mol*s)'),
Expand Down Expand Up @@ -3450,7 +3450,7 @@

entry(
index = 126,
label = "CH4N + CH3 <=> CH4 + CH3N-2",
label = "CH4N + CH3_r1 <=> CH4 + CH3N-2",
degeneracy = 3.0,
kinetics = Arrhenius(
A = (7.2e+06, 'cm^3/(mol*s)'),
Expand Down Expand Up @@ -3538,7 +3538,7 @@

entry(
index = 130,
label = "CH3 + CH4N-2 <=> CH4 + CH3N-3",
label = "CH3_r1 + CH4N-2 <=> CH4 + CH3N-3",
degeneracy = 2.0,
kinetics = Arrhenius(
A = (3.2e+06, 'cm^3/(mol*s)'),
Expand Down Expand Up @@ -3607,7 +3607,7 @@

entry(
index = 133,
label = "CH3 + CH2NO <=> CH4 + CHNO",
label = "CH3_r1 + CH2NO <=> CH4 + CHNO",
degeneracy = 2.0,
kinetics = Arrhenius(
A = (3.2e+06, 'cm^3/(mol*s)'),
Expand Down
9 changes: 7 additions & 2 deletions input/kinetics/families/H_Abstraction/groups.py
Original file line number Diff line number Diff line change
Expand Up @@ -21,14 +21,19 @@
entry(
index = 0,
label = "X_H_or_Xrad_H_Xbirad_H_Xtrirad_H",
group = "OR{Xtrirad_H, Xbirad_H, Xrad_H, X_H}",
group = """
1 *1 R u[0,1,2,3] {2,S}
2 *2 H u0 {1,S}
""",
kinetics = None,
)

entry(
index = 1,
label = "Y_rad_birad_trirad_quadrad",
group = "OR{Y_rad, Y_1centerbirad, Y_1centertrirad, Y_1centerquadrad}",
group = """
1 *3 R u[1,2,3,4]
""",
kinetics = None,
)

Expand Down
119 changes: 78 additions & 41 deletions input/kinetics/families/H_Abstraction/training/dictionary.txt
Original file line number Diff line number Diff line change
Expand Up @@ -251,20 +251,6 @@ CH3N-2
4 H u0 p0 c0 {2,S}
5 H u0 p0 c0 {1,S}

CH3_p1
multiplicity 2
1 *1 C u1 p0 c0 {2,S} {3,S} {4,S}
2 H u0 p0 c0 {1,S}
3 H u0 p0 c0 {1,S}
4 H u0 p0 c0 {1,S}

CH3_r3
multiplicity 2
1 *3 C u1 p0 c0 {2,S} {3,S} {4,S}
2 H u0 p0 c0 {1,S}
3 H u0 p0 c0 {1,S}
4 H u0 p0 c0 {1,S}

C5H11
multiplicity 2
1 C u0 p0 c0 {2,S} {3,S} {4,S} {5,S}
Expand Down Expand Up @@ -1164,14 +1150,21 @@ H3NO-2
4 H u0 p0 c0 {2,S}
5 *2 H u0 p0 c0 {1,S}

CH4b
CH4
1 C u0 p0 c0 {2,S} {3,S} {4,S} {5,S}
2 H u0 p0 c0 {1,S}
3 H u0 p0 c0 {1,S}
4 H u0 p0 c0 {1,S}
5 H u0 p0 c0 {1,S}

CH4_r12
1 *1 C u0 p0 c0 {2,S} {3,S} {4,S} {5,S}
2 *2 H u0 p0 c0 {1,S}
3 H u0 p0 c0 {1,S}
4 H u0 p0 c0 {1,S}
5 H u0 p0 c0 {1,S}

CH4p
CH4_p23
1 *3 C u0 p0 c0 {2,S} {3,S} {4,S} {5,S}
2 *2 H u0 p0 c0 {1,S}
3 H u0 p0 c0 {1,S}
Expand Down Expand Up @@ -1285,13 +1278,6 @@ multiplicity 2
4 H u0 p0 c0 {1,S}
5 H u0 p0 c0 {2,S}

CH4
1 C u0 p0 c0 {2,S} {3,S} {4,S} {5,S}
2 H u0 p0 c0 {1,S}
3 H u0 p0 c0 {1,S}
4 H u0 p0 c0 {1,S}
5 H u0 p0 c0 {1,S}

C3H7O2
multiplicity 2
1 O u0 p2 c0 {2,S} {3,S}
Expand Down Expand Up @@ -2075,11 +2061,16 @@ H2S
2 *2 H u0 p0 c0 {1,S}
3 H u0 p0 c0 {1,S}

H2O
H2O_r12
1 *1 O u0 p2 c0 {2,S} {3,S}
2 *2 H u0 p0 c0 {1,S}
3 H u0 p0 c0 {1,S}

H2O_p23
1 *3 O u0 p2 c0 {2,S} {3,S}
2 *2 H u0 p0 c0 {1,S}
3 H u0 p0 c0 {1,S}

C4H10b
1 *1 C u0 p0 c0 {2,S} {3,S} {5,S} {6,S}
2 C u0 p0 c0 {1,S} {4,S} {7,S} {8,S}
Expand Down Expand Up @@ -2413,6 +2404,14 @@ multiplicity 2
5 H u0 p0 c0 {3,S}
6 H u0 p0 c0 {2,S}

H3N2
multiplicity 2
1 N u0 p1 c0 {2,S} {4,S} {5,S}
2 *3 N u1 p1 c0 {1,S} {3,S}
3 H u0 p0 c0 {2,S}
4 H u0 p0 c0 {1,S}
5 H u0 p0 c0 {1,S}

C6H6-2
1 *1 C u0 p0 c0 {2,S} {3,S} {4,S} {7,S}
2 C u0 p0 c0 {1,S} {8,S} {9,S} {10,S}
Expand Down Expand Up @@ -2594,6 +2593,27 @@ C4H8-8
11 H u0 p0 c0 {4,S}
12 H u0 p0 c0 {4,S}

CH3_r12
multiplicity 2
1 *1 C u1 p0 c0 {2,S} {3,S} {4,S}
2 H u0 p0 c0 {1,S}
3 H u0 p0 c0 {1,S}
4 *2 H u0 p0 c0 {1,S}

CH3_r3
multiplicity 2
1 *3 C u1 p0 c0 {2,S} {3,S} {4,S}
2 H u0 p0 c0 {1,S}
3 H u0 p0 c0 {1,S}
4 H u0 p0 c0 {1,S}

CH3_p1
multiplicity 2
1 *1 C u1 p0 c0 {2,S} {3,S} {4,S}
2 H u0 p0 c0 {1,S}
3 H u0 p0 c0 {1,S}
4 H u0 p0 c0 {1,S}

CH3_p23
multiplicity 2
1 *3 C u1 p0 c0 {2,S} {3,S} {4,S}
Expand Down Expand Up @@ -3243,12 +3263,6 @@ NH3_p23
3 H u0 p0 c0 {1,S}
4 H u0 p0 c0 {1,S}

NH3_r
1 *1 N u0 p1 c0 {2,S} {3,S} {4,S}
2 *2 H u0 p0 c0 {1,S}
3 H u0 p0 c0 {1,S}
4 H u0 p0 c0 {1,S}

N2H2(T)_r3
multiplicity 3
1 *3 N u1 p1 c0 {2,S} {3,S}
Expand Down Expand Up @@ -3332,13 +3346,6 @@ multiplicity 3
3 H u0 p0 c0 {1,S}
4 H u0 p0 c0 {1,S}

H2NN(T)_p1
multiplicity 3
1 N u0 p1 c0 {2,S} {3,S} {4,S}
2 *1 N u2 p1 c0 {1,S}
3 H u0 p0 c0 {1,S}
4 H u0 p0 c0 {1,S}

CH3
multiplicity 2
1 *1 C u1 p0 c0 {2,S} {3,S} {4,S}
Expand Down Expand Up @@ -3858,6 +3865,11 @@ multiplicity 2
1 *3 O u1 p2 c0 {2,S}
2 *2 H u0 p0 c0 {1,S}

OH_r3
multiplicity 2
1 *3 O u1 p2 c0 {2,S}
2 H u0 p0 c0 {1,S}

SH_p1
multiplicity 2
1 *1 S u1 p2 c0 {2,S}
Expand Down Expand Up @@ -4000,10 +4012,17 @@ H2N2
3 *2 H u0 p0 c0 {1,S}
4 H u0 p0 c0 {2,S}

OH
multiplicity 2
1 *3 O u1 p2 c0 {2,S}
2 H u0 p0 c0 {1,S}
H2N2-2
1 *3 N u0 p1 c0 {2,D} {3,S}
2 N u0 p1 c0 {1,D} {4,S}
3 *2 H u0 p0 c0 {1,S}
4 H u0 p0 c0 {2,S}

H2N2-3
1 N u0 p1 c0 {2,S} {3,S} {4,S}
2 *1 N u2 p1 c0 {1,S}
3 H u0 p0 c0 {1,S}
4 H u0 p0 c0 {1,S}

C2H5OS
multiplicity 2
Expand Down Expand Up @@ -4475,6 +4494,24 @@ C6H12O2-4
19 H u0 p0 c0 {7,S}
20 H u0 p0 c0 {7,S}

NH2
multiplicity 2
1 N u1 p1 c0 {2,S} {3,S}
2 H u0 p0 c0 {1,S}
3 H u0 p0 c0 {1,S}

NH3
1 N u0 p1 c0 {2,S} {3,S} {4,S}
2 H u0 p0 c0 {1,S}
3 H u0 p0 c0 {1,S}
4 H u0 p0 c0 {1,S}

NH3_r12
1 *1 N u0 p1 c0 {2,S} {3,S} {4,S}
2 *2 H u0 p0 c0 {1,S}
3 H u0 p0 c0 {1,S}
4 H u0 p0 c0 {1,S}

C3H6S
1 S u0 p2 c0 {4,D}
2 *1 C u0 p0 c0 {3,S} {4,S} {5,S} {6,S}
Expand Down
Loading