{"id":"304ef4de-04fc-4fc1-90ba-0bbbead0407b","arxiv_id":"2412.08143","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Symbolic regression rediscovers the exact von Weizsäcker kinetic-energy functional for one electron and Thomas-Fermi-like variants for many electrons in 1D, while showing semi-local forms fail for few-electron systems.","lead":"This paper uses symbolic regression to fit analytical kinetic-energy density functionals to exact one-dimensional quantum data, recovering the known von Weizsäcker limit for one electron and Thomas-Fermi-like forms for many electrons. It concludes that simple semi-local functionals built from the density and its derivatives are not accurate enough for few-electron systems, which matters for the design of orbital-free density functional theory.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Failure to find accurate semi-local KEDFs for few electrons may reflect the restricted SR operator set (+,×,÷, square only) and in-sample loss, not an intrinsic limit of semi-local functionals.","rationale":"The reader's weakest assumption concerns dataset representativeness and lack of a train/test split, which is a valid and important issue. I agree with that concern, but I identify an additional, more fundamental weakness: the restricted SR operator set and complexity penalty undermine the inference from 'SR did not find a good semi-local functional' to 'no good semi-local functional exists.' This is central because the paper's headline conclusion is precisely that semi-local forms are insufficient for few-electron systems. Both concerns are addressable, and the paper's positive results (exact recovery of vW for Ne=1, TF-like trend for larger Ne) remain valuable as an exploratory demonstration. The verdict of CONDITIONAL is therefore still appropriate, but the conditions should include both a held-out test set and evidence that the search space is expressive enough to represent accurate semi-local functionals. No verdict change is needed; the paper should be revised to qualify the negative claim and provide the missing experimental details.","tokens_in":15050,"tokens_out":5849,"duration_ms":62202,"concrete_test":"For Ne=2 and Ne=3, generate a new dataset of 200 random three-Gaussian potentials with explicitly reported parameter ranges, split into 150 training and 50 test systems. Train (a) the same PySR setup but with an expanded operator set including ρ^(1/3), ρ^(2/3), exp, log, and a larger complexity budget, and (b) a flexible neural-network-based semi-local KEDF taking ρ, ρ′, ρ″ as inputs. Compare test-set percentage loss. If either achieves test loss substantially below the reported in-sample values (e.g., <0.5%), then the paper's negative conclusion is an artifact of search restrictions rather than a genuine property of semi-local functionals.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central negative claim—that a semi-local expression of the KEDF is insufficient for few-electron systems (Introduction, Section IV)—rests on the SR search failing to find a low-loss functional. But the search space in Section III.A is extremely narrow: only the operators +, ×, ÷, and squaring are allowed, with inputs limited to ρ, ρ′, and ρ″. No fractional powers, exponentials, logarithms, or other nonlinear building blocks are available, so the search cannot represent many plausible semi-local forms. The paper itself notes in Section III.B that the second-best functionals often have substantially lower loss than the Pareto-front choice (Fig. 6), but they are discarded because of the complexity penalty—evidence that the search's complexity trade-off, not the intrinsic accuracy of semi-local forms, may be the limiting factor. Furthermore, Section II.D provides no train/test split, no sample count, and no explicit parameter ranges for the random Gaussian potentials, so the reported percentage losses (e.g., ~2% for Ne=2 in Fig. 4) are in-sample errors that could reflect overfitting to an unrepresentative ensemble. Thus the conclusion that semi-local functionals are insufficient for few electrons is an absence-of-evidence claim; without a broader search or a held-out test, it is not established.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses symbolic regression (SR), with a modified PySR implementation, to search for one-dimensional kinetic-energy density functionals (KEDFs) from exact-diagonalization ground states of non-interacting spinless electrons in random three-Gaussian external potentials. Three search strategies are presented: an unconstrained search over densities and their first two derivatives, a search seeded with the von Weizsäcker (vW) and Thomas-Fermi (TF) functionals as input features, and a search constrained to the form of a TF enhancement factor. For Ne=1 the SR recovers the vW functional with a fitted prefactor 0.12499846, very close to the exact 1/8; for Ne>=2 the best-scoring expressions are variants of the TF functional, with in-sample percentage losses decreasing with Ne and largest around 2% for Ne=2. The authors conclude that semi-local expressions are insufficient for the few-electron KEDF in 1D.","tokens_in":15363,"tokens_out":3537,"duration_ms":40077,"significance":"If the conclusions are robust, the paper provides a useful benchmark showing that symbolic regression can rediscover exact KEDF limits (vW for one electron, TF for many electrons) and it highlights the difficulty of constructing accurate semi-local KEDFs for few-electron systems. The exact-diagonalization data generation is simple and the Ne=1 recovery is a clean sanity check. However, the central negative claim rests on a search with a narrow operator set and on in-sample loss only; the manuscript does not report held-out tests, dataset sizes, or parameter ranges, and it does not validate functional derivatives. These omissions currently prevent the negative conclusion from being established as a general statement about semi-local KEDFs.","major_comments":[{"comment":"The dataset description in Section II.D gives no number of random potentials per Ne, no parameter ranges for Ai, bi, ci (aside from a vague statement that the Ai range is chosen according to Ne), and no train/test split. All reported percentage losses in Figs. 4, 5, 6, and 8 are therefore in-sample fit errors. This is load-bearing because the central negative claim, that semi-local expressions are insufficient for few-electron systems, depends on the losses not being an artifact of fitting a small or unrepresentative set of potentials. The authors should report the dataset size, the exact parameter sampling ranges, and a held-out test (e.g., cross-validation over potentials) before drawing conclusions about generalizability.","section":"II.D, Figs. 4-8"},{"comment":"The unconstrained search allows only the operators +, ×, ÷, and squaring, with inputs restricted to ρ, ρ', and ρ''. Many plausible semi-local KEDF forms, including fractional powers of the density such as ρ^α with non-integer α, cannot be represented exactly in this search space. The claim that 'a semi-local expression of the functional appears not sufficient' (Introduction) is therefore only supported with respect to this restricted operator set. Moreover, Fig. 6 shows that second-best functionals often have substantially lower loss than the Pareto-front choice but are discarded by the complexity penalty; this indicates that the parsimony criterion, not an intrinsic limitation of semi-local forms, may be responsible for the reported losses. The conclusions should be rephrased as applying to the specific SR search space and complexity trade-off used here.","section":"III.A"},{"comment":"The Introduction and Abstract motivate the work by the need for KEDFs with accurate functional derivatives for OFDFT minimization, and the paper states that analytical SR expressions make derivatives easy to compute. However, no functional derivative is ever computed or tested. Given that the manuscript explicitly frames derivative accuracy as a central advantage of SR, the authors should at least verify δTs/δρ for the Ne=1 and Ne=2 fitted functionals against the exact Kohn-Sham potential, or otherwise clearly delimit derivative accuracy as future work. As written, the claimed advantage for OFDFT is unsubstantiated.","section":"II.A and IV"},{"comment":"Symbolic regression is a stochastic global optimization procedure, yet the paper reports only a single Pareto front for each Ne and gives no statistics over independent SR runs or random seeds. The reported losses and functional forms could be non-reproducible if the genetic search is sensitive to initialization. The authors should provide repeated-run statistics (e.g., median and spread of losses across seeds) or justify why a single run is representative. This concern applies to all three searches and to the specific functional forms quoted in Figs. 4, 5, and 8.","section":"III.A-III.C"}],"minor_comments":[{"comment":"Equation (5) defines ρ(r) = Σ_i |φ_i(r)|, but the square is missing; it should be ρ(r) = Σ_i |φ_i(r)|^2.","section":"II.A, Eq. (5)"},{"comment":"The text refers to 'Born-von Karman boundary conditions' but states ψ(x0)=ψ(x_N−1)=0, which are hard-wall Dirichlet boundary conditions, not periodic Born-von Karman conditions.","section":"II.D"},{"comment":"The abstract and introduction contain the typo 'von von Weizsäcker' and later 'Plack constant' should be 'Planck constant'.","section":"Abstract and I"},{"comment":"The caption of Fig. 8 has 'grey dashed line os to guide the eye'; 'os' should be 'is'.","section":"III.C, Fig. 8"},{"comment":"The discussion of the λ prefactor in 'TF - λvW' says the standard deviation is small but 'statistically insignificant'; this phrasing is confusing because the standard deviation is a descriptive statistic, not a significance test. Clarify what is meant.","section":"III.B"},{"comment":"Reference [14] has the same title as reference [13] but appears to be a different paper (by del Mazo-Sevillano and Hermann); please verify the title and citation details.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper is an exploratory study and the central negative claim is currently stronger than the evidence supports. The missing dataset details and absent functional-derivative validation are fixable within the manuscript's scope, so I recommend major revision rather than rejection. I would also encourage the editor to ask the authors to provide data/code availability or at least the dataset construction parameters, since the paper's conclusions depend on reproducibility of the SR runs."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis is a useful paper to know about, but read it as a methods proof-of-concept, not as a settled negative result. The genuinely new piece is the first symbolic-regression search for KEDFs, and they had to modify PySR to evaluate candidate expressions as integrands rather than pointwise functions. That is a real technical contribution, and they show it works: for Ne=1 the unconstrained search recovers the von Weizsäcker functional with a prefactor of 0.12499846, essentially 1/8. The systematic scan over Ne=1..20 and the clean rediscovery of the TF limit at large Ne are well executed. I believe the exact-diagonalization data generation is sound.\n\nThe soft spots are not fatal to the demonstration, but they are important. First, there is no held-out test set, no reported dataset size or parameter ranges for the random Gaussian potentials, and no repeated-run statistics. The percentage losses in Figs. 4–8 are in-sample fit errors, so we do not know how the discovered functionals generalize. Second, and more damaging to the central conclusion, the unconstrained search only permits +, ×, ÷, and squaring, with inputs limited to ρ, ρ′, and ρ″. That search space cannot represent many plausible semi-local forms, e.g., fractional powers, exponentials, or logarithms. The paper itself shows in Fig. 6 that second-best functionals often have substantially lower loss than the Pareto-front choices but are discarded by the complexity penalty. So the claim that 'a semi-local expression of the functional appears not sufficient' for few electrons is an absence-of-evidence claim. The SR search failed within a narrow expression space; it does not establish a fundamental limitation of semi-local KEDFs. The authors do hedge by calling it a first experiment, but the abstract and conclusion state the insufficiency more strongly than the evidence supports.\n\nThe citation pattern looks fair. The related work on machine-learned KEDFs and SR for XC functionals is covered. I would send this to a serious referee: the method contribution is valuable, the vW rediscovery is clean, and the weaknesses are addressable with a test set, broader operators, and more runs. It deserves review time, but the referee should push for a more careful negative claim.","headline":"First symbolic-regression KEDF search, clean vW rediscovery; the negative semi-local claim is undercut by a narrow operator set and in-sample losses.","tokens_in":15866,"tokens_out":3257,"would_cite":true,"duration_ms":30300,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["71.15.Mb"],"model":"deepseek-v4-flash","headline":"The paper claims that symbolic regression can rediscover the exact von Weizsäcker kinetic-energy functional for a single 1D electron, but for two or more electrons the best semi-local fits are Thomas-Fermi-like and leave the largest error…","keywords":["symbolic regression","kinetic energy density functional","orbital-free DFT","Thomas-Fermi functional","von Weizsäcker functional","one-dimensional electron gas","enhancement factor"],"falsifier":"Train or test the same symbolic-regression searches on a different ensemble of 1D potentials, such as shallow well potentials, periodic potentials, or potentials producing multi-peaked densities for Ne = 2 to 5, and check whether any semi-local expression built from ρ, ρ′, and ρ″ reaches near-zero loss on a held-out set of exact ground states. A held-out test yielding near-zero error for a semi-local expression at intermediate electron counts would falsify the paper's negative claim, while a demonstration that the reported losses grow substantially out-of-sample would support it.","tokens_in":14858,"feed_emoji":"⚛️","tokens_out":6655,"duration_ms":61307,"temperature":0.7,"pith_summary":"This paper tries to establish that symbolic regression, a search over closed-form expressions, can rediscover exact kinetic-energy density functionals in one-dimensional systems, but that semi-local functional forms are not enough to represent the non-interacting kinetic energy for few-electron systems. Using exact ground-state densities from non-interacting electrons in random three-Gaussian wells, the search recovers the von Weizsäcker functional for one electron and variants of the Thomas-Fermi functional for every higher electron count. The largest errors appear at two electrons, which the authors attribute to the non-local nodal structure of orbitals that a semi-local expression built from density and its first two derivatives cannot capture. If correct, the result sharpens the known difficulty of orbital-free DFT: accurate kinetic-energy functionals likely need non-local information, not just local density and derivative terms.","feed_headline":"Symbolic regression finds exact 1D kinetic functional for one electron","feed_subtitle":"For two or more electrons, only Thomas-Fermi-like forms emerge, and semi-local kinetic functionals fall short.","key_machinery":"The machinery is symbolic regression adapted to functionals: candidate expressions are built from the density, its first and second derivatives, and, in constrained runs, from the dimensionless reduced variables s(x), q(x), and k(x); each candidate kinetic-energy density is integrated over space and compared to the exact Ts, with complexity penalized. The core identity is the uniform-scaling law Ts[ρλ] = λ³Ts[ρ], which is enforced in the enhancement-factor search by writing Ts[ρ] = (π²/6)∫ρ³ F(s,q,k) dx. The Pareto front of loss versus expression complexity is the selection device that decides which functional wins at each electron number.","core_discovery":"The central claim is that, for one-dimensional non-interacting electron systems, the only accurately learnable semi-local kinetic-energy density functionals are the two exact limits: the von Weizsäcker functional for one electron and Thomas-Fermi-like expressions for many electrons. The search across electron numbers 1 to 20 shows a smooth transition in which the vW functional is recovered exactly at one electron and the TF functional, with a fitted prefactor below the exact 1D value π²/6, dominates the Pareto front at all other electron counts. The authors interpret the elevated losses at intermediate electron counts, especially two electrons, as evidence that the semi-local paradigm is fundamentally limited: orbital orthogonality imposes nodal constraints that cannot be expressed locally through the density and its derivatives.","pith_inferences":["A natural extension the paper leaves implicit is to split the functional into a fixed non-local kernel plus a learned semi-local correction; the authors mention this avenue but do not test it, and the negative result here is what makes that split worth trying.","Since the losses are reported in-sample, the gap between semi-local and exact might grow when measured on densities outside the training set; a train/test split could make the few-electron failure more severe than the figures suggest.","The paper does not evaluate the variational accuracy of the functional derivatives of the discovered expressions; if those derivatives are inaccurate, even the recovered vW and TF forms may not be directly usable in an orbital-free minimization.","The same symbolic-regression setup could be applied to 3D confined systems or to electrons with spin, where the density has nodal surfaces; a failure there would extend the paper's conclusion beyond 1D."],"forward_implications":["At one electron, symbolic regression returns the exact von Weizsäcker functional with a percentage loss of about 0.00027%, so the method can rediscover a known exact limit without being told it.","At every other electron count, the best semi-local form is a Thomas-Fermi-like functional whose fitted prefactor stays below π²/6 and approaches about 1.635 for large Ne, so the 1D TF functional systematically overestimates the kinetic energy of confined non-uniform systems.","The largest percentage loss occurs at Ne = 2, identifying the few-electron regime as the hardest case for semi-local kinetic-energy density functionals.","When the search is seeded with vW and TF, a 'TF − λvW' linear combination becomes the second-best expression for intermediate Ne with λ ≈ 0.265 ± 0.008 in 1D, while in 3D the optimal combination reported in the literature is TF + λvW with λ ≈ 0.2.","Constraining the search to an enhancement factor that respects uniform scaling does not change the picture: the winning forms are still vW for one electron and TF-like variants otherwise, with the extra dimensionless terms integrating to essentially zero."],"supporting_citations":[{"why":"Supplies the Thomas-Fermi functional, the target that the many-electron searches converge to.","marker":"[28]"},{"why":"Supplies the Thomas-Fermi functional alongside Thomas, the target of the many-electron searches.","marker":"[29]"},{"why":"Defines the von Weizsäcker functional, which the search recovers exactly for one electron.","marker":"[31]"},{"why":"Provides the 1D Thomas-Fermi kinetic-energy density π²/6 ρ(x)³ that the symbolic-regression results are compared against.","marker":"[18]"},{"why":"Provides the 3D TF + λvW optimal coefficient (about 0.2) that the paper contrasts with its 1D TF − λvW result.","marker":"[32]"},{"why":"Supplies Sturm's separation theorem, which underpins the argument that orbital nodal structure defeats semi-local forms.","marker":"[63]"},{"why":"Establishes the enhancement-factor symbolic-regression strategy for density functionals that the paper adapts to the kinetic-energy density functional.","marker":"[47]"},{"why":"Is the symbolic-regression implementation the authors modified to work with functionals rather than plain functions.","marker":"[58]"}],"fun_headline_variants":["Symbolic regression recovers exact 1D kinetic limits","1D kinetic functional search lands on vW and Thomas-Fermi","Machine learning 1D KE: only exact limits emerge","Semi-local 1D kinetic functionals stop at exact limits"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The conclusion that semi-local functionals are insufficient for few-electron systems rests on treating the random three-Gaussian external potentials used to generate the training densities as representative of all one-dimensional non-interacting electron systems; if that potential ensemble does not sample the density space where semi-local functionals could succeed, the negative result is a property of the dataset rather than of semi-local functionals in general.","fun_headline_variants_meta":{"raw":{"variants":["Symbolic regression recovers exact 1D kinetic limits","1D kinetic functional search lands on vW and Thomas-Fermi","Machine learning 1D KE: only exact limits emerge","Semi-local 1D kinetic functionals stop at exact limits"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000227,"raw_usage":{"total_tokens":1434,"prompt_tokens":868,"completion_tokens":566,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":484,"completion_tokens_details":{"reasoning_tokens":494}},"tokens_in":484,"tokens_out":566,"duration_ms":5824,"temperature":1.0,"reasoning_tokens":494,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T18:09:53.858620+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Train or test the same symbolic-regression searches on a different ensemble of 1D potentials, such as shallow well potentials, periodic potentials, or potentials producing multi-peaked densities for Ne = 2 to 5, and check whether any semi-local expression built from ρ, ρ′, and ρ″ reaches near-zero loss on a held-out set of exact ground states. A held-out test yielding near-zero error for a semi-local expression at intermediate electron counts would falsify the paper's negative claim, while a demonstration that the reported losses grow substantially out-of-sample would support it.","supporting_citations":[{"cited_title":"The calculation of atomic fields,","cited_arxiv_id":null,"evidence_quote":"Supplies the Thomas-Fermi functional, the target that the many-electron searches converge to."},{"cited_title":"Un metodo statistico per la determinazione di alcune priorietá dell’atomo,","cited_arxiv_id":null,"evidence_quote":"Supplies the Thomas-Fermi functional alongside Thomas, the target of the many-electron searches."},{"cited_title":"Zur theorie der kernmassen,","cited_arxiv_id":null,"evidence_quote":"Defines the von Weizsäcker functional, which the search recovers exactly for one electron."},{"cited_title":"Understanding machine-learned density func- tionals,","cited_arxiv_id":null,"evidence_quote":"Provides the 1D Thomas-Fermi kinetic-energy density π²/6 ρ(x)³ that the symbolic-regression results are compared against."},{"cited_title":"On the Weizsäcker Correc- tion to the Thomas-Fermi Theory of the Atom,","cited_arxiv_id":null,"evidence_quote":"Provides the 3D TF + λvW optimal coefficient (about 0.2) that the paper contrasts with its 1D TF − λvW result."},{"cited_title":"Teschl,Ordinary Differential Equations and Dynami- cal Systems, Graduate studies in mathematics (American Mathematical Soc.)","cited_arxiv_id":null,"evidence_quote":"Supplies Sturm's separation theorem, which underpins the argument that orbital nodal structure defeats semi-local forms."},{"cited_title":"Evolv- ing symbolic density functionals,","cited_arxiv_id":null,"evidence_quote":"Establishes the enhancement-factor symbolic-regression strategy for density functionals that the paper adapts to the kinetic-energy density functional."}],"review_version":1}