{"id":"be4caf15-e1af-43c3-a40a-74ce0f9ef96d","arxiv_id":"2501.05117","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Treating lattice site count as an extra thermodynamic variable yields a generalized Gibbs-Duhem equation and new response functions for vacancy formation in compressible crystals.","lead":"This paper derives a generalized Gibbs-Duhem equation for compressible crystal lattices by treating the number of lattice sites as an independent thermodynamic variable. It uses this framework to compare how crystals respond to pressure and chemical potential changes at constant volume versus constant pressure.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The generalized Gibbs-Duhem formalism is internally consistent, but the model's pressure decomposition (Eq. 80) is admitted to lack microscopic justification, leaving the quantitative adsorption predictions unsubstantiated.","rationale":"The paper's formal core—the extended Euler equation and the generalized Gibbs-Duhem relation (Eqs. 5, 11)—is internally consistent and follows from the model's manifest extensivity. The reader's weakest assumption (homogeneous variation of M) is a standard postulate of Larché-Cahn theory; in the thermodynamic limit, M scales with the crystal size and the surface mechanism of site addition contributes negligibly. I do not see a genuine flaw there. The truly load-bearing weakness is the constitutive decomposition of pressure into a molecular component p_c(N/V) and an elastic component σ_j(V/M) (Eqs. 80-87). The paper explicitly states (Section I, Section IX) that it cannot supply a statistical-mechanical justification for this decomposition. Every quantitative prediction—the lattice site potential (Eq. 74), the open-system bulk modulus (Eq. 141), the adsorption susceptibilities (Eqs. 149, 156, 160)—depends on this separation. If in a real crystal the hydrostatic pressure is not separable in this way, the advertised quantification of vacancy creation under (μ,p,T) conditions is unsupported, even though the formal Gibbs-Duhem identity remains true. The paper is therefore best read as a formal demonstration plus a toy model; the condition for acceptance is an independent microscopic justification of the pressure decomposition or a substantial caveat that the results are constitutive-model-dependent. This matches the reader's CONDITIONAL verdict, so no change is needed.","tokens_in":114,"tokens_out":18525,"duration_ms":304836,"concrete_test":"Run canonical Monte Carlo simulations of a Lennard-Jones (or hard-sphere) crystal at fixed number of sites M over a grid of (N,V) values, computing the hydrostatic pressure from the virial. Fit the scaled pressure p v_c^R to the functional form γ_b(c/J)^2 - α_j(J-1) (Eq. 87). If the fit requires cross-terms coupling c and J that cannot be absorbed into a single density-dependent p_c and a single stretch-dependent elastic stress, the pressure decomposition underlying all quantitative results fails. A simpler analytical check: re-derive the pressure from the partition function of a one-component lattice gas with binding energy I(ρ) without postulating the additive decomposition; if cross-derivatives ∂^2 p/∂c∂J are nonzero, the decomposition is not generally valid.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central formal claim—that including M as an extensive variable yields the Euler equation (5) and the generalized Gibbs-Duhem relation (11)—is internally consistent for the model, given the Gibbs equation (1). The load-bearing weakness is the constitutive decomposition of pressure into a molecular part p_c(N/V) and an elastic part σ_j(V/M) (Eq. 80). The paper itself states in Section I and again in Section IX that it cannot provide a statistical-mechanical justification ('requires justification in terms of structural correlations. We are not yet able to provide such a microscopic underpinning'). All model-specific predictions—the lattice site potential (Eq. 74), the response functions (Eqs. 149, 156, 160), and the claimed quantification of vacancy creation—depend on this separation. If the true pressure of a compressible lattice gas does not decompose this way, the quantitative content of the paper is unsubstantiated, even though the formal Gibbs-Duhem identities (Eqs. 12-13) remain true. The paper's own limitation statement thus undermines the advertised adsorption application.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops an extended thermodynamics for a one-component compressible lattice-gas crystal by promoting the number of lattice sites M to an independent extensive variable. It derives the Euler relation U = TS - pV + mu N + nu M (Eq. 5), identifies the grand free enthalpy E = G - mu N = nu M (Eq. 6), and obtains a generalized Gibbs-Duhem relation M dnu = V dp - N dmu (Eq. 11) together with the adsorption relations Eqs. (12)-(13). After a rigid-lattice warm-up, it introduces a compressible-lattice free energy (Eq. 56) with a molecular binding term depending on N/V and an elastic term depending on V/M, and derives the chemical potential, lattice-site potential, pressure, and a set of response functions for the (mu,V,T) and (mu,p,T) ensembles. The formal identities are internally consistent; the quantitative model predictions are conditional on a constitutive pressure decomposition that the author explicitly states lacks microscopic justification.","tokens_in":26494,"tokens_out":13066,"duration_ms":121574,"significance":"If the constitutive assumption could be justified, or if the paper is read as a deliberately minimal model, the work would be a useful contribution: it gives a concrete thermodynamic meaning to the lattice-site potential nu, shows how the generalized Gibbs-Duhem relation yields an adsorption route to occupation and cell volume, and clarifies why the usual liquid-state Gibbs-Duhem relation fails for crystalline solids. The formal part is derived analytically from a stated free energy with no parameter fitting, and the author is unusually candid about the missing microscopic underpinning and about the unclosed nature of the response functions. The main value is therefore in the thermodynamic framework; the quantitative vacancy-creation predictions are not yet established.","major_comments":[{"comment":"The quantitative content of the model, including Eqs. (74), (149), (156), and (160), rests on the additive decomposition of pressure into a molecular component p_c^*(rho) and an elastic component sigma_j(V/M) stated in Eqs. (80)-(87), together with the linearized binding energy in Eq. (82). The manuscript itself declares in Section I and again in Section IX that it cannot yet provide a statistical-mechanical justification for this decomposition, which it says 'requires justification in terms of structural correlations.' Because this split is the only mechanism coupling occupation to strain, all adsorption and vacancy-creation predictions are conditional on an unvalidated constitutive assumption. The formal generalized Gibbs-Duhem identities (Eqs. (11)-(13)) do not depend on this assumption and are on solid ground, but the advertised quantitative application is not. The authors should either supply a microscopic derivation or explicitly re-scope the paper as a formal theory illustrated by a toy model.","section":"Section IV, Eqs. (56), (80)-(87); Section IX"},{"comment":"The Euler integration leading to Eq. (5) treats M as a homogeneous extensive variable, while Section III B correctly notes that lattice sites 'are not inserted but added at the periphery.' These statements can be reconciled in the thermodynamic limit, where surface contributions are negligible, but the paper never says this explicitly. It also does not address how changes in M are realized in the bulk, for example through dislocation climb or grain-boundary sources, within a homogeneous equilibrium model. Please add an explicit statement of the intended domain of validity: for macroscopic crystals with negligible surface-to-volume ratio, peripheral addition is thermodynamically equivalent to homogeneous scaling of M; for finite systems Eq. (5) is an approximation.","section":"Section II A; Section III B"}],"minor_comments":[{"comment":"In Eq. (129), the coefficient of dc should be k_B T/(h c) + Gamma_b J/c, not k_B T/(h c) + B_c. With B_c defined as Gamma_b c/J in Eq. (123), the printed equation does not follow from Eq. (128). The final result Eq. (130) is nevertheless correct if the coefficient is replaced by Gamma_b J/c.","section":"Section VII C, Eq. (129)"},{"comment":"The sentence claiming that the relation kappa_mu = 1/B_mu is satisfied by Eqs. (160) and (141) is not correct as printed: Eq. (160) is a small-h approximation, and the product of Eqs. (160) and (141) differs from unity for typical parameter values. The reciprocal relation holds for the exact kappa_mu obtained from Eqs. (144) and (149). Please correct Eq. (160) or qualify the statement.","section":"Section VII D, after Eq. (160)"},{"comment":"The phrase 'For attractive interactions' is left dangling; no statement is made about the sign of rho_0/rho_R for gamma_b < 0. Please complete the sentence.","section":"Section VI D"},{"comment":"The heading uses 'Absorption' but the text and standard terminology is 'adsorption'; Eq. (38) is an adsorption isotherm, not an absorption isotherm.","section":"Section III C"},{"comment":"Please correct typographical errors: 'expnsion' in Section I and Section VII, 'pessure' in Section VI D, 'vlaue' in Section VIII C, 'Uderstanding' in Ref. 57, and 'Spinger' in Ref. 37.","section":"Throughout"},{"comment":"The model name 'Blume-Emmery-Grifiths' should be 'Blume-Emery-Griffiths'.","section":"Section VIII A"}],"recommendation":"major_revision","confidential_remarks":"This is a single-author, programmatic theoretical paper. Its main obstruction to acceptance is the admitted lack of microscopic justification for the pressure decomposition; if the editor believes that a formally correct but intentionally minimal model is publishable, the paper could be accepted after the authors re-scope the claims. The relation to the author's previous work (Refs. 13 and 14) is openly acknowledged, and I see no citation-pattern concern. The paper is somewhat over-long for its content, with several sections of qualitative discussion that could be shortened."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know up front. First, the core claim—that treating lattice-site number M as an additional extensive variable restores the Euler equation and a generalized Gibbs-Duhem relation, with ν obeying the adsorption relations (∂ν/∂μ)_{p,T} = -c and (∂ν/∂p)_{μ,T} = v_c—is internally consistent and correctly worked out. The formal identities are not in question. Second, every quantitative prediction of the model, including the response functions and the vacancy-creation analysis, rests on the pressure decomposition p⋆ = p⋆_c − σ_j, with the molecular part a function of N/V and the elastic part a function of V/M. The author states plainly in Sections I and IX that he cannot yet provide a microscopic justification for this split. That is an explicit limitation, not a hidden flaw. But it does mean the abstract's promise to 'quantify' vacancy creation is stronger than what the paper actually delivers—the equation of state is never determined, as Section IX.B admits.\n\nThe genuinely new content is the specific compressible lattice-gas model and the systematic comparison of (μ,p,T) and (μ,V,T) response functions. The M/ν formalism itself is present in Larché-Cahn and the Voorhees–Johnson review; the paper says so and builds on it. The relation ξ_p = κ_N ξ_V (Eq. 153) is new and is the kind of concrete, checkable result that could be useful in simulation design.\n\nSoft spots, in order of importance. (1) The Euler integration (Eq. 5) treats M as a variable that can be scaled homogeneously. In a real crystal, sites are only added or removed at surfaces, grain boundaries, or dislocations—the paper makes this point itself in Section III.B. So the extensivity and generalized Gibbs-Duhem relation are postulates of the LC network picture, not consequences of ordinary bulk thermodynamics. I don't think this is fatal; it is the standard LC assumption and the paper is explicit. But a referee should not let it pass without discussion. (2) The pressure decomposition is load-bearing and unvalidated. This qualifies the quantitative content severely, and the paper deserves credit for saying so. (3) Two minor issues: Eq. (129) has a coefficient typo (B_c should be Γ_b J/c, not Γ_b c/J; the following Eq. (130) is correct), and the assertion that κ_μ = 1/B_μ is made without proof. The author should supply the proof or a citation.\n\nThe paper is clear, honest, and technically careful. It is written for readers who work on solid thermodynamics, Larché-Cahn theory, or open-system simulation ensembles, and it deserves a serious referee. I would send it to review with a request for minor revisions only: fix the typo, prove or reference the κ_μ–B_μ identity, and tone down the abstract's 'quantify' language to match the admitted limitations.","headline":"A careful, honest theory paper that gets the formal Larché-Cahn thermodynamics right and is upfront about the one constitutive assumption its quantitative conclusions depend on.","tokens_in":26959,"tokens_out":3146,"would_cite":true,"duration_ms":31370,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["82B20","80A10","82D25"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper claims that the extensivity failure of Larché–Cahn crystals is cured by promoting the number of lattice sites $M$ to an independent extensive variable, which introduces a conjugate force $\\nu$ and restores a generalized…","keywords":["compressible lattice gas","Larché-Cahn theory","lattice site potential","generalized Gibbs-Duhem equation","adsorption isotherm","vacancy formation","extensivity","chemical potential of solids"],"falsifier":"Simulate a finite crystal at fixed $M$ under $(\\mu,p,T)$ control and measure $E=G-\\mu N$ by thermodynamic integration; the claim fails if $E/M$ is not an intensive function of occupation and stretch alone, or if the value of $\\nu$ inferred from $(\\partial\\nu/\\partial\\mu)_{p,T}=-c$ disagrees with the directly computed free-energy derivative with respect to $M$.","tokens_in":26012,"feed_emoji":"🧊","tokens_out":6084,"duration_ms":51940,"temperature":0.7,"pith_summary":"Larché–Cahn theory holds that the number of lattice sites in a crystal is conserved under deformation, so composition and strain couple through two distinct densities: occupation $N/M$ and cell volume $V/M$. The paper argues that because of this, the standard Gibbs–Duhem relation of liquids fails for crystals, and extensivity is restored only by treating $M$ as an additional extensive variable. The new conjugate force $\\nu$, called the lattice site potential, makes the Euler equation read $U = TS - pV + \\mu N + \\nu M$ and turns the Gibbs–Duhem relation into an adsorption equation whose derivatives are the occupation and the deformed cell volume. A nonzero $\\nu$ is what thermodynamically separates a solid from a liquid under open isobaric–isothermal conditions. Working out a one-component compressible lattice gas, the paper derives the $(\\mu,p,T)$ response functions and shows that their deviation from $(\\mu,V,T)$ behavior is controlled by the closed-system compressibility.","feed_headline":"One extra variable restores crystal thermodynamics","feed_subtitle":"Treating the number of lattice sites as extensive yields a new force ν and an adsorption equation for solids.","key_machinery":"The central object is the lattice site potential $\\nu$, defined as the thermodynamic force conjugate to the conserved number of lattice sites $M$ in the generalized Gibbs equation $dU = T\\,dS - p\\,dV + \\mu\\,dN + \\nu\\,dM$. It carries the argument: with $M$ included in the system-size scaling, the Euler–Gibbs integration is admissible, giving $U = TS - pV + \\mu N + \\nu M$; differentiating $E = \\nu M$ produces the generalized Gibbs–Duhem equation, and converting that equation into derivatives of $\\nu$ with respect to $\\mu$ and $p$ yields the adsorption isotherms. In the model, $\\nu$ is evaluated from the free energy per site $f_e(c,J) = f_s(c) - cI(\\rho) + g_j(J)$, which gives $\\nu = k_B T \\ln(1-c) - \\tfrac{1}{2}(J+1)\\sigma_j$, coupling vacancy population to elastic stress.","core_discovery":"The central claim is that extensivity of a Larché–Cahn crystal is restored by treating the number of lattice sites $M$ as an additional extensive variable, so that the fundamental relation integrates to $U = TS - pV + \\mu N + \\nu M$ (Eq. 5), the grand Gibbs free energy satisfies $E = G - \\mu N = \\nu M$ (Eq. 6), and the Gibbs–Duhem relation becomes $M\\,d\\nu = V\\,dp - N\\,d\\mu$ (Eq. 11). This yields the adsorption isotherms $(\\partial\\nu/\\partial\\mu)_{p,T}=-c$ and $(\\partial\\nu/\\partial p)_{\\mu,T}=v_c$ (Eqs. 12–13). The paper works this out for a uniform one-component compressible lattice gas in which pressure is decomposed into a molecular component depending on deformed density $N/V$ and an elastic term linear in the volume strain set by $V/M$; in the harmonic approximation the lattice site potential is $\\nu = k_B T \\ln(1-c) - \\tfrac{1}{2}(J+1)\\sigma_j$. The main conclusion is that $\\nu$ is generally nonzero, so crystals under open constant-pressure, constant-temperature conditions do not obey the liquid relations $G=\\mu N$ and $\\Omega=-pV$, and the difference is a measurable solid-state quantity.","pith_inferences":["The same $M$-scaling logic could plausibly extend to shear deformations, where the Gibbs-prism paradox blocks a unique chemical potential; a tensorial analogue of $\\nu$ might resolve that ambiguity for non-hydrostatic stress states.","If $\\nu$ is accessible through thermodynamic integration in experiments, mapping it across the $(\\mu,p)$ plane would yield a new solid-state equation of state that supplements the usual pressure–density data.","The paper's pressure decomposition suggests a microscopic test: in a lattice-fixed ensemble with $M$ held constant, the molecular pressure $p_c$ should equal the virial of the deformed-density interactions while the elastic stress $\\sigma_j$ is the derivative of the reference-cell energy; a simulation that measures these pieces separately would check the constitutive split.","The formal analogy to Hill's nanothermodynamics implies that finite crystals should show replica-energy-like corrections scaling as $1/M$, which could be tested by comparing thermodynamic functions of small and large clusters."],"forward_implications":["If $\\nu$ is nonzero, a crystal in equilibrium with reservoirs of fixed $\\mu$ and $p$ does not satisfy the liquid identities $G=\\mu N$ and $\\Omega=-pV$; the mismatch $E = G - \\mu N = \\nu M$ becomes a measurable thermodynamic property of the solid.","The generalized Gibbs–Duhem relation provides integration paths: measuring occupation $c$ and deformed cell volume $v_c$ over the $(\\mu,p)$ plane determines relative values of the lattice site potential $\\nu$.","The adsorption form of the Gibbs–Duhem relation quantifies vacancy creation under isothermal–isobaric conditions through the derivatives $(\\partial\\nu/\\partial\\mu)_{p,T}=-c$ and $(\\partial\\nu/\\partial p)_{\\mu,T}=v_c$.","Open-system response functions in $(\\mu,p,T)$ differ from their $(\\mu,V,T)$ counterparts, and the difference is controlled by the closed-system compressibility; the two ensembles become equivalent only in the limit of an ideal crystal with no vacancies.","Accretion—enlarging the crystal by scaling $M$, $N$, and $V$ together—is thermodynamically consistent only if the $\\nu$ term is included in the energy balance."],"supporting_citations":[{"why":"Supplies the network principle that the number of lattice sites is conserved under deformation, the foundational postulate from which $M$ becomes an extensive variable.","marker":"[3]"},{"why":"Establishes the Larché–Cahn theory of composition–stress coupling and the diffusion-potential framework that the paper extends.","marker":"[2]"},{"why":"Provides the review context for the lattice-site potential $\\nu$ and its role in elastically stressed crystals.","marker":"[9]"},{"why":"Shows that surface thermodynamics of crystals generates nonzero $\\nu$ through the difference between interface energy and surface stress, supporting the claim that $\\nu$ is generally nonvanishing.","marker":"[5]"},{"why":"Serves as the formal template for the extended thermodynamics; the replica energy in Hill's nanothermodynamics is analogous to $\\nu M$.","marker":"[26]"},{"why":"Provides the contrasting simulation scheme that treats $M$ as an unconstrained internal degree of freedom with $\\nu=0$, which the paper argues differs from a Larché–Cahn crystal.","marker":"[15]"}],"fun_headline_variants":["Lattice sites: the missing variable in crystal thermodynamics","Solids break G=μN and Ω=−pV—new force ν appears","Adsorption equation for crystals from lattice-site force","Counting lattice sites fixes crystal thermodynamics"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the number of lattice sites $M$ can be scaled like an ordinary extensive variable, so that the fundamental relation integrates to $U = TS - pV + \\mu N + \\nu M$; in real crystals, sites are added only at the periphery (surfaces, grain boundaries, or dislocations), so this homogeneous scaling may fail for finite crystals.","fun_headline_variants_meta":{"raw":{"variants":["Lattice sites: the missing variable in crystal thermodynamics","Solids break G=μN and Ω=−pV—new force ν appears","Adsorption equation for crystals from lattice-site force","Counting lattice sites fixes crystal thermodynamics"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00109,"raw_usage":{"total_tokens":4659,"prompt_tokens":1156,"completion_tokens":3503,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":772,"completion_tokens_details":{"reasoning_tokens":3438}},"tokens_in":772,"tokens_out":3503,"duration_ms":27630,"temperature":1.0,"reasoning_tokens":3438,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T21:19:09.696642+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Simulate a finite crystal at fixed $M$ under $(\\mu,p,T)$ control and measure $E=G-\\mu N$ by thermodynamic integration; the claim fails if $E/M$ is not an intensive function of occupation and stretch alone, or if the value of $\\nu$ inferred from $(\\partial\\nu/\\partial\\mu)_{p,T}=-c$ disagrees with the directly computed free-energy derivative with respect to $M$.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the network principle that the number of lattice sites is conserved under deformation, the foundational postulate from which $M$ becomes an extensive variable."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the Larché–Cahn theory of composition–stress coupling and the diffusion-potential framework that the paper extends."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the review context for the lattice-site potential $\\nu$ and its role in elastically stressed crystals."},{"cited_title":"Mullins ,\\ title title T hermodynamic equilibrium of a crystalline sphere in a fluid , \\ https://doi.org/10.1063/1.447779 journal journal J","cited_arxiv_id":null,"evidence_quote":"Shows that surface thermodynamics of crystals generates nonzero $\\nu$ through the difference between interface energy and surface stress, supporting the claim that $\\nu$ is generally nonvanishing."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Serves as the formal template for the extended thermodynamics; the replica energy in Hill's nanothermodynamics is analogous to $\\nu M$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the contrasting simulation scheme that treats $M$ as an unconstrained internal degree of freedom with $\\nu=0$, which the paper argues differs from a Larché–Cahn crystal."}],"review_version":1}