The Physics of a Gas Mixture: Analyzing the Air We Breathe

Today, I’m using modern code to solve fundamental chemistry problems, and I hope this series of posts will help students and enthusiasts bridge the gap between abstract equations and practical, logical solutions. Problem: We know that the atmosphere isn't just a mix of 78% Nitrogen and 21% Oxygen. It is a complex mixture of at least 12 significant gases. Our goal is to use this Volume Percentage to calculate the Average Molecular Weight (\bar{M}) and the Density (\rho) of dry air at Standard Temperature and Pressure (STP). To do this, we rely on a fundamental constant from physical chemistry: the molar volume of an ideal gas at STP, which is 22,400 cc/mol (cm^3/mol). The Logical Steps (The "Math" before the Code) Instead of just asking an AI for the answer, we will build a logical script that mirrors a scientific calculation: * Define the knowns: We create a dictionary to store the Volume Percentage and Molecular Weight of each gas (Nitrogen, Oxygen, CO_2, etc.). * Calculate the masses: We use the 22,400 cc/mol constant to determine how many grams of each gas would exist in a 100 cc sample. * Sum it up: We find the total mass of the 100 cc sample. * Derive the answers: We find the average molecular weight from this total mass and the density by dividing the mass by the volume. The Python Solution

Python: The Physics of a Gas Mixture
Author: [Ranjit Singh]
Molar Volume at STP (Standard Temp and Pressure) in cc/mol
MOLAR_VOLUME = 22400
{ "Name": (Volume %, Molecular Weight) }
COMPOSITION = {
"Nitrogen, N2": (78.08, 28.013),
"Oxygen, O2": (20.95, 31.998),
"Argon, Ar": (0.93, 39.948),
"Carbon dioxide, CO2": (0.031, 44.009),
"Hydrogen, H2": (5.9e-03, 2.016),
"Neon, Ne": (1.8e-03, 20.180),
"Helium, He": (5.2e-04, 4.003),
"Methane, CH4": (2.0e-04, 16.043),
"Krypton, Kr": (1.1e-04, 83.798),
"Nitric oxide, NO": (5.0e-05, 30.006),
"Xenon, Xe": (8.7e-06, 131.29),
"Ozone, O3": (7.0e-06, 47.998)
}
sum_vol = 0
sum_mass = 0
mass_dict = {}
Header for the detailed output
print(f"{'Component':<22} {'Vol %':<10} {'Mass (g/100cc)'}")
print("-" * 52)
for gas, (vol, mlwt) in COMPOSITION.items():
sum_vol += vol
# Calculate mass of each component in 100cc of air
# Logic: (Molar_Mass / Molar_Volume) * Vol %
mass = (mlwt / MOLAR_VOLUME) * vol
sum_mass += mass
mass_dict[gas] = mass
print(f"{gas:<22} {vol:<10} {mass:.8f}")

Average Molecular Weight calculation
Avg_MLWT = (Total Mass * Molar_Volume) / Volume Sample
avg_mlwt = (sum_mass * MOLAR_VOLUME) / 100
Density (mass / volume)
density = sum_mass / sum_vol
print("\n" + "="*45)
print(f"Total Volume Checked:  {sum_vol:.2f}%")
print(f"Avg Mol Weight of Air: {avg_mlwt:.2f} g/mol")
print(f"Density of Air:        {density:.5e} g/cc")
print("="*45 + "\n")
Percentage by Mass calculation
print(f"{'Component':<22} {'% By Mass'}")
print("*" * 35)
for g, m in mass_dict.items():
m_percent = (m / sum_mass) * 100
print(f"{g:<22} {m_percent:>10.5f}%")

For any students reading, the constant 22,400 is derived from the Ideal Gas Law (PV = nRT) for a single mole (n=1) at 273.15K and 1 atm. In real-world exploration at ONGC, we knew that conditions were rarely "standard," but this script provides the necessary logical framework. We could easily modify this code to calculate the density at a different pressure (like at the bottom of a well) or even account for humidity. If you found this useful, my next post will look at solving the Ideal Gas Law using Python and handling common unit conversions automatically.