for reactant in reactants:
        obj_reaction = obj_reactants[reactant]
        conv_model.objective_terms = [conv_opt.LinearTerm(obj_reaction, 1.),]
        result = conv_model.solve()
        if result.value == 0.0: