{"id":"afee02ab-dccd-49a5-a5aa-365bf48308ee","arxiv_id":"2506.13483","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A modified extended Enskog theory with a Fermi-like free-volume multiplier and two adjustable parameters reproduces simulated self-diffusion coefficients of hard-sphere fluids in hard-sphere porous media.","lead":"The authors modify Enskog theory to describe how colloidal hard spheres diffuse through a disordered matrix of hard-sphere obstacles. They add a Fermi-like, free-volume factor with two adjustable parameters to improve agreement with 2004 computer simulations of self-diffusion.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equation (3.1) cannot reproduce the claimed agreement at η0=0.2: it reduces the common multiplier by a factor of ~2.6 (α=1) or ~4.9 (α=3), shifting D upward, contrary to Fig. 2 and the text.","rationale":"The reader's weakest assumption was the ad hoc, fitted nature of the Fermi-like multiplier and the absence of out-of-sample tests. That is a real concern, and the authors partly disclose it: Section 3 says 'We consider this parameter as the fitting parameter', and the Conclusions call both α and η*0 adjustable. My stress-test pass, however, found a more elementary problem that precedes the extrapolation question. As printed, Eq. (3.1) cannot be the formula behind Fig. 2: at η0=0.2 it replaces the multiplier φ0/φ ≈ 9.64 by 3.76 (α=1) or 1.97 (α=3), which shifts D1 upward by the inverse factors, while the text says the original EET already agreed at η0=0.2. One of the two claims is wrong, or the figure was computed from a formula not in the manuscript. This is an internal consistency failure in the central construction, not a disagreement with prevailing opinion. It also means the conditional acceptance proposed by the reader is too weak for the current text: the in-sample agreement itself must be reproducible from Eqs. (3.1)–(3.5) before fitting and generalization can be assessed. I recommend REJECT for the manuscript as submitted; if the numerical check shows a typographical error and the corrected curves match the simulation data, a revised manuscript could be reconsidered. No judgment about author intent is implied; the discrepancy is directly checkable from the printed equations.","tokens_in":9515,"tokens_out":17522,"duration_ms":165338,"concrete_test":"Recompute the NEET self-diffusion curves with the printed equations: for τ=1, η0=0.2, α=1 and α=3, η*0=0.4, evaluate M from Eq. (3.1), G10 and G11 from Eqs. (3.4)–(3.5), and D1 from Eq. (2.11). Compare with the original EET result using the same bracket but the multiplier φ0/φ from Eq. (2.6). The printed formula forces D_NEET = (φ0/φ)/M × D_EET at η0=0.2, i.e. a uniform upward shift of factor ≈2.6 (α=1) or ≈4.9 (α=3) at all η1. If the reproduced figure does not show this shift, the published figure was made with an unstated expression; if it does show the shift, the claimed agreement with simulation [12] at η0=0.2 is false. Either result settles whether the central claim is reproducible.","verdict_should_be":"REJECT","load_bearing_attack":"At τ=1 and η0=0.2, Eq. (2.6) gives φ0/φ ≈ 9.64 (φ/φ0 ≈ 0.104). In Eq. (3.1), for η0 < η*0 = 0.4 the exponential is less than 1, so the denominator 1 + (φ/φ0 − 1) exp[α(η0 − η*0)] is larger than φ/φ0. The new multiplier is therefore smaller than φ0/φ: with α=1 it is ≈ 3.76, and with α=3 it is ≈ 1.97, instead of 9.64. Because both G10 and G11 in Eqs. (3.4)–(3.5) carry this common multiplier, D1 from Eq. (2.11) is larger by the inverse factor at every η1: about 2.6-fold for α=1 and 4.9-fold for α=3 at η0=0.2. The paper states that the original EET (with φ0/φ) already agrees with simulation at η0=0.2 (Fig. 1a). If that is true, the NEET curves in Fig. 2 cannot also agree at η0=0.2; if the figure was generated from a different expression than the printed one, the manuscript must say so. This is an internal inconsistency in the central construction, independent of whether the two parameters are fitted. The fitting issue raised by the reader remains valid as well: α is called a fitting parameter in Section 3, η*0 is set after inspecting Fig. 1a, and no out-of-sample test beyond τ=1 is reported.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a new extended Enskog theory (NEET) for the self-diffusion coefficient of a hard-sphere fluid in a disordered hard-sphere matrix. The key ingredient is the replacement of the porosity ratio φ0/φ, which multiplies the modified contact values G10 and G11, by a Fermi-like function of the matrix packing fraction η0 with two adjustable parameters α and η0* (Eq. 3.1). The authors report that with α=3 and η0*=0.4 the theoretical prediction is in the best agreement with the simulation data of Ref. [12] for all considered values of the fluid packing fraction η1 and matrix packing fraction η0 at size ratio τ=1.","tokens_in":9853,"tokens_out":10776,"duration_ms":90182,"significance":"If the central claim held, the paper would offer a simple semi-empirical extension of Enskog theory for diffusion in model porous media, and it would be a useful reference for follow-up work. The manuscript is clearly structured and connects to the authors' earlier developments in scaled particle theory and Enskog theory. However, the claimed agreement is not reproducible from the printed equations: at η0=0.2 the proposed multiplier is far smaller than φ0/φ, moving the theory in the wrong direction. The two-parameter fit to a single data set also limits the significance of the reported agreement.","major_comments":[{"comment":"At τ=1 and η0=0.2, Eq. (2.6) gives φ0/φ≈9.64. The multiplier in Eq. (3.1) for η0=0.2<η0*=0.4 is about 3.76 for α=1 and 1.97 for α=3, i.e., smaller than φ0/φ by factors of 2.6 and 4.9, respectively. Since this multiplier is common to G10 and G11 in Eqs. (3.4)–(3.5), the NEET self-diffusion coefficient from Eq. (2.11) is larger than the EET value by the inverse factors at all η1. The paper states that the original EET at η0=0.2 is in 'more or less correct agreement' with simulation (Section 3, discussion of Fig. 1a). The NEET curves in Fig. 2 therefore cannot also agree with simulation at η0=0.2. Either the figure was generated from a different expression than the one printed, or the text description is in error. This contradiction directly undermines the central claim of the paper.","section":"Section 3, Eqs. (3.1), (3.4), (3.5), (2.11)"},{"comment":"The two parameters α and η0* are not derived or independently motivated: α is called a 'fitting parameter' in Section 3, and η0* is fixed by inspecting Fig. 1a to align the theoretical curves with the simulation data of Ref. [12]. The paper reports no out-of-sample test of the resulting theory, for example on different size ratios τ, different matrix structures, or densities outside the fitted range. Thus the 'best agreement' is a demonstration of interpolation within the fitted data set, not a falsifiable prediction of the proposed theory.","section":"Section 3, Eq. (3.1)"}],"minor_comments":[{"comment":"The caption states that G10 and G11 are given by expressions (3.2) and (3.3), but those equations are limiting forms of the multiplier, not the full functions. The correct references should be Eqs. (3.4) and (3.5).","section":"Figure 2 caption"},{"comment":"The phrase 'slightly reestimates the value of D1' appears to be a typo; the intended verb is likely 'underestimates' or 'overestimates'. The direction of the deviation matters for the motivation of the modification, so it should be stated clearly.","section":"Section 3, first paragraph"},{"comment":"The notation 2πσ^2_{11} in the prefactor is unusual for a binary-mixture expression; please clarify whether σ11 denotes σ1 and make the notation consistent with the rest of the paper.","section":"Section 2, Eq. (2.3)"},{"comment":"There are minor language issues, such as 'Fermi-like for dependence on η0', that could be polished.","section":"Section 3, text near Eq. (3.1)"}],"recommendation":"reject","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: the paper is a transparent, compact extension of the authors' earlier EET, and the equations are explicit enough to reimplement. But the central formula as printed cannot produce the figures, so the main claim is not supported until that is sorted out.\n\nThe new thing is Eq. (3.1): replacing the porosity ratio phi0/phi in both G10 and G11 by a Fermi-like function with two parameters, alpha and eta*0. That is new relative to [10], which modified only G10. Credit where due: the authors say plainly that alpha is a fitting parameter and eta*0 is fixed after viewing Fig. 1a; they show alpha=1 and alpha=3; and the comparison set is the standard Chang/Yethiraj simulation data. The derivation is not first-principles — G10 uses the truncated expansion (2.8), G11 follows by \"similar manipulation\", and (3.1) is borrowed from Poisson-Fermi theory by analogy — but the disclosure is honest and the formulas are easy to test.\n\nThe soft spot is bigger than a fitted parameter. As written, (3.1) makes the new multiplier smaller than phi0/phi for every eta0 below eta*0, because phi/phi0 - 1 is negative. Numerically at tau=1, eta0=0.2, phi0/phi ~ 9.64; Eq. (3.1) gives ~3.8 for alpha=1 and ~2.0 for alpha=3. Since G10 and G11 both carry this common factor, D1 from (2.11) is raised by a factor of ~2.6 or ~4.9 relative to the original EET. But the text says the original EET already agrees with simulation at eta0=0.2 and overestimates D1 at eta0=0.05-0.15. A correction that increases D1 cannot improve those curves. So either (3.1) has a sign/notation error (phi and phi0 swapped, or a missing negative in the exponent), or Fig. 2 was made with different expressions than the printed ones. The manuscript needs a correction and a re-plot before the \"best agreement\" claim can be evaluated.\n\nThe reader's concern about fitting is valid but, by itself, would only justify CONDITIONAL. The equation/figure mismatch is the deciding issue. No out-of-sample test is reported: tau=1 only, and the same simulation data used to set eta*0 is the data used for validation. Citation practice is fine; self-citations are to the authors' own prior equations, which are necessary context.\n\nWho is this for? A narrow soft-matter readership. If the formula is fixed, it is a small useful improvement with an explicit ansatz. As it stands, I would not cite it. I would send it to a competent referee rather than desk-reject, because the inconsistency looks fixable and the paper is transparent; the referee should check Eq. (3.1) against Fig. 2 and the original EET curves. Major revision.","headline":"An honest but fitted Fermi-like ansatz, and as printed Eq. (3.1) moves D in the wrong direction for the matrix densities where the paper claims improvement.","tokens_in":10441,"tokens_out":8480,"would_cite":false,"duration_ms":79865,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A new extended Enskog theory with a Fermi-like free-volume factor reproduces the simulated self-diffusion of hard-sphere fluids in disordered hard-sphere matrices across all densities tested.","keywords":["hard-sphere fluid","disordered porous media","scaled particle theory","extended Enskog theory","self-diffusion coefficient","Fermi-like distribution","free volume","trapping"],"falsifier":"Take the same formula (2.11) with the Fermi-like substitution (3.1), $\\alpha=3$ and $\\eta_0^*=0.4$, and compare its predictions to new simulations with fluid particles smaller or larger than the matrix particles (for example size ratios 0.5 or 2), or with matrix packing fractions above 0.2; a systematic deviation at any fluid density would show the ansatz is fitted rather than general.","tokens_in":9229,"feed_emoji":"🔬","tokens_out":12562,"duration_ms":110208,"temperature":0.7,"pith_summary":"This paper tries to fix a known failure: Enskog theory, the standard kinetic description of dense hard-sphere fluids, systematically overestimates how fast a hard-sphere fluid diffuses when it is confined in a disordered matrix of hard-sphere obstacles. The authors' explanation is that the old theory treats the immobile matrix only as extra scatterers and misses the volume lost because fluid particles get trapped against matrix particles. Their new extended Enskog theory replaces the porosity ratio $\\phi_0/\\phi$ in the collision terms by a Fermi-like distribution with two adjustable parameters, $\\alpha$ and $\\eta_0^*$. With $\\alpha=3$ and $\\eta_0^*=0.4$, the predicted self-diffusion coefficient agrees with the computer-simulation data [12] for equal-sized particles at every fluid and matrix packing fraction considered. If the ansatz holds beyond the fitted case, the payoff is a cheap analytic formula for transport in porous media that captures trapping physics without expensive simulation.","feed_headline":"Fermi-like factor aligns hard-sphere diffusion theory with simulations","feed_subtitle":"The two-parameter factor makes the theory match simulated diffusion at every tested density.","key_machinery":"The load-bearing mechanism is the Fermi-like multiplier (3.1), $[1+(\\phi/\\phi_0-1)\\exp(\\alpha(\\eta_0-\\eta_0^*))]^{-1}$, inserted in place of $\\phi_0/\\phi$ in the definitions of $G_{10}(\\sigma_{10})$ and $G_{11}(\\sigma_{11})$ (equations (3.4) and (3.5)). Near the bulk limit $\\eta_0\\to 0$ the multiplier tends to unity, recovering the earlier extended Enskog theory; at $\\eta_0=\\eta_0^*$ it equals $\\phi_0/\\phi$, recovering the authors' previous modified contact values. The parameter $\\alpha$ controls how sharply the matrix packing fraction suppresses diffusion, while $\\eta_0^*$ is interpreted as the saturated matrix packing fraction beyond which percolation-driven changes would require a modified description. These generalized functions are not true pair contact values, which is what lets them absorb the free-volume effects of fluid trapping.","core_discovery":"The central claim is that the earlier extended Enskog treatment overestimates self-diffusion because it represents the matrix only as infinitely massive scatterers; it does not account for the volume made unavailable when mobile fluid particles are trapped by immobile matrix particles. The fix is to replace the contact values $g_{11}(\\sigma_{11})$ and $g_{10}(\\sigma_{10})$ by effective functions $G_{11}$ and $G_{10}$ built from scaled-particle-theory expressions, multiplied by the porosity ratio $\\phi_0/\\phi$, and then to replace that multiplier by the Fermi-like expression (3.1), $[1+(\\phi/\\phi_0-1)\\exp(\\alpha(\\eta_0-\\eta_0^*))]^{-1}$. With $\\alpha=3$ and $\\eta_0^*=0.4$, equation (2.11) reproduces the simulation data [12] for $\\tau=1$ across all studied fluid and matrix packing fractions, which the paper reports as the best agreement between theory and simulation. The new functions are explicitly not contact values of pair distribution functions; they are effective inputs carrying the free-volume and trapping physics.","pith_inferences":["An inference not drawn in the paper: if this Fermi ansatz is a true free-volume effect, the fitted parameters $\\alpha$ and $\\eta_0^*$ should vary smoothly with the size ratio $\\tau$; measuring them at $\\tau\\neq1$ would be a direct test.","The same Fermi-like multiplier could be applied to other Enskog-level transport coefficients of confined fluids, such as shear viscosity or thermal conductivity, which the paper does not compute, yielding a full transport description from the same two parameters.","The interpolation between the bulk limit and the saturated porosity limit has the same shape as saturation laws familiar from adsorption phenomena, which suggests that $\\alpha$ might be derivable from a thermodynamic argument rather than treated as adjustable; the paper does not attempt such a derivation.","A further consequence of the paper's own logic is that in very dense matrices the theory should break down near $\\eta_0^*$ because the trapped-fluid fraction itself becomes a collective, percolation-controlled quantity; simulations scanning $\\eta_0$ up to that threshold would map where the Fermi fit begins to fail."],"forward_implications":["With $\\alpha=3$ and $\\eta_0^*=0.4$, equation (2.11) matches the simulated self-diffusion of [12] for equal-sized hard spheres at all the matrix and fluid packing fractions shown.","At small matrix packing fraction the Fermi multiplier is essentially unity, so the new theory continuously reduces to the standard Enskog result for bulk hard spheres, and at $\\eta_0=\\eta_0^*$ it reduces to the authors' previous extended Enskog theory.","Because $\\alpha$ is the explicit control on the influence of matrix packing, comparing the theory with simulation at one matrix density effectively fixes the strength of matrix-induced slowing for other densities.","The separate effective functions $G_{10}$ and $G_{11}$ let the theory attribute deviations to fluid-matrix versus fluid-fluid trapping, and the paper states that above $\\eta_0^*$ percolation makes diffusion anomalous so the present description applies below that threshold."],"supporting_citations":[{"why":"Supplies the computer-simulation self-diffusion data for $\\tau=1$ that the Fermi-like theory is fitted to and compared against.","marker":"[12]"},{"why":"Provides the previous extended Enskog theory for hard spheres in porous media, including the modified fluid-matrix function that this paper generalizes.","marker":"[10]"},{"why":"Gives the scaled-particle-theory contact-value expressions for $g_{11}$ and $g_{10}$ that the new effective functions replace.","marker":"[9]"},{"why":"Establishes the SPT2 approach that fixes the SPT1 inconsistency and yields the porosity coefficients $\\phi$ and $\\phi^*$ used here.","marker":"[5]"},{"why":"Introduces the scaled-particle-theory description of hard spheres in porous media and the geometric porosity parameter $\\phi_0$.","marker":"[2]"},{"why":"Sets out the bulk Enskog theory to which the new formulas reduce at zero matrix packing.","marker":"[11]"},{"why":"Supplies the Poisson-Fermi analogy that motivates replacing the porosity ratio by a Fermi-like distribution.","marker":"[16]"},{"why":"Underlies the interpretation of $\\eta_0^*$ as a saturation packing fraction where percolation sharply changes transport.","marker":"[17]"}],"fun_headline_variants":["Fermi factor sharpens hard-sphere diffusion in porous media","Trapped volume rescues Enskog theory for porous fluids","New Enskog theory nails diffusion in disordered matrices","Free-volume factor aligns diffusion theory with simulations"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that one fixed Fermi-like factor, with $\\alpha=3$ and $\\eta_0^*=0.4$ chosen by hand after inspecting the simulation data, captures the trapping physics of hard-sphere fluids in any disordered hard-sphere matrix; the paper tests it only for equal-sized particles and matrix packing fractions up to 0.2.","fun_headline_variants_meta":{"raw":{"variants":["Fermi factor sharpens hard-sphere diffusion in porous media","Trapped volume rescues Enskog theory for porous fluids","New Enskog theory nails diffusion in disordered matrices","Free-volume factor aligns diffusion theory with simulations"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000154,"raw_usage":{"total_tokens":1191,"prompt_tokens":904,"completion_tokens":287,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":520,"completion_tokens_details":{"reasoning_tokens":221}},"tokens_in":520,"tokens_out":287,"duration_ms":3042,"temperature":1.0,"reasoning_tokens":221,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T20:00:01.640228+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the same formula (2.11) with the Fermi-like substitution (3.1), $\\alpha=3$ and $\\eta_0^*=0.4$, and compare its predictions to new simulations with fluid particles smaller or larger than the matrix particles (for example size ratios 0.5 or 2), or with matrix packing fractions above 0.2; a systematic deviation at any fluid density would show the ansatz is fitted rather than general.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the computer-simulation self-diffusion data for $\\tau=1$ that the Fermi-like theory is fitted to and compared against."},{"cited_title":"F., Korvatska M","cited_arxiv_id":null,"evidence_quote":"Provides the previous extended Enskog theory for hard spheres in porous media, including the modified fluid-matrix function that this paper generalizes."},{"cited_title":"V., Holovko M., Patsahan T., Cummings P","cited_arxiv_id":null,"evidence_quote":"Gives the scaled-particle-theory contact-value expressions for $g_{11}$ and $g_{10}$ that the new effective functions replace."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the SPT2 approach that fixes the SPT1 inconsistency and yields the porosity coefficients $\\phi$ and $\\phi^*$ used here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the scaled-particle-theory description of hard spheres in porous media and the geometric porosity parameter $\\phi_0$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Sets out the bulk Enskog theory to which the new formulas reduce at zero matrix packing."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Poisson-Fermi analogy that motivates replacing the porosity ratio by a Fermi-like distribution."},{"cited_title":"(Eds.), Fractals and Disordered Systems, Springer, Berlin, Heidelberg, 1991, doi:10.1007/978-3-642-51435-7","cited_arxiv_id":null,"evidence_quote":"Underlies the interpretation of $\\eta_0^*$ as a saturation packing fraction where percolation sharply changes transport."}],"review_version":2}