{"id":"ac569c15-3b50-4b52-920b-b06460553477","arxiv_id":"2501.04278","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"The authors benchmark LDA, LDA+U, QSGW, and DFT+HI methods for RECo5 magnets and propose a penalty functional to make Hund's rule-constrained DFT total energies variational.","lead":"This paper compares several computer methods that calculate the electronic structure and magnetic anisotropy of rare-earth cobalt magnets (RECo5). It suggests a small technical improvement to enforce Hund's rule in density functional calculations, aiming to make such calculations reliable enough for screening new permanent magnet materials.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The accuracy claim rests on a qualitative CF sign trend that the Hund's-rule constraint itself imposes, with no quantitative benchmark against experimental anisotropy or CF parameters.","rationale":"The reader's weakest assumption identifies 4f localization as the key premise for open-core/Hubbard-I treatments. That premise is necessary, but even if localization is granted, the present manuscript does not provide evidence that the computed anisotropy is quantitatively accurate. The HRC penalty explicitly fixes the 4f spin and orbital moments and their orientation, so the oblate/prolate density shape responsible for the CF anisotropy is inserted by hand. The paper's own statement that the sign variation follows the Hund's-rule shape of the tripositive RE ion confirms that Figs. 3(b) and 4 primarily reproduce an imposed trend. A variational penalty functional improves the constrained-energy formalism but cannot by itself establish that the constrained state is the true ground state. The strongest claim in Sec. IV is therefore not tested by anything in the paper: no experimental anisotropy energy or CF parameter is compared with the computed number. The proposed check, benchmarking computed anisotropy magnitudes against measured values for several RECo5 compounds while verifying λ-independence, would settle whether the method is genuinely accurate or merely restates the input CF model. Since this concern is closely related to but more specific than the reader's localization-based concern, agreement is partial. The verdict should remain CONDITIONAL because the concern is serious but the central claim is not yet disproved; a quantitative benchmark could support it, so a conditional accept with this test as a condition is appropriate.","tokens_in":6482,"tokens_out":6458,"duration_ms":69061,"concrete_test":"Compute the X/Z anisotropy energy for SmCo5, PrCo5, NdCo5, TbCo5, and ErCo5 with the HRC functional, report the numerical values in meV per RE site, and compare against experimental anisotropy constants or anisotropy energies from magnetometry and neutron data. Require both correct sign and agreement within about 30% for all non-Ce compounds. In the same runs, sweep the penalty strength λ over at least two orders of magnitude to verify that the X−Z energy difference is converged and independent of λ; if the difference changes with λ or fails to match experiment, the Sec. IV accuracy claim is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section III introduces the Hund's-rule constraint (HRC) penalty functional E = E0 + Σ_i λ_i (M_i − e_i M_i)^2 and then reports the X/Z magnetic-anisotropy differences shown in Fig. 4. The central problem is that the constraint prescribes both the spin and orbital 4f moments and their spatial direction e_i for each field orientation; hence the oriented, oblate/prolate 4f charge cloud that produces the CF anisotropy is an input, not an output of the calculation. The authors acknowledge this in Sec. III: the sign trend 'follows oblate or prolate shapes of 4f electron cloud ... in their ground states of the Hund's rule.' Consequently, seeing the expected CF-like sign pattern in Figs. 3(b) and 4 validates the constraint implementation but not the accuracy of the predicted anisotropy. The paper provides no quantitative comparison of computed anisotropy energies or CF splittings with experimental values for any RECo5 compound; the cited Brooks et al. agreement is from Ref. 18 and is not reproduced here. The new penalty restores variationality of E0+Epen, but variationality alone does not show that the constrained state is the physical ground state. Without a benchmark, the Sec. IV conclusion that this combination is 'the most flexible and accurate approach' is unsupported. The method may be a useful re-derivation of CF results, but the accuracy and predictive-power claim is the load-bearing assertion that has not been tested.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper surveys electronic-structure approaches for RECo5 permanent-magnet materials, focusing on the 4f moments and magnetic anisotropy of the rare-earth sublattice. It compares open-core LDA, QSGW, DFT+HI, and a constrained DFT scheme in which a penalty functional enforces Hund's-rule spin and orbital moments on the 4f site. The authors propose that adding such a penalty functional restores variationality to the total energy, and they decompose the SmCo5 anisotropy using crystal-field parameters obtained from DFT+HI. The central claim, stated in the Conclusion, is that the proposed combination of methods is 'the most flexible and accurate approach for material science progress in this area.'","tokens_in":6799,"tokens_out":2246,"duration_ms":23608,"significance":"If the central claim were substantiated, the work would be useful for high-throughput screening of rare-earth permanent magnets: it would provide a practical way to obtain 4f moments and crystal-field anisotropies with variational total energies. The paper has several genuine strengths: it explicitly compares QSGW and LDA densities of states for SmCo5, it proposes a physically motivated penalty functional for Hund's-rule constraints, and it connects constrained DFT results to the standard crystal-field picture. However, the significance is limited by the absence of quantitative benchmarks: no computed anisotropy energy or crystal-field splitting is compared with experimental values for any RECo5 compound, and the load-bearing accuracy claim is therefore unsupported. The paper is better read as a methodological scoping study than as a validated predictive method.","major_comments":[{"comment":"The penalty functional is written as E = E0 + Σ_i λ_i (M_i − e_i M_i)^2, but no Lagrange multiplier values, convergence criteria, or demonstration that the constraints are actually satisfied are given. Without this information the reader cannot assess whether the constrained solutions are the intended Hund's-rule states or whether the penalty introduces large uncontrolled energy shifts. This is load-bearing because the anisotropy differences in Fig. 4 are computed from total-energy differences that include the penalty term.","section":"Section III, Eq. (1)"},{"comment":"The Hund's-rule constraint prescribes both the target spin moment M_s and orbital moment M_l and the spatial direction e_i for each field orientation. The resulting oriented, oblate/prolate 4f charge density is therefore an input to the calculation, not an independent prediction. The authors themselves state in Section III that the sign trend 'follows oblate or prolate shapes of 4f electron cloud ... in their ground states of the Hund's rule.' Consequently, agreement of Fig. 4 with the crystal-field sign pattern validates the constraint implementation but does not validate the magnitude or predictive accuracy of the computed anisotropy.","section":"Section III, Figs. 3 and 4"},{"comment":"The paper reports no quantitative comparison of computed anisotropy energies with experimental magnetocrystalline anisotropy for any RECo5 compound, nor a comparison of computed crystal-field splittings with experimental inelastic neutron scattering or specific-heat data. The only quantitative agreement cited is from Ref. 18 (Brooks et al.), which is not reproduced here. The Conclusion's claim that the proposed combination is 'the most flexible and accurate approach' is therefore unsupported by the evidence presented in this manuscript.","section":"Section III, Fig. 4 and Conclusion"},{"comment":"The LDA+U calculation uses Hubbard U = 6.7 eV, but no justification for this value, no U-dependence study, and no comparison with open-core (U = ∞) results is provided. Since Fig. 4 reports both 'open core' and LDA+U results, the difference between the two sets of points depends on an ad hoc parameter. A sensitivity scan or a derivation of U from first principles is needed before the magnitude of the predicted anisotropy can be trusted.","section":"Section III, Fig. 4"},{"comment":"The SmCo5 anisotropy decomposition in Fig. 5 relies on crystal-field and exchange-field parameters taken from Table VI of Ref. 23, obtained by a non-spin-polarized DFT+HI calculation. The transferability of those parameters to the single-site Hamiltonian used here is asserted but not tested, and no comparison with the constrained-DFT anisotropy of Fig. 4 is made. The text would need either a direct reproduction of the Ref. 23 calculation or an explicit error estimate for the parameter transfer before the 15% spin-orbit contribution and the 10% non-CF contribution can be considered quantitative.","section":"Section II and Fig. 5"}],"minor_comments":[{"comment":"The equation for the penalty functional is garbled in the typeset text ('𝐸=𝐸!+𝐸\"=𝐸!+%𝜆##'(𝑴#(−𝒆#𝑴#'); please restore proper math notation and define all symbols, including what e_i multiplies.","section":"Section III, Eq. (1)"},{"comment":"The phrase 'exited properties' should be 'excited properties'.","section":"Section IV, first paragraph"},{"comment":"The statement that 'nonspherical components of the 4f charge density were ignored' is clear, but the subsequent sentence about reducing the effective potential to spherical form 'only when evaluating the CF energies' is confusing; please clarify which potential is used in which calculation step.","section":"Section II"},{"comment":"The caption reports U = 6.7 eV but does not state whether the same U is used for all rare-earth atoms or whether the open-core results are U = ∞ in the same code; please specify.","section":"Fig. 4 caption"},{"comment":"Fig. 1 shows the hierarchy of spin-orbit, exchange, and crystal-field energies, but the manuscript does not state how these values were computed; please add a sentence describing the method or cite the source.","section":"Section I, Fig. 1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a useful comparison of existing methods, but it overclaims accuracy without quantitative benchmarks. The referee's main concern is that the constrained-DFT anisotropy is largely determined by the imposed Hund's-rule state, so the qualitative agreement with crystal-field theory is not a validation of predictive power. The authors should either add experimental comparisons (e.g., anisotropy fields or CF splittings) or substantially soften the Conclusion. The paper also needs to provide the missing technical details of the penalty functional and the U-dependence analysis. With those additions, the manuscript could become a solid methodological reference; in its present form it is not ready for publication in a top-tier journal."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read it. The useful part is methodological: they implement a penalty functional after Ma-Dudarev to make Hund's-rule-constrained DFT total energies variational, and they compare open-core LDA, LDA+U, QSGW, and DFT+HI on the same RECo5 family. The QSGW density of states for Sm is a nice sanity check that the 4f manifold is atomic-like. The paper is also honest: it says HRC should not be enforced for Ce, that DMFT total energies are not reliable enough for anisotropy, and that the sign trend follows known CF theory.\n\nThe soft spot is the one the stress test flags. The anisotropy sign pattern in Fig. 4 is substantially a consequence of the imposed Hund's rule state. Constraining both spin and orbital moments along a chosen field direction fixes the oriented 4f charge cloud, so recovering the oblate/prolate sign rule validates the constraint implementation, not the predictive accuracy of the anisotropy. There is no quantitative benchmark of computed anisotropy energies or CF splittings against experiment for any RECo5 compound; the Brooks et al. agreement is cited but not reproduced. The penalty functional is only sketched—no Lagrange multiplier values, no convergence details—and U = 6.7 eV enters as an ad hoc input. Data availability says the data are in the article, but no raw values or code are provided.\n\nThese are real gaps, but they are not fatal for the paper's narrower purpose. As a comparison of available approaches and a proposal for a variational constraint, it is useful for the RE-magnet computational community. What does not hold is the conclusion that the combination is 'the most flexible and accurate approach' in the field. That claim needs a benchmark, or at least a direct comparison with measured anisotropy constants. The self-acknowledged limitations soften the issue: the authors themselves note limited applicability and uncontrolled electronic structure for HRC methods.\n\nI would send it to peer review, because the methodological point is worth airing and the comparison data are not available elsewhere. But I would insist the referee asks for a quantitative benchmark, a fuller specification of the penalty parameters, and error estimates before publication. If the authors add those, it becomes a solid contribution; in its current form it is a reasonable preprint with an overbroad abstract.\n\nWho this is for: computationalists working on rare-earth permanent magnets or constrained DFT. I would not bring it to a general reading group, but if the group does correlated-electron methods, it is worth one slot.","headline":"A useful methods comparison for RECo5 that restores variationality to Hund's-rule-constrained DFT, but the headline accuracy claim is not yet backed by a quantitative benchmark.","tokens_in":7286,"tokens_out":2757,"would_cite":false,"duration_ms":28809,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["71.15.Mb","75.30.Gw","71.70.Ch","75.50.Ww"],"model":"deepseek-v4-flash","headline":"The paper argues that reliable RECo5 electronic-structure calculations require treating 4f electrons as localized atomic shells, adding a Hubbard U, removing the anisotropic part of the 4f charge density, and finishing with crystal-field…","keywords":["rare-earth permanent magnets","RECo5","crystal-field theory","4f electron localization","magnetic anisotropy","Hubbard-I approximation","constrained DFT","Hund's rules"],"falsifier":"Measure the occupied 4f bandwidth in a RECo5 compound with angle-resolved photoemission: if the 4f states disperse by an energy comparable to the crystal-field splittings rather than forming flat atomic levels, or if inelastic neutron scattering resolves crystal-field levels that disagree with the DFT+HI spherical-average prediction by more than experimental error, the localization premise fails.","tokens_in":6251,"feed_emoji":"🧲","tokens_out":9985,"duration_ms":92464,"temperature":0.7,"pith_summary":"This paper seeks to establish a reliable computational recipe for the electronic structure and magnetic anisotropy of RECo5 permanent magnets, where rare-earth 4f electrons are the source of the anisotropy. The authors argue that standard DFT-based treatments fail because 4f levels sit at the Fermi energy, hybridize too strongly, and develop a spurious self-interaction from the anisotropic part of the 4f charge density, which inflates crystal-field splittings by 10 to 30 times. Their recipe combines a localized treatment of 4f electrons (open-core or Hubbard-I), a Hubbard U that reproduces the experimental spin splitting, and removal of the nonspherical 4f charge density, followed by a crystal-field-theory step to obtain anisotropy. They also introduce a penalty functional that makes Hund's-rule-constrained DFT total energies variational. If this recipe is right, it makes high-throughput searches for new rare-earth-containing magnets possible because the difficult secondary magnetic properties become computable without fitting parameters.","feed_headline":"Rare-earth magnet anisotropy now computable with atomic 4f shells","feed_subtitle":"Treating 4f electrons as localized and stripping their nonspherical density reproduces crystal fields and anisotropies.","key_machinery":"The load-bearing object is the anisotropic 4f charge density and its removal. In the constrained-DFT branch, the 4f shell is treated as an open core whose nonspherical charge components are excluded from the self-consistent density and potential, eliminating the self-interaction that otherwise exaggerates crystal-field splittings by an order of magnitude. To make Hund's-rule-constrained calculations variational, the paper adds a penalty functional to the DFT energy that penalizes deviation of each site's spin and orbital moments from the wanted Hund's values. In the DFT+HI branch, localized 4f electrons are solved with a Hubbard-I atomistic impurity solver and the f orbitals are represented by extended Wannier functions formed from a narrow f-band window, which keeps the hybridization contribution to the crystal field while preventing artificial delocalization. Crystal-field theory then converts the computed CF and exchange-field parameters into magnetic anisotropy energies.","core_discovery":"The authors' central claim is that no single existing method is enough: LDA puts 4f states at the Fermi level, LDA+U keeps them too hybridized, open-core and Hubbard-I localize them but still see a wrong crystal field unless the anisotropic 4f charge density is removed, and QSGW produces an atomistic 4f structure but lacks a total energy accurate enough for anisotropy. The paper concludes that the correct description must include a localized 4f shell (open core or Hubbard I), a Hubbard U chosen to match experimental spin splittings, and elimination of the anisotropic part of the 4f charge density; when a reliable total energy is unavailable, crystal-field theory with Wannier-function f orbitals supplies the required CF splittings and anisotropies. Hund's-rule-constrained DFT is made variational by adding a penalty functional that drives spin and orbital moments to their Hund's values, and such constraints should not be enforced for Ce systems, whose 4f electrons are not sufficiently localized. In SmCo5 the calculation decomposes the 4f anisotropy into crystal-field, spin-orbit, and other terms, finding the crystal field dominant with spin-orbit contributing roughly 15 percent with the opposite sign.","pith_inferences":["A testable consequence the paper does not spell out: compounds with the same rare-earth ground multiplet should show crystal-field splittings that scale in a predictable way with the lattice environment, since the spherical-averaging step makes the 4f density a rigid shape that samples the ligand potential.","The penalty-functional construction could be extended to map total-energy surfaces across coupling schemes such as LS, jj, and Jj, turning the method into a diagnostic for how well Hund's rules hold in a given compound.","A high-throughput screen using this pipeline could rank substitutions such as Fe or Ni on the Co sites by their predicted anisotropy before synthesis, which is the concrete route to new permanent magnets the paper points toward but does not carry out.","QSGW could serve as a localization check: where QSGW 4f densities of states lose their atomic multiplet structure, the open-core or Hubbard-I assumption should be abandoned, as in CeCo5."],"forward_implications":["For any RECo5 compound whose 4f electrons are localized, the combination yields 4f spin and orbital moments plus crystal-field anisotropy without adjustable parameters.","The same pipeline can be transplanted to other rare-earth-transition-metal magnets, where the localized-4f and crystal-field steps are independent of the transition-metal sublattice.","Ce-based magnets must be treated differently: Hund's-rule constraints, open-core, and Hubbard-I treatments all lose validity when 4f states are not localized.","With the penalty functional, constrained Hund's-rule calculations return a variational total energy, so total-energy comparisons of magnetic anisotropy become meaningful.","In SmCo5 the crystal-field contribution dominates 4f anisotropy, with spin-orbit anisotropy about 15 percent of it and opposite in sign, making the small non-CF terms a controlled approximation."],"supporting_citations":[{"why":"Supplies the constrained DFT total-energy method that removes the nonspherical 4f self-interaction and yields crystal-field energies.","marker":"[18]"},{"why":"Establishes the DFT+HI method with spherical averaging and extended Wannier f orbitals that produces crystal-field parameters in good agreement with experiment.","marker":"[11]"},{"why":"Provides the higher-order crystal-field and exchange-field parameters for SmCo5 used in the paper's anisotropy decomposition.","marker":"[23]"},{"why":"Introduced spherical averaging of the f-charge in LDA+U to remove the parasitic 4f self-interaction.","marker":"[6]"},{"why":"The Hubbard-I standard model of rare earths that DFT+HI builds on for treating 4f electrons as localized.","marker":"[13]"},{"why":"Supplies the penalty-functional recipe in constrained DFT that the paper adapts for Hund's-rule constraints.","marker":"[29]"},{"why":"Describes the QSGW method used to show that 4f states preserve an atomistic multiplet structure.","marker":"[27]"},{"why":"A prior application of orbital-polarization-type constraint to SmCo5 that motivates the total-energy-corrected treatment.","marker":"[31]"}],"fun_headline_variants":["Localized 4f shells plus charge stripping compute RECo5 anisotropy","Rare-earth magnet anisotropy from localized 4f electrons","Penalty functional makes Hund's rule DFT variational for rare-earth magnets","Combined DFT methods compute RECo5 crystal fields and anisotropy"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The method's success depends on the 4f electrons staying tightly bound to the rare-earth atom, so that treating them as a rigid atomic shell and discarding the non-spherical part of their charge still reproduces the crystal field the atom actually feels.","fun_headline_variants_meta":{"raw":{"variants":["Localized 4f shells plus charge stripping compute RECo5 anisotropy","Rare-earth magnet anisotropy from localized 4f electrons","Penalty functional makes Hund's rule DFT variational for rare-earth magnets","Combined DFT methods compute RECo5 crystal fields and anisotropy"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000557,"raw_usage":{"total_tokens":2639,"prompt_tokens":920,"completion_tokens":1719,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":536,"completion_tokens_details":{"reasoning_tokens":1646}},"tokens_in":536,"tokens_out":1719,"duration_ms":12874,"temperature":1.0,"reasoning_tokens":1646,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T21:37:56.518564+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the occupied 4f bandwidth in a RECo5 compound with angle-resolved photoemission: if the 4f states disperse by an energy comparable to the crystal-field splittings rather than forming flat atomic levels, or if inelastic neutron scattering resolves crystal-field levels that disagree with the DFT+HI spherical-average prediction by more than experimental error, the localization premise fails.","supporting_citations":[],"review_version":1}