Figuring out the limiting reactant in a chemical reaction integrates many of the basic chemistry concepts including: unit conversions from moles to mass and vice versa; the meaning of chemical formulas; and understanding the stoichiometry of chemical reactions. So, since we’ll need a number of these, I wrote a python program to help me design the questions (and figure out the answers).
Program examples come zipped because they require the program file and the elements_database.py library:
Baking Powder and Vinegar (Common Molecules)
Limiting_component-Common.py: This has the baking powder and vinegar reaction limited by 5 g of baking soda. It’s nice because it uses a few pre-defined “common molecules” (which are defined in the elements_database.py library.
You enter the reactants and products and the program checks if the reaction is balanced, then calculates the moles and masses based on the limiting component, and finally double checks to make sure the reaction is mass balanced.
Limiting_component-Common.py
from elements_database import *
import sys
print "LIMITING REACTANT PROGRAM"
print
print " Determines the needed mass and moles of reactants and products if reaction is limited by one of the components"
c = common_molecules()
'''Create Reaction'''
rxn = reaction()
# Add reactants (and products)
# Use rxn.add_reactant(molecule[, stoichiometry])
# molecule: from molecule class in elements_database
# stoichiometry: integer number
rxn.add_reactant(c.baking_soda,1)
rxn.add_reactant(c.hydrochloric_acid,1)
rxn.add_product(c.carbon_dioxide)
rxn.add_product(c.water, 1)
rxn.add_product(c.salt)
'''Print out the reaction'''
print
print "Chemical Formula"
print " " + rxn.print_reaction()
print
'''Check if reaction is balanced'''
balanced = rxn.check_for_balance(printout=True)
'''Calculate limits of reaction'''
if balanced:
rxn.limited_by_mass(c.baking_soda, 5, printout=True)
Outputs results in the Results table (using scientific notation):
LIMITING REACTANT PROGRAM
Determines the needed mass and moles of reactants and products if reaction is limited by one of the components
Chemical Formula
NaHCO3 + HCl --> CO2 + H2O + NaCl
Check for balance
---------------------------------------------
| Element | Reactants | Products | Difference |
---------------------------------------------
| Na | 1 | -1 | 0 |
| H | 2 | -2 | 0 |
| C | 1 | -1 | 0 |
| O | 3 | -3 | 0 |
| Cl | 1 | -1 | 0 |
---------------------------------------------
Balance is: True
Given: Limiting component is 5 g of NaHCO3.
Molar mass = 84.00676928
Moles of NaHCO3 = 0.0595190130849
Results
------------------------------------------------------------------
| Molecule | Stoich.*| Molar Mass (g) | Moles req. | Mass (g) |
------------------------------------------------------------------
|NaHCO3 | 1 | 84.0068 | 5.952e-02 | 5.000e+00 |
|HCl | 1 | 36.4602 | 5.952e-02 | 2.170e+00 |
|CO2 | -1 | 44.0096 | 5.952e-02 | 2.619e+00 |
|H2O | -1 | 18.0156 | 5.952e-02 | 1.072e+00 |
|NaCl | -1 | 58.4418 | 5.952e-02 | 3.478e+00 |
------------------------------------------------------------------
* negative stoichiometry means the component is a product
Final Check: Confirm Mass balance:
Reactants: 7.1701 g
Products: -7.1701 g
--------------------------
= 0.0000 g
--------------------------
General Example
If not using the common molecules database, you need to define the components in the reaction as molecules yourself. This example reacts magnesium sulfate and sodium hydroxide, and limits the reaction with 20 g of magnesium sulfate.
Limiting_component-General.zip
The main file is:
Chem_Exam-limiting.py
from elements_database import *
import sys
print "LIMITING REACTANT PROGRAM"
print
print " Determines the needed mass and moles of reactants and products if reaction is limited by one of the components"
# create reaction
rxn2 = reaction()
# Add reactants (and products)
# Use rxn.add_reactant(molecule[, stoichiometry])
# molecule: from molecule class in elements_database
# stoichiometry: integer number
rxn2.add_reactant(molecule("Mg:1,S:1,O:4"))
rxn2.add_reactant(molecule("Na:1,O:1,H:1"), 2)
rxn2.add_product(molecule("Mg:1,O:2,H:2"))
rxn2.add_product(molecule("Na:2,S:1,O:4"))
'''Print out the reaction'''
print
print "Chemical Formula"
print " " + rxn2.print_reaction()
print
'''Check if reaction is balanced'''
balanced = rxn2.check_for_balance(printout=True)
'''Calculate limits of reaction'''
if balanced:
rxn2.limited_by_mass(molecule("Mg:1,S:1,O:4"), 20, True)
# print out masses
print "Masses Involved in Reaction"
print " Reactants:"
for i in rxn2.reactants:
#print "rxn", i
print " {m}: {g} g".format(m=i.molecule.print_formula().ljust(10), g=i.mass)
print " Products:"
for i in rxn2.products:
#print "rxn", i
print " {m}: {g} g".format(m=i.molecule.print_formula().ljust(10), g=-i.mass)
This program is slightly different from the common molecules example in that, at the end, it prints out masses calculated in a more easily readable format in addition to the other data.
Masses Involved in Reaction
Reactants:
MgSO4 : 20.0 g
NaOH : 13.2919931774 g
Products:
MgO2H2 : 9.69066512026 g
Na2SO4 : 23.6013280571 g
When I have some time I’ll convert this to JavaScript like the molecular mass example.