for product in products:
        obj_reaction = obj_products[product]
        conv_model.objective_terms = [conv_opt.LinearTerm(obj_reaction, 1.),]
        result = conv_model.solve()
        if result.value == 0.0: