def METHOD_NAME(mechanism, reaction_info, qssa_species):
    """Identify the reactions that should be removed."""
    reactions_to_remove = {}
    for orig_idx, _ in reaction_info.idxmap.items():
        reaction = mechanism.reaction(orig_idx)
        # Check for quadratic relations
        qssa_species_involved = cu.intersection(
            list(reaction.reactants.keys()), qssa_species
        )
        sum_coeff = sum(
            v for k, v in reaction.reactants.items() if k in qssa_species_involved
        )
        if sum_coeff > 1:
            if orig_idx not in reactions_to_remove:
                reactions_to_remove[orig_idx] = []
            reactions_to_remove[orig_idx].append("f")
        if reaction.reversible:
            qssa_species_involved = cu.intersection(
                list(reaction.products.keys()), qssa_species
            )
            sum_coeff = sum(
                v for k, v in reaction.products.items() if k in qssa_species_involved
            )
            if sum_coeff > 1:
                if orig_idx not in reactions_to_remove:
                    reactions_to_remove[orig_idx] = []
                reactions_to_remove[orig_idx].append("r")
    return reactions_to_remove