{"id":"f78785ee-f326-450b-865c-9df1e093f0a0","arxiv_id":"2506.03578","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The LAK meta-GGA predicts the alpha and gamma phases of cerium with lattice constants close to experiment, unlike other semi-local functionals.","lead":"A new type of density functional calculation, the LAK meta-GGA, correctly predicts the two phases of cerium metal and the transition between them. If confirmed, this could make affordable computer simulations of strongly localized f-electron materials much more reliable.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central claim rests on the spin-polarized gamma solution; the paper never shows it survives independent initialization or magnetic ordering, so the double minimum could be a metastable artifact.","rationale":"The strongest claim would be true if LAK's spin-polarized solution at large volume is the correct zero-temperature representation of the local-moment gamma phase. This is exactly the least supported step in the paper. The double minimum appears only with a 3.0 μB spin-polarized initialization, and the gamma phase is experimentally paramagnetic with local moments, not a collinear ferromagnet. The text even contradicts itself by calling gamma-Ce ferromagnetic after citing Curie-Weiss paramagnetism. Without convergence against initialization and magnetic order, the gamma minimum could be a metastable broken-symmetry solution rather than the functional's true ground state. This is not an internal inconsistency in the LAK construction, but an external validity gap in the identification of the computed state with the physical phase. The concern is therefore load-bearing, but it is also checkable at modest cost. A positive outcome of the proposed test would substantially strengthen the paper; a negative outcome would reduce the result to a metastable-solution report. I agree with the reader's assessment, and the CONDITIONAL verdict remains appropriate, so no change is needed.","tokens_in":6422,"tokens_out":3452,"duration_ms":36743,"concrete_test":"Recompute the LAK energy-volume curve for bulk Ce at 520 eV cutoff and 24x24x24 k-mesh using: (a) spin-unpolarized initialization; (b) initial moments of 0.5, 1.0, 2.0, 3.0, 4.0 μB; and (c) a two-atom supercell with +3.0/−3.0 μB AFM initialization, reporting final total energy, magnetization, and lattice constants. Also run a noncollinear spin-spiral at the gamma volume. If the large-volume minimum disappears, shifts by more than ~0.05 eV, or is replaced by an AFM/noncollinear state, the claim that LAK predicts the correct gamma phase and ordering is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is the identification of the large-volume spin-polarized LAK state with the gamma phase. In Fig. 1 the double minimum exists only in the spin-polarized run, which was initialized with 3.0 μB; the spin-unpolarized curve has only the alpha minimum. The gamma phase is therefore entirely carried by a broken-symmetry magnetic solution. Yet the paper itself cites Curie-Weiss paramagnetism in gamma-Ce, and later calls gamma-Ce 'ferromagnetic in nature' (Sec. II, Fig. 2 discussion), a contradiction that is never resolved. No calculation is reported with other initial moments, with an antiferromagnetic two-atom cell, or with noncollinear moments, so we do not know whether the gamma minimum is the ground state of the LAK functional or a metastable artifact of the 3.0 μB initialization. If it is an artifact, the central claim reduces to 'LAK has an alpha phase only,' like OFR2 and r2SCAN. The magnetic moment (~1.2 μB) is also not compared with the local-moment magnitude expected for a Kondo/paramagnetic gamma-Ce, so the physical correspondence is unvalidated.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript applies the LAK meta-GGA to bulk cerium and reports, from a single self-consistent spin-polarized calculation, a double minimum in the energy-volume curve corresponding to the alpha and gamma phases, with lattice parameters 4.87 Å and 5.14 Å, a magnetic moment of about 1.2 μB, a critical pressure of -0.4 GPa, and a volume collapse of 15.8%. Two other meta-GGAs, OFR2 and r2SCAN, are shown to yield only the alpha phase. The paper argues that LAK's ultranonlocality is responsible for capturing localization/delocalization of f electrons and that this constitutes a parameter-free first-principles description of the isostructural alpha-gamma transition.","tokens_in":6606,"tokens_out":3775,"duration_ms":45366,"significance":"If the central claim holds, the result is significant: LAK would be the first semilocal-level functional without a Hubbard U or exact exchange to produce both cerium phases with the correct ordering and a realistic volume collapse, at a fraction of the cost of hybrid or RPA approaches. The study also provides a useful benchmark comparison of modern meta-GGAs on a correlated f-electron material. The paper's strengths are that it uses an external functional without fitting to cerium data, reports lattice constants close to experiment, and includes a computed volume collapse in good agreement with the measured value. However, the central claim depends entirely on a broken-symmetry spin-polarized solution whose stability and physical interpretation are not established, so the result is promising but not yet fully validated.","major_comments":[{"comment":"The double minimum appears only in the spin-polarized calculation initialized with an initial magnetic moment of 3.0 μB; the spin-unpolarized LAK curve has only the alpha minimum. Because the gamma phase is carried entirely by this particular broken-symmetry magnetic solution, the manuscript must demonstrate that this minimum is a stable stationary point of the LAK functional rather than a metastable artifact. Concretely, the authors should scan initial magnetic moments (for example 0.5, 1.0, 2.0, 3.0, 4.0 μB), test an antiferromagnetic two-atom cell, and test a noncollinear ordering. Without this evidence, the central claim reduces to the weaker statement that LAK has only an alpha phase unless artificially initialized into a ferromagnetic-like state.","section":"Figure 1 and the paragraph describing the spin-polarized LAK calculation"},{"comment":"The text states that the gamma phase is 'ferromagnetic in nature' and reports a magnetic moment of about 1.2 μB, but the Introduction correctly describes gamma-Ce as paramagnetic with Curie-Weiss susceptibility and localized moments. This contradiction is load-bearing because the physical identification of the calculated broken-symmetry state with the gamma phase depends on the local-moment picture. The paper should either provide evidence that the computed ferromagnetic-like solution represents the paramagnetic local-moment state (for example, by comparing with the expected Ce3+ effective moment or with the reduced moment expected from Kondo screening) or explicitly discuss the discrepancy.","section":"Figure 2 discussion and the magnetic-moment paragraph"},{"comment":"LAK's predicted critical pressure of -0.4 GPa is compared only with the value -0.8 GPa taken from Ref. [22], which is a theoretical study, while the experimental room-temperature values listed in Table II are +1.8 and +1.44 GPa. The sign of the predicted pressure is opposite to the room-temperature experimental sign, and the manuscript does not provide a temperature or pressure argument justifying the comparison to the extrapolated or theoretical value. This weakens the claim that LAK predicts the transition pressure 'reasonably close' to experiment; the authors should clarify what quantity is being compared and why a negative pressure is consistent with the measured phase transition.","section":"Table II and the critical-pressure discussion"}],"minor_comments":[{"comment":"There are duplicated words in the text, for example 'the the total energy' and 'the the double minimum'; these should be corrected.","section":"Abstract and text near Figure 1"},{"comment":"The projected density difference is defined as n_alpha - n_gamma, and the text says there is more charge localization in the gamma phase at the atomic sites. With a conventional color scale, this would make the atomic sites negative (blue) rather than deep red. The caption and color scale need to be clarified so the sign convention is unambiguous.","section":"Figure 3 caption and the corresponding paragraph"},{"comment":"The density of states figure is described as the sum of spin-up and spin-down states, but the figure does not show the spin-resolved components. Since the magnetic moment and the f-electron localization are central to the argument, showing the spin channels or at least reporting the f-occupancy in each phase would make the comparison more informative.","section":"Figure 4 discussion"},{"comment":"The same reference [22] is used for both the lattice parameters and the extrapolated pressure; the authors should explicitly state which entries in the tables are experimental and which are theoretical, since the text sometimes calls the -0.8 GPa value 'extrapolated to 0 K' without making clear that it is a theoretical estimate.","section":"Table II and Ref. [22]"}],"recommendation":"major_revision","confidential_remarks":"The paper has no circularity problem: LAK is an externally developed functional, and the only fitting is the equation-of-state fit. The main risk is purely technical but decisive: the gamma-phase minimum may be an initialization artifact of the spin-polarized calculation. I would encourage the editor to request the additional magnetic-order tests before considering publication, because the claim as stated is otherwise not robust."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper's one new result is clear: with the LAK meta-GGA, a spin-polarized self-consistent calculation gives a double minimum in cerium's energy-volume curve, with α and γ lattice constants close to experiment, where LDA, PBE, SCAN, r2SCAN, and OFR2 do not. That is worth taking seriously. The comparison with OFR2 and r2SCAN is useful, and the density-difference and DOS plots give some texture to the claim. The paper is also honest that the spin-unpolarized calculation yields only the α minimum.\n\nThe soft spot is the load-bearing one. The γ minimum exists only in the spin-polarized run, initialized with a 3.0 μB moment. The paper calls γ-Ce 'ferromagnetic in nature' in the figure discussion, while the introduction correctly states that γ is paramagnetic with Curie-Weiss local moments. That contradiction is never resolved. No other initial magnetization, no antiferromagnetic cell, no noncollinear calculation is reported. So we don't know whether the large-volume spin-polarized solution is the ground state of the LAK functional or a metastable artifact of the initialization. If it is an artifact, the central claim shrinks to 'LAK has only the α phase,' like the other semilocal functionals. The ~1.2 μB moment is also not compared with the expected local moment for γ-Ce.\n\nThere are smaller issues: the predicted critical pressure is -0.4 GPa, negative; the paper compares it to an extrapolated -0.8 GPa at 0K rather than the positive room-temperature experimental values, and the close agreement is a stretch. Zero-temperature lattice constants are compared directly to room-temperature experiment without comment. The language ('spectacular,' 'remarkable result') oversells.\n\nThese caveats are fixable in a revision. The core question—whether a semilocal functional without exact exchange can capture the α–γ transition—is important, and LAK's apparent success deserves a careful referee. But the magnetic solution needs to be validated before the claim can be accepted.","headline":"A potentially important meta-GGA result for cerium, but the γ phase rests on an unvalidated spin-polarized solution that could be a metastable artifact.","tokens_in":7154,"tokens_out":3215,"would_cite":false,"duration_ms":34079,"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 single DFT calculation reproduces both phases of cerium's volume collapse.","keywords":["cerium","alpha-gamma phase transition","meta-GGA","LAK functional","ultranonlocality","f-electron localization","volume collapse","density functional theory"],"falsifier":"Repeat the LAK calculation with initial magnetic moments of 0, 1, 2, and 4 μB, and with antiferromagnetic starting configurations. If the γ-phase minimum disappears or the energy ordering reverses in any of these runs, the claim that LAK captures the phase transition from a single self-consistent calculation fails. A second check is to compare the predicted ferromagnetic moment of ~1.2 μB with magnetic-susceptibility or neutron-scattering data that show γ-Ce has no long-range magnetic order.","tokens_in":6205,"feed_emoji":"⚛️","tokens_out":6662,"duration_ms":64423,"temperature":0.7,"pith_summary":"The paper claims that the LAK meta-GGA — a density functional that uses the orbital kinetic energy density to go beyond standard semilocal approximations — reproduces both the α and γ phases of cerium and their correct energetic ordering in a single spin-polarized calculation. This matters because cerium's isostructural α–γ transition, driven by the localization and delocalization of f electrons, has resisted ordinary semilocal functionals, which see only the α phase. The same calculation yields lattice parameters, a magnetic-moment jump, a critical pressure, and a volume collapse close to experiment, suggesting that a carefully constructed meta-GGA can handle localized f electrons without exact exchange or Hubbard corrections.","feed_headline":"Cerium's two phases from one self-consistent calculation","feed_subtitle":"The LAK meta-GGA gets the lattice constants, magnetic jump, and volume collapse right without exact exchange.","key_machinery":"The load-bearing object is the LAK meta-GGA. It is a semilocal functional built from the density, its gradient, and the orbital kinetic energy density τ, entering through the variable α = (τ − τ_W)/τ_unif that separates one-electron regions from the uniform electron gas. LAK adjusts the relative weights of gradient and τ contributions so that its exchange enhancement factor satisfies ∂F_x(s, α)/∂α < 0, the ultranonlocality condition that mimics some effects of exact exchange. The double minimum appears only in the spin-polarized solution: allowing a magnetic moment at large volume localizes the f electrons and stabilizes γ-Ce, so the mechanism is the combination of ultranonlocality with spin polarization.","core_discovery":"The central discovery is that the LAK meta-GGA produces a double minimum in the total energy as a function of lattice constant for bulk cerium, but only when the calculation is spin-polarized and starts from a magnetic moment of 3.0 μB. The minimum at larger volume is γ-Ce, with a magnetic moment of about 1.2 μB per atom and f electrons localized at atomic sites; the minimum at smaller volume is α-Ce, with zero moment and delocalized f electrons. From that single self-consistent calculation the authors obtain lattice constants of 5.14 Å and 4.87 Å for the γ and α phases, a critical pressure of -0.4 GPa, and a volume collapse of about 16%, all in reasonable agreement with experiment. The projected electron-density difference shows more charge buildup at atomic sites in the γ phase and more interstitial charge in the α phase, consistent with the localization picture.","pith_inferences":["I infer that the same mechanism may apply to other f-electron metals with volume-collapse transitions, such as plutonium or americium, where a spin-polarized LAK calculation could be tested for a double minimum.","The sensitivity to the initial magnetic moment suggests a systematic scan over starting magnetizations is needed to determine whether the γ minimum is the ground state or a metastable ferromagnetic state.","A natural extension is to compute phonons and finite-temperature free energies with LAK to verify the experimental finding that phonon entropy does not drive the transition.","The ultranonlocality of LAK may also improve charge-transfer and band-gap descriptions in other metals or Mott insulators, not just cerium; this could be tested on simple oxides."],"forward_implications":["LAK captures both cerium phases at zero temperature without exact exchange, Hubbard U, or random-phase approximation, unlike LDA, PBE, SCAN, r2SCAN, and OFR2.","The zero-temperature lattice constants from LAK (4.87 Å for α, 5.14 Å for γ) are close to room-temperature experimental values.","The predicted volume collapse of about 16% is much closer to the experimental 14–15% than the roughly 30% collapse from the hybrid-based (EX+cRPA)@PBE0 calculation.","The magnetic moment jumps from 0 in α-Ce to about 1.2 μB in γ-Ce, marking the localization of f electrons at the transition.","If these results hold, ultranonlocal meta-GGAs could serve as a low-cost alternative to fourth- and fifth-rung methods for f-electron materials."],"supporting_citations":[{"why":"Defines the LAK meta-GGA, the functional whose enhanced nonlocality is the central object of the paper.","marker":"[12]"},{"why":"Previous (EX+cRPA)@PBE0 calculation that found a double minimum for cerium and supplies the main comparison for lattice parameters, critical pressure, and volume collapse.","marker":"[4]"},{"why":"Experimental lattice constants for α- and γ-cerium used to benchmark the LAK results.","marker":"[23]"},{"why":"Experimental critical pressure and volume collapse used to benchmark the LAK predictions.","marker":"[25]"},{"why":"Theoretical reference for the critical pressure extrapolated to zero temperature, against which LAK's -0.4 GPa is compared.","marker":"[22]"},{"why":"Self-interaction-corrected LDA calculation of cerium whose density of states is used as a comparison and whose earlier attempt at the transition is referenced.","marker":"[5]"}],"fun_headline_variants":["One self-consistent run gives both cerium phases","LAK meta-GGA reproduces cerium's α–γ transition in one shot","Cerium's two phases from one LAK calculation","LAK delivers both cerium phases without exact exchange"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The γ-phase minimum is obtained from a spin-polarized calculation that assumes an initial magnetic moment of 3.0 μB, and the paper treats this ferromagnetic-like solution as physical even though γ-Ce is experimentally paramagnetic with local moments.","fun_headline_variants_meta":{"raw":{"variants":["One self-consistent run gives both cerium phases","LAK meta-GGA reproduces cerium's α–γ transition in one shot","Cerium's two phases from one LAK calculation","LAK delivers both cerium phases without exact exchange"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000619,"raw_usage":{"total_tokens":2909,"prompt_tokens":1022,"completion_tokens":1887,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":638,"completion_tokens_details":{"reasoning_tokens":1826}},"tokens_in":638,"tokens_out":1887,"duration_ms":15833,"temperature":1.0,"reasoning_tokens":1826,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T10:59:02.335936+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Repeat the LAK calculation with initial magnetic moments of 0, 1, 2, and 4 μB, and with antiferromagnetic starting configurations. If the γ-phase minimum disappears or the energy ordering reverses in any of these runs, the claim that LAK captures the phase transition from a single self-consistent calculation fails. A second check is to compare the predicted ferromagnetic moment of ~1.2 μB with magnetic-susceptibility or neutron-scattering data that show γ-Ce has no long-range magnetic order.","supporting_citations":[{"cited_title":"Lebeda, T","cited_arxiv_id":null,"evidence_quote":"Defines the LAK meta-GGA, the functional whose enhanced nonlocality is the central object of the paper."},{"cited_title":"Casadei, X","cited_arxiv_id":null,"evidence_quote":"Previous (EX+cRPA)@PBE0 calculation that found a double minimum for cerium and supplies the main comparison for lattice parameters, critical pressure, and volume collapse."},{"cited_title":"Olsen, L","cited_arxiv_id":null,"evidence_quote":"Experimental lattice constants for α- and γ-cerium used to benchmark the LAK results."},{"cited_title":"Decremps, L","cited_arxiv_id":null,"evidence_quote":"Experimental critical pressure and volume collapse used to benchmark the LAK predictions."},{"cited_title":"Casadei, X","cited_arxiv_id":null,"evidence_quote":"Theoretical reference for the critical pressure extrapolated to zero temperature, against which LAK's -0.4 GPa is compared."},{"cited_title":"Szotek, W","cited_arxiv_id":null,"evidence_quote":"Self-interaction-corrected LDA calculation of cerium whose density of states is used as a comparison and whose earlier attempt at the transition is referenced."}],"review_version":1}