for reactant in reactants:            
        obj_reaction = conv_opt.Variable(name=reactant + '_test_reaction', 
                                type=conv_opt.VariableType.continuous,
                                lower_bound=0.0)
        obj_reactants[reactant] = obj_reaction