Skip to content

Uncertainty Part 1.5: Adding in quadrature#2933

Draft
sevyharris wants to merge 8 commits intoReactionMechanismGenerator:mainfrom
sevyharris:uncertainty_tool_quadrature
Draft

Uncertainty Part 1.5: Adding in quadrature#2933
sevyharris wants to merge 8 commits intoReactionMechanismGenerator:mainfrom
sevyharris:uncertainty_tool_quadrature

Conversation

@sevyharris
Copy link
Copy Markdown
Contributor

Dependencies

Motivation or Problem

This overhauls the uncertainty tool to add uncertainties in quadrature (variances add instead of standard deviations)

Description of Changes (each change corresponds to a commit)

  • Make uncertainty contributions add in quadrature.
  • Set dG_group=0.586 kcal/mol as determined by the fitting in Demo notebook (fits dG_group to best match of GAV error for a set of 564 species with high quality DFT to compare to)
  • TODO: update math explanation in local uncertainty notebook to reflect the changes

Testing

  • Demo notebook comparing quadrature uncertainty estimation to previous method, and why dG_group was set to 0.7159 kcal/mol
  • Here's a little explainer pdf I made for myself why uncertainty should be added in quadrature
    1. Theory of Adding in Quadrature.pdf

Up until now the uncertainty tool has assumed the kinetics trees are the original hand-built type. This enables the uncertainty tool to grab the calculated standard deviation and number of training reactions from the specific autogen tree node used, instead of incorrectly treating it as an exact match of a single rate rule with a std dev of 0.5
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

Successfully merging this pull request may close these issues.

1 participant