{"id":"766ed045-9633-4053-9f0a-adbfef730c92","arxiv_id":"2411.17742","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"The paper derives a composition compressibility k = beta(T, mu_t) - beta(T, y) from the Gibbs function and shows how it distinguishes multicomponent incompressible fluids from constant-density fluids.","lead":"This paper introduces a new coefficient, k, that measures how much a multicomponent fluid's volume changes with composition, expressed as a compressibility comparable to the usual isothermal compressibility. It shows that for incompressible mixtures, k can be large even when pressure-driven compression vanishes, and illustrates this with water-salt, water-sucrose, and water-ethanol solutions.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The exact decomposition (3) is sound only where the compositional Hessian is positive definite, and the paper never verifies this for the fitted ternary water-ethanol model, so the reported k ratios may be artifacts rather than validated compressibilities.","rationale":"The reader's weakest assumption identifies the need for a positive-definite compositional Hessian, which is indeed a necessary condition for the central decomposition. I agree that this is the theoretical hinge of the paper, but the more acute form of the same concern is that the paper never verifies this condition for its own fitted models, especially the ternary water-ethanol model in Appendix C, where chemical equilibrium is imposed and only two parameters are fitted. If the effective Hessian loses positive definiteness somewhere, eqs. (54)–(55) fail exactly where the paper draws quantitative conclusions. Even in the stable regime, the ideal-model fits are not independently validated for k: the seawater comparison shows sizable disagreement at high salinity, and no equivalent validation is offered for sucrose or ethanol. The theoretical claim itself, eq. (3), appears mathematically clean and is supported by the general derivation in §3.2 and by the ideal-mixture formula (57), so I do not see a reason to reject or reframe the paper's core identity. The weakness is in the quantitative application to real fluids, which is precisely why the reader's CONDITIONAL verdict is appropriate. My read does not change that verdict; it sharpens the condition: before the ratio k/β(T,y) is used as a criterion, the fitted models must be checked for stability and validated against full thermodynamic data.","tokens_in":18474,"tokens_out":11628,"duration_ms":118702,"concrete_test":"For the ternary water-ethanol model of Appendix C, construct the effective Gibbs function g(T, p, y) with the fitted parameters, and compute the smallest eigenvalue of D²_yy g over y ∈ [0,1] at T = 25°C and p = 1 atm. Separately, evaluate k from eq. (55) using numerical derivatives of a validated non-ideal Gibbs function for ethanol-water (or from experimental partial-molar-volume data) across the same composition range and compare with the values behind Fig. 6. If the minimal eigenvalue is non-positive or the independent k differs by more than ~50% at any composition, the reported ratio k/β(T,y) is not a reliable quantitative measure for this mixture.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central identity β(T, μ_t) = β(T, y) + k, eq. (3), is a calculus result that requires the map (T, p, y) ↦ (T, p, μ_t) to be a local diffeomorphism, i.e. it requires D²_yy g to be positive definite. The paper states this stability condition in §2.1, but it never checks it for the fitted models used in §4, and the water-ethanol ternary model of Appendix C is especially exposed. There, g is an effective Gibbs function obtained after imposing the chemical equilibrium condition μ3 = (1−θ)μ1 + θμ2 and solving the algebraic system (79). The effective reduced Hessian in y is not computed; the parameters D(T), D'(T) are fitted to one reference density and one compressibility, so nothing guarantees convexity on the plotted composition range. If an eigenvalue of D²_yy g approaches zero, the inverse in eqs. (54)–(55) diverges, k ceases to be a physical compressibility, and the ratios shown in Figs. 5–6 become artefacts. Even where the Hessian is positive, the quantitative estimates of k come from ideal models fitted to limited data; the one direct comparison against a full Gibbs library (seawater, Fig. 2) shows order-level disagreement at high salinity. The paper itself acknowledges in the conclusion that wider data sets and full Gibbs functions are needed. Thus the theoretical decomposition is internally consistent, but the paper's practical claim that k/β(T,y) provides a quantitative criterion for real multicomponent fluids is not yet established.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces a coefficient k defined as the difference between the isothermal compressibility at fixed chemical potentials, β(T,μ_t), and the usual isothermal, isocompositional compressibility β(T,y). The central identity is β(T,μ_t) = β(T,y) + k, where k is expressed via the Gibbs function and the inverse of its compositional Hessian (Eqs. (2), (54), (55)). The author proves this identity by transforming variables from (T,p,y) to (T,p,μ_t) and argues that in the incompressible limit β(T,y)→0, k can remain nonzero, so an incompressible multicomponent fluid need not have constant density. The ratio k/β(T,y) is proposed as a quantitative criterion for the validity of the constant-density approximation. Three aqueous solutions (seawater, sucrose, and ethanol) are used as illustrations, with the ethanol case treated by both a binary ideal model and a ternary ideal model with a chemical equilibrium constraint (Appendix C).","tokens_in":18844,"tokens_out":4790,"duration_ms":42840,"significance":"If the decomposition is valid and can be computed reliably, the paper provides a useful conceptual distinction between 'incompressible' and 'constant density' for multicomponent fluids, with a concrete criterion for when the constant-density approximation fails. The calculus derivation in §3.2 is transparent and internally consistent, and the explicit ideal-mixture formula (57) is a clean closed-form result. The paper is also honest in acknowledging that wider data sets and full Gibbs functions are needed (conclusion). However, the numerical support is not yet sufficient to establish the quantitative claims: the seawater ideal model has a 27.9% error in the very compressibility used in the ratio, the sucrose and ethanol ideal models have errors of 7.2% and 35% respectively, and the ternary water-ethanol model is fitted with only two parameters and its thermodynamic stability is never verified. Thus the theoretical decomposition is sound, but the practical claim that k/β(T,y) provides a quantitative criterion for real multicomponent fluids is not established.","major_comments":[{"comment":"The decomposition β(T,μ_t) = β(T,y)+k requires the compositional Hessian D²_yy g to be positive definite (hence invertible) on the relevant domain. The paper states this stability condition in §2.1 but never verifies it for the fitted models used in §4. In the ternary water-ethanol model of Appendix C, g is an effective Gibbs function obtained after imposing chemical equilibrium (Eq. (79)); the effective reduced Hessian is not computed, and the parameters D(T), D′(T) are fitted to one reference density and one reference compressibility (Eqs. (81)–(82)). If an eigenvalue of D²_yy g approaches zero or the effective model is not convex on the plotted composition range, the inverse in Eqs. (54)–(55) fails or becomes indefinite, so k loses its interpretation as a physical compressibility and the ratios shown in Figs. 5–6 become artefacts. This issue is load-bearing because the proposed criterion k/β(T,y) presupposes that k is well defined.","section":"§3.2, Eqs. (53)–(55)"},{"comment":"The ideal binary model for seawater has relative error Err(β(T,y)) = 0.279 in the isothermal compressibility, yet this model is used to compute k = β(q)−β(T,y) and to compare it with β(T,y). Figure 2 shows that for high salinity the ideal model yields a k about six times larger than β(T,y), while the Gibbs seawater library attains a significantly smaller maximum. With a 27.9% error in the very quantity that defines the ratio k/β(T,y), the quantitative estimates for seawater are not reliable. The paper should either compute k from the full Gibbs function provided by the gsw library (as it already does for comparison) or provide uncertainty bounds before drawing conclusions about the incompressible-limit criterion.","section":"§4.1, Eq. (65) and following"},{"comment":"The same reliability issue affects the sucrose and ethanol examples. For sucrose, Err(β(T,y)) = 0.072 (§4.2); for the binary ethanol model, Err(β(T,y)) = 0.35 and Err(ϱ) = 0.14 (§4.3). The ternary ethanol model improves the fits but is calibrated with only two parameters D(T) and D′(T) to one reference density and one reference compressibility (Eqs. (81)–(82)). Consequently, the reported maximum ratios (more than 10 for sucrose, slightly above 1 for ethanol) are model outcomes, not validated physical predictions. The paper's own caveats in §4.1 ('Further investigations will be necessary') and in the conclusion ('wider data sets and ... full Gibbs functions are available') are appropriate, but they are in tension with the abstract's statement that the phenomenon becomes quantitatively comparable through k.","section":"§4.2 and §4.3"}],"minor_comments":[{"comment":"The notation D²_{y,y}g is used in Eq. (2) before the tangential derivative is defined in §2.1; please define the notation at first use.","section":"Eq. (2) and §2.1"},{"comment":"The phrase 'completely equivalent ot (8)' contains a typo: 'ot' should be 'to'.","section":"Page 5, first paragraph"},{"comment":"The sentence 'We let ϱs be the bulk density of the mixture entropy possesses of the special form' is grammatically awkward and should be rephrased.","section":"§2.1, paragraph on entropy"},{"comment":"The error metric Err(·) is defined as the maximum absolute deviation normalized by the data range; this makes the values Err(ϱ)=0.029 and Err(β)=0.279 difficult to interpret. Please also report a typical or root-mean-square error, which would give a clearer picture of the model's accuracy.","section":"§4.1, Eq. (65)"},{"comment":"The wave-equation analysis is only linearized and heuristic; the statement that diffusion processes 'accelerate' compression waves (around Eq. (73)) should be qualified as a linearized, near-equilibrium result.","section":"Appendix B"}],"recommendation":"major_revision","confidential_remarks":"The paper's theoretical core is a calculus identity that is correctly derived under the stated stability assumptions, and the ideal-mixture formula (57) is a valuable explicit result. The main weakness is that the numerical illustrations are used to support quantitative claims despite large model errors and without verifying the positive definiteness of the compositional Hessian for the effective ternary model. These issues are fixable within the scope of the manuscript: the author can add a stability check for the fitted models, compute k from the Gibbs seawater library for the seawater case, and explicitly frame the other examples as model-based illustrations with uncertainty estimates. If the author does not add such support, the paper should be weakened to a purely theoretical contribution. I recommend major revision rather than rejection because the central derivation is sound and the limitations are acknowledged in the conclusion."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this paper defines a coefficient k = β(T, μt) − β(T, y) that measures composition-driven volume change as a compressibility, shows it arises naturally when you hold chemical potentials rather than composition fixed, and gives an explicit ideal-mixture formula. The theory is clean and the core identity is new as far as the cited literature goes. The numerics are weaker and honestly labeled as illustrative.\n\nWhat's new: the decomposition β(T, μt) = β(T, y) + k, derived from the Gibbs function, and the interpretation that for incompressible mixtures (β(T,y)→0) the density can still vary through composition. That is a real conceptual clarification. The ideal-mixture formula (64) is explicit and useful. The derivation in Section 3.2 is direct and does not assume the target result. The paper also distinguishes two incompressible limits, which is worth taking seriously.\n\nSoft spots: all quantitative claims about real solutions rest on fitted ideal models, not full Gibbs free energies. The seawater ideal model gives 27.9% error in β(T,y); the ethanol binary model 35% error; the ternary ethanol model is better but is never checked for thermodynamic stability. The stress-test note is fair: the central identity requires D²_yy g positive definite, and for the water-ethanol ternary model of Appendix C the effective Hessian is never computed or checked. If an eigenvalue approaches zero, k diverges and the plotted ratios in Figs. 5–6 become artifacts. The paper itself acknowledges in the conclusion that wider data sets and full Gibbs functions are needed. So the theoretical framework holds up; the numbers for real mixtures are model-dependent. Also, no error bars are given, which is a fair criticism for a paper reporting quantitative ratios.\n\nThe self-citations for the thermodynamic framework are heavy, but the derivations here are self-contained enough, and the framework is standard in this niche. Not a flaw.\n\nWho it's for: people working on multicomponent incompressible flow models, or on thermodynamic closures for mixtures. It deserves a serious referee. My recommendation: send it to peer review. The core definition is likely to be reused, and the numerical section can be tightened with a stability check for the ternary model and an explicit statement of the fitted models' domain of validity.","headline":"A clean derivation of a new composition-compressibility coefficient, with honest but model-dependent numerics that need stability checks.","tokens_in":4,"tokens_out":3196,"would_cite":true,"duration_ms":50190,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["76T30","76N99","80A17","35Q35"],"pacs":[],"model":"deepseek-v4-flash","headline":"For a multicomponent fluid, zero isothermal compressibility does not force a constant density: the paper defines $k=\\beta(T,\\mu_t)-\\beta(T,y)$, built from derivatives of the Gibbs function, which measures composition-driven volume change…","keywords":["multicomponent fluid","compressibility","Gibbs function","incompressible limit","constant density approximation","volume of mixing","chemical potential","aqueous solutions"],"falsifier":"Take a binary mixture such as sucrose-water, measure density and isothermal compressibility over a range of temperature, pressure, and composition, construct the Gibbs function from these data, then compare $\\beta(T,\\mu_t)-\\beta(T,y)$ from eq. (55) with the measured difference between compressibility at fixed relative chemical potential and at fixed composition; a systematic mismatch would show the claimed identity does not hold.","tokens_in":18219,"feed_emoji":"🧪","tokens_out":11302,"duration_ms":88961,"temperature":0.7,"pith_summary":"For a single-component fluid, zero isothermal compressibility forces a constant density. The paper shows that the same conclusion does not hold for mixtures: even when $\\beta(T,y)=0$, the specific volume can still change because composition changes. To make that effect quantitative, the paper defines a coefficient $k=\\beta(T,\\mu_t)-\\beta(T,y)$, expressed through derivatives of the Gibbs function, which treats volume change on mixing as a true compressibility. The ratio $k/\\beta(T,y)$ then gives a criterion for when the constant-density approximation is justified. Three aqueous solutions, salt water, sucrose water, and ethanol-water, illustrate the sizes of this ratio in practice.","feed_headline":"Zero compressibility no longer forces constant density in mixtures","feed_subtitle":"The new coefficient ratio tells when a mixture can be treated as constant density, and when not.","key_machinery":"The load-bearing object is the coefficient $k=\\beta(T,\\mu_t)-\\beta(T,y)$, constructed from the Gibbs function $g(T,p,y)$ and its tangential compositional Hessian $\\mathrm{D}^2_{yy}g$. Explicitly, $k=\\partial_p g\\,(\\mathrm{D}^2_{yy}g)^{-1}(\\partial_y v/v)\\cdot(\\partial_y v/v)$, with derivatives taken along the hypersurface $\\sum_i y_i=1$. Because $\\beta(T,\\mu_t)$ is the pressure derivative of volume at fixed relative chemical potentials, $k$ has the same dimension and meaning as a compressibility; the inverse Hessian converts the composition response into this pressure-like measure. The derivation also uses the entropic-variable transformation of the full mixture equations to identify $\\beta(T,\\mu_t)$, and the ideal-mixture limit yields the closed form $k=\\frac{1}{RT}\\sum_i M_i\\rho_i(1/\\varrho-v_i)^2$, which makes the vanishing cases, a pure fluid or equal specific volumes, explicit.","core_discovery":"The central claim is that composition-induced volume changes in a multicomponent fluid can and should be quantified as a compressibility. If $v(T,p,y)$ is the specific volume and $\\beta(T,y)=-(1/v)\\partial_p v$ is the isothermal compressibility at fixed composition, the paper proves that the compressibility measured at fixed relative chemical potentials is $\\beta(T,\\mu_t)=\\beta(T,y)+k$, where $k$ is defined from the Gibbs function by $k=\\partial_p g\\,(\\mathrm{D}^2_{yy}g)^{-1}(\\partial_y v/v)\\cdot(\\partial_y v/v)$. Thus $k$ is the compressibility of mixing. In the incompressible limit $\\beta(T,y)\\to 0$ the coefficient $k$ need not vanish, so an incompressible mixture can still have a density that varies with composition; the density is constant only when $k$ also vanishes. For ideal mixtures the paper derives the explicit formula $\\beta(T,\\mu_t)-\\beta(T,y)=\\frac{1}{RT}\\sum_i M_i\\rho_i(1/\\varrho-v_i)^2$, which vanishes only for a pure fluid or when all species have the same specific volume. The ratio $k/\\beta(T,y)$ is proposed as the quantitative criterion for the validity of the constant-density approximation.","pith_inferences":["[Editorial inference] The same coefficient $k$ could be extracted from tabulated mixture-density data by finite differencing, even without a full Gibbs function, which would make the constant-density test routinely available for engineering mixtures.","[Editorial inference] Because the appendix shows the compressive part of the pressure gradient is controlled by $\\beta(T,\\mu_t)$ rather than $\\beta(T,y)$, sound-speed measurements in a mixture could provide a direct experimental route to $k$.","[Editorial inference] In a ternary reactive model such as the ethanol-water one, $k$ absorbs the volume-of-mixing defect into an effective species, so $k$ could serve as a general measure of non-ideal volume behaviour in mixtures with chemical association.","[Editorial inference] Since $k$ involves the inverse compositional Hessian, it should diverge as a mixture approaches a spinodal, potentially making $k$ an early-warning indicator of phase instability; the paper does not pursue this."],"forward_implications":["In a multicomponent fluid with $\\beta(T,y)=0$, the density can still vary with composition through $k$; only $k=0$ gives a truly constant density.","The ratio $k/\\beta(T,y)$ is a quantitative test for the constant-density approximation: when it is large, composition effects dominate and the approximation fails.","For ideal mixtures, $k$ is nonnegative and vanishes only for a pure fluid or when all species have identical specific volumes, so mixing alone always contributes a positive compressibility-like term.","The incompressible limit splits into two distinct regimes: $\\beta(T,y)\\to 0$ with the constraint $\\sum_i \\bar v_i\\rho_i=1$ (variable density), versus $\\beta(T,\\mu_t)\\to 0$ with $\\varrho=$ constant.","In the aqueous examples, the ratio exceeds about 6 for concentrated salt water and 10 for sucrose water, while for ethanol-water it stays near order one, suggesting the constant-density approximation is unsafe for the first two but may be acceptable for the third."],"supporting_citations":[{"why":"Defines incompressibility for multicomponent fluids and proves that zero isothermal compressibility still allows composition-dependent volume.","marker":"[4]"},{"why":"Supplies the entropic-variable transformation and the ideal-mixture Hessian formula from which the explicit form of k is derived.","marker":"[9]"},{"why":"Provides the thermodynamic framework, stability conditions, and the incompressibility definition used in Section 2.","marker":"[3]"},{"why":"Supplies the generalized incompressibility constraint sum of bar-v_i rho_i = 1 used to interpret the beta(T,y) tending to 0 limit.","marker":"[5]"},{"why":"Provides the incompressible-limit analysis separating the beta(T,y) tending to 0 and beta(T,mu_t) tending to 0 regimes.","marker":"[10]"},{"why":"Gives the association model for alcohol-water mixtures used to build the ternary ideal model in the ethanol illustration.","marker":"[21]"},{"why":"Supplies measured specific volumes of ethanol-water mixtures under pressure, used to compute isothermal compressibility for comparison.","marker":"[23]"},{"why":"Supplies tabulated densities for water-sucrose and water-ethanol at ambient conditions, used in the numerical illustrations.","marker":"[15]"}],"fun_headline_variants":["Mixtures can compress even when incompressible","Zero compressibility doesn't fix mixture density","New coefficient quantifies mixing-induced compression","Incompressible mixtures may still vary in density","Volume change from composition becomes a compressibility"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derivation assumes the mixture is thermodynamically stable, so the curvature of the Gibbs function with respect to composition is positive and invertible; near a spinodal or unstable state this inverse fails, $k$ diverges, and the split into $\\beta(T,y)$ and $k$ no longer describes a physical compressibility.","fun_headline_variants_meta":{"raw":{"variants":["Mixtures can compress even when incompressible","Zero compressibility doesn't fix mixture density","New coefficient quantifies mixing-induced compression","Incompressible mixtures may still vary in density","Volume change from composition becomes a compressibility"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000121,"raw_usage":{"total_tokens":1075,"prompt_tokens":907,"completion_tokens":168,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":523,"completion_tokens_details":{"reasoning_tokens":115}},"tokens_in":523,"tokens_out":168,"duration_ms":2591,"temperature":1.0,"reasoning_tokens":115,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T13:41:42.014139+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a binary mixture such as sucrose-water, measure density and isothermal compressibility over a range of temperature, pressure, and composition, construct the Gibbs function from these data, then compare $\\beta(T,\\mu_t)-\\beta(T,y)$ from eq. (55) with the measured difference between compressibility at fixed relative chemical potential and at fixed composition; a systematic mismatch would show the claimed identity does not hold.","supporting_citations":[{"cited_title":"Bothe, W","cited_arxiv_id":null,"evidence_quote":"Defines incompressibility for multicomponent fluids and proves that zero isothermal compressibility still allows composition-dependent volume."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the entropic-variable transformation and the ideal-mixture Hessian formula from which the explicit form of k is derived."},{"cited_title":"Bothe and W","cited_arxiv_id":null,"evidence_quote":"Provides the thermodynamic framework, stability conditions, and the incompressibility definition used in Section 2."},{"cited_title":"Bothe and P .-E","cited_arxiv_id":null,"evidence_quote":"Supplies the generalized incompressibility constraint sum of bar-v_i rho_i = 1 used to interpret the beta(T,y) tending to 0 limit."},{"cited_title":"Roux and J","cited_arxiv_id":null,"evidence_quote":"Gives the association model for alcohol-water mixtures used to build the ternary ideal model in the ethanol illustration."},{"cited_title":"Tanaka, T","cited_arxiv_id":null,"evidence_quote":"Supplies measured specific volumes of ethanol-water mixtures under pressure, used to compute isothermal compressibility for comparison."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies tabulated densities for water-sucrose and water-ethanol at ambient conditions, used in the numerical illustrations."}],"review_version":1}