{"id":"6adf67f8-2895-46d1-b39b-be7bac52ad6c","arxiv_id":"1908.09063","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":3.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Fitting a power-law correction to two noble-gas atoms recovers the LDA exchange functional's exponent within 0.3%, but with no held-out test the benchmark is more a consistency check than a prediction.","lead":"This paper benchmarks a method that improves approximate density functional theory by fitting a power-law correction to two noble-gas atoms and then checking whether it recovers the known LDA exchange functional. The recovered exponent is within 0.3% of the true value, but the test is a self-consistency check because the same atoms both set and evaluate the fitted parameters.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (9) is an unvalidated power-law ansatz, and the benchmark fits and evaluates on the same atoms, so the reported improvement is an in-sample consistency check rather than an out-of-sample validation.","rationale":"The reader's weakest assumption matches my own reading. The derivation of Eq. (7) is internally plausible, and the numerical implementation appears consistent, so I do not see a mathematical contradiction in the formalism. The soft spot is external validity: Eq. (9) is arbitrary, and the benchmark is constructed so that the target lies inside the ansatz and the evaluation set overlaps the fit set. This does not falsify the method, but it leaves the central generalization claim conditional. A held-out test with a non-power-law target is the decisive check. I therefore keep the reader's CONDITIONAL verdict.","tokens_in":5567,"tokens_out":6163,"duration_ms":61052,"concrete_test":"Fit A and α on two noble-gas atoms (e.g., He and Xe) using Eq. (10), with \\tilde E_Hxc pure Hartree and E_target changed to Hartree + LDA exchange + VWN correlation, whose density dependence is not a single power of ρ. Then predict ground-state energies and densities for Ne, Ar, Kr, and Rn and compare with target values. If the transfer errors remain at the 1-5% level seen in Tables 1-2, Eq. (9) is adequate; if errors are large or the fitted parameters shift strongly with the choice of the fitting pair, the power-law ansatz is the limiting assumption and the central claim is not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing assumption is Eq. (9), E_Hxc^(1)[ρ] = A ∫ ρ^α dr, with A and α fixed by two systems. Eq. (7) is a linear functional equation for E_Hxc^(1), and Eq. (9) is a tractability ansatz, not a consequence of DFT or of the perturbation expansion. In the benchmark this is exactly the form of the target: \\tilde E_Hxc is pure Hartree and E_target adds only the Dirac term -C∫ρ^{4/3} (Eq. 11), so the solution space of the ansatz contains the target by construction. Moreover, the parameters are fitted to Ar and Kr (or another pair) and then Tables 1-2 and Fig. 2 report energies and densities for those same atoms, including the Ar-Kr case. The quoted 2-3 orders of magnitude improvement is therefore an in-sample consistency check of the inversion, not an out-of-sample prediction. If the true correction is not of the single-power-law form over the relevant density range, the fitted A and α will absorb projection error and the improved functional need not transfer to other systems.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This proceedings paper presents the IKS-DFPT method for improving approximate energy density functionals. The idea is to combine the inverse Kohn-Sham method with first-order density functional perturbation theory: given a known approximate functional and an exact ground-state density, the first-order correction to the Hartree-exchange-correlation functional is determined from a functional equation, Eq. (7). To make the equation tractable, the correction is assumed to have the power-law form E_Hxc^(1)[ρ] = A ∫ ρ^α dr in Eq. (9), and the two parameters λA and α are fixed using two noble-gas atoms. The benchmark uses the Hartree functional as the approximate starting point and the Hartree plus LDA exchange functional as the target, and reports that the fitted coefficients and ground-state energies are close to the target for Ar and Kr.","tokens_in":5768,"tokens_out":3893,"duration_ms":44841,"significance":"If the method is validated, it would offer a systematic route to improving approximate functionals from exact densities, and the extension to nuclear DFT mentioned in the conclusion could be of interest. The paper is honest about the central assumption in Eq. (9), and the derivation of Eq. (7) is explicit and internally consistent. However, the benchmark as presented is not an out-of-sample validation: the correction ansatz has exactly the same functional form as the target exchange term, and the fitted parameters are evaluated on the same systems used for the fit. The reported two-to-three-order improvements in energy are therefore an in-sample consistency check rather than evidence that the method improves conventional functionals in general. The strength of the formal idea is real, but the central claim of the paper currently rests on a benchmark that is partly guaranteed by construction.","major_comments":[{"comment":"The benchmark cannot distinguish the method from a two-parameter fit because the assumed correction in Eq. (9) has exactly the same power-law form as the target exchange correction in Eq. (11). The target differs from the starting Hartree functional by -C∫ρ^{4/3}dr, and the ansatz is A∫ρ^αdr; for α=4/3 and λA=-C, the target is inside the ansatz class by construction. A meaningful test would use a target functional that is not of this form, for example a gradient-dependent functional or a functional with a different density dependence, and quantify how well the two-parameter ansatz can represent it. Without such a test, the statement that the method is promising for improving conventional functionals remains unsupported.","section":"Sec. 2, Eq. (9) and Sec. 3, Eq. (11)"},{"comment":"The parameters λA and α are determined from the same atoms on which the energy and density improvements are then reported: Table 1 reports Ar and Kr after fitting to Ar and Kr, and Table 2 reports each pair after fitting to that same pair. The quoted improvement by two to three orders of magnitude in ground-state energies is therefore an in-sample residual of the inversion, not a predictive result. Please provide out-of-sample tests, such as fitting on Ar and Kr and then reporting errors for Ne, Xe, or another atom, including density errors before and after improvement. This is necessary to support the transferability claim in Sec. 4.","section":"Sec. 3, Tables 1-2 and Fig. 2"},{"comment":"The conclusion states that the accuracy of ground-state energies is improved by two to three orders of magnitude and that the IKS-DFPT is promising to improve conventional functionals. This overstates what the benchmark shows: Table 2 reports errors in α of about 0.2–1.0% and errors in λA of about 2.3–7.6%, and these are errors in recovering a functional that lies exactly inside the assumed ansatz. The two-to-three-order improvement in energies is a consequence of using the same target-derived quantities in the fit and in the evaluation. The conclusion should be limited to the statement that the inversion works when the correction has the assumed form, and the open question of transferability to other functional forms should be stated explicitly.","section":"Sec. 4, conclusion"}],"minor_comments":[{"comment":"The phrase 'their errors with respect to the target valued are shown' contains a typo; it should be 'target values'.","section":"Table 2 caption"},{"comment":"The mass m in the Kohn-Sham equation is not specified; since the benchmark uses Hartree atomic units, the notation would be clearer if m is explicitly set to the electron mass or the atomic-unit convention is stated.","section":"Sec. 2, Eq. (2)"},{"comment":"The sentence 'the KS potential V_KS(r) is unique concerning the system' is imprecise; the Kohn-Sham potential is unique up to an additive constant for a given density, and the statement could be clarified to refer to the Hohenberg-Kohn and Kohn-Sham uniqueness theorems.","section":"Sec. 2, after Eq. (7)"},{"comment":"The inset in Fig. 2 is very small and the ratio rs/r_s^target is hard to read; a separate panel or larger inset would improve readability.","section":"Fig. 2"}],"recommendation":"major_revision","confidential_remarks":"The formal derivation in Sec. 2 is sound, but the benchmark design makes the headline improvement an in-sample fit. I would encourage the editor to request out-of-sample tests or a target functional outside the ansatz class; if those cannot be provided, the paper should be reframed as a proof-of-principle of the inversion step rather than as evidence that conventional functionals can be improved generally. This is a proceedings contribution, so the scope of the requested revision could reasonably be limited to reframing the claims and adding at least one held-out atom."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this is a follow-up proceedings from the same group. The IKS-DFPT formalism already appears in their earlier preprint (Ref. [5]); the new material here is a benchmark for all noble-gas pairs, with Ar and Kr as the focus. If you come looking for new theory, you won't find it. What the paper does well: the derivation in Sec. 2 is clear and the comparison of the two energy expressions gives a legitimate functional equation for the first-order correction, Eq. (7). The authors are upfront that Eq. (9) is an assumed power-law form and that two reference systems are needed to fix the parameters. The error tables are systematic and the trend—heavier atoms, wider density range, better parameter accuracy—makes physical sense. The density improvement shown in Fig. 2 is real, for the cases tested.\n\nThe soft spot is the one that matters: the target in the benchmark is LDA exchange, which is exactly a single power of the density, ρ^{4/3}, and Eq. (9) is the same functional family. So the solution space of the ansatz contains the target by construction. Worse, λA and α are fitted to the same atoms whose energies and densities are then reported as improved in Tables 1 and 2. That makes the benchmark an in-sample consistency check, not an out-of-sample validation. The order-of-magnitude energy gains also look partly like a baseline effect: the Hartree starting point is very crude, so adding a fitted exchange term has a lot of room to help. The conclusion that IKS-DFPT is 'promising to improve the conventional functionals' overreaches the evidence. A held-out test, fitting on one pair and evaluating on another, or using a target whose form is not known to sit inside the ansatz, would substantially change my confidence.\n\nRead as a verification that the inversion machinery works when the correction has the assumed form, the paper is acceptable. Read as a demonstration of general functional improvement, it falls short. The novelty is low and the claims are stronger than the benchmark supports. People working on functional construction—especially in nuclear DFT—will find it a compact example of how easy it is to mistake an in-sample fit for a predictive method. I'd use it in a reading group as a cautionary tale, and I'd send it to peer review if it came across my desk, because the derivation is sound and the assumptions are transparent, but I'd ask the authors to add a cross-pair test or soften the conclusion.","headline":"A narrow, clean benchmark of a previously proposed functional-improvement scheme; the reported gains are in-sample and the power-law ansatz matches the target exactly, so the broader claim of general improvement is not supported.","tokens_in":6346,"tokens_out":4759,"would_cite":false,"duration_ms":42946,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A two-system fit using exact densities restores the LDA exchange functional from the Hartree baseline.","keywords":["density functional theory","inverse Kohn-Sham method","density functional perturbation theory","exchange-correlation functional","LDA exchange","noble-gas atoms","ground-state energy","functional improvement"],"falsifier":"Apply the same two-system fitting to a target functional that is not a single power law, for example Hartree plus LDA exchange plus LDA correlation; because Equation (10) cannot then be satisfied exactly, the residual error in the reproduced energy density across a wide density range would show whether the power-law ansatz is adequate.","tokens_in":5308,"feed_emoji":"⚛️","tokens_out":6881,"duration_ms":67582,"temperature":0.7,"pith_summary":"This paper tests a scheme, IKS-DFPT, for improving approximate energy density functionals by using the exact ground-state density. The idea is to treat the unknown correction to the Hartree-exchange-correlation functional as a small perturbation, compute the right-hand side of the first-order equation from known quantities, and fit a simple power-law ansatz to two systems. In benchmark calculations on noble-gas atoms, starting from the Hartree functional and targeting Hartree plus LDA exchange, the fitted exponent agrees with the target to within about 1 percent and the prefactor to within a few percent. The resulting IKS-DFPT energies are two to three orders of magnitude closer to the target, and densities one to two orders closer. The authors conclude that the method is promising for improving conventional functionals, including in nuclear density functional theory.","feed_headline":"Two noble-gas densities recover LDA exchange to ~1%","feed_subtitle":"Starting from the Hartree functional, the IKS-DFPT fit cuts ground-state energy errors by two to three orders of magnitude.","key_machinery":"The load-bearing object is the functional equation for the first-order correction. By comparing two first-order expressions for the exact ground-state energy, the paper obtains an equation whose right-hand side, $C[\\rho^{\\mathrm{exact}}_{\\mathrm{gs}}]$, depends only on known quantities: the inverse Kohn-Sham orbital energies, the known functional $\\tilde{E}_{\\mathrm{Hxc}}$, and the exact density. The ansatz $E^{(1)}_{\\mathrm{Hxc}}[\\rho] = A\\int\\rho(\\mathbf{r})^{\\alpha}\\,d\\mathbf{r}$ converts that functional equation into two algebraic equations for $\\lambda A$ and $\\alpha$ using two calibrating systems. The improved functional is then $E_{\\mathrm{Hxc}}[\\rho] = \\tilde{E}_{\\mathrm{Hxc}}[\\rho] + \\lambda E^{(1)}_{\\mathrm{Hxc}}[\\rho]$.","core_discovery":"The central claim is that the first-order correction functional $E^{(1)}_{\\mathrm{Hxc}}[\\rho]$ can be extracted without solving the full many-body problem. Equating two first-order expressions for the exact ground-state energy gives a functional equation whose right-hand side, $C[\\rho^{\\mathrm{exact}}_{\\mathrm{gs}}]$, is computable from the known functional, the exact density, and inverse Kohn-Sham orbital energies. With the ansatz $E^{(1)}_{\\mathrm{Hxc}}[\\rho] = A\\int \\rho(\\mathbf{r})^{\\alpha}\\,d\\mathbf{r}$, two systems fix the constants $\\lambda A$ and $\\alpha$. In the benchmark where the base functional is Hartree and the target is Hartree plus LDA exchange, the Ar-Kr pair yields $\\alpha = 1.3290958$ against a target of $1.3333333$, and $\\lambda A$ within 3.7 percent; across all noble-gas pairs $\\alpha$ is within about 1 percent and $\\lambda A$ within about 5 percent, with heavier pairs giving more accurate coefficients. The resulting functional improves ground-state energies by two to three orders of magnitude and densities by one to two orders of magnitude.","pith_inferences":["Because the benchmark target, LDA exchange, has exactly the same power-law form as the ansatz, the benchmark verifies the coefficient-extraction machinery rather than the generality of the ansatz; a target with different density dependence would be a harder test.","The persistent few-percent error in the prefactor across pairs suggests that two densities do not fully determine the correction; adding a third calibrating system or matching energy-density moments could tighten the fit.","In practical use, the exact density must come from experiment or high-accuracy methods, so the method's reliability will depend on the quality and range of that input density."],"forward_implications":["If the exact density of any two systems is known, the IKS-DFPT1 procedure yields the two constants in the corrected functional, so no direct minimization over functional space is needed.","For noble-gas benchmarks, the improved Hartree-based functional reproduces the LDA exchange target with the exponent within about 1 percent and the prefactor within a few percent, with heavier atoms giving more accurate fits.","The improved functional yields ground-state energies two to three orders of magnitude closer to the target and densities one to two orders closer than the original Hartree functional.","Because the improved functional is constructed to be system independent, the pair-specific calibration is meant to transfer to other electron systems.","The authors identify nuclear energy density functionals as a promising application of the same correction scheme."],"supporting_citations":[{"why":"Proposes the IKS-DFPT method whose benchmark this paper reports.","marker":"[5]"},{"why":"Defines the LDA exchange functional used as the target.","marker":"[13]"},{"why":"Provide the inverse Kohn-Sham method that yields exact orbital energies from a given density.","marker":"[6, 7]"},{"why":"Provide the density functional perturbation theory framework used for the first-order correction.","marker":"[8, 9, 10, 11]"},{"why":"Supplies the first-order perturbation theorem used in deriving the correction equation.","marker":"[12]"},{"why":"Establish density functional theory and the Kohn-Sham equations underlying the energy expressions.","marker":"[1, 2]"}],"fun_headline_variants":["IKS-DFPT: functional correction from two noble-gas densities","Ar and Kr fix exchange functional to within 1%","No full many-body: functional improved via IKS-DFPT","Two atoms, one functional: energy error drops 100x","Inverse KS meets DFPT to build better functionals"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The correction $E^{(1)}_{\\mathrm{Hxc}}[\\rho]$ is exactly $A\\int \\rho(\\mathbf{r})^{\\alpha}\\,d\\mathbf{r}$, with $A$ and $\\alpha$ fixed by two systems; if the true correction is not of this exact power-law form, the two fitted numbers cannot represent it.","fun_headline_variants_meta":{"raw":{"variants":["IKS-DFPT: functional correction from two noble-gas densities","Ar and Kr fix exchange functional to within 1%","No full many-body: functional improved via IKS-DFPT","Two atoms, one functional: energy error drops 100x","Inverse KS meets DFPT to build better functionals"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000838,"raw_usage":{"total_tokens":3606,"prompt_tokens":850,"completion_tokens":2756,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":466,"completion_tokens_details":{"reasoning_tokens":2670}},"tokens_in":466,"tokens_out":2756,"duration_ms":21127,"temperature":1.0,"reasoning_tokens":2670,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:23:07.826599+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Apply the same two-system fitting to a target functional that is not a single power law, for example Hartree plus LDA exchange plus LDA correlation; because Equation (10) cannot then be satisfied exactly, the residual error in the reproduced energy density across a wide density range would show whether the power-law ansatz is adequate.","supporting_citations":[{"cited_title":"Improvement of functionals in density functional theory by the inverse Kohn--Sham method and density functional perturbation theory","cited_arxiv_id":"1812.09285","evidence_quote":"Proposes the IKS-DFPT method whose benchmark this paper reports."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the LDA exchange functional used as the target."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the first-order perturbation theorem used in deriving the correction equation."}],"review_version":1}