Skip to content

Commit

Permalink
Revert "remove non-spin conserving reactions in explorer"
Browse files Browse the repository at this point in the history
This reverts commit 1ac1dd5.
  • Loading branch information
donerancl committed Jun 27, 2024
1 parent c7b1679 commit 9fdc594
Showing 1 changed file with 0 additions and 10 deletions.
10 changes: 0 additions & 10 deletions arkane/explorer.py
Original file line number Diff line number Diff line change
Expand Up @@ -264,16 +264,6 @@ def execute(self, output_file, plot, file_format='pdf', print_summary=True, spec
for r in rxn.reactants + rxn.products:
if forbidden_structures.is_molecule_forbidden(r.molecule[0]):
rm_rxns.append(rxn)

print("In ExplorerJob.execute\n",flush=True)
print(f"Num reactants {len(rxn.reactants)}, num products {len(rxn.products)}\n",flush=True)
if len(rxn.reactants) > 1 and len(rxn.products) > 1:
num_unpaired_electrons_reactants = sum([r.multiplicity-1 for r in rxn.reactants])
num_unpaired_electrons_products = sum([r.multiplicity-1 for r in rxn.products])
print(f"Unpaired electrons {num_unpaired_electrons_reactants} => {num_unpaired_electrons_products}\n",flush=True)
if num_unpaired_electrons_products!=num_unpaired_electrons_reactants:
print("Removing reaction\n",flush=True)
rm_rxns.append(rxn)

for rxn in rm_rxns:
logging.info('Removing forbidden reaction: {0}'.format(rxn))
Expand Down

0 comments on commit 9fdc594

Please sign in to comment.