Skip to content

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

@kfir4444

Description

@kfir4444

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">

Metadata

Metadata

Assignees

No one assigned

    Labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions