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

Could not determine reaction family for reaction C3H6-2 + O <=> C3H6O #738

Open
kfir4444 opened this issue Apr 15, 2024 · 1 comment
Open

Comments

@kfir4444
Copy link
Collaborator

Describe the bug
A reaction family could not be found for the reaction C3H6-2 + O <=> C3H6O.

How to reproduce

from arc.rmgdb import make_rmg_database_object, load_families_only
from arc.reaction import ARCReaction
from arc.species import ARCSpecies

db = make_rmg_database_object()
load_families_only(db, "all")

rxn = ARCReaction(r_species=[ARCSpecies(label="C3H6", smiles = "CC=C"),
                             ARCSpecies(label="O", smiles="[O]", multiplicity=1)],
                  p_species=[ARCSpecies(label="C3H6O", smiles = "CC1OC1")])

rxn.determine_family(db)
print(rxn.family) # result is None with errors.

This remains true even if multiplicity=3 in the [O] species. This reaction originates from the training depository of 1+2_Cycloaddition reaction family originally.
Expected behavior

print(rxn.family) -> <ReactionFamily "1+2_Cycloaddition">
@kfir4444
Copy link
Collaborator Author

kfir4444 commented Apr 15, 2024

Also true for:

'C=C(C)C + O <=> CC1(C)CO1'
'CC=CC + O <=> CC1OC1C'
'CC(C)=C(C)C + O <=> CC1(C)OC1(C)C'

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

No branches or pull requests

1 participant