{"id":"80837298-01da-4848-8f4a-25fc7eaec86c","arxiv_id":"2506.01082","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"New density measurements for a high-calorific natural gas with 9% ethane and 3% propane, and its 10% and 20% hydrogen blends, reveal equation-of-state deviations up to 0.44%.","lead":"This paper reports new, precisely measured densities of a hydrogen-enriched natural gas mixture with significant ethane and propane content, at temperatures from 260 to 350 K and pressures up to 20 MPa. The data show that standard equations of state predict these densities less accurately than they do for methane-rich natural gas, with deviations up to 0.44%.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified: the density-based central claim is robust; only a secondary compressibility inconsistency needs correcting.","rationale":"I read the paper as an experimental data contribution whose central comparative claim is well supported. The mixtures were gravimetrically prepared under ISO 6142-1 with independent GC validation, and the densimeter uncertainty budget is detailed. Tables 3 to 5 provide state-by-state densities and deviations, and Table 6 shows consistently larger AARD and MaxRD values for the ethane/propane-rich G432/G455/G456 set than for the methane-rich G431/G453/G454 set across all three equations of state. The largest deviations occur at high pressure and low temperature, where the experimental relative uncertainty is roughly 0.02% to 0.03%; a 0.44% deviation is therefore a genuine model-experiment discrepancy, not an artifact of low-pressure sorption. The sorption limitation flagged by the reader is real for the lowest-pressure points, but it is explicitly disclosed and is not the basis for the central conclusion. The only notable internal inconsistency I found is in Section 4.2: the text reports MaxRD values up to 1.9% for the derived compressibility comparison while claiming that all deviations fall within Ur(kappa_T) = 0.7%. These statements are mutually incompatible. This is a secondary presentation or analysis error in the derived-property section and does not affect the density comparisons that form the paper's main conclusion. The manuscript should be accepted; the compressibility section should be corrected before final publication.","tokens_in":31139,"tokens_out":14096,"duration_ms":157140,"concrete_test":"Recompute the Table 6 AARD/MaxRD statistics for G432, G455, and G456 with the 260 K isotherm removed, and independently rederive the Section 4.2 partial derivatives and kappa_T deviations from Tables 3 to 5; if the qualitative ordering against the methane-rich mixtures persists and the reported MaxRD values are reconciled with Ur(kappa_T) = 0.7%, the paper's claims are fully secure.","verdict_should_be":"UNCHANGED","load_bearing_attack":"No load-bearing concern identified. The central claim is directly supported by the AARD and MaxRD statistics in Table 6, and the largest deviations (0.19% to 0.44%) are many times the experimental expanded uncertainty at the same state points, so they cannot be dismissed as measurement noise or as the low-pressure sorption effect acknowledged in Section 4.1. The reader's sorption concern is a plausible limitation for low-pressure data, but it does not explain the high-pressure, low-temperature deviations on which the main conclusion rests.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports new (p, ρ, T) density measurements for a high-calorific natural gas mixture containing about 9 % ethane and 3 % propane, and for two hydrogen-enriched derivatives with nominal 10 % and 20 % hydrogen content, over the range (260–350) K and pressures up to 20 MPa, using a single-sinker magnetic suspension densimeter. The experimental densities are compared with three equations of state: AGA8-DC92, GERG-2008, and an improved GERG-2008 variant. The central finding is that all three models reproduce the densities of the methane-rich mixtures from the authors' prior work better than those of the ethane/propane-rich mixtures studied here, with maximum deviations reaching 0.44 % for GERG-2008 on the 10 % hydrogen blend. The data tables and uncertainty budgets are presented in detail, and the statistical comparison with earlier literature is systematic.","tokens_in":31171,"tokens_out":5844,"duration_ms":55588,"significance":"The paper provides valuable reference-quality density data for hydrogen-enriched natural gas with realistic heavier-hydrocarbon content, a topic of direct relevance to EoS validation for gas-grid decarbonization scenarios. The mixture preparation follows ISO 6142-1 with independent GC validation; the uncertainty budget (Eqs. 2–3) is detailed; complete data tables are given; and the comparison with the prior methane-rich study is quantitative. The main conclusion is robust: the largest deviations occur at high pressure and low temperature, where they exceed the stated experimental expanded uncertainty by a wide margin, so they cannot be dismissed as measurement noise or as the low-pressure sorption effect discussed in Section 4.1. The compressibility section contains an internal inconsistency (see minor comments) that does not affect the density-based central claim.","major_comments":[],"minor_comments":[{"comment":"The statement that sorption-induced deviations 'are well below the experimental uncertainty' is not backed by a quantitative estimate; the paper correctly states that the single-sinker densimeter cannot quantify this influence, so this sentence should be reworded as an explicit assumption or supported by an order-of-magnitude estimate based on the Richter and Kleinrahm analysis.","section":"Section 4.1"},{"comment":"The text reports MaxRD values of 1.9 %, 1.6 %, and 1.3 % for the derivative comparisons and then states that 'all deviations are located within the estimated expanded (k = 2) uncertainty of κT, Ur(κT) = 0.7 %'; since the MaxRD values exceed 0.7 %, this is internally inconsistent and needs to be corrected or rephrased.","section":"Section 4.2 and Table 7"},{"comment":"The sentence 'The κT values span from (0.0361 to 0.5427) MPa-1 at 250 K' appears to refer to Table 7, which starts at 260 K; the temperature should be 260 K or the table should include 250 K data.","section":"Section 4.2"},{"comment":"The text uses the abbreviation 'MarRD' where 'MaxRD' is meant, and in Table 6 the GERG-2008 BiasRD entry for G456 appears as '0.003 9' with an awkward line break; these should be fixed.","section":"Section 4.1 and Table 6"},{"comment":"The 'improved GERG-2008' is described as combining updates from references [53] and [54], but the manuscript does not state whether the calculations were performed with REFPROP 10 or a modified in-house code, nor which binary parameters were changed; a short implementation note or reference to the supplement would improve reproducibility.","section":"Section 3.2"},{"comment":"In the normalized composition block for G456, the rows for nitrogen, carbon dioxide, ethane, and propane appear in an order that does not match the original table; please verify that the values and uncertainties are aligned with the correct components.","section":"Table 1"}],"recommendation":"minor_revision","confidential_remarks":"The paper is a solid experimental contribution well within the journal's scope. The density data and their comparison with equations of state are thoroughly documented, and the central claim is well supported. The requested changes are local: fix the compressibility inconsistency, clarify the sorption assumption, and correct a few presentation details. I see no novelty or attribution concerns."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Straight to it: this is a useful, well-executed measurement paper. The new contribution is concrete: (p,ρ,T) data for a high-calorific 11-component natural gas with 9% ethane and 3% propane, and for the 10% and 20% H2 derivatives, measured with a calibrated single-sinker densimeter over 260–350 K and up to 20 MPa. Mixtures are gravimetrically prepared per ISO 6142-1, compositions independently validated by GC, and the uncertainty budget in Eqs. 2–3 is transparent. The data fill a gap: most prior work, including the authors' own earlier study [32], used methane-rich (>97%) base gas.\n\nThe central claim—that AGA8-DC92, GERG-2008, and the improved GERG-2008 all lose accuracy when the base gas has significant ethane and propane, with deviations up to 0.44% against a stated 0.1% uncertainty—is supported by Table 6 and the pointwise tables. The high-pressure, low-temperature deviations are many times the experimental expanded uncertainty, so they cannot be dismissed as the low-pressure sorption effect the authors acknowledge in Section 4.1. I agree with the stress-test note: the sorption concern is a real limitation for the 1–2 MPa region, but it does not touch the main conclusion.\n\nSoft spots, in proportion. The sorption influence is discussed only qualitatively; a bound, even approximate, would make the low-pressure comparison more useful. The \"improved GERG-2008\" is a composite of two published modifications (EOS-LNG and Beckmüller hydrogen-mixture updates), so it is not a single testable equation; the comparison is still legitimate, but readers should not treat it as an off-the-shelf model. There is also a small inconsistency in Section 4.2 / Table 7: the text says κT values span 0.0361–0.5427 MPa⁻¹ at 250 K for G432, but this study's data start at 260 K and the table is for 260–350 K. That looks like a leftover from the 250 K baseline of [32] and should be corrected. Minor, not load-bearing.\n\nThe citation pattern is fine. The paper cites its own force-transmission calibration work [45] appropriately, and the REFPROP-based derivatives in Eq. 3 are standard uncertainty propagation, not fitting to the data being tested.\n\nBottom line: this is reference-data work that deserves a serious referee. It is not glamorous, but it is the kind of careful experimental work that EoS development and custody-transfer standardization rely on. I would cite it, and I would bring it to a reading group focused on hydrogen-blend metrology.","headline":"Solid new density data for H2-enriched high-calorific natural gas; the central EoS-degradation claim holds, with a minor compressibility inconsistency to fix.","tokens_in":31745,"tokens_out":2747,"would_cite":true,"duration_ms":23338,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["51.30.+i"],"model":"deepseek-v4-flash","headline":"High-precision density data on ethane-rich natural gas with 10% and 20% hydrogen show that all three reference equations of state predict densities less accurately than for methane-dominant gas, deviating up to 0.44% at low temperatures…","keywords":["hydrogen-enriched natural gas","single-sinker densimeter","natural gas density","equation of state","GERG-2008","AGA8-DC92","gravimetric mixture preparation","isothermal compressibility"],"falsifier":"Measure the same three gravimetric mixtures in a two-sinker densimeter, which compensates adsorption effects, or after repeated evacuation and flushing at 260 K and pressures from 1 to 20 MPa; if the 0.2 to 0.44% deviations persist, they are true equation-of-state errors, while if they shrink, part of the reported deviation is a sorption artifact.","tokens_in":30940,"feed_emoji":"🧪","tokens_out":5453,"duration_ms":47973,"temperature":0.7,"pith_summary":"This paper reports new high-precision density measurements for a synthetic high-calorific natural gas containing 85% methane, 9% ethane, and 3% propane, plus two hydrogen-enriched derivatives with 10% and 20% hydrogen. The measurements, made with a single-sinker densimeter at 260 to 350 K and pressures up to 20 MPa, are compared with three reference equations of state: AGA8-DC92, GERG-2008, and an improved GERG-2008. The central finding is that all three models predict densities less accurately for these ethane- and propane-rich mixtures than for the methane-dominant mixtures studied earlier, with deviations growing at low temperature and high pressure and reaching 0.44% for GERG-2008 on the 10% hydrogen mixture. This matters because custody transfer, pipeline operation, and hydrogen-grid blending rely on these equations to convert pressure and temperature into density and energy content; if the models overstate their accuracy on realistic heavier natural gas, billing and capacity calculations inherit the error.","feed_headline":"Hydrogen-rich gas data: equations of state off by up to 0.44%","feed_subtitle":"Ethane-rich hydrogen blends reveal up to 0.44% density errors in all three reference models.","key_machinery":"The load-bearing instrument is a single-sinker magnetic-suspension densimeter, whose density equation combines sinker buoyancy with corrections for a force-transmission error split into an apparatus-specific factor and a fluid-specific magnetic-susceptibility term. The mixtures themselves are the other carrier of the argument: they were prepared gravimetrically under ISO 6142-1 and validated by gas chromatography, so the composition is known to roughly 0.001 to 0.005 mol% and the density deviations cannot be blamed on preparation. The three equations of state — AGA8-DC92, GERG-2008, and an improved GERG-2008 built from updated pure-fluid and binary departure functions — supply the predictions against which the measured densities are compared. The statistical engine is a set of relative-deviation metrics (AARD, BiasRD, RMSRD, MaxRD) computed at every state point.","core_discovery":"On the paper's own terms, the discovery is a documented accuracy boundary: all three reference equations of state perform worse for a high-calorific natural gas with significant ethane (9%) and propane (3%) content than for the methane-rich (>97%) gas tested in the companion study, both with and without hydrogen addition. For the hydrogen-free mixture G 432, AGA8-DC92 stays within its claimed 0.1% uncertainty except near 260 K and 2 to 10 MPa, while GERG-2008 and the improved GERG-2008 deviate up to 0.21% and 0.28%. For the 10% hydrogen mixture G 455, GERG-2008 reaches -0.44% and the improved version -0.38%; for the 20% hydrogen mixture G 456, maximum deviations are 0.28% (AGA8-DC92), 0.34% (GERG-2008), and 0.30% (improved GERG-2008). The paper concludes that equation-of-state accuracy decreases for mixtures that deviate from simple methane-rich compositions, especially at lower temperatures and higher pressures.","pith_inferences":["If the error scales with ethane and propane content, real grid gas (which often contains 5-15% ethane) may be systematically miscalculated; testing intermediate ethane levels would show whether the trend is monotonic.","The nonzero low-pressure deviations suggest sorption may contaminate low-density data points; a quantitative sorption correction could shrink the apparent equation-of-state errors.","The non-monotonic hydrogen effect (10% hydrogen worse than 20% for the GERG models) hints at cancellation between opposing binary-model errors, which could guide future parameterization of hydrogen and hydrocarbon interactions.","Re-analyzing the cell gas composition by gas chromatography after the density runs would reveal whether adsorption or desorption altered the effective mixture, separating experimental artifacts from genuine model deficiencies."],"forward_implications":["For custody transfer and pipeline metering, using AGA8-DC92 or GERG-2008 on ethane-rich hydrogen-enriched gas at low temperature and high pressure can produce density errors several times larger than the models' claimed 0.1% uncertainty.","The improved GERG-2008 does not fix the problem: on the 10% hydrogen mixture it is among the worse performers, deviating by -0.38% at 260 K near 13 MPa.","The new data give reference-quality targets for refitting binary departure functions, particularly for methane-ethane, methane-propane, and hydrogen-heavy-hydrocarbon interactions.","AGA8-DC92 remains the most consistent of the three models on these mixtures, but its deviations still grow with hydrogen content, from 0.15% maximum on the hydrogen-free gas to 0.28% on the 20% hydrogen gas."],"supporting_citations":[{"why":"Defines the AGA8-DC92 equation of state whose density predictions are tested against the new data.","marker":"[19]"},{"why":"Defines the GERG-2008 equation of state whose density predictions are tested.","marker":"[20,21]"},{"why":"Supplies the methane-rich natural gas and hydrogen-enriched density data from the companion study that the present deviations are compared against.","marker":"[32]"},{"why":"The ISO standard governing the gravimetric preparation of the reference mixtures, central to the claimed composition accuracy.","marker":"[33]"},{"why":"ISO 20765-2 formulation that extends GERG-2008 and anchors the improved-GERG comparison.","marker":"[52]"},{"why":"Provides the LNG-range departure-function update that is one of the two modifications comprising the improved GERG-2008.","marker":"[53]"},{"why":"Provides updated pure-fluid and hydrogen-binary departure functions, the other modification comprising the improved GERG-2008.","marker":"[54]"},{"why":"Examines adsorption and desorption influence on gas-mixture density measurements, the study the authors cite to frame the sorption caveat.","marker":"[58]"}],"fun_headline_variants":["Density errors up to 0.44% in ethane-rich hydrogen gas blends","Equations of state miss hydrogen-rich gas density by 0.44%","High-calorific natural gas exposes EoS limits up to 0.44%","Ethane-heavy gas blends strain gas equations of state"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument assumes that gas adsorption and desorption on the densimeter cell walls do not measurably change the sample composition or the measured density, especially at low pressures where the largest relative uncertainties occur.","fun_headline_variants_meta":{"raw":{"variants":["Density errors up to 0.44% in ethane-rich hydrogen gas blends","Equations of state miss hydrogen-rich gas density by 0.44%","High-calorific natural gas exposes EoS limits up to 0.44%","Ethane-heavy gas blends strain gas equations of state"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000327,"raw_usage":{"total_tokens":1856,"prompt_tokens":998,"completion_tokens":858,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":614,"completion_tokens_details":{"reasoning_tokens":776}},"tokens_in":614,"tokens_out":858,"duration_ms":6915,"temperature":1.0,"reasoning_tokens":776,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T11:51:05.037886+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the same three gravimetric mixtures in a two-sinker densimeter, which compensates adsorption effects, or after repeated evacuation and flushing at 260 K and pressures from 1 to 20 MPa; if the 0.2 to 0.44% deviations persist, they are true equation-of-state errors, while if they shrink, part of the reported deviation is a sorption artifact.","supporting_citations":[{"cited_title":"AGA Report No","cited_arxiv_id":null,"evidence_quote":"Defines the AGA8-DC92 equation of state whose density predictions are tested against the new data."},{"cited_title":"Thermodynamic (p, ρ, T) characterization of a reference high-calorific NG mixture when hydrogen is added up to 20 % (mol/mol)","cited_arxiv_id":null,"evidence_quote":"Supplies the methane-rich natural gas and hydrogen-enriched density data from the companion study that the present deviations are compared against."},{"cited_title":"ISO 6142-1 Gas analysis — Preparation of calibration gas mixtures — Part 1: Gravimetric method for Class I mixtures","cited_arxiv_id":null,"evidence_quote":"The ISO standard governing the gravimetric preparation of the reference mixtures, central to the claimed composition accuracy."},{"cited_title":"ISO 20765-2 NG — Calculation of thermodynamic properties — Part 2: Single-phase properties (gas, liquid, and dense fluid) for extended ranges of application","cited_arxiv_id":null,"evidence_quote":"ISO 20765-2 formulation that extends GERG-2008 and anchors the improved-GERG comparison."},{"cited_title":"EOS-LNG: A Fundamental Equation of State for the Calculation of Thermodynamic Properties of Liquefied NGes","cited_arxiv_id":null,"evidence_quote":"Provides the LNG-range departure-function update that is one of the two modifications comprising the improved GERG-2008."},{"cited_title":"Influence of adsorption and desorption on accurate density measurements of gas mixtures","cited_arxiv_id":null,"evidence_quote":"Examines adsorption and desorption influence on gas-mixture density measurements, the study the authors cite to frame the sorption caveat."}],"review_version":1}