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

Thermo inconsistencies in training reactions #538

Open
davidfarinajr opened this issue Oct 21, 2021 · 0 comments
Open

Thermo inconsistencies in training reactions #538

davidfarinajr opened this issue Oct 21, 2021 · 0 comments

Comments

@davidfarinajr
Copy link
Contributor

In H abstraction training, there is at least on example where a reaction is specified in both the forward and reverse directions, but the kinetics are inconsistent with the thermo. This can cause a lot of trouble if using the the thermo to generate reverse rate constants when training the tree since the reverse kinetics of one direction will not match the forward kinetics of the other direction. This is happening for OH + N <=> HN + O, which results in a node with very high uncertainty.

entry(
    index = 625,
    label = "HO-4 + N <=> HN + O-2",
    degeneracy = 1.0,
    kinetics = Arrhenius(A=(6.4e+12,'cm^3/(mol*s)'), n=0.1, Ea=(88.9518,'kJ/mol'), T0=(1,'K'), Tmin=(300,'K'), Tmax=(2000,'K')),
    rank = 2,
    shortDesc = """Added by Beat Buesser""",
    longDesc = 
"""
Added by Beat Buesser, value for reaction: N + OH = NH + O (B&D #26b) in 'Gas-Phase Combustion Chemistry' (ISBN: 978-1-4612-7088-1), chapter 2, 'Combustion Chemistry of Nitrogen', Anthony M. Dean, Joseph W. Bozzelli",

Converted to training reaction from rate rule: OH_rad_H;N_atom_quartet
""",
)
entry(
    index = 629,
    label = "HN-4 + O <=> HO-2 + N-2",
    degeneracy = 1.0,
    kinetics = Arrhenius(A=(1.7e+08,'cm^3/(mol*s)'), n=1.5, Ea=(217.878,'kJ/mol'), T0=(1,'K'), Tmin=(300,'K'), Tmax=(2000,'K')),
    rank = 2,
    shortDesc = """Added by Beat Buesser""",
    longDesc = 
"""
Added by Beat Buesser, value for reaction: NH + O = N + OH (B&D #27e2) in 'Gas-Phase Combustion Chemistry' (ISBN: 978-1-4612-7088-1), chapter 2, 'Combustion Chemistry of Nitrogen', Anthony M. Dean, Joseph W. Bozzelli",

Converted to training reaction from rate rule: NH_triplet_H;O_atom_triplet
""",
)

note: the training indices are different then above because I am using a different branch
Screen Shot 2021-10-21 at 12 19 23 PM
Screen Shot 2021-10-21 at 12 23 44 PM

tree_cross_val.csv
tree_cross_val.pdf

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant