Skip to content

Commit 3d8fb22

Browse files
committed
Added sulfur group to ketoenol Family
1 parent caf4827 commit 3d8fb22

2 files changed

Lines changed: 33 additions & 6 deletions

File tree

input/kinetics/families/ketoenol/groups.py

Lines changed: 32 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -4,10 +4,12 @@
44
name = "ketoenol/groups"
55
shortDesc = u""
66
longDesc = u"""
7-
7+
Sulfur was added to this family, and is treated the same as oxygen.
8+
Ideally we would like to branch this into a new family "R=RSR <=> RRR=S"
9+
once relevant kinetic data is available
810
"""
911

10-
template(reactants=["R_ROR"], products=["keton"], ownReverse=False)
12+
template(reactants=["R_ROSR"], products=["keton"], ownReverse=False)
1113

1214
reverse = "none"
1315

@@ -50,6 +52,15 @@
5052
kinetics = None,
5153
)
5254

55+
56+
entry(
57+
index = 500,
58+
label = "R_ROSR",
59+
group = "OR{R_ROR, R_RSR}",
60+
kinetics = None,
61+
)
62+
63+
5364
entry(
5465
index = 1,
5566
label = "R_ROR",
@@ -63,6 +74,20 @@
6374
kinetics = None,
6475
)
6576

77+
78+
entry(
79+
index = 501,
80+
label = "R_RSR",
81+
group =
82+
"""
83+
1 *1 R!H u0 {2,D}
84+
2 *2 R!H u0 {1,D} {3,S}
85+
3 *3 S2s u0 {2,S} {4,S}
86+
4 *4 R u0 {3,S}
87+
""",
88+
kinetics = None,
89+
)
90+
6691
entry(
6792
index = 22,
6893
label = "R_O_H",
@@ -110,9 +135,9 @@
110135
label = "R1_doublebond_CHR",
111136
group =
112137
"""
113-
1 *1 C u0 {2,S} {3,S}
138+
1 *1 C u0 {2,S} {3,S}
114139
2 R!H u0 {1,S}
115-
3 H u0 {1,S}
140+
3 H u0 {1,S}
116141
""",
117142
kinetics = None,
118143
)
@@ -209,7 +234,9 @@
209234

210235
tree(
211236
"""
212-
L1: R_ROR
237+
L1: R_ROSR
238+
L2: R_ROR
239+
L2: R_RSR
213240
L1: R1_doublebond
214241
L2: R1_doublebond_CHR
215242
L3: R1_doublebond_CHCH3

input/kinetics/families/ketoenol/rules.py

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -8,7 +8,7 @@
88
"""
99
entry(
1010
index = 1,
11-
label = "R_ROR;R1_doublebond;R2_doublebond;R_O",
11+
label = "R_ROSR;R1_doublebond;R2_doublebond;R_O",
1212
kinetics = ArrheniusEP(
1313
A = (100000, 's^-1'),
1414
n = 2,

0 commit comments

Comments
 (0)