Renamed Oa_R_Recombination family as Birad_R_Recombination

Elaborated with [S] and [NH]
Added training reactions for [O], [S], [NH]

Author: alongd

…s/Oa_R_Recombination/NIST/dictionary.txt …irad_R_Recombination/NIST/dictionary.txtinput/kinetics/families/Oa_R_Recombination/NIST/dictionary.txt renamed to input/kinetics/families/Birad_R_Recombination/NIST/dictionary.txt input/kinetics/families/Oa_R_Recombination/NIST/dictionary.txt renamed to input/kinetics/families/Birad_R_Recombination/NIST/dictionary.txt

@@ -1,16 +1,16 @@
 #!/usr/bin/env python
 # encoding: utf-8
 
-name = "Oa_R_Recombination/groups"
+name = "Birad_R_Recombination/groups"
 shortDesc = u""
 longDesc = u"""
-This reaction family is reserved for O_atom (triplet only). The forbidden groups
-at the bottom prevent it from reacting with other forms of O.
+This reaction family is reserved for recombination of O_atom, S_atom, N_R_birad (triplets only).
+The forbidden groups at the bottom prevent it from reacting with other forms of O, S, NH.
 """
 
-template(reactants=["Y_rad", "Oa"], products=["YO."], ownReverse=False)
+template(reactants=["Y_rad", "Birad"], products=["YOS."], ownReverse=False)
 
-reverse = "RO_Bond_Dissociation"
+reverse = "ROS_Bond_Dissociation"
 
 recipe(actions=[
     ['FORM_BOND', '*1', 1, '*2'],
@@ -30,10 +30,95 @@
 
 entry(
     index = 2,
+    label = "Birad",
+    group = 
+"""
+1 *2 R!H u2
+""",
+    kinetics = None,
+)
+
+entry(
+    index = 50,
     label = "Oa",
     group = 
 """
-1 *2 O u2
+1 *2 O u2 p2
+""",
+    kinetics = None,
+)
+
+entry(
+    index = 51,
+    label = "Sa",
+    group = 
+"""
+1 *2 S u2 p2
+""",
+    kinetics = None,
+)
+
+entry(
+    index = 52,
+    label = "N_R_birad",
+    group = 
+"""
+1 *2 N u2 p1
+""",
+    kinetics = None,
+)
+
+entry(
+    index = 53,
+    label = "N_birad/H",
+    group = 
+"""
+1 *2 N u2 p1 {2,S}
+2    H u0 p0 {1,S}
+""",
+    kinetics = None,
+)
+
+entry(
+    index = 54,
+    label = "N_birad/C",
+    group = 
+"""
+1 *2 N u2 p1 {2,S}
+2    C ux    {1,S}
+""",
+    kinetics = None,
+)
+
+entry(
+    index = 55,
+    label = "N_birad/O",
+    group = 
+"""
+1 *2 N u2 p1 {2,S}
+2    O ux    {1,S}
+""",
+    kinetics = None,
+)
+
+entry(
+    index = 56,
+    label = "N_birad/N",
+    group = 
+"""
+1 *2 N u2 p1 {2,S}
+2    N ux    {1,S}
+""",
+    kinetics = None,
+)
+
+entry(
+    index = 57,
+    label = "N_birad/S",
+    group = 
+"""
+1 *2 N u2 p1 {2,S}
+2    S ux    {1,S}
 """,
     kinetics = None,
 )
@@ -673,7 +758,15 @@
                 L5: C_rad/Cs
                 L5: C_rad/TDMustO
             L4: C_rad/ThreeDe
-L1: Oa
+L1: Birad
+    L2: Oa
+    L2: Sa
+    L2: N_R_birad
+        L3: N_birad/H
+        L3: N_birad/C
+        L3: N_birad/O
+        L3: N_birad/N
+        L3: N_birad/S
 """
 )
 
@@ -688,6 +781,20 @@
     longDesc = 
 u"""
 
+""",
+)
+ 
+forbidden(
+    label = "S2_1centeredBirad",
+    group = 
+"""
+1 *1 S u2 p1 {2,S}
+2    R u0 {1,S}
+""",
+    shortDesc = u"""""",
+    longDesc = 
+u"""
+
 """,
 )
 
@@ -702,6 +809,20 @@
     longDesc = 
 u"""
 
+""",
+)
+ 
+forbidden(
+    label = "S2_birad",
+    group = 
+"""
+1 *1 S u1 p2 {2,S}
+2    S u1 p2 {1,S}
+""",
+    shortDesc = u"""""",
+    longDesc = 
+u"""
+
 """,
 )
 
@@ -718,3 +839,29 @@
 """,
 )
 
+forbidden(
+    label = "S_atom_singlet",
+    group = 
+"""
+1 *1 S u0 p3
+""",
+    shortDesc = u"""""",
+    longDesc = 
+u"""
+
+""",
+)
+
+forbidden(
+    label = "N_R_singlet",
+    group = 
+"""
+1 *1 N u0 p2
+""",
+    shortDesc = u"""""",
+    longDesc = 
+u"""
+
+""",
+)
+

…ies/Oa_R_Recombination/NIST/reactions.py …/Birad_R_Recombination/NIST/reactions.pyinput/kinetics/families/Oa_R_Recombination/NIST/reactions.py renamed to input/kinetics/families/Birad_R_Recombination/NIST/reactions.py input/kinetics/families/Oa_R_Recombination/NIST/reactions.py renamed to input/kinetics/families/Birad_R_Recombination/NIST/reactions.py

@@ -0,0 +1,50 @@
+#!/usr/bin/env python
+# encoding: utf-8
+
+name = "Birad_R_Recombination/rules"
+shortDesc = u""
+longDesc = u"""
+
+"""
+entry(
+    index = 1000,
+    label = "Y_rad;Oa",
+    kinetics = ArrheniusEP(
+        A = (1e+13, 'cm^3/(mol*s)'),
+        n = 0,
+        alpha = 0,
+        E0 = (0, 'kcal/mol'),
+        Tmin = (300, 'K'),
+        Tmax = (1500, 'K'),
+    ),
+    rank = 2,
+)
+
+entry(
+    index = 2000,
+    label = "Y_rad;Sa",
+    kinetics = ArrheniusEP(
+        A = (1e+13, 'cm^3/(mol*s)'),
+        n = 0,
+        alpha = 0,
+        E0 = (0, 'kcal/mol'),
+        Tmin = (300, 'K'),
+        Tmax = (1500, 'K'),
+    ),
+    rank = 2,
+)
+
+entry(
+    index = 2000,
+    label = "Y_rad;N_R_birad",
+    kinetics = ArrheniusEP(
+        A = (1e+13, 'cm^3/(mol*s)'),
+        n = 0,
+        alpha = 0,
+        E0 = (0, 'kcal/mol'),
+        Tmin = (300, 'K'),
+        Tmax = (1500, 'K'),
+    ),
+    rank = 2,
+)
+

…cs/families/Oa_R_Recombination/groups.py …families/Birad_R_Recombination/groups.pyinput/kinetics/families/Oa_R_Recombination/groups.py renamed to input/kinetics/families/Birad_R_Recombination/groups.py input/kinetics/families/Oa_R_Recombination/groups.py renamed to input/kinetics/families/Birad_R_Recombination/groups.py

@@ -0,0 +1,57 @@
+O
+multiplicity 3
+1 *2 O u2 p2 c0
+
+S
+multiplicity 3
+1 *2 S u2 p2 c0
+
+NH
+multiplicity 3
+1 *2 N u2 p1 c0 {2,S}
+2    H u0 p0 c0 {1,S}
+
+HSO
+multiplicity 2
+1    S u0 p2 c0 {2,S} {3,S}
+2 *1 O u1 p2 c0 {1,S}
+3    H u0 p0 c0 {1,S}
+
+HSOO
+multiplicity 2
+1    S u0 p2 c0 {2,S} {4,S}
+2 *1 O u0 p2 c0 {1,S} {3,S}
+3 *2 O u1 p2 c0 {2,S}
+4    H u0 p0 c0 {1,S}
+
+NO
+multiplicity 2
+1 *1 N u1 p1 c0 {2,D}
+2    O u0 p2 c0 {1,D}
+
+NO2_p
+multiplicity 2
+1 *1 N u0 p1 c0 {2,D} {3,S}
+2    O u0 p2 c0 {1,D}
+3 *2 O u1 p2 c0 {1,S}
+
+NO2_r
+multiplicity 2
+1    O u0 p2 c0 {2,D}
+2 *1 N u1 p0 c+1 {1,D} {3,S}
+3    O u0 p3 c-1 {2,S}
+
+HNNO2
+multiplicity 2
+1 *2 N u1 p1 c0 {2,S} {3,S}
+2    H u0 p0 c0 {1,S}
+3 *1 N u0 p0 c+1 {1,S} {4,S} {5,D}
+4    O u0 p3 c-1 {3,S}
+5    O u0 p2 c0 {3,D}
+
+SNO
+multiplicity 2
+1 *2 S u1 p2 c0 {2,S}
+2 *1 N u0 p1 c0 {1,S} {3,D}
+3    O u0 p2 c0 {2,D}
+