{"id":"2955130b-264d-4747-903e-c45d3e441da1","arxiv_id":"2508.19496","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A multipole-ratio correction fixes DFT's overestimation of light rare-earth 4f charge asphericity, bringing the calculated magnetic anisotropy energy of SmCo5 into the experimental range.","lead":"Rare-earth magnets are hard to model with standard computer simulations because a single electron wavefunction distorts the 4f electron cloud too much. The authors add a many-body correction based on the true quantum state of the rare-earth ion, and their corrected calculation matches experiment for SmCo5.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"SmCo5 agreement may be partly fortuitous because Eq. (4) scales the total fitted anisotropy, including Co-sublattice and exchange terms, not just the Sm 4f multipoles.","rationale":"The central claim has two parts: (i) single-Slater-determinant DFT overestimates the 4f charge asphericity for light REs, and (ii) the MAE error can be removed by multiplying the DFT anisotropy by ion-only ratios Q_true/Q_DFT. Part (i) is well supported by the angular-momentum argument and Table I. Part (ii) is the delicate step. In a TbV6Sn6-like case with a nonmagnetic V sublattice, the DFT total-energy anisotropy is essentially the RE single-ion contribution, so fitting Eq. (3) and later scaling κ_l may be legitimate. But in RCo5 the total E(θ,φ) also contains Co sublattice and exchange anisotropy. Since the fit to Eq. (3) is a symmetry-adapted expansion, the fitted κ_l absorb all such terms; Eq. (4) then scales all of them by 4f angular-momentum ratios. This is not a controlled test of the correction. The reader's weakest assumption about configuration independence of A_q^k is related, but the failure identified here can occur even if all A_q^k are perfectly configuration independent: the non-RE part of κ_l is simply not proportional to Q_l of the RE ion. The proposed YCo5 subtraction experiment is decisive because it isolates the non-4f anisotropy that should not be multiplied by Table I ratios. If the corrected SmCo5 MAE is stable to that subtraction, the central claim survives; if not, the apparent agreement is partly accidental and the ion-only correction needs reformulation. I therefore do not change the reader's CONDITIONAL verdict: the SmCo5 demonstration is promising but not yet a rigorous test of the broad generality claim.","tokens_in":9301,"tokens_out":12364,"duration_ms":119490,"concrete_test":"Compute the same constrained DFT+U E(θ,φ) for YCo5 and fit Eq. (3) to obtain the non-4f κ_l^m; then recompute the SmCo5 corrected MAE as (κ_l^Sm - κ_l^YCo5)·(Q_true/Q_DFT)_l + κ_l^YCo5. If the result shifts by more than about 1-2 meV/f.u. or leaves the 13-16 meV/f.u. range, the total-energy scaling in Eq. (4) is not valid and the SmCo5 agreement is not a clean confirmation of the ion-only correction.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing step is the application of Eq. (4) to the κ_l^m coefficients obtained by fitting the total DFT energy surface E(θ,φ) to Eq. (3). In RCo5 the fitted κ_l^m are not pure single-ion 4f crystal-field terms: they absorb the magnetocrystalline anisotropy of the Co sublattice, 4f-3d exchange anisotropy, and hybridization contributions, all of which have the same P_l(cosθ) angular forms. Eq. (4) is derived from angular-momentum algebra of the 4f shell alone (Table I), so it should multiply only the 4f part of κ_l^m. Scaling the entire fitted κ_l^m by the Sm ratios (13/21 for k=2, 13/84 for k=4, 0 for k=6) also scales those non-4f contributions. The paper does not quantify the non-Sm part, for example by computing the same anisotropy for YCo5. The validations on TbV6Sn6 and TbCo5 are performed for a heavy RE where |JJ> = |LL> and no many-body correction is needed, so they do not test the correction procedure in a system with a magnetic Co sublattice. If the Co contribution is non-negligible, the reported 13-16 meV/f.u. agreement for SmCo5 could be the result of a partial cancellation rather than a controlled test of the ion-only correction, directly threatening the paper's generality claim.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript addresses the long-standing failure of DFT-based methods to reproduce the magnetocrystalline anisotropy (MA) of light rare-earth (RE) compounds. The authors argue that the single Slater-determinant representation in DFT overestimates the 4f charge asphericity, leading to an overestimated MA. They propose a correction that multiplies the crystal-field anisotropy coefficients extracted from constrained DFT+U total-energy surfaces by ratios of multipole moments computed from angular momentum algebra for the true |J,J> and the DFT |L,L> states. They validate the crystal-field-parameter extraction by showing consistency across different 4f configurations in TbV6Sn6 and TbCo5, and apply the correction to SmCo5, obtaining a corrected MAE in the experimental range of 13-16 meV/f.u.","tokens_in":9612,"tokens_out":11783,"duration_ms":100004,"significance":"The central idea is elegant and the angular-momentum algebra leading to the correction factors in Table I is on firm ground. If the approach is valid, it offers a practical, efficient DFT-based route to light-RE anisotropy, which is important for permanent-magnet design. The internal consistency check across m_l configurations is a valuable test of the assumption that DFT total-energy surfaces can be described by a crystal-field model. However, the only external validation of the many-body correction is SmCo5, where the correction is applied to the total anisotropy without separating the Co-sublattice contribution. The paper's conclusions are therefore plausible but not yet fully established.","major_comments":[{"comment":"The correction ratios in Table I are derived from the angular momentum algebra of the 4f shell alone, yet they are applied to the κ_l^m coefficients obtained by fitting the total DFT energy E(θ,φ) of SmCo5 to Eq. (3). In RCo5, the total anisotropy also contains Co-sublattice magnetocrystalline anisotropy and 4f-3d exchange anisotropy, which have the same P_l(cosθ) angular forms. The manuscript does not separate these contributions, so scaling the entire fitted κ's by the Sm ratios in Table I also scales the non-4f terms. The validation on TbCo5 does not test the many-body correction because Tb is a heavy RE with |JJ>=|LL>, so no correction is applied. To support the SmCo5 result, the authors should either compute the corresponding anisotropy for YCo5 and subtract it from the SmCo5 total before applying the correction, or otherwise demonstrate that the Co contribution is negligible. Without this decomposition, the reported agreement with the experimental MAE of 13-16 meV/f.u. could be fortuitous.","section":"Eq. (3), Eq. (4), and the paragraph 'Now we use these corrections'"},{"comment":"The correction formula assumes that the DFT constrained 4f state can be identified with the maximally polarized |L,L> Slater determinant for the purpose of computing Q_DFT_k. This assumption is not directly tested for light RE ions. The consistency check across |m_l> states is performed for Tb (a heavy RE) and does not establish that the Sm 4f state in the DFT calculation is exactly |L,L>; hybridization or self-interaction errors could distort the occupation matrix away from this ideal configuration. The authors should provide a quantitative comparison of the DFT 4f density matrix with the ideal |L,L> projection, or validate the correction on a light-RE compound with a non-magnetic transition-metal sublattice, where the 4f contribution can be isolated.","section":"Eq. (4) and the paragraph 'Our procedure relies on the consistency of CF theory within DFT'"},{"comment":"The SmCo5 result is reported for a single value of U (10 eV) and without uncertainty estimates. Since the extracted κ_l^m depend on U and other computational parameters while the correction factors in Table I do not, a sensitivity analysis is needed to show that the agreement with experiment is robust rather than parameter-dependent. The Supplemental Material is cited for U-dependence but is not available in the preprint.","section":"SmCo5 application paragraph and Supplemental Material note"}],"minor_comments":[{"comment":"The abstract states that the method is 'confirmed' on TbV6Sn6 and TbCo5, but for these heavy-RE compounds only the crystal-field-parameter extraction is validated; the many-body correction is not applied. Please phrase this distinction explicitly.","section":"Abstract"},{"comment":"The spread of the extracted CFPs across m_l states is described only qualitatively ('within ~14%'). Providing numerical values or error bars would strengthen the consistency claim.","section":"Figure 2(e,f)"},{"comment":"The symbols θ_k, O0_k, and ⟨r^k⟩ are not fully defined in the text. Please define them and state that the ratios in Table I are independent of ⟨r^k⟩.","section":"Eq. (4)"},{"comment":"The paper refers to the Supplemental Material for U-dependent results, but the preprint does not include it. Please ensure the supplemental file is provided with the submission.","section":"Supplemental Material"},{"comment":"The statement that the overestimation 'generally increases with both atomic number and multipole order k' is not strictly monotonic in Z for k=2 (e.g., Pr has a larger Q2/Q_DFT_2 ratio than Ce). Consider rephrasing to avoid overstatement.","section":"Figure 1 caption"},{"comment":"The limitations paragraph mentions that DFT may probe bare CFPs that could be renormalized by hybridization and interactions. This is an important caveat for the interpretation of the 'close agreement with experiment' and should be discussed in more detail.","section":"Limitations paragraph"},{"comment":"The constrained-DFT protocol is cited to Ref. [12], but a brief description of the constraint (e.g., the occupation matrix setting) in the text would make the paper more self-contained.","section":"Methods description"}],"recommendation":"major_revision","confidential_remarks":"The paper fits the scope of the journal. The main concern is the SmCo5 validation; the authors should be encouraged to provide a YCo5 calculation or equivalent decomposition of the anisotropy into RE and TM contributions. Also, the preprint lacks the Supplemental Material referenced for U-dependence; it should be provided to the reviewers."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe headline: this paper has a genuinely good idea—DFT's single Slater determinant overestimates the 4f charge asphericity for light rare-earths, and an ion-dependent multipole ratio can fix that. The SmCo5 result falls in the experimental range, which is encouraging. But there's a soft spot the paper doesn't address: the correction is applied to the total fitted anisotropy, which in RCo5 includes Co-sublattice and exchange terms, not just the Sm 4f crystal field. So the agreement could be partly accidental.\n\nWhat's new: the multipole-ratio table (Table I) is clean angular momentum algebra. The consistency check across different m_l states in TbV6Sn6 is a nice validation that CF theory works inside DFT—nontrivial and worth credit. The paper also openly lists limitations (integer valence, Ce-Hill, actinides), which I appreciate.\n\nWhat's soft: first, the Co contamination issue. If the fitted κ_l^m absorb Co anisotropy, multiplying the whole thing by Q_l^true/Q_l^DFT scales those non-4f parts too. A YCo5 calculation would quantify this. Without it, the SmCo5 agreement is a single point with an unquantified systematic risk. Second, configurational independence of the CFPs is only tested on Tb (heavy RE), where no correction is needed; for Sm it's an assumption. Third, no error bars, no U-dependence scan, and the supplemental details are missing. These are fixable.\n\nThe core idea is solid. The multipole algebra is exact; the question is whether the DFT fit isolates the 4f part. That's an empirical question the paper doesn't fully nail.\n\nFor who: anyone doing DFT for RE permanent magnets will want this table and the method. I'd bring it to group meeting.\n\nRecommendation: send it to peer review. A good referee will ask for the YCo5 control, a clearer statement of what the correction multiplies, and the missing supplement. The paper deserves that.","headline":"A physically motivated correction for light-RE MAE in DFT, with a promising SmCo5 result that is less controlled than it looks because the correction scales the whole fitted anisotropy, not just the 4f contribution.","tokens_in":10154,"tokens_out":3627,"would_cite":true,"duration_ms":34295,"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":"Density functional theory systematically overestimates the magnetic anisotropy of light rare-earths because its single Slater determinant distorts the 4f charge cloud; the paper shows an ion-only multipole correction fixes this, bringing…","keywords":["rare-earth magnets","magnetic anisotropy energy","4f charge asphericity","density functional theory","crystal field theory","many-body correction","SmCo5","light rare-earths"],"falsifier":"Measure the full magnetocrystalline anisotropy of NdCo5 at low temperature: the correction implies the in-plane term should be about three times smaller than uncorrected DFT predicts, so an experimentally robust sixfold term close to the uncorrected DFT value would disprove the ion-only ratios.","tokens_in":9030,"feed_emoji":"🧲","tokens_out":10832,"duration_ms":94367,"temperature":0.7,"pith_summary":"This paper identifies why density functional theory (DFT) fails for light rare-earth magnets: a single Slater determinant cannot represent the Hund's-rule $|J=L-S,J\\rangle$ ground state, so the $4f$ charge cloud appears more aspherical than it really is, and the magnetic anisotropy energy (MAE) comes out too large. The fix is a systematic many-body correction: each charge multipole moment $Q_k$ of the DFT density is rescaled by a ratio that depends only on the rare-earth ion, not on the material. The authors validate the underlying assumption that crystal-field parameters can be probed consistently inside DFT on TbV$_6$Sn$_6$ and TbCo$_5$, then apply the correction to SmCo$_5$, where the uniaxial MAE drops from roughly 25 meV/f.u. into the experimental range of 13–16 meV/f.u. This matters because light rare-earth permanent magnets such as SmCo$_5$ and Nd$_2$Fe$_{14}$B are the strongest permanent magnets known, and a computationally cheap correction opens them to systematic theory-driven design.","feed_headline":"One ion-only ratio brings SmCo5 magnetism in line with experiment","feed_subtitle":"DFT overstates the 4f charge shape; a per-ion rescaling recovers measured anisotropy.","key_machinery":"The load-bearing object is the set of ion-specific correction ratios $R_k = Q_k/Q_k^{\\mathrm{DFT}} = \\theta_J\\langle JJ|O_0^k|JJ\\rangle/(\\theta_L\\langle LL|O_0^k|LL\\rangle)$ for $k=2,4,6$, where $O_0^k$ are the operator-equivalent angular-momentum polynomials used in crystal-field theory and $\\theta_J$, $\\theta_L$ are the corresponding operator-equivalent coefficients. These ratios rescale the DFT-derived anisotropy coefficients $\\kappa_\\ell^m = A_\\ell^m Q_\\ell$ before the MAE is evaluated, so all material dependence is carried by the crystal-field parameters $A_\\ell^m$ extracted from constrained DFT+U calculations. The companion machinery is the rotation of the magnetization axis in constrained DFT+U with a fixed Hund's-rule occupancy, producing the energy surface $E(\\theta,\\phi)$ that is fit to the hexagonal $C_{3v}$ form of Eq. (3).","core_discovery":"The central claim is that the long-standing failure of DFT-based methods for light rare-earths has a specific, correctable origin: the single-Slater-determinant restriction of DFT represents the $4f$ shell by its maximally polarized $|L,L\\rangle$ component rather than the true $|J,J\\rangle$ ground state, so the angular anisotropy of the $4f$ charge density is systematically exaggerated for the $J=L-S$ light ions. The exaggeration is quantified by the multipole moments $Q_k$ ($k=2,4,6$), which are always overestimated in DFT; the correction factor is the ratio $Q_k/Q_k^{\\mathrm{DFT}}$, a pure number for each ion and multipole order. The paper shows that crystal-field parameters extracted from constrained DFT+U total-energy surfaces are essentially independent of which $4f$ configuration probes them in TbV$_6$Sn$_6$ and TbCo$_5$, giving a self-consistency check for combining DFT with crystal-field theory. With the correction, SmCo$_5$'s calculated uniaxial MAE falls from about 25 meV/f.u. to the experimentally measured 13–16 meV/f.u., and a spurious in-plane anisotropy disappears because $Q_6=0$ for the $J=5/2$ Sm ground state.","pith_inferences":["If the configuration-independence of the crystal-field parameters holds, the same correction ratios should predict how the anisotropy changes as the 4f occupancy or the J multiplet changes, allowing finite-temperature and doping trends to be computed from one set of crystal-field parameters.","The same overestimated charge asphericity should affect other orbital-sensitive observables in light rare-earths, such as magnetic form factors or X-ray magnetic circular dichroism; comparing those against experiment would test the correction beyond magnetocrystalline anisotropy.","The method's boundary is likely set by hybridization: for mixed-valent Ce or 5f actinides the assumption of a fixed crystal field probed by an integer-valence DFT state can break down, so a natural test is whether the extracted crystal-field parameters remain stable for a Ce compound with sizable f-f hopping."],"forward_implications":["DFT+U with the ion-only correction reproduces the measured MAE of SmCo5, making light rare-earth permanent magnets tractable with standard DFT codes.","The extracted crystal-field parameters can be reused for low-temperature properties and for excited spin-orbit multiplets, opening a route to optical-transition calculations.","Because the correction ratios depend only on the ion, the same tabulated factors apply across the isostructural rare-earth series without material-specific fitting.","In cubic Sm compounds the correction cuts the 4f anisotropy to about 15% of the uncorrected DFT value, which is where the leading anisotropy first appears.","For lighter light rare-earths such as Nd, the in-plane anisotropy is expected to be overestimated by about a factor of three before correction."],"supporting_citations":[{"why":"supplies the constrained DFT+U protocol that keeps the 4f shell in the correct Hund's-rule occupancy while the magnetization is rotated.","marker":"[12]"},{"why":"reports the uniaxial ferromagnetism of TbV6Sn6 used as a validation case.","marker":"[19]"},{"why":"provides the operator-equivalent formalism and multipole definitions underlying Eq. (4).","marker":"[23]"},{"why":"gives experimental anisotropic-magnetism data for TbV6Sn6 used for comparison during validation.","marker":"[33]"},{"why":"analyzes the rare-earth contribution to anisotropy in RCo5, the reference needed to interpret the TbCo5 sign change.","marker":"[34]"},{"why":"documents the magnetic properties of R ions in RCo5 used as a TbCo5 comparison.","marker":"[35]"},{"why":"provides the experimental SmCo5 MAE range of 13-16 meV/f.u. that the corrected calculation must match.","marker":"[36]"}],"fun_headline_variants":["DFT overestimates 4f charge shape; per-ion ratio fixes it","SmCo5 MAE corrected: 25 meV to 13-16 meV with simple ratio","Light rare-earth magnetism: DFT's 4f asphericity error fixed","Slater-determinant limit causes DFT error; ratio correction works","Single ratio aligns SmCo5 anisotropy with experiment"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The correction works only if the crystal-field potential seen by the 4f electrons is unchanged no matter which 4f orbital configuration is used to probe it in DFT, so that rescaling the charge asphericity alone captures the full many-body effect.","fun_headline_variants_meta":{"raw":{"variants":["DFT overestimates 4f charge shape; per-ion ratio fixes it","SmCo5 MAE corrected: 25 meV to 13-16 meV with simple ratio","Light rare-earth magnetism: DFT's 4f asphericity error fixed","Slater-determinant limit causes DFT error; ratio correction works","Single ratio aligns SmCo5 anisotropy with experiment"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001148,"raw_usage":{"total_tokens":4760,"prompt_tokens":946,"completion_tokens":3814,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":562,"completion_tokens_details":{"reasoning_tokens":3714}},"tokens_in":562,"tokens_out":3814,"duration_ms":26597,"temperature":1.0,"reasoning_tokens":3714,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T16:51:49.172149+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the full magnetocrystalline anisotropy of NdCo5 at low temperature: the correction implies the in-plane term should be about three times smaller than uncorrected DFT predicts, so an experimentally robust sixfold term close to the uncorrected DFT value would disprove the ion-only ratios.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the constrained DFT+U protocol that keeps the 4f shell in the correct Hund's-rule occupancy while the magnetization is rotated."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"provides the operator-equivalent formalism and multipole definitions underlying Eq. (4)."},{"cited_title":"Skomski, Simple Models of Magnetism , Oxford Grad- uate Texts (OUP Oxford, 2008)","cited_arxiv_id":null,"evidence_quote":"gives experimental anisotropic-magnetism data for TbV6Sn6 used for comparison during validation."},{"cited_title":"Hoffmann, Generalization of stevens’ operator- equivalent method, Journal of Physics A: Mathematical and General 24, 35 (1991)","cited_arxiv_id":null,"evidence_quote":"analyzes the rare-earth contribution to anisotropy in RCo5, the reference needed to interpret the TbCo5 sign change."},{"cited_title":"Duros, A","cited_arxiv_id":null,"evidence_quote":"documents the magnetic properties of R ions in RCo5 used as a TbCo5 comparison."}],"review_version":1}