Article — Langmuir Isotherm Calculator (θ, qₑ, R_L)
Langmuir Isotherm Calculator: θ, qₑ, and Separation Factor
The Langmuir isotherm describes monolayer adsorption on a uniform surface. Its core equation θ = KP/(1+KP) (gases) or qₑ = q_max·K_L·Cₑ/(1+K_L·Cₑ) (solutions) captures how surface coverage grows from zero at empty conditions to saturation at high concentration. Irving Langmuir won the 1932 Nobel Prize in Chemistry for the underlying surface kinetics.
The calculator handles three tasks: equilibrium loading qₑ from concentration, fractional coverage θ from KP, and the separation factor R_L that diagnoses whether the isotherm is favorable.
What is the Langmuir isotherm?
The Langmuir isotherm is a two-parameter mathematical model for adsorption equilibrium. At low concentration the model predicts linear behavior (Henry's law region); at high concentration it saturates at q_max. The crossover happens around K_L·Cₑ = 1, where θ = 0.5.
The model is derived kinetically. At equilibrium, the rate of adsorption (proportional to vacant sites and concentration) equals the rate of desorption (proportional to occupied sites). Setting the two equal and solving gives the famous hyperbolic form.
The Langmuir isotherm formula
θ = KP / (1 + KP)q_e = q_max · K_L · C_e / (1 + K_L · C_e)R_L = 1 / (1 + K_L · C_0)C_e/q_e = 1/(K_L·q_max) + C_e/q_maxFor activated carbon with q_max = 200 mg/g and K_L = 0.05 L/mg at Cₑ = 20 mg/L: qₑ = 200 · 0.05 · 20 / (1 + 0.05 · 20) = 200/2 = 100 mg/g. Half-saturation. The separation factor with C₀ = 50 mg/L gives R_L = 1/(1 + 2.5) = 0.286 — a favorable isotherm.
Langmuir isotherm assumptions
The model rests on five idealizations:
- Monolayer only — adsorbate forms a single molecular layer; no stacking
- Uniform surface — every site has the same binding energy
- No interactions — adsorbed molecules do not influence neighbors
- Reversibility — adsorption and desorption are both possible
- Localized binding — molecules occupy discrete sites, not a mobile film
Real surfaces violate at least one of these. The Langmuir fit often works anyway because two parameters give enough flexibility to match a moderate range of data. Outside the fitted range, predictions can be unreliable.
Irving Langmuir worked at General Electric for 41 years and made fundamental contributions to gas-discharge physics, plasma chemistry, surface chemistry, and meteorology (cloud seeding). His 1918 derivation of the adsorption isotherm was originally meant to explain the operation of tungsten filaments in vacuum tubes.
Langmuir vs Freundlich vs BET
Three isotherms dominate textbook coverage of adsorption:
Langmuir — saturating monolayer on uniform sites. Best for chemisorption and low-coverage physisorption.
Freundlich — empirical q = K_F · C^(1/n), no saturation. Handles heterogeneous surfaces and multilayers approximately. Common for dye and heavy metal adsorption on natural sorbents.
BET — Brunauer-Emmett-Teller extends Langmuir to multilayer gas adsorption. The standard for measuring solid surface area by N₂ adsorption at 77 K (the "BET surface area" reported for catalysts and porous materials).
Practical Langmuir isotherm applications
The Langmuir framework appears in nearly every field that deals with surface chemistry:
- Water treatment — activated carbon removal of dyes, phenols, pesticides
- Heterogeneous catalysis — modeling site occupancy for catalyst design
- Gas separation — zeolite and MOF capacity for CO₂, CH₄, H₂
- Biotechnology — receptor-ligand binding (mathematically identical form)
- Pharmacology — drug binding curves, dose-response
- Soil science — phosphate and trace metal retention
- Sensor design — surface plasmon resonance, biosensor calibration
A linearized fit can show R² > 0.99 even when the Langmuir model is fundamentally wrong. The linearization transforms errors unevenly — small absolute errors at high concentration become large reciprocal errors. Always check residuals against the non-linear fit, or run a non-linear least squares regression instead.
Fitting Langmuir isotherm data
Three common linearizations:
Hanes-Woolf: Cₑ/qₑ = 1/(K_L·q_max) + Cₑ/q_max. Plot Cₑ/qₑ vs Cₑ. Slope = 1/q_max, intercept = 1/(K_L·q_max). Often the cleanest linearization.
Lineweaver-Burk: 1/qₑ = 1/q_max + 1/(K_L·q_max·Cₑ). Plot 1/qₑ vs 1/Cₑ. Familiar from enzyme kinetics; weights low concentrations heavily.
Eadie-Hofstee: qₑ = q_max − qₑ/(K_L·Cₑ). Plot qₑ vs qₑ/Cₑ. Errors appear on both axes; not recommended for primary fitting.
Modern practice is non-linear regression on the raw qₑ vs Cₑ data. Most spreadsheet tools (Excel solver, Python scipy.optimize, GraphPad Prism) handle the two-parameter fit in seconds.
Design adsorption experiments around the half-saturation concentration C_50 = 1/K_L. Collect data points spaced logarithmically from 0.1·C_50 to 10·C_50. Below that range you only see linear behavior; above it you only see saturation. The mid-range is where Langmuir and competitor models diverge most clearly.
Common Langmuir isotherm mistakes
The frequent errors:
- Extrapolating beyond data — Langmuir parameters fit to 0 to 50 mg/L may fail at 500 mg/L
- Wrong units — K_L in L/mg vs L/mol changes the numerical value by molar mass
- Ignoring equilibrium time — sampling before equilibrium gives apparent capacities below true q_max
- Mixing pH — both K_L and q_max usually shift with pH; report and control it
- Forgetting temperature — physisorption K_L drops at higher temperature
- Confusing R_L with R² — R_L is the separation factor (0 to 1), R² is the regression goodness of fit
Irving Langmuir published the monolayer model in 1918, working on tungsten filaments at General Electric. Twenty years later, in 1938, Stephen Brunauer, Paul Emmett, and Edward Teller extended the analysis to multilayer adsorption, giving the BET isotherm. The BET surface-area determination — using nitrogen adsorption at 77 K — remains the global standard for characterizing porous solids. Catalysts, activated carbons, and pharmaceutical excipients are all spec'd with BET areas, and every BET measurement starts from Langmuir's 1918 mass-balance argument.
For day-to-day adsorption problems, Langmuir is the right starting point. If the fit is poor, move to Freundlich; if the isotherm shows multilayer behavior, move to BET; if cooperative effects appear, move to Sips or Toth. But always start with Langmuir — two parameters and a transparent physical picture.
The model also generalizes beautifully outside chemistry. Receptor-ligand binding curves in pharmacology are mathematically identical to the Langmuir form, with K_d playing the role of 1/K_L and B_max playing the role of q_max. Substrate binding in enzyme kinetics (Michaelis-Menten) shares the same hyperbolic shape. Once you recognize the Langmuir form, you see it everywhere a finite number of binding sites are filling up reversibly.