Volatile Liquid Molar Mass theory

The Dumas method is rooted in the fundamental gas equation of state PV = nRT, which establishes that at a specified thermodynamic state defined by pressure P and absolute temperature T, the mass density of a vapor ρ_vapor = m / V is directly proportional to its molecular weight M: ρ_vapor = (P · M) / (R · T). Jean-Baptiste Dumas recognized that by maintaining a closed vessel at constant temperature in a boiling water bath (T = 373.15 K) and open to the atmosphere through a pinhole (P = P_atm), the condition inside the vessel becomes invariant. When an excess of volatile liquid boils inside the vessel, the volume of vapor produced is typically 200 to 500 times greater than the original liquid volume. This continuous outward flow of expanding vapor mechanically purges all atmospheric nitrogen and oxygen from the container. At the exact instant the last liquid droplet evaporates, the container is filled exclusively with pure vapor at T_bath and P_atm. The number of moles of vapor trapped inside is given by n = (P_atm · V_flask) / (R · T_bath). When the vessel is cooled, this exact quantity of vapor condenses into liquid, whose mass m_vapor can be weighed with fourth-decimal precision on an analytical balance. Rearranging the ideal gas law gives the classical Dumas expression: M = (m_vapor · R · T_bath) / (P_atm · V_flask). Because R is a universal constant, determining M requires only three macroscopic measurements: the mass of condensed vapor m_vapor, the internal flask volume V_flask, and the thermodynamic state variables (P_atm, T_bath).

The physical mechanism that enables the Dumas bulb to trap pure vapor without contamination by ambient air relies upon continuous convective displacement and fluid flow through the micro-pinhole. Consider a Dumas bulb containing 5.0 mL of liquid cyclohexane (density ρ_liq = 0.778 g/mL at 20 °C, net mass ~3.89 g). At 100 °C and 1.000 atm, the ideal vapor volume generated by 3.89 g of cyclohexane is V_vapor = n · R · T / P = (3.89 / 84.16) · 0.08206 · 373.15 = 1.414 L (1414 mL). For a 275 mL Dumas flask, the boiling liquid generates more than five full container volumes of vapor. As cyclohexane boils, the rapid generation of vapor increases internal pressure slightly above atmospheric pressure (P_inside = P_atm + ΔP), driving a vigorous, unidirectional jet of vapor through the 0.5 mm pinhole. This convective mass flux sweeps out all initial air molecules via plug flow and turbulent mixing. The Peclet number for mass transfer Pe = (u · L) / D_air-vapor is very large, meaning convective outflow completely overwhelms inward molecular diffusion of atmospheric air. However, once the last droplet of liquid vaporizes, the generation of vapor ceases abruptly, ΔP drops to zero, and the system transitions from convective outflow to diffusive exchange. If the flask is left in the bath too long, ambient air will slowly back-diffuse through the pinhole. Conversely, if the flask is removed too early, residual unvaporized liquid droplets (possessing 300 times the density of vapor) will remain. Precise timing—removing the flask approximately 60 to 90 seconds after visible liquid disappears—is critical for achieving accurate vapor density.

A subtle but significant systematic error in classical Dumas gravimetry arises from air buoyancy. When the clean, dry flask capped with aluminum foil is weighed initially to record the tare mass m_tare, it is not empty; it is completely filled with ambient laboratory air at room temperature T_room and barometric pressure P_atm. The analytical balance reading reflects the mass of the glass and foil minus the buoyant force of the surrounding air, plus the mass of air trapped inside: m_tare,app = m_glass + m_foil + m_air(inside) − B_outside. Because the exterior volume of the flask displaces an equal volume of ambient air, the buoyant force B_outside exactly cancels the mass of air inside the flask if the walls have negligible thickness. However, when the flask is removed from the bath, cooled, and weighed as m_gross, the internal space contains condensed liquid plus room air that rushed in through the pinhole to replace the condensed vapor. If the condensed liquid mass is calculated simply as m_gross − m_tare, an error occurs because the initial tare included a full flask volume of air (m_air = ρ_air · V_flask), whereas the gross weighing contains liquid plus air occupying V_flask − V_liquid. Because air density at 20 °C and 1 atm is ρ_air ≈ 0.00120 g/mL, for a 275 mL flask, the mass of air displaced is approximately 275 · 0.00120 = 0.330 g. If uncorrected, this systematic bias shifts the calculated vapor mass by over 40%! The rigorous buoyancy-corrected mass of condensed vapor is m_vapor,true = (m_gross − m_tare) + (ρ_air · V_flask). Multi-point regression across varying volumes provides an elegant alternative: by plotting m_vapor versus V_flask, the air buoyancy term scales linearly with volume and is absorbed directly into the slope, providing an unbiased estimator.

The ideal gas law assumes that molecules exert no intermolecular forces and occupy negligible physical volume. In reality, cyclohexane (C₆H₁₂) vapor at 100 °C (373.15 K) and 1 atm is only 19 °C above its normal boiling point (80.7 °C), operating in a thermodynamic regime where intermolecular attractive London dispersion forces are substantial. In the van der Waals equation of state, [P + a·(n/V)²] · [V - n·b] = n·R·T, the parameter a for cyclohexane is large: a = 23.3 L²·bar/mol² (compared to 0.25 for H₂ and 3.6 for CO₂), reflecting the high polarizability of the 84 g/mol hydrocarbon ring. The excluded volume parameter is b = 0.144 L/mol. The compressibility factor Z = PV / (nRT) evaluates to Z ≈ 0.985 at 1 atm and 373 K. Because Z < 1, attractive intermolecular forces dominate, causing the actual pressure exerted by the vapor to be lower than an ideal gas, or equivalently, causing a given volume at 1 atm to contain more molecules (higher mass density) than predicted by the ideal gas law. Specifically, the real vapor density is ρ_real = ρ_ideal / Z = ρ_ideal / 0.985 ≈ 1.015 · ρ_ideal. If an analyst uses the uncorrected ideal gas formula M = (ρ · R · T) / P, the experimental molar mass will be systematically overestimated by approximately 1.5% (yielding ~85.4 g/mol instead of 84.16 g/mol). Applying the Berthelot equation of state or van der Waals correction M_true = M_ideal · Z eliminates this non-ideality bias, restoring accuracy to within 0.1%.

In Dumas experiments, two experimental errors relating to phase equilibrium produce catastrophic positive errors in calculated molar mass: incomplete vaporization and premature condensation. (1) Incomplete Vaporization: Liquid cyclohexane has a density of 0.778 g/mL, whereas its vapor at 100 °C has a density of 0.00275 g/mL. The liquid is therefore approximately 283 times more dense than the vapor. If an analyst mistakenly removes the flask from the boiling water bath while even a single tiny droplet of unvaporized liquid remains—for example, a droplet with volume of 0.05 mL (mass ~0.039 g)—that droplet alone contributes as much mass as 14.2 mL of pure vapor! In a 250 mL flask, this represents a +5.7% positive error in m_vapor, falsely inflating the calculated molar mass from 84.16 g/mol to ~89.0 g/mol. (2) Premature Condensation (Refluxing): If the flask is not submerged deeply enough in the boiling water bath, the upper neck and shoulder of the flask will be exposed to ambient air and remain at temperatures below 80.7 °C. Cyclohexane vapor rising into this cooler neck condenses on the glass walls and drips back down (refluxing). Furthermore, droplets clinging to the neck at the moment of removal are retained. Ensuring that the boiling water level reaches within 1 cm of the aluminum foil cap, using a tall beaker, and providing bright backlighting to confirm complete transparency inside the bulb are critical procedural safeguards.

A standard laboratory Erlenmeyer flask has nominal volumetric graduations with manufacturing tolerances of ±5% to ±10%, which is completely unacceptable for high-precision physical chemistry metrology. In the Dumas method, the internal volume V_flask must be determined individually for every flask via gravimetric water pycnometry. Following the vapor density determination, the foil cap is removed and the flask is filled to overflowing with degassed, deionized water. A flat, polished glass cover slip is slid horizontally across the ground rim of the flask to shear off excess liquid without trapping air bubbles. The exterior glass is dried thoroughly with a microfiber towel, and the gross mass m_water+flask is recorded on an analytical balance. Subtracting the clean, dry empty flask mass m_empty isolates the net mass of water filling the vessel: m_water = m_water+flask − m_empty. The internal volume is then calculated from the thermodynamic density of liquid water at the measured water temperature T_water: V_flask = m_water / ρ_water(T_water). At 20.0 °C, ρ_water = 0.998203 g/mL; at 22.0 °C, ρ_water = 0.997770 g/mL. Because analytical balances measure mass to ±0.0001 g, gravimetric pycnometry determines flask volume with an extraordinary accuracy of ±0.02 mL (relative uncertainty < 0.01% in a 250 mL bulb), completely eliminating volumetric uncertainty as a limiting factor in the Dumas method.

A rigorous sensitivity analysis of the Dumas equation M = (m_vapor · R · T_bath) / (P_atm · V_flask) evaluates how experimental uncertainties in individual measured variables propagate into the final molar mass. Taking the natural logarithm and total differential yields: dM / M = dm / m + dT / T − dP / P − dV / V. In accordance with the GUM (Guide to the Expression of Uncertainty in Measurement), the combined relative standard uncertainty is (u_c(M) / M)² = (u(m) / m)² + (u(T) / T)² + (u(P) / P)² + (u(V) / V)². Evaluating typical laboratory magnitudes for a 275 mL bulb containing ~0.7559 g of cyclohexane vapor: (1) Mass Uncertainty: Four-decimal balance weighing error of ±0.0003 g on ~0.75 g gives u(m)/m = 0.0003 / 0.7559 = 0.040%. (2) Temperature Uncertainty: Immersion thermometer calibrated to ±0.1 °C in 373.15 K gives u(T)/T = 0.1 / 373.15 = 0.027%. (3) Pressure Uncertainty: Mercury barometer vernier reading to ±0.1 mmHg in 760.0 mmHg gives u(P)/P = 0.1 / 760.0 = 0.013%. (4) Volume Uncertainty: Pycnometric water calibration to ±0.05 mL in 275 mL gives u(V)/V = 0.05 / 275.0 = 0.018%. Combining these in quadrature yields a theoretical instrumental uncertainty of u_c(M) / M = sqrt(0.040² + 0.027² + 0.013² + 0.018²) ≈ 0.053% (corresponding to ±0.045 g/mol). In practice, non-instrumental systematic effects—primarily residual water droplets on the exterior foil crimp and vapor non-ideality—dominate. Multi-point regression across six calibrated flasks eliminates random transfer biases and verifies stoichiometric linearity with R² = 1.00000.