{"id":"e09d6b8a-e11e-4d16-ae63-2e22c8046712","arxiv_id":"1908.04762","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":3.0,"correctness_risk":"high","formal_verification":"none","parameter_count":2,"one_line_summary":"LASSO regression on DFT data suggests a linear relation between MCA energy and orbital moment anisotropy in Fe-Rh/MgO films, but the paper's own analysis admits the prediction is poor.","lead":"This paper applies a machine learning method called LASSO to look for a link between magnetic anisotropy energy and orbital magnetic moments in iron-rhodium thin films. It claims to have found a linear relation, but its own results section says the prediction is not accurate.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The paper self-refutes its central claim: Section 3 states predicted EMCAs are 'not satisfying' and that orbital-moment descriptors are insufficient, so the claimed linear Bruno-type relation is unsupported.","rationale":"The reader's weakest-assumption diagnosis matches the load-bearing flaw: the sufficiency of orbital-moment descriptors is assumed in Eq. (1), and the paper's own Section 3 disowns that assumption. The strongest claim ('successfully found a linear behavior') is contradicted by the same section's admission that predictions are 'not satisfying' and that the breakdown of the Bruno relation may occur. Because the manuscript itself contains the decisive counter-evidence, and because no quantitative metrics or reproducible artifacts are provided, the REJECT verdict is appropriate. I agree with the reader's assessment; the only adjustment is to emphasize that the concern is not a subtle external risk but a direct internal inconsistency between the abstract/conclusion and the results section.","tokens_in":3693,"tokens_out":2126,"duration_ms":22823,"concrete_test":"Obtain the 64 DFT EMCA values and the corresponding layer-resolved orbital-moment anisotropies/averages used in Fig. 3; refit Eq. (1) by LASSO with the same leave-one-out cross-validation; report LOO RMSE, R2, and the RMSE of a null model that predicts the training mean for every test point. If the LOO RMSE is not substantially below the standard deviation of EMCA (or below roughly 0.1 meV/unit-area), the claimed linear behavior is not a valid predictor. If the data are not released, the claim remains unverifiable.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that LASSO finds a linear behavior between MCA energy and the anisotropy of orbital moments, implying that the Bruno relation holds for Fe-Rh/MgO films. For that claim to stand, the descriptor set in Eq. (1)---anisotropies and averages of orbital moments only---must be sufficient to determine EMCA. The manuscript's own Section 3 states the opposite: 'those predicted using the LASSO coefficients are not satisfying'; 'we considered only the orbital moment anisotropy and averages as the data descriptors for the MCA energy. The discrepancy between the first-principles calculated and the predicted EMCAs reveals that such descriptors might not be sufficient'; and 'resulting in a breakdown of the Bruno relation.' This is not an external objection but an explicit internal limitation. No quantitative error metric (RMSE, R2, MAE) is reported; Fig. 3c is described only narratively, with both 'agree well' and 'not satisfying' appearing in the same section. Leave-one-out cross-validation on 64 samples with roughly 28 descriptors cannot establish predictive power without error bars or comparison to a null model. The paper also provides no fitted coefficients, data, or code, so the claimed linear behavior is neither reproducible nor quantitatively supported. The load-bearing condition--descriptor sufficiency--is thus contradicted by the authors' own results, making the central claim unsupported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript applies LASSO regression to DFT-computed magnetocrystalline anisotropy (MCA) energies and orbital moments of Fe-Rh thin films on MgO(001), claiming that a linear relationship exists between the MCA energy and the anisotropy of orbital moments. The authors model 128 seven-layer binary Fe-Rh configurations, fit four coefficient sets (A, B, C, D) in Eq. (1) using leave-one-out cross-validation, and compare predicted with calculated EMCA values. The abstract and conclusions assert success, while the Results section states that the predicted values are 'not satisfying' and that the descriptors 'might not be sufficient', explicitly acknowledging a breakdown of the Bruno relation.","tokens_in":4009,"tokens_out":2551,"duration_ms":26619,"significance":"If the claimed linear relation were established with validated predictive accuracy, the work could provide a computationally cheap route to estimate MCA in Fe-Rh/MgO films from orbital moments alone. However, the manuscript reports no quantitative goodness-of-fit metrics, no fitted coefficients, no data, and no code, and its own Results section undermines the central claim. As presented, the paper does not constitute a verifiable scientific contribution beyond the underlying DFT data.","major_comments":[{"comment":"The central claim in the Abstract and Conclusions that 'we have successfully found a linear behavior' is directly contradicted by the same Section's statements: the predicted EMCAs 'are not satisfying', the orbital-moment descriptors 'might not be sufficient', and the analysis results in 'a breakdown of the Bruno relation'. The paper thus self-refutes its main assertion, and no evidence is offered to reconcile these contradictions.","section":"Section 3 (Results)"},{"comment":"No quantitative error metric (RMSE, MAE, R²) is reported for the LASSO predictions, and the same paragraph describes the comparison as both 'agree well' and 'not satisfying'. Without error bars or a comparison against a null model, leave-one-out cross-validation on 64 samples with roughly 28 descriptors cannot establish predictive power.","section":"Section 3, Fig. 3c"},{"comment":"The descriptor set consisting only of anisotropies and averages of orbital moments is explicitly admitted to be insufficient for systems with strong spin-orbit coupling. Because the linear relation is the paper's central assertion, this admission invalidates the load-bearing assumption of Eq. (1). A concrete test would be to report the fitted coefficients and evaluate predictions on a held-out set of configurations that were not used in fitting; neither is provided.","section":"Eq. (1) and following text"},{"comment":"The manuscript does not disclose the LASSO coefficients, the regularization parameter, the data, or the code. The claim that a large coefficient indicates physical importance cannot be checked, and the regression is not reproducible from the information given.","section":"Reproducibility"}],"minor_comments":[{"comment":"The title contains typos ('anisotripy', 'leaning' for 'anisotropy' and 'learning'), and the abstract has 'expectially' instead of 'especially'.","section":"Title and Abstract"},{"comment":"Figure numbering is duplicated and inconsistent: both Section 2 and Section 3 contain 'Figure 1' and the captions are reused, which makes the narrative difficult to follow.","section":"Figures"},{"comment":"Equation (1) is referenced but never displayed in the text, so the exact form of the regression model is unclear.","section":"Eq. (1)"},{"comment":"The text refers to 'AAu,1' in the coefficient explanation, but the system contains only Fe and Rh atoms, not Au.","section":"Notation"},{"comment":"The number of configurations is written as '27 = 128', which should be '2^7 = 128'.","section":"Section 2"},{"comment":"The reference list contains formatting errors, including missing closing brackets in [13] and [14], and a misspelled author name ('Tibshirani, Rober' should be 'Robert').","section":"References"}],"recommendation":"reject","confidential_remarks":"This manuscript is not in a publishable state: the authors' own Results section contradicts the abstract's central claim, and the absence of quantitative validation and missing regression coefficients preclude any independent assessment. The paper conflates in-sample fitting with prediction, and the internal admission of a breakdown of the Bruno relation is a fatal inconsistency. I recommend rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper's abstract says it 'successfully found a linear behavior between the MCA energy and the anisotropy of orbital moments' for Fe-Rh/MgO films, but Section 3 says the opposite: the LASSO-predicted EMCAs are 'not satisfying' and the chosen orbital-moment descriptors 'might not be sufficient.' The conclusion nevertheless claims 'reasonable agreement.' That internal contradiction alone sinks the central claim.\n\nWhat is actually worthwhile is the dataset. The authors computed MCA and orbital-moment anisotropies for all 128 possible atomic-layer configurations of a seven-layer Fe-Rh slab on MgO(001), using a well-converged FLAPW setup (3600 k-points, <0.01 meV numerical noise). That systematic enumeration is a legitimate contribution, and if the raw results were published, they could be a useful benchmark.\n\nThe soft spots are serious. There is no quantitative error metric anywhere—no R2, RMSE, or MAE—just a scatter plot with a red line. The LASSO coefficients are never given, so the model is not reproducible. LOOCV on 64 samples with roughly 28 descriptors cannot support a predictive claim without error bars or a null-model comparison. The paper also asserts a 'breakdown of the Bruno relation' in one paragraph and a linear Bruno-type relation in the abstract; those cannot both be true.\n\nThe citation pattern is not the main problem, though the heavy self-citation and the title typos ('anisotripy', 'leaning') suggest the manuscript is a rough draft. The real issue is that the authors' own results undercut the headline. This is not an external objection; it is written in the paper.\n\nWho is this for? A researcher working on Fe-Rh/MgO might appreciate the systematic DFT data, but the modeling conclusion is unsupported. I would not cite it in its current form. It does not deserve a serious referee; it deserves a desk reject with an invitation to resubmit if the authors provide the data, the coefficients, proper metrics, and a conclusion that matches their results.\n\nMy recommendation: reject, and don't spend referee time on this version.","headline":"The paper's own results section refutes its headline claim of a linear Bruno-type relation for Fe-Rh/MgO; only the systematic DFT dataset has value.","tokens_in":4470,"tokens_out":2473,"would_cite":false,"duration_ms":24675,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A LASSO linear model using orbital-moment anisotropies and averages predicts the magnetocrystalline anisotropy of Fe-Rh thin films on MgO(001).","keywords":["magnetocrystalline anisotropy","LASSO","orbital moment","Fe-Rh thin films","density functional theory","machine learning","MgO substrate","Bruno relation"],"falsifier":"Apply the identical LASSO descriptor set (orbital-moment anisotropies and averages) to a different strong-SOC binary thin film, such as Co-Pt or Fe-Pd on MgO, and compare the predicted EMCA against first-principles calculations; if the linear relation fails there, the descriptor set is not sufficient.","tokens_in":3534,"feed_emoji":"🧲","tokens_out":8463,"duration_ms":76038,"temperature":0.7,"pith_summary":"The paper applies LASSO regression to first-principles density-functional data for 128 atomic-layer configurations of seven-layer Fe-Rh films on MgO(001). It claims that the magnetocrystalline anisotropy energy (EMCA) is linearly correlated with the anisotropies and averages of the Fe and Rh orbital moments. The fitted model is tested with leave-one-out cross-validation and reported to agree with the DFT-computed EMCA values. The paper also includes a passage noting that the same descriptors may not be sufficient for systems with strong spin-orbit coupling, and that the discrepancy between calculated and predicted EMCAs suggests a breakdown of the Bruno relation.","feed_headline":"Orbital moments linearly predict Fe-Rh film anisotropy","feed_subtitle":"A LASSO model trained on density-functional data claims a simple linear rule for magnetocrystalline anisotropy in Fe-Rh/MgO films.","key_machinery":"The machinery is the linear regression model of Eq. (1), in which EMCA is expanded in the orbital-moment anisotropies and averages of Fe and Rh at each atomic layer, together with LASSO's L1 regularization for coefficient selection and leave-one-out cross-validation. The supporting DFT calculations use the FLAPW method with the force theorem to obtain EMCA and orbital moments for the 128 layer configurations. The 'Bruno relation' — the proportionality between EMCA and the anisotropy of the orbital moment — is the physical hypothesis the model tests.","core_discovery":"The paper's central claim is that LASSO regression on orbital-moment descriptors yields a linear relationship between EMCA and the anisotropy of orbital moments for Fe-Rh thin films, so that EMCA can be predicted from orbital moments alone. The authors assert this confirms the applicability of the Bruno relation to these binary films, where the strong spin-orbit coupling of Rh enhances the perpendicular anisotropy. The paper also contains an explicit caveat: the discrepancy between the first-principles and predicted EMCAs indicates that orbital-moment anisotropies and averages alone may not be sufficient descriptors, and that spin-flip terms can break the Bruno relation in strong-SOC systems.","pith_inferences":["A natural next step is to test the same descriptor set on other strong-SOC film systems, such as Co-Pt or Fe-Pd; if the linear relation fails there, the boundary of the Bruno relation's validity would be mapped.","The paper's own admission that the descriptors may be insufficient suggests that augmenting the model with spin-flip terms or interface-specific orbital-moment components could recover predictive power.","Because the dataset exhaustively covers all 128 layer orderings, it could be reused to benchmark richer descriptor sets and settle whether the observed linear trend is an artifact of collinearity in the orbital-moment features."],"forward_implications":["If the linear relation is correct, magnetocrystalline anisotropy energies for Fe-Rh/MgO films could be estimated directly from orbital-moment calculations, avoiding the more expensive force-theorem computations.","The fitted LASSO coefficients identify which atomic layers and which species contribute most to MCA, offering a possible design rule for tuning perpendicular anisotropy in Fe-Rh films.","The result supports the Bruno relation as a valid approximation for Fe-Rh binary films, extending its applicability to systems where 5d elements provide strong spin-orbit coupling.","The same regression procedure could be reused to screen other binary magnetic film configurations for high perpendicular MCA before performing full DFT calculations."],"supporting_citations":[{"why":"Introduces the LASSO estimator that the paper uses to select and fit the linear MCA descriptors.","marker":"[23]"},{"why":"Provides statistical guarantees for LASSO coefficient estimates, justifying its use on the small 128-configuration dataset.","marker":"[24]"},{"why":"Establishes that 4d and 5d elements with strong spin-orbit coupling can yield perpendicular magnetocrystalline anisotropy, motivating the Fe-Rh study.","marker":"[17]"},{"why":"Gives first-principles predictions of magnetocrystalline anisotropy in related ferromagnetic monolayers, informing the calculation setup.","marker":"[19]"},{"why":"Documents the strong volume magnetostriction and magnetic properties of FeRh, providing the material background for the film system.","marker":"[20]"}],"fun_headline_variants":["LASSO links orbital moments to Fe-Rh anisotropy","Fe-Rh MCA predicted from orbital moments via LASSO","Linear MCA rule from orbital moments in Fe-Rh films","Orbital-moment predictor for Fe-Rh anisotropy: linear but incomplete"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The model assumes that magnetocrystalline anisotropy energy is fully determined by a linear combination of the orbital-moment anisotropies and averages of Fe and Rh at each atomic layer, with no other physical ingredients.","fun_headline_variants_meta":{"raw":{"variants":["LASSO links orbital moments to Fe-Rh anisotropy","Fe-Rh MCA predicted from orbital moments via LASSO","Linear MCA rule from orbital moments in Fe-Rh films","Orbital-moment predictor for Fe-Rh anisotropy: linear but incomplete"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000237,"raw_usage":{"total_tokens":1436,"prompt_tokens":804,"completion_tokens":632,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":420,"completion_tokens_details":{"reasoning_tokens":560}},"tokens_in":420,"tokens_out":632,"duration_ms":6787,"temperature":1.0,"reasoning_tokens":560,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:56:14.805304+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Apply the identical LASSO descriptor set (orbital-moment anisotropies and averages) to a different strong-SOC binary thin film, such as Co-Pt or Fe-Pd on MgO, and compare the predicted EMCA against first-principles calculations; if the linear relation fails there, the descriptor set is not sufficient.","supporting_citations":[],"review_version":1}