A 1913 oscillator argument gets used to explain why cubic equations of state work

The Korea Institute of Civil Engineering and Building Technology has published a cubic equation of state that, it says, supplies a physical reason for a mathematical structure the chemical and petroleum industries have used for more than half a century without a first-principles explanation ^[1]. For anyone sizing a distillation column or a refrigeration cycle, the practical content is not a new correlation but a boundary: the workhorse form is now argued to be the minimum that satisfies three constraints, which tells you what it is doing and what it is not ^[2][8].

Since the 1970s, SRK, Peng-Robinson and Patel-Teja have improved accuracy by adopting a particular quadratic form for reshaping the attractive-force term ^[3]. That structure was largely arrived at by trial and error, and why it works so well has stayed open ^[4]. According to KICT, Jai-Yeop Lee of its Environmental Research Division traced the form to a little-known 1913 argument by the Dutch physicist Hugo Tetrode, who treated fluids as collections of vibrating oscillators rather than freely moving particles ^[5]. The work appears in Chemical Engineering Science ^[6].

Adding Tetrode's vibrational correction through a new parameter, d, the paper argues the quadratic structure is the minimal form that at once reduces correctly to the ideal-gas law at low density, remains solvable as a cubic, and keeps enough flexibility to reproduce each substance's critical compressibility ^[7][8]. That is the useful part. Three requirements, one surviving shape.

The numbers are modest and worth reading carefully. Against NIST REFPROP reference data for 76 fluids, from argon and methane to water and ammonia ^[9], the model in fully predictive mode, using only basic critical properties and no adjustable volume-translation correction, gave the lowest average error in saturated-liquid volume at 4.0%, against 4.6% for VPT, 5.5% for PT, 7.2% for PR and 13.7% for SRK ^[10]. So SRK's average error is about 3.4 times the new model's, and PR's about 1.8 times ^[11][12], but the winner is still off by 4% on average ^[10]. This is a structural result dressed in a benchmark, not a step change in liquid-density accuracy.

The parameter d also behaves like a physical quantity rather than a fudge factor. Its values correlate at R-squared of about 0.93 with an empirical constant in vapor-pressure equations whose theoretical basis had been unclear, and they sort the 76 fluids into four chemical families ^[13]. Scaled by molecular size, d rises steadily from weakly interacting argon to strongly hydrogen-bonded water, which KICT reads as a measure of interaction strength ^[14]. Lee's own framing is unusually restrained: modern cubic equations are extraordinarily useful, he said, but part of their success has rested on empirical mathematical structure rather than physical understanding ^[15].

Because the model predicts liquid and vapor behavior from critical properties without the extra empirical corrections comparable methods need, KICT presents it as a more transparent basis for process and equipment design ^[16], with validation spanning hydrogen, carbon dioxide and ammonia and suggested relevance to hydrogen energy, carbon capture and utilization, and clean-ammonia fuels ^[17].

Two things to watch. First, whether d can be assigned to a fluid nobody has measured, using the four-family grouping, rather than back-fitted. Second, mixtures: the released material reports pure-component saturated-liquid volume, and mixture performance is where cubic equations earn their keep in plant design ^[10][18].