Bridging the gap in statistical thermodynamics: A rigorous justification of the cartesian quantum harmonic oscillator model for diatomic molecules

In advanced undergraduate and graduate mechanical engineering courses on statistical thermodynamics, the one-dimensional Cartesian quantum harmonic oscillator (QHO) is a standard model used to calculate energy levels used in the vibrational partition functions of diatomic gases. However, when inquisitive students compare this model to the full radial Schrödinger equation, which allows for motion in three dimensions, they encounter a confusing discrepancy: the radial equation contains a first-derivative term, <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline" overflow="scroll"> <mml:mrow> <mml:mfrac> <mml:mn>2</mml:mn> <mml:mi>r</mml:mi> </mml:mfrac> </mml:mrow> <mml:mrow> <mml:mfrac> <mml:mrow> <mml:mi>d</mml:mi> <mml:mi>R</mml:mi> </mml:mrow> <mml:mrow> <mml:mi>d</mml:mi> <mml:mi>r</mml:mi> </mml:mrow> </mml:mfrac> </mml:mrow> </mml:math> , which is absent in the Cartesian approximation. This leads to the pedagogical question: ”Where did the first derivative go?” To the best knowledge of the author, the answer to this question is not available in the literature. This paper provides engineering educators with a derivation to answer this question by using matched asymptotic expansions. We demonstrate that the Cartesian QHO equation corresponds to the leading-order inner-region problem of the full radial equation in the case where angular momentum effects are negligible. This approach provides a mathematical foundation for a core engineering thermodynamics concept and introduces students to perturbation methods. The results clarify a commonly overlooked theoretical issue and offer a framework for integrating perturbation methods into engineering curricula. This work strengthens the conceptual foundation of a core topic in statistical thermodynamics while enhancing its pedagogical presentation.