Version v1.2 · 2026-10-14
This is a living document. It is revised as further method sets are verified in the underlying reading programme; each revision is logged at the end. Statements that the underlying sources leave unresolved are marked [to be verified]. The review records methods side by side — their forms, validity ranges, reference states and verification status — and draws no ranking among them.
Every absorption refrigerator or heat-pump calculation passes through the same narrow gate: the thermophysical properties of the working pair. A cycle model needs the solution vapour pressure Psol(T, ξ), the bubble-point and dew-point temperatures, and the specific enthalpy of the solution hsol at every state point; the generator, absorber and solution heat-exchanger balances are only as definite as the correlations behind these quantities.
Two practical difficulties recur throughout the literature surveyed here. The first is dispersion. For a single pair such as LiBr–H₂O, method sets coexist from journal papers, theses, handbooks and manufacturers' booklets spanning 1971–2025, fitted to different data batches, in different functional forms (per-property polynomials, Dühring-type relations, activity-coefficient models, single Gibbs-energy equations), and with three concentration conventions in simultaneous use (weight-percent numbers, mass fractions, mole fractions). The second is degradation in transmission: correlations are reprinted from paper to thesis to textbook, and along these chains coefficients acquire sign errors, lost exponents and altered values (Section 5); enthalpy reference states also differ between sets, so enthalpies from two sets cannot be compared — or combined — until their zero points are aligned.
This review compiles the method sets that an ongoing survey has located, transcribed and, where the sources permit, checked by recomputation. The pairs covered are listed in Section 2.1; the methods are organised by refrigerant and, within each refrigerant, by absorbent (Section 3), followed by the pure-fluid fits on which several pair models rest (Section 4), an account of the provenance and transcription-error auditing that accompanies the survey (Section 5), and a list of open gaps (Section 6).
Census. A census of the published literature (begun 2026-09-28) admits a source as a method only if it prints explicit equations with numerical constants; sources giving data tables without fitted equations are registered as data-type leads. Each admitted pair is graded: A — bubble point and enthalpy both programmable; B — vapour–liquid equilibrium (VLE) only; C — candidates, gaps and exclusions.
Verification. Four forms of checking are used, and the review states which form stands behind each entry: (i) recomputation against the source's own measured points or printed tables; (ii) cross-implementation agreement between independent Python and C implementations of the same transcription; (iii) design-point and endpoint back-checks (for example, an enthalpy polynomial whose pure-refrigerant endpoint is structurally required to vanish); (iv) same-point comparison of two method sets against a third reference (CoolProp 8.0 in the cases reported here). A distinction is maintained throughout between a source's self-reported fit residuals and residuals recomputed in this survey.
Validity ranges and reference states. Every method is recorded with the temperature, pressure and composition range of its underlying data; implementations in the associated library refuse evaluation outside these ranges rather than extrapolating silently. Enthalpy reference states are recorded per set, because they differ (Section 3, passim).
Provenance. Each set is labelled original fit, compilation, transcription or re-fit, with the basis of the judgement recorded: equation-by-equation comparison, coefficient fingerprints, endpoint structure tests, or numerical back-checks (Section 5).
The associated library. Verified sets are built into an open-source library (MIT licence) with three implementations per pair — Python (reference), C, and Excel VBA worksheet functions — under a unified naming system: Psol, Pref, TdewByTsolXi, XiByTsolTdew, hgas, hliq, hsol, hlat = hgas − hliq, TsolByHsol. The VBA modules are ported line-by-line from the verified Python/C code but have not been compiled in a real Excel session; this caveat applies to every "built" status below.
Classification. This review organises the pairs by a two-level scheme: the first level is the refrigerant, the second the absorbent. Section 2.1 lists every pair covered, in that order, and Section 3 follows the same order. A derivative complex absorbent is registered under its own name within its parent absorbent family, and a binary refrigerant blend counts as one refrigerant entry.
The pairs covered are listed below in the classification order (refrigerant, then absorbent). The last column states, in one line, the method sets on record for the pair and their verification status; the full accounts are in Section 3. Pairs for which the census found no source printing explicit equations with constants are not listed; they appear as gaps in Section 6.
Second-level pages: NH₃ · H₂O · CH₃OH · Halocarbons and hydrocarbons · Pure fluids
| Refrigerant | Absorbent | Pair | Method sets on record and verification status |
|---|---|---|---|
| NH₃ | H₂O | NH₃–H₂O | 3 sets; the principal set (Pátek & Klomfar 1995) recomputed and built |
| NH₃ | LiNO₃ | NH₃–LiNO₃ | 3 sets + the 1984 source entry; the Castilla (2014) compilation built, two self-fits in notes only |
| NH₃ | NaSCN | NH₃–NaSCN | 4 sets of one Ferreira lineage; bubble point verified at the source's design points; enthalpy coefficient sign unresolved (Section 3.1.3) |
| NH₃ | [Na(TX-7)]SCN complex | NH₃–[Na(TX-7)]SCN | Derivative system, 1 set (Sun Qiang 2020); its NRTL α is flagged in Section 3.1.3 |
| NH₃ | [emim][SCN] | NH₃/[emim][SCN] | 1 set; census grade A; not yet built |
| NH₃ | [emim][EtSO₄] | NH₃/[emim][EtSO₄] | 1 set; census grade A; not yet built |
| NH₃ | [emim]Cu₂Cl₅ | NH₃/[emim]Cu₂Cl₅ | 1 set (Chen & Bai 2016); census grade A; not yet built |
| NH₃ | [bmim]Zn₂Cl₅ | NH₃/[bmim]Zn₂Cl₅ | 1 set (Chen et al. 2013); grade B — VLE only |
| H₂O | LiBr | LiBr–H₂O | 10 sets; two built, one verified at the source's own points (not yet built) |
| H₂O | [emim][DMP] | H₂O/[emim][DMP] | 1 set (2018); grade B — VLE only |
| CH₃OH | LiBr | CH₃OH–LiBr | Appendix compilation (Critoph & Rivero Pacho 2022), after Saravanan & Maiya 1998 and Iyoki et al. 1991; misprints TE-020–TE-022 |
| CH₃OH | LiBr + ZnCl₂ (1:1) | CH₃OH–(LiBr + ZnCl₂, 1:1) | Appendix compilation (2022) |
| CH₃OH | LiBr + ZnBr₂ (2:1) | CH₃OH–(LiBr + ZnBr₂, 2:1) | Appendix compilation (2022); the same pair also built from El-Shamarka (1981) and Iedema (1984) |
| CH₃OH | LiI + ZnBr₂ (2:1) | CH₃OH–(LiI + ZnBr₂, 2:1) | Appendix compilation (2022) |
| CH₃OH | ZnBr₂ | CH₃OH–ZnBr₂ | Appendix compilation (2022), after Iedema (1984) |
| R134a | DMF | R134a–DMF | Zehioua et al. (2009) + the library's γ–φ re-fit; a second published set is excluded (TE-024) |
| R134a | DMA | R134a–DMA | Nezu et al. (2002); the Table 2 B-coefficient signs are suspect (TE-028, Section 3.4.1) |
| R124 | NMP | R124–NMP | Xu et al. (2017); VLE only — no enthalpy published (Section 6) |
| R124 | DMF | R124–DMF | Xu et al. (2017); VLE only — no enthalpy published (Section 6) |
| R124 | DMAC | R124–DMAC | Borde et al. (1997), a complete set; transcription verified against the paper's endpoints and figures |
| R22 | DMF | R22–DMF | Agarwal (1982) VLE fit + Fatouh (1993) enthalpy; built |
| R22 | DEGDME | R22–DEGDME | Ando & Takeshita (1984); the paper-reported deviations are reproduced by the library implementation |
| R161 | DMETrEG | R161–DMETrEG | Feng et al. (2016) five-parameter NRTL; one coefficient coincidence is flagged in Section 3.4.4 |
| R152a | DMETrEG | R152a–DMETrEG | Feng et al. (2016) five-parameter NRTL |
| R152a | NMP | R152a–NMP | Feng et al. (2016); α = −37.1878 flagged in Section 3.4.4 |
| (R32 + R227ea) | DMF | (R32 + R227ea)/DMF | Gao Wanli et al. (2010); no measured VLE published and no solution-enthalpy model — no experimental anchor |
| R600a | POE ISO 7 oil | R600a + POE oil | Neto & Barbosa (2008); Raoult-form reduction AAD ≈ 5% |
| R290 | Oil (unspecified) | R290–oil | Wang, Jia & Wang (2021) p–h–X chart; full text not yet obtained |
| R290 | NM100P oil | R290–NM100P | Li Keqiao (2018) thesis fit; printed form and units await checking (Section 3.4.6) |
| — | — | Pure fluids (nine fluids) | Saturation fits for H₂O, NH₃, CH₃OH, R134a, R22, R124, R32, R218 and CO₂ (Section 4) |
Overview. Table 1 collects the principal sets. It is maintained as an extensible table: newly verified sets are added as rows in the same columns, in the classification order of Section 2.1, and a row's verification status is upgraded as checking progresses.
Table 1. Principal method sets by refrigerant and absorbent.
| Pair | Method source | Correlation type | Validity range | Verification status |
|---|---|---|---|---|
| NH₃–H₂O | Pátek & Klomfar 1995 | Five explicit polynomial correlations | 0.002–2 MPa, full composition | Recomputed; built |
| NH₃–H₂O | Conde 2004 (transport) | Corresponding states + mixing rules | Saturation line | Spot-checked vs. CoolProp; built |
| NH₃–H₂O | Schulz 1971 | Gibbs-function equation of state | 0.001–2.5 MPa, 200–450 K | Coefficients collated; faces [to be verified] |
| NH₃–LiNO₃ | Castilla 2014, App. B (Libotean et al.) | Per-property polynomials | 293.15–353.15 K, x 0.35–0.65 | Built; C/Python 1715 pts; viscosity excluded |
| NH₃–LiNO₃ | Cai Dehua 2015 | ln p form + linear cp, self-fit | 293.15–353.15 K, 0.2 < w < 0.6 | Source-reported residuals; notes only |
| NH₃–LiNO₃ / NH₃–NaSCN | Liang Xiao 2022 | Unified cubic template, re-fit | Libotean / Chaudhari data domains | Source-reported RMSE; notes only |
| NH₃–NaSCN | Ferreira lineage, via Yin Jiayuan 2015 | ln P = A + B/T; enthalpy cubic in t | Verified band X ≈ 0.45–0.52, 310–355 K | Bubble verified at design points; enthalpy C-term sign [to be verified] |
| NH₃–NaSCN | Cai Dehua 2015 | Self-fit + double-reciprocal heat-of-solution model | 253.1–393.1 K, 0.35 < w < 0.9 (vapour pressure) | Source-reported; notes only |
| NH₃/[emim]Cu₂Cl₅ | Chen & Bai 2016 | Modified UNIFAC (Dortmund) + enthalpy equations | xNH₃ 0.89–0.96 | Census grade A; not yet built |
| LiBr–H₂O | Patterson & Perez-Blanco 1988 | Unified double polynomials, 7 functions | X ≤ 70 wt%, T ≤ 180 °C (per function) | Built; C/Python point-identical |
| LiBr–H₂O | Chua et al. 2000 | Dühring equation + density polynomial | p: 0–190 °C, x ≤ 60 wt%; ρ: 0–200 °C, ≤ 70 wt% | Built; checked against paper's figures |
| LiBr–H₂O | Pátek & Klomfar 2006 | Mole-fraction correlations | 273–500 K, 0–75 wt% | Verified at six paper points; not yet built |
| LiBr–H₂O | Kaita set, via Li Jing 2009 | Dew-point polynomial + Antoine-type p | 20–210 °C, 40–65 wt% | Cross-checked; attribution [to be verified] |
| LiBr–H₂O | Kim & Ferreira 2006 | Single Gibbs free-energy equation | Curve family 0–210 °C, 0–75 wt% | Constants collated; faces/units [to be verified] |
| CH₃OH–(LiBr–ZnBr₂ 2:1) | El-Shamarka 1981; Iedema 1984 | Antoine-type + enthalpy polynomials | As per the two theses | Built; sources differ ≈20% on p at w ≈ 0.35–0.40 |
| R134a–DMF | Zehioua et al. 2009; re-fit in this library | PR + MHV1 + NRTL (paper); NRTL γ–φ (re-fit) | 303.30–353.24 K | Recomputed, AARD 1.49% (58 pts) |
| R124–NMP / R124–DMF | Xu et al. 2017 | Five-parameter NRTL | 303.15–363.15 K | Recomputed, AARD 1.25% / 1.17%; VLE only |
| R124–DMAC | Borde et al. 1997 | Full polynomial correlation set | As published (see §3.5) | Transcription verified vs. paper's endpoints/figures |
| R22–DMF | Agarwal 1982 (VLE fit); Fatouh 1993 (enthalpy) | NRTL-class fit + enthalpy correlation | Lower boundary −25 °C | Fit AARD 2.38% (132 pts); built |
Because these correlations circulate by reprinting, the survey maintains a standing register of transcription and printing errors (entries TE-001 to TE-076 at the time of writing), each with the wrong and correct values, the detection method, the quantitative impact, and the disposition. Four detection methods recur: coefficient-by-coefficient comparison of two chains; endpoint structure tests; back-calculation from the source's own tables or figures; and numerical fingerprinting. Three confirmed cases and one falsification illustrate the scale of the effects.
TE-001 (transcription chain). In the SCIRP chain of the Ferreira NH₃–NaSCN set, the X³ coefficient of the enthalpy polynomial A is printed as −3.5137; the parallel Yin chain gives −295.67. The endpoint test decides: the Yin chain gives A(1) = −0.05 ≈ 0, the SCIRP chain A(1) = +292.11, which is physically impossible; and −3.5137 is identical, digit for digit, to the SCIRP chain's own B-polynomial X³ coefficient — the coefficient was copied from the wrong polynomial. The chain's enthalpies run 33–52 kJ/kg high at six comparison points; its enthalpy equation is excluded, while its bubble-point equation (within +0.014 to +0.016% of the Yin chain at six points) remains usable as corroboration.
TE-014 (original-paper misprint). In Pátek & Klomfar (2006), the 22nd of the 30 enthalpy coefficients is printed as 2.33910E1; back-calculation at the paper's own six Table 9 points reproduces both p and h only with 2.33910E−1 — a factor-of-100 exponent error in the original printing.
TE-015 (original-paper misprint). In Chua et al. (2000), Table 1, coefficient B₁₂ is printed as +3.12628E−16 without its minus sign (the glyph is clean at 800 dpi, so this is not a reading error). With the positive value, BD(30) = 336.8 against the ≈4.5 implied by the paper's Fig. 2; with the negative value Fig. 2 is reproduced. The library implements −3.12628E−16.
TE-050 (non-self-consistent parameter set). A 2014 thesis reports new VLE measurements for R161 + DMAC, R161 + NMP, R134a + DMAC and R134a + NMP with five-parameter NRTL correlations. Substituting the printed equations and parameters at one of its own points (R161 + DMAC, T = 303.15 K, x₁ = 0.4058; measured p = 0.3315 MPa) returns a bubble pressure of ≈0.000078 MPa — a factor of about 4000 from the measurement; the printed α values (5.935–18.97) also lie far outside the conventional NRTL span of 0.2–0.5. The parameter sets are marked [to be verified] and excluded from building; the measurement tables remain usable. The register exists so that such sets are documented and excluded rather than silently "corrected".
Bibliographic fields not present in the survey materials are marked [bibliographic detail pending] rather than completed from memory.
| Version | Date | Additions | Source batch |
|---|---|---|---|
| v1 | 2026-10-14 | Initial compilation: NH₃–H₂O (3 sets + endpoint account), LiBr–H₂O (10 sets), NH₃–LiNO₃ (3 sets + source entry), NH₃–NaSCN (4 sets + 1 derivative system), refrigerant–solvent sets, ionic-liquid and methanol–salt sets, nine pure fluids, transcription-error register summary | Compiled from the existing method-comparison drafts (NH₃/H₂O and ammonia salts; LiBr–H₂O; nine pure-fluid drafts), the 2026-09-28 working-pair census, and the open-source working-pair library's verification records |
| v1.1 | 2026-10-14 | Method sections re-organised by refrigerant (first level) and absorbent (second level); list of pairs covered added as Section 2.1; Table 1 re-ordered to the same scheme; HTML edition split into an overview page plus one page per refrigerant (pure fluids on their own page). No method content, figure or status changed | User direction, 2026-10-14 |
| v1.2 | 2026-10-14 | Pressure units unified to MPa throughout: every pressure value converted exactly from its previous unit (kPa/bar), matching the working-pair library's own MPa unification of the same date. Values explicitly quoted as printed in a source keep the source unit and are marked as printed. No method content, figure or status changed | Library-wide MPa unification, 2026-10-14 |