{"id":"f897926f-b2e5-4c17-b2c7-50fb51357c09","arxiv_id":"2502.02378","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A new excited-state uniform electron gas with a Fermi-surface gap yields closed-form kinetic and exchange energies and a gap-dependent leading correlation coefficient in the high-density limit.","lead":"This paper proposes a model of an excited electron gas by moving electrons near the Fermi surface across a gap, and derives closed-form formulas for its kinetic and exchange energy plus the leading correlation correction. The goal is to provide a reference system for building density functionals that can describe excited states in molecules and materials.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equation (25) fails its own zero-gap limit by a factor of 2, so the central correlation coefficient is not correct as printed; the six-channel derivation is also only sketched.","rationale":"The reader's verdict is CONDITIONAL, and my independent check of Eq. (25) confirms the specific factor-of-2 inconsistency the reader noted in the rationale: at Δ=0 the printed formula yields 2λ0 rather than λ0. This is the most load-bearing concern because it directly undermines the central new result, the claimed closed-form correlation coefficient. The kinetic (Eq. (15)) and exchange (Eq. (18)) expressions are derived from direct integrals over the occupation pattern and their limiting values are consistent, so they are not the weak point. I did not make the transferability issue (mapping a global Δ to a local variable) the load-bearing concern: the conclusion explicitly acknowledges that this mapping 'requires further investigation,' so it is an honest, stated limitation rather than a hidden flaw in the closed-form results. The paper's own text asserts the zero-gap limit, so the factor-of-2 error is a demonstrable internal inconsistency, not a disagreement with consensus. The derivation of the six-term formula is only sketched, so the Δ-dependence must be independently verified, which is what the concrete numerical test addresses. The verdict should remain CONDITIONAL because the error appears correctable and the overall model plus kinetic/exchange results are plausible, but the correlation formula must be corrected and validated before the framework is usable.","tokens_in":11133,"tokens_out":11464,"duration_ms":107276,"concrete_test":"Numerically extract the logarithmic coefficient of the direct second-order energy by integrating Eq. (21) with the excited-state occupation (10) for several values of Δ (e.g., 0, 0.25, 0.5, 0.75, 1). Use a small-k cutoff δ, integrate over the allowed p,q and k regions, and fit the coefficient of ln δ for each Δ. Compare these coefficients with Eq. (25) using both the printed prefactor 1/π^2 and the corrected prefactor 1/(2π^2). If the numerical coefficient at Δ=0 is (1−ln2)/π^2 and matches Eq. (25) after the prefactor correction at all sampled Δ, the correlation formula is vindicated; otherwise the six-term sum or its prefactor needs revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central new result is the leading high-density correlation coefficient λ0(Δ) in Eq. (25). As printed, Eq. (25) is internally inconsistent with its stated ground-state limit. Setting Δ=0 forces κ=1, so the six terms reduce to F(1,1)=2(1−ln2) and the sum Σλ0^(k)=F(1,1). With the printed prefactor 1/π^2, this gives λ0(0)=2(1−ln2)/π^2≈0.0621814, exactly twice the Macke coefficient λ0=(1−ln2)/π^2≈0.0310907 that the very next sentence asserts. This is not a rounding discrepancy but a factor-of-2 error, forcing either a correction of the prefactor to 1/(2π^2) or a revision of the six-term enumeration. All downstream numbers—Fig. 4, the stated value λ0(1)≈0.00578826, and any functional built on Λ0(Δ)—change accordingly. Moreover, the derivation preceding Eq. (25) is only a sketch: the six channels are enumerated verbally, but the explicit reductions of Eq. (21) to the functions F(α,β) are not shown. Consequently, the Δ-dependence of the claimed closed form is not verifiable from the text as it stands. Because the correlation coefficient is the load-bearing new result, the paper's central claim is not supported as printed.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces a model of an excited-state uniform electron gas (UEG) in which a fraction of electrons near the Fermi surface is promoted across a gap, with occupation given by the step function in Eq. (10). The density-matching condition fixes the parameter κ(Δ), and the author derives closed-form expressions for the reduced kinetic and exchange energies of this model as functions of density and gap, Eqs. (14)-(18). The central new result is the leading high-density direct correlation coefficient λ0(Δ), Eq. (25), obtained from second-order perturbation theory and a six-channel enumeration of logarithmically divergent momentum-transfer processes. The paper argues that such excited-state UEGs can serve as a reference system for constructing state-specific local and semilocal density functionals.","tokens_in":11384,"tokens_out":7477,"duration_ms":69924,"significance":"If the results are correct, the kinetic and exchange closed forms provide a clean, parameter-free extension of the UEG paradigm, and the correlation coefficient λ0(Δ) would be a genuinely new quantity for excited-state DFT. The model has no fitted parameters, with κ fixed by density conservation, and the exchange-hole expression in Eq. (20) is a simple and appealing combination of ground-state holes. However, the central correlation result as printed is internally inconsistent by a factor of 2 in both stated limits, and the derivation leading to Eq. (25) is only sketched. Because λ0(Δ) is the load-bearing new result, the paper in its current form does not support the claimed correlation-energy coefficient. The kinetic and exchange parts appear sound and are not affected by this issue.","major_comments":[{"comment":"Equation (25) is inconsistent with the stated ground-state limit by a factor of 2. Setting Δ=0 forces κ=1, and the six terms reduce to F(1,1) + F(1,1) + F(1,1) - 2F(1,1) - 2F(1,1) + 2F(1,1) = F(1,1). Since F(1,1)=2(1-ln2), the printed formula gives λ0(0)=2(1-ln2)/π^2 ≈ 0.06218, exactly twice the value λ0=(1-ln2)/π^2 ≈ 0.03109 asserted in the following sentence. Similarly, evaluating Eq. (25) at Δ=1 gives λ0(1) ≈ 0.0116, not the stated 0.00578826. Either the prefactor 1/π^2 should be 1/(2π^2), or the six-term enumeration must be revised; all downstream values and Fig. 4 depend on this correction.","section":"III, Eq. (25)"},{"comment":"The derivation of the six-channel result is only a verbal enumeration of the divergent processes; the reduction of Eq. (21) to the functions F(α,β) in Eq. (26) is not shown. In particular, the text does not display how the admissible regions of p, q, and k produce the specific prefactors (1-Δ)^3, (1+κΔ)^3, and the mixed terms -2F(1-Δ,1), -2F(1,1+κΔ), and 2F(1-Δ,1+κΔ). Without this algebra, the claimed Δ-dependence cannot be independently verified. Please provide the full derivation, for example in an appendix or supplemental material.","section":"III, derivation of Eq. (25)"},{"comment":"The proposed construction of state-specific semilocal functionals relies on mapping the global UEG parameter Δ onto a local degree-of-excitation variable at each point of an inhomogeneous system. This mapping is not constructed or tested in the manuscript, and the conclusion acknowledges that it 'requires further investigation.' Since the abstract claims the paper 'introduces a new framework' for building such functionals, the unsupported nature of this transferability claim should be stated more prominently, or the claim should be tempered to refer to the formal UEG results only.","section":"IV, Conclusion"}],"minor_comments":[{"comment":"The notation k·(p−q+k) in the denominator of Eq. (21) is ambiguous; it should be written explicitly as k·(p−q)+|k|^2, or with vector arrows, so that the reader can follow the angular integrations.","section":"III, Eq. (21)"},{"comment":"The passage 'k·p = kpx and k·q = −kqy' introduces coordinates x and y without defining the reference frame; please clarify the coordinate choice and the resulting integration limits.","section":"III, Eq. (22)"},{"comment":"The statement that the upper limit '√rs < k < 1' is arbitrary should be clarified: the coefficient of ln rs is independent of the cutoff, but the finite part is not; this distinction is worth stating to avoid confusion.","section":"III, Eq. (23)"},{"comment":"After correcting Eq. (25), the value λ0(1), the curve in Fig. 4, and the normalization Λ0(Δ)=λ0(Δ)/λ0 must be recomputed and reported consistently.","section":"III, Fig. 4 and text after Eq. (25)"},{"comment":"The exchange coefficient Ξx is given by a long expression; at Δ=0 it correctly reduces to 1, but the reader would benefit from a brief indication of how the logarithmic terms arise from the angular integrations over the two-shell Fermi hole.","section":"III, Eq. (18)"}],"recommendation":"major_revision","confidential_remarks":"The factor-of-2 discrepancy in Eq. (25) is very likely a typo in the prefactor (1/(2π^2) instead of 1/π^2), but because the derivation is not shown in full, the referee cannot fully exclude a deeper issue in the six-channel enumeration. The kinetic and exchange parts are straightforward and appear correct. The paper is within the scope of the journal and the model is interesting, but the central correlation coefficient must be corrected and the derivation completed before publication. A major revision rather than rejection seems appropriate, since the error appears fixable within the manuscript's scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is worth a look: it defines a density-preserving excited-state UEG with a Fermi-surface gap, and the kinetic and exchange closed forms (Eqs. 15 and 18) are clean, direct integrals that check out. That part is solid and useful. The six-channel decomposition for the leading correlation coefficient is a real attempt to go beyond the jellium-with-gap models, and the idea of a gap-dependent Λ0 is a sensible extension of the Macke/Gell-Mann-Brueckner machinery.\n\nThe soft spot is exactly where the reader and the stress-test put it. Equation (25) as printed fails its own zero-gap limit by a factor of 2: setting Δ=0 gives Λ0=2 with the stated prefactor, so λ0(0) comes out twice the Macke coefficient that the next sentence asserts. This is not a rounding issue. A missing 1/2 in the prefactor, or a revised channel enumeration, fixes it, but as written the central new result is numerically wrong and all downstream numbers shift. The derivation of the six terms is also only sketched; the text lists the channels verbally but does not show how Eq. (21) reduces to the F(α,β) functions. That matters because the Δ-dependence is the whole point.\n\nThe transferability concern is real but honestly flagged by the author in the conclusion. The mapping from a global gap parameter to a local variable is the foundation for any state-specific LDA, and the paper admits it is unvalidated. That does not sink the model as a reference system; it just means the functional part is speculative, which the author says plainly.\n\nOn balance: the kinetic and exchange results are publishable as they stand, and the correlation part is likely correct once the prefactor is fixed and the derivation is expanded. The paper deserves a serious referee, but not in its current form. I would not cite it in the next twelve months, because the load-bearing correlation coefficient cannot be used until corrected. If the author fixes Eq. (25) and shows the six-channel reductions explicitly, it becomes a useful reference for excited-state DFT work. I'd bring it to reading group, mostly to see whether the six-channel enumeration holds up under close inspection.","headline":"Genuinely new excited-state UEG model with clean kinetic/exchange results, but the central correlation coefficient has a factor-of-2 error and the derivation is under-sketched.","tokens_in":11969,"tokens_out":1112,"would_cite":false,"duration_ms":12499,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Gapped electron gas yields exact kinetic and exchange energies","keywords":["uniform electron gas","excited states","density functional theory","local density approximation","exchange energy","correlation energy","kinetic energy","Fermi surface gap"],"falsifier":"Evaluate the second-order direct correlation integral of Eq. (21) numerically with the excited-state occupation factors for several values of Δ; if the extracted coefficient of ln r_s disagrees with the six-term formula in Eq. (25), the list of divergent channels is incomplete. The kinetic and exchange formulas can likewise be checked by direct numerical integration of Eqs. (14) and (17).","tokens_in":10852,"feed_emoji":"⚛️","tokens_out":5838,"duration_ms":56489,"temperature":0.7,"pith_summary":"This paper extends the uniform electron gas (UEG), the model behind the local-density approximation in density-functional theory, to excited states. It defines an excited-state UEG by opening a gap at the Fermi surface: electrons in a fraction Δ of the occupied band are excited to states just above the Fermi level instead of sitting just below it. For this model it derives closed-form expressions for the reduced kinetic and exchange energies as functions of density and Δ, and the leading ln r_s term of the correlation energy in the high-density limit. These formulas provide a reference system for building local and semilocal density functionals that are specific to individual excited states.","feed_headline":"Gapped electron gas yields exact kinetic and exchange energies","feed_subtitle":"Leading correlation term also derived, opening a route to state-specific density functionals.","key_machinery":"The machinery is the occupation pattern of Eq. (10): occupied plane-wave states up to k_F(1−Δ) and between k_F and k_F(1+κΔ), with the factor κ chosen so that the spin density matches the ground-state value. All energy expressions are integrals over these occupied states. The gap-dependent scaling functions Ξ_s(Δ) and Ξ_x(Δ) encode the change in kinetic and exchange energy relative to the ground state at fixed density, while the correlation coefficient λ_0(Δ) is assembled from six divergent second-order processes, each contributing a term of the form F(α,β)/π² with F(α,β) = α²β + αβ² + α³ ln α + β³ ln β − (α³+β³) ln(α+β). Together these functions convert a global excitation fraction into modifications of the standard UEG energy scalings.","core_discovery":"The central result is that an excited-state UEG with the step occupation pattern of Eq. (10) has a reduced kinetic energy t_sσ = Ξ_s(Δ_σ) C_F $ρ_σ^{{2/3}}$, a reduced exchange energy ε_xσ = Ξ_x(Δ_σ) C_x $ρ_σ^{{1/3}}$, and a high-density direct correlation energy $ε^{{(2d)}}$ ~ λ_0(Δ) ln r_s. The dimensionless factors Ξ_s and Ξ_x are closed-form functions of the excitation fraction Δ, with Ξ_x also containing the density-matching factor κ_σ(Δ) and logarithmic terms; λ_0(Δ) is a sum over six momentum-transfer channels, each giving a logarithmic divergence, and it reduces to the known ground-state coefficient (1 − ln 2)/π² when Δ → 0. The formulas recover the Thomas-Fermi, Dirac, and Gell-Mann-Brueckner results in the ground-state limit and provide a pure-state excited-state generalization of the UEG.","pith_inferences":["The Δ parameter could be linked to local descriptors such as the fraction of excited electrons in a region or natural-orbital occupation numbers, providing a concrete path from the uniform model to molecules.","The six-channel logarithmic decomposition may transfer to coupled-cluster theory, where the infrared catastrophe in metals is handled by similar resummations; the excited-state UEG offers a testbed for those resummations.","Spin-forbidden transitions would require two coupled gaps, one per spin channel; the formulas already allow Δ_↑ ≠ Δ_↓, so an extension to singlet-triplet excitations may be straightforward.","The low-density limit is expected to be degenerate with the ground state based on cited thermodynamic arguments, which the paper does not itself prove; a numerical check in the Wigner-crystal regime would settle that expectation."],"forward_implications":["The closed-form Ξ_s and Ξ_x allow immediate construction of excited-state local exchange and kinetic functionals by replacing the constants C_F and C_x with Δ-dependent coefficients in an LDA-type expression.","The high-density λ_0(Δ) ln r_s term fixes the exact leading correlation behavior of a gapped electron gas, giving a target that any state-specific correlation functional must reproduce in the high-density limit.","Because the exchange hole remains normalized and tightens as Δ increases, ground-state LDA machinery such as hole-based corrections can be adapted to excited states with the same formal guarantees.","The model reduces continuously to the ground-state UEG at Δ = 0, so any functional built on it interpolates between known ground-state and new excited-state limits."],"supporting_citations":[{"why":"Gives the ground-state direct correlation energy and the original logarithmic divergence that the excited-state analysis extends.","marker":"[56]"},{"why":"Establishes the ln r_s high-density expansion and the screening cutoff used in the correlation integral.","marker":"[57]"},{"why":"Supplies the explicit second-order direct-term integral (Eq. 21) from which the six divergent channels are extracted.","marker":"[72]"},{"why":"Introduces the jellium-with-a-gap model that the excited-state UEG generalizes by shifting unoccupied states.","marker":"[35]"},{"why":"Presents the constant-occupation-factor ensemble UEG for excited states, the alternative approach this pure-state model is contrasted with.","marker":"[11]"},{"why":"Provides the low-density degeneracy argument used to justify focusing on the high-density limit.","marker":"[71]"},{"why":"The standard review of ground-state UEG formulas that the kinetic and exchange derivations build on.","marker":"[33]"}],"fun_headline_variants":["Excited electron gas yields exact kinetic and exchange energies","Exact kinetic and exchange energies for excited electron gas","Exact kinetic and exchange energies for gapped electron gas","Excited electron gas: exact kinetic and exchange energies","Gapped electron gas: exact kinetic and exchange energies"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole program rests on the assumption that the single global fraction Δ of the uniform gas can be replaced by a local 'degree of excitation' at each point of an inhomogeneous molecule; the paper states that this mapping remains to be constructed.","fun_headline_variants_meta":{"raw":{"variants":["Excited electron gas yields exact kinetic and exchange energies","Exact kinetic and exchange energies for excited electron gas","Exact kinetic and exchange energies for gapped electron gas","Excited electron gas: exact kinetic and exchange energies","Gapped electron gas: exact kinetic and exchange energies"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001789,"raw_usage":{"total_tokens":7017,"prompt_tokens":881,"completion_tokens":6136,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":497,"completion_tokens_details":{"reasoning_tokens":6059}},"tokens_in":497,"tokens_out":6136,"duration_ms":42433,"temperature":1.0,"reasoning_tokens":6059,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T12:22:50.344755+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate the second-order direct correlation integral of Eq. (21) numerically with the excited-state occupation factors for several values of Δ; if the extracted coefficient of ln r_s disagrees with the six-term formula in Eq. (25), the list of divergent channels is incomplete. The kinetic and exchange formulas can likewise be checked by direct numerical integration of Eqs. (14) and (17).","supporting_citations":[{"cited_title":"Raimes ,\\ @noop title Many-Electron Theory \\ ( publisher Amsterdam, North-Holland ,\\ year 1972 ) NoStop","cited_arxiv_id":null,"evidence_quote":"Supplies the explicit second-order direct-term integral (Eq. 21) from which the six divergent channels are extracted."},{"cited_title":"Callaway ,\\ title title Correlation energy in a model semiconductor , \\ https://doi.org/10.1103/PhysRev.116.1368 journal journal Phys","cited_arxiv_id":null,"evidence_quote":"Introduces the jellium-with-a-gap model that the excited-state UEG generalizes by shifting unoccupied states."}],"review_version":1}